An embedded mesh method using piecewise constant multipliers with stabilization: mathematical and numerical aspects

Summary An embedded mesh method using piecewise constant multipliers originally proposed by Puso et al. (CMAME, 2012) is analyzed here to determine effects of the pressure stabilization term and small cut cells. The approach is implemented for transient dynamics using the central difference scheme for the time discretization. It is shown that the resulting equations of motion are a stable linear system with a condition number independent of mesh size. Next, it is shown that the constraints and the stabilization terms can be recast as non-proportional damping such that the time integration of the scheme is provably stable with a critical time step computed from the undamped equations of motion. Effects of small cuts are discussed throughout the presentation. A mesh study is conducted to evaluate the effects of the stabilization on the discretization error and conditioning and is used to recommend an optimal value for stabilization scaling parameter. Several nonlinear problems are also analyzed and compared with comparable conforming mesh results. Finally, several demanding problems highlighting the robustness of the proposed approach are shown. Copyright © 2014 John Wiley & Sons, Ltd.

[1]  P. Hansbo,et al.  A FINITE ELEMENT METHOD ON COMPOSITE GRIDS BASED ON NITSCHE'S METHOD , 2003 .

[2]  T. Hughes.  Generalization of selective integration procedures to anisotropic and nonlinear media , 1980 .

[3]  Peter Hansbo,et al.  Fictitious domain finite element methods using cut elements: I. A stabilized Lagrange multiplier method , 2010 .

[4]  Michael A. Puso,et al.  An embedded mesh method in a multiple material ALE , 2012 .

[5]  R. Glowinski,et al.  Error analysis of a fictitious domain method applied to a Dirichlet problem , 1995 .

[6]  Nicolas Moës,et al.  A stable Lagrange multiplier space for stiff interface conditions within the extended finite element method , 2009 .

[7]  Ralf Deiterding,et al.  A virtual test facility for the efficient simulation of solid material response under strong shock and detonation wave loading , 2006, Engineering with Computers.

[8]  C. Peskin.  Flow patterns around heart valves: A numerical method , 1972 .

[9]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[10]  R. Glowinski,et al.  A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow , 2001 .

[11]  Ted Belytschko,et al.  Fluid–structure interaction by the discontinuous‐Galerkin method for large deformations , 2009 .

[12]  Wolfgang A. Wall,et al.  An embedded Dirichlet formulation for 3D continua , 2010 .

[13]  Michael A. Puso,et al.  An embedded mesh method for treating overlapping finite element meshes , 2011 .

[14]  Kevin G. Wang,et al.  Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems , 2011 .

[15]  A. Lew,et al.  A discontinuous‐Galerkin‐based immersed boundary method , 2008 .

[16]  R. Taylor,et al.  Lagrange constraints for transient finite element surface contact , 1991 .

[17]  John E. Dolbow,et al.  Stable imposition of stiff constraints in explicit dynamics for embedded finite element methods , 2012 .