Difference between revisions of "Step 1"

From CFD Benchmark
Jump to: navigation, search
m (Step 1 Description)
(Step 1 Description)
Line 20: Line 20:
 
</math>
 
</math>
  
with <math>u_0 = 1 \; \mathrm{m/s}$.
+
with <math>u_0 = 1 \; \mathrm{m/s}</math>.
 
The flow is thus initially composed of four vortices, one in each quarter of the box.
 
The flow is thus initially composed of four vortices, one in each quarter of the box.
 
The vortex turnover time is defined as <math>\tau_\mathrm{ref} = L_0/u_0 = 1\;\mathrm{s}</math>.
 
The vortex turnover time is defined as <math>\tau_\mathrm{ref} = L_0/u_0 = 1\;\mathrm{s}</math>.

Revision as of 05:42, 27 July 2020

Step 1 Description

The 2-D TGV test case has been already well documented in the scientific literature, e.g. [1][2]. Therefore, only a brief description is given here. This configuration is used as verification step, since an analytic solution exists~\cite{Taylor1937}.

An incompressible flow with constant density and viscosity is simulated.

The fluid is contained in a square box of dimension with periodic boundary conditions in both directions. The value with has been arbitrarily chosen in the present benchmark for dimensional codes. The numerical value for the kinematic viscosity is chosen as to get suitable conditions, as discussed below. The initial velocity components are prescribed as follows:

with . The flow is thus initially composed of four vortices, one in each quarter of the box. The vortex turnover time is defined as . For non-dimensional codes, the relevant parameter for this simulation is the vortex Reynolds number, in the present setup. It should be noted that the theoretical solution is only stable for sufficiently small Reynolds numbers, explaining why the particular value has been chosen here, as done in many previous studies using the TGV configuration.

A major interest of this test-case is that an analytical solution has been derived for such conditions~\cite{Taylor1937}:

All physical conditions are now completely defined and only numerical parameters remain to be set. Based on previously published studies, the recommended parameters for time and space discretization are and grid points in each direction. The flow must be simulated for a physical time of .

Aside suggestions

Notes

  1. A. Abdelsamie and G. Fru and F. Dietzsch and G. Janiga and D. Thévenin , Towards direct numerical simulations of low-Mach number turbulent reacting and two-phase flows using immersed boundaries, Comput. Fluids, 2016 (131-5, pp123--141)
  2. S. Laizet and E. Lamballais , High-order compact schemes for incompressible flows: A simple and efficient method with quasi-spectral accuracy, J. Comput. Phys., 2009 (228, pp5989--6015)