CFD_WIND_TUNNEL

CFD

Beta command

This command is in the beta stage and the format may change over time.

*CFD_WIND_TUNNEL
"Optional title"
fid${}_v$

Parameter definition

Variable
Description
fid${}_v$
FUNCTION defining the inlet air flow velocity

Description

This command is used to convert the CFD domain into a virtual wind tunnel. The ambient air conditions are assumed be at pressure $1 \mathrm{atm}$ and temperature $293 \mathrm{K}$. The flow inlet is at $x_0$ and the outlet is at $x_1$ (see CFD_DOMAIN).

Note that this command is only useful at supersonic speeds, as boundary layer and viscous drag effects are not accounted for.

Example

Retardation of an irregular object

In this example we start with the air inlet flow velocity $v_{inlet}(0) = v_0$. The drag force $F(t)$ acting on the object results in a virtual retardation $a(t) = -F(t)/m$, where $m$ is object mass. Instead of actually changing the velocity of the object, this retardation if used to dynamically adjust the air speed at the inlet:

$\displaystyle{\frac{\mathrm{d} v_{inlet}}{\mathrm{d} t} = a(t) = -\frac{F(t)}{m}}$

In this example the inlet velocity is defined through the FUNCTION with ID=12. The function uses the flag derivative=1, which means that the functional expression defines its time derivative and not its current value.

*UNIT_SYSTEM SI *PARAMETER m = 0.006586, "object mass" v0 = 1000.0, "initial velocity" tend = 0.01, "termination time" # # --- TIME --- # *TIME [%tend] # # --- MESH --- # *INCLUDE mesh.k # # --- MATERIAL --- # *MAT_RIGID 1, 60.0 # # --- PART --- # *PART "object" 1, 1 # # --- CONSTRAIN OBJECT IN X-DIRECTION --- # *BC_MOTION "object" 1 P, 1, X # # --- INITIAL OBJECT SPIN --- # *INITIAL_VELOCITY P, 1, 0, 0, 0, 10000, 5000, 3000 # # --- CFD --- # *CFD_DOMAIN 1 P, 1, 0.005, 1, 1 -0.3, -0.3, -0.3, 0.5, 0.3, 0.3 0, 0, 0, 0, 0, 0 *CFD_WIND_TUNNEL 12 *FUNCTION 12, 1, [%v0] -fxr(1)/%m *END