3,475 research outputs found
Convergence of a cell-centered finite volume discretization for linear elasticity
We show convergence of a cell-centered finite volume discretization for
linear elasticity. The discretization, termed the MPSA method, was recently
proposed in the context of geological applications, where cell-centered
variables are often preferred. Our analysis utilizes a hybrid variational
formulation, which has previously been used to analyze finite volume
discretizations for the scalar diffusion equation. The current analysis
deviates significantly from previous in three respects. First, additional
stabilization leads to a more complex saddle-point problem. Secondly, a
discrete Korn's inequality has to be established for the global discretization.
Finally, robustness with respect to the Poisson ratio is analyzed. The
stability and convergence results presented herein provide the first rigorous
justification of the applicability of cell-centered finite volume methods to
problems in linear elasticity
Stable cell-centered finite volume discretization for Biot equations
In this paper we discuss a new discretization for the Biot equations. The
discretization treats the coupled system of deformation and flow directly, as
opposed to combining discretizations for the two separate sub-problems. The
coupled discretization has the following key properties, the combination of
which is novel: 1) The variables for the pressure and displacement are
co-located, and are as sparse as possible (e.g. one displacement vector and one
scalar pressure per cell center). 2) With locally computable restrictions on
grid types, the discretization is stable with respect to the limits of
incompressible fluid and small time-steps. 3) No artificial stabilization term
has been introduced. Furthermore, due to the finite volume structure embedded
in the discretization, explicit local expressions for both momentum-balancing
forces as well as mass-conservative fluid fluxes are available.
We prove stability of the proposed method with respect to all relevant
limits. Together with consistency, this proves convergence of the method.
Finally, we give numerical examples verifying both the analysis and convergence
of the method
Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity
Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and devise the hybridizable discontinuous Galerkin (HDG) method. Exact geometry described by non-uniform rational B-splines (NURBS) is integrated into HDG using the framework of the NURBS-enhanced finite element method (NEFEM). Moreover, optimal convergence and superconvergence properties of HDG-Voigt formulation in presence of symmetric second-order tensors are exploited to construct inexpensive error indicators and drive degree adaptive procedures. Applications involving the numerical simulation of problems in electrostatics, linear elasticity and incompressible viscous flows are presented. Moreover, this is done for both high-order HDG approximations and the lowest-order framework of face-centered finite volumes (FCFV).Peer ReviewedPostprint (author's final draft
Three-dimensional CFD simulations with large displacement of the geometries using a connectivity-change moving mesh approach
This paper deals with three-dimensional (3D) numerical simulations involving 3D moving geometries with large displacements on unstructured meshes. Such simulations are of great value to industry, but remain very time-consuming. A robust moving mesh algorithm coupling an elasticity-like mesh deformation solution and mesh optimizations was proposed in previous works, which removes the need for global remeshing when performing large displacements. The optimizations, and in particular generalized edge/face swapping, preserve the initial quality of the mesh throughout the simulation. We propose to integrate an Arbitrary Lagrangian Eulerian compressible flow solver into this process to demonstrate its capabilities in a full CFD computation context. This solver relies on a local enforcement of the discrete geometric conservation law to preserve the order of accuracy of the time integration. The displacement of the geometries is either imposed, or driven by fluid–structure interaction (FSI). In the latter case, the six degrees of freedom approach for rigid bodies is considered. Finally, several 3D imposed-motion and FSI examples are given to validate the proposed approach, both in academic and industrial configurations
- …