1,950 research outputs found
Strong and auxiliary forms of the semi-Lagrangian method for incompressible flows
We present a review of the semi-Lagrangian method for advection-diusion and incompressible Navier-Stokes equations discretized with high-order methods. In particular, we compare the strong form where the departure points are computed directly via backwards integration with the auxiliary form where an auxiliary advection equation is solved instead; the latter is also referred to as Operator Integration Factor Splitting (OIFS) scheme. For intermediate size of time steps the auxiliary form is preferrable but for large time steps only the strong form is stable
Simulating water-entry/exit problems using Eulerian-Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library
In this paper we employ two implementations of the fictitious domain (FD)
method to simulate water-entry and water-exit problems and demonstrate their
ability to simulate practical marine engineering problems. In FD methods, the
fluid momentum equation is extended within the solid domain using an additional
body force that constrains the structure velocity to be that of a rigid body.
Using this formulation, a single set of equations is solved over the entire
computational domain. The constraint force is calculated in two distinct ways:
one using an Eulerian-Lagrangian framework of the immersed boundary (IB) method
and another using a fully-Eulerian approach of the Brinkman penalization (BP)
method. Both FSI strategies use the same multiphase flow algorithm that solves
the discrete incompressible Navier-Stokes system in conservative form. A
consistent transport scheme is employed to advect mass and momentum in the
domain, which ensures numerical stability of high density ratio multiphase
flows involved in practical marine engineering applications. Example cases of a
free falling wedge (straight and inclined) and cylinder are simulated, and the
numerical results are compared against benchmark cases in literature.Comment: The current paper builds on arXiv:1901.07892 and re-explains some
parts of it for the reader's convenienc
A conservative, optimization-based semi-lagrangian spectral element method for passive tracer transport
We present a new optimization-based, conservative, and quasi-monotone
method for passive tracer transport. The scheme combines high-order spectral element
discretization in space with semi-Lagrangian time stepping. Solution of a singly linearly
constrained quadratic program with simple bounds enforces conservation and physically
motivated solution bounds. The scheme can handle efficiently a large number of passive
tracers because the semi-Lagrangian time stepping only needs to evolve the grid
points where the primitive variables are stored and allows for larger time steps than a
conventional explicit spectral element method. Numerical examples show that the use
of optimization to enforce physical properties does not affect significantly the spectral
accuracy for smooth solutions. Performance studies reveal the benefits of high-order approximations,
including for discontinuous solutions
Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes
We present a new family of high order accurate fully discrete one-step
Discontinuous Galerkin (DG) finite element schemes on moving unstructured
meshes for the solution of nonlinear hyperbolic PDE in multiple space
dimensions, which may also include parabolic terms in order to model
dissipative transport processes. High order piecewise polynomials are adopted
to represent the discrete solution at each time level and within each spatial
control volume of the computational grid, while high order of accuracy in time
is achieved by the ADER approach. In our algorithm the spatial mesh
configuration can be defined in two different ways: either by an isoparametric
approach that generates curved control volumes, or by a piecewise linear
decomposition of each spatial control volume into simplex sub-elements. Our
numerical method belongs to the category of direct
Arbitrary-Lagrangian-Eulerian (ALE) schemes, where a space-time conservation
formulation of the governing PDE system is considered and which already takes
into account the new grid geometry directly during the computation of the
numerical fluxes. Our new Lagrangian-type DG scheme adopts the novel a
posteriori sub-cell finite volume limiter method, in which the validity of the
candidate solution produced in each cell by an unlimited ADER-DG scheme is
verified against a set of physical and numerical detection criteria. Those
cells which do not satisfy all of the above criteria are flagged as troubled
cells and are recomputed with a second order TVD finite volume scheme. The
numerical convergence rates of the new ALE ADER-DG schemes are studied up to
fourth order in space and time and several test problems are simulated.
Finally, an application inspired by Inertial Confinement Fusion (ICF) type
flows is considered by solving the Euler equations and the PDE of viscous and
resistive magnetohydrodynamics (VRMHD).Comment: 39 pages, 21 figure
- …