ACU-T: 3204 Radiation Heat Transfer in a Simple Headlamp using the Discrete Ordinate Model

Prerequisites

This tutorial introduces you to setting up a radiation heat transfer problem using the Discrete Ordinate radiation model in HyperWorks CFD. Prior to starting this tutorial, you should have already run through the introductory HyperWorks tutorial, ACU-T: 1000 HyperWorks UI Introduction, and have a basic understanding of HyperWorks CFD, AcuSolve, and HyperView. To run this simulation, you will need access to a licensed version of HyperWorks CFD and AcuSolve.

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

Problem Description

The problem to be solved is shown schematically in Figure 1 and Figure 2. It consists of a simple headlamp with a housing, lens, and a bulb. The inner cavity of the bulb is filled with air and the wattage of the bulb is 1W, which is modeled as a volumetric heat source. The Boussinesq density model is used for the air to consider the natural convection effects in the fluid volume. The heat generated in the bulb is transferred by three means: conduction from the bulb to the housing, natural convection in the air volume, and radiation from the bulb to the air and lens volume. The outer surfaces of the headlamp are assumed to be at a constant temperature of 300 K.


Figure 1.


Figure 2.
The air and the lens are modeled as participating media for radiation heat transfer. For the purpose of this tutorial, the absorption coefficient of air is considered to be zero. The following radiation material properties are used for the air and the lens volumes.
Table 1.
  Absorption Coefficient Refractive Index
Air 0 1.0
Lens 900 1.57
For simulations involving participating media, the following radiation surface types must be defined at the boundaries of each participating medium.
Table 2.
Surface Type Radiation Surface Type
Participating medium - Participating medium Interface Radiation Interface - Internal
Participating medium - (Non- Participating) medium Interface Wall
External boundaries of Participating medium Radiation Interface - External or Wall
The discrete ordinate radiation model implementation in AcuSolve allows for modeling specular, diffuse, and partially specular interfaces. This can be done by specifying the appropriate value of diffused fraction while defining the radiation surfaces. For the diffused fraction, a value of 1 means the surface is purely diffuse, a value of 0 means it is purely specular. Any value between 0 and 1 means the surface is partially specular. In this tutorial, since both the air and lens are modeled as participating media, the lens-air interface and the lens-outer surfaces will be modeled as internal and external radiation interface respectively. The internal interface will be defined as a specular interface whereas the external interface will be defined as a diffuse interface.


Figure 3.

Start HyperWorks CFD and Create the HyperMesh Model Database

  1. Start HyperWorks CFD from the Windows Start menu by clicking Start > Altair <version> > HyperWorks CFD.
    When HyperWorks CFD is loaded, the Geometry ribbon is open by default.
  2. Create a new .hm database in one of the following ways:
    • From the menu bar, click File > Save.
    • From the Home tools, Files tool group, click the Save As tool.


      Figure 4.
  3. In the Save File As dialog, navigate to the directory where you would like to save the database.
  4. Enter Headlamp_DO as the name for the database then click Save.
    This will be your problem directory and all the files related to the simulation will be stored in this location.

Import and Validate the Geometry

Import the Geometry

  1. From the menu bar, click File > Import > Geometry Model.
  2. In the Import File dialog, browse to your working directory then select ACU-T3204_headlamp.x_t and click Open.
  3. In the Geometry Import Options dialog, leave all the default options unchanged then click Import.


    Figure 5.


    Figure 6.

Validate the Geometry

  1. From the Geometry ribbon, click the Validate tool.


    Figure 7.
    The Validate tool scans through the entire model, performs checks on the surfaces and solids, and flags any defects in the geometry, such as free edges, closed shells, intersections, duplicates, and slivers.
    The surface and solid errors display in the list below the tool.


    Figure 8.
  2. Click SolidChecks.
    The Solid Repair tool opens, which you can use to fix the geometric errors in the model.
    From the SolidChecks legend, you can see the model's solids have five intersections.


    Figure 9.
  3. Click Intersections.
    A guide bar used to fix intersecting solids displays.
  4. Optional: Click and to review each error.
  5. Activate the Keep common interface option then click Combine All.
    The SolidChecks legend should now display zero for all errors.
  6. Click the Validate tool once again.
    Observe that a blue check mark now appears on the top-left corner of the tool icon. This indicates that no issues are detected and you are ready to continue.


    Figure 10.

Set Up Flow

Set the General Simulation Parameters

  1. From the Flow ribbon, click the Physics tool.


    Figure 11.
    The Setup dialog opens.
  2. Click the Time setting and ensure that Steady is selected.


    Figure 12.
  3. Click the Flow setting and change the option to Laminar.
  4. Activate the Include gravitational acceleration option and set the gravity in the y-direction to -9.81.


    Figure 13.
  5. Click the Heat Transfer setting and activate the Heat transfer option.


    Figure 14.
  6. Click the Solver controls setting and activate the Thermal flow equation.


    Figure 15.
  7. Close the dialog and save the model.

Define the Material Models

  1. From the Flow ribbon, click the Material Library tool.


    Figure 16.
    The Material Library dialog opens.
  2. Click the My Materials tab.
  3. Click to add a new fluid material model.
  4. In the material creation dialog, click the name in the top-left corner and rename the material to Air_Boussinesq.
  5. In the Density tab,
    1. Set the Type to Boussinesq.
    2. Set the Density value to 1.225.
    3. Set the Expansivity value to 0.00347222.
    4. Reference temperature value to 288.


    Figure 17.
  6. Click the Specific Heat tab and set the Specific heat value to 1005.


    Figure 18.
  7. Click the Viscosity tab and set the Viscosity value to 1.781e-05.


    Figure 19.
  8. Click the Conductivity tab and set the Conductivity value to 0.02521.


    Figure 20.
  9. Close the material creation dialog to return to the Material Library dialog.
  10. Select Solid in the Settings menu, click the My Materials tab, then click to create a new solid material model.
  11. Name the material Plastic and set the following values.
    The Type should be Constant for each property.
    • Density: 1270
    • Specific Heat: 1900
    • Conductivity: 0.2
  12. Close the material creation dialog to return to the Material Library dialog.
  13. Similarly, create new solid material models named Arnite and LED with the following properties.
    The Type should be Constant for each property.
    Arnite:
    • Density: 1670
    • Specific Heat: 2050
    • Conductivity: 1.6
    LED:
    • Density: 5500
    • Specific Heat: 0.3
    • Conductivity: 5.0


    Figure 21.
  14. Close all dialogs and save the model.

Assign Material Properties

  1. From the Flow ribbon, click the Material tool.


    Figure 22.
  2. Click the lens volume highlighted in the figure below and select Arnite from the Material drop-down menu.


    Figure 23.
  3. On the guide bar, click to execute the command and remain in the tool.
  4. Click the housing volume and assign the Plastic material model.


    Figure 24.
  5. On the guide bar, click to execute the command and remain in the tool.
  6. In the Materials legend, right-click on Air and select Isolate.
  7. Click the air volume and assign the Air_Boussinesq material model.


    Figure 25.
  8. On the guide bar, click to execute the command and remain in the tool.
  9. In the Materials legend, right-click on Air and select Isolate.
  10. Click the bulb volume and assign the LED material model.


    Figure 26.
  11. On the guide bar, click to execute the command and exit the tool.
  12. Save the model.

Define the Heat Source

  1. From the Flow ribbon, click the Sources > Heat tool.


    Figure 27.
  2. In the modeling window, select the bulb volume.
  3. In the Heat Source dialog, set the heat source value to 2049180 W/m3.


    Figure 28.
  4. On the guide bar, click to execute the command and exit the tool.
  5. Press Esc to exit the Sources tool then press the A key to turn on the display of all solids.
  6. Save the model.

Define Flow Boundary Conditions

In this problem, all the surfaces are walls and will therefore be assigned the default wall boundary condition. The outer walls of the headlamp will be given a no-slip wall boundary condition with a constant temperature.

  1. From the Flow ribbon, click the No Slip tool.


    Figure 29.
  2. In the modeling window, select the surfaces highlighted in the figure below.


    Figure 30.
  3. In the microdialog, enter the values shown in the figure below.


    Figure 31.
  4. In the Boundaries legend, double-click Wall, rename it to Outerwalls, then press Enter.
  5. On the guide bar, click to execute the command and exit the tool.
  6. Save the model.

Set Up Radiation

In this step, you will specify the parameters related to the thermal radiation setup.

Define Radiation Model Settings

  1. From the Radiation ribbon, Thermal Radiation tools, click the Physics tool.


    Figure 32.
    The Radiation Settings dialog opens.
  2. Activate Thermal radiation and set the Radiation model to Discrete Ordinate.
  3. Verify that the Radiation quadrature is set to S4.


    Figure 33.
  4. Close the dialog.

Define the Emissivity Models

  1. From the Radiation ribbon, click the Surface Finish Library tool.


    Figure 34.
    The Surface finish library opens.
    Click to add a new emissivity model.
  2. Set the Name of the model to Inner and the Emissivity value to 0.05 by double-clicking on the entity fields.


    Figure 35.
  3. Close the dialog.

Define the Participating Media Radiation Models

  1. From the Radiation ribbon, Participating Media tools, click the Model tool.


    Figure 36.
    The Participating media model library opens.
  2. Click to add a new model.
  3. Name the new model Air and assign the following property values:
    • Absorption coefficient - 0
    • Scattering coefficient - 0
    • Phase coefficient - 0
    • Refractive index - 1
  4. Similarly, create another model named Lens and assign the following property values:
    • Absorption coefficient - 900
    • Scattering coefficient - 0
    • Phase coefficient - 0
    • Refractive index - 1.57


    Figure 37.
  5. Close the dialog.

Assign Participating Media Models

  1. From the Radiation ribbon, click the Assign tool.


    Figure 38.
  2. Click the lens volume and assign the Lens model from the microdialog drop-down.


    Figure 39.
  3. On the guide bar, click to execute the command and remain in the tool.
  4. In the Participating Media legend, right-click on Lens and select Hide.
  5. Select the air volume and assign the Air model.


    Figure 40.
  6. On the guide bar, click to execute the command and remain in the tool.
  7. In the Participating Media legend, right-click on Air and select Hide.
  8. Save the model.

Assign Surface Finish Models

  1. From the Radiation ribbon, click the Surface Finish tool.


    Figure 41.
  2. Select only the inner surfaces of the housing - i.e. the housing-air interface.


    Figure 42.

    A total of nine surfaces should be selected.

  3. Assign the Inner model from the microdialog drop-down
  4. On the guide bar, click to execute the command and exit the tool.
  5. Turn on the display of all surfaces in the model by pressing the A key.

Define the Radiation Interfaces

Since the lens and air are defined as semi-transparent media, the interface between the two will be modeled as a specular interface of type: internal. Accordingly, the outer surfaces of the lens will be modeled as external interface radiation surface.

  1. From the Radiation ribbon, click the Interface tool.


    Figure 43.
    Notice that all the surfaces except the Lens-Air interface surface are made transparent. Since the Air and Lens volumes are defined as participating media, HyperWorks CFD automatically detects the surfaces that can be defined as an internal radiation interface and only those surfaces will be available for selection when this tool is active.
  2. Using the box selection method, draw a window around the model.
    Observe that the lens-air interface is the only surface that is selected.
  3. In the microdialog, set the Radiation diffused fraction to 0.


    Figure 44.
  4. On the guide bar, click to execute the command and exit the tool.
  5. Click the External Boundary tool.


    Figure 45.
  6. Select the three outer surfaces of the lens volume.
  7. In the microdialog, set the Radiation diffused fraction to 1, the External refractive index to 1, and the External temperature to 300 K.


    Figure 46.
  8. On the guide bar, click to execute the command and exit the tool.
  9. Save the model.

Generate the Mesh

In this step, you will define the mesh controls and then generate the mesh.

Define the Surface Mesh Controls

  1. From the Mesh ribbon, click the Surface tool.


    Figure 47.
  2. Using the box selection method, select all the surfaces in the model.
  3. In the microdialog, set the Average element size to 0.002.


    Figure 48.
  4. On the guide bar, click to execute the command and remain in the tool.
  5. Select the surface highlighted in the figure below and hide it by pressing the H key or by right-clicking and select Hide from the context menu.


    Figure 49.
  6. Select the three surfaces highlighted in the figure below.
  7. In the microdialog, set the Average element size to 0.001 and the Mesh growth rate to 1.


    Figure 50.
  8. On the guide bar, click to execute the command and exit the tool.

Define the Boundary Layer Controls

  1. From the Mesh ribbon, click the Boundary Layer tool.


    Figure 51.
  2. Right-click in the modeling window and go to Select > Advanced Select > By Material > Air_Boussinesq.
    All the fluid wall surfaces should be selected and a microdialog for BL specification appears.
  3. Enter the following values in the microdialog:
    • First layer thickness definition: Constant
    • First layer thickness: 0.0005
    • Total number of layers: 4
    • Growth method: Constant
    • Initial growth rate: 1.3
    • Termination policy: Truncate
    • Activate the Enable surface mesh modification option


    Figure 52.
  4. On the guide bar, click to execute the command and exit the tool.
  5. Turn on the display of all the surfaces.

Define the Volume Mesh Controls

Since the thickness of the housing and the lens solids are small, you will use the thin layer meshing tool so that when the volume mesh is generated, there will be two layers across the thickness of those solids.

  1. From the Mesh ribbon, click the Volume Mesh tool.


    Figure 53.
  2. Select the housing and the lens solids.
  3. In the microdialog,
    1. Set the Average size to 0.001.
    2. Set the Growth rate to 1.0.
    3. Activate the Thin layer meshing option and set the Number of layers to 2.


    Figure 54.
  4. On the guide bar, click to execute the command and exit the tool.

Generate the Mesh

  1. From the Mesh ribbon, click the Batch tool.


    Figure 55.
    The Meshing Operations dialog opens.
  2. Set the Mesh growth rate to 1.
  3. Click Mesh.
    The Run Status dialog opens. Once the run is complete, the status is updated and you can close the dialog.
    Tip: Right-click on the mesh job and select View log file to view a summary of the meshing process.
  4. Save the model.

Run AcuSolve

  1. From the Solution ribbon, click the Run tool.


    Figure 56.
    The Launch AcuSolve dialog opens.
  2. Set the Parallel processing option to Intel MPI.
  3. Optional: Set the number of processors to 4 or 8 based on availability.
  4. Expand the Default initial conditions menu and deactivate the Pre-compute flow option.
  5. Set the x-velocity to 0 and the Temperature to 300.
  6. Leave the remaining options as default and click Run to launch AcuSolve.


    Figure 57.
    The Run Status dialog opens. Once the run is complete, the status is updated and you can close the dialog.
    Tip: While AcuSolve is running, right-click on the AcuSolve job in the Run Status dialog and select View Log File to monitor the AcuSolve solution process.

Post-Process the Results With HyperView

Open HyperView and Load the Model and Results

  1. Start HyperView from the Windows Start menu by clicking Start > All Programs > Altair <version> > HyperView.
    Once the HyperView window is loaded, the Load model and results panel should be open by default. If you do not see the panel, click File > Open > Model.
  2. In the Load model and results panel, click next to Load model.
  3. In the Load Model File dialog, navigate to your working directory and select the AcuSolve .Log file for the solution run that you want to post-process. In this example, the file to be selected is Headlamp_DO.1.Log.
  4. Click Open.
  5. Click Apply in the panel area to load the model and results.
    The model is colored by geometry after loading.

Create Temperature and Incident Radiation Contours on a Section Cut

  1. Orient the display to the xy-plane by clicking on the Standard Views toolbar.
  2. On the 3DViewControls toolbar, right-click on multiple times until the orientation of the model is as shown in the figure below.


    Figure 58.
  3. Click the Section cut icon on the HV-Display toolbar.
  4. In the panel area, click Add to create a new section cut named Section 1.
  5. In the Define plane section, set the axis to X Axis then click Apply.
  6. Change the Display options from Clipping plane to Cross section.


    Figure 59.
  7. Click Gridline. In the Gridline Options dialog, deactivate the Show check box under Grid line then click OK.
  8. Click on the Results toolbar to open the Contour panel.
  9. In the panel area, set the Result type to Temperature (s).
  10. Click the Components entity selector. In the Extended Entity Selection dialog, select All.
  11. Click Apply.
  12. In the panel area, under the Display tab, turn off the Discrete color option.


    Figure 60.
  13. Click the Legend tab then click Edit Legend. In the dialog, change the Numeric precision to 2 then click OK.
    The contour plot should look like the one shown in the figure below.


    Figure 61.
  14. In the panel area, change the Result type to Incident_radiation then click Apply.


    Figure 62.

    As seen in the figure above, the rays travel from the bulb through the air and lens before getting radiated to the atmosphere.

Summary

In this tutorial, you learned how to set up and solve a radiation heat transfer problem in a headlamp using the discrete ordinate radiation model in HyperWorks CFD. You started by importing and cleaning up the geometry and then set up the flow, thermal, and radiation boundary conditions. Once the solution was computed, you processed the results using HyperView, where you created contour plots of temperature and incident radiation.