*HEADING
Test for fluid exchange bewteen fluid cavities filled with hydraulic fluids.

This input deck tests the following six types of fluid exchange behaviors for hydraulic fluids:
1. MASS FLUX
2. VOLUME FLUX
3. BULK VISCOSITY
4. FABRIC LEAKAGE
5. MASS RATE LEAKAGE
6. VOLUME RATE LEAKAGE

All the above six fluid exchange behaviors are tested with single cavity and dual cavity.

The surface based fluid cavities in all the test cases are considered to be enclosed within a cubic block of fixed volume=0.125 m^3. Each face of the cubic block is modeled using single SFM3D4R element.

Effective area of 0.02 is considered for bulk viscosity type exchange behavior
Effective area based on surface is considered for fabric leakage type fluid exchange behaviors.

Mass flow rate dependence on average absolute pressure is considered for mass rate leakage fluid exchange behavior.
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16
**#############################################################################################################
**      TEST CASE 1: Single Cavity-Mass Flux Type Exchange Property
**#############################################################################################################
*SYSTEM
0, 0, 0
*NODE, NSET=SINGLECAVITY_M_FLUX
11, 0., 0., 0.
12, 0.5, 0., 0.
13, 0.5, 0.5, 0.
14, 0., 0.5, 0.
15, 0., 0., 0.5
16, 0.5, 0., 0.5
17, 0.5, 0.5, 0.5
18, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SINGLECAVITY_M_FLUX_REF
199, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SINGLECAVITY_M_FLUX
11, 12, 13, 14, 11
12, 18, 17, 16, 15
13, 15, 16, 12, 11
14, 16, 17, 13, 12
15, 17, 18, 14, 13
16, 18, 15, 11, 14
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SINGLECAVITY_M_FLUX, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SINGLECAVITY_M_FLUX_SURF, TYPE=ELEMENT
SINGLECAVITY_M_FLUX, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SINGLECAVITY_M_FLUX_CAV , REF NODE=SINGLECAVITY_M_FLUX_REF, SURFACE=SINGLECAVITY_M_FLUX_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=SINGLECAVITY_M_FLUX, PROPERTY=M_FLUX
199
**#############################################################################################################
**      TEST CASE 2: Single Cavity-Volume Flux Type Exchange Property
**#############################################################################################################
*SYSTEM
0, 1, 0
*NODE, NSET=SINGLECAVITY_V_FLUX
21, 0., 0., 0.
22, 0.5, 0., 0.
23, 0.5, 0.5, 0.
24, 0., 0.5, 0.
25, 0., 0., 0.5
26, 0.5, 0., 0.5
27, 0.5, 0.5, 0.5
28, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SINGLECAVITY_V_FLUX_REF
299, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SINGLECAVITY_V_FLUX
21, 22, 23, 24, 21
22, 28, 27, 26, 25
23, 25, 26, 22, 21
24, 26, 27, 23, 22
25, 27, 28, 24, 23
26, 28, 25, 21, 24
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SINGLECAVITY_V_FLUX, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SINGLECAVITY_V_FLUX_SURF, TYPE=ELEMENT
SINGLECAVITY_V_FLUX, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SINGLECAVITY_V_FLUX_CAV , REF NODE=SINGLECAVITY_V_FLUX_REF, SURFACE=SINGLECAVITY_V_FLUX_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=SINGLECAVITY_V_FLUX, PROPERTY=V_FLUX
299
**#############################################################################################################
**      TEST CASE 3: Single Cavity-Bulk Viscosity Type Exchange Property
**#############################################################################################################
*SYSTEM
0, 2, 0
*NODE, NSET=SINGLECAVITY_B_VISCOSITY
31, 0., 0., 0.
32, 0.5, 0., 0.
33, 0.5, 0.5, 0.
34, 0., 0.5, 0.
35, 0., 0., 0.5
36, 0.5, 0., 0.5
37, 0.5, 0.5, 0.5
38, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SINGLECAVITY_B_VISCOSITY_REF
399, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SINGLECAVITY_B_VISCOSITY
31, 32, 33, 34, 31
32, 38, 37, 36, 35
33, 35, 36, 32, 31
34, 36, 37, 33, 32
35, 37, 38, 34, 33
36, 38, 35, 31, 34
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SINGLECAVITY_B_VISCOSITY, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SINGLECAVITY_B_VISCOSITY_SURF, TYPE=ELEMENT
SINGLECAVITY_B_VISCOSITY, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SINGLECAVITY_B_VISCOSITY_CAV , REF NODE=SINGLECAVITY_B_VISCOSITY_REF, SURFACE=SINGLECAVITY_B_VISCOSITY_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

**
**Effective area for bulk viscosity type fluid exchange is considered to be 0.02
**
*FLUID EXCHANGE, NAME=SINGLECAVITY_B_VISCOSITY, PROPERTY=B_VISCOSITY, EFFECTIVE AREA=0.02
399
**#############################################################################################################
**      TEST CASE 4: Single Cavity-Fabric Leakage Type Exchange Property
**#############################################################################################################
*SYSTEM
0, 3, 0
*NODE, NSET=SINGLECAVITY_F_LEAKAGE
41, 0., 0., 0.
42, 0.5, 0., 0.
43, 0.5, 0.5, 0.
44, 0., 0.5, 0.
45, 0., 0., 0.5
46, 0.5, 0., 0.5
47, 0.5, 0.5, 0.5
48, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SINGLECAVITY_F_LEAKAGE_REF
499, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SINGLECAVITY_F_LEAKAGE
41, 42, 43, 44, 41
42, 48, 47, 46, 45
43, 45, 46, 42, 41
44, 46, 47, 43, 42
45, 47, 48, 44, 43
46, 48, 45, 41, 44
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SINGLECAVITY_F_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SINGLECAVITY_F_LEAKAGE_SURF, TYPE=ELEMENT
SINGLECAVITY_F_LEAKAGE, SPOS
**
**Create surface on primary cavity to define effective area with fabric leakage type of exchange behavior
**
*SURFACE, NAME=SINGLECAVITY_F_LEAKAGE_EFFECTIVE_A_SURF, TYPE=ELEMENT
41, SPOS
44, SPOS
42, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SINGLECAVITY_F_LEAKAGE_CAV , REF NODE=SINGLECAVITY_F_LEAKAGE_REF, SURFACE=SINGLECAVITY_F_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=SINGLECAVITY_F_LEAKAGE, PROPERTY=F_LEAKAGE, SURFACE=SINGLECAVITY_F_LEAKAGE_EFFECTIVE_A_SURF
499
**#############################################################################################################
**      TEST CASE 5: Single Cavity-Mass Rate Leakage Type Exchange Property
**#############################################################################################################
*SYSTEM
0, 4, 0
*NODE, NSET=SINGLECAVITY_MR_LEAKAGE
51, 0., 0., 0.
52, 0.5, 0., 0.
53, 0.5, 0.5, 0.
54, 0., 0.5, 0.
55, 0., 0., 0.5
56, 0.5, 0., 0.5
57, 0.5, 0.5, 0.5
58, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SINGLECAVITY_MR_LEAKAGE_REF
599, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SINGLECAVITY_MR_LEAKAGE
51, 52, 53, 54, 51
52, 58, 57, 56, 55
53, 55, 56, 52, 51
54, 56, 57, 53, 52
55, 57, 58, 54, 53
56, 58, 55, 51, 54
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SINGLECAVITY_MR_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SINGLECAVITY_MR_LEAKAGE_SURF, TYPE=ELEMENT
SINGLECAVITY_MR_LEAKAGE, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SINGLECAVITY_MR_LEAKAGE_CAV , REF NODE=SINGLECAVITY_MR_LEAKAGE_REF, SURFACE=SINGLECAVITY_MR_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=SINGLECAVITY_MR_LEAKAGE, PROPERTY=MR_LEAKAGE
599
**#############################################################################################################
**      TEST CASE 6: Single Cavity-Volume Rate Leakage Type Exchange Property
**#############################################################################################################
*SYSTEM
0, 5, 0
*NODE, NSET=SINGLECAVITY_VR_LEAKAGE
61, 0., 0., 0.
62, 0.5, 0., 0.
63, 0.5, 0.5, 0.
64, 0., 0.5, 0.
65, 0., 0., 0.5
66, 0.5, 0., 0.5
67, 0.5, 0.5, 0.5
68, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SINGLECAVITY_VR_LEAKAGE_REF
699, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SINGLECAVITY_VR_LEAKAGE
61, 62, 63, 64, 61
62, 68, 67, 66, 65
63, 65, 66, 62, 61
64, 66, 67, 63, 62
65, 67, 68, 64, 63
66, 68, 65, 61, 64
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SINGLECAVITY_VR_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SINGLECAVITY_VR_LEAKAGE_SURF, TYPE=ELEMENT
SINGLECAVITY_VR_LEAKAGE, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SINGLECAVITY_VR_LEAKAGE_CAV , REF NODE=SINGLECAVITY_VR_LEAKAGE_REF, SURFACE=SINGLECAVITY_VR_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=SINGLECAVITY_VR_LEAKAGE, PROPERTY=VR_LEAKAGE
699
**#############################################################################################################
**      TEST CASE 7: Dual Cavity-Mass Flux Type Exchange Property
**#############################################################################################################
*SYSTEM
2, 0, 0
*NODE, NSET=FIRSTCAVITY_M_FLUX
71, 0., 0., 0.
72, 0.5, 0., 0.
73, 0.5, 0.5, 0.
74, 0., 0.5, 0.
75, 0., 0., 0.5
76, 0.5, 0., 0.5
77, 0.5, 0.5, 0.5
78, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=FIRSTCAVITY_M_FLUX_REF
799, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=FIRSTCAVITY_M_FLUX
71, 72, 73, 74, 71
72, 78, 77, 76, 75
73, 75, 76, 72, 71
74, 76, 77, 73, 72
75, 77, 78, 74, 73
76, 78, 75, 71, 74
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
FIRSTCAVITY_M_FLUX, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=FIRSTCAVITY_M_FLUX_SURF, TYPE=ELEMENT
FIRSTCAVITY_M_FLUX, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=FIRSTCAVITY_M_FLUX_CAV , REF NODE=FIRSTCAVITY_M_FLUX_REF, SURFACE=FIRSTCAVITY_M_FLUX_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*SYSTEM
4, 0, 0
*NODE, NSET=SECONDCAVITY_M_FLUX
81, 0., 0., 0.
82, 0.5, 0., 0.
83, 0.5, 0.5, 0.
84, 0., 0.5, 0.
85, 0., 0., 0.5
86, 0.5, 0., 0.5
87, 0.5, 0.5, 0.5
88, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SECONDCAVITY_M_FLUX_REF
899, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SECONDCAVITY_M_FLUX
81, 82, 83, 84, 81
82, 88, 87, 86, 85
83, 85, 86, 82, 81
84, 86, 87, 83, 82
85, 87, 88, 84, 83
86, 88, 85, 81, 84
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SECONDCAVITY_M_FLUX, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SECONDCAVITY_M_FLUX_SURF, TYPE=ELEMENT
SECONDCAVITY_M_FLUX, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SECONDCAVITY_M_FLUX_CAV , REF NODE=SECONDCAVITY_M_FLUX_REF, SURFACE=SECONDCAVITY_M_FLUX_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=DUALCAVITY_M_FLUX, PROPERTY=M_FLUX
799, 899
**#############################################################################################################
**      TEST CASE 8: Dual Cavity-Volume Flux Type Exchange Property
**#############################################################################################################
*SYSTEM
2, 1, 0
*NODE, NSET=FIRSTCAVITY_V_FLUX
91, 0., 0., 0.
92, 0.5, 0., 0.
93, 0.5, 0.5, 0.
94, 0., 0.5, 0.
95, 0., 0., 0.5
96, 0.5, 0., 0.5
97, 0.5, 0.5, 0.5
98, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=FIRSTCAVITY_V_FLUX_REF
999, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=FIRSTCAVITY_V_FLUX
91, 92, 93, 94, 91
92, 98, 97, 96, 95
93, 95, 96, 92, 91
94, 96, 97, 93, 92
95, 97, 98, 94, 93
96, 98, 95, 91, 94
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
FIRSTCAVITY_V_FLUX, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=FIRSTCAVITY_V_FLUX_SURF, TYPE=ELEMENT
FIRSTCAVITY_V_FLUX, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=FIRSTCAVITY_V_FLUX_CAV , REF NODE=FIRSTCAVITY_V_FLUX_REF, SURFACE=FIRSTCAVITY_V_FLUX_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*SYSTEM
4, 1, 0
*NODE, NSET=SECONDCAVITY_V_FLUX
101, 0., 0., 0.
102, 0.5, 0., 0.
103, 0.5, 0.5, 0.
104, 0., 0.5, 0.
105, 0., 0., 0.5
106, 0.5, 0., 0.5
107, 0.5, 0.5, 0.5
108, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SECONDCAVITY_V_FLUX_REF
1099, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SECONDCAVITY_V_FLUX
101, 102, 103, 104, 101
102, 108, 107, 106, 105
103, 105, 106, 102, 101
104, 106, 107, 103, 102
105, 107, 108, 104, 103
106, 108, 105, 101, 104
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SECONDCAVITY_V_FLUX, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SECONDCAVITY_V_FLUX_SURF, TYPE=ELEMENT
SECONDCAVITY_V_FLUX, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SECONDCAVITY_V_FLUX_CAV , REF NODE=SECONDCAVITY_V_FLUX_REF, SURFACE=SECONDCAVITY_V_FLUX_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=DUALCAVITY_V_FLUX, PROPERTY=V_FLUX
999, 1099
**#############################################################################################################
**      TEST CASE 9: Dual Cavity-Bulk Viscosity Type Exchange Property
**#############################################################################################################
*SYSTEM
2, 2, 0
*NODE, NSET=FIRSTCAVITY_B_VISCOSITY
111, 0., 0., 0.
112, 0.5, 0., 0.
113, 0.5, 0.5, 0.
114, 0., 0.5, 0.
115, 0., 0., 0.5
116, 0.5, 0., 0.5
117, 0.5, 0.5, 0.5
118, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=FIRSTCAVITY_B_VISCOSITY_REF
1199, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=FIRSTCAVITY_B_VISCOSITY
111, 112, 113, 114, 111
112, 118, 117, 116, 115
113, 115, 116, 112, 111
114, 116, 117, 113, 112
115, 117, 118, 114, 113
116, 118, 115, 111, 114
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
FIRSTCAVITY_B_VISCOSITY, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=FIRSTCAVITY_B_VISCOSITY_SURF, TYPE=ELEMENT
FIRSTCAVITY_B_VISCOSITY, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=FIRSTCAVITY_B_VISCOSITY_CAV , REF NODE=FIRSTCAVITY_B_VISCOSITY_REF, SURFACE=FIRSTCAVITY_B_VISCOSITY_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*SYSTEM
4, 2, 0
*NODE, NSET=SECONDCAVITY_B_VISCOSITY
121, 0., 0., 0.
122, 0.5, 0., 0.
123, 0.5, 0.5, 0.
124, 0., 0.5, 0.
125, 0., 0., 0.5
126, 0.5, 0., 0.5
127, 0.5, 0.5, 0.5
128, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SECONDCAVITY_B_VISCOSITY_REF
1299, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SECONDCAVITY_B_VISCOSITY
121, 122, 123, 124, 121
122, 128, 127, 126, 125
123, 125, 126, 122, 121
124, 126, 127, 123, 122
125, 127, 128, 124, 123
126, 128, 125, 121, 124
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SECONDCAVITY_B_VISCOSITY, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SECONDCAVITY_B_VISCOSITY_SURF, TYPE=ELEMENT
SECONDCAVITY_B_VISCOSITY, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SECONDCAVITY_B_VISCOSITY_CAV , REF NODE=SECONDCAVITY_B_VISCOSITY_REF, SURFACE=SECONDCAVITY_B_VISCOSITY_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

**
**Effective area for bulk viscosity type fluid exchange is considered to be 0.02
**
*FLUID EXCHANGE, NAME=DUALCAVITY_B_VISCOSITY, PROPERTY=B_VISCOSITY, EFFECTIVE AREA=0.02
1199, 1299
**#############################################################################################################
**      TEST CASE 10: Dual Cavity-Fabric Leakage Type Exchange Property
**#############################################################################################################
*SYSTEM
2, 3, 0
*NODE, NSET=FIRSTCAVITY_F_LEAKAGE
131, 0., 0., 0.
132, 0.5, 0., 0.
133, 0.5, 0.5, 0.
134, 0., 0.5, 0.
135, 0., 0., 0.5
136, 0.5, 0., 0.5
137, 0.5, 0.5, 0.5
138, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=FIRSTCAVITY_F_LEAKAGE_REF
1399, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=FIRSTCAVITY_F_LEAKAGE
131, 132, 133, 134, 131
132, 138, 137, 136, 135
133, 135, 136, 132, 131
134, 136, 137, 133, 132
135, 137, 138, 134, 133
136, 138, 135, 131, 134
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
FIRSTCAVITY_F_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=FIRSTCAVITY_F_LEAKAGE_SURF, TYPE=ELEMENT
FIRSTCAVITY_F_LEAKAGE, SPOS
**
**Create surface on primary cavity to define effective area with fabric leakage type of exchange behavior
**
*SURFACE, NAME=FIRSTCAVITY_F_LEAKAGE_EFFECTIVE_A_SURF, TYPE=ELEMENT
131, SPOS
134, SPOS
132, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=FIRSTCAVITY_F_LEAKAGE_CAV , REF NODE=FIRSTCAVITY_F_LEAKAGE_REF, SURFACE=FIRSTCAVITY_F_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*SYSTEM
4, 3, 0
*NODE, NSET=SECONDCAVITY_F_LEAKAGE
141, 0., 0., 0.
142, 0.5, 0., 0.
143, 0.5, 0.5, 0.
144, 0., 0.5, 0.
145, 0., 0., 0.5
146, 0.5, 0., 0.5
147, 0.5, 0.5, 0.5
148, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SECONDCAVITY_F_LEAKAGE_REF
1499, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SECONDCAVITY_F_LEAKAGE
141, 142, 143, 144, 141
142, 148, 147, 146, 145
143, 145, 146, 142, 141
144, 146, 147, 143, 142
145, 147, 148, 144, 143
146, 148, 145, 141, 144
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SECONDCAVITY_F_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SECONDCAVITY_F_LEAKAGE_SURF, TYPE=ELEMENT
SECONDCAVITY_F_LEAKAGE, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SECONDCAVITY_F_LEAKAGE_CAV , REF NODE=SECONDCAVITY_F_LEAKAGE_REF, SURFACE=SECONDCAVITY_F_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=DUALCAVITY_F_LEAKAGE, PROPERTY=F_LEAKAGE, SURFACE=FIRSTCAVITY_F_LEAKAGE_EFFECTIVE_A_SURF
1399, 1499
**#############################################################################################################
**      TEST CASE 11: Dual Cavity-Mass Rate Leakage Type Exchange Property
**#############################################################################################################
*SYSTEM
2, 4, 0
*NODE, NSET=FIRSTCAVITY_MR_LEAKAGE
151, 0., 0., 0.
152, 0.5, 0., 0.
153, 0.5, 0.5, 0.
154, 0., 0.5, 0.
155, 0., 0., 0.5
156, 0.5, 0., 0.5
157, 0.5, 0.5, 0.5
158, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=FIRSTCAVITY_MR_LEAKAGE_REF
1599, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=FIRSTCAVITY_MR_LEAKAGE
151, 152, 153, 154, 151
152, 158, 157, 156, 155
153, 155, 156, 152, 151
154, 156, 157, 153, 152
155, 157, 158, 154, 153
156, 158, 155, 151, 154
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
FIRSTCAVITY_MR_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=FIRSTCAVITY_MR_LEAKAGE_SURF, TYPE=ELEMENT
FIRSTCAVITY_MR_LEAKAGE, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=FIRSTCAVITY_MR_LEAKAGE_CAV , REF NODE=FIRSTCAVITY_MR_LEAKAGE_REF, SURFACE=FIRSTCAVITY_MR_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*SYSTEM
4, 4, 0
*NODE, NSET=SECONDCAVITY_MR_LEAKAGE
161, 0., 0., 0.
162, 0.5, 0., 0.
163, 0.5, 0.5, 0.
164, 0., 0.5, 0.
165, 0., 0., 0.5
166, 0.5, 0., 0.5
167, 0.5, 0.5, 0.5
168, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SECONDCAVITY_MR_LEAKAGE_REF
1699, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SECONDCAVITY_MR_LEAKAGE
161, 162, 163, 164, 161
162, 168, 167, 166, 165
163, 165, 166, 162, 161
164, 166, 167, 163, 162
165, 167, 168, 164, 163
166, 168, 165, 161, 164
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SECONDCAVITY_MR_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SECONDCAVITY_MR_LEAKAGE_SURF, TYPE=ELEMENT
SECONDCAVITY_MR_LEAKAGE, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SECONDCAVITY_MR_LEAKAGE_CAV , REF NODE=SECONDCAVITY_MR_LEAKAGE_REF, SURFACE=SECONDCAVITY_MR_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=DUALCAVITY_MR_LEAKAGE, PROPERTY=MR_LEAKAGE
1599, 1699
**#############################################################################################################
**      TEST CASE 12: Dual Cavity-Volume Rate Leakage Type Exchange Property
**#############################################################################################################
*SYSTEM
2, 5, 0
*NODE, NSET=FIRSTCAVITY_VR_LEAKAGE
171, 0., 0., 0.
172, 0.5, 0., 0.
173, 0.5, 0.5, 0.
174, 0., 0.5, 0.
175, 0., 0., 0.5
176, 0.5, 0., 0.5
177, 0.5, 0.5, 0.5
178, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=FIRSTCAVITY_VR_LEAKAGE_REF
1799, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=FIRSTCAVITY_VR_LEAKAGE
171, 172, 173, 174, 171
172, 178, 177, 176, 175
173, 175, 176, 172, 171
174, 176, 177, 173, 172
175, 177, 178, 174, 173
176, 178, 175, 171, 174
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
FIRSTCAVITY_VR_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=FIRSTCAVITY_VR_LEAKAGE_SURF, TYPE=ELEMENT
FIRSTCAVITY_VR_LEAKAGE, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=FIRSTCAVITY_VR_LEAKAGE_CAV , REF NODE=FIRSTCAVITY_VR_LEAKAGE_REF, SURFACE=FIRSTCAVITY_VR_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*SYSTEM
4, 5, 0
*NODE, NSET=SECONDCAVITY_VR_LEAKAGE
181, 0., 0., 0.
182, 0.5, 0., 0.
183, 0.5, 0.5, 0.
184, 0., 0.5, 0.
185, 0., 0., 0.5
186, 0.5, 0., 0.5
187, 0.5, 0.5, 0.5
188, 0., 0.5, 0.5
**
** Cavity Ref. node placed at the center of the cubic cavity
**
*NODE, NSET=SECONDCAVITY_VR_LEAKAGE_REF
1899, 0.25, 0.25, 0.25
**
** Elements defined such that all normals point inward
**
*ELEMENT, TYPE=SFM3D4R, ELSET=SECONDCAVITY_VR_LEAKAGE
181, 182, 183, 184, 181
182, 188, 187, 186, 185
183, 185, 186, 182, 181
184, 186, 187, 183, 182
185, 187, 188, 184, 183
186, 188, 185, 181, 184
**
**Fix all DOF's of all the nodes of the cubic structure
**
*BOUNDARY
SECONDCAVITY_VR_LEAKAGE, 1, 3
**
**Create surface to define surface based fluid cavity
**
*SURFACE, NAME=SECONDCAVITY_VR_LEAKAGE_SURF, TYPE=ELEMENT
SECONDCAVITY_VR_LEAKAGE, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=SECONDCAVITY_VR_LEAKAGE_CAV , REF NODE=SECONDCAVITY_VR_LEAKAGE_REF, SURFACE=SECONDCAVITY_VR_LEAKAGE_SURF, BEHAVIOR=LIQUID, AMBIENT PRESSURE=1.0E5

*FLUID EXCHANGE, NAME=DUALCAVITY_VR_LEAKAGE, PROPERTY=VR_LEAKAGE
1799, 1899
**
**Element set for section assignment
**
*ELSET, ELSET=ELSET_ALL
SINGLECAVITY_M_FLUX
SINGLECAVITY_V_FLUX
SINGLECAVITY_B_VISCOSITY
SINGLECAVITY_F_LEAKAGE
SINGLECAVITY_MR_LEAKAGE
SINGLECAVITY_VR_LEAKAGE
FIRSTCAVITY_M_FLUX
SECONDCAVITY_M_FLUX
FIRSTCAVITY_V_FLUX
SECONDCAVITY_V_FLUX
FIRSTCAVITY_B_VISCOSITY
SECONDCAVITY_B_VISCOSITY
FIRSTCAVITY_F_LEAKAGE
SECONDCAVITY_F_LEAKAGE
FIRSTCAVITY_MR_LEAKAGE
SECONDCAVITY_MR_LEAKAGE
FIRSTCAVITY_VR_LEAKAGE
SECONDCAVITY_VR_LEAKAGE
**
**Node set for prescribing initial temperature and pressure initial conditions
**
*NSET, NSET=PRIM_CAV_REF
SINGLECAVITY_M_FLUX_REF
SINGLECAVITY_V_FLUX_REF
SINGLECAVITY_B_VISCOSITY_REF
SINGLECAVITY_F_LEAKAGE_REF
SINGLECAVITY_MR_LEAKAGE_REF
SINGLECAVITY_VR_LEAKAGE_REF
FIRSTCAVITY_M_FLUX_REF
FIRSTCAVITY_V_FLUX_REF
FIRSTCAVITY_B_VISCOSITY_REF
FIRSTCAVITY_F_LEAKAGE_REF
FIRSTCAVITY_MR_LEAKAGE_REF
FIRSTCAVITY_VR_LEAKAGE_REF
**
**Node set for prescribing initial temperature and pressure initial conditions
**
*NSET, NSET=SECOND_CAV_REF
SECONDCAVITY_M_FLUX_REF
SECONDCAVITY_V_FLUX_REF
SECONDCAVITY_B_VISCOSITY_REF
SECONDCAVITY_F_LEAKAGE_REF
SECONDCAVITY_MR_LEAKAGE_REF
SECONDCAVITY_VR_LEAKAGE_REF
**
*NSET, NSET=qa_testPRIM_CAV_REF
PRIM_CAV_REF
*NSET, NSET=qa_testSECOND_CAV_REF
SECOND_CAV_REF
**
**Section assignment
**
*SURFACE SECTION, ELSET=ELSET_ALL
**
**Fluid behavior definition
**
*FLUID BEHAVIOR, NAME=LIQUID
*FLUID DENSITY
1000
*FLUID BULK MODULUS
4.0E6
**
**Initial conditions
**
*INITIAL CONDITIONS, TYPE=TEMPERATURE
PRIM_CAV_REF, 20
SECOND_CAV_REF, 20
*INITIAL CONDITIONS, TYPE=FLUID PRESSURE
PRIM_CAV_REF, 4.0E5
**
**Fluid exchange property definitions
**
*FLUID EXCHANGE PROPERTY, NAME=M_FLUX, TYPE=MASS FLUX
7.00
*FLUID EXCHANGE PROPERTY, NAME=V_FLUX, TYPE=VOLUME FLUX
7.00E-3
**
**Effective area, A=0.02 m^2
**
*FLUID EXCHANGE PROPERTY, NAME=B_VISCOSITY, TYPE=BULK VISCOSITY
400, 200
**
**Effective leakage area of the fabric = 3*Area of one face of the cubic primary cavity = 3*0.25 = 0.75 m^2
**
*FLUID EXCHANGE PROPERTY, NAME=F_LEAKAGE, TYPE=FABRIC LEAKAGE
2.0E-4
**
**Mass flow rate is considered a function of pressure difference and average absolute pressure
**
*FLUID EXCHANGE PROPERTY, NAME=MR_LEAKAGE, TYPE=MASS RATE LEAKAGE
0, 0, 0
2.0, 0.025E5, 1.0E5
5.0, 3.000E5, 1.0E5
3.0, 0.025E5, 5.0E5
7.0, 3.000E5, 5.0E5
**
**Volume flow rate is considered only a function of pressure difference
**
*FLUID EXCHANGE PROPERTY, NAME=VR_LEAKAGE, TYPE=VOLUME RATE LEAKAGE
0, 0
7.00E-3, 4.0E5
**
**Step 1: Activate fluid exchange between single cavity and its environment and between two fluid cavities.
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT USER CONTROL
0.001, 1.0
**
**Activate fluid exchange interaction
**
*FLUID EXCHANGE ACTIVATION
SINGLECAVITY_M_FLUX
SINGLECAVITY_V_FLUX
SINGLECAVITY_B_VISCOSITY
SINGLECAVITY_F_LEAKAGE
SINGLECAVITY_MR_LEAKAGE
SINGLECAVITY_VR_LEAKAGE
DUALCAVITY_M_FLUX
DUALCAVITY_V_FLUX
DUALCAVITY_B_VISCOSITY
DUALCAVITY_F_LEAKAGE
DUALCAVITY_MR_LEAKAGE
DUALCAVITY_VR_LEAKAGE
**
**Time interval for history output request should be changed to smaller value for smooth XY curves.
**
*OUTPUT, HISTORY, TIME INTERVAL=0.25
*NODE OUTPUT, NSET=qa_testPRIM_CAV_REF
PCAV, CMASS, CMFL, CMFLT
*NODE OUTPUT, NSET=qa_testSECOND_CAV_REF
PCAV, CMASS, CMFL, CMFLT
*END STEP