MOSQUITO: An efficient finite difference scheme for numerical simulation of 2D advection

An explicit finite difference method for the treatment of the advective terms in the 2D equation of unsteady scalar transport is presented. The scheme is a conditionally stable extension to two dimensions of the popular QUICKEST scheme. It is deduced imposing the vanishing of selected components of the truncation error for the case of steady uniform flow. The method is then extended to solve the conservative form of the depth‐averaged transport equation. Details of the accuracy and stability analysis of the numerical scheme with test case results are given, together with a comparison with other existing schemes suitable for the long‐term computations needed in environmental modelling. Although with a truncation error of formal order 0(ΔxΔt, ΔyΔt, Δt²), the present scheme is shown actually to be of an accuracy comparable with schemes of third‐order in space, while requiring a smaller computational effort and/or having better stability properties. In principle, the method can be easily extended to the 3D case. Copyright © 1999 John Wiley & Sons, Ltd.     

High‐order compact finite difference schemes for the vorticity–divergence representation of the spherical shallow water equations

This work is devoted to the application of the super compact finite difference method (SCFDM) and the combined compact finite difference method (CCFDM) for spatial differencing of the spherical shallow water equations in terms of vorticity, divergence, and height. The fourth‐order compact, the sixth‐order and eighth‐order SCFDM, and the sixth‐order and eighth‐order CCFDM schemes are used for the spatial differencing. To advance the solution in time, a semi‐implicit Runge–Kutta method is used. In addition, to control the nonlinear instability, an eighth‐order compact spatial filter is employed. For the numerical solution of the elliptic equations in the problem, a direct hybrid method, which consists of a high‐order compact scheme for spatial differencing in the latitude coordinate and a fast Fourier transform in longitude coordinate, is utilized. The accuracy and convergence rate for all methods are verified against exact analytical solutions. Qualitative and quantitative assessments of the results for an unstable barotropic mid‐latitude zonal jet employed as an initial condition are addressed. It is revealed that the sixth‐order and eighth‐order CCFDMs and SCFDMs lead to a remarkable improvement of the solution over the fourth‐order compact method. Copyright © 2015 John Wiley & Sons, Ltd.     

浅水方程组合型超紧致差分格式

提出一族组合型超紧致差分格式(CSCD),对CSCD的数值特性作了分析,并同其他中心型差分格式进行比较。从定性角度,得出同阶中心差分格式中,CSCD格式的截断误差系数最小的结论。从定量角度,利用Fou-rier分析方法分析了CSCD格式的分辨率,并同其他中心型差分格式比较,得出CSCD格式有较高的分辨率的结论。把10阶CSCD格式应用于KdV-Burgers方程和浅水方程的数值模拟,给出两个应用算例。数值实验表明CSCD格式不仅有理论上的高精度,而且有良好的稳定性和收敛性。     

Compact finite difference schemes on non‐uniform meshes. Application to direct numerical simulations of compressible flows

In this paper, the development of a fourth‐ (respectively third‐) order compact scheme for the approximation of first (respectively second) derivatives on non‐uniform meshes is studied. A full inclusion of metrics in the coefficients of the compact scheme is proposed, instead of methods using Jacobian transformation. In the second part, an analysis of the numerical scheme is presented. A numerical analysis of truncation errors, a Fourier analysis completed by stability calculations in terms of both semi‐ and fully discrete eigenvalue problems are presented. In those eigenvalue problems, the pure convection equation for the first derivative, and the pure diffusion equation for the second derivative are considered. The last part of this paper is dedicated to an application of the numerical method to the simulation of a compressible flow requiring variable mesh size: the direct numerical simulation of compressible turbulent channel flow. Present results are compared with both experimental and other numerical (DNS) data in the literature. The effects of compressibility and acoustic waves on the turbulent flow structure are discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

Symmetry techniques for the numerical solution of the 2D Euler equations at impermeable boundaries

The implementation of boundary conditions at rigid, fixed wall boundaries in inviscid Euler solutions by upwind, finite volume methods is considered. Some current methods are reviewed. Two new boundary condition procedures, denoted as the symmetry technique and the curvature-corrected symmetry technique are then presented. Their behaviour in relation to the problem of the subsonic flow about blunt and slender elliptic bodies is analysed. The subsonic flow inside the Stanitz elbow is then computed. The symmetry technique is proven to be as accurate as one of the current methods, second-order pressure extrapolation technique. Finally, for arbitrary curved geometries, dramatic advantages of the curvature-corrected symmetry technique over the other methods are shown. © 1998 John Wiley & Sons, Ltd.     

Numerical solution of the incompressible Navier–Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier–Stokes equations is presented. The scheme is constructed on a staggered‐mesh grid system. The convection terms are discretized with a fifth‐order‐accurate upwind compact difference approximation, the viscous terms are discretized with a sixth‐order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth‐order difference approximation on a cell‐centered mesh. Time advancement uses a three‐stage Runge–Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented. Copyright © 1999 John Wiley & Sons, Ltd.     

A spectral–finite difference solution of the Navier–Stokes equations in three dimensions

A new computational code for the numerical integration of the three-dimensional Navier–Stokes equations in their non-dimensional velocity–pressure formulation is presented. The system of non-linear partial differential equations governing the time-dependent flow of a viscous incompressible fluid in a channel is managed by means of a mixed spectral–finite difference method, in which different numerical techniques are applied: Fourier decomposition is used along the homogeneous directions, second-order Crank–Nicolson algorithms are employed for the spatial derivatives in the direction orthogonal to the solid walls and a fourth-order Runge–Kutta procedure is implemented for both the calculation of the convective term and the time advancement. The pressure problem, cast in the Helmholtz form, is solved with the use of a cyclic reduction procedure. No-slip boundary conditions are used at the walls of the channel and cyclic conditions are imposed at the other boundaries of the computing domain. Results are provided for different values of the Reynolds number at several time steps of integration and are compared with results obtained by other authors. © 1998 John Wiley & Sons, Ltd.     

Goal-oriented error estimation applied to direct solution of steady-state analysis with frequency-domain finite element method

Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems,a goal-oriented error estimation method with extended degrees of freedom is developed.It leads to the high quality local error bounds in the problem of the direct-solution steady-state dynamic analysis with a frequency-domain finite element,which involves the enrichments with plural variable basis functions.The solution of the steady-state dynamic procedure calculates the harmonic response directly in terms of the physical degrees of freedom in the model,which uses the mass,damping,and stiffness matrices of the system.A three-dimensional finite element example is carried out to illustrate the computational procedures.     

Analysis of super compact finite difference method and application to simulation of vortex–shock interaction

Turbulence and aeroacoustic noise high‐order accurate schemes are required, and preferred, for solving complex flow fields with multi‐scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth‐order traditional and compact finite difference approximation. The comparison shows that the sixth‐order accurate super compact method has higher resolving efficiency. The sixth‐order super compact method, with a three‐stage Runge–Kutta method for approximation of the compressible Navier–Stokes equations, is used to solve the complex flow structures induced by vortex–shock interactions. The basic nature of the near‐field sound generated by interaction is studied. Copyright © 2001 John Wiley & Sons, Ltd.     

Analysis of hybrid algorithms for the Navier–Stokes equations with respect to hydrodynamic stability theory

Hybrid three‐dimensional algorithms for the numerical integration of the incompressible Navier–Stokes equations are analyzed with respect to hydrodynamic stability in both linear and nonlinear fields. The computational schemes are mixed—spectral and finite differences—and are applied to the case of the channel flow driven by constant pressure gradient; time marching is handled with the fractional step method. Different formulations—fully explicit convective term, partially and fully implicit viscous term combined with uniform, stretched, staggered and non‐staggered meshes, x‐velocity splitted and non‐splitted in average and perturbation component – are analyzed by monitoring the evolution in time of both small and finite amplitude perturbations of the mean flow. The results in the linear field are compared with correspondent solutions of the Orr–Sommerfeld equation; in the nonlinear field, the comparison is made with results obtained by other authors. Copyright © 2002 John Wiley & Sons, Ltd.