Fill out the form to download

Required field
Required field
Not a valid email address
Required field

Documentation

Validation Case: Heat Transfer in a Perforated Plate

This project belongs to the heat transfer analysis type. The aim of this test case is to validate the following parameter using a transient thermal solver, with a surface and a convective heat flux boundary condition on a perforated plate:

  • Temperature \([K]\) on a single node for 12 consecutive timesteps.

The simulation results of SimScale were compared to the results presented in case B of [TTLP301]\(^1\).

Geometry

The geometry used for the analysis with the highlighted node N1 can be seen below. It is a random thick plate with straight and curved edges.

perforated plate geometry points
Figure 1: Geometry of a perforated plate with highlighted node N1 (red)

The coordinates for each vertex of the plate are displayed in the following table:

 ABCDEFGHIJ
x \([m]\)0.6350.93950.93950.54270.31750.6350.93950.93950.54280.3175
y \([m]\)000.83380.94010.5499000.83380.94010.5499
z \([m]\)000000.20.20.20.20.2
Table 1: The coordinates of the main geometry points

Analysis Type and Domain

Tool type : Code_Aster

Analysis type : Heat transfer, linear

Time dependency: Transient

Mesh and element types: Four mesh cases were considered for the analysis of the perforated plate; with 1st order and 2nd order tetrahedral, as well as 1st order and 2nd order hexahedral elements. The hexahedral mesh was locally computed and uploaded to the SimScale Workbench.

CaseMesh typeNumber of nodesNumber of 3D elementsElement type
(A)1st order tetrahedral534025676Standard
(B)2nd order tetrahedral1680710976Standard
(C)1st order hexahedral1210840Standard
(D)2nd order hexahedral4433840Standard
Table 2: The mesh characteristics for all the cases

Below, the 1st order tetrahedral mesh for case A is visualized:

mesh standard first order tetrahedral elements
Figure 2: The mesh used for case A, created with SimScale’s standard meshing algorithm

And the mesh that was uploaded for case C is presented below:

mesh standard first order hexahedral elements
Figure 3: The mesh used for case C, created externally with 1st order hexahedral elements

Simulation Setup

Material/Solid

  • Isotropic:
    • Density \(ρ\) = 1 \( kg \over \ m³ \),
    • Thermal conductivity \(\kappa\) = 0.1 \( W \over \ m \ K\),
    • Specific heat = 1 \( J \over \ kg \ K\)  

Initial conditions

  • Initial Temperature = 273.15 \(K\)  

Boundary conditions

  • Surface heat flux on face AFJE:
    • Heat flux value = 1 \(W \over \ m²\)  
  • Convective heat flux on face CHID
    • Reference temperature \(T_{0}\) = 273.15 \(K\)
    • Heat transfer coefficient = 1 \(W \over \ m² \ K\)  

Results Comparison

The temperature values at node N1 for 12 time steps obtained with SimScale are compared against the results presented in case B of [TTLP301]\(^1\).

Timesteps \([s]\)[TTLP301]
\([K]\)
Case A \([K]\)Error [%]Case B \([K]\)Error [%]Case C \([K]\)Error [%]Case D \([K]\)Error [%]
0.1274.195274.2012.19E-03274.2032.92E-03274.1981.09E-03274.2043.28E-03
0.2274.597274.6011.46E-03274.6042.55E-03274.61.09E-03274.6042.55E-03
0.3274.892274.8951.09E-03274.8982.18E-03274.8947.28E-04274.8982.18E-03
0.4275.132275.1351.09E-03275.1382.18E-03275.1347.27E-04275.1382.18E-03
0.5275.339275.3421.09E-03275.3452.18E-03275.3417.26E-04275.3452.18E-03
0.6275.523275.5261.09E-03275.5292.18E-03275.5257.26E-04275.5292.18E-03
0.7275.691275.6951.45E-03275.6972.18E-03275.6941.09E-03275.6972.18E-03
0.8275.848275.8521.45E-03275.8542.18E-03275.8511.09E-03275.8542.18E-03
0.9276.996275.9991.09E-03276.0011.81E-03275.9987.25E-04276.0011.81E-03
1.0276.136276.141.45E-03276.1422.17E-03276.1391.09E-03276.1422.17E-03
1.1276.27276.2741.45E-03276.2762.17E-03276.2731.09E-03276.2762.17E-03
1.2276.398276.4031.81E-03276.4042.17E-03276.4021.45E-03276.4042.17E-03
Table 3: Temperature results for all cases at different time steps

All of the cases are in a good agreement with the reference results.

You can see the temperature distribution on the plate for the last time step of Case A below:

temperature distribution perforated plate
Figure 4: Temperature distribution on the plate for Case A at time 1.2s

Last updated: August 24th, 2020

Contents
Data Privacy