Fictitious domain method

Short description: Numerical analysis method

In numerical analysis, the fictitious domain method (FD) is a numerical technique designed to solve partial differential equations on complex geometries by embedding the physical domain into a larger and simpler computational domain. The method consists of extending the equations beyond the physical boundaries and enforcing the interface conditions through a distributed Lagrange multiplier in order to recover the correct solution within the original domain of interest.[1]

Example of a fitted mesh for a circular obstacle Example of a background unfitted mesh with an overlapping circular obstacle

This method belongs to the more general family of unfitted methods (also known as embedded or immersed), which allow solving interface problems on complex or evolving domains without generating a mesh that conforms to the domain’s boundaries.[2] For this reason, the construction of two independent meshes is considered in the fictitious domain method: one for the unfitted background domain and another overlapping mesh that conforms to the interface, but which can move freely during the simulation. As with the immersed boundary method[3] and the level-set method,[4] the fictitious domain method is particularly effective for problems with evolving geometries since there is no need to rebuild the mesh but just to move the overlapping one.[5]

The Dirichlet boundary conditions on immersed closed interfaces are imposed weakly[6] by the introduction of the distributed Lagrange multiplier which enforces the boundary condition over the whole overlapping region.[1]

A particular application of the fictitious domain method is in fluid–structure interaction (FSI) problems, where a fluid and an immersed solid interact one each other, exchanging forces through the interface.[7] In this context, the solid is considered to be embedded within the computational fluid dynamics domain and the coupling conditions at the interface are enforced weakly over the whole fictitious fluid region.[8]

## Contents

* 1 Weak imposition of Dirichlet boundary conditions
* 2 Fluid-structure interaction application
* 3 Numerical aspects
* 4 See also
* 5 References
* 6 External links

## Weak imposition of Dirichlet boundary conditions

Let Ωf⊂ℝd be a fluid domain in which we want to solve the Stokes equation for incompressible Newtonian fluids. We want to impose the Dirichlet boundary condition u=g on the boundary of an immersed circle described by the interface Γ=∂Ωc. In this context, the fictitious domain method allows us to solve the Stokes equation over the whole domain Ω=Ωf∪Ωc, and to impose weakly the constraint u=g by means of the Lagrange multiplier λ.

Velocity field around a fixed circular obstacle

The fictitious domain weak formulation reads:[9]{∫Ωμ∇u:∇v−∫Ωp∇⋅v+⟨λ,v⟩Ωc=0∀v∈H01(Ω)∫Ωq ∇⋅uf=0∀q∈L02(Ω)⟨u−g,η⟩Ωc=0∀η∈Λ

The last equation of the system guarantees that, in a weak sense, the interface condition u=g is satisfied over the whole domain Ωc, and in particular on the interface Γ, while in Ωf the equation is left untouched.[6]

The Lagrange multiplier λ acts on the fluid as an external force defined only in Ωc that enforces the fictitious fluid velocity u to fulfil the constraint inside the circle.[8]

## Fluid-structure interaction application

As it is typically done in fluid-structure interaction problems, we consider a fluid domain Ωf⊂ℝd, where the Navier–Stokes equations are usually solved, and ℬ⊂ℝd representing the region occupied by the structure.[7] The fluid is described in an Eulerian framework while the solid is described in a Lagrangian formulation. This means that the solid domain ℬ can be expressed as the image of a fixed reference configuration ℬ* through the Lagrangian map x=X(s,t).[7]

The solid structure interacts with the fluid through its interface Γ=∂ℬ, and from the mathematical point of view, it is required that the following kinematic condition is satisfied at the interface Γ:u(x,t)=∂X(s,t)∂t|x=X(s,t)where u is the fluid velocity and x=X(s,t) is the position of the Lagrangian coordinate s at time t. The latter kinematic condition can be seen as a specific Dirichlet boundary condition at the interface Γ which involves both the unknowns u and X of the problem.[7]

Therefore, in the spirit of the fictitious domain method, we can extend the fluid equation to the whole fixed computational domain Ω=Ωf∪ℬ and impose weakly the kinematic coupling condition through the Lagrange multiplier λ, requiring that the fictitious velocity u matches the time derivative of the Lagrangian map X in a weak sense inside the solid domain ℬ.[7]

Numerical simulation of a heart valve

This requirement leads to the following weak formulation:[7]{af(u,v)−b(p,v)+⟨λ,v(X)⟩ℬ*=Ff(v)∀v∈H01(Ω)b(q,u)=0∀q∈L02(Ω)as(X,Y)−⟨λ,Y⟩ℬ*=Fs(Y)∀Y∈H1(ℬ*)⟨u(X)−∂X∂t,η⟩ℬ*=0∀η∈Λ

where af(⋅,⋅) and b(⋅,⋅) are the classical bilinear forms coming from the Navier–Stokes equations, while as(⋅,⋅) is the bilinear form associated with the elastic model describing the solid.[7]

The kinematic condition is represented in a weak sense by the last equation, and the Lagrange multiplier λ acts as a body force on both the fluid and solid equations in order to fulfil this condition.[8]

## Numerical aspects

The fictitious domain method proves to be effective for the numerical simulation of moving obstacles, as it employs two independent meshes: a fixed background grid and an independent, overlapping mesh that describes the obstacle geometry.[5] As we have seen in the previous applications, in order to weakly enforce the interface condition, we have to compute the duality pairing ⟨λ,v⟩Ωob over the obstacle domain Ωob.[7][8][9]

The Lagrange multiplier λ is defined over Ωob while the test function 𝐯 is defined over the whole computational domain Ω. Therefore, upon finite element method discretization, the discrete Lagrange multiplier λh is defined on the overlapping mesh, whereas the discrete test function 𝐯h lives on the background mesh.[10] As a result, computing the above coupling term requires the evaluation of basis functions defined on the background mesh at locations that belong to the overlapping mesh, which may be change over time. This task is computationally non-trivial and may introduce numerical errors and possibly wrong solutions.[10]

This issue can be handled with different strategies:[11]

* Augmented quadrature rule: sub-triangulate the overlapping region in order to adopt an enhanced quadrature formula capable of accurately computing the coupling term. While exact, this approach can be computationally expensive.[11]
* Inexact quadrature rule: maintain the original quadrature used for the overlapping mesh, accepting a quadrature error that is controlled by the mesh size h.[10]
* Collocation method: project λh over the discrete space of 𝐯h in order to have both the functions defined on the same mesh.[8]

## See also

* Computational fluid dynamics
* Finite element method
* Fluid structure interaction
* Immersed boundary method
* Level set method
* Meshfree method

## References

1. ↑ 1.0 1.1 Glowinski, Roland; Pan, Tsorng-Whay; Periaux, Jacques (1994-01-01). "A fictitious domain method for Dirichlet problem and applications". Computer Methods in Applied Mechanics and Engineering 111 (3): 283–303. doi:10.1016/0045-7825(94)90135-X. ISSN 0045-7825. Bibcode: 1994CMAME.111..283G. https://dx.doi.org/10.1016/0045-7825(94)90135-X.
2. ↑ Bordas, Stéphane P. A.; Burman, Erik; Larson, Mats G. et al., eds (2017) (in en). Geometrically Unfitted Finite Element Methods and Applications. Lecture Notes in Computational Science and Engineering. 121. doi:10.1007/978-3-319-71431-8. ISBN 978-3-319-71430-1. https://link.springer.com/book/10.1007/978-3-319-71431-8.
3. ↑ Mittal, Rajat; Iaccarino, Gianluca (2005-01-21). "Immersed Boundary Methods" (in en). Annual Review of Fluid Mechanics 37 (1): 239–261. doi:10.1146/annurev.fluid.37.061903.175743. ISSN 0066-4189. Bibcode: 2005AnRFM..37..239M. https://www.annualreviews.org/content/journals/10.1146/annurev.fluid.37.061903.175743.
4. ↑ Gibou, Frederic; Fedkiw, Ronald; Osher, Stanley (2018-01-15). "A review of level-set methods and some recent applications". Journal of Computational Physics 353: 82–109. doi:10.1016/j.jcp.2017.10.006. ISSN 0021-9991. Bibcode: 2018JCoPh.353...82G. https://www.sciencedirect.com/science/article/pii/S0021999117307441.
5. ↑ 5.0 5.1 Boffi, Daniele; Cavallini, Nicola; Gastaldi, Lucia (2015-03-04). "The Finite Element Immersed Boundary Method with Distributed Lagrange Multiplier" (in en). SIAM Journal on Numerical Analysis 53 (6): 2584–2604. doi:10.1137/140978399. ISSN 0036-1429. http://epubs.siam.org/doi/10.1137/140978399.
6. ↑ 6.0 6.1 Lu, Kaizhou; Augarde, Charles E.; Coombs, William M.; Hu, Zhendong (2019-05-01). "Weak impositions of Dirichlet boundary conditions in solid mechanics: A critique of current approaches and extension to partially prescribed boundaries". Computer Methods in Applied Mechanics and Engineering 348: 632–659. doi:10.1016/j.cma.2019.01.035. ISSN 0045-7825. Bibcode: 2019CMAME.348..632L. https://www.sciencedirect.com/science/article/pii/S004578251930060X.
7. ↑ 7.0 7.1 7.2 7.3 7.4 7.5 7.6 7.7 Boffi, Daniele; Gastaldi, Lucia (2017). "A fictitious domain approach with Lagrange multiplier for fluid-structure interactions". Numerische Mathematik 135 (3): 711–732. doi:10.1007/s00211-016-0814-1.
8. ↑ 8.0 8.1 8.2 8.3 8.4 Yu, Zhaosheng (2005-07-20). "A DLM/FD method for fluid/flexible-body interactions". Journal of Computational Physics 207 (1): 1–27. doi:10.1016/j.jcp.2004.12.026. ISSN 0021-9991. Bibcode: 2005JCoPh.207....1Y. https://www.sciencedirect.com/science/article/pii/S0021999105000148.
9. ↑ 9.0 9.1 Bertrand, F.; Tanguy, P. A.; Thibault, F. (1997). "A three-dimensional fictitious domain method for incompressible fluid flow problems" (in en). International Journal for Numerical Methods in Fluids 25 (6): 719–736. doi:10.1002/(SICI)1097-0363(19970930)25:6<719::AID-FLD585>3.0.CO;2-K. ISSN 1097-0363. Bibcode: 1997IJNMF..25..719B. https://onlinelibrary.wiley.com/doi/abs/10.1002/(SICI)1097-0363(19970930)25:6<719::AID-FLD585>3.0.CO;2-K.
10. ↑ 10.0 10.1 10.2 Boffi, Daniele; Credali, Fabio; Gastaldi, Lucia (2024). "Quadrature error estimates on non-matching grids in a fictitious domain framework for fluid-structure interaction problems". arXiv:2406.03981 [math.NA].
11. ↑ 11.0 11.1 Boffi, Daniele; Credali, Fabio; Gastaldi, Lucia (2022). "On the interface matrix for fluid–structure interaction problems with fictitious domain approach". Computer Methods in Applied Mechanics and Engineering 401. doi:10.1016/j.cma.2022.115650.

## External links

* vanDANA – FSI solver based on the fictitious domain method
* Daniele Boffi's home page
* Lucia Gastaldi's home page

* v
* t
* e

Numerical partial differential equations by method

Finite difference|

| Parabolic|

* Forward-time central-space (FTCS)
* Crank–Nicolson

|

Hyperbolic|

* Lax–Friedrichs
* Lax–Wendroff
* MacCormack
* Upwind
* Method of characteristics

Others|

* Alternating direction-implicit (ADI)
* Finite-difference time-domain (FDTD)

Finite volume|

* Godunov
* High-resolution
* Monotonic upstream-centered (MUSCL)
* Advection upstream-splitting (AUSM)
* Riemann solver
* essentially non-oscillatory (ENO)
* weighted essentially non-oscillatory (WENO)

Finite element|

* hp-FEM
* Extended (XFEM)
* Discontinuous Galerkin (DG)
* Spectral element (SEM)
* Mortar
* Gradient discretisation (GDM)
* Loubignac iteration
* Smoothed (S-FEM)

Meshless/Meshfree|

* Smoothed-particle hydrodynamics (SPH)
* Moving particle semi-implicit method (MPS)
* Material point method (MPM)
* Particle-in-cell (PIC)

Domain decomposition|

* Schur complement
* Fictitious domain
* Schwarz alternating
* additive
* abstract additive
* Neumann–Dirichlet
* Neumann–Neumann
* Poincaré–Steklov operator
* Balancing (BDD)
* Balancing by constraints (BDDC)
* Tearing and interconnect (FETI)
* FETI-DP

Others|

* Spectral
* Pseudospectral (DVR)
* Method of lines
* Multigrid
* Collocation
* Level-set
* Boundary element
* Immersed boundary
* Analytic element
* Particle-in-cell
* Isogeometric analysis
* Infinite difference method
* Infinite element method
* Galerkin method
* Petrov–Galerkin method
* Validated numerics
* Computer-assisted proof

0.00

(0 votes)

Original source: https://en.wikipedia.org/wiki/Fictitious domain method. Read more

|

Retrieved from "https://handwiki.org/wiki/index.php?title=Fictitious_domain_method&oldid=3983128"
*[v]: View this template
*[t]: Discuss this template
*[e]: Edit this template