8,956 research outputs found

    Hydrodynamic Flows on Curved Surfaces: Spectral Numerical Methods for Radial Manifold Shapes

    Full text link
    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 L2L^2-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

    Full text link
    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

    Get PDF
    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

    Full text link
    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

    Get PDF
    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
    • …
    corecore