How to estimate the transport properties: dynamic viscosity, thermal conductivity and diffusivity of a gas, if measurements are not available?
Tagged: 16, cfx, fluid-dynamics, General, General - CFX, material-properties
-
-
January 25, 2023 at 7:16 am
FAQ
ParticipantSome gases are rare or have low vapour pressure and so measurements of transport properties may not be available from the literature. Dynamic Viscosity ============== Transport properties can still be estimated based on the kinetic theory of gases. For example, Poling. Prausnitz and O’Connell, “The Properties of Gases and Liquids” lists critical constants for many substances. Then from Bird, Stewart and Lightfoot, “Transport Phenomena”, Wiley NY 1960, p22-23, for spherical non-polar molecules, we can estimate the Lennard-Jones parameters: epsilon/k = 0.77*Tc where the critical temperature, Tc is in K, the characteristic diameter, sigma in Angstroms, = 0.841*Vc^(1/3) with the critical volume, Vc in cm3/mol (or sigma = 2.44(Tc/pc)^(1/3) with the critical pressure, pc, in atm). The dynamic viscosity for a gas at low density is then given by: Mu = 2.6693×10^-6*(MT)^(1/2)/sigma^2*Omega in kg/m/s, where M is the molecular weight in kg/kmol and Omega can be found as a function of epsilon /k in Table B-2, p746. Specific Heat Capacity ================= The specific heat at constant pressure, Cp, is usually easily found in the literature, for example, Alexander Burcat and Branko Ruscic, Ideal Gas Thermochemical Database with updates from Active Thermochemical Tables in ftp://ftp.technion.ac.il/pub/supported/aetdd/thermodynamics mirrored at http://garfield.chem.elte.hu/Burcat/burcat.html has NASA polynomials for many materials in the form of two-range NASA polynomials (S. Gordon and B.J. McBride, “Computer Program for Calculation of Complex Chemical Equilibrium Composition, Rocket Performance, Incident and Reflected Shocks and Chapman-Jouguet Detonations”, NASA SP-273 (1971)),. So for example, THERMO 200.0 1000.0 6000.0 Br J6/82BR 1 0 0 0G 200.000 6000.000 1 2.08851053E+00 7.12118611E-04-2.70003073E-07 4.14986299E-11-2.31188294E-15 2 1.28568767E+04 9.07351144E+00 2.48571711E+00 1.50647525E-04-5.37267333E-07 3 7.20921065E-10-2.50205558E-13 1.27092168E+04 6.86030804E+00 1.34535890E+04 4 END The lower polynomial is valid for 200K < T < 1000K: Cp = R/M*(a1+a2*T+a3*T^2+a4*T^3+a5*T^4), where a1 = 2.48571711E+00 a2 = 1.50647525E-04 a3 = -5.37267333E-07 a4 = 7.20921065E-10 a5 = -2.50205558E-13 and the upper for 1000K < T < 6000K a1 = 2.08851053E+00 a2 = 7.12118611E-04 a3 = -2.70003073E-07 a4 = 4.14986299E-11 a5 = -2.31188294E-15 and R is the universal gas constant (8314.41 J/kmol/K) and M is the molecular weight. Thermal Conductivity ================ You can estimate k from Bird, Stewart and Lightfoot, p257. the Eucken formula gives k = (Cp + 5/4*R/M)*Mu, which for monotonic gases can be simplified to k = 15/4*(R/M)*Mu.
-

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 High-Mounted Stop Lamp (CHMSL)
This video article demonstrates how to create a reflector for a center high-mounted 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 K-Omega) turbulence model offers a flexible, robust, general-purpose 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
- How can I create a Cell Register from a Cell Zone?
- ANSYS System Coupling: Two Way Fluid Structure Interaction – Part 1
- Delete or Deactivate Zone in Fluent
- Left-handed faces troubleshooting
- Running Python Script from Workbench
- ANSYS Polyflow: Adaptive Meshing Based on Contact
- How to overcome the model information incompatible with incoming mesh error?
- Check CPU Time in ANSYS FLUENT
- Apply Custom Material Properties in Fluent
© 2023 Copyright ANSYS, Inc. All rights reserved.