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
|
|
lmix
|
|
|
|
|
|
node_ID
|
|
|
|
|
|
|
|
|
|
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
|
|
|
Turbulence coefficient for honeycomb substrate
Default = 0.1 (Real)
|
|
lmix
|
Turbulence mixing length
(Real)
|
|
node_ID
|
Master node identifier of rigid body modeling rigid substrate
Default = 0 (Integer)
|
|
|
#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----|
|
#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----|
|
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)