Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

The rest you can keep as default. Click Ok.

FIGURE 11Image Added

blades

Highlight blade_top and click "Edit...". 

...

Go ahead and input that for X and Y of "Rotation-Axis Origin (Relative)".

FIGURE 12Image Added

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.

...

Locate and highlight "farfield1". Change its Type to "velocity-inlet".

FIGURE 13Image Added

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 14Image Added

Pressure at the 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 15Image Added

Note: is this window does not automatically pop up, click "Edit...", next to where you specified the Type of the boundary condition.

...

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 16Image Added

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 17Image Added

Click Ok on the window that popped.

...

For this problem, just change the first option, "Scheme" to "Coupled". Leave the rest as default.

FIGURE 18Image Added

Monitors

Here is where we define parameters to help checking for convergence. Residuals, the first one already embedded, is a measure of the linearization error. We define some threshold value and when the Residuals are below that value, the solution has converged.

...

The right most box containing the 0.001 number is the threshold we just discussed. There are 5 of them, continuity, x-velocity, y-velocity, k and epsilon (you might have to scroll down to see the last one). Change all values to "1e-6". (Note, "Ctrl C" does not work here...). Click Ok.

FIGURE 19Image Added

Now, let's create a monitor for Moment.

Click "Create > Moment...".

FIGURE 20Image Added

Under "Options", check the boxes to "Print to Console" and "Plot". Under "Wall Zones", select "wall_blade_right". Change the "Moment Center" to X=0.04m and Y=0. Click Ok.

FIGURE 21Image Added

Similarly, create a similar monitor for "wall_blade_top". Remember to change the Moment Center to X=-0.02m and Y=0.034641m.

...

Highlight  "Run Calculation". Change the Number of Iterations to about 400 and click Calculate600 and click Calculate.

Note that on the first iterations, a message will be printed saying "Turbulent viscosity limited to viscosity ratio of 1.000000e+05 in number cells". Don't worry about this, as it will stop happening after 70 iterations or so.

The solution will converge in about 470 iterations.

Image Added

You can change what is displayed by clicking on the drop-down menu right above the plot (upper-left corner). Select between Residuals or one of the two moment monitors created. 

Note: if you want to display them side by side as the above picture, click on the "Arrange graphics window layout" and select the desired layout.

Image Added

Save your project (if a window pop up, you can leave the default option selected and click Ok).

 

 

 

MAYBE it's to big? Video...? (sad)

...