HyperWorks Solvers

/PROP/TYPE15 (POROUS)

/PROP/TYPE15 (POROUS)

Previous topic Next topic Expand/collapse all hidden text  

/PROP/TYPE15 (POROUS)

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

Block Format Keyword

/PROP/TYPE15 - Porous Solid Element Property Set

Description

Describes the porous solid element property set (extended Darcy's law).

Format

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

/PROP/TYPE15/prop_ID or /PROP/POROUS/prop_ID

prop_title

Blank Format

qa

qb

h

 

 

 

 

Por

 

 

 

 

 

 

 

 

R1

R2

R3

 

 

 

 

skew_IDr

Ihon

 

 

 

 

 

 

 

 

Itu

symbol_a_14

lmix

 

 

 

 

 

node_ID

 

 

 

 

 

 

 

 

 

hmtoggle_plus1Flag Definition

Field

Contents

SI Unit Example

prop_ID

Property identifier

(Integer, maximum 10 digits)

 

prop_title

Property title

(Character, maximum 100 characters)

 

qa

Quadratic bulk viscosity

Default = 10-20  (Real)

 

qb

Linear bulk viscosity

Default = 10-20  (Real)

 

h

Hourglass viscosity coefficient

Default = 0.10  (Real)

 

Por

Porosity

Default = 1.0  (Real)

 

R1

Specific resistance factor in direction 1

Default = 0.0  (Real)

 

R2

Specific resistance factor in direction 2

Default = 0.0  (Real)

 

R3

Specific resistance factor in direction 3

Default = 0.0  (Real)

 

skew_IDr

Skew identifier for resistance orthotropy

(Integer)

 

Ihon

Honeycomb substrate flag

Default = 0.0  (Integer)

=1: resistance applies only in direction 1 of skew_IDr and velocities are constrained in directions 2 and 3.

 

Itu

Turbulence imposed by medium flag (Comment 7)

Default = 0.0  (Integer)

= 0: turbulence is done according to model

= 1: turbulence is imposed by porous model

 

symbol_a_14

Turbulence coefficient for honeycomb substrate

Default = 0.1  (Real)

 

lmix

Turbulence mixing length

(Real)

symbol_m

node_ID

Master node identifier of rigid body modeling rigid substrate

Default = 0  (Integer)

 

hmtoggle_plus1Example (Air)

#RADIOSS STARTER

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

/PROP/POROUS/4

AIR (unit kg_m_s)

 

#                q_a                 q_b                   h

               1E-20               1E-20                   0

#                Por

                   1

#                 R1                  R2                  R3

                 400                 400                 400

# skew_IDr      Ihon

         0         0

#      Itu               alpha               I_mix

         0                   0                   0

#  node_ID

         0

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

#enddata

/END

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

hmtoggle_plus1Example (Gas)

#RADIOSS STARTER

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

/PROP/POROUS/13

GAS (unit kg_m_s)

 

#                q_a                 q_b                   h

                 1.1                 .05                   0

#                Por

                   0

#                 R1                  R2                  R3

                3500                   0                   0

# skew_IDr      Ihon

         3         1

#      Itu               alpha               I_mix

         1                 .05                .003

#  node_ID

     43218

#---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.The qa and qb default values are equal to 1.10 and 0.05 in non-CFD versions.
2.Porosity (relative volume of pores) is taken into account for 3D Euler materials only. Porosity is otherwise set to 1.
3.Force Fi = -mRij Vj is added to nodal force vector, Rij being the tension representation of the specific resistance factor in the global skew system.
4.If all resistance factors are set to 0, no computation is performed but the grid velocity can be tied to rigid body node_ID (Line 9).
5.Specific resistances along orthotropic axis are defined in Line 7 with respect to skew_IDr (by default global system).
6.Relevant only if material is material law 6 with turbulence.
7.If Itu =1, turbulence is imposed by porous model, and .
8.Fluid reaction force vector and moment are transmitted to specified rigid body; this option is available in 3D only.
9.Grid velocity is computed according to node_ID rigid body movement. This option is useful to impose such a movement to a grid, even for non-porous material.

See Also:

Skew and Frame (/SKEW & /FRAME)