Eigensolution analysis of immersed boundary method based on volume penalization: applications to high-order schemes

This paper presents eigensolution and non-modal analyses for immersed boundary methods (IBMs) based on volume penalization for the linear advection equation. This approach is used to analyze the behavior of flux reconstruction (FR) discretization, including the influence of polynomial order and penalization parameter on numerical errors and stability. Through a semi-discrete analysis, we find that the inclusion of IBM adds additional dissipation without changing significantly the dispersion of the original numerical discretization. This agrees with the physical intuition that in this type of approach, the solid wall is modelled as a porous medium with vanishing viscosity. From a stability point view, the selection of penalty parameter can be analyzed based on a fully-discrete analysis, which leads to practical guidelines on the selection of penalization parameter. Numerical experiments indicate that the penalization term needs to be increased to damp oscillations inside the solid (i.e. porous region), which leads to undesirable time step restrictions. As an alternative, we propose to include a second-order term in the solid for the no-slip wall boundary condition. Results show that by adding a second-order term we improve the overall accuracy with relaxed time step restriction. This indicates that the optimal value of the penalization parameter and the second-order damping can be carefully chosen to obtain a more accurate scheme. Finally, the approximated relationship between these two parameters is obtained and used as a guideline to select the optimum penalty terms in a Navier-Stokes solver, to simulate flow past a cylinder. © 2021 Elsevier Inc. All rights reserved.

[1]  Philippe Angot,et al.  A penalization method to take into account obstacles in incompressible viscous flows , 1999, Numerische Mathematik.

[2]  Thomas de Quincey.  [C] , 2000, The Works of Thomas De Quincey, Vol. 1: Writings, 1799–1820.

[3]  J. Wu,et al.  Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications , 2009 .

[4]  B. Wetton,et al.  Analysis of Stiffness in the Immersed Boundary Method and Implications for Time-Stepping Schemes , 1999 .

[5]  Antony Jameson,et al.  Energy stable flux reconstruction schemes for advection-diffusion problems on triangles , 2013, J. Comput. Phys..

[6]  Spencer J. Sherwin,et al.  Spatial eigensolution analysis of discontinuous Galerkin schemes with practical insights for under-resolved computations and implicit LES , 2017, Computers & Fluids.

[7]  Gianmarco Mengaldo,et al.  Non-modal analysis of spectral element methods: Towards accurate and robust large-eddy simulations , 2018, Computer Methods in Applied Mechanics and Engineering.

[8]  K. Schneider,et al.  Analysis and discretization of the volume penalized Laplace operator with Neumann boundary conditions , 2014, 1403.5948.

[9]  Do Y. Kwak,et al.  Immersed finite element method for eigenvalue problem , 2014, J. Comput. Appl. Math..

[10]  X. Yao,et al.  A coupled Volume Penalization-Thermal Lattice Boltzmann method for thermal flows , 2018, International Journal of Heat and Mass Transfer.

[11]  H. T. Huynh,et al.  A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods , 2007 .

[12]  Gregor Gassner,et al.  A Comparison of the Dispersion and Dissipation Errors of Gauss and Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods , 2011, SIAM J. Sci. Comput..

[13]  Huaxiong Huang,et al.  Stability analysis of the immersed boundary method for a two-dimensional membrane with bending rigidity , 2008 .

[14]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[15]  Spencer J. Sherwin,et al.  Linear dispersion-diffusion analysis and its application to under-resolved turbulence simulations using discontinuous Galerkin spectral/hp methods , 2015, J. Comput. Phys..

[16]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[17]  Sylvain Laizet,et al.  A DNS study of jet control with microjets using an immersed boundary method , 2014 .

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

[19]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[20]  M. Carpenter,et al.  Fourth-order 2N-storage Runge-Kutta schemes , 1994 .

[21]  Gorjan Alagic,et al.  #p , 2019, Quantum information & computation.

[22]  Frederick Stern,et al.  A simple and efficient direct forcing immersed boundary framework for fluid-structure interactions , 2012, J. Comput. Phys..

[23]  P. Alam.  ‘A’ , 2021, Composites Engineering: An A–Z Guide.

[24]  He Yang,et al.  Dispersion and Dissipation Errors of Two Fully Discrete Discontinuous Galerkin Methods , 2012, Journal of Scientific Computing.

[25]  G. Iaccarino,et al.  Immersed boundary technique for turbulent flow simulations , 2003 .

[26]  B. Yin,et al.  On the numerical oscillation of the direct-forcing immersed-boundary method for moving boundaries , 2012 .

[27]  Bernardo Cockburn,et al.  The Runge-Kutta local projection discontinous Galerkin finite element method for conservation laws , 1990 .

[28]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[29]  E. Ferrer,et al.  Design of a Smagorinsky spectral Vanishing Viscosity turbulence model for discontinuous Galerkin methods , 2020 .

[30]  Li Wang,et al.  An immersed boundary method for fluid-structure interaction with compressible multiphase flows , 2017, J. Comput. Phys..

[31]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[32]  Philippe Angot,et al.  A volume penalization method for incompressible flows and scalar advection-diffusion with moving obstacles , 2011, J. Comput. Phys..

[33]  Hyung Jin Sung,et al.  Simulation of flexible filaments in a uniform flow by the immersed boundary method , 2007, J. Comput. Phys..

[34]  Gilles Carbou,et al.  Boundary layer for a penalization method for viscous incompressible flow , 2003, Advances in Differential Equations.

[35]  John H. Kolias,et al.  A CONSERVATIVE STAGGERED-GRID CHEBYSHEV MULTIDOMAIN METHOD FOR COMPRESSIBLE FLOWS , 1995 .

[36]  Dmitry Kolomenskiy,et al.  Numerical simulation of fluid-structure interaction with the volume penalization method , 2015, J. Comput. Phys..

[37]  V. Guinot.  Approximate Riemann Solvers , 2010 .

[38]  Pietro De Palma,et al.  An immersed boundary method for compressible flows using local grid refinement , 2007, J. Comput. Phys..

[39]  Will Trojak,et al.  Effect of Mesh Quality on Flux Reconstruction in Multi-dimensions , 2018, J. Sci. Comput..

[40]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[41]  Rémi Abgrall,et al.  An immersed boundary method using unstructured anisotropic mesh adaptation combined with level-sets and penalization techniques , 2014, J. Comput. Phys..

[42]  W. Shyy,et al.  Regular Article: An Accurate Cartesian Grid Method for Viscous Incompressible Flows with Complex Immersed Boundaries , 1999 .

[43]  Zhi Jian Wang,et al.  Fourier analysis and evaluation of DG, FD and compact difference methods for conservation laws , 2018, J. Comput. Phys..

[44]  Dmitry Kolomenskiy,et al.  A Fourier spectral method for the Navier-Stokes equations with volume penalization for moving solid obstacles , 2009, J. Comput. Phys..

[45]  Wei Guo,et al.  Superconvergence of discontinuous Galerkin and local discontinuous Galerkin methods: Eigen-structure analysis based on Fourier approach , 2013, J. Comput. Phys..

[46]  M. Quintard,et al.  A penalization technique applied to the “Volume-Of-Fluid” method: wettability condition on immersed boundaries , 2014 .

[47]  Antony Jameson,et al.  A New Class of High-Order Energy Stable Flux Reconstruction Schemes , 2011, J. Sci. Comput..

[48]  Lilia Krivodonova,et al.  Relaxing the CFL Number of the Discontinuous Galerkin Method , 2014, SIAM J. Sci. Comput..

[49]  Tim Colonius,et al.  The immersed boundary method: A projection approach , 2007, J. Comput. Phys..

[50]  Spencer J. Sherwin,et al.  Spatial eigensolution analysis of energy-stable flux reconstruction schemes and influence of the numerical flux on accuracy and robustness , 2018, J. Comput. Phys..

[51]  Yoshiharu Tamaki,et al.  Near-Wall Modification of Spalart–Allmaras Turbulence Model for Immersed Boundary Method , 2016 .

[52]  H. S. Udaykumar,et al.  A Sharp Interface Cartesian Grid Methodfor Simulating Flows with ComplexMoving Boundaries , 2001 .

[53]  Tapan K. Sengupta,et al.  Global spectral analysis for convection-diffusion-reaction equation in one and two-dimensions: Effects of numerical anti-diffusion and dispersion , 2020, J. Comput. Phys..

[54]  Kai Schneider,et al.  Volume penalization for inhomogeneous Neumann boundary conditions modeling scalar flux in complicated geometry , 2019, J. Comput. Phys..

[55]  G. Karniadakis,et al.  Spectral/hp Element Methods for Computational Fluid Dynamics , 2005 .

[56]  Antony Jameson,et al.  Insights from von Neumann analysis of high-order flux reconstruction schemes , 2011, J. Comput. Phys..

[57]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[58]  S. Balachandar,et al.  An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing , 2021, J. Comput. Phys..

[59]  Krishnan Mahesh,et al.  High order finite difference schemes with good spectral resolution , 1997 .

[60]  Kai Schneider,et al.  Immersed boundary methods for numerical simulation of confined fluid and plasma turbulence in complex geometries: a review , 2015, 1508.04593.

[61]  M. Y. Hussaini,et al.  An Analysis of the Discontinuous Galerkin Method for Wave Propagation Problems , 1999 .

[62]  Xiaodong Jing,et al.  An immersed boundary computational model for acoustic scattering problems with complex geometries. , 2012, The Journal of the Acoustical Society of America.

[63]  P. Alam.  ‘K’ , 2021, Composites Engineering.

[64]  M. Bussmann,et al.  A volume-of-fluid ghost-cell immersed boundary method for multiphase flows with contact line dynamics , 2018 .

[65]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[66]  L. Sirovich,et al.  Modeling a no-slip flow boundary with an external force field , 1993 .

[67]  Jung Hee Seo,et al.  A high-order immersed boundary method for acoustic wave scattering and low-Mach number flow-induced sound in complex geometries , 2011, J. Comput. Phys..

[68]  Zhi Jian Wang,et al.  A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids , 2009, J. Comput. Phys..

[69]  Soledad Le Clainche,et al.  Data-driven eigensolution analysis based on a spatio-temporal Koopman decomposition, with applications to high-order methods , 2022, J. Comput. Phys..

[70]  M. Chávez-Modena,et al.  Optimizing free parameters in the D3Q19 Multiple-Relaxation lattice Boltzmann methods to simulate under-resolved turbulent flows , 2020, J. Comput. Sci..

[71]  Marcel Vinokur,et al.  Spectral difference method for unstructured grids I: Basic formulation , 2006, J. Comput. Phys..

[72]  Spencer J. Sherwin,et al.  Eigensolution analysis of spectral/hp continuous Galerkin approximations to advection-diffusion problems: Insights into spectral vanishing viscosity , 2016, J. Comput. Phys..

[73]  P. Alam,et al.  H , 1887, High Explosives, Propellants, Pyrotechnics.

[74]  Esteban Ferrer,et al.  Dispersion-Dissipation Analysis for Advection Problems with Nonconstant Coefficients: Applications to Discontinuous Galerkin Formulations , 2018, SIAM J. Sci. Comput..

[75]  Xijun He,et al.  Dispersion-dissipation analysis of triangular numerical-flux-based discontinuous Galerkin method for elastic wave equations , 2020, J. Comput. Phys..

[76]  Eric Brown-Dymkoski,et al.  A characteristic based volume penalization method for general evolution problems applied to compressible viscous flows , 2014, J. Comput. Phys..

[77]  Hu Dai,et al.  Fluid-structure interaction involving large deformations: 3D simulations and applications to biological systems , 2014, J. Comput. Phys..

[78]  Charles S. Peskin,et al.  Stability and Instability in the Computation of Flows with Moving Immersed Boundaries: A Comparison of Three Methods , 1992, SIAM J. Sci. Comput..

[79]  Brian C. Vermeire,et al.  On the behaviour of fully-discrete flux reconstruction schemes , 2017 .

[80]  High-Order Flux Reconstruction Based on Immersed Boundary Method , 2021, 14th WCCM-ECCOMAS Congress.

[81]  F. Sotiropoulos,et al.  Immersed boundary methods for simulating fluid-structure interaction , 2014 .

[82]  Kai Schneider,et al.  Approximation of the Laplace and Stokes operators with Dirichlet boundary conditions through volume penalization: a spectral viewpoint , 2012, Numerische Mathematik.

[83]  G. Rubio,et al.  Improving the stability of multiple-relaxation lattice Boltzmann methods with central moments , 2018, Computers & Fluids.

[84]  Boyce E. Griffith,et al.  Immersed Methods for Fluid-Structure Interaction. , 2020, Annual review of fluid mechanics.

[85]  Adrián Navas-Montilla,et al.  Application of approximate dispersion-diffusion analyses to under-resolved Burgers turbulence using high resolution WENO and UWC schemes , 2021, J. Comput. Phys..

[86]  Chris Lacor,et al.  Dispersion and dissipation properties of the 1D spectral volume method and application to a p-multigrid algorithm , 2007, J. Comput. Phys..

[87]  Robert D. Guy,et al.  Immersed boundary smooth extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods , 2015, J. Comput. Phys..

[88]  Esteban Ferrer,et al.  Immersed boundary method for high-order flux reconstruction based on volume penalization , 2022, J. Comput. Phys..