Include Page | ||||
---|---|---|---|---|
|
Include Page | ||||
---|---|---|---|---|
|
Physics Setup
Launch Fluent
To launch FLUENT, double click Setup under Project Schematic.
Check the option "Double Precision". This will make the simulation more accurate, since Fluent will use twice more as many bits in the calculations.
...
You can leave the rest as Default. Select Parallel if desired. Now press OK to launch FLUENT.
The FLUENT window layout is basically divided into 4 zones:
1) The left most column contain all the possible operations that one may want to perform to setup a simulation, like setting the boundary conditions, selecting the fluid used, and which simplifications to make.
We usually go from the first "General", one by one, until "Run Calculation". However in this tutorial we will have to break a little bit this symmetry to specify what is called "Interface Zones" (not (yet) shown in the list of options).
Here are just the main options; the actual "knob tweak" is made right next to it. Note that once you click in one item from the menu, a grey column appear right next to it showing more detailed options.
2) Sub-options for each main option selected from the left most column. Here we actually modify the way Fluent will behave.
3) The graphical window, showing at first launch the mesh. You can navigate on that using the left mouse button to Pan the mesh.
Drawing a box using the middle mouse button from left to right (either top or bottom) will zoom in, while drawing a box from right to left will zoom out.
The right mouse button will not often be used for our relative simple simulations. It displays the "surface group" the "id" and the "name" of the surface containing the point you clicked.
4) A mesh box where Fluent will print its operations, messages, warnings, iterations, results (as requested), etc.
Enough of talking! Let's start the CFD Simulation!
Check Mesh
It is always a good practice to check the mesh, specially especially if you are importing a mesh that wasn't was not created by you. It is not frequent rare to get matching errors, ; particularly when with meshes with that have multiple Cell-Zones like our casemesh.
In order to do thatcheck the mesh, you can click on Mesh > Check under the , which is located under the upper tab with options, next to File.
You can also note that , when Fluent is first launched, the second menu box (described previously) contains the options for "General". There, on the first option you can see under "Mesh", the option for Checkan option to check the mesh is available under "General". You can see this highlighted in the above figure.
Lastly, you can type "mesh/check" and hit Enter in the messages box (number 4 from previous description).
Create mesh interfaces
As briefly discussed before, we need to start by creating the interfaces between the different zones of the mesh.
To account for the rotation of the mesh, we couldn't simply have created cannot simply create one single mesh. Instead, we had to create sub regions , or Zones of meshing.
The mesh created is "non-conformal", i.e. the nodes does do not match across the interface), so we need to tell Fluent that the adjacent cell across that interface share information.
To make the "Interfaces" option appear, first go to "Boundary Conditions". Select the first Zone on the list, "blade_bot_in," and change the type to "Interface". Click OK in the popped window to keep the default name.
Note that the "Mesh Interfaces" option appear appears on the main menu to the left.
Go ahead and set the other interfaces Boundary Conditions. There is no quick way to do that , as Fluent does not allow you to Copy Interfaces Boundary Conditions. So, proceed as before and assign the following Zones to Interface:
...
In our case, each pair of Interface Boundary Condition will be an Interface Zone. So, there will be 4 zones in total.
...
In the first box, "Mesh Interface", write "blade_bot". That's the name of the interface. No select as Under Interface Zone 1, select blade_bot_in. For Interface Zone 2, select blade_bot_outer. Click "Create". The window will not close by it self, itself so go ahead and close it.
Repeat this and create the remaining 3 interface zones.
Mesh Interface | Interface Zone 1 | Interface Zone 2 |
---|---|---|
blade_bot | blade_bot_in | blade_bot_outer |
blade_right | blade_right_in | blade_right_outer |
blade_top | blade_top_in | blade_top_outouter |
hub | hub_inner | hub_outer |
Be careful to select the correct in-outer pairs!!
Note: you don't do not have to close the Mesh Interface window all the time. But be aware that after selecting the first zone, a warning will pop saying that the second zone is already assigned to other interface. Just close the message and select the correct Interface Zone 2After clicking "Create" you can go ahead and give the name for the new Interface and create it.
By the end of this step, you your window should look like this.
Set Solver informations
Now, resuming...
Select the first option under "Setup": "General".
Here we can change key information about the solver, like whether it's steady or transient, planar or axisymmetric, etc.
For our problem, you can leave the default options , i.e., Pressure-Based, Steady and Planar.
Lastly, click in "Units" and change the "angular-velocity" unit to rpm.
Set the Model
Here is where we tell FLUENT all the simplifications for the model it can assume. For instance, here you specify if the model is Inviscid, Laminar or Turbulent, if you should consider the Energy Equation (for supersonic flows), and other options.
For us, we are will be considering turbulence, and we will use the k-epsilon Realizable model.
Highlight "ModelModels" and double click the third item in the list, "Viscous - Laminar". Select "k-epsilon (2 eqn)". Change "k-epsilon Model" to Realizable. Retain the rest as default. Click Ok.
Materials
We will use the default properties for air. Go ahead and check if they are correct.
...
Viscosity: 1.7894e-05 kg/m-s
Cell Zone Conditions
Here we specify to fluent the material of each meshed zone (usually correct by default, unless we create a new material).
Also, we set this if we 're are solving a moving mesh problem or moving frame of reference problem (our case!). We can also set here the centroid of each zone and the angular velocity of that zone here.
Highlight "Cell Zone Conditions". Note that our problem is made out of 5 zones:
- blade_bot: corresponds to the mesh around the bottom blade
- blade_right: similar as above, but for the right blade
- blade_top: againsimilar as above, but for the top blade
- fluid-surface_body: corresponds to the big circular mesh around the main geometry of the turbine.
- inner: corresponds to the zone between the blades , inside the hub.
These are the same zones described previously , when creating the mesh interfaces. We will have to edit parameters for each of them.
fluid-surface_body
Highlight fluid-surface_body and click "Edit...".
No options should be selected here, and the material must be set to air. This is the default. You can verify that and click Ok.
inner
Highlight inner and click "Edit...". Here we will set that the inner portion of the hub is spinning, adding more terms to the equations to account for local acceleration, even though the mesh is not actually moving.
...
From the geometry, the centroid is at the global origin (0,0). Verify that that 's is inputted.
From the problem statement, the turbine is spinning at 40rpm, so go ahead and input 40 to "Rotational Velocity". This is the absolute velocity. Note that "absolute" is selected under "Relative Specification".
The rest you can keep as default. Click Ok.
FIGURE 11
blades
Highlight blade_top and click "Edit...".
...
Go ahead and input that for X and Y of "Rotation-Axis Origin (Relative)".
FIGURE 12
Do the same for the right and bottom blades. Select "Frame Motion" and make it relative to the "inner" cell zone. The needed values for centroid and angular velocity are summarized below.
...
We are now ready to set the other boundary conditions!
Boundary Conditions
Remember that we had already set some Boundary Conditions, before doing the Mesh Interfaces. However we still need to tell FLUENT what is wall, and specify some pressures and velocities.
Velocity at the inlet
First, specify the velocity at the inlet.
Locate and highlight "farfield1". Change its Type to "velocity-inlet".
FIGURE 13
Click "Edit...". In "Velocity Specification Method", change to "Components". Set "X-Velocity (m/s)" to 10.
Under "Turbulence", change the "Specification Method" to "Intensity and Length Scale". Set the "Turbulence Intensity (%)" to 5 and the "Turbulent Length Scale (m)" to 1.FIGURE 14
Note |
---|
It's very important to change the Specification Method for the Turbulence to "Intensity and Length Scale", otherwise your model won't converge. |
Pressure at the outlet
Locate and highlight "farfield2". Change its Type to "pressure-outlet".
A window will pop up. Leave "Gauge Pressure (pascal)" as default (zero). Apply the same turbulence conditions as describe above for Velocity at inlet boundary condition.
FIGURE 15
Note: is this window does not automatically pop up, click "Edit...", next to where you specified the Type of the boundary condition.
Wall
We need to tell FLUENT that the blades of the turbine are walls (i.e., no-slip condition, or no velocity normal or tangential to there), that are rotating together with the mesh around it.
...
Under "Wall Motion" change to "Moving Wall". On the new options that appeared, change the Motion to "Rotational". Now, specify the Rotation-Axis Origin with the same corresponding values as before. X=-0.02m, Y=-0.034641m. Click Ok.
FIGURE 16
Now, let's copy that to the other blades. Click "Copy", near to "Edit...". On the left list, select the boundary condition you just made,"wall_blade_bot". The two othe wall_blade boundary condition appeared on the right column. Go ahead, click on both (the two must be highlighted), then click "Copy".
FIGURE 17
Click Ok on the window that popped.
MAYBE it's to big? Video...?
ADD MATCH CONTROL TO MESH? (see wind blade tutorial)
Import mesh as a wbpz file instead o msh?
To explain a little bit about the MFR method, see Numerical Results part, where I wrote something very short.
Boundary Conditions:
-farfield 1: velocity inlet....
-farfield 2: pressure outlet...
etc
Note |
---|
Under Construction |
Lastly, we need to change the Origin of rotation of the copied wall boundary conditions. To do so, select each of them ("wall_blade_top" and "wall_blade_right"), click "Edit..." and change the "Rotational-Axis Origin" to the respective centroid.
Boundary Zone | Centroid (X,Y) |
---|---|
wall_blade_top | (-0.02, 0.034641) |
wall_blade_right | (0.04, 0) |
Observe that the options selected for the first wall boundary condition were copied to each new one.
Now that we're done with Boundary Conditions, we're almost ready to run the simulation!
Save your project.
Go to Step 5: Numerical Solution
Go to all (ANSYS or FLUENT ) Learning Modules