ACU-T: 2200 Transition Flow over an Airfoil using the SST Transition Models

This tutorial provides the instructions for setting up, solving and viewing results for a steady simulation of transition flow over a S809 airfoil using the SST (Shear Stress Transport k-ω) turbulence model with transition models (Gamma or Gamma-ReTheta). AcuSolve is used to compute the intermittency and predict the point where the boundary layer transitions from the laminar mode to turbulence mode. This tutorial is designed to introduce you to the modeling concepts necessary to perform simulations using the transition models coupled with the SST Turbulence model.

The basic steps in any CFD simulation are shown in ACU-T: 2000 Turbulent Flow in a Mixing Elbow and the steps to set up external turbulent flow are shown in ACU-T: 2100 Turbulent Flow Over an Airfoil Using the SST Turbulence Model. The following additional capabilities of AcuSolve are introduced in this tutorial:
  • Use of the SST turbulence model with the Gamma transition model
  • Use of the SST turbulence model with Gamma-ReTheta transition model
In this tutorial you will do the following:
  • Analyze the problem
  • Start AcuConsole and create a simulation database
  • Set general problem parameters
  • Set solution strategy parameters
  • Import the geometry for the simulation
  • Create a volume group and apply volume parameters
  • Create surface groups and apply volume parameters
  • Set global and local meshing parameters
  • Generate the mesh
  • Run AcuSolve
  • Monitor the solution with AcuProbe
  • Post-processing the nodal output with AcuFieldView

Prerequisites

You should have already run through the introductory tutorials, ACU-T: 2000 Turbulent Flow in a Mixing Elbow and ACU-T: 2100 Turbulent Flow Over an Airfoil Using the SST Turbulence Model. It is assumed that you have some familiarity with AcuConsole, AcuSolve, and AcuFieldView. You will also need access to a licensed version of AcuSolve.

Prior to running through this tutorial, copy AcuConsole_tutorial_inputs.zip from <Altair_installation_directory>\hwcfdsolvers\acusolve\win64\model_files\tutorials\AcuSolve to a local directory. Extract s809_blunt.x_t from AcuConsole_tutorial_inputs.zip.

The color of objects shown in the modeling window in this tutorial and those displayed on your screen may differ. The default color scheme in AcuConsole is "random," in which colors are randomly assigned to groups as they are created. In addition, this tutorial was developed on Windows. If you are running this tutorial on a different operating system, you may notice a slight difference between the images displayed on your screen and the images shown in the tutorial.

Analyze the Problem

An important step in any CFD simulation is to examine the engineering problem and determine the important parameters that need to be provided to AcuSolve. Parameters can be based on geometrical components (such as volumes, inlets, outlets, or walls) and on flow conditions (such as fluid properties, velocity, or whether the flow should be modeled as turbulent or as laminar).

The problem to be addressed in this tutorial is shown schematically in Figure 1. It consists of a cylindrical bounding region containing air that flows past a S809 airfoil profile. The simulation is performed as the two-dimensional problem by including only a single layer of extruded elements in the airfoil span direction. The velocity vector at the far field boundary of the domain is specified to yield an angle of attack of one degree and a Reynolds Number of ~250000. The airfoil chord is 1 m, and standard air material properties are used for the simulation.


Figure 1. S809 Airfoil in the Flow Domain

The diameter of the cylindrical bounding volume for the airfoil is set to 500 times the airfoil chord. This large bounding volume is selected to ensure that the farfield boundaries are sufficiently far from the airfoil to prevent any influence of blockage of the domain on the solution.

The simulation will be carried out by activating the turbulence transition models. The underlying turbulence model used will be the two-equation SST turbulence model. The problem will be solved with the transition models available in AcuSolve – the one-equation Gamma and the two-equation Gamma-ReTheta transition models. As the name suggests, the transition models predict the point where the boundary layer transitions from the laminar mode to turbulence mode. When in the turbulence regime, the SST turbulence model will be used to determine the flow characteristics.

Introduction to Theory

Transition and Transition Models

In the real world, laminar flow and turbulent flow coexist when obstacles are located inside of a fluid flow. Transition flows from laminar to turbulent flow regimes can be found in many industrial applications including turbomachinery, vehicles, drones, and wind turbines.

A transition from laminar flow to turbulent flow occurs due to various external factors such as freestream turbulence level (vertical disturbances), sound waves (acoustic disturbances), temperature fluctuations (entropy disturbances), streamwise pressure gradients, surface roughness, surface curvature, and vibration. Depending on the external factor or factors, the development of various transition (instability) mechanisms cause the flow to trip from laminar to turbulent. Here are three main mechanisms responsible for transition:
  1. Natural Transition: In the laminar flow regime, viscous forces usually damp out the disturbances. However, when the free stream turbulence is low (below 1 percent) and the Reynolds number is higher than the critical Reynolds number, viscous forces destabilize the shear layer, causing the fluid to transition to a turbulent regime via the development of the three initial disturbances (waves, vorticity and vortex breakdown). This mechanism is called the Natural Transition. It is usually very subtle and progresses very slowly, over a long distance scale.
  2. Bypass Transition: When the initial disturbances are high (due to surface roughness or freestream turbulence levels higher than 1%), turbulence are generated without development of the three initial disturbances (T-S waves, spanwise vorticity, and vortex breakdown) observed in the natural transition.
  3. Separated-Inducted Transition: When a laminar flow experiences adverse pressure gradients (such as airfoil suction surfaces and flow over a sphere), the fluid flow detaches from the wall surface. If the disturbances are low, separation can cause the generation of structures found in natural transition. Larger disturbances generate Kelvin-Helmholtz instabilities, where vortices roll up before breaking down into turbulence. The process involved in separation-induced transition depends on the size of the adverse pressure gradient and the presence of additional disturbances (for example, obstacles).
Since transition affects the flow development, lift, drag, heat transfer, and others, it is important to predict the transition phenomena properly. It is known that conventional turbulence models over-predict the wall-shear stress for transitional flows. Thus, transition models can be used to improve the accuracy of CFD solutions when flows encounter transition in the boundary layer.

Transition Models in AcuSolve

AcuSolve has two transition models (Gamma-ReTheta and Gamma) coupled with the Spalart-Allmaras (SA) and the Shear Stress Transport (SST) turbulence models, resulting in total of four choices of the transition modelling, including: SA-Gamma, SA-Gamma-ReTheta, SST-Gamma, and SST-Gamma-ReTheta.
  1. Gamma-ReTheta Transition Model: The Gamma-ReTheta ( γ R e θ ) model is a correlation-based intermittency model that predicts natural, bypass, and separation-induced transition mechanisms. It is based on two transport equations for intermittency ( γ ) and transition momentum thickness Reynolds number ( R e θ ). The intermittency is a measure of the flow regime that varies between zero (laminar) and one (fully turbulent). It is used to turn on the turbulent kinetic energy production in the turbulent kinetic energy equation, while the transition momentum thickness Reynolds number is used as the transition onset criteria. The Gamma-ReTheta model is more suitable for cases when freestream turbulence intensity may be high or adverse pressure gradients are present in the flow. Most internal flows and some external flows fall in this category.
  2. Gamma Transition Model: The Gamma ( γ ) transition model is a one-equation transition model and follows Galilean invariance by modifying the correlations. The Gamma transition model is well suited for the external aerodynamic cases where freestream turbulence intensity is low.

Transition Model Usage Guidelines

  1. Meshing Guidelines
    1. Generate boundary layer meshes on surfaces where transition onset occurs
    2. The first layer height of the boundary layer mesh should be y+ < 5
    3. The stretch ratio of the boundary layer mesh should be in the range between 1.1 and 1.2
    4. The transition between the boundary layer mesh and the freestream mesh should be smooth
    5. The streamwise surface mesh spacing should be small enough to capture separation bubbles
  2. Convergence Monitoring
    1. The residual and solution ratios of the transition simulations are sometimes stagnated (commonly seen when the bubble forms near the leading-edge of an airfoil)
    2. Despite the appearance of residual and/or solution ratio stagnations, the integrated force and moment over the airfoil could be converged to the satisfactory level. Thus, the monitoring of quantities such as lift, drag, thrust, power, torque, etc. is highly suggested.

The purpose of this tutorial is to provide the instructions for setting up the steady transition simulation of the S809 airfoil using the SST transition models. To keep the total computing time within 5 minutes, some of the meshing guidelines described above have not been followed. However, you must follow these guidelines when solving a problem if a high accuracy in results is desired.

This tutorial consists of two parts. In the first part, you will setup and solve the problem, using the Gamma transition model. After successfully having run the problem with the Gamma model, you will modify the database setup to use the Gamma-ReTheta transition model, and generate the solution again.

Define the Simulation Parameters

Start AcuConsole and Create the Simulation Database

In this tutorial, you will begin by creating a database, populating the geometry-independent settings, loading the geometry, creating volume and surface groups, setting group attributes, adding geometry components to groups and assigning mesh controls and boundary conditions to the groups. Next you will generate a mesh and run AcuSolve to solve for the number of time steps specified. Finally, you will review some characteristics using AcuFieldView.

In the next steps you will start AcuConsole and create the database for storage of the simulation settings.

  1. Start AcuConsole from the Windows Start menu by clicking Start > All Programs > Altair <version> > AcuConsole.
  2. Click the File menu, then click New to open the New data base dialog.
    Note: You can also open the New data base dialog by clicking on the toolbar.
  3. Browse to the location that you would like to use as your working directory.
    This directory is where all files related to the simulation will be stored. The AcuConsole database file (.acs) is stored in this directory. Once the mesh and solution are created, additional files and directories will be created within this directory.
  4. Create a new folder named SST_Transition_Models and open this folder.
  5. Enter sst_transition as the File name for the database.
    Note: In order for other applications to be able to read the files written by AcuConsole, the database path and name should not include spaces.
  6. Save the database to create a backup of your settings.

Define Expressions and Variables Using the Variable Manager

In this step, you will use the Variable Manager in AcuConsole to create a list of expressions that will be used during the model setup process.

  1. Click the Variable List icon from the main toolbar:


    Figure 2.

    The Variable Manager opens.

  2. To create new variables in the Variable Manager, click Add.
    To create new variables in the Variable Manager, you will need to click Add, then populate the columns of the table with appropriate inputs. The columns in the table represent the name of the variable, the expression that defines it, the current value, and a description of the expression.
  3. For Name, enter Uinf.
  4. For Expression, enter 4.0.
  5. For Description, enter Free Stream Velocity of Air.
  6. Repeat this process for the remaining variables shown in the table below:
    Name Expression Description
    Uinf 4.0 Free Stream Velocity of Air
    alpha 1.0 Angle of attack (degrees)
    xVelocity :=Uinf * cos(alpha*deg) X component of velocity
    yVelocity :=Uinf * sin(alpha*deg) Y component of velocity
    Once the expressions are entered, the Variable Manager should appear similar to what is shown below:


    Figure 3.
  7. At this point all variables and expressions are defined and you can close the Variable Manager by selecting Close.

Set General Simulation Parameters

In next steps you will set parameters that apply globally to the simulation. To make this simple, the basic settings applicable for any simulation can be filtered using the BAS filter in the Data Tree Manager. This filter enables display of only a small subset of the available items in the data tree and makes navigation of the entries easier.

  1. Click BAS in the Data Tree Manager to switch to basic view in the Data Tree.


    Figure 4.
  2. Double-click the Global Data Tree item to expand it.
  3. Double-click Problem Description to open the Problem Description detail panel.
  4. Enter S809 Airfoil as the Title for this case.
  5. Enter SST Transition as the Sub title for this case.
  6. Ensure that the Analysis type is set to Steady State.
  7. Change the Turbulence equation from Laminar to SST.
  8. Change the Turbulence transition model from None to Gamma.


    Figure 5.

Set Solution Strategy Parameters

  1. Double-click Auto Solution Strategy to open the Auto Solution Strategy detail panel.
  2. Check that the Analysis type is set to Steady State.
  3. Set the Max time steps to 100.
  4. Set the Convergence tolerance to 0.0001.
  5. Enter 0.5 for Relaxation factor.
  6. Ensure that both the Turbulence and Turbulence transition flags are set to On.


    Figure 6.

Set Material Model Parameters

AcuConsole has three pre-defined standard materials, Air, Aluminium, and Water, with standard parameters defined. In the next steps you will check and if needed modify the material characteristics of the predefined "Air" model to match the desired properties for this problem.

  1. Double-click Material Model in the Data Tree to expand it.


    Figure 7.
  2. Double-click Air in the Data Tree to open the Air detail panel.
    The Material type for air is Fluid. Fluid is the default material type for any new material created in AcuConsole.
  3. In the Density tab check the following.
    1. The Type is set to Constant
    2. The Density value is set to 1.225 kg/m3
  4. Click the Viscosity tab.
    The viscosity of air is 1.781 x 10-5 kg/m - sec

The remaining thermal and other material properties are not critical to this simulation. However, you may browse through the tabs to check the complete material specification.

  1. Save the database to create a backup of your settings. This can be achieved with any of the following methods.
    • Click the File menu, then click Save.
    • Click on the toolbar.
    • Click Ctrl+S.
    Note: Changes made in AcuConsole are saved into the database file (.acs) as they are made. A save operation copies the database to a backup file, which can be used to reload the database from that saved state in the event that you do not want to commit future changes.

Import the Geometry and Define the Model

Import Airfoil Geometry

You will import the geometry in the next part of this tutorial. You will need to know the location of s809_blunt.x_t in order to complete these steps. This file contains information about the geometry in Parasolid ASCII format.
  1. Click File > Import.
  2. Browse to the directory containing s809_blunt.x_t.
  3. Change the file name filter to Parasolid File (*.x_t *.xmt *X_T …).
  4. Select s809_blunt.x_t and click Open to open the Import Geometry dialog.


    Figure 8.

    For this tutorial, the default values for the Import Geometry dialog are used to load the geometry. If you have previously used AcuConsole, be sure that any settings that you might have altered are manually changed to match the default values shown in the figure. With the default settings, volumes from the CAD model are added to a default volume group. Surfaces from the CAD model are added to a default surface group. You will work with groups later in this tutorial to create new groups, set flow parameters, add geometric components, and set meshing parameters.

  5. Click OK to complete the geometry import.
  6. Rotate and zoom in the visualization to view the entire model.


    Figure 9.

Apply Volume Parameters

Volume groups are containers used for storing information about a volume region. This information includes solution and meshing parameters applied to the volume and the geometric regions that these settings are applied to.

When a new geometry is imported, by default AcuConsole will place all volumes from a geometry in a single volume container named "default". You should be able to see it in the data tree upon successful import of your model in the last step, under Model > Volumes > default.

Since the model for this tutorial has only a single volume, it will be the only volume in the default volume group when the geometry is imported. Even when there is a single volume in the model, it is advisable to rename the volume for ease of identification in future. In the next steps you will rename the default volume group container, and set the material and other properties for it.

  1. Minimize Global in the Data Tree Manager and expand the Model tree item by clicking .
  2. Right-click on Surfaces and select Display off to turn off the display of surfaces
  3. Expand the Volumes tree item.
  4. Toggle the display of the default volume container by clicking and next to the volume name.
    Note: You may not see any change when toggling the display if Surfaces are being displayed, as surfaces and volumes may overlap.
  5. Rename the default volume group to Fluid.
    1. Right-click default.
    2. Click Rename from the context menu.
    3. Enter Fluid as the new name, and click Enter.
      Note: When an item in the Data Tree is renamed, the change is not saved until you press the Enter key on your keyboard. If you move the input focus away from the item without entering it, your changes will be lost.
  6. Set up the fluid volume element set.
    1. Expand the Fluid volume group in the tree.
    2. Double-click Element Set under fluid to open it in the detail panel.
    3. Check that the Medium for the volume is set to Fluid. If not, click on the drop-down selector next to Medium and select Fluid.
      Note: The default Medium for a volume in AcuConsole is Fluid.
    4. Click on the drop-down selector next to Material model and select Air.


    Figure 10.

Create Surface Groups

Surface groups are containers used for storing information about a surface, including solution and meshing parameters, and the corresponding surface in the geometry that the parameters will apply to.

In the next steps you will define surface groups, assign the appropriate attributes for each group in the problem and add surfaces to the groups.

In the process of setting up a simulation, you need to move into different panels for setting up the boundary conditions, mesh parameters, etc. which can sometimes be cumbersome (especially for models with too many surfaces). To make it easier, less error prone, and for saving time two new dialogs are provided in AcuConsole which you can use to verify and provide the information for all surface or volume entities at once. They are the Volume Manager and Surface Manager. In this section some features of Surface Manager are exploited.

  1. Right-click on Volumes and select Display off to turn off the display of Volumes.
  2. Right-click Surfaces and select Surface Manager.
  3. In the Surface Manager dialog, click New 3 times to create 3 new surface groups.
  4. Turn display off for all surfaces except for default surface.
  5. Rename Surface Names (column 1) for Surface 1 to Surface 3, and set the Simple BC Active and Simple BC Type columns as per the table shown below.


    Figure 11.
  6. Assign the surfaces to the z_pos and z_neg surface groups.
    1. In the z-pos row, click Add to.
    2. Select the planar surface with the maximum z-coordinate as shown in figure below and click Done.
    3. Similarly assign the surface with the minimum z-coordinate to the z_neg surface group.


      Figure 12.
  7. Assign the outer peripheral surface of the domain to the farfield surface group. Use the following figure as the reference for selecting the required surfaces.


    Figure 13.

    When the geometry was loaded into AcuConsole, all geometry surfaces were placed in the default surface group container. This default surface group was renamed to airfoil. In the previous steps, you assigned some surfaces to various other surface groups that you created. At this point, all that is left in the airfoil surface group are the surfaces which make up the airfoil.

  8. Close the Surface Manager dialog.

Assign the Surface Parameters

  1. Set up the airfoil surface parameters:
    1. In the Data Tree, under Surfaces, expand the airfoil surface group.
    2. Double-click Simple Boundary Condition to open the detail panel.
    3. Ensure the Type is set to Wall.
  2. Set up the z_pos and z_neg surface parameters:
    1. In the Data Tree, under Surfaces, expand the z_pos surface group.
    2. Double-click Simple Boundary Condition to open the detail panel.
    3. Ensure the Type is set to Symmetry.
    4. Similarly, check that the Simple Boundary Condition type for z_neg is also set to Symmetry.
  3. Set up the farfield surface parameters:
    1. In the Data Tree, under Surfaces, expand the farfield surface group.
    2. Double-click Simple Boundary Condition to open the detail panel.
    3. Ensure the Type is set to Far Field.
      When using the Far field boundary condition, you need to enter the vector direction of the flow velocity. Based on this direction, AcuSolve determines whether each element face should be acting as an inlet or an outlet. Using this direction vector, in conjunction with the outward facing normal at each element, it applies the appropriate boundary conditions to the faces of the model. For faces where the flow is entering the model, a nodal boundary condition for velocity and turbulence variables is assigned. For faces where the flow is exiting the model, an element boundary condition for pressure is assigned (that is, an outflow condition).

      To facilitate ease in changing the angle of attack for this simulation, you will set the inlet velocity vector at the far field boundary using the variables that were previously defined. Set the values using the following inputs.

    4. For X Velocity, enter :=xVelocity.
    5. For Y Velocity, enter :=yVelocity.
    6. Set Turbulence input type to Auto.
    7. Set the Turbulence intensity type to Auto.


    Figure 14.

    You can see the value of Percent turbulence intensity is set to 0.5. This value is automatically selected by AcuConsole based on parameters like the turbulence model selected, flow type, etc.

Set Initial Conditions

  1. Under Global in the Data Tree, double-click Nodal Initial Condition to open the dialog in the detail panel.
  2. Set the X velocity to :=xVelocity.
  3. Set the Y velocity to :=yVelocity.
    This will set the X and Y velocities equal to the values specified in the Variable Manager.
  4. Ensure all initial condition types in this panel are set to Constant.

For the SST turbulence model, you need to provide the initial value for kinetic energy and eddy frequency. For the Gamma transition model, you need to provide the initial value for Intermittency, or γ . If you have a reasonable estimate of these values, you can enter them directly in the nodal initial condition fields. One option is to use the same values that are assigned at the inlet boundary. In the absence of good estimates for the initial conditions, it is also possible to let AcuSolve perform an automatic initialization of the turbulence and transition variables. By setting these values to zero, AcuSolve will trigger an automatic initialization of these variables.

  1. Set the Kinetic energy, Eddy frequency and Intermittency to 0.0 to trigger the automatic initialization.


    Figure 15.
  2. Save the database to create a backup of your settings.

Assign Mesh Controls

Set Global Meshing Parameters

Now that the flow characteristics have been set for the whole problem, a sufficiently refined mesh has to be generated.

Global mesh attributes are the meshing parameters applied to the model as a whole without reference to a specific geometric volume, surface, edge, or point. Local mesh attributes are used to create mesh generation controls for specific geometry components of the model.

In the next steps you will set the global mesh attributes.

  1. Click MSH in the Data Tree Manager to filter the settings in the Data Tree to show only the controls related to meshing.
  2. Double-click the Global Data Tree item to expand it.
  3. Double-click Global Mesh Attributes to open the Global Mesh Attributes detail panel.
  4. Change the Mesh size type to Absolute.
    1. Click the drop-down control next to Mesh size type.
    2. Click Absolute.
  5. Set Absolute mesh size to 200.0 m.
  6. Set the Mesh growth rate to 1.4.
  7. Set Maximum sweep angle to 5.0 degrees.


    Figure 16.

Set Mesh Process Parameters

For this simulation, the mesher will locally reduce the height of the boundary layer stack but will keep the total number of layers constant. Using this approach, the height of each layer is scaled by a constant factor to reduce the total height of the stack and avoid the creation of the poor quality boundary layer elements.

  1. Under the Global Tree, double-click on Mesh Process Attributes.
  2. Set Shrink mixed flag to On.
  3. Set the 2D Boundary layer adjacent height smoothing ratio to 0.25.
    This parameter controls how smoothly the local boundary layer heights vary from one element to the next after the layers height are adjusted locally to resolve poor quality elements. A low value of this parameter smooths the variation in height over a large distance, while a value closer to 1.0 enforces a more abrupt change in height. Note that there are separate values of this setting for 2D and 3D boundary layers. For this application, you will be creating a 2D mesh and extruding it in the third direction to create the volume. Therefore, the 2D setting will control the behavior of the mesh in this case.


    Figure 17.

Set Surface Mesh Parameters

Surface mesh attributes are applied to a specific surface in the model. It is a type of local meshing parameter used to create targeted mesh controls for one or more specific surfaces.

Setting local mesh attributes, such as surface mesh attributes, is not mandatory. When a local mesh attribute is not found for a component, the global attributes are used as the mesh generation control for that component. If a local mesh attribute is present, it will take precedence over the global setting.

In the next steps you will set the surface meshing attributes.

  1. Expand Model, Surfaces and then Airfoil.
  2. Double click Surface Mesh Attributes to open the detail panel.
  3. Click the Surface Mesh Attributes check box to activate it.
    The detail panel is populated with more options.
  4. Set the Mesh size type to Absolute.
  5. For Absolute mesh size, enter 0.01 m.
  6. Set the Region of influence parameters to On.
    Mesh controls related to influence region from the surface are now visible.

    Region of influence is a size control that allows you to control the size and growth rate of the surface and volume mesh surrounding a surface based on the distance from the surface.

  7. Set the Influence parameters as follows:
    1. Influence type: Simple
    2. Influence size factor: 2.0
    3. Influence distance: 0.1
  8. Set Boundary layer flag to On.
    Mesh controls related to boundary layer meshing will be visible now.
  9. Set the Boundary layer type to Full Control.
  10. For Resolve, select Number of Layers.
  11. Set the remaining parameters as follows:
    1. First element height: 0.0004
    2. Total layer height: 0.04
    3. Growth rate: 1.2
    4. Boundary layer elements type: Tetrahedron


    Figure 18.

Set Edge Mesh Parameters

To create an optimum mesh on the surface of the airfoil, it is necessary to have high levels of refinement near the leading and trailing edges and a large element size near the mid chord. Since our surface mesh size was set to constant to serve as the size that is propagated into the volume for the region of influence refinement, you will use an edge mesh attribute to control the placement of nodes along the airfoil surface. To accomplish this, you will first need to create an edge group that contains the perimeter edges of the airfoil.

  1. Add a new edge by right-clicking Edges > New.
    A new entity, Edge 1, is created under Edges in the data tree.
  2. Rename Edge 1 to Perimeter Edges.
  3. Select the two perimeter edges of the airfoil to add them to this group.
    1. Right-click Perimeter Edges and click Add to.
    2. Select the two perimeter edges of the airfoil shown below.
    3. Click Done.


      Figure 19.
  4. Activate Edge Mesh Attributes and then double-click on it.
  5. Set the Mesh size type to Absolute Expression.
  6. Click Absolute mesh size expression to open a text editor where you can write an expression for the mesh size.
    This expression will be evaluated by the mesh engine to determine the local mesh size that should be used at each location on the edges in the edge set.
  7. Enter the following expression in the text area editor:
    sizeMin = 0.002
    sizeMax = 0.05
    
    xMin = (0.0)
    xMax = (1.0)
    
    xCenter = ((xMin+xMax)/2.0)
    dx = ((xMax-xMin)/2.0)
    
    fac = ((dx - abs( x-xCenter) )/dx )
    
    power = (ln( sizeMax/sizeMin))
    expVal = sizeMin*exp(power*fac)
    
    linVal = sizeMin+(sizeMax-sizeMin)*fac
    
    value = expVal

    This expression takes the min and max surface mesh size (sizeMin and sizeMax) along with the location of the leading and trailing edge (xMin and xMax) and computes a logarithmic expansion of the surface mesh size as a function of distance from the leading and trailing edges. The mesh size at the leading and trailing edge corresponds to sizeMin and the size at the mid chord location corresponds to sizeMax.

  8. Select Ok to close this text editor area and accept the expression.

Define Mesh Extrusion

The present simulation is equivalent to a representation of a 2D cross section of the model. In AcuSolve 2D models are simulated by having just one element across the faces of the cross section. Thus when these faces are set up with a similar boundary condition, it coerces the corresponding nodes across the faces to have same results. In this problem, these faces are the negative and positive z-surfaces. This kind of mesh is achieved in AcuSolve with mesh extrusion process. In the following steps you will define the process of extrusion of the mesh between these surfaces.

  1. Right-click Mesh Extrusions and select New.
  2. Right-click Mesh Extrusion 1, select Rename and enter Z extrusion.
  3. Double-click z extrusion to open the settings dialog.
  4. Verify that Geometry type is set to surface.
  5. For Side 1, from the drop down menu, select z_neg.
  6. For Side 2, from the drop down menu, select z_pos.
  7. Check that the Extrusion type is set to Number of layers.
  8. Set Number of layers to 1.
  9. For Extrusion options, select All tets.
  10. Click OK to accept these settings.


    Figure 20.

Generate the Mesh

In the next steps you will generate the mesh that will be used when computing a solution for the problem.

  1. Click on the toolbar to open the Launch AcuMeshSim dialog.
  2. Leave the default settings and select OK.

    During meshing an AcuTail window opens. Meshing progress is reported in this window. A summary of the meshing process indicates that the mesh has been generated.



    Figure 21.
    Note: The actual number of nodes and elements, and memory usage may vary slightly from machine to machine.
  3. Visualize the mesh in the modeling window. Turn on the display of surfaces and set the display type to solid and wire.
  4. You can rotate and zoom in the model to analyse the various mesh regions.


    Figure 22.

Compute the Solution and Review the Results

Run AcuSolve

In the next steps, you will launch AcuSolve to compute the solution for this case.

  1. Click on the toolbar to open the Launch AcuSolve dialog.
    For this case, the default settings will be used. AcuSolve will run using four processors (if available, higher number of processors may be specified) and AcuConsole will generate AcuSolve input files and will launch AcuSolve. AcuSolve will calculate the steady state solution for this problem.
  2. Select Ok to start the solution process.

    While computing the solution, an AcuTail window opens. Solution progress is reported in this window. A summary of the solution process indicates that the run has been completed.

    The information provided in the summary is based on the number of processors used by AcuSolve. If you use a different number of processors than indicated in this tutorial, the summary for your run may be slightly different than the summary shown.



    Figure 23.
  3. Close the AcuTail window and save the database to create a backup of your settings.

Post-Process with AcuProbe

AcuProbe can be used to monitor various variables over solution time.

  1. Open AcuProbe by clicking on the toolbar.
  2. Double-click Residual Ratio to expand the tree.
  3. Right-click Final and click Plot All.
    This will plot the residuals for all the variables, pressure, velocity, eddy viscosity, and intermittency, in the plot area.
    Note: You might need to click on the toolbar in order to properly display the plot.


    Figure 24.
  4. Right-click Final and click Plot None.
  5. Expand Surface Output > airfoil > Forces and Moments.
  6. Right-click on x-traction and select Plot then right-click on y-traction and select plot Plot.
    Note: You might need to click on the toolbar in order to properly display the plot.


    Figure 25.

    The traction values on the airfoil surface have nearly converged. When using correlation-based transition models, it is always a good practice to examine not only the residuals but also the actual solution quantities of interest for convergence before accepting the solution. In some other cases, it is also possible that the flow field has converged even while the residuals show minor oscillations. The user thus should observe both in tandem before taking an informed decision about the validity of the solution.

  7. You can also save the plots as an image.
    1. From the AcuProbe dialog, click File > Save.
    2. Enter a name for the image and click Save.

Post-Process with AcuFieldView

The tutorial has been written with the assumption that you have become familiar with the AcuFieldView interface and basic operations. In general, it will be helpful to understand the following basics:
  • How to find the data readers in the File menu and open up the desired reader panel for data input.
  • How to find the visualization panels either from the Side toolbar or the Visualization panel menu to create and modify surfaces in AcuFieldView
  • How to move the data around the graphics window using mouse actions to translate, rotate and zoom in to the data.
This tutorial shows you how to work with steady state analysis data.
  1. Click on the AcuConsole toolbar to open the Launch AcuFieldView dialog.
  2. Click OK to start AcuConsole.
    You will see that the pressure contours have already been displayed on all the boundary surfaces with mesh.


    Figure 26.

Set Up AcuFieldView

  1. Click Viewer Options.


    Figure 27.
  2. In the Viewer Options dialog:
    1. Deselect Perspective to turn off the perspective view.
    2. Click Axis Markers to disable the axis markers.
    3. Click Close.
  3. On the toolbar, click the Colormap icon .
  4. In the Scalar Colormap Specification dialog, click Background and select White.
  5. Close the Scalar Colormap Specification dialog.
  6. Click the Toggle Outline icon on the toolbar to turn off the outline display.
    Your display should look similar to Figure 2.


    Figure 28.

Coordinate the Surface Showing Turbulence Viscosity on the Mid Coordinate Surface

  1. From the Surface tab in the Boundary Surface dialog, turn off the Visibility for the active boundary surfaces.
  2. Click to open the Coordinate Surface dialog.
  3. Click Create to create a new coordinate surface.
  4. Select Z as the Coordinate Plane.
    This will create a coordinate surface at the mid – Z coordinate plane. The coordinate surface created is the mid plane between the two side surfaces in the model.
  5. Change the Coloring to Scalar.
  6. Change the Display Type to Smooth.
  7. For Scalar Function, select kinetic_energy.
  8. Close the Coordinate Surface dialog.
  9. From the Defined Views menu, select +Z as the viewing direction.
    Your view should be similar to Figure 1.


    Figure 29.

    You can clearly see the turbulent flow developing at about halfway through the chord of the airfoil. Before the onset of turbulence, the boundary layer on the airfoil surface is laminar.

  10. Open the Coordinate Surface dialog, click Select next to Scalar Function and select Intermittency.


    Figure 30.

    A closer look at the contour plot of intermittency will show that the value of intermittency transitions to one at the point where the flow transitions from laminar to turbulent. In the region where the flow is laminar, from the leading edge of the airfoil to the halfway, the intermittency is zero.

Post-Process to Calculate Flow Coefficients

AcuSolve is shipped with a number of utility scripts to facilitate the pre and post-processing of a problem solved using the solver. You will be introduced to two of these scripts, AcuLiftDrag and AcuGetCpCf in this section, and their usage. These two scripts are focused on aerodynamic simulations as the ones solved in this tutorial.

Run AcuLiftDrag

AcuLiftDrag is a utility script used to calculate the lift and drag coefficients for an airfoil. The following options are required as input to run this script.

Option Name Description
-pb problem name
-dir working directory. The default working directory is ACUSIM.DIR, which contains all the solution files for the problem.
-run run number (0, use last)
-osis comma seperated list of surface outputs to use for force computation (all for all)
-angle_of_attack / aoa angle of attack of the wing (in degrees)
-reference_velocity / ref_vel reference velocity
-wing_chord / chord chord of the wing
-wing_span / span span of the wing
-density / rho denisty of the air
-drag_direction / drag_dir direction of drag at zero angle of attack: x, y, z. The default drag direction is x.
-lift_direction / lift_dir direction of lift at angle of attack: x, y, z. The default drag direction is y.
-x_axis_vars / x_vars x axis variable for the exported data: step, time, both. The default is step.

In the Analyze the Problem section, it was described that the simulation is performed as 2D by including only a single layer of extruded elements in the airfoil span direction. When solving a problem in such a way, the span of the airfoil should be set equal to the thickness of the domain in the extrusion direction. When solving a 3D problem, the actual span of the airfoil should be used.

To execute the AcuLiftDrag script for this position, follow the steps below:

  1. Start AcuSolve Command Prompt from the Start menu by clicking Start > All Programs > Altair HyperWorks <version> > AcuSolve > AcuSolve Cmd Prompt.
  2. Change the directory to the present problem directory using the 'cd' command.
  3. Enter the following command at the prompt:
    acuLiftDrag -osis airfoil -aoa 1 -ref_vel 4 -chord 1 -rho 1.225 
    -span 50  
    
    The output of the command should look like the image below:


    Figure 31.

    The value C_l corresponds to the lift coefficient, and C_d corresponds to the drag coefficient. The final keyword refers to the respective coefficient values at the last time step of the simulation. The rest of the values provide the basic statistics for these coefficients over all the time steps of the simulation. These statistics are more meaningful if the simulation is transient.

    The script also creates the file Turbulent_Airfoil_SST.liftDrag.dat in the problem directory. The file contains the lift and drag data for all the available time steps in a determined tabular arrangement. The first column is the time step, the second column is the lift coefficient, and the third column is the drag coefficient.

Run AcuGetCpCf

AcuGetCpCf is another utility script used to calculate the pressure coefficient (Cp) for the airfoil. The following options are required as input to run this script.

Option Name Description
-pb problem name
-dir working directory. The default working directory is ACUSIM.DIR, which contains all the solution files for the problem.
-run run number (0, use last)
-ts timestep to extract (0=use last)
-osis comma seperated list of surface outputs to use for force computation (all for all)
-type type of data to generate: cp, cf
-point_type method used to generate the points for the curve: file, auto
-points_file / pts comma seperated list of files to read points from (used with type=file)
-radial_locations / rad_locs list of radial locations
-z_direction / z_dir direction of iso-cut: x, y, z. Default is z.
-reference_pressure / ref_pres reference pressure
-density / rho density of the air
-reference_velocity / ref_vel reference velocity
-wing_chord / chord chord of the wing
-normalize_chord / nc normalize the local chord to unity. Default is True.
-chord_scale_fac / csf chord scale factor
-cp_scale_fac / cpsf c_p / c_f scale factor
-pressure_type / ptype pressure field to use: standard, time_ave, running_ave. Default is standard.

To execute the AcuGetCpCf script for this problem, follow the steps below:

  1. Start AcuSolve Command Prompt from the Start menu by clicking Start > All Programs > Altair HyperWorks <version> > AcuSolve > AcuSolve Cmd Prompt.
  2. Change the directory to the present problem directory using the 'cd' command.
  3. Enter the following command at the prompt:
    acuGetCpCf -osis airfoil -type cp -ref_vel 4 -rho 
    1.225 -no_nc

    The output of the command should look like the image below:



    Figure 32.

    The script prints the minimum and maximum vlaues for the presure coefficient. It also creates a file, cp.1.dat, in the problem directory. The file contains the pressure coefficient data along the chord of the airfoil. The first column is the x-coordinate along the chord, and the second column is the pressure coefficient. You can use an external plotting utility to the plot the data. The resulting plot is shown below.



    Figure 33.

Set Up the Gamma-ReTheta Transition Model

At this stage you have successfully setup and ran the S809 airfoil problem with the Gamma transition model and SST turbulence model. In this part of the tutorial you will modify your open database to setup the problem so as to use the Gamma-ReTheta transition model.

  1. Close the open AcuFieldView window and return to the open AcuConsole window.
  2. Save the database to retain the setup for the Gamma transition model.
  3. Create a new directory within your existing working directory, or at any other location of your choice, and name it SST_Gamma_ReTheta.
  4. Click File > Save As
  5. Navigate into the SST_Gamma_Re_Theta directory. Enter sst_transition_gammaretheta as the File name for the database, or choose any name of your preference.
  6. Save the database to create a backup of your settings.

Update General Simulation Parameters

  1. Click BAS in the Data Tree Manager to switch to basic view in the Data Tree.


    Figure 34.
  2. Double-click the Global Data Tree item to expand it.
  3. Double-click Problem Description to open the Problem Description detail panel.
  4. Enter SST Transition Gamma ReTheta as the new Sub title for this case.
  5. Change the Turbulence transition model from Gamma to Gamma-ReTheta.


    Figure 35.

Update the Nodal Initial Conditions

The Gamma-ReTheta transition model is a two equation model and introduces a new variable, transition ReTheta, or R e θ t . Like other variables, an initial value for this variable also needs to be provided. As before, you will set it to zero to trigger the automatic initialization of R e θ t by AcuSolve.

  1. Double-click Nodal Initial Condition in the data tree to open the detail panel.
  2. Set Transition re theta to 0.0 to trigger the automatic initialization.


    Figure 36.

    The remaining settings in the setup need not be modified. You can now launch AcuSolve to get the solution of the S809 airfoil problem with the Gamma-ReTheta transition case. Follow the same steps as in the previous case to post-process the results.



    Figure 37.








    Results from running AcuLiftDrag and AcuGetCpCf on the Gamma-Re Theta database are shown below:



    Figure 38.


    Figure 39.


    Figure 40.

Summary

In this AcuSolve tutorial, you successfully set up and solved a turbulence transition problem. The underlying turbulence model employed was the SST model. The problem simulated a S809 wind turbine airfoil in an external flow field. You started the tutorial by creating a database in AcuConsole, importing and meshing the geometry, and setting up the simulation parameters. The database was initially set up with the one-equation Gamma transition model. Once the case was setup, the solution was generated with AcuSolve. Results were post-processed in AcuProbe and AcuFieldView. In AcuFieldView you observed the inter-relation between onset of turbulence viscosity and intermittency. After successfully getting a solution for the Gamma transition model, you modified the database to use the two-equation Gamma-ReTheta as the transition model.