1,507 research outputs found
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
Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods
In this paper, we develop a new tensor-product based preconditioner for
discontinuous Galerkin methods with polynomial degrees higher than those
typically employed. This preconditioner uses an automatic, purely algebraic
method to approximate the exact block Jacobi preconditioner by Kronecker
products of several small, one-dimensional matrices. Traditional matrix-based
preconditioners require storage and
computational work, where is the degree of basis polynomials used, and
is the spatial dimension. Our SVD-based tensor-product preconditioner requires
storage, work in two spatial
dimensions, and work in three spatial dimensions.
Combined with a matrix-free Newton-Krylov solver, these preconditioners allow
for the solution of DG systems in linear time in per degree of freedom in
2D, and reduce the computational complexity from to
in 3D. Numerical results are shown in 2D and 3D for the
advection and Euler equations, using polynomials of degree up to . For
many test cases, the preconditioner results in similar iteration counts when
compared with the exact block Jacobi preconditioner, and performance is
significantly improved for high polynomial degrees .Comment: 40 pages, 15 figure
- …