Search for: All records

The paper discusses a reuse of matrix factorization as a building block in the Augmented Lagrangian (AL) and modified AL preconditioners for nonsymmetric saddle point linear algebraic systems. The strategy is applied to solve two‐dimensional incompressible fluid problems with efficiency rates independent of the Reynolds number. The solver is then tested to simulate motion of a surface fluid, an example of a two‐dimensional flow motivated by an interest in lateral fluidity of inextensible viscous membranes. Numerical examples include the Kelvin–Helmholtz instability problem posed on the sphere and on the torus. Some new eigenvalue estimates for the AL preconditioner are derived.

 

The paper studies the equilibrium configurations of inextensible elastic membranes exhibiting lateral fluidity. Using a continuum description of the membrane's motions based on the surface Navier–Stokes equations with bending forces, the paper derives differential equations governing the mechanical equilibrium. The equilibrium conditions are found to be independent of lateral viscosity and relate tension, pressure, and tangential velocity of the fluid. These conditions suggest that either the lateral fluid motion ceases or non-decaying stationary flow of mass can only be supported by surfaces with Killing vector fields, such as axisymmetric shapes. A shape equation is derived that extends the classical Helfrich model with an area constraint to membranes of non-negligible mass. Furthermore, the paper suggests a simple numerical method to compute solutions of the shape equation. Numerical experiments conducted reveal a diverse family of equilibrium configurations. The stability of equilibrium states involving lateral flow of mass remains an unresolved question. 

We construct and analyze a CutFEM discretization for the Stokes problem based on the Scott–Vogelius pair. The discrete piecewise polynomial spaces are defined on macro-element triangulations which are not fitted to the smooth physical domain. Boundary conditions are imposed via penalization through the help of a Nitsche-type discretization, whereas stability with respect to small and anisotropic cuts of the bulk elements is ensured by adding local ghost penalty stabilization terms. We show stability of the scheme as well as a divergence–free property of the discrete velocity outside an O ( h ) neighborhood of the boundary. To mitigate the error caused by the violation of the divergence–free condition, we introduce local grad–div stabilization. The error analysis shows that the grad–div parameter can scale like O ( h −1 ), allowing a rather heavy penalty for the violation of mass conservation, while still ensuring optimal order error estimates. 

The paper considers a system of equations that models a lateral flow of a Boussinesq–Scriven fluid on a passively evolving surface embedded in [Formula: see text]. For the resulting Navier–Stokes type system, posed on a smooth closed time-dependent surface, we introduce a weak formulation in terms of functional spaces on a space-time manifold defined by the surface evolution. The weak formulation is shown to be well-posed for any finite final time and without smallness conditions on data. We further extend an unfitted finite element method, known as TraceFEM, to compute solutions to the fluid system. Convergence of the method is demonstrated numerically. In another series of experiments we visualize lateral flows induced by smooth deformations of a material surface. 

Abstract This paper studies a model of two-phase flow with an immersed material viscous interface and a finite element method for the numerical solution of the resulting system of PDEs.The interaction between the bulk and surface media is characterized by no-penetration and slip with friction interface conditions.The system is shown to be dissipative, and a model stationary problem is proved to be well-posed.The finite element method applied in this paper belongs to a family of unfitted discretizations.The performance of the method when model and discretization parameters vary is assessed.Moreover, an iterative procedure based on the splitting of the system into bulk and surface problems is introduced and studied numerically.