href="#fb3_img_img_49fa2839-833b-55bb-a3d9-39882f7cf217.png" alt="images"/> uniquely, the differential equation is complemented by boundary conditions imposed at the boundary points and . A common example is the homogeneous Dirichlet conditions, corresponding to keeping the temperature at zero at the endpoints of the wire. The result is a two‐point BVP:
(1.5.7)
The boundary condition may be replaced by , corresponding to prescribing zero heat flux, or insulating the wire, at . Later, we also consider nonhomogeneous boundary conditions of the form or where and may be different from zero. For other types of boundary conditions, see Trinities (Section 1.2).
The time‐dependent heat equation in (1.5.2) describes the diffusion of thermal energy in a homogeneous material, where is the temperature at a position at time and is called thermal diffusivity or heat conductivity (corresponding to in (1.5.5)–(1.5.7)) of the material.
Remark 1.3
The heat equation can be used to model the heat flow in solids and fluids, in the latter case, however, it does not take into account the convection phenomenon; and provides a reasonable model only if phenomena such as macroscopic currents in the fluid are not present (or negligible). Further, the heat equation is not a fundamental law of physics, and it does not give reliable answers at very low or very high temperatures.
Since temperature is related to heat, which is a form of energy, the basic idea in deriving the heat equation is to use the law of conservation of energy. Below we derive the general form of the heat equation in arbitrary dimension.
1.5.2.2 Fourier's Law of Heat Conduction, Derivation of the Heat Equation
Let be a fixed spatial domain with boundary . The rate of change of thermal energy with respect to time in is equal to the net flow of energy across the boundary of plus the rate at which heat is generated within .
Let denote the temperature at the position and at time . We assume that the solid is at rest and that it is rigid, so that the only energy present is thermal energy and the density is independent of the time and temperature . Let denote the energy per unit mass. Then the amount of thermal energy in is given by
and the time rate (time derivative) of change of thermal energy in is:
Let denote the heat flux vector and denote the outward unit normal to the boundary , at the point . Then represents the flow of heat per unit cross‐sectional area per unit time crossing a surface element. Thus,
is the amount of heat per unit time flowing into across the boundary . Here, represents the element of surface area. The minus sign reflects the fact that if more heat flows out of the domain than in, the energy in decreases. Finally, in general, the heat production is determined by external sources that are independent of the temperature. In some cases, (such as an air conditioner controlled by a thermostat), it depends on temperature itself, but not on its derivatives. Hence, in the presence of a source (or sink), we denote the corresponding rate at which heat is produced per unit volume by so that the source term becomes
Now, the law of conservation of energy takes the form
(1.5.8)
Applying the Gauss divergence theorem to the integral over
