Implementation of high‐order compact finite‐difference method to parabolized Navier–Stokes schemes

The numerical solution to the parabolized Navier–Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth‐order compact finite‐difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite‐difference algorithm of Beam and Warming type with a high‐order compact accuracy. A shock‐fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously computed with the fourth‐order accuracy. The computations are carried out for a benchmark case: hypersonic axisymmetric flow over a blunt cone at Mach 8. A sensitivity study is performed for the basic flowfield, including profiles and their derivatives obtained from the fourth‐order compact PNS and IPNS solutions, and the effects of grid size and numerical dissipation term used are discussed. The present results for the flowfield variables and also their derivatives are compared with those of other basic flow models to demonstrate the accuracy and efficiency of the proposed method. The present work represents the first known application of a high‐order compact finite‐difference method to the PNS schemes, which are computationally more efficient than Navier–Stokes solutions. Copyright © 2008 John Wiley & Sons, Ltd.     

Finite‐difference 4th‐order compact scheme for the direct numerical simulation of instabilities of shear layers

An algorithm based on the 4th‐order finite‐difference compact scheme is developed and applied in the direct numerical simulations of instabilities of channel flow. The algorithm is illustrated in the context of stream function formulation that leads to field equation involving 4th‐order spatial derivatives. The finite‐difference discretization in the wall‐normal direction uses five arbitrarily spaced points. The discretization coefficients are determined numerically, providing a large degree of flexibility for grid selection. The Fourier expansions are used in the streamwise direction. A hybrid Runge–Kutta/Crank–Nicholson low‐storage scheme is applied for the time discretization. Accuracy tests demonstrate that the algorithm does deliver the 4th‐order accuracy. The algorithm has been used to simulate the natural instability processes in channel flow as well as processes occurring when the flow is spatially modulated using wall transpiration. Extensions to three‐dimensional situations are suggested. Copyright © 2005 John Wiley & Sons, Ltd.     

Velocity–pressure coupling in finite difference formulations for the Navier–Stokes equations

A new numerical procedure for solving the two‐dimensional, steady, incompressible, viscous flow equations on a staggered Cartesian grid is presented in this paper. The proposed methodology is finite difference based, but essentially takes advantage of the best features of two well‐established numerical formulations, the finite difference and finite volume methods. Some weaknesses of the finite difference approach are removed by exploiting the strengths of the finite volume method. In particular, the issue of velocity–pressure coupling is dealt with in the proposed finite difference formulation by developing a pressure correction equation using the SIMPLE approach commonly used in finite volume formulations. However, since this is purely a finite difference formulation, numerical approximation of fluxes is not required. Results presented in this paper are based on first‐ and second‐order upwind schemes for the convective terms. This new formulation is validated against experimental and other numerical data for well‐known benchmark problems, namely developing laminar flow in a straight duct, flow over a backward‐facing step, and lid‐driven cavity flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Verification of a mixed high‐order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high‐order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier–Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity–velocity formulation for a two‐dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high‐order compact and non‐compact finite‐differences from fourth‐order to sixth‐order of accuracy. The other scheme used spectral methods instead of finite‐difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed‐order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright © 2009 John Wiley & Sons, Ltd.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

The advantages of using high‐order finite differences schemes in laminar–turbulent transition studies

This paper presents various finite difference schemes and compare their ability to simulate instability waves in a given flow field. The governing equations for two‐dimensional, incompressible flows were solved in vorticity–velocity formulation. Four different space discretization schemes were tested, namely, a second‐order central differences, a fourth‐order central differences, a fourth‐order compact scheme and a sixth‐order compact scheme. A classic fourth‐order Runge–Kutta scheme was used in time. The influence of grid refinement in the streamwise and wall normal directions were evaluated. The results were compared with linear stability theory for the evolution of small‐amplitude Tollmien–Schlichting waves in a plane Poiseuille flow. Both the amplification rate and the wavenumber were considered as verification parameters, showing the degree of dissipation and dispersion introduced by the different numerical schemes. The results confirmed that high‐order schemes are necessary for studying hydrodynamic instability problems by direct numerical simulation. Copyright © 2005 John Wiley & Sons, Ltd.     

The improved surface gradient method for flows simulation in variable bed topography channel using TVD‐MacCormack scheme

This paper reports four different approaches to discretize the source terms for the simulation of one‐dimensional open‐channel flows with rapidly varied bottom topography using TVD‐MacCormack scheme. Compared with other high‐resolution shock‐capturing schemes, MacCormack‐type predictor–corrector method is easy to implement and does not present any additional difficulty in dealing with the source terms. To avoid the generation of artificial numerical waves, if the bottom topography shows strong variation, special treatment of the source terms is still required to eliminate or reduce the artificial numerical error caused by adding TVD corrections to the method. The computed results demonstrated that the improved surface gradient method is more suitable for simulating open‐channel flow with highly irregular bed topography by using the surface gradient instead of the depth gradient for TVD corrections and considering the balancing of the source terms and the flux gradients. Copyright © 2003 John Wiley & Sons, Ltd.     

Large Eddy Simulation of Flow around a Rectangular Cylinder Using the Finite Element Method

In previous papers, we proposed finite element schemes based on the Petrov-Galerkin weak formulation using exponential weighting functions for solving accurately, and in a stable manner, the flow field of an incompressible viscous fluid. In this paper, we present the Petrov-Galerkin finite element scheme for turbulent flow fields based on large eddy simulation using the standard Smagorinsky model with the Van Driest damping function. The filtered incompressible Navier-Stokes equations are numerically integrated in time by using a fractional step strategy with second-order accurate Adams-Bashforth explicit differencing for both convection an diffusion terms. In order to evaluate more accurately a mass matrix, the well-known multi-pass algorithm was also adopted in this study. Numerical results obtained are compared through flow around a rectangular cylinder at Re = 22,000 with the experimental data and other existing numerical data.     

Central schemes for open‐channel flow

The resolution of the Saint‐Venant equations for modelling shock phenomena in open‐channel flow by using the second‐order central schemes of Nessyahu and Tadmor (NT) and Kurganov and Tadmor (KT) is presented. The performances of the two schemes that we have extended to the non‐homogeneous case and that of the classical first‐order Lax–Friedrichs (LF) scheme in predicting dam‐break and hydraulic jumps in rectangular open channels are investigated on the basis of different numerical and physical conditions. The efficiency and robustness of the schemes are tested by comparing model results with analytical or experimental solutions. Copyright © 2003 John Wiley & Sons, Ltd.     

An upstream flux‐splitting finite‐volume scheme for 2D shallow water equations

An upstream flux‐splitting finite‐volume (UFF) scheme is proposed for the solutions of the 2D shallow water equations. In the framework of the finite‐volume method, the artificially upstream flux vector splitting method is employed to establish the numerical flux function for the local Riemann problem. Based on this algorithm, an UFF scheme without Jacobian matrix operation is developed. The proposed scheme satisfying entropy condition is extended to be second‐order‐accurate using the MUSCL approach. The proposed UFF scheme and its second‐order extension are verified through the simulations of four shallow water problems, including the 1D idealized dam breaking, the oblique hydraulic jump, the circular dam breaking, and the dam‐break experiment with 45° bend channel. Meanwhile, the numerical performance of the UFF scheme is compared with those of three well‐known upwind schemes, namely the Osher, Roe, and HLL schemes. It is demonstrated that the proposed scheme performs remarkably well for shallow water flows. The simulated results also show that the UFF scheme has superior overall numerical performances among the schemes tested. Copyright © 2005 John Wiley & Sons, Ltd.     

A finite element thermally coupled flow formulation for phase‐change problems

A finite element, thermally coupled incompressible flow formulation considering phase‐change effects is presented. This formulation accounts for natural convection, temperature‐dependent material properties and isothermal and non‐isothermal phase‐change models. In this context, the full Navier–Stokes equations are solved using a generalized streamline operator (GSO) technique. The highly non‐linear phase‐change effects are treated with a temperature‐based algorithm, which provides stability and convergence of the numerical solution. The Boussinesq approximation is used in order to consider the temperature‐dependent density variation. Furthermore, the numerical solution of the coupled problem is approached with a staggered incremental‐iterative solution scheme, such that the convergence criteria are written in terms of the residual vectors. Finally, this formulation is used for the solutions of solidification and melting problems validating some numerical results with other existing solutions obtained with different methodologies. Copyright © 2000 John Wiley & Sons, Ltd.     

A simple strategy for constructing bounded convection schemes for unstructured grids

Due to the great geometrical flexibility, popularity for unstructured grid methods in fluid dynamics has been increasing in recent years. In parallel with this interest there is a need for bounded second or higher order convection schemes which can be implemented easily in the unstructured setting. In the present work a simple strategy for achieving convective boundedness in the context of a vertex‐centered unstructured finite volume algorithm is demonstrated. Testing is carried out on an inviscid oblique step problem using both structured and unstructured grid arrangements. Further testing for numerical diffusion is done using a distorted grid in a two dimensional channel. The proposed scheme is straightforward to implement and is found to perform well for the cases considered. The overall algorithm converges well and the limiter appears to introduce little extra numerical diffusion beyond that inherently present in the base scheme. Copyright © 2004 John Wiley & Sons, Ltd.     

Implicit weighted essentially non‐oscillatory schemes for the incompressible Navier–Stokes equations

A class of lower–upper/approximate factorization (LUAF) implicit weighted essentially non‐oscillatory (ENO; WENO) schemes for solving the two‐dimensional incompressible Navier–Stokes equations in a generalized co‐ordinate system is presented. The algorithm is based on the artificial compressibility formulation, and symmetric Gauss–Seidel relaxation is used for computing steady state solutions while symmetric successive overrelaxation is used for treating time‐dependent flows. WENO spatial operators are employed for inviscid fluxes and central differencing for viscous fluxes. Internal and external viscous flow test problems are presented to verify the numerical schemes. The use of a WENO spatial operator not only enhances the accuracy of solutions but also improves the convergence rate for the steady state computation as compared with using the ENO counterpart. It is found that the present solutions compare well with exact solutions, experimental data and other numerical results. Copyright © 1999 John Wiley & Sons, Ltd.     

Explicit finite volume non‐oscillatory schemes for 2D transient free‐surface flows

A class of high‐resolution non‐oscillatory shock‐capturing Roe, TVD and ENO explicit schemes in finite volume approach are presented for the computation of 2D unsteady rapidly varied open channel flows. In order to apply these schemes to simulate the hydraulic phenomena in field, the Strang‐type operator splitting technique is adopted to treat the flow with bottom slope and friction terms. Verifications of the proposed schemes are made by comparison with analytical solutions or experimental data, and very good agreements are obtained. To illustrate the efficiency and stability of the present algorithms, four typical problems of rapidly varied flows are solved and the results of different schemes are compared. It is demonstrated that the proposed method is accurate, robust and highly stable even in the flows with very strong discontinuites, which need no tuning of any adjustable parameter, such as artificial viscosity coefficient, as other methods do, and is a reliable mathematical modeling for 2D practical hydraulic engineering applications. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

A semi‐Lagrangian multi‐moment finite volume method with fourth‐order WENO projection

Numerical oscillation has been an open problem for high‐order numerical methods with increased local degrees of freedom (DOFs). Current strategies mainly follow the limiting projections derived originally for conventional finite volume methods and thus are not able to make full use of the sub‐cell information available in the local high‐order reconstructions. This paper presents a novel algorithm that introduces a nodal value‐based weighted essentially non‐oscillatory limiter for constrained interpolation profile/multi‐moment finite volume method (CIP/MM FVM) (Ii and Xiao, J. Comput. Phys., 222 (2007), 849–871) as an effort to pursue a better suited formulation to implement the limiting projection in schemes with local DOFs. The new scheme, CIP‐CSL‐WENO4 scheme, extends the CIP/MM FVM method by limiting the slope constraint in the interpolation function using the weighted essentially non‐oscillatory (WENO) reconstruction that makes use of the sub‐cell information available from the local DOFs and is built from the point values at the solution points within three neighboring cells, thus resulting a more compact WENO stencil. The proposed WENO limiter matches well the original CIP/MM FVM, which leads to a new scheme of high accuracy, algorithmic simplicity, and computational efficiency. We present the numerical results of benchmark tests for both scalar and Euler conservation laws to manifest the fourth‐order accuracy and oscillation‐suppressing property of the proposed scheme. Copyright © 2016 John Wiley & Sons, Ltd.     

A dual-porosity model for gas reservoir flow incorporating adsorption behaviour—Part II. Numerical algorithm and example applications

In the present study, an algorithm is presented for the dual-porosity model formulated in Part I of this series. The resultant flow equation with the dual-porosity formulation is of an integro-(partial) differential equation involving differential terms for the Darcy flow in large fractures and integrals in time for diffusion within matrix blocks. The algorithm developed here to solve this equation involves a step-by-step finite difference procedure combined with a quadrature scheme. The quadrature scheme, used for the integral terms, is based on the trapezoidal method which is of second-order precision. This order of accuracy is consistent with the first- and second-order finite difference approximations used here to solve the differential terms in the derived flow equation. In an approach consistent with many petroleum reservoir and groundwater numerical flow models, the example formulation presented uses a first-order implicit algorithm. A two-dimensional example is also demonstrated, with the proposed model and numerical scheme being directly incorporated into the commercial gas reservoir simulator SIMED II that is based on a fully implicit finite difference approach. The solution procedure is applied to several problems to demonstrate its performance. Results from the derived dual-porosity formulation are also compared to the classic Warren–Root model. Whilst some of this work confirmed previous findings regarding Warren–Root inaccuracies at early times, it was also found that inaccuracy can re-enter the Warren–Root results whenever there are changes in boundary conditions leading to transient variation within the domain.     

High‐order implicit Runge–Kutta time integrators for fluid‐structure interactions

This paper presents an approach to develop high‐order, temporally accurate, finite element approximations of fluid‐structure interaction (FSI) problems. The proposed numerical method uses an implicit monolithic formulation in which the same implicit Runge–Kutta (IRK) temporal integrator is used for the incompressible flow, the structural equations undergoing large displacements, and the coupling terms at the fluid‐solid interface. In this context of stiff interaction problems, the fully implicit one‐step approach presented is an original alternative to traditional multistep or explicit one‐step finite element approaches. The numerical scheme takes advantage of an arbitrary Lagrangian–Eulerian formulation of the equations designed to satisfy the geometric conservation law and to guarantee that the high‐order temporal accuracy of the IRK time integrators observed on fixed meshes is preserved on arbitrary Lagrangian–Eulerian deforming meshes. A thorough review of the literature reveals that in most previous works, high‐order time accuracy (higher than second order) is seldom achieved for FSI problems. We present thorough time‐step refinement studies for a rigid oscillating‐airfoil on deforming meshes to confirm the time accuracy on the extracted aerodynamics reactions of IRK time integrators up to fifth order. Efficiency of the proposed approach is then tested on a stiff FSI problem of flow‐induced vibrations of a flexible strip. The time‐step refinement studies indicate the following: stability of the proposed approach is always observed even with large time step and spurious oscillations on the structure are avoided without added damping. While higher order IRK schemes require more memory than classical schemes (implicit Euler), they are faster for a given level of temporal accuracy in two dimensions. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of fourth‒order difference method to prediction of turbulence control in two‒dimensional channel with surface roughness

The fourth‒order finite difference method is combined with the vorticity–streamfunction formulation in generalized co‒ordinates. Direct numerical simulations are performed for channel flows with and without surface roughness at a Reynolds number of 10⁴. The present results are in good agreement with those of the pseudospectral method with respect to the flow in a smooth channel. It is shown that the present method predicts well the precise change in the flow with the channel length and roughness height. The turbulence is generally weakened by the roughness. Laminarization is also accomplished under the appropriate condition. © 1997 John Wiley & Sons, Ltd.     

An extension of the SIMPLE based discontinuous Galerkin solver to unsteady incompressible flows

In this paper, we present a SIMPLE based algorithm in the context of the discontinuous Galerkin method for unsteady incompressible flows. Time discretization is done fully implicit using backward differentiation formulae (BDF) of varying order from 1 to 4. We show that the original equation for the pressure correction can be modified by using an equivalent operator stemming from the symmetric interior penalty (SIP) method leading to a reduced stencil size. To assess the accuracy as well as the stability and the performance of the scheme, three different test cases are carried out: the Taylor vortex flow, the Orr‐Sommerfeld stability problem for plane Poiseuille flow and the flow past a square cylinder. (1) Simulating the Taylor vortex flow, we verify the temporal accuracy for the different BDF schemes. Using the mixed‐order formulation, a spatial convergence study yields convergence rates of k + 1 and k in the L²‐norm for velocity and pressure, respectively. For the equal‐order formulation, we obtain approximately the same convergence rates, while the absolute error is smaller. (2) The stability of our method is examined by simulating the Orr–Sommerfeld stability problem. Using the mixed‐order formulation and adjusting the penalty parameter of the symmetric interior penalty method for the discretization of the viscous part, we can demonstrate the long‐term stability of the algorithm. Using pressure stabilization the equal‐order formulation is stable without changing the penalty parameter. (3) Finally, the results for the flow past a square cylinder show excellent agreement with numerical reference solutions as well as experiments. Copyright © 2015 John Wiley & Sons, Ltd.