Dynamical modeling of sub-grid scales in 2D turbulence

被引:18
作者
Laval, JP
Dubrulle, B
Nazarenko, S
机构
[1] CEA, DAPNIA, F-91191 Gif Sur Yvette, France
[2] CEA, DAPNIA, CNRS, URA 285, F-91191 Gif Sur Yvette, France
[3] CNRS, URA 285, Observ Midi Pyrenees, F-31400 Toulouse, France
[4] Univ Warwick, Inst Math, Coventry CV4 7AL, W Midlands, England
[5] Univ Arizona, Dept Math, Tucson, AZ 85721 USA
来源
PHYSICA D | 2000年 / 142卷 / 3-4期
关键词
sub-grid scales; two-fluid simulation; pseudo-spectral simulation;
D O I
10.1016/S0167-2789(00)00060-9
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We develop a new numerical method which treats resolved and sub-grid scales as two different fluid components evolving according to their own dynamical equations. These two fluids are nonlinearly interacting and can be transformed one into another when their scale becomes comparable to the grid size. Equations describing the two-fluid dynamics were rigorously derived from Euler equations [B. Dubrulle, S. Nazarenko, Physica D 110 (1997) 123-138] and they do not involve any adjustable parameters. The main assumption of such a derivation is that the large-scale vortices are so strong that they advect the sub-grid scales as a passive scalar, and the interactions of small scales with small and intermediate scales can be neglected. As a test for our numerical method, we performed numerical simulations of 2D turbulence with a spectral gap, and we found a good agreement with analytical results obtained for this case by Nazarenko and Laval [Non-local 2D turbulence and passive scalars in Batchelor's regime, J. Fluid Mech., in press]. We used the two-fluid method to study three typical problems in 2D dynamics of incompressible fluids: decaying turbulence, vortex merger and forced turbulence. The two-fluid simulations performed on at 128(2) and 256(2) resolution were compared with pseudo-spectral simulations using hyperviscosity performed at the same and at much higher resolution. This comparison shows that performance of the two-fluid method is much better than one of the pseudo-spectral method at the same resolution and comparable computational cost. The most significant improvement is observed in modeling of the small-scale component, so that effective inertial interval increases by about two decades compared to the high-resolution pseudo-spectral method. Using the two-fluid method, we demonstrated that the k(-3) tail always exists for the energy spectrum, although its amplitude is slowly decreasing in decaying turbulence. (C) 2000 Elsevier Science B.V. All rights reserved.
引用
收藏
页码:231 / 253
页数:23
相关论文
共 33 条
[1]   Some properties of two-dimensional inverse energy cascade dynamics [J].
Babiano, A ;
Dubrulle, B ;
Frick, P .
PHYSICAL REVIEW E, 1997, 55 (03) :2693-2706
[2]   Self-similarity of decaying two-dimensional turbulence [J].
Bartello, P ;
Warn, T .
JOURNAL OF FLUID MECHANICS, 1996, 326 :357-372
[3]  
Batchelor G K., 1969, Phys. Fluids, V12, pII, DOI [10.1063/1.1692443, DOI 10.1063/1.1692443]
[4]  
Birdsall C., 1985, PLASMA PHYSICS VIA C
[5]   SPECTRAL EXPONENTS OF ENSTROPHY CASCADE IN STATIONARY 2-DIMENSIONAL HOMOGENEOUS TURBULENCE [J].
BORUE, V .
PHYSICAL REVIEW LETTERS, 1993, 71 (24) :3967-3970
[6]   INVERSE ENERGY CASCADE IN STATIONARY 2-DIMENSIONAL HOMOGENEOUS TURBULENCE [J].
BORUE, V .
PHYSICAL REVIEW LETTERS, 1994, 72 (10) :1475-1478
[7]   THE DYNAMICS OF FREELY DECAYING TWO-DIMENSIONAL TURBULENCE [J].
BRACHET, ME ;
MENEGUZZI, M ;
POLITANO, H ;
SULEM, PL .
JOURNAL OF FLUID MECHANICS, 1988, 194 :333-349
[9]   Incremental unknowns, multilevel methods and the numerical simulation of turbulence [J].
Dubois, T ;
Jauberteau, F ;
Temam, R .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 1998, 159 (1-2) :123-189
[10]   Interaction of turbulence and large-scale vortices in incompressible 2D fluids [J].
Dubrulle, B ;
Nazarenko, S .
PHYSICA D, 1997, 110 (1-2) :123-138