skip to main content


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
NSF-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. In this paper, we consider Maxwell’s equations in linear dispersive media described by a single-pole Lorentz model for electronic polarization. We study two classes of commonly used spatial discretizations: finite difference methods (FD) with arbitrary even order accuracy in space and high spatial order discontinuous Galerkin (DG) finite element methods. Both types of spatial discretizations are coupled with second order semi-implicit leap-frog and implicit trapezoidal temporal schemes. By performing detailed dispersion analysis for the semi-discrete and fully discrete schemes, we obtain rigorous quantification of the dispersion error for Lorentz dispersive dielectrics. In particular, comparisons of dispersion error can be made taking into account the model parameters, and mesh sizes in the design of the two types of schemes. This work is a continuation of our previous research on energy-stable numerical schemes for nonlinear dispersive optical media [6,7]. The results for the numerical dispersion analysis of the reduced linear model, considered in the present paper, can guide us in the optimal choice of discretization parameters for the more complicated and nonlinear models. The numerical dispersion analysis of the fully discrete FD and DG schemes, for the dispersive Maxwell model considered in this paper, clearly indicate the dependence of the numerical dispersion errors on spatial and temporal discretizations, their order of accuracy, mesh discretization parameters and model parameters. The results obtained here cannot be arrived at by considering discretizations of Maxwell’s equations in free space. In particular, our results contrast the advantages and disadvantages of using high order FD or DG schemes and leap-frog or trapezoidal time integrators over different frequency ranges using a variety of measures 
    more » « less
  2. Abstract

    Understanding spatial development of a turbulent mixing layer is essential for many engineering applications. However, the flow development is difficult to replicate in physical or numerical experiments. For this reason, the most attractive method for the mixing layer analysis is the direct numerical simulation (DNS), with the most control over the simulation inputs and free from modeling assumptions. On the other hand, the DNS cost often prevents conducting the sensitivity analysis of the simulation results to variations in the numerical procedure and thus, separating numerical and physical effects. In this paper, effects of the computational domain dimensions on statistics collected from DNS of a spatially developing incompressible turbulent mixing layer are analyzed with the focus on determining the domain dimensions suitable for studying the flow asymptotic state. In the simulations, the mixing layer develops between two coflowing laminar boundary layers formed on two sides of a sharp-ended splitter plate of a finite thickness with characteristics close to those of the untripped boundary layers in the experiments by Bell and Mehta, AIAA J., 28(12), 2034 (1990). The simulations were conducted using the spectral-element code Nek5000.

     
    more » « less
  3. Abstract

    Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations (PDEs). The Gauss-Seidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady-state solutions. A nice property of fixed-point fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems. Hence, they are easy to be applied to a general hyperbolic system. To deal with the difficulties associated with numerical boundary treatment when high-order finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains, inverse Lax-Wendroff (ILW) procedures were developed as a very effective approach in the literature. In this paper, we combine a fifth-order fixed-point fast sweeping WENO method with an ILW procedure to solve steady-state solution of hyperbolic conservation laws on complex computing regions. Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids. Numerical results show high-order accuracy and good performance of the method. Furthermore, the method is compared with the popular third-order total variation diminishing Runge-Kutta (TVD-RK3) time-marching method for steady-state computations. Numerical examples show that for most of examples, the fixed-point fast sweeping method saves more than half CPU time costs than TVD-RK3 to converge to steady-state solutions.

     
    more » « less
  4. This work is centered on high-fidelity modeling, analysis, and rigorous experiments of vibrations and guided (Lamb) waves in a human skull in two connected tracks: (1) layered modeling of the cranial bone structure (with cortical tables and diploë) and its vibration-based elastic parameter identification (and validation); (2) transcranial leaky Lamb wave characterization experiments and radiation analyses using the identified elastic parameters in a layered semi analytical finite element framework, followed by time transient simulations that consider the inner porosity as is. In the first track, non-contact vibration experiments are conducted to extract the first handful of modal frequencies in the auditory frequency regime, along with the associated damping ratios and mode shapes, of dry cranial bone segments extracted from the parietal and frontal regions of a human skull. Numerical models of the bone segments are built with a novel image reconstruction scheme that employs microcomputed tomographic scans to build a layered bone geometry with separate homogenized domains for the cortical tables and the diploë. These numerical models and the experimental modal frequencies are then used in an iterative parameter identification scheme that yields the cortical and diploic isotropic elastic moduli of each domain, whereas the corresponding densities are estimated using the total experimental mass and layer mass ratios obtained from the scans. With the identified elastic parameters, the average error between experimental and numerical modal frequencies is less than 1.5% and the modal assurance criterion values for most modes are above 0.90. Furthermore, the extracted parameters are in the range of the results reported in the literature. In the second track, the focus is placed on the subject of leaky Lamb waves, which has received growing attention as a promising alternative to conventional ultrasound techniques for transcranial transmission, especially to access the brain periphery. Experiments are conducted on the same cranial bone segment set for leaky Lamb wave excitation and radiation characterization. The degassed skull bone segments are used in submersed experiments with an ultrasonic transducer and needle hydrophone setup for radiation pressure field scanning. Elastic parameters obtained from the first track are used in guided wave dispersion simulations, and the radiation angles are accurately predicted using the aforementioned layered model in the presence of fluid loading. The dominant radiation angles are shown to correspond to guided wave modes with low attenuation and a significant out-of-plane polarization. The experimental radiation spectra are finally compared against those obtained from time transient finite element simulations that leverage geometric models reconstructed from microcomputed tomographic scans. 
    more » « less
  5. 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