skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: An Eulerian finite element method for PDEs in time-dependent domains
The paper introduces a new finite element numerical method for the solution of partial differential equations on evolving domains. The approach uses a completely Eulerian description of the domain motion.The physical domain is embedded in a triangulated computational domain and can overlap the time-independent background mesh in an arbitrary way. The numerical method is based on finite difference discretizations of time derivatives and a standard geometrically unfitted finite element method with an additional stabilization term in the spatial domain.The performance and analysis of the method rely on the fundamental extension result in Sobolev spaces for functions defined on bounded domains. This paper includes a complete stability and error analysis, which accounts for discretization errors resulting from finite difference and finite element approximations as well as for geometric errors coming from a possible approximate recovery of the physical domain. Several numerical examples illustrate the theory and demonstrate the practical efficiency of the method.  more » « less
Award ID(s):
1717516
PAR ID:
10090137
Author(s) / Creator(s):
;
Date Published:
Journal Name:
ESAIM: Mathematical Modelling and Numerical Analysis
ISSN:
0764-583X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Computational Methods in Applied Mathematics)
    Abstract The paper introduces a finite element method for an Eulerian formulation of partial differential equations governing the transport and diffusion of a scalar quantity in a time-dependent domain. The method follows the idea from[C. Lehrenfeld and M. Olshanskii,An Eulerian finite element method for PDEs in time-dependent domains,ESAIM Math. Model. Numer. Anal. 53 2019, 2, 585–614]of a solution extension to realise the Eulerian time-stepping scheme. However, a reformulation of the partial differential equation is suggested to derive a scheme which conserves the quantity under consideration exactly on the discrete level. For the spatial discretisation, the paper considers an unfitted finite element method. Ghost-penalty stabilisation is used to realise the discrete solution extension and gives a scheme robust against arbitrary intersections between the mesh and geometry interface. The stability is analysed for both first- and second-order backward differentiation formula versions of the scheme. Several numerical examples in two and three spatial dimensions are included to illustrate the potential of this method. 
    more » « less
  2. ; ; ; ; ; (, Geophysical Journal International)
    SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency. 
    more » « less
  3. ; ; (, Advances in applied mathematics and mechanics)
    Interface problems have wide applications in modern scientific research. Obtaining accurate numerical solutions of multi-domain problems involving triple junction conditions remains a significant challenge. In this paper, we develop an efficient finite element method based on non-body-fitting meshes for solving multi-domain elliptic interface problems. We follow the idea of immersed finite element by modifying local basis functions to accommodate interface conditions. We enrich the local finite element space by adding new basis functions for handling non-homogeneous flux jump. The numerical scheme is symmetric and positive definite. Numerical experiments are provided to demonstrate the features of our method. 
    more » « less
  4. ; ; (, IMA Journal of Numerical Analysis)
    Abstract In this paper we study the biharmonic equation with Navier boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the fourth-order problem as a system of Poisson equations. Our method differs from the naive mixed method that leads to two Poisson problems but only applies to convex domains; our decomposition involves a third Poisson equation to confine the solution in the correct function space, and therefore can be used in both convex and nonconvex domains. A $C^0$ finite element algorithm is in turn proposed to solve the resulting system. In addition, we derive optimal error estimates for the numerical solution on both quasi-uniform meshes and graded meshes. Numerical test results are presented to justify the theoretical findings. 
    more » « less
  5. ; (, IMA Journal of Numerical Analysis)
    Abstract The paper addresses an error analysis of an Eulerian finite element method used for solving a linearized Navier–Stokes problem in a time-dependent domain. In this study, the domain’s evolution is assumed to be known and independent of the solution to the problem at hand. The numerical method employed in the study combines a standard backward differentiation formula-type time-stepping procedure with a geometrically unfitted finite element discretization technique. Additionally, Nitsche’s method is utilized to enforce the boundary conditions. The paper presents a convergence estimate for several velocity–pressure elements that are inf-sup stable. The estimate demonstrates optimal order convergence in the energy norm for the velocity component and a scaled $$L^{2}(H^{1})$$-type norm for the pressure component. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email