An inverse problem formulation of the immersed‐boundary method

We formulate the immersed‐boundary method (IBM) as an inverse problem. A control variable is introduced on the boundary of a larger domain that encompasses the target domain. The optimal control is the one that minimizes the mismatch between the state and the desired boundary value along the immersed target‐domain boundary. We begin by investigating a naïve problem formulation that we show is ill‐posed: in the case of the Laplace equation, we prove that the solution is unique, but it fails to depend continuously on the data; for the linear advection equation, even solution uniqueness fails to hold. These issues are addressed by two complimentary strategies. The first strategy is to ensure that the enclosing domain tends to the true domain, as the mesh is refined. The second strategy is to include a specialized parameter‐free regularization that is based on penalizing the difference between the control and the state on the boundary. The proposed inverse IBM is applied to the diffusion, advection, and advection‐diffusion equations using a high‐order discontinuous Galerkin discretization. The numerical experiments demonstrate that the regularized scheme achieves optimal rates of convergence and that the reduced Hessian of the optimization problem has a bounded condition number, as the mesh is refined.

[1]  Jason E. Hicken,et al.  Immersed Boundary Method as an Inverse Problem , 2018, 2018 Fluid Dynamics Conference.

[2]  Pavel B. Bochev,et al.  An Optimization-Based Approach for Elliptic Problems with Interfaces , 2017, SIAM J. Sci. Comput..

[3]  Ken Mattsson,et al.  A high-order accurate embedded boundary method for first order hyperbolic equations , 2017, J. Comput. Phys..

[4]  Peter Hansbo,et al.  CutFEM: Discretizing geometry and partial differential equations , 2015 .

[5]  William Gropp,et al.  CFD Vision 2030 Study: A Path to Revolutionary Computational Aerosciences , 2014 .

[6]  P. G. Ciarlet,et al.  Linear and Nonlinear Functional Analysis with Applications , 2013 .

[7]  Hermann F. Fasel,et al.  A novel concept for the design of immersed interface methods , 2013, J. Comput. Phys..

[8]  Huafei Sun,et al.  A robust simplex cut-cell method for adaptive high-order discretizations of aerodynamics and multi-physics problems , 2013 .

[9]  John E. Dolbow,et al.  A robust Nitsche’s formulation for interface problems , 2012 .

[10]  Alfio Borzì,et al.  Computational Optimization of Systems Governed by Partial Differential Equations , 2012, Computational science and engineering.

[11]  Alexandre Noll Marques,et al.  A correction function method to solve incompressible fluid flows to high accuracy with immersed geometries , 2012 .

[12]  R. Glowinski,et al.  A Least-Squares/Fictitious Domain Method for Linear Elliptic Problems with Robin Boundary Conditions , 2011 .

[13]  Jean-Christophe Nave,et al.  A Correction Function Method for Poisson problems with interface jump conditions , 2010, J. Comput. Phys..

[14]  Fredi Tröltzsch,et al.  On Finite Element Error Estimates for Optimal Control Problems with Elliptic PDEs , 2009, LSSC.

[15]  Isaac Harari,et al.  An efficient finite element method for embedded interface problems , 2009 .

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

[17]  David L. Darmofal,et al.  A triangular cut-cell adaptive method for high-order discretizations of the compressible Navier-Stokes equations , 2007, J. Comput. Phys..

[18]  David E. Keyes,et al.  Parallel Algorithms for PDE-Constrained Optimization , 2006, Parallel Processing for Scientific Computing.

[19]  George Biros,et al.  Parallel Lagrange-Newton-Krylov-Schur Methods for PDE-Constrained Optimization. Part I: The Krylov-Schur Solver , 2005, SIAM J. Sci. Comput..

[20]  Khosro Shahbazi,et al.  An explicit expression for the penalty parameter of the interior penalty method , 2022 .

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

[22]  Otmar Scherzer,et al.  The use of Morozov's discrepancy principle for Tikhonov regularization for solving nonlinear ill-posed problems , 1993, Computing.

[23]  P. Hansbo,et al.  An unfitted finite element method, based on Nitsche's method, for elliptic interface problems , 2002 .

[24]  P. Roache.  Code Verification by the Method of Manufactured Solutions , 2002 .

[25]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[26]  Paul Houston,et al.  Discontinuous hp-Finite Element Methods for Advection-Diffusion-Reaction Problems , 2001, SIAM J. Numer. Anal..

[27]  A. Polyanin.  Handbook of Linear Partial Differential Equations for Engineers and Scientists , 2001 .

[28]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[29]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[30]  W. McLean.  Strongly Elliptic Systems and Boundary Integral Equations , 2000 .

[31]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

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

[33]  M. Berger,et al.  Robust and efficient Cartesian mesh generation for component-based geometry , 1998 .

[34]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[35]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[36]  S. Axler,et al.  Harmonic Function Theory , 1992 .

[37]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[38]  Anita Mayo,et al.  The rapid evaluation of volume integrals of potential theory on general regions , 1992 .

[39]  R. LeVeque,et al.  An Adaptive Cartesian Mesh Algorithm for the Euler Equations in Arbitrary Geometries , 1989 .

[40]  P. Grisvard.  Elliptic Problems in Nonsmooth Domains , 1985 .

[41]  A. Mayo.  The Fast Solution of Poisson’s and the Biharmonic Equations on Irregular Regions , 1984 .

[42]  D. Arnold.  An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[43]  J. W. Purvis,et al.  Prediction of critical Mach number for store configurations , 1979 .

[44]  C. Peskin.  Numerical analysis of blood flow in the heart , 1977 .

[45]  Hermann Weyl,et al.  The method of orthogonal projection in potential theory , 1940 .