OS-T: 1085 Linear Steady-state Heat Convection Analysis

This tutorial performs a heat transfer analysis on a steel pipe.

The temperature on the inside surface of the pipe is 60°C. The outside surface is exposed to the surrounding air, which is at 20°C. The temperature distribution within the pipe can be determined by solving the linear steady-state heat conduction and convection solution.

rd1085_cylinder
Figure 1. Model Review

Launch HyperMesh and Set the OptiStruct User Profile

  1. Launch HyperMesh.
    The User Profile dialog opens.
  2. Select OptiStruct and click OK.
    This loads the user profile. It includes the appropriate template, macro menu, and import reader, paring down the functionality of HyperMesh to what is relevant for generating models for OptiStruct.

Import the Model

  1. Click File > Import > Solver Deck.
    An Import tab is added to your tab menu.
  2. For the File type, select OptiStruct.
  3. Select the Files icon files_panel.
    A Select OptiStruct file browser opens.
  4. Select the thermal.fem file you saved to your working directory. Refer to Access the Model Files.
  5. Click Open.
  6. Click Import, then click Close to close the Import tab.

Set Up the Model

Create Thermal Material and Properties

Create the material and property collectors before creating the component collectors.
  1. In the Model Browser, right-click and select Create > Material.
    A default MAT1 material displays in the Entity Editor.
  2. For Name, enter steel.
  3. Check the box in front of MAT4.
    MAT4 card image appears below MAT1 in the Entity Editor. The MAT1 card defines the isotropic structural material. MAT4 card is for the constant thermal material. MAT4 uses the same material ID as MAT1.
  4. Enter the following values for the material, steel, in the Entity Editor.
    [E] Young's modulus
    2.1 x 1011 Pa
    [NU] Poisson's ratio
    0.3
    [RHO] Material density
    7.9 x 103 Kg/m3
    [A] Thermal expansion coefficient
    1.0 x 10-5 / °C
    [K] Thermal conductivity
    73W / m °C
    [H] Heat transfer coefficient
    40W / m2 °C

    os_1080_steel
    Figure 2. Material Entity Editor

    A new material, steel, is created with both structural and thermal properties.

  5. In the Model Browser, right-click and select Create > Property.
    A default PSHELL property displays in the Entity Editor.
  6. For Name, enter solid.
  7. For Card Image, select PSOLID and click Yes to confirm.
  8. For Material, click Unspecified > Material.
  9. In the Select Material dialog, select steel and click OK.
    The property of the solid steel pipe has been created as 3D PSOLID. Material information is linked to this property.

Apply Thermal Loads and Boundary Conditions

In this exercise the thermal boundary conditions are applied on the model and saved in a predefined load collector spc_temp. A predefined node 4679 specifies the ambient temperature. A predefined node set node_temp contains the nodes on the inside surface of the pipe.

Create Temperature on the Inner Surface of the Pipe

  1. From the Analysis page, click constraints.
  2. Go to the create subpanel.
  3. Make sure the current selection field is set to nodes.
  4. Click nodes > by sets.
  5. Select node_temp and click select.
  6. Uncheck the box in front of dof1, dof2, dof3, dof4, dof5, and dof6 and verify that the entry fields are set to 0.0.
  7. Set load types = to SPC.
  8. Click create.
    This applies the temperature 0.0 on the inside nodes. In the next step, the temperature value is updated to 60.
  9. Click the Card edit icon entityCards-24.
  10. Click loads > by collector.
  11. Check the box in front of spc_temp and click select.
  12. Click config= and select const.
  13. Click type= and select SPC.
  14. Click edit.
  15. In the field of D, enter 60.0.
  16. Click return three times to go back to the Analysis page.

Create Ambient Temperature

Make sure spc_temp is the current load collector.
  1. From the Analysis page, select the Constraints panel.
  2. Go to the create subpanel.
  3. Click nodes > by id.
  4. Input the ID of the predefined node 4679.
    Node 4679 should be highlighted.
  5. Uncheck the box in front of dof1, dof2, dof3, dof4, dof5, and dof6 and verify that the entry fields are set to 0.0.
  6. Click create.
  7. Click the Card edit icon entityCards-24.
  8. Select loads entry.
  9. Select the ambient spc just created on the screen.
  10. Click config= and select const.
  11. Click type= and select SPC.
  12. Click edit.
  13. In the field of D, enter 20.0.
    The temperature boundary conditions are created, as shown below.

    rd1085_thermal_bc
    Figure 3. Thermal Boundary Conditions
  14. Click return three times to go back to the Analysis page.

Create CHBDYE Surface Elements

Surface elements are to be created to simulate the heat exchange between the solid pipe and the surrounding air. A predefined element set elem_convec, which contains the solid elements on the outer surface of the pipe, is used to define the surface elements.
  1. Click BCs > Create > Interfaces.
  2. For Name, enter convection.
  3. For Card Image, select CONVECTION from the drop-down menu.
  4. Click Color and select a color from the palette.
  5. Click MID to activate it.
  6. For Material, click Unspecified > Material.
  7. In the Select Material dialog, select steel and click OK.
    An element group convection and a free convection property PCONV are created.
  8. For Secondary Entity IDs, select Elements.
    A panel appears under the modeling window.
  9. Click on the switch button beside elems and select faces from the list.
  10. Click the highlighted solid elems and select by sets from the selection menu.
  11. Select element set elem_convec and click select.
  12. Click nodes in face nodes field.
  13. Select 4 nodes on the surface face of a solid element, as shown below.

    rd1085_4nodes
    Figure 4. Selected Surface Nodes on the Solid Element Outside the Pipe
  14. In the break angle = field, enter 89.0.
  15. Click add.
    This adds the CHBDYE surface elements to the solid elements on the outer surface following the same side convention, as shown below.

    rd1085_4nodes
    Figure 5. Surface Elements on the Outer Layer of the Pipe
  16. Click return to go back to the Create group window.

Define Convection Boundary Condition to Surface Elements

  1. Click the Card Edit icon entityCards-24.
  2. Select elems.
  3. Click elems > by group.
  4. Check the box in front of CONVECTION and click select.
  5. Click config= and select second4.
  6. Click type= and select CHBDYE4.
  7. Click edit and go to the CHBDYE Card Image panel.
  8. Check the box in front of CONV.
  9. Click TA1 and input the ambient node ID 4679, as shown below.

    rd1085_4679
    Figure 6. Define the Convection Boundary Condition
  10. Click return three times to go back to the Analysis page.

Create Heat Transfer Load Step

An OptiStruct steady-state heat convection loadstep is created, which references the thermal boundary conditions in the load collector spc_temp. The gradient, flux, and temperature output for the heat transfer analysis is also requested in the Loadsteps panel.
  1. In the Model Browser, right-click and select Create > Load Step.
  2. A default loadstep displays in the Entity Editor.
  3. For Name, enter heat_transfer.
  4. Click on the Analysis type field and select Heat transfer (steady-state) from the drop-down menu.
  5. For SPC, click Unspecified > Loadcol.
  6. In the Select Loadcol dialog, select spc_temp and click OK.
  7. Check the box next to Output.
  8. Activate the options of FLUX and THERMAL on the sub-list.
  9. Activate the FORMAT fields for both outputs and select H3D format.
  10. Activate the OPTION fields for both outputs and select ALL.
    The FORMAT and OUTPUT fields for THERMAL output may open up a new window. Click on the first field in the window to select the corresponding values.

    FLUX and THERMAL output can also be requested in the Control cards panel on the Analysis page.

Submit the Job

  1. From the Analysis page, click the OptiStruct panel.

    OS_1000_13_17
    Figure 7. Accessing the OptiStruct Panel
  2. Click save as.
  3. In the Save As dialog, specify location to write the OptiStruct model file and enter thermal_complete for filename.
    For OptiStruct input decks, .fem is the recommended extension.
  4. Click Save.
    The input file field displays the filename and location specified in the Save As dialog.
  5. Set the export options toggle to all.
  6. Set the run options toggle to analysis.
  7. Set the memory options toggle to memory default.
  8. Click OptiStruct to launch the OptiStruct job.
If the job is successful, new results files should be in the directory where the thermal_complete.fem was written. The thermal_complete.out file is a good place to look for error messages that could help debug the input deck if any errors are present.

View the Results

Gradient temperatures and flux contour results for the steady-state heat conduction analysis and the stress and displacement results for the structural analysis are computed from OptiStruct. HyperView will be used to post process the results.

  1. From the OptiStruct panel, click HyperView.
    HyperView is launched and the results are loaded. A message window appears to inform of the successful model and result files loading into HyperView.
  2. Click Close to close the message window, if one appears.
  3. On the Results toolbar, click resultsContour-16 to open the Contour panel.
  4. Select the first pull-down menu below Result type and select Grid Temperatures(s).
  5. Click Apply.
    You may have to use Edit Legend in the Contour panel to get the contour, as shown in Figure 8.
    A contour plot of grid temperatures is created.
  6. Select the first pull-down menu below Result type and select Element Fluxes (V).
  7. Click Apply.
    You may have to use Edit Legend in the Contour panel to get the contour. Both temperature and flux contour plots are shown in Figure 8.

    rd1085_results
    Figure 8. Results of Heat Transfer Analysis