14.55. PLANE55 - 2-D Thermal Solid
| Matrix or Vector | Geometry | Shape Functions | Integration Points |
|---|---|---|---|
| Conductivity Matrix and Heat Generation Load Vector | Quad | (Equation 12–108) | 2 x 2 |
| Triangle | (Equation 12–89) | 1 if planar | |
| Specific Heat Matrix | Same as conductivity matrix. Matrix is diagonalized as described in Lumped Matrices. | Same as conductivity matrix | |
| Convection Surface Matrix and Load Vector | Same as conductivity matrix evaluated at the face | 2 | |
14.55.1. Other Applicable Sections
"Heat Flow" describes the derivation of the element matrices and load vectors as well as heat flux evaluations. SOLID70 - 3-D Thermal Solid describes fluid flow in a porous medium, accessed in PLANE55 with KEYOPT(9) = 1.
14.55.2. Mass Transport Option
If KEYOPT(8) > 0, the mass transport option is included as described
in Heat Flow Fundamentals with (Equation 6–1)
and by 
of (Equation 6–21).
The solution accuracy is dependent on the element size. The accuracy is
measured in terms of the non-dimensional criteria called the element Peclet
number (Gresho(58)):
| (14–382) |
where:
| V = magnitude of the velocity vector |
| L = element length dimension along the velocity vector direction |
| ρ = density of the fluid (input as DENS on MP command) |
| Cp = specific heat of the fluid (input as C on MP command) |
| K = equivalent thermal conductivity along the velocity vector direction |
The terms V, L, and K are explained more thoroughly below:
| (14–383) |
where:
| Vx = fluid velocity (mass transport) in x direction (input as VX on R command) |
| Vy = fluid velocity (mass transport) in y direction (input as VY on R command) |
Length L is calculated by finding the intersection points of the velocity vector which passes through the element origin and intersects at the element boundaries.
For orthotropic materials, the equivalent thermal conductivity K is given by:
| (14–384) |
where:
| Kx, Ky = thermal conductivities in the x and y directions (input as KXX and KYY on MP command) |
 |
For the solution to be physically valid, the following condition has to be satisfied (Gresho(58)):
| (14–385) |
This check is carried out during the element formulation and an error message is printed out if equation (14.431) is not satisfied. When this error occurs, the problem should be rerun after reducing the element size in the direction of the velocity vector.