Analysis for a molten carbonate fuel cell.
We present the convergence analysis of an efficient numerical method for the solution of an initial-boundary value problem for a scalar nonlinear conservation law equation with a diffusion term. Nonlinear convective terms are approximated with the aid of a monotone finite volume scheme considered over the finite volume barycentric mesh, whereas the diffusion term is discretized by piecewise linear nonconforming triangular finite elements. Under the assumption that the triangulations are of weakly...
A coupled finite/boundary element method to approximate the free vibration modes of an elastic structure containing an incompressible fluid is analyzed in this paper. The effect of the fluid is taken into account by means of one of the most usual procedures in engineering practice: an added mass formulation, which is posed in terms of boundary integral equations. Piecewise linear continuous elements are used to discretize the solid displacements and the fluid-solid interface variables. Spectral...
A coupled finite/boundary element method to approximate the free vibration modes of an elastic structure containing an incompressible fluid is analyzed in this paper. The effect of the fluid is taken into account by means of one of the most usual procedures in engineering practice: an added mass formulation, which is posed in terms of boundary integral equations. Piecewise linear continuous elements are used to discretize the solid displacements and the fluid-solid interface variables....
We analyze a force-based quasicontinuum approximation to a one-dimensional system of atoms that interact by a classical atomistic potential. This force-based quasicontinuum approximation can be derived as the modification of an energy-based quasicontinuum approximation by the addition of nonconservative forces to correct nonphysical “ghost” forces that occur in the atomistic to continuum interface during constant strain. The algorithmic simplicity and consistency with the purely atomistic model at...
We present a new stabilized mixed finite element method for the linear elasticity problem in . The approach is based on the introduction of Galerkin least-squares terms arising from the constitutive and equilibrium equations, and from the relation defining the rotation in terms of the displacement. We show that the resulting augmented variational formulation and the associated Galerkin scheme are well posed, and that the latter becomes locking-free and asymptotically locking-free for Dirichlet...
The smoothing-type algorithm is a powerful tool for solving the second-order cone programming (SOCP), which is in general designed based on a monotone line search. In this paper, we propose a smoothing-type algorithm for solving the SOCP with a non-monotone line search. By using the theory of Euclidean Jordan algebras, we prove that the proposed algorithm is globally and locally quadratically convergent under suitable assumptions. The preliminary numerical results are also reported which indicate...
We show that a non-standard mixed finite element method proposed by Barrios and Gatica in 2007, is a higher order perturbation of the least-squares mixed finite element method. Therefore, it is also superconvergent whenever the least-squares mixed finite element method is superconvergent. Superconvergence of the latter was earlier investigated by Brandts, Chen and Yang between 2004 and 2006. Since the new method leads to a non-symmetric system matrix, its application seems however more expensive...
In order to describe a solid which deforms smoothly in some region, but non smoothly in some other region, many multiscale methods have recently been proposed. They aim at coupling an atomistic model (discrete mechanics) with a macroscopic model (continuum mechanics). We provide here a theoretical ground for such a coupling in a one-dimensional setting. We briefly study the general case of a convex energy, and next concentrate on a specific example of a nonconvex energy, the Lennard-Jones case....
In order to describe a solid which deforms smoothly in some region, but non smoothly in some other region, many multiscale methods have recently been proposed. They aim at coupling an atomistic model (discrete mechanics) with a macroscopic model (continuum mechanics). We provide here a theoretical ground for such a coupling in a one-dimensional setting. We briefly study the general case of a convex energy, and next concentrate on a specific example of a nonconvex energy, the Lennard-Jones case....
The quasicontinuum method is a coarse-graining technique for reducing the complexity of atomistic simulations in a static and quasistatic setting. In this paper we aim to give a detailed a priori and a posteriori error analysis for a quasicontinuum method in one dimension. We consider atomistic models with Lennard–Jones type long-range interactions and a QC formulation which incorporates several important aspects of practical QC methods. First, we prove the existence, the local uniqueness...
We consider the spherically symmetric Vlasov-Einstein system in the case of asymptotically flat spacetimes. From the physical point of view this system of equations can model the formation of a spherical black hole by gravitational collapse or describe the evolution of galaxies and globular clusters. We present high-order numerical schemes based on semi-Lagrangian techniques. The convergence of the solution of the discretized problem to the exact solution is proven and high-order error estimates...
In this paper we propose a time discretization of a system of two parabolic equations describing diffusion-driven atom rearrangement in crystalline matter. The equations express the balances of microforces and microenergy; the two phase fields are the order parameter and the chemical potential. The initial and boundary-value problem for the evolutionary system is known to be well posed. Convergence of the discrete scheme to the solution of the continuous problem is proved by a careful development...