Relativistic heat conduction

Last updated

Relativistic heat conduction refers to the modelling of heat conduction (and similar diffusion processes) in a way compatible with special relativity. In special (and general) relativity, the usual heat equation for non-relativistic heat conduction must be modified, as it leads to faster-than-light signal propagation. [1] [2] Relativistic heat conduction, therefore, encompasses a set of models for heat propagation in continuous media (solids, fluids, gases) that are consistent with relativistic causality, namely the principle that an effect must be within the light-cone associated to its cause. Any reasonable relativistic model for heat conduction must also be stable, in the sense that differences in temperature propagate both slower than light and are damped over time (this stability property is intimately intertwined with relativistic causality [3] ).

Contents

Parabolic model (non-relativistic)

Heat conduction in a Newtonian context is modelled by the Fourier equation, [4] namely a parabolic partial differential equation of the kind: where θ is temperature, [5] t is time, α = k/(ρc) is thermal diffusivity, k is thermal conductivity, ρ is density, and c is specific heat capacity. The Laplace operator, , is defined in Cartesian coordinates as

This Fourier equation can be derived by substituting Fourier’s linear approximation of the heat flux vector, q, as a function of temperature gradient, into the first law of thermodynamics where the del operator, , is defined in 3D as

It can be shown that this definition of the heat flux vector also satisfies the second law of thermodynamics, [6] where s is specific entropy and σ is entropy production. This mathematical model is inconsistent with special relativity: the Green function associated to the heat equation (also known as heat kernel) has support that extends outside the light-cone, leading to faster-than-light propagation of information. For example, consider a pulse of heat at the origin; then according to Fourier equation, it is felt (i.e. temperature changes) at any distant point, instantaneously. The speed of propagation of heat is faster than the speed of light in vacuum, which is inadmissible within the framework of relativity.

Hyperbolic model (relativistic)

The parabolic model for heat conduction discussed above shows that the Fourier equation (and the more general Fick's law of diffusion) is incompatible with the theory of relativity [7] for at least one reason: it admits infinite speed of propagation of the continuum field (in this case: heat, or temperature gradients). To overcome this contradiction, workers such as Carlo Cattaneo, [2] Vernotte, [8] Chester, [9] and others [10] proposed that Fourier equation should be upgraded from the parabolic to a hyperbolic form, where the n, the temperature field is governed by:

In this equation, C is called the speed of second sound (that is related to excitations and quasiparticles, like phonons). The equation is known as the "hyperbolic heat conduction" (HHC) equation. [11] Mathematically, the above equation is called "telegraph equation", as it is formally equivalent to the telegrapher's equations, which can be derived from Maxwell’s equations of electrodynamics.

For the HHC equation to remain compatible with the first law of thermodynamics, it is necessary to modify the definition of heat flux vector, q, to where is a relaxation time, such that This equation for the heat flux is often referred to as "Maxwell-Cattaneo equation". The most important implication of the hyperbolic equation is that by switching from a parabolic (dissipative) to a hyperbolic (includes a conservative term) partial differential equation, there is the possibility of phenomena such as thermal resonance [12] [13] [14] and thermal shock waves. [15]

Notes

  1. van Kampen, N. G. (1970-03-02). "A model for relativistic heat transport". Physica. 46 (2): 315–332. Bibcode:1970Phy....46..315V. doi:10.1016/0031-8914(70)90231-4. ISSN   0031-8914.
  2. 1 2 Cattaneo, C. R. (1958). "Sur une forme de l'équation de la chaleur éliminant le paradoxe d'une propagation instantanée". Comptes Rendus. 247 (4): 431.
  3. Gavassino, L.; Antonelli, M.; Haskell, B. (2022-01-06). "Thermodynamic Stability Implies Causality". Physical Review Letters. 128 (1): 010606. arXiv: 2105.14621 . Bibcode:2022PhRvL.128a0606G. doi:10.1103/PhysRevLett.128.010606. PMID   35061457. S2CID   235254457.
  4. Carslaw, H. S.; Jaeger, J. C. (1959). Conduction of Heat in Solids (Second ed.). Oxford: University Press.
  5. Some authors also use T, φ,...
  6. Barletta, A.; Zanchini, E. (1997). "Hyperbolic heat conduction and local equilibrium: a second law analysis". International Journal of Heat and Mass Transfer. 40 (5): 1007–1016. doi:10.1016/0017-9310(96)00211-6.
  7. Eckert, E. R. G.; Drake, R. M. (1972). Analysis of Heat and Mass Transfer. Tokyo: McGraw-Hill, Kogakusha.
  8. Vernotte, P. (1958). "Les paradoxes de la theorie continue de l'équation de la chaleur". Comptes Rendus. 246 (22): 3154.
  9. Chester, M. (1963). "Second sound in solids". Physical Review . 131 (15): 2013–2015. Bibcode:1963PhRv..131.2013C. doi:10.1103/PhysRev.131.2013.
  10. Morse, P. M.; Feshbach, H. (1953). Methods of Theoretical Physics. New York: McGraw-Hill.
  11. Joseph, D. D.; Preziosi, L. (1989). "Heat waves". Reviews of Modern Physics. 61 (1): 47–71. Bibcode:1989RvMP...61...41J. doi:10.1103/RevModPhys.61.41.
  12. Mandrusiak, G. D. (1997). "Analysis of non-Fourier conduction waves from a reciprocating heat source". Journal of Thermophysics and Heat Transfer. 11 (1): 82–89. doi:10.2514/2.6204.
  13. Xu, M.; Wang, L. (2002). "Thermal oscillation and resonance in dual-phase-lagging heat conduction". International Journal of Heat and Mass Transfer. 45 (5): 1055–1061. doi:10.1016/S0017-9310(01)00199-5.
  14. Barletta, A.; Zanchini, E. (1996). "Hyperbolic heat conduction and thermal resonances in a cylindrical solid carrying a steady periodic electric field". International Journal of Heat and Mass Transfer. 39 (6): 1307–1315. doi:10.1016/0017-9310(95)00202-2.
  15. Tzou, D. Y. (1989). "Shock wave formation around a moving heat source in a solid with finite speed of heat propagation". International Journal of Heat and Mass Transfer. 32 (10): 1979–1987. doi:10.1016/0017-9310(89)90166-X.

Related Research Articles

<span class="mw-page-title-main">Laplace's equation</span> Second-order partial differential equation

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace, who first studied its properties. This is often written as or where is the Laplace operator, is the divergence operator, is the gradient operator, and is a twice-differentiable real-valued function. The Laplace operator therefore maps a scalar function to another scalar function.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

The thermal conductivity of a material is a measure of its ability to conduct heat. It is commonly denoted by , , or and is measured in W·m−1·K−1.

<span class="mw-page-title-main">Schrödinger equation</span> Description of a quantum-mechanical system

The Schrödinger equation is a partial differential equation that governs the wave function of a non-relativistic quantum-mechanical system. Its discovery was a significant landmark in the development of quantum mechanics. It is named after Erwin Schrödinger, who postulated the equation in 1925 and published it in 1926, forming the basis for the work that resulted in his Nobel Prize in Physics in 1933.

Thermal conduction is the diffusion of thermal energy (heat) within one material or between materials in contact. The higher temperature object has molecules with more kinetic energy; collisions between molecules distributes this kinetic energy until an object has the same kinetic energy throughout. Thermal conductivity, frequently represented by k, is a property that relates the rate of heat loss per unit area of a material to its rate of change of temperature. Essentially, it is a value that accounts for any property of the material that could change the way it conducts heat. Heat spontaneously flows along a temperature gradient. For example, heat is conducted from the hotplate of an electric stove to the bottom of a saucepan in contact with it. In the absence of an opposing external driving energy source, within a body or between bodies, temperature differences decay over time, and thermal equilibrium is approached, temperature becoming more uniform.

<span class="mw-page-title-main">Heat equation</span> Partial differential equation describing the evolution of temperature in a region

In mathematics and physics, the heat equation is a certain partial differential equation. Solutions of the heat equation are sometimes known as caloric functions. The theory of the heat equation was first developed by Joseph Fourier in 1822 for the purpose of modeling how a quantity such as heat diffuses through a given region.

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols , (where is the nabla operator), or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p).

The Klein–Gordon equation is a relativistic wave equation, related to the Schrödinger equation. It is second-order in space and time and manifestly Lorentz-covariant. It is a differential equation version of the relativistic energy–momentum relation .

<span class="mw-page-title-main">Four-vector</span> 4-dimensional vector in relativity

In special relativity, a four-vector is an object with four components, which transform in a specific way under Lorentz transformations. Specifically, a four-vector is an element of a four-dimensional vector space considered as a representation space of the standard representation of the Lorentz group, the representation. It differs from a Euclidean vector in how its magnitude is determined. The transformations that preserve this magnitude are the Lorentz transformations, which include spatial rotations and boosts.

Linear elasticity is a mathematical model as to how solid objects deform and become internally stressed by prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

A continuity equation or transport equation is an equation that describes the transport of some quantity. It is particularly simple and powerful when applied to a conserved quantity, but it can be generalized to apply to any extensive quantity. Since mass, energy, momentum, electric charge and other natural quantities are conserved under their respective appropriate conditions, a variety of physical phenomena may be described using continuity equations.

<span class="mw-page-title-main">Onsager reciprocal relations</span> Relations between flows and forces, or gradients, in thermodynamic systems

In thermodynamics, the Onsager reciprocal relations express the equality of certain ratios between flows and forces in thermodynamic systems out of equilibrium, but where a notion of local equilibrium exists.

In mathematics, the Helmholtz equation is the eigenvalue problem for the Laplace operator. It corresponds to the elliptic partial differential equation: where 2 is the Laplace operator, k2 is the eigenvalue, and f is the (eigen)function. When the equation is applied to waves, k is known as the wave number. The Helmholtz equation has a variety of applications in physics and other sciences, including the wave equation, the diffusion equation, and the Schrödinger equation for a free particle.

In fluid mechanics, potential vorticity (PV) is a quantity which is proportional to the dot product of vorticity and stratification. This quantity, following a parcel of air or water, can only be changed by diabatic or frictional processes. It is a useful concept for understanding the generation of vorticity in cyclogenesis, especially along the polar front, and in analyzing flow in the ocean.

The Richards equation represents the movement of water in unsaturated soils, and is attributed to Lorenzo A. Richards who published the equation in 1931. It is a quasilinear partial differential equation; its analytical solution is often limited to specific initial and boundary conditions. Proof of the existence and uniqueness of solution was given only in 1983 by Alt and Luckhaus. The equation is based on Darcy-Buckingham law representing flow in porous media under variably saturated conditions, which is stated as

<span class="mw-page-title-main">HydroGeoSphere</span>

HydroGeoSphere (HGS) is a 3D control-volume finite element groundwater model, and is based on a rigorous conceptualization of the hydrologic system consisting of surface and subsurface flow regimes. The model is designed to take into account all key components of the hydrologic cycle. For each time step, the model solves surface and subsurface flow, solute and energy transport equations simultaneously, and provides a complete water and solute balance.

<span class="mw-page-title-main">Heat flux</span> Vector representing the energy passing through a given area per unit time

In physics and engineering, heat flux or thermal flux, sometimes also referred to as heat flux density, heat-flow density or heat-flow rate intensity, is a flow of energy per unit area per unit time. Its SI units are watts per square metre (W/m2). It has both a direction and a magnitude, and so it is a vector quantity. To define the heat flux at a certain point in space, one takes the limiting case where the size of the surface becomes infinitesimally small.

In fluid dynamics, the Oseen equations describe the flow of a viscous and incompressible fluid at small Reynolds numbers, as formulated by Carl Wilhelm Oseen in 1910. Oseen flow is an improved description of these flows, as compared to Stokes flow, with the (partial) inclusion of convective acceleration.

Miniaturizing components has always been a primary goal in the semiconductor industry because it cuts production cost and lets companies build smaller computers and other devices. Miniaturization, however, has increased dissipated power per unit area and made it a key limiting factor in integrated circuit performance. Temperature increase becomes relevant for relatively small-cross-sections wires, where it may affect normal semiconductor behavior. Besides, since the generation of heat is proportional to the frequency of operation for switching circuits, fast computers have larger heat generation than slow ones, an undesired effect for chips manufacturers. This article summaries physical concepts that describe the generation and conduction of heat in an integrated circuit, and presents numerical methods that model heat transfer from a macroscopic point of view.

Conservative temperature is a thermodynamic property of seawater. It is derived from the potential enthalpy and is recommended under the TEOS-10 standard as a replacement for potential temperature as it more accurately represents the heat content in the ocean.