Top Qs
Timeline
Chat
Perspective
Fictitious domain method
Numerical analysis method From Wikipedia, the free encyclopedia
Remove ads
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]


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]
Remove ads
Weak imposition of Dirichlet boundary conditions
Let 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 on the boundary of an immersed circle described by the interface . In this context, the fictitious domain method allows us to solve the Stokes equation over the whole domain , and to impose weakly the constraint by means of the Lagrange multiplier .

The fictitious domain weak formulation reads:[9]
The last equation of the system guarantees that, in a weak sense, the interface condition is satisfied over the whole domain , and in particular on the interface , while in the equation is left untouched.[6]
The Lagrange multiplier acts on the fluid as an external force defined only in that enforces the fictitious fluid velocity to fulfil the constraint inside the circle.[8]
Remove ads
Fluid-structure interaction application
As it is typically done in fluid-structure interaction problems, we consider a fluid domain , where the Navier–Stokes equations are usually solved, and 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 .[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 :where is the fluid velocity and is the position of the Lagrangian coordinate at time . The latter kinematic condition can be seen as a specific Dirichlet boundary condition at the interface which involves both the unknowns and 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 and impose weakly the kinematic coupling condition through the Lagrange multiplier , requiring that the fictitious velocity matches the time derivative of the Lagrangian map in a weak sense inside the solid domain .[7]

This requirement leads to the following weak formulation:[7]
where and are the classical bilinear forms coming from the Navier–Stokes equations, while 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]
Remove ads
Numerical aspects
Summarize
Perspective
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 over the obstacle domain .[7][8][9]
The Lagrange multiplier is defined over while the test function is defined over the whole computational domain . Therefore, upon finite element method discretization, the discrete Lagrange multiplier is defined on the overlapping mesh, whereas the discrete test function 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 .[10]
- Collocation method: project over the discrete space of in order to have both the functions defined on the same mesh.[8]
Remove ads
See also
References
External links
Wikiwand - on
Seamless Wikipedia browsing. On steroids.
Remove ads