VM-LSDYNA-SOLVE-001

VM-LSDYNA-SOLVE-001
Steady State Thermal Analysis of 3D Hollow Spheres with Convection Boundary Condition

Overview

Reference: Kreith, F. (1959). Principles of Heat Transfer (2nd ed.). International Textbook Co.
Analysis Type(s): Steady State Thermal
Element Type(s): Solid
Input Files:Link to Input Files Download Page

Test Case

A hollow sphere with an inner radius ra and an outer radious rb has an inner surface temp of 100°C. The material of the hollow sphere is isotropic. A convective boundary condition is applied at the outer surface rb. The external ambient temperature T is 0°C, and the convective heat transfer coefficient h is 1 W/(m2°C). Compute the Steady state temperature at rc = 3.51419 m.

Figure 14: Problem Sketch

Problem Sketch

Material PropertiesGeometric PropertiesLoading
k = 1 W/(m°C)ra = 2 mTa = 100°C
h = 1 W/(m2°C)rb = 5.02839 mT = 0°C

Analysis Assumptions and Modeling Notes

LS-DYNA Thermal Solver 11 is used for this example. Three hollow spheres are modeled using ELFORM =  1, 4, 16. The temperature and convection boundary conditions are applied using the keywords *BOUNDARY_TEMPERATURE and *BOUNDARY_CONVECTION respectively. One eighth symmetry is used, and the heat flow is limited to the radial direction only.

Figure 15: Heat Flow Schematic

Heat Flow Schematic

At steady state, the rate of heat dissipation Q is

where the thermal resistances Rac, Rcb, and Rb∞ can be expressed as follows:

From the above equations, the analytical solution to Tc is 36.7617°C.

Figure 16: Hollow sphere models with hex, tet4 and tet10 element meshes

Hollow sphere models with hex, tet4 and tet10 element meshes

Results Comparison

Result Target LS-DYNA Error (%)
Temperature at rc (hex)36.7617°C36.9823°C.6%
Temperature at rc (tet4)36.7617°C36.9371°C.48%
Temperature at rc (tet10)36.7617°C36.3149°C-1.2%