Transient Thermal Analysis in ANSYS® Mechanical (Workbench): Dealing with Non-physical Temperature Results
|Figure 1: |
Temperatures at First Recorded Substep - Results Outside Bounds
When Transient Thermal Analysis is performed in ANSYS, whether via the APDL interface or Mechanical (Workbench), there are circumstances in which non-physical results can occur. An example is a temperature result that is outside any temperature applied to a model. This may be seen with extreme Biot numbers (high convection coefficients) or with inappropriate time substep sizes, and is more common with high-order thermal elements.
This article shows an example of this difficulty, and gives suggestions on reducing its occurrence.
A Transient Thermal Model
In a test model, a rectangular block with an initial temperature of 300C has had a convection load of a fluid at 22C and a convection coefficient of 0.01 W/mm2 /C applied to five of its six faces:
Figure 2: Convection Load on 5 Faces of a Block
After a short substep, temperatures that are greater than 300C are seen... this is not a physically realistic response:
Figure 3: Non-physical Temperature Response
The above temperature result is affected by the use of high-order elements, tetrahedral shapes, and high Biot numbers. Means of reducing the problem include smaller elements, low-order elements, and hex meshing.
The same model meshed with hex element of the same general size produces a not quite perfect, but much better result at the same substep:
Figure 4: Result with High-Order Hex Elements
Temperatures just slightly above 300C are seen. The non-physical temperatures are not seen if low-order elements are employed:
Figure 5: Result with Low-Order Elements
Thermal results are often moved into structural models. In ANSYS Workbench, the same geometry can have a different mesh applied, and the temperatures can be mapped from the thermal to the structural system.
Figure 6: Mapping Temperatures between Thermal and Structural Systems with Differing Meshes of the Same Geometry
A thermal system can use a mesh of low-order hex and tet elements in order to avoid non-physical results, while the associated structural system can mesh the same geometry with high-order elements in order to give good structural results, particularly with tet elements. The mapping includes controls for the imported body temperature in the structural system:
Figure 7: Imported Temperatures Mapped from a Different Thermal Mesh
In thermal FEA models, choices of elements size, shape and order, as well as high Biot number convective loads, can sometimes result in non-physical temperature results such as temperatures that are higher or lower than any applied temperature. In transient models, the use of small time substeps can amplify the effect with high-order elements.
Fewer problems of this sort are seen in thermal models that use low-order elements such as 4-node tet elements and 8-node brick elements. Related structural FEA models of the same geometry can use high-order structural elements, and recent versions of ANSYS Mechanical (Workbench), such as v14.0, can map temperatures between the non-matching meshes. There are a number of user-set controls for how the mapping is performed.
Difficulties in thermal responses will still be occasionally seen with low-order thermal elements. The use of layers of elements that are thin at exterior surfaces can be used in attempts to address this (reducing the Biot number). The difficulty may also sometimes be seen in thermal elements that have different convective loads on two faces of one element. Work-around methods could include small elements on edges, or a small strip on one side of an edge with no convective load applied.
The ANSYS Mechanical APDL interface has mapping methods for temperatures using the BFINT command. The mapping methods available in the most recent Workbench Mechanical interface are more advanced, and offer higher success rates for mapping between non-matching meshes. In this way, many non-physical thermal response difficulties can be avoided, while still performing related structural analyses with high-order elements.