Direct Numerical Simulation of Turbulence at Lower Costs

Direct Numerical Simulation (DNS) is the most accurate, but also the most expensive, way of computing turbulent flow. To cut the costs of DNS we consider a family of second-order, explicit one-leg time-integration methods and look for the method with the best linear stability properties. It turns out that this method requires about two times less computational effort than Adams–Bashforth. Next, we discuss a fourth-order finite-volume method that is constructed as the Richardson extrapolate of a classical second-order method. We compare the results of this fourth-order method and the underlying second-order method for a DNS of the flow in a cubical driven cavity at Re= 104. Experimental results are available for comparison. For this example, the fourth-order results are clearly superior to the second-order results, whereas their computational effort is about twenty times less. With the improved simulation method, a DNS of a turbulent flow in a cubical lid-driven flow at Re = 50,000 and a DNS of a turbulent flow past a square cylinder at Re = 22,000 are performed.

[1]  Ronald D. Joslin,et al.  Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate , 1992 .

[2]  Shmuel Einav,et al.  A laser-Doppler velocimetry study of ensemble-averaged characteristics of the turbulent near wake of a square cylinder , 1995, Journal of Fluid Mechanics.

[3]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[4]  A. Veldman,et al.  Playing with nonuniform grids , 1992 .

[5]  I. S. Gartshore,et al.  Spanwise correlations of pressure on a rigid square section cylinder , 1992 .

[6]  Joel H. Ferziger,et al.  Status of Large Eddy Simulation: Results of a Workshop , 1997 .

[7]  Parviz Moin,et al.  Direct simulations of turbulent flow using finite-difference schemes , 1991 .

[8]  B. Lee.  The effect of turbulence on the surface pressure field of a square prism , 1975, Journal of Fluid Mechanics.

[9]  I. Gustafsson.  A class of first order factorization methods , 1978 .

[10]  H. Fasel,et al.  Numerical Simulation of Instability and Transition in Boundary Layer Flows , 1990 .

[11]  S. Eisenstat.  Efficient Implementation of a Class of Preconditioned Conjugate Gradient Methods , 1981 .

[12]  Fourth order finite difference and multigrid methods for modeling instabilities in flat plate boundary layers—2-D and 3-D approaches , 1994 .

[13]  Peter R. Voke,et al.  Flow Past a Square Cylinder: Test Case LES2 , 1997 .

[14]  J. Koseff,et al.  Reynolds number and end‐wall effects on a lid‐driven cavity flow , 1989 .

[15]  Parviz Moin,et al.  Zonal Embedded Grids for Numerical Simulations of Wall-Bounded Turbulent Flows , 1996 .

[16]  Parviz Moin,et al.  An improvement of fractional step methods for the incompressible Navier-Stokes equations , 1991 .

[17]  Werner Liniger,et al.  Contractive methods for stiff differential equations part I , 1978 .

[18]  Arthur Veldman,et al.  Data-Parallel Solution of the Incompressible Navier-Stokes Equations , 1996 .

[19]  Arthur Veldman,et al.  Direct Numerical Simulation of Driven Cavity Flows , 1993 .