How is the dissipation rate calculated in the initialization panel from the turbulent intensity and viscosity ratio specified at the boundary?
Tagged: 16, BCs & Interfaces, fluent, fluiddynamics, General, materials, turbulence


January 25, 2023 at 7:16 amFAQParticipant
Let’s say you have specified the following conditions at a pressureinlet boundary: (Operating Pressure = 0 Pa) Total Pressure, Pt = 51056 Pa Supersonic/Initial Pressure, Ps = 27364 Pa Total Temperature, Tt = 793.2 K Turbulent Intensity, I = 0.05 (5%) Turbulent Viscosity Ratio, Mut/Mu, Rv = 10 The material used is Air, with Density specified as ‘Ideal Gas’ and Viscosity specified using ‘Sutherland’ method. Standard kepsilon turbulence model is being used. When you go to the initialization panels and use that boundary as the source (Compute From), the following values are given by Fluent: Temperature, T = 663.9422 K X Velocity, Vx = 510.0763 m/s Turbulent Kinetic Energy, k = 975.6668 (m/s)^2 Turbulent Dissipation Rate, e = 6.874611e+07 m2/s3 From display contour, the density and viscosity are as follows: Density, rho = 0.1435854 kg/m3 Molecular Viscosity, mu = 3.22212e05 kg/(ms) Using ideal gas, isentropic relations, and the above equations, we can confirm that all the initialization values given by Fluent are correct, except for the value of Turbulent Dissipation Rate, e. The Turbulent Kinetic Energy, k, and Turbulent Dissipation Rate, e, are given by the following equations: k = 1.5*(Uavg*I)^2 e = rho*Cmu*(k*k/Mu)*(1/Rv) For Standard kepsilon turbulence model the model constant, Cmu = 0.09 The value of e computed using the above equation is 3.817808E+07 m2/s3. Why is the value of e in the initialization panel different? Explanation For initialization only, the value of mu used in the calculation of epsilon above is the constant value of mu for standard air and NOT the value computed using the Sutherland formula above. It is true that for consistencies, the value of mu computed by the Sutherland formula should be used instead, but most likely the flow is insensitive to the initial value of e, hence the procedure used by Fluent is acceptable. As verification, using mu=1.7894e5 kg/(ms), the value of Turbulent Dissipation Rate is given as e = 6.874616E+07 m2/s3 (which is very close to the value given by the initialization panel). Once the iteration is started, the correct value of mu computed using the Sutherland formula is used.

Introducing Ansys Electronics Desktop on Ansys Cloud
The Watch & Learn video article provides an overview of cloud computing from Electronics Desktop and details the product licenses and subscriptions to ANSYS Cloud Service that are...
How to Create a Reflector for a Center HighMounted Stop Lamp (CHMSL)
This video article demonstrates how to create a reflector for a center highmounted stop lamp. Optical Part design in Ansys SPEOS enables the design and validation of multiple...
Introducing the GEKO Turbulence Model in Ansys Fluent
The GEKO (GEneralized KOmega) turbulence model offers a flexible, robust, generalpurpose approach to RANS turbulence modeling. Introducing 2 videos: Part 1 provides background information on the model and a...
Postprocessing on Ansys EnSight
This video demonstrates exporting data from Fluent in EnSight Case Gold format, and it reviews the basic postprocessing capabilities of EnSight.
 ANSYS Fluent: Introduction to the GEKO Turbulence Model Part I
 ANSYS Polyflow: Adaptive Meshing Based on Contact
 Delete or Deactivate Zone in Fluent
 Apply Custom Material Properties in Fluent
 AeroMechanical Simulation of Turbomachinery Blading
 What is meant by Warning: Flow boundary zone 18 is adjacent to a solid zone?
 ANSYS System Coupling: Two Way Fluid Structure Interaction – Part 1
 Check CPU Time in ANSYS FLUENT
 Predict Gearbox Lubrication, Oil Temperature and Churning Losses using CFD Simulation
 Running Python Script from Workbench
© 2023 Copyright ANSYS, Inc. All rights reserved.