8,956 research outputs found
Hydrodynamic Flows on Curved Surfaces: Spectral Numerical Methods for Radial Manifold Shapes
We formulate hydrodynamic equations and spectrally accurate numerical methods
for investigating the role of geometry in flows within two-dimensional fluid
interfaces. To achieve numerical approximations having high precision and level
of symmetry for radial manifold shapes, we develop spectral Galerkin methods
based on hyperinterpolation with Lebedev quadratures for -projection to
spherical harmonics. We demonstrate our methods by investigating hydrodynamic
responses as the surface geometry is varied. Relative to the case of a sphere,
we find significant changes can occur in the observed hydrodynamic flow
responses as exhibited by quantitative and topological transitions in the
structure of the flow. We present numerical results based on the
Rayleigh-Dissipation principle to gain further insights into these flow
responses. We investigate the roles played by the geometry especially
concerning the positive and negative Gaussian curvature of the interface. We
provide general approaches for taking geometric effects into account for
investigations of hydrodynamic phenomena within curved fluid interfaces.Comment: 14 figure
High-Order Unstructured Lagrangian One-Step WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems: Applications to Compressible Multi-Phase Flows
In this article we present the first better than second order accurate
unstructured Lagrangian-type one-step WENO finite volume scheme for the
solution of hyperbolic partial differential equations with non-conservative
products. The method achieves high order of accuracy in space together with
essentially non-oscillatory behavior using a nonlinear WENO reconstruction
operator on unstructured triangular meshes. High order accuracy in time is
obtained via a local Lagrangian space-time Galerkin predictor method that
evolves the spatial reconstruction polynomials in time within each element. The
final one-step finite volume scheme is derived by integration over a moving
space-time control volume, where the non-conservative products are treated by a
path-conservative approach that defines the jump terms on the element
boundaries. The entire method is formulated as an Arbitrary-Lagrangian-Eulerian
(ALE) method, where the mesh velocity can be chosen independently of the fluid
velocity.
The new scheme is applied to the full seven-equation Baer-Nunziato model of
compressible multi-phase flows in two space dimensions. The use of a Lagrangian
approach allows an excellent resolution of the solid contact and the resolution
of jumps in the volume fraction. The high order of accuracy of the scheme in
space and time is confirmed via a numerical convergence study. Finally, the
proposed method is also applied to a reduced version of the compressible
Baer-Nunziato model for the simulation of free surface water waves in moving
domains. In particular, the phenomenon of sloshing is studied in a moving water
tank and comparisons with experimental data are provided
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
Relaxation of a dewetting contact line Part 1: A full-scale hydrodynamic calculation
The relaxation of a dewetting contact line is investigated theoretically in
the so-called "Landau-Levich" geometry in which a vertical solid plate is
withdrawn from a bath of partially wetting liquid. The study is performed in
the framework of lubrication theory, in which the hydrodynamics is resolved at
all length scales (from molecular to macroscopic). We investigate the
bifurcation diagram for unperturbed contact lines, which turns out to be more
complex than expected from simplified 'quasi-static' theories based upon an
apparent contact angle. Linear stability analysis reveals that below the
critical capillary number of entrainment, Ca_c, the contact line is linearly
stable at all wavenumbers. Away from the critical point the dispersion relation
has an asymptotic behaviour sigma~|q| and compares well to a quasi-static
approach. Approaching Ca_c, however, a different mechanism takes over and the
dispersion evolves from |q| to the more common q^2. These findings imply that
contact lines can not be treated as universal objects governed by some
effective law for the macroscopic contact angle, but viscous effects have to be
treated explicitly.Comment: 21 pages, 9 figure
Spectral/hp element methods: recent developments, applications, and perspectives
The spectral/hp element method combines the geometric flexibility of the
classical h-type finite element technique with the desirable numerical
properties of spectral methods, employing high-degree piecewise polynomial
basis functions on coarse finite element-type meshes. The spatial approximation
is based upon orthogonal polynomials, such as Legendre or Chebychev
polynomials, modified to accommodate C0-continuous expansions. Computationally
and theoretically, by increasing the polynomial order p, high-precision
solutions and fast convergence can be obtained and, in particular, under
certain regularity assumptions an exponential reduction in approximation error
between numerical and exact solutions can be achieved. This method has now been
applied in many simulation studies of both fundamental and practical
engineering flows. This paper briefly describes the formulation of the
spectral/hp element method and provides an overview of its application to
computational fluid dynamics. In particular, it focuses on the use the
spectral/hp element method in transitional flows and ocean engineering.
Finally, some of the major challenges to be overcome in order to use the
spectral/hp element method in more complex science and engineering applications
are discussed
- …