HyperWorks Solvers

/MAT/LAW6 (K-EPS)

/MAT/LAW6 (K-EPS)

 Previous topic Next topic Expand/collapse all hidden text  

/MAT/LAW6 (K-EPS)

 Previous topic Next topic JavaScript is required for expanding text JavaScript is required for the print function  

Block Format Keyword

/MAT/LAW6 - k - E1 Turbulence Material (with k - E1)

Description

Describes the k - E1 turbulence viscous material for fluid.

Format

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

/MAT/LAW6/mat_ID/unit_ID or /MAT/K-EPS/mat_ID/unit_ID

mat_title

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C0

C1

C2

C3

 

 

Pmin

Psh

 

 

 

 

 

 

C4

C5

E0

 

 

 

 

SSL

 

 

 

 

 

 

csymbol_u

symbolk

symbole

Pr / Prt

 

 

c1ε

c2ε

c3ε

 

 

 

 

k

E

symbol_a_14

Xt

 

 
hmtoggle_plus1Flag Definition

Field

Contents

SI Unit Example

mat_ID

Material identifier

(Integer, maximum 10 digits)

 

unit_ID

Optional unit identifier

(Integer, maximum 10 digits)

 

mat_title

Material title

(Character, maximum 100 characters)

 

Initial density

(Real)

Reference density used in E.O.S (equation of state)

Default (Real)

Kinematic viscosity

(Real)

symbol_m2

C0

Constant parameters coefficient

(Real)

symbol_Pa

C1

Constant parameters coefficient

(Real)

symbol_Pa

C2

Constant parameters coefficient

(Real)

symbol_Pa

C3

Constant parameters coefficient

(Real)

symbol_Pa

Pmin

Pressure cutoff ( < 0 )

Default = -1030  (Real)

symbol_Pa

Psh

Pressure shift

(Real)

symbol_Pa

C4

Energy coefficient

(Real)

 

C5

Energy coefficient

(Real)

 

E0

Initial energy per unit volume

(Real)

fail_energy

Initial turbulent energy (first part)

(Real)

symbol_J

SSL

Subgrid scale length (first part)

Default = 1e+10  (Real)

symbol_m

csymbol_u

Turbulent viscosity coefficient (second part)

Default = 0.09  (Real)

 

symbolk

k diffusion coefficient (second part)

Default = 1.00  (Real)

 

symbole

Prandtl number of dissipation (second part)

Default = 1.30  (Real)

 

Pr / Prt

Laminar/turbulent Prandtl ratio (second part)

Default = 0.7/0.9  (Real)

 

c1ε

E1 equation coefficient 1 (third part)

Default = 1.440  (Real)

 

c2ε

E1 equation coefficient 2 (third part)

Default = 1.920  (Real)

 

c3ε

E1 equation coefficient 3 (third part)

Default = -0.375  (Real)

 

k

Kappa wall constant (fourth part)

Default = 0.4187  (Real)

 

E

E wall constant (fourth part)

Default = 9.7930  (Real)

 

symbol_a_14

excentration (fourth part)

Default = 0.5000  (Real)

 

Xt

Source term factor (fourth part)

(Real)

 
hmtoggle_plus1Example (Gas)

#RADIOSS STARTER

/UNIT/1

unit for mat

                 kg                   m                   s

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

/MAT/K-EPS/4/1

GAS

#              RHO_I

               .3828

#                KNU

             1.05E-4

#                 C0                  C1                  C2                  C3

                   0                   0                   0                   0

#               Pmin                 Psh

                   0              101320

#                 C4                  C5                  E0

                 0.4                 0.4              253300

#            RHO0_K0                 SSL

                  20                   0

#               C_MU               SIG_k             SIG_EPS         P_R_ON_P_RT

                   0                   0                   0

#                C_1                 C_2                 C_3

                   0                   0                   0

#              KAPPA                   E              ALPHA                GSI_T

                   0                   0                   0                   0

/ALE/MAT/4

#     Modif. factor.

                   0

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

#enddata

/END

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
hmtoggle_plus1Comments
1.If

Where, Sij is the deviatoric stress tensor and eij is the deviatoric strain tensor.

If the element is connected to a boundary condition, a turbulent boundary layer model is used:

Where, κ is the turbulent kinetic energy.

If the ratio between the laminar and the turbulent Prantl numbers is higher than Pr / Prt, then:

For laminar flow:

For turbulent flow:

Where, E1 is the turbulent dissipation and it is calculated using the following equations:

With,

(turbulent viscosity)

Where, is the material velocity and w is the grid velocity.

2.In case of a perfect gas:

and

with symbol_y being the perfect gas constant.

3.In case of an incompressible gas:

  and

with c being the sound velocity.

4.In case of a linear material with a volumetric dilatation:

and

and

then:

If P = cst = 0,  then

so

Where, symbol_u is the dilatation coefficient, symbol_u < 0 means dilatation, and in this case, the parameters C2 and C3 will not be taken into account.

5.If using Law 6 coupled with Law 37 for liquid phase (without gas phase), the compatibility of the liquid EOS is as follows:
•  for Law 37
•law6_eq for Law 6

with

then,

6.If using Law 6 coupled with Law 37 for gas phase (without liquid phase), the compatibility of the gas EOS is as follows:

for Law 37

for Law 6

with and

then,

Where, E is the energy per unit volume and e is the energy per unit mass.

7.All thermal data () can be defined with keyword /HEAT.

See Also:

Material Compatibility