12.9.4 Inlet Boundary Conditions for the LES Model

This section describes the three algorithms available in FLUENT to model the fluctuating velocity at velocity inlet boundaries.

No Perturbations

The stochastic components of the flow at the velocity-specified inlet boundaries are neglected if the No Perturbations option is used. In such cases, individual instantaneous velocity components are simply set equal to their mean velocity counterparts. This option is suitable only when the level of turbulence at the inflow boundaries is negligible or does not play a major role in the accuracy of the overall solution.

Vortex Method

To generate a time-dependent inlet condition, a random 2D vortex method is considered. With this approach, a perturbation is added on a specified mean velocity profile via a fluctuating vorticity field (i.e. two-dimensional in the plane normal to the streamwise direction). The vortex method is based on the Lagrangian form of the 2D evolution equation of the vorticity and the Biot-Savart law. A particle discretization is used to solve this equation. These particles, or "vortex points'' are convected randomly and carry information about the vorticity field. If is the number of vortex points and is the area of the inlet section, the amount of vorticity carried by a given particle is represented by the circulation and an assumed spatial distribution :

 (12.9-24) (12.9-25)

where is the turbulence kinetic energy. The parameter provides control over the size of a vortex particle. The resulting discretization for the velocity field is given by
 (12.9-26)

Where is the unit vector in the streamwise direction. Originally [ 327], the size of the vortex was fixed by an ad hoc value of . To make the vortex method generally applicable, a local vortex size is specified through a turbulent mixing length hypothesis. is calculated from a known profile of mean turbulence kinetic energy and mean dissipation rate at the inlet according to the following:
 (12.9-27)

where . To ensure that the vortex will always belong to resolved scales, the minimum value of in Equation  12.9-27 is bounded by the local grid size. The sign of the circulation of each vortex is changed randomly each characteristic time scale . In the general implementation of the vortex method, this time scale represents the time necessary for a 2D vortex convected by the bulk velocity in the boundary normal direction to travel along times its mean characteristic 2D size ( ), where is fixed equal to 100 from numerical testing. The vortex method considers only velocity fluctuations in the plane normal to the streamwise direction.

In FLUENT however, a simplified linear kinematic model (LKM) for the streamwise velocity fluctuations is used [ 231]. It is derived from a linear model that mimics the influence of the two-dimensional vortex in the streamwise mean velocity field. If the mean streamwise velocity is considered as a passive scalar, the fluctuation resulting from the transport of by the planar fluctuating velocity field is modeled by
 (12.9-28)

where is the unit vector aligned with the mean velocity gradient . When this mean velocity gradient is equal to zero, a random perturbation can be considered instead.

 Since the vortex method theory is based on the modification of the velocity field normal to the streamwise direction, it is imperative that the user creates an inlet plane normal (or as close as possible) to the streamwise velocity direction.

Spectral Synthesizer

The spectral synthesizer provides an alternative method of generating fluctuating velocity components. It is based on the random flow generation technique originally proposed by Kraichnan [ 186] and modified by Smirnov et al. [ 338]. In this method, fluctuating velocity components are computed by synthesizing a divergence-free velocity-vector field from the summation of Fourier harmonics. In the implementation in FLUENT, the number of Fourier harmonics is fixed to 100.

Previous: 12.9.3 Subgrid-Scale Models
Up: 12.9 Large Eddy Simulation
Next: 12.10 Near-Wall Treatments for