Skip to content

Commit 5e0478d

Browse files
authored
Merge pull request #47 from Nicola-Fonzi/master
Typos in unsteady fsi
2 parents efb0ee1 + 068f803 commit 5e0478d

File tree

1 file changed

+15
-9
lines changed

1 file changed

+15
-9
lines changed

_tutorials/multiphysics/Unsteady_FSI_Python/Dynamic_FSI_Python.md

Lines changed: 15 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -46,25 +46,31 @@ Here, the difference is due to the fact that the simulation is unsteady. Thus, t
4646
The aerodynamic model is based on the compressible Reynolds-averaged Navier-Stokes equations. A central JST scheme is used for the convective fluxes, and a weighted least square
4747
scheme is used for the gradients. The turbulence model is the SST and a CFL number of 20, for the psuedo time step, is used.
4848

49-
Different Mach numbers will be considered, namely $$M=[0.1, 0.2, 0.3, 0.357, 0.364]$$. The Reynolds number is fixed at 4 millions, and the temperature is equal to 273K.
49+
Different Mach numbers will be considered, namely $$Ma=[0.1, 0.2, 0.3, 0.357, 0.364]$$. The Reynolds number is fixed at 4 millions, and the temperature is equal to 273K.
5050

5151
The structural model is made by a single point, positioned at the rotation axis, with two degrees of freedom, pitch and plunge.
5252
Inertia and mass of the airfoil are concentrated at the center of mass of the profile, at a certain distance from the rotation axis. The equations of motions are available
5353
analytically and read:
5454

55-
$$m\ddot{h} + S\ddot{\alpha} + C_{h}\dot{h} + K_{h}h = -L$$
56-
$$S\ddot{h} + I\ddot{\alpha} + C_{\alpha}\dot{\alpha} + K_{\alpha}\alpha = M$$
55+
$$
56+
\begin{cases}
57+
\begin{array}{ll}
58+
m\ddot{h} + S\ddot{\alpha} + C_{h}\dot{h} + K_{h}h = -L \\
59+
S\ddot{h} + I\ddot{\alpha} + C_{\alpha}\dot{\alpha} + K_{\alpha}\alpha = M \\
60+
\end{array}
61+
\end{cases}
62+
$$
5763

5864
Where $$m$$ is the mass of the airfoil, $$I$$ the inertia around the center of mass, $$S$$ the static moment of inertia at the rotation axis, $$C$$ and $$K$$ the dampings and stiffnesses respectively. $$L$$ and $$M$$ are the lift and pitching up moment.
5965

6066
These equations are usually adimensionalised to obtain results independent from the free-stream density of the flow.
6167
Indeed, we can define the following parameters:
6268

63-
$$\Csi=\frac{S}{mb}$$, $$r_{\alpha}^2=\frac{I_f}{mb^2}$$, $$\bar{\omega}=\frac{\omega_h}{\omega_{\alpha}}$$, $$\mu=\frac{m}{\pi \rho_{\inf} b^2}$$
69+
$$\chi=\frac{S}{mb}$$, $$r_{\alpha}^2=\frac{I_f}{mb^2}$$, $$\bar{\omega}=\frac{\omega_h}{\omega_{\alpha}}$$, $$\mu=\frac{m}{\pi \rho_{\infty} b^2}$$
6470

65-
Where $$b$$ is the semi chord of the airfoil, $$\omega_h = \sqrt{\frac{K_h}{m}}$$ $$\omega_{\alpha} = \sqrt{\frac{K_{\alpha}}{I_f}}$$. If we fix them, the structure will behave always the same regardless of $$\rho_{\inf}$$.
71+
Where $$b$$ is the semi chord of the airfoil, $$\omega_h = \sqrt{\frac{K_h}{m}}$$ $$\omega_{\alpha} = \sqrt{\frac{K_{\alpha}}{I_f}}$$. If we fix them, the structure will behave always the same regardless of $$\rho_{\infty}$$.
6672

67-
In this context $$\Csi=0.25$$, $$r_{\alpha}=0.5$$, $$\omega_{\alpha} = 45 rad/s$$ ,$$\bar{\omega}=0.3185$$ and $$\mu=100$$.
73+
In this context $$\chi=0.25$$, $$r_{\alpha}=0.5$$, $$\omega_{\alpha} = 45 rad/s$$, $$\bar{\omega}=0.3185$$ and $$\mu=100$$.
6874

6975
Note that, as we will vary the Mach number, the density will also change accordingly. Thus, with given nondimensional parameters, the inertias and stiffnesses must be
7076
varied accordingly.
@@ -78,8 +84,8 @@ integrate a set of ODEs for the modes of the structure.
7884
To perform the preprocessing step, in the case control section of the Nastran model (i.e. at the very beginning of
7985
the file), the following lines must be added:
8086

81-
ECHO = SORT
82-
DISPLACEMENT(PRINT,PUNCH)=ALL
87+
* __ECHO = SORT__
88+
* __DISPLACEMENT(PRINT,PUNCH)=ALL__
8389

8490
A real egeinvalue analysis will then be performed.
8591
This will produce, in the f06 file, an equivalent, ordered, model that will
@@ -187,7 +193,7 @@ The solver can work in two ways:
187193

188194
2. It can integrate in time the modal equations of motion to study the linearised structural deformations when the body is surrounded by the flow
189195

190-
Available keyword for the config file:
196+
Available keywords for the config file:
191197

192198
* __NMODES__ (int): number of modes to use in the analysis. If n modes are available in the punch file, but only the first m<n are required, set this to m
193199

0 commit comments

Comments
 (0)