
LARGE EDDY SIMULATIONS

Researcher: 
Lars Davidson lada@chalmers.se 
Supervisor:   
Cooperation:   
Sponsors:   
Publications:  [112], see references below 
Start of project:  spring 1996 
THE PROJECT The standard dynamic subgrid model (the Germano model) has numerical stability problems. The remedy is to average in some homogeneous flow direction(s) or to introduce some artificial clipping. Thus this type of models does not seem to be applicable to real threedimensional flow without introducing ad hoc user modifications. In the present study a new oneequation subgrid model is presented which eliminates the need of this type of user modification. The idea is to include all local dynamic information into the source terms of the transport equation of k_{sgs}. In this way the numerical algorithm can accept local, wildly fluctuating dynamic coefficients. Negative values of the dynamic coefficient in the production term are permitted, thus allowing for back scatter. In the momentum equations, a homogeneously constant value <C>_{xyz} (varying in time) is used, which is computed by requiring that, when integrating over the whole computational domain, <C>_{xyz} should give the same production as the local value C. The model is applied to recirculating flow in an enclosure. The new model is shown to give better agreement with experimental data, to be considerably more stable numerically, and, as a consequence, to be computationally cheaper than the traditional, spanwiseaveraging Germano model. The idea is to include all local dynamic information through the source terms of the transport equation for k_{sgs}. This is probably physically more sound since large local variations in C appear only in the source term, and the effect of the large fluctuations in the dynamic coefficients will be smoothed out in a natural way. In this way, it turns out that the need to restrict or limit the dynamic coefficient is eliminated altogether. The spatial variation of C is included via the production term in the modelled k_{sgs} equation. In this way, back scatter is taken into account in an indirect way. Although it is not fed directly back to the resolved flow, it influences the resolved flow via the kinetic subgrid energy in two ways. First, a negative production reduces k_{sgs} and thereby the subgrid viscosity. Second, a negative production also reduces <C>_{xyz}. THE NUMERICAL METHOD An implicit, twostep time advancement method is used. When the filtered NavierStokes equation a velocity field is obtained which does not satisfy continuity. An intermediate velocity field is computed by subtracting the implicit part of the pressure gradient, and the resulting Poisson equation is solved employing an efficient multigrid method. Central differencing is used for the convective terms in all equations, including the k_{sgs} equation. For more details, see [4,5]. A 64processor Origin 2000 Silicon Graphics machine was used. The code has been parallelized by Zacharov [6]. For a backwardfacing step (maximum CFL approximately 2), the achieved speedup on 4 (8) processors is 3.5 (6.3) on a 1.3 million node mesh Each time step requires 28 seconds of CPU time on eight processors. REFERENCES  click to view postscript file

This page, LARGE EDDY SIMULATIONS: DEVELOPMENT OF ONEEQUATION SUBGRID MODELS, should be part of a frames system at www.tfd.chalmers.se/~lada/projects/proind.html by Webbmaster: IngaLena Ljungström ilj@flowsim.se 