5,882 research outputs found
High-order numerical methods for 2D parabolic problems in single and composite domains
In this work, we discuss and compare three methods for the numerical
approximation of constant- and variable-coefficient diffusion equations in both
single and composite domains with possible discontinuity in the solution/flux
at interfaces, considering (i) the Cut Finite Element Method; (ii) the
Difference Potentials Method; and (iii) the summation-by-parts Finite
Difference Method. First we give a brief introduction for each of the three
methods. Next, we propose benchmark problems, and consider numerical tests-with
respect to accuracy and convergence-for linear parabolic problems on a single
domain, and continue with similar tests for linear parabolic problems on a
composite domain (with the interface defined either explicitly or implicitly).
Lastly, a comparative discussion of the methods and numerical results will be
given.Comment: 45 pages, 12 figures, in revision for Journal of Scientific Computin
Accurate gradient computations at interfaces using finite element methods
New finite element methods are proposed for elliptic interface problems in
one and two dimensions. The main motivation is not only to get an accurate
solution but also an accurate first order derivative at the interface (from
each side). The key in 1D is to use the idea from \cite{wheeler1974galerkin}.
For 2D interface problems, the idea is to introduce a small tube near the
interface and introduce the gradient as part of unknowns, which is similar to a
mixed finite element method, except only at the interface. Thus the
computational cost is just slightly higher than the standard finite element
method. We present rigorous one dimensional analysis, which show second order
convergence order for both of the solution and the gradient in 1D. For two
dimensional problems, we present numerical results and observe second order
convergence for the solution, and super-convergence for the gradient at the
interface
Finite Element Analysis of an Arbitrary Lagrangian–Eulerian Method for Stokes/Parabolic Moving Interface Problem With Jump Coefficients
In this paper, a type of arbitrary Lagrangian–Eulerian (ALE) finite element method in the monolithic frame is developed for a linearized fluid–structure interaction (FSI) problem — an unsteady Stokes/parabolic interface problem with jump coefficients and moving interface, where, the corresponding mixed finite element approximation in both semi- and fully discrete scheme are developed and analyzed based upon one type of ALE formulation and a novel H1- projection technique associated with a moving interface problem, and the stability and optimal convergence properties in the energy norm are obtained for both discretizations to approximate the solution of a transient Stokes/parabolic interface problem that is equipped with a low regularity. Numerical experiments further validate all theoretical results. The developed analytical approaches and numerical implementations can be similarly extended to a realistic FSI problem in the future
Parallel numerical modeling of hybrid-dimensional compositional non-isothermal Darcy flows in fractured porous media
This paper introduces a new discrete fracture model accounting for
non-isothermal compositional multiphase Darcy flows and complex networks of
fractures with intersecting, immersed and non immersed fractures. The so called
hybrid-dimensional model using a 2D model in the fractures coupled with a 3D
model in the matrix is first derived rigorously starting from the
equi-dimensional matrix fracture model. Then, it is dis-cretized using a fully
implicit time integration combined with the Vertex Approximate Gradient (VAG)
finite volume scheme which is adapted to polyhedral meshes and anisotropic
heterogeneous media. The fully coupled systems are assembled and solved in
parallel using the Single Program Multiple Data (SPMD) paradigm with one layer
of ghost cells. This strategy allows for a local assembly of the discrete
systems. An efficient preconditioner is implemented to solve the linear systems
at each time step and each Newton type iteration of the simulation. The
numerical efficiency of our approach is assessed on different meshes, fracture
networks, and physical settings in terms of parallel scalability, nonlinear
convergence and linear convergence
- …