20,990 research outputs found
Efficient and robust monolithic finite element multilevel Krylov subspace solvers for the solution of stationary incompressible Navier-Stokes equations
Multigrid methods belong to the best-known methods for solving linear systems
arising from the discretization of elliptic partial differential equations. The
main attraction of multigrid methods is that they have an asymptotically meshindependent
convergence behavior. Multigrid with Vanka (or local multilevel
pressure Schur complement method) as smoother have been frequently used for
the construction of very effcient coupled monolithic solvers for the solution of
the stationary incompressible Navier-Stokes equations in 2D and 3D. However,
due to its innate Gauß-Seidel/Jacobi character, Vanka has a strong influence
of the underlying mesh, and therefore, coupled multigrid solvers with Vanka
smoothing very frequently face convergence issues on meshes with high aspect
ratios. Moreover, even on very nice regular grids, these solvers may fail when
the anisotropies are introduced from the differential operator.
In this thesis, we develop a new class of robust and efficient monolithic finite
element multilevel Krylov subspace methods (MLKM) for the solution of the
stationary incompressible Navier-Stokes equations as an alternative to the coupled
multigrid-based solvers. Different from multigrid, the MLKM utilizes a
Krylov method as the basis in the error reduction process. The solver is based
on the multilevel projection-based method of Erlangga and Nabben, which accelerates
the convergence of the Krylov subspace methods by shifting the small
eigenvalues of the system matrix, responsible for the slow convergence of the
Krylov iteration, to the largest eigenvalue.
Before embarking on the Navier-Stokes equations, we first test our implementation
of the MLKM solver by solving scalar model problems, namely the
convection-diffusion problem and the anisotropic diffusion problem. We validate
the method by solving several standard benchmark problems. Next, we
present the numerical results for the solution of the incompressible Navier-Stokes
equations in two dimensions. The results show that the MLKM solvers produce
asymptotically mesh-size independent, as well as Reynolds number independent
convergence rates, for a moderate range of Reynolds numbers. Moreover, numerical
simulations also show that the coupled MLKM solvers can handle (both
mesh and operator based) anisotropies better than the coupled multigrid solvers
A Two-Level Finite Element Discretization of the Streamfunction Formulation of the Stationary Quasi-Geostrophic Equations of the Ocean
In this paper we proposed a two-level finite element discretization of the
nonlinear stationary quasi-geostrophic equations, which model the wind driven
large scale ocean circulation. Optimal error estimates for the two-level finite
element discretization were derived. Numerical experiments for the two-level
algorithm with the Argyris finite element were also carried out. The numerical
results verified the theoretical error estimates and showed that, for the
appropriate scaling between the coarse and fine mesh sizes, the two-level
algorithm significantly decreases the computational time of the standard
one-level algorithm.Comment: Computers and Mathematics with Applications 66 201
Finite element formulation of general boundary conditions for incompressible flows
We study the finite element formulation of general boundary conditions for
incompressible flow problems. Distinguishing between the contributions from the
inviscid and viscid parts of the equations, we use Nitsche's method to develop
a discrete weighted weak formulation valid for all values of the viscosity
parameter, including the limit case of the Euler equations. In order to control
the discrete kinetic energy, additional consistent terms are introduced. We
treat the limit case as a (degenerate) system of hyperbolic equations, using a
balanced spectral decomposition of the flux Jacobian matrix, in analogy with
compressible flows. Then, following the theory of Friedrich's systems, the
natural characteristic boundary condition is generalized to the considered
physical boundary conditions. Several numerical experiments, including standard
benchmarks for viscous flows as well as inviscid flows are presented
A time dependent Stokes interface problem: well-posedness and space-time finite element discretization
In this paper a time dependent Stokes problem that is motivated by a standard
sharp interface model for the fluid dynamics of two-phase flows is studied.
This Stokes interface problem has discontinuous density and viscosity
coefficients and a pressure solution that is discontinuous across an evolving
interface. This strongly simplified two-phase Stokes equation is considered to
be a good model problem for the development and analysis of finite element
discretization methods for two-phase flow problems. In view of the unfitted
finite element methods that are often used for two-phase flow simulations, we
are particularly interested in a well-posed variational formulation of this
Stokes interface problem in a Euclidean setting. Such well-posed weak
formulations, which are not known in the literature, are the main results of
this paper. Different variants are considered, namely one with suitable spaces
of divergence free functions, a discrete-in-time version of it, and variants in
which the divergence free constraint in the solution space is treated by a
pressure Lagrange multiplier. The discrete-in-time variational formulation
involving the pressure variable for the divergence free constraint is a natural
starting point for a space-time finite element discretization. Such a method is
introduced and results of numerical experiments with this method are presented
Numerical approximation of phase field based shape and topology optimization for fluids
We consider the problem of finding optimal shapes of fluid domains. The fluid
obeys the Navier--Stokes equations. Inside a holdall container we use a phase
field approach using diffuse interfaces to describe the domain of free flow. We
formulate a corresponding optimization problem where flow outside the fluid
domain is penalized. The resulting formulation of the shape optimization
problem is shown to be well-posed, hence there exists a minimizer, and first
order optimality conditions are derived.
For the numerical realization we introduce a mass conserving gradient flow
and obtain a Cahn--Hilliard type system, which is integrated numerically using
the finite element method. An adaptive concept using reliable, residual based
error estimation is exploited for the resolution of the spatial mesh.
The overall concept is numerically investigated and comparison values are
provided
- …