14.7. Blade Film Cooling for Gas Turbines

The cooling of the turbine blades within a gas turbine engine is usually achieved via the film-cooling approach. This is where cool gas is injected into the boundary layer of the flow over the blade from lots of small holes on its surface. This presents a problem for CFD modeling, as it can be both difficult and time-consuming to create and mesh a geometry which adequately resolves all the holes.

Figure 14.19: Schematic Diagrams of a Cooled Turbine Blade

Schematic Diagrams of a Cooled Turbine Blade

Within Ansys Fluent, you can easily create models for Blade Film Cooling (BFC) using the Boundary Interface feature. This involves taking a standard blade model setup, and creating virtual holes on an existing boundary mesh. The holes can then be treated as an inlet or outlet, and with the associated boundary condition settings, the simulation can be solved in the usual way. The holes can also be moved, or have their size and shape changed, without the need for costly remeshing.

14.7.1. Specifying the Virtual Hole Geometry

Steps for creating a virtual Boundary Interface to perform Blade Film Cooling are outlined below:

  1. Enable Turbo Models for the domain.

     Domain Turbomachinery Turbo Models

    Once enabled, Boundary Interfaces appears in the outline view.

  2. Import hole center locations from a profile in .csv format. See Profiles for details.

    Each data row in the profile represents an individual hole (x,y,z coordinates) and additional columns can be used to provide flow data for each hole. The Interpolation Method should remain the default constant.

    The Reference Frame specified for the profile controls how the hole positions behave during the simulation. For steady cases, the frame should be set to global. For transient cases, the chosen reference frame depends on the motion of the cell zone to which the holes are connected.

    • If the cell zone is stationary (for example, a stator row), then the Reference Frame should be set to global.

    • If the cell zone is moving (for example, a rotor), then selecting global for Reference Frame keeps the holes stationary as the rotor moves past them. Such an arrangement would be appropriate for shroud cooling. However, if the holes are intended to cool the moving rotor blade, then the Reference Frame should be set such that it that tracks the rotor Cell Zone (for details on creating a Reference Frame, see Creating and Using Reference Frames). In this configuration, the holes will be stationary within the rotating frame and move position with the rotor.


    Note:  The hole coordinates do not have to be defined precisely on the mesh surface, and a tolerance is used to locate the holes on the boundary. However, if the hole location is too far from the surface, it will not be located and is ignored.


    An example of a hole definition profile is shown below:

    [Name]								
    Rotor Pressure Surface Cooling								
    								
    [Parameters]								
    InjectionTemperature = 420 [K]								
    MyParameter = 333 [s^-1]								
    MyParameter2 = 44.44 []								
    Diameter = 3.4 [mm]								
    								
    [Spatial Fields]								
    x	 y	 z						
    								
    [Data]			     Velocity in Stn  Velocity in Stn  Velocity in Stn   Mass Flow Some Variable  Sector
    x [m]	 y [m]     z [m]   Frame u [m s^-1] Frame v [m s^-1] Frame w [m s^-1]  [kg s^-1]  [kg s^-1]     Tag []
    -2.50E-01  1.08E-03  2.22E-01     2.63E+01        -9.87E+01	  5.09E+01	 3.71E-04     12.12         1
    -2.61E-01  1.06E-03  2.22E-01     1.08E+01        -1.03E+02	  5.32E+01	 3.71E-04     12.12         1
    -2.72E-01  1.04E-03  2.22E-01    -6.59E+00        -1.06E+02	  4.60E+01	 3.71E-04     12.12         1
    -2.83E-01  1.03E-03  2.22E-01    -2.97E+01        -1.10E+02	  4.41E+01	 3.71E-04     12.12         1
    -2.94E-01  1.01E-03  2.22E-01     2.07E+00        -1.12E+02	  3.32E+01	 3.87E-04     12.12         1
    #-- End of profile Rotor Pressure Surface Cooling--
    
    [Name]								
    Rotor Suction Surface Cooling								
    								
    [Parameters]								
    InjectionTemperature = 420 [K]								
    MyParameter = 333 [s^-1]								
    MyParameter2 = 44.44 []								
    Diameter = 3.4 [mm]								
    								
    [Spatial Fields]								
    x	 y	 z						
    								
    [Data]			     Velocity in Stn  Velocity in Stn  Velocity in Stn   Mass Flow Some Variable  Sector
    x [m]	 y [m]     z [m]   Frame u [m s^-1] Frame v [m s^-1] Frame w [m s^-1]  [kg s^-1]  [kg s^-1]     Tag []
    -2.97E-01  -1.22E-03  2.25E-01     2.63E+01        -9.87E+01	  5.09E+01	 3.71E-04     12.12         1
    -2.86E-01  -1.24E-03  2.22E-01     1.08E+01        -1.03E+02	  5.32E+01	 3.71E-04     12.12         1
    -2.75E-01  -1.26E-03  2.22E-01    -6.59E+00        -1.06E+02	  4.60E+01	 3.71E-04     12.12         1
    -2.63E-01  -1.29E-03  2.22E-01    -2.97E+01        -1.10E+02	  4.41E+01	 3.71E-04     12.12         1
    -2.52E-01  -1.32E-03  2.22E-01     2.07E+00        -1.12E+02	  3.32E+01	 3.87E-04     12.12         1
    #-- End of profile Rotor Pressure Surface Cooling--								
    

    Note:  Multiple hole definition profiles can included in one .csv file as shown in the example.


  3. Create a virtual geometry to represent the holes from the Create/Edit Boundary Geometry dialog box.

     Setup Boundary Interfaces Create/Edit Geometry

    1. Under Boundary Geometry, enter a name for the virtual hole geometry.

    2. Under Coordinate Profile, select a profile to define the hole center locations. The hole center profile must be defined in .csv file format.

    3. Under Geometry BC-type, select the type of boundary for the holes. The options are mass-flow-inlet, mass-flow-outlet, pressure-inlet, or wall. In most cases for blade film cooling, you will specify the holes as mass flow inlets, although for a bleed you would require an outlet.

    4. Under Hole Direction Specification, select one of the following options:

      • normal-to-boundary

        The holes are oriented normal to the surface based on their location.

      • cartesian

        The holes are oriented by specifying a X, Y, and Z Direction Component. Using a constantdirection component orients all holes in the same direction. If you want individual holes to be oriented in different directions, select profile and choose the appropriate Profile Variable to assign the direction components for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to a Direction Component for all holes.

      • cylindrical

        The holes are oriented by specifying a Radial, Tangential, and Axial Direction Component. The axial direction will be defined by the Reference Frame properties of the parent cell zone. Using a constant direction component orients all holes in the same direction. If you want individual holes to be oriented in different directions, select profile and choose the appropriate Profile Variable to assign the direction components for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to a Direction Component for all holes.

    5. Under Flow Direction Specification, you have the option of specifying the flow direction independently from the hole direction. You have the following options:

      • hole-direction

        The flow is aligned with the hole direction.

      • from-boundary-panel

        The flow direction is taken from the overlap boundary condition specification. For details, see Specifying Overlap Boundary Conditions.

      • normal-to-boundary

        The flow is oriented normal to the surface regardless of hole direction.

      • cartesian

        The flow is oriented by specifying a X, Y, and Z Component of Flow. Using a constant Component of Flow orients the flow in the same direction for each hole. If you want individual holes to have different flow directions, select profile and choose the appropriate Profile Variable to assign the flow direction for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to a Component of Flow for all holes.

      • cylindrical

        The flow is oriented by specifying a Radial, Tangential, and Axial Component of Flow. The axial direction will be defined by the Reference Frame properties of the parent cell zone. Using a constantComponent of Flow orients the flow in the same direction for each hole. If you want individual holes to have different flow directions, select profile and choose the appropriate Profile Variable to assign the flow direction for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to a Component of Flow for all holes.

    6. Under Hole Shape, select cylindrical or rectangle and specify the relevant dimensions.

      • For cylindrical holes, specify the Diameter as constant if each hole is the same size or select profile and choose the appropriate Profile Variable to assign the diameter for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to all holes.

      • For rectangular holes, specify the Length and Width as constant if each hole is the same size or select profile and choose the appropriate Profile Variable to assign the width or length for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to all holes.


      Note:  For both cylindrical and rectangular hole shapes, the orientation of the hole is based on the projection of the specified hole direction onto the local underlying boundary surface. If the specified hole direction is not normal to the boundary, then the hole shape is modified to account for this as if the hole were "drilled" into the boundary in the specified direction. As such, a cylindrical hole definition can result in an elliptical shape on the boundary, and similarly a rectangle defined with equal width and length can result in a non-square hole, as shown in Figure 14.20: Hole Shape Dependency from Hole Direction.


      Figure 14.20: Hole Shape Dependency from Hole Direction

      Hole Shape Dependency from Hole Direction

    7. Under Flow Variables, there exist additional settings depending on the physics present in the simulation. If Energy is enabled, and the Geometry BC-type is an inlet or wall, then a Temperature specification is required. If the Geometry BC-type is an inlet or an outlet, then a Mass Flow rate or Total Pressure specification is required.

      For Mass Flow, the following options are available:

      • from-boundary-panel

        The mass flow is taken from the overlap boundary condition specification. For details, see Specifying Overlap Boundary Conditions.

      • profile

        If you want individual holes to have different mass flow rates, select profile and choose the appropriate Profile Variable to assign the mass flow rate for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to all holes.

      • constant-per-hole

        A single value for the mass flow rate is applied to each hole.

      • total-mass-flow

        This value corresponds to the total mass flow for all the located holes, which is distributed equally between those holes.

      • total-for-all-sectors

        This value corresponds to the total mass flow of all the located holes, assuming a full 360 degree wheel were being modeled. The specified mass flow is distributed proportionately to any modeled sectors, and then equally between the holes within those sectors.

      For Temperature, the following options are available:

      • from-boundary-panel

        The temperature is taken from the overlap boundary condition specification. For details, see Specifying Overlap Boundary Conditions.

      • profile

        If you want individual holes to have different temperatures, select profile and choose the appropriate Profile Variable to assign the temperature for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to all holes.

      • radial-distribution

        Allows for a temperature variation within each individual hole. The temperature for each face within a hole is defined by the equation:

        (14–14)

        Where A and B are user-defined constants and is the Reference Temperature, which can be defined as constant under Reference Value or select profile to use hole-specific values for Reference Temperature as defined from the input profile.

        The resultant temperature variation is quadratic, based on the radial position of the hole face centers from the input hole center, compared to the input hole radius. The value of k is evaluated automatically by Fluent, to ensure the enthalpy flux for each hole is approximately the same as if the hole were specified with a uniform temperature equal to the Reference Temperature. This helps to reduce the sensitivity of the model to the mesh resolution of the holes themselves.

      For a pressure-inlet, a Total Pressure and Loss Coefficient specification is required.

      For total pressure, the following options are available:

      • from-boundary-panel

        The total pressure is taken from the overlap boundary condition specification. For details, see Specifying Overlap Boundary Conditions.

      • profile

        If you want individual holes to have different total pressure specifications, select profile and choose the appropriate Profile Variable to assign the total pressure for each individual hole based on data from the profile. You can also select a profile Parameter which will apply the constant value to all holes.

      The Loss Coefficient reproduces the pressure head loss due to the hole drilling itself. For details on the loss coefficient, see Specifying the Loss Coefficient. The following options are available:

      • from-boundary-panel

        The loss coefficient is taken from the overlap boundary condition specification. For details, see Specifying Overlap Boundary Conditions.

      • constant

        A single, constant value is applied to all holes.

      • factor-and-length-ratio

        Two further values are required for an internal calculation of the loss coefficient. Both are constant for all holes: A Darcy Friction Factor and an assumed Length-to-Diameter Ratio for the hole drilling.

    8. Once fully defined, click Create/Edit. Note that existing geometry definitions can be modified and you can define multiple geometries if needed. It is often useful to group holes into separate geometries according to their location: for example, blade pressure side, blade suction side, hub, shroud, and so on.

14.7.2. Specifying the Boundary Interface

After the hole geometry has been defined, you must create the Boundary Interface between the geometry and an existing boundary. To do this, right-click Boundary Interfaces in the outline view and select New to open the Create/Edit Boundary Interfaces dialog box.

 Setup Boundary Interfaces New

  1. Under Boundary Interface, enter a name for the interface.

  2. Under Boundary Zones Side 1, select the boundary that the holes will be located on.

  3. Under Geometry Zones Side 2, select the boundary geometry that was created in Specifying the Virtual Hole Geometry.

  4. Click Preview and the selected Boundary Zones will be displayed, along with symbols indicating the proposed locations of the holes as implied by the selected Geometry Zones.

    Note that the symbols do not indicate the actual size and shape of the holes; they are only markers. The default symbols and colors are used but these can be changed using the following console command for each Geometry: /define/virtual-boundary/hole-geometry/edit/[geometry-name]/preview.

  5. Click Create/Edit and the boundary interface will appear in the list. You can also select previous defined boundary interfaces to edit or delete.

Once created, Fluent intersects the selected boundaries with the selected geometries, resulting in one overlap boundary for each geometry and separate non-overlap boundaries for each of the underlying boundary zone locations, as shown in Figure 14.21: Creation of Overlap and Non-overlap Surfaces. The new overlap boundaries will automatically be created as mass-flow inlets or outlets, or an inlet vent (for pressure), depending on the type defined within the geometry.


Note:   If a profile option, or any option other than constant or from-boundary-panel is selected for Flow Direction Specification or within Flow Variables for the relevant Geometry Zone, then Fluent creates a special internal profile that is only valid for this overlap boundary. It has a name beginning with vbr and also contains the overlap boundary name. Fluent then assigns this internal profile to the relevant boundary conditions on the overlap boundary, corresponding to the Boundary Geometry specification. The internal profile will not be visible or usable anywhere but on the appropriate overlap boundary.


Figure 14.21: Creation of Overlap and Non-overlap Surfaces

Creation of Overlap and Non-overlap Surfaces

14.7.2.1. Hole Locations for Rotationally Periodic Interface Pairs

In the configuration where you have a rotationally periodic model with hole locations specified for a Boundary Zone that fall outside of that Boundary Zone and the periodic boundaries, then Fluent will attempt to rotate the hole locations until they overlap with the Boundary Zone as shown in Figure 14.22: Holes Outside the Parent Cell Zone Rotated to Overlap with Mesh

This is triggered automatically when a Boundary Interface is created using a Boundary Zone that is attached to a parent Cell Zone that has an axis vector defined and at least on rotationally periodic interface pair. Any hole locations specified within the Boundary Geometry that are not initially overlapping with the parent Cell Zone mesh will be rotated (using the pitch angle defined by the periodic pair) until they do overlap. This rotation moves the hole location from just outside one periodic boundary, to just inside the other periodic boundary, as if it were located in an adjacent passage. If, even after successive rotation, the hole cannot be located within the mesh, then that hole will remain undefined.

This automated feature is designed to lower the sensitivity of the positioning of the rotationally periodic boundary pairs, in relation to hole locations.

Figure 14.22: Holes Outside the Parent Cell Zone Rotated to Overlap with Mesh

Holes Outside the Parent Cell Zone Rotated to Overlap with Mesh

Additionally, if a hole is large enough to touch or cross one of the periodic boundaries, then the portion of the hole outside that boundary will be transformed so that it appears inside the other periodic boundary. The complete hole will then appear to span the periodic pair.

14.7.3. Boundary Interface Definitions

After a Boundary Interface has been created, you can expand its entry in the outline view to see all of the components that make up that boundary interface.

 Setup Boundary Interfaces[Name]

In this example, expanding the Boundary Interface int-cool-hub reveals the following:

  • Boundary zone (rotor-hub-passage-2)

  • Non-overlap zone (rotor-hub-passage-2-non-overlapping)

  • Geometry (cool-hub)

  • Overlap zone (int-cool-hub-overlap-1-1)

You can right-click the overlap or non-overlap zone and select Display to view the zone in the graphics viewer. The non-overlap zone displays the non-intersected faces and the overlap zone displays the intersected mesh, revealing the exact hole shapes and sizes.

14.7.4. Specifying Overlap Boundary Conditions

Once the Boundary Interface is created, you may need to specify additional boundary conditions at the new overlap boundary.

  1. The boundary type will already be set based on the original geometry definition (Mass-Flow Inlet, Mass-Flow Outlet, or Wall).

  2. The overlap boundary will have some of its options automatically set by Fluent, depending on the associated Boundary Geometry specification. For example, if you selected profile or constant-per-hole for the relevant Boundary Geometry, then Fluent will automatically assign the relevant internal profile to the Mass Flux specification. However, if the from-boundary-panel option is chosen then you must manually specify the mass flow (or mass flux) directly. This value then corresponds to the total mass flow for all the located holes which is then distributed between those holes based on area-weighting.

  3. Verify the flow direction for each hole. The Direction Specification Method, Coordinate System, and the Flow Direction components should all have been set automatically according to the Boundary Geometry specification. Note that if Fluent has not set these options, and instead you specify the flow direction components as constant values, then each hole will have the same flow direction.

  4. Set appropriate values for other inflow quantities, where necessary, such as turbulence or temperature (shown below). For a quantity such as temperature, hole-specific data can be used by selecting the profile option for Temperature on the relevant Boundary Geometry, and choosing an additional data column from the original geometry definition profile.


    Important:  You are still able to manually edit the overlap Boundary Condition and change any automatically generated settings. However, this can lead to inconsistency with the Boundary Geometry specification. It is generally recommended that you change any overlap Boundary Condition indirectly by editing the relevant geometry (see Specifying the Virtual Hole Geometry).


14.7.5. Limitations for Boundary Interfaces and their Geometries

  • Only available for boundaries, and not volumes. Flow injection is supported at a Boundary Interface only, which must refer to an existing boundary definition.

  • Holes must have their size and shape defined such that they do not overlap with another hole.

  • Hole centers must be located within, or close to, a Cell Zone mesh. If any are not, then a warning will be printed that those holes will not be included and they will not contribute to fluid flow.

  • Flow boundary conditions are currently specified as constant for each hole, and cannot vary across an individual hole.

  • Reported boundary fluxes will show the specified mass flow for each hole overlap boundary, as well as the total mass flow for all geometries intersecting the original underlying boundary. However, the Net Results will not include this duplication and should still tend to zero across all boundaries for a converged solution.

  • When modeling Blade Film Cooling on a moving rotor blade in a transient simulation, where the Boundary Interface Geometry references a profile defined in a local reference frame rotating with the rotor, an incorrect Mass Flow will be applied to the holes unless the from-boundary-panel option has been chosen for Mass Flow on the geometry. Also, Fluent may become unstable if the rotation angle becomes greater than one pitch.

    Even if the from-boundary-panel option has been selected, the hole positions may not be updated correctly with the rotor position if more than one Boundary Interface is defined.