Applications of Mathematics最新文献

We propose some approaches for the generation of conforming simplicial partitions with various regularity properties for polytopic domains that are products or a union of products, thus generalizing our earlier results. The techniques presented can be used for finite element simulations of higher-dimensional problems.

Conjugate gradient methods are widely used for solving large-scale unconstrained optimization problems, because they do not need the storage of matrices. Based on the self-scaling memoryless Broyden-Fletcher-Goldfarb-Shanno (SSML-BFGS) method, new conjugate gradient algorithms CG-DESCENT and CGOPT have been proposed by W. Hager, H. Zhang (2005) and Y. Dai, C. Kou (2013), respectively. It is noted that the two conjugate gradient methods perform more efficiently than the SSML-BFGS method. Therefore, C. Kou, Y. Dai (2015) proposed some suitable modifications of the SSML-BFGS method such that the sufficient descent condition holds. For the sake of improvement of modified SSML-BFGS method, in this paper, we present an efficient SSML-BFGS-type three-term conjugate gradient method for solving unconstrained minimization using Ford-Moghrabi secant equation instead of the usual secant equations. The method is shown to be globally convergent under certain assumptions. Numerical results compared with methods using the usual secant equations are reported.

We introduce a new scaling parameter for the Dai-Kou family of conjugate gradient algorithms (2013), which is one of the most numerically efficient methods for unconstrained optimization. The suggested parameter is based on eigenvalue analysis of the search direction matrix and minimizing the measure function defined by Dennis and Wolkowicz (1993). The corresponding search direction of conjugate gradient method has the sufficient descent property and the extended conjugacy condition. The global convergence of the proposed algorithm is given for both uniformly convex and general nonlinear objective functions. Also, numerical experiments on a set of test functions of the CUTER collections and the practical problem of the manipulator of robot movement control show that the proposed method is effective.

We consider pressure-driven flow between adjacent surfaces, where the fluid is assumed to have constant density. The main novelty lies in using implicit algebraic constitutive relations to describe the fluid’s response to external stimuli, enabling the modeling of fluids whose responses cannot be accurately captured by conventional methods. When the implicit algebraic constitutive relations cannot be solved for the Cauchy stress in terms of the symmetric part of the velocity gradient, the traditional approach of inserting the expression for the Cauchy stress into the equation for the balance of linear momentum to derive the governing equation for the velocity becomes inapplicable. Instead, a non-standard system of first-order equations governs the flow. This system is highly complex, making it important to develop simplified models. Our primary contribution is the development of a framework for achieving this. Additionally, we apply our findings to a fluid that exhibits an S-shaped curve in the shear stress versus shear rate plot, as observed in some colloidal solutions.

We present a precise anisotropic interpolation error estimate for the Morley finite element method (FEM) and apply it to fourth-order elliptic equations. We do not impose the shape-regularity mesh condition in the analysis. Anisotropic meshes can be used for this purpose. The main contributions of this study include providing a new proof of the term consistency. This enables us to obtain an anisotropic consistency error estimate. The core idea of the proof involves using the relationship between the Raviart-Thomas and Morley finite-element spaces. Our results indicate optimal convergence rates and imply that the modified Morley FEM may be effective for errors.

In an error estimation of finite element solutions to the Poisson equation, we usually impose the shape regularity assumption on the meshes to be used. In this paper, we show that even if the shape regularity condition is violated, the standard error estimation can be obtained if “bad” elements that violate the shape regularity or maximum angle condition are covered virtually by simplices that satisfy the minimum angle condition. A numerical experiment illustrates the theoretical result.

We construct a new initial data to prove the ill-posedness of both Navier-Stokes and Euler equations in weaker Besov spaces in the sense that the solution maps to these equations starting from u₀ are discontinuous at t = 0.

A model two-dimensional acoustic waveguide with lateral impedance boundary conditions (and outgoing boundary conditions at the waveguide outlet) is considered. The governing operator is proved to be bounded below with a stability constant inversely proportional to the length of the waveguide. The presence of impedance boundary conditions leads to a non self-adjoint operator which considerably complicates the analysis. The goal of this paper is to elucidate these complications and tools that are applicable, as simply as possible. This work is a continuation of prior waveguide studies (where self-adjoint operators arose) by J. M. Melenk et al. (2023), and L. Demkowicz et al. (2024).

We deal with the a-posteriori estimation of the error for finite element solutions of nonstationary heat conduction problems with mixed boundary conditions on bounded polygonal domains. The a-posteriori error estimates are constucted by solving stationary “reconstruction” problems, obtained by replacing the time derivative of the exact solution by the time derivative of the finite element solution. The main result is that the reconstructed solutions, or reconstructions, are superconvergent approximations of the exact solution (they are more accurate than the finite element solution) when the error is measured in the gradient or the energy-norm. Because of this, the error in the gradient of the finite element solution can be estimated reliably, by computing its difference from the gradient of its reconstructions. Numerical examples show that “reconstruction estimates” are reliable for the most general classes of solutions which can occur in practical computations.