A Least Squares Method for Finding the Preisach Hysteresis Operator from Measurements.
We address in this article the computation of the convex solutions of the Dirichlet problem for the real elliptic Monge − Ampère equation for general convex domains in two dimensions. The method we discuss combines a least-squares formulation with a relaxation method. This approach leads to a sequence of Poisson − Dirichlet problems and another sequence of low dimensional algebraic eigenvalue problems of a new type. Mixed finite element approximations with a smoothing procedure are used for the...
A Legendre spectral collocation method is presented for the solution of the biharmonic Dirichlet problem on a square. The solution and its Laplacian are approximated using the set of basis functions suggested by Shen, which are linear combinations of Legendre polynomials. A Schur complement approach is used to reduce the resulting linear system to one involving the approximation of the Laplacian of the solution on the two vertical sides of the square. The Schur complement system is solved by a...
This paper examines a linear extrapolation time-discretization of a phase relaxation model with temperature dependent convection and reaction. The model consists of a diffusion-advection PDE for temperature and an ODE with double obstacle for phase variable. Under a stability constraint, this scheme is shown to converge with optimal orders for temperature and enthalpy, and for heat flux as time-step .
In this work we study a fully discrete mixed scheme, based on continuous finite elements in space and a linear semi-implicit first-order integration in time, approximating an Ericksen–Leslie nematic liquid crystal model by means of a Ginzburg–Landau penalized problem. Conditional stability of this scheme is proved via a discrete version of the energy law satisfied by the continuous problem, and conditional convergence towards generalized Young measure-valued solutions to the Ericksen–Leslie problem...
This paper proposes a linear discrete-time scheme for general nonlinear cross-diffusion systems. The scheme can be regarded as an extension of a linear scheme based on the nonlinear Chernoff formula for the degenerate parabolic equations, which proposed by Berger et al. [RAIRO Anal. Numer.13 (1979) 297–312]. We analyze stability and convergence of the linear scheme. To this end, we apply the theory of reaction-diffusion system approximation. After discretizing the scheme in space, we obtain a versatile,...
This paper proposes a linear discrete-time scheme for general nonlinear cross-diffusion systems. The scheme can be regarded as an extension of a linear scheme based on the nonlinear Chernoff formula for the degenerate parabolic equations, which proposed by Berger et al. [RAIRO Anal. Numer.13 (1979) 297–312]. We analyze stability and convergence of the linear scheme. To this end, we apply the theory of reaction-diffusion system approximation. After discretizing the scheme in space, we obtain a versatile,...
An extension of the local projection stabilization (LPS) finite element method for convection-diffusion-reaction equations is presented and analyzed, both in the steady-state and the transient setting. In addition to the standard LPS method, a nonlinear crosswind diffusion term is introduced that accounts for the reduction of spurious oscillations. The existence of a solution can be proved and, depending on the choice of the stabilization parameter, also its uniqueness. Error estimates are derived...
In this paper we propose and analyze a localized orthogonal decomposition (LOD) method for solving semi-linear elliptic problems with heterogeneous and highly variable coefficient functions. This Galerkin-type method is based on a generalized finite element basis that spans a low dimensional multiscale space. The basis is assembled by performing localized linear fine-scale computations on small patches that have a diameter of order H | log (H) | where H is the coarse mesh size. Without any assumptions...
The aim of this paper is to develop a finite element method which allows computing the buckling coefficients and modes of a non-homogeneous Timoshenko beam. Studying the spectral properties of a non-compact operator, we show that the relevant buckling coefficients correspond to isolated eigenvalues of finite multiplicity. Optimal order error estimates are proved for the eigenfunctions as well as a double order of convergence for the eigenvalues using classical abstract spectral approximation theory...
The aim of this paper is to develop a finite element method which allows computing the buckling coefficients and modes of a non-homogeneous Timoshenko beam. Studying the spectral properties of a non-compact operator, we show that the relevant buckling coefficients correspond to isolated eigenvalues of finite multiplicity. Optimal order error estimates are proved for the eigenfunctions as well as a double order of convergence for the eigenvalues using classical abstract spectral approximation theory...