What number can i use for Turbulent Kinetic Energy and Turb. Dissipation Rate, what is the criteria for chosing it. In my case i have solar radiation raising the temperature of the wall to 75 deg.C in a completely calm day and am simulating for buoyancy flows of air.

The following descriptions may help you: In some situations, it is appropriate to specify a uniform value of the turbulence quantity at the boundary where inflow occurs. Examples are fluid entering a duct, far-field boundaries, or even fully-developed duct flows where accurate profiles of turbulence quantities are unknown.

In most turbulent flows, higher levels of turbulence are generated within shear layers than enter the domain at flow boundaries, making the result of the calculation relatively insensitive to the inflow boundary values. Nevertheless, caution must be used to ensure that boundary values are not so unphysical as to contaminate your solution or impede convergence. This is particularly true of external flows where unphysically large values of effective viscosity in the free stream can "swamp'' the boundary layers.

You can use the turbulence specification methods described above to enter uniform constant values instead of profiles. Alternatively, you can specify the turbulence quantities in terms of more convenient quantities such as turbulence intensity, turbulent viscosity ratio, hydraulic diameter, and turbulence length scale. These quantities are discussed further in the following sections.

Turbulence Intensity

The turbulence intensity, , is defined as the ratio of the root-mean-square of the velocity fluctuations, , to the mean flow velocity, .

A turbulence intensity of 1% or less is generally considered low and turbulence intensities greater than 10% are considered high. Ideally, you will have a good estimate of the turbulence intensity at the inlet boundary from external, measured data. For example, if you are simulating a wind-tunnel experiment, the turbulence intensity in the free stream is usually available from the tunnel characteristics. In modern low-turbulence wind tunnels, the free-stream turbulence intensity may be as low as 0.05%.

For internal flows, the turbulence intensity at the inlets is totally dependent on the upstream history of the flow. If the flow upstream is under-developed and undisturbed, you can use a low turbulence intensity. If the flow is fully developed, the turbulence intensity may be as high as a few percent. The turbulence intensity at the core of a fully-developed duct flow can be estimated from the following formula derived from an empirical correlation for pipe flows:

(7.2-1)

At a Reynolds number of 50,000, for example, the turbulence intensity will be 4%, according to this formula.

Turbulence Length Scale and Hydraulic Diameter

The turbulence length scale, , is a physical quantity related to the size of the large eddies that contain the energy in turbulent flows.

In fully-developed duct flows, is restricted by the size of the duct, since the turbulent eddies cannot be larger than the duct. An approximate relationship between and the physical size of the duct is

(7.2-2)

where is the relevant dimension of the duct. The factor of 0.07 is based on the maximum value of the mixing length in fully-developed turbulent pipe flow, where is the diameter of the pipe. In a channel of non-circular cross-section, you can base on the hydraulic diameter.

If the turbulence derives its characteristic length from an obstacle in the flow, such as a perforated plate, it is more appropriate to base the turbulence length scale on the characteristic length of the obstacle rather than on the duct size.

It should be noted that the relationship of Equation 7.2-2, which relates a physical dimension ( ) to the turbulence length scale ( ), is not necessarily applicable to all situations. For most cases, however, it is a suitable approximation.

Guidelines for choosing the characteristic length or the turbulence length scale for selected flow types are listed below:

For fully-developed internal flows, choose the Intensity and Hydraulic Diameter specification method and specify the hydraulic diameter in the Hydraulic Diameter field.

For flows downstream of turning vanes, perforated plates, etc., choose the Intensity and Hydraulic Diameter method and specify the characteristic length of the flow opening for in the Hydraulic Diameter field.

For wall-bounded flows in which the inlets involve a turbulent boundary layer, choose the Intensity and Length Scale method and use the boundary-layer thickness, , to compute the turbulence length scale, , from . Enter this value for in the Turbulence Length Scale field.

Turbulent Viscosity Ratio

The turbulent viscosity ratio, , is directly proportional to the turbulent Reynolds number ( ). is large (on the order of 100 to 1000) in high-Reynolds-number boundary layers, shear layers, and fully-developed duct flows. However, at the free-stream boundaries of most external flows, is fairly small. Typically, the turbulence parameters are set so that .

To specify quantities in terms of the turbulent viscosity ratio, you can choose Turbulent Viscosity Ratio (for the Spalart-Allmaras model) or Intensity and Viscosity Ratio (for the - models, the - models, or the RSM).

Relationships for Deriving Turbulence Quantities

--------------------------------------------------------------------------------

To obtain the values of transported turbulence quantities from more convenient quantities such as , , or , you must typically resort to an empirical relation. Several useful relations, most of which are used within FLUENT, are presented below.

Estimating Modified Turbulent Viscosity from Turbulence Intensity and Length Scale

To obtain the modified turbulent viscosity, , for the Spalart-Allmaras model from the turbulence intensity, , and length scale, , the following equation can be used:

(7.2-3)

This formula is used in FLUENT if you select the Intensity and Hydraulic Diameter specification method with the Spalart-Allmaras model. is obtained from Equation 7.2-2.

Estimating Turbulent Kinetic Energy from Turbulence Intensity

The relationship between the turbulent kinetic energy, , and turbulence intensity, , is

(7.2-4)

where is the mean flow velocity.

This relationship is used in FLUENT whenever the Intensity and Hydraulic Diameter, Intensity and Length Scale, or Intensity and Viscosity Ratio method is used instead of specifying explicit values for and .

Estimating Turbulent Dissipation Rate from a Length Scale

If you know the turbulence length scale, , you can determine from the relationship

(7.2-5)

where is an empirical constant specified in the turbulence model (approximately 0.09). The determination of was discussed previously.

This relationship is used in FLUENT whenever the Intensity and Hydraulic Diameter or Intensity and Length Scale method is used instead of specifying explicit values for and .

Estimating Turbulent Dissipation Rate from Turbulent Viscosity Ratio

The value of can be obtained from the turbulent viscosity ratio and using the following relationship:

(7.2-6)

where is an empirical constant specified in the turbulence model (approximately 0.09).

This relationship is used in FLUENT whenever the Intensity and Viscosity Ratio method is used instead of specifying explicit values for and .

Estimating Turbulent Dissipation Rate for Decaying Turbulence

If you are simulating a wind-tunnel situation in which the model is mounted in the test section downstream of a grid and/or wire mesh screens, you can choose a value of such that

(7.2-7)

where is the approximate decay of you wish to have across the flow domain (say, 10% of the inlet value of ), is the free-stream velocity, and is the streamwise length of the flow domain. Equation 7.2-7 is a linear approximation to the power-law decay observed in high-Reynolds-number isotropic turbulence. Its basis is the exact equation for in decaying turbulence, .

If you use this method to estimate , you should also check the resulting turbulent viscosity ratio to make sure that it is not too large, using Equation 7.2-6.

Although this method is not used internally by FLUENT, you can use it to derive a constant free-stream value of that you can then specify directly by choosing K and Epsilon in the Turbulence Specification Method drop-down list. In this situation, you will typically determine from using Equation 7.2-4.

[ohmsree]

is there any significant difference in results when the mesh size of a turbulent pipe is uniform through out than with more meshes at the wall and less mesh at the axis??

[adriant]

Thank you very much, could you detail formulas (7.2-1) an so on?