Finite element method

The finite element method computes the temperature values at discrete points, solid nodes, by solving the heat conduction equation.

The three-dimensional governing equation for the heat conduction in a solid domain Ω [1, 2] is given as:

where:

  • ρ is the material density.
  • c is the specific heat capacity.
  • k is the thermal conductivity, where subscripts x, y, and z are the anisotropic specifications in space directions.
  • qv is the heat generation per unit volume.

The thermal solver calculates temperature distribution using the assigned initial condition and Dirichlet and Neumann boundary conditions.

The required initial condition is given as:

The Dirichlet boundary condition that applies a constant temperature TD at the boundary ΩD is given as:

The Neumann boundary condition that applies a constant heat flux qN at the boundary region ΩN is given as:

where n is the unit normal vector directed to the exterior of the boundary surface.

The example shows both of the assigned boundary conditions for conduction heat transfer on a solid domain.