Pre-Analysis & Start-Up
Pre-Analysis
In the Pre-Analysis step, we'll review the following:
- Mathematical model: We'll look at the governing equations + boundary conditions and the assumptions contained within the mathematical model.
- Numerical solution procedure in ANSYS: We'll briefly overview the solution strategy used by ANSYS and contrast it to the hand calculation approach.
- Hand-calculations of expected results: We'll use an analytical solution of the mathematical model to predict the expected particle velocity at the blade tip. We'll pay close attention to additional assumptions that have to be made in order to obtain an analytical solution.
Mathematical Model
Governing Equations
The governing equations are the continuity and Navier-Stokes equations. These equations are written in a steady rotating frame of reference. This has the advantage of making our simulation not require a moving mesh.
This form of the Navier Stokes equations has additional terms, namely the centripetal acceleration term and the Coriolis term. The equations that we will use looks as follow:
Conservation of mass:
\begin
\frac
+\nabla \cdot \rho \vec
^
_r =0
\end
Conservation of Momentum (Navier-Stokes):
\begin
\nabla \cdot (\rho \vec
^
_r \vec
^
_r)+\rho(2 \vec
^
\times \vec
^
_r+\vec
^
\times \vec
^
\times \vec
^
)=-\nabla p +\nabla \cdot \overline{\overline{\tau}}_r
\end
Where
$vec
^
_r)$
is the relative velocity (the velocity viewed from the moving frame) and
$vec
^
)$
Under Construction
Numerical Solution Procedure in ANSYS
Under Construction
Hand-Calculations of Expected Results
One simple hand-calculation that we can do now before even starting our simulation is to find theoretical wind velocity at the tip. We can then later compare our answer with what we get from our simulation to verify that they agree.
The velocity, v, on the blade should follow the formula
\begin
v=r \times \omega_{}
\end
\begin
\mathbf{m \ddot
+ k x =0}
\end
\begin
m \ddot
+ k x =0
\end
\begin
\boldsymbol{m \ddot
+ k x =0}
\end
Plugging in our angular velocity of -2.22 rad/s and using the blade length of 43.2 meters plus 1 meter to account for the distance from the root to the hub, we get
\begin
v=-2.22\ \mathrm
\ \mathbf{\hat{k}} \times -44.2\ \mathrm
\ \mathbf{\hat{i}}
\end
\begin
v=98.1\ \mathrm
\ \mathbf{\hat{j}}
\end
Start-Up
*Insert video
It will cover:
-Opening ANSYS
-Setting up the project schematic for the geometry and the Fluid flow, naming them.
Under Construction