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.