Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA fictitious domain approach for the Stokes problem based on the extended finite element method
146
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In the present work, we propose to extend to the Stokes problem a fictitious domain approach inspired by eXtended Finite Element Method and studied for Poisson problem in [Renard]. The method allows computations in domains whose boundaries do not match. A mixed finite element method is used for fluid flow. The interface between the fluid and the structure is localized by a level-set function. Dirichlet boundary conditions are taken into account using Lagrange multiplier. A stabilization term is introduced to improve the approximation of the normal trace of the Cauchy stress tensor at the interface and avoid the inf-sup condition between the spaces for velocity and the Lagrange multiplier. Convergence analysis is given and several numerical tests are performed to illustrate the capabilities of the method.
rate research
Read More
In this work we develop a fictitious domain method for the Stokes problem which allows computations in domains whose boundaries do not depend on the mesh. The method is based on the ideas of Xfem and has been first introduced for the Poisson problem. The fluid part is treated by a mixed finite element method, and a Dirichlet condition is imposed by a Lagrange multiplier on an immersed structure localized by a level-set function. A stabilization technique is carried out in order to get the convergence for this multiplier. The latter represents the forces that the fluid applies on the structure. The aim is to perform fluid-structure simulations for which these forces have a central role. We illustrate the capacities of the method by extending it to the incompressible Navier-Stokes equations coupled with a moving rigid solid.
In this article we study a mixed finite element formulation for solving the Stokes problem with general surface forces that induce a jump of the normal trace of the stress tensor, on an interface that splits the domain into two subdomains. Equality of velocities is assumed at the interface. The interface conditions are taken into account with multipliers. A suitable Lagrangian functional is introduced for deriving a weak formulation of the problem. A specificity of this work is the consideration of the interface with a fictitious domain approach. The latter is inspired by the XFEM approach in the sense that cut-off functions are used, but it is simpler to implement since no enrichment by singular functions is provided. In that context, getting convergence for the dual variables defined on the interface is non-trivial. For that purpose, an augmented Lagrangian technique stabilizes the convergence of the multipliers, which is important because their value would determine the dynamics of the interface in an unsteady framework. Theoretical analysis is provided, where we show that a discrete inf-sup condition, independent of the mesh size, is satisfied for the stabilized formulation. This guarantees optimal convergence rates, that we observe with numerical tests. The capacity of the method is demonstrated with robustness tests, and with an unsteady model tested for deformations of the interface that correspond to ellipsoidal shapes in dimension 2.
The paper extends a stabilized fictitious domain finite element method initially developed for the Stokes problem to the incompressible Navier-Stokes equations coupled with a moving solid. This method presents the advantage to predict an optimal approximation of the normal stress tensor at the interface. The dynamics of the solid is governed by the Newtons laws and the interface between the fluid and the structure is materialized by a level-set which cuts the elements of the mesh. An algorithm is proposed in order to treat the time evolution of the geometry and numerical results are presented on a classical benchmark of the motion of a disk falling in a channel.
This paper constructs and analyzes a boundary correction finite element method for the Stokes problem based on the Scott-Vogelius pair on Clough-Tocher splits. The velocity space consists of continuous piecewise quadratic polynomials, and the pressure space consists of piecewise linear polynomials without continuity constraints. A Lagrange multiplier space that consists of continuous piecewise quadratic polynomials with respect to boundary partition is introduced to enforce boundary conditions as well as to mitigate the lack of pressure-robustness. We prove several inf-sup conditions, leading to the well-posedness of the method. In addition, we show that the method converges with optimal order and the velocity approximation is divergence free.
In this paper, a stabilized extended finite element method is proposed for Stokes interface problems on unfitted triangulation elements which do not require the interface align with the triangulation. The velocity solution and pressure solution on each side of the interface are separately expanded in the standard nonconforming piecewise linear polynomials and the piecewise constant polynomials, respectively. Harmonic weighted fluxes and arithmetic fluxes are used across the interface and cut edges (segment of the edges cut by the interface), respectively. Extra stabilization terms involving velocity and pressure are added to ensure the stable inf-sup condition. We show a priori error estimates under additional regularity hypothesis. Moreover, the errors {in energy and $L^2$ norms for velocity and the error in $L^2$ norm for pressure} are robust with respect to the viscosity {and independent of the location of the interface}. Results of numerical experiments are presented to {support} the theoretical analysis.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
