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

A "poor man's" approach to topology optimization of natural convection problems

Topology optimization of natural convection problems is computationally expensive, due to the large number of degrees of freedom (DOFs) in the model and its two-way coupled nature. Herein, a method is presented to reduce the computational effort by use of a reduced-order model governed by simplified physics. The proposed method models the fluid flow using a potential flow model, which introduces an additional fluid property. This material property currently requires tuning of the model by comparison to numerical Navier-Stokes based solutions. Topology optimization based on the reduced-order model is shown to provide qualitatively similar designs, as those obtained using a full Navier-Stokes based model. The number of DOFs is reduced by 50% in two dimensions and the computational complexity is evaluated to be approximately 12.5% of the full model. We further compare to optimized designs obtained utilizing Newton's convection law.Comment: Preprint version. Please refer to final version in Structural Multidisciplinary Optimization https://doi.org/10.1007/s00158-019-02215-

Numerical integration of the contravariant integral form of the Navier–Stokes equations in time-dependent curvilinear coordinate systems for three-dimensional free surface flows

We propose a three-dimensional non-hydrostatic shock-capturing numerical model for the simulation of wave propagation, transformation and breaking, which is based on an original integral formulation of the contravariant Navier–Stokes equations, devoid of Christoffel symbols, in general time-dependent curvilinear coordinates. A coordinate transformation maps the time-varying irregular physical domain that reproduces the complex geometries of coastal regions to a fixed uniform computational one. The advancing of the solution is performed by a second-order accurate strong stability preserving Runge–Kutta fractional-step method in which, at every stage of the method, a predictor velocity field is obtained by the shock-capturing scheme and a corrector velocity field is added to the previous one, to produce a non-hydrostatic divergence-free velocity field and update the water depth. The corrector velocity field is obtained by numerically solving a Poisson equation, expressed in integral contravariant form, by a multigrid technique which uses a four-colour Zebra Gauss–Seidel line-by-line method as smoother. Several test cases are used to verify the dispersion and shock-capturing properties of the proposed model in time-dependent curvilinear grids

Brownian dynamics of rigid particles in an incompressible fluctuating fluid by a meshfree method

A meshfree Lagrangian method for the fluctuating hydrodynamic equations (FHEs) with fluid-structure interactions is presented. Brownian motion of the particle is investigated by direct numerical simulation of the fluctuating hydrodynamic equations. In this framework a bidirectional coupling has been introduced between the fluctuating fluid and the solid object. The force governing the motion of the solid object is solely due to the surrounding fluid particles. Since a meshfree formulation is used, the method can be extended to many real applications involving complex fluid flows. A three-dimensional implementation is presented. In particular, we observe the short and long-time behaviour of the velocity autocorrelation function (VACF) of Brownian particles and compare it with the analytical expression. Moreover, the Stokes-Einstein relation is reproduced to ensure the correct long-time behaviour of Brownian dynamics.Comment: 24 pages, 2 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