In this paper, we construct a combined upwinding and mixed finite element method for the numerical solution of a two-dimensional mean field model of superconducting vortices. An advantage of our method is that it works for any unstructured regular triangulation. A simple convergence analysis is given without resorting to the discrete maximum principle. Numerical examples are also presented.
In this work, we address the problem of fluid-structure interaction (FSI) with moving structures that may come into contact. We propose a penalization contact algorithm implemented in an unfitted numerical framework designed to treat large displacements. In the proposed method, the fluid mesh is fixed and the structure meshes are superimposed to it without any constraint on the conformity. Thanks to the Extended Finite Element Method (XFEM), we can treat discontinuities of the fluid solution on...
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 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...
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...
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...
The authors study problems of existence and uniqueness of solutions of various variational formulations of the coupled problem of dynamical thermoelasticity and of the convergence of approximate solutions of these problems. First, the semidiscrete approximate solutions is defined, which is obtained by time discretization of the original variational problem by Euler’s backward formula. Under certain smoothness assumptions on the date authors prove existence and uniqueness of the solution and establish...
Compatible schemes localize degrees of freedom according to the physical nature of the underlying fields and operate a clear distinction between topological laws and closure relations. For elliptic problems, the cornerstone in the scheme design is the discrete Hodge operator linking gradients to fluxes by means of a dual mesh, while a structure-preserving discretization is employed for the gradient and divergence operators. The discrete Hodge operator is sparse, symmetric positive definite and is...
We analyze and suggest improvements to a recently developed approximate continuum-electrostatic model for proteins. The model, called BIBEE/I (boundary-integral based electrostatics estimation with interpolation), was able to estimate electrostatic solvation free energies to within a mean unsigned error of 4% on a test set of more than 600 proteins¶a significant improvement over previous BIBEE models. In this work, we tested the BIBEE/I model for its capability to predict residue-by-residue interactions...