/SPH/INOUT
Block Format Keyword Describes the SPH inlet/outlet conditions.
Format
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

/SPH/INOUT/condition_ID  
condition_name  
Ityp  part_ID  surf_ID  Dist  node_ID_{1}  node_ID_{2}  node_ID_{3}  F_{cut} 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

X_{M}  Y_{M}  Z_{M}  
X_{M1}  Y_{M1}  Z_{M1}  
X_{M2}  Y_{M2}  Z_{M2} 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

fct_ID_{r}  Fscale_{r}  fct_ID_{E}  Fscale_{E}  
fct_ID_{Vn} 
Definition
Field  Contents  SI Unit Example 

condition_ID  Inlet/Outlet condition
identifier. (Integer, maximum 10 digits) 

condition_name  Inlet/Outlet condition
name. (Character, maximum 100 characters) 

Ityp  Condition type. 22
24
(Integer) 

part_ID  Part identifier is used in
order to define the SPH particles concerned by the
condition. (Integer) 

surf_ID  Surface
identifier. (Integer) 

Dist  Distance from the surface
for particle control. (Real) 
$\left[\text{m}\right]$ 
node_ID_{1}  Optional ID of node 1
(Itype ≠ 1). (Integer) 

node_ID_{2}  Optional ID of node 2
(Itype ≠ 1). (Integer) 

node_ID_{3}  Optional ID of node 3
(Itype ≠ 1). (Integer) 

Fcut  Optional Cutoff frequency
(Itype ≠ 1). Default = 0 (Real) 

X_{M}  Optional X coordinate of M
(Itype ≠ 1). (Real) 

Y_{M}  Optional Y coordinate of M
(Itype ≠ 1). (Real) 

Z_{M}  Optional Z coordinate of M
(Itype ≠ 1). (Real) 

X_{M1}  Optional X coordinate of
M1 (Itype ≠ 1). (Real) 

Y_{M1}  Optional Y coordinate of
M1 (Itype ≠ 1). (Real) 

Z_{M1}  Optional Z coordinate of
M1 (Itype ≠ 1). (Real) 

X_{M2}  Optional X coordinate of
M2 (Itype ≠ 1). (Real) 

Y_{M2}  Optional Y coordinate of
M2 (Itype ≠ 1). (Real) 

Z_{M2}  Optional Z coordinate of
M2 (Itype ≠ 1). (Real) 

fct_ID_{r}  Function
f_{r}(t)
identifier for density. (Integer) 

Fscale_{r}  Density scale
factor. Default = 0.0 (Real) 
$\left[\frac{\text{kg}}{{\text{m}}^{\text{3}}}\right]$ 
fct_ID_{E}  Function
f_{E}(t)
identifier for energy. (Integer) 

Fscale_{E}  Energy per volume unit
scale factor. Default = 0.0 (Real) 
$\left[\frac{\text{J}}{{\text{m}}^{\text{3}}}\right]$ 
fct_ID_{Vn}  Function
f_{Vn}(t)
identifier for velocity in normal direction. (Integer) 

fct_ID_{P}  Function
f_{P}(t)
identifier for pressure. (Integer) 

Fscale_{P}  Pressure scale
factor. Default = 1.0 (Real) 
$\left[\text{Pa}\right]$ 
${l}_{c}$  Characteristic length.
21 (Real) 
$\left[\text{m}\right]$ 
Comments
 The surface segments must be orientated so that their normal vectors point towards the interior of the domain.
 The surface must be fixed.
 In case of an inlet condition, the condition enters particles belonging to its related part, so long as inactive particles are available for this part. The behavior of the particles belonging to the part which is related to the condition is set with respect to the condition characteristics for all particles lying on the positive side of the surface, within the distance Dist from the inlet surface.
 In case of an outlet condition, the behavior of the particles belonging to the part which is related to the condition is set with respect to the condition characteristics for all particles lying on the negative side of the surface, within the distance Dist from the outlet surface. Such a particle is deactivated if it does not interact with any nonoutgoing particle.
 A particle deactivated by an outlet condition can be reused by an inlet condition acting on the same part for incoming.
 If using outlets, order = 1 is recommended in the relative SPH property.
 In case of an outlet, the initial
net must be provided up to the distance of 2h down to the outlet surface (where h is
the smoothing length into the relative property).In case of inlet or outlet, the distance must be large enough, in order to control incoming or outgoing particles within at least a distance of 2h.
 The domains defined by two inlet/outlet surfaces and distances must not overlap.
 It is recommended for both, inlets and outlets, particles to be initially defined and controlled within more than twice the smoothing length of the particles.
 Inlet/outlet conditions option is allowed for SPMD parallel version. Parallel Arithmetic (same numerical results obtained whatever the number of processors) is not guaranteed for inlet conditions.
 Each incoming particle belonging to
the part which is related to the condition gets the same mass
m_{p} (defined into the
geometrical property which is attached to the part).A particle belonging to this part is entered at the center of a surface segment each time t such that:
(1) $${m}_{p}\ge {\displaystyle \underset{{t}_{last}}{\overset{t}{\int}}\rho (t)\cdot {S}_{i}\cdot v(t)dt}$$Where, S_{i}
 Area of the segment
 $\rho (t)$ and $v(t)$
 The density and velocity of the incoming matier (Lines 2 and 3)
 t_{last}
 Time at last incoming through this segment
It is recommended to use a regular surface mesh.
 If no inactive particle belonging to this part is available for incoming, the program stops and you should provide a larger set of inactive particles for this part.
 If a particle belonging to the part which is related to the condition lies on the positive side of the surface within the Dist, its velocity is set with respect to the data given at Line 5.
 If
fct_ID_{r} = 0, density of the
incoming particles is set to:
(2) $${\rho}_{a}=Fscal{e}_{r}$$else(3) $${\rho}_{a}=Fscal{e}_{{}_{r}}\cdot {\mathrm{f}}_{r}(t)$$  If fct_ID_{E} = 0, energy per volume unit of the incoming particles is set to ${E}_{a}=Fscal{e}_{E}$ , else ${E}_{a}=Fscal{e}_{E}\cdot {\mathrm{f}}_{E}(t)$ .
 If a particle belonging to the part which is related to the condition lies on the negative side of the surface within the Dist, its internal pressure is set with respect to the data given at Line 6.
 If the particle does not interact with any nonoutgoing particle, the particle is deactivated.
 If fct_ID_{P} = 0, internal pressure of the outgoing particles is set to the internal pressure of the closest particle lying above the outlet surface, else it is set to $Fscal{e}_{P}\cdot {\mathrm{f}}_{P}(t)$ .
 If a particle belonging to the part
which is related to the condition lies on the negative side of the surface within
the Dist, its internal pressure is set with respect to the equation:
(4) $$\frac{\partial P}{\partial t}=\rho c\frac{\partial {V}_{n}}{\partial t}+c\frac{({P}_{\infty}P)}{2{l}_{c}}$$  If fct_ID_{P} = 0, pressure in the far field ${P}_{\infty}$ is set to Fscale_{P}, else it is set to Fscale_{P} f_{P}(t).

${l}_{c}$
is the characteristic length, it allows to compute
cutoff frequency
${f}_{c}$
as:
(5) $${f}_{c}=\frac{c}{4\cdot \pi \cdot {l}_{c}}$$  If
Itype = 2, 3 or 4, the surface can be defined by a meshed
surface (surf_ID > 0), by 3 nodes (node_ID_{1} > 0, node_ID_{2} > 0 and node_ID_{3} > 0) or by the coordinates of 3 nodes
(M, M1 and M2).
 If Itype = 2, 3 or 4 and if the surface is defined by 3 coordinates, then the surface will be fixed. If the surface is defined by a surface ID or by 3 nodes, the surface will move according to the displacement of the shell elements or nodes.
 If Itype = 2, 3 or 4, a computation of the total mass crossing the surface is automatically performed and can be plotted using /TH/SPH_FLOW.