Adaptive time-stepping for incompressible flow part I: scalar advection-diffusion

Even the simplest advection-diffusion problems can exhibit multiple time scales. This means that robust variable step time integrators are a prerequisite if such problems are to be efficiently solved computationally. The performance of the second order trapezoid rule using an explicit Adams–Bashforth method for error control is assessed in this work. This combination is particularly well suited to long time integration of advection-dominated problems. Herein it is shown that a stabilized implementation of the trapezoid rule leads to a very effective integrator in other situations: specifically diffusion problems with rough initial data; and general advection-diffusion problems with different physical time scales governing the system evolution

A sequential regularization method for time-dependent incompressible Navier--Stokes equations

The objective of the paper is to present a method, called sequential regularization method (SRM), for the nonstationary incompressible Navier-Stokes equations from the viewpoint of regularization of differential-algebraic equations (DAEs) , and to provide a way to apply a DAE method to partial differential-algebraic equations (PDAEs). The SRM is a functional iterative procedure. It is proved that its convergence rate is O(ffl m ), where m is the number of the SRM iterations and ffl is the regularization parameter. The discretization and implementation issues of the method are considered. In particular, a simple explicit difference scheme is analyzed and its stability is proved under the usual step size condition of explicit schemes. It appears that the SRM formulation is new in the Navier-Stokes context. Unlike other regularizations or pseudo-compressibility methods in the Navier-Stokes context, the regularization parameter ffl in the SRM need not be very small, and the regularized..

On the calculation of normals in free-surface flow problems

The use of boundary-conforming finite element methods is considered for the solution of surface-tension-dominated free-surface flow problems in three dimensions. This class of method is based upon the use of a moving mesh whose velocity is driven by the motion of the free surface, which is in turn determined via a kinematic boundary condition for the normal velocity. The significance of the method used to compute the normal direction at the finite element node points for a C0 piecewise-polynomial free surface is investigated. In particular, it is demonstrated that the concept of mass-consistent normals on an isoparametric quadratic tetrahedral mesh is flawed. In this case an alternative, purely geometric, normal is shown to lead to a far more robust numerical algorithm

Numerical wave propagation for the triangular P1DGP1_{DG}-P2P2 finite element pair

Inertia-gravity mode and Rossby mode dispersion properties are examined for discretisations of the linearized rotating shallow-water equations using the $P1_{DG}$-$P2$ finite element pair on arbitrary triangulations in planar geometry. A discrete Helmholtz decomposition of the functions in the velocity space based on potentials taken from the pressure space is used to provide a complete description of the numerical wave propagation for the discretised equations. In the $f$-plane case, this decomposition is used to obtain decoupled equations for the geostrophic modes, the inertia-gravity modes, and the inertial oscillations. As has been noticed previously, the geostrophic modes are steady. The Helmholtz decomposition is used to show that the resulting inertia-gravity wave equation is third-order accurate in space. In general the \pdgp finite element pair is second-order accurate, so this leads to very accurate wave propagation. It is further shown that the only spurious modes supported by this discretisation are spurious inertial oscillations which have frequency $f$, and which do not propagate. The Helmholtz decomposition also allows a simple derivation of the quasi-geostrophic limit of the discretised $P1_{DG}$-$P2$ equations in the $\beta$-plane case, resulting in a Rossby wave equation which is also third-order accurate.Comment: Revised version prior to final journal submissio

Comparison of Wide and Compact Fourth Order Formulations of the Navier-Stokes Equations

In this study the numerical performances of wide and compact fourth order formulation of the steady 2-D incompressible Navier-Stokes equations will be investigated and compared with each other. The benchmark driven cavity flow problem will be solved using both wide and compact fourth order formulations and the numerical performances of both formulations will be presented and also the advantages and disadvantages of both formulations will be discussed

On the stability of bubble functions and a stabilized mixed finite element formulation for the Stokes problem

In this paper we investigate the relationship between stabilized and enriched finite element formulations for the Stokes problem. We also present a new stabilized mixed formulation for which the stability parameter is derived purely by the method of weighted residuals. This new formulation allows equal order interpolation for the velocity and pressure fields. Finally, we show by counterexample that a direct equivalence between subgrid-based stabilized finite element methods and Galerkin methods enriched by bubble functions cannot be constructed for quadrilateral and hexahedral elements using standard bubble functions.Comment: 25 pages, 13 figures (The previous version was compiled by mistake with the wrong style file, the current one uses amsart, and there is no difference in the text or the figures

Improving convergence in smoothed particle hydrodynamics simulations without pairing instability

The numerical convergence of smoothed particle hydrodynamics (SPH) can be severely restricted by random force errors induced by particle disorder, especially in shear flows, which are ubiquitous in astrophysics. The increase in the number NH of neighbours when switching to more extended smoothing kernels at fixed resolution (using an appropriate definition for the SPH resolution scale) is insufficient to combat these errors. Consequently, trading resolution for better convergence is necessary, but for traditional smoothing kernels this option is limited by the pairing (or clumping) instability. Therefore, we investigate the suitability of the Wendland functions as smoothing kernels and compare them with the traditional B-splines. Linear stability analysis in three dimensions and test simulations demonstrate that the Wendland kernels avoid the pairing instability for all NH, despite having vanishing derivative at the origin (disproving traditional ideas about the origin of this instability; instead, we uncover a relation with the kernel Fourier transform and give an explanation in terms of the SPH density estimator). The Wendland kernels are computationally more convenient than the higher-order B-splines, allowing large NH and hence better numerical convergence (note that computational costs rise sub-linear with NH). Our analysis also shows that at low NH the quartic spline kernel with NH ~= 60 obtains much better convergence then the standard cubic spline.Comment: substantially revised version, accepted for publication in MNRAS, 15 pages, 13 figure

SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch

We present a novel implementation of Smoothed Particle Hydrodynamics (SPHS) that uses the spatial derivative of the velocity divergence as a higher order dissipation switch. Our switch -- which is second order accurate -- detects flow convergence before it occurs. If particle trajectories are going to cross, we switch on the usual SPH artificial viscosity, as well as conservative dissipation in all advected fluid quantities (for example, the entropy). The viscosity and dissipation terms (that are numerical errors) are designed to ensure that all fluid quantities remain single-valued as particles approach one another, to respect conservation laws, and to vanish on a given physical scale as the resolution is increased. SPHS alleviates a number of known problems with `classic' SPH, successfully resolving mixing, and recovering numerical convergence with increasing resolution. An additional key advantage is that -- treating the particle mass similarly to the entropy -- we are able to use multimass particles, giving significantly improved control over the refinement strategy. We present a wide range of code tests including the Sod shock tube, Sedov-Taylor blast wave, Kelvin-Helmholtz Instability, the `blob test', and some convergence tests. Our method performs well on all tests, giving good agreement with analytic expectations.Comment: 21 pages; 15 Figures. Submitted to MNRAS. Comments welcom

Dynamically Correct Formulations of the Linearised Navier-Stokes Equations

Motivated by the need to efficiently obtain low-order models of fluid flows around complex geometries for the purpose of feedback control system design, this paper considers the effect on system dynamics of basing plant models on different formulations of the linearised Navier-Stokes equations. We consider the dynamics of a single computational node formed by spatial discretisation of the governing equations in both primitive variables (momentum equation & continuity equation) and pressure Poisson equation (PPE) formulations. This reveals fundamental numerical differences at the nodal level, whose effects on the system dynamics at the full system level are exemplified by considering the corresponding formulations of a two-dimensional (2D) channel flow, subjected to a variety of different boundary conditions

Singular and regular solutions of a non-linear parabolic system

We study a dissipative nonlinear equation modelling certain features of the Navier-Stokes equations. We prove that the evolution of radially symmetric compactly supported initial data does not lead to singularities in dimensions $n\leq 4$. For dimensions $n>4$ we present strong numerical evidence supporting existence of blow-up solutions. Moreover, using the same techniques we numerically confirm a conjecture of Lepin regarding existence of self-similar singular solutions to a semi-linear heat equation.Comment: 16 page