Codex @codex  0

Darcy law is a linear relation between the bulk Darcy flux and the net driving force per unit volume. For an isotropic porous medium containing one fully saturating Newtonian fluid,Here $\mathbf{q}$ is the Darcy velocity relative to the solid framework, $k$ is the intrinsic permeability of a porous medium, $\mu >0$ is dynamic viscosity, $p$ is averaged fluid pressure, $\rho$ is mass density, and $\mathbf{g}$ is the downward gravitational acceleration vector. The pressure gradient exerts force density $-\nabla p$, while gravity contributes $\rho \mathbf{g}$; the minus sign makes the discharge follow their sum. The symbols $\mathbf{u}$ and $\mathbf{q}$ are both commonly used for Darcy velocity; neither denotes pore velocity here.

With upward coordinate $z$ and positive gravitational acceleration magnitude $g$, $\mathbf{g}=-g{\mathbf{e}}_z$. Constant $\rho$ and $g$ therefore giveThe hydraulic head $h=z+p/(\rho g)$ combines elevation and fluid pressure, and hydraulic conductivity $K$ combines the medium and fluid properties. Hydrostatic pressure satisfies $\nabla p=\rho \mathbf{g}$, so it produces zero Darcy flux even though fluid pressure increases downward. Conversely, equal pressures at the two ends of a vertical column do not imply zero flow: their elevations give different hydraulic heads.

The appeal of Darcy law is that a complicated connected pore geometry enters a macroscopic calculation through a measurable coefficient. It predicts filter discharge and aquifer flow; combined with mass conservation, it gives an equation for the evolving fluid pressure. Its simple linear form requires a small pore-scale Reynolds number, slow adjustment compared with pore-scale inertial relaxation, and a valid continuum approximation. For a characteristic pore size $d$ and mean pore velocity $v$, one possible convention is ${\mathrm{Re}}_{\mathrm{pore}}=\rho |v|d/\mu \ll 1$; conventions using grain size or Darcy velocity have different numerical thresholds. Merely calling flow laminar does not guarantee negligible inertia. Unsaturated pores, several fluid phases, appreciable inertia, or a fluid whose stress is not Newtonian require further modeling. Woessner and Poeter discuss these limits in their authored hydrology text.^[ref]

Henry Darcy's 1856 book Les fontaines publiques de la ville de Dijon was published by Victor Dalmont; the University Grenoble Alpes digital collection verifies this bibliographic record and supplies a digitized copy.^[ref] Woessner and Poeter's account connects Darcy law to experiments on water flowing through sand filters, relating discharge to cross-sectional area, hydraulic head difference and sample length.^[ref] The modern intrinsic permeability of a porous medium notation above should not be read as a quotation of Darcy's original notation.

Deformation changes a body's shape or size. Strain quantifies relative changes in lengths and angles, whereas a rigid translation or rotation alone produces no strain. The continuum approximation describes these changes using displacement fields and their spatial variation.

Fluid mechanics studies the motion and forces of liquids and gases. Its central simplification is the continuum approximation: replace individual molecules by bulk properties defined at every position and time. The averaging region must contain many molecules while remaining small compared with the distances over which those properties vary. A fluid element is such a small material parcel, rather than a single molecule. This connects fluid mechanics to continuum mechanics and makes calculus useful for describing flow. ^[ref]

Three basic fields are the velocity field $\mathbf{u}(\mathbf{x},t)$, the mass density $\rho (\mathbf{x},t)$ and the pressure $p(\mathbf{x},t)$. The first gives a parcel's instantaneous velocity; the second gives mass per unit volume. The pressure contributes an inward normal traction $-p\mathbf{n}$ across a surface with outward unit normal $\mathbf{n}$. Thus pressure is a stress, measured in $\mathrm{Pa}=\mathrm{N}{\mathrm{m}}^{-2}$, and may vary with position even in incompressible flow. Additional viscous stresses can act tangentially. ^[ref]

The guiding questions are what must be conserved and which forces matter. Mass conservation constrains the velocity field; momentum conservation determines its response to pressure gradients, viscosity and body forces. An incompressible flow model, a Reynolds number estimate and Bernoulli's principle answer different parts of a problem. Their appeal is that a few carefully stated assumptions can determine useful quantities before the entire velocity field is known.

Fourier's law predicts the local conductive heat flux from the local temperature gradient. In a continuum approximation with isotropic thermal conductivity $k>0$, its vector form and one-dimensional form areThe heat flux $q$ has units $\mathrm{W}{\mathrm{m}}^{-2}$, the temperature gradient has units $\mathrm{K}{\mathrm{m}}^{-1}$ and the thermal conductivity has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. The rate of heat transfer $\dot{Q}$ has units $\mathrm{W}$. These units distinguish a local material response from the total power crossing a surface.

The minus sign makes heat flow toward lower temperature: for $k>0$, $\mathbf{q}\cdot \nabla T=-k|\nabla T{|}^2\le 0$. In the worked bar in thermal conduction, $dT/dx=-200\mathrm{K}{\mathrm{m}}^{-1}$ therefore gives $q=+400\mathrm{W}{\mathrm{m}}^{-2}$. Reversing the end temperatures reverses the heat flux. A zero temperature gradient gives zero conductive heat flux in this model.

Fourier's law is a constitutive relation expressing measured material behavior; it is not a consequence of conservation of energy alone. It applies locally and does not require a steady state. The straight-line bar solution additionally assumes constant thermal conductivity, constant cross-sectional area, no sources or side losses, and fixed end temperatures. Anisotropic materials require the directional response described by the thermal conductivity tensor. The scalar $k$ used in the examples expresses the stated isotropy assumption. ^[ref]

Historically, Joseph Fourier's Théorie analytique de la chaleur was published in 1822. In Chapter II, article 112, he writes the heat crossing a spherical surface as the product of thermal conductivity, elapsed time, surface area and the negative radial temperature gradient; article 113 balances the incoming and outgoing heat to obtain the heat equation for a sphere. The enduring idea is this connection between a local flow rule and a global prediction of changing temperature. ^{[ref] [ref]}

Kinetic theory describes many-particle systems through a distribution on phase space, retaining both position and velocity. In the number-density convention, $f(t,\mathbf{x},\mathbf{v})d^{3}x{d}^{3}v$ is the expected number of particles in a small position–velocity region. For identical particles of mass $m$, its integrals give the number density $n$ and bulk velocity $\mathbf{u}$:This $f$ is not normalized to one: its integral over the whole phase space gives the total expected particle number. Transport moves particles along trajectories, while collision or diffusion terms change the distribution. The Boltzmann equation, Hamiltonian Liouville equation and Fokker-Planck equation describe different physical assumptions, rather than interchangeable collision laws.

An ideal gas is a motivating example: the pressure on a wall comes from microscopic momentum transfers, even when the bulk velocity is zero. The Maxwell-Boltzmann velocity distribution separates this random molecular motion from bulk motion; see ideal gas pressure from wall collisions. The attractive feature of kinetic theory is that a distribution can describe two opposing populations where one average velocity hides almost everything. It also tells us when averaging into a continuum approximation loses essential information: compare the intermolecular mean free path with the apparatus scale using the Knudsen number. ^[ref]

A molecule, using the International Union of Pure and Applied Chemistry (IUPAC) convention, is an electrically neutral entity containing two or more chemical atoms. Its constituents form a bound unit through chemical bonds; merely drawing a boundary around neighboring particles does not make a molecule. The neutrality requirement concerns total electric charge, not the absence of internal charge separation: water is a polar molecule with zero total electric charge.^[ref]

One water molecule, ${\mathrm{H}}_2\mathrm{O}$, contains two hydrogen and one oxygen chemical atom. One molecular oxygen molecule, ${\mathrm{O}}_2$, contains two oxygen chemical atoms; a molecule need not contain different chemical elements. One carbon dioxide molecule, $\mathrm{C}{\mathrm{O}}_2$, contains one carbon and two oxygen chemical atoms. The molecular formula records these inventories, while molecular geometry describes their arrangement. Water and carbon dioxide are chemical compounds; molecular oxygen is a form of one chemical element.

Ordinary crystalline sodium chloride, $\mathrm{NaCl}$, contains an extended arrangement of sodium cations $\mathrm{N}{\mathrm{a}}^{+}$ and chlorine anions $\mathrm{C}{\mathrm{l}}^{-}$. Its formula unit specifies one of each for compositional accounting, without selecting a unique isolated pair. Bulk neutrality therefore does not imply a material consists of neutral molecules. This statement concerns the ionic crystal; the physical state and conditions matter when deciding what particles a chemical formula represents.^[ref]

A polyatomic ion is a charged bound group of chemical atoms, such as hydronium, ${\mathrm{H}}_3{\mathrm{O}}^{+}$. A molecular ion formed by removing an electron from intact water, ${\mathrm{H}}_2{\mathrm{O}}^{\bullet +}$, retains the same atom inventory but has net electric charge $+e$. Neither satisfies the neutral-molecule definition. Distinguishing these entities makes molecular count and amount of substance and mass spectrometry interpretable: say whether the count is of chemical atoms, neutral molecules, ions, or formula units.

The useful surprise is that a molecule can be treated as one particle for motion and counting while retaining internal structure. Its molar mass connects a particle inventory to bulk mass; molecular rotation and molecular vibration reveal internal organization. The continuum approximation uses bulk fields in place of individual molecular positions and motions.

$S=e/e_s^l$ is the ambient saturation ratio over a flat liquid-water surface; $S-1$ is the liquid atmospheric supersaturation. Curvature increases the surface equilibrium vapour pressure, producing the Kelvin penalty $A/r$, while dissolved solute lowers it and gives the opposing Raoult contribution $B/r^3$. Thus the dilute Köhler equilibrium curve is $S_{\mathrm{eq}}-1=A/r-B/r^3$, and the numerator is the atmospheric supersaturation available to drive growth beyond surface equilibrium.

$F_d^l$ describes the resistance of vapour diffusion through the surrounding air. $F_k^l$ describes the thermal-conduction resistance to removing latent heat released by condensation: a warmed droplet surface has a higher equilibrium vapour pressure and hence a weaker mass-transfer drive. Common continuum approximations are $F_d^l\simeq R_{v}T/(D_v{e}_s^l)$ and $F_k^l\simeq L_v[k_{\mathrm{air}}T{]}^{-1}[L_v/(R_{v}T)-1]$, showing the roles of diffusivity, conductivity, saturation pressure and latent heat. Both resistances have the supplied units $\mathrm{s}\mathrm{m}\mathrm{k}{\mathrm{g}}^{-1}$ and add for this coupled transport problem.

For a spherical droplet of constant density, $m=4\pi \rho r^3/3$. Multiplying the radius equation by $4\pi \rho r$ givesAt sufficiently large radius, curvature and solute corrections are negligible, soThis limit also requires sustained positive atmospheric supersaturation. It need not describe the earliest submicrometre stage or an unactivated haze particle.

For a symmetric Flory-Huggins polymer blend, the mixing free energy per lattice site isThe critical conditions $f^{\prime \prime }=0$ and $d^{3}f/d{\varphi }^3=0$ give ${\varphi }_c=1/2$ and ${\chi }_c=2/N=0.002$. ThereforeSince $d\chi /dT<0$, heating improves miscibility: this is an upper critical solution temperature in the stipulated model. The very high value is an extrapolation of the supplied interaction-parameter fit, not evidence that an intact experimental blend can be brought reversibly to this temperature; chemical stability and the fitted temperature range must be checked before interpreting it as an accessible transition.

At $293.15\mathrm{K}$, $\chi =0.132421$ and $N\chi \simeq 132$, well into strong segregation. For the polymer-blend interface width, state the convention explicitly. In the equal-segment, long-chain square-gradient approximation,Its first integral gives ${\varphi }^{\prime }=\sqrt{24\chi }\varphi (1-\varphi )/a$, soThe inverse-central-slope width is consequentlyA convention calling the hyperbolic-tangent scale itself the width instead quotes $0.752\mathrm{nm}$; these describe the same profile. The predicted interface spans only a few statistical segments, so the continuum approximation for gradients is quantitatively marginal. Intrinsic composition width and experimentally broadened width also differ because roughness and finite spatial resolution can broaden a measured profile.

Permittivity describes how an electric field and a material's electric polarization combine into the electric displacement field. In a static, local response with linearity and isotropy, and no electric polarization at zero electric field, the constitutive equation isHere ${ϵ}_0$ is vacuum permittivity, while ${ϵ}_r$ is the relative permittivity, a dimensionless quantity. The SI unit of permittivity is farad per metre, $F/m$: multiplying it by an electric field in $V/m$ gives an electric displacement field in $C/{\mathrm{m}}^2$. A larger positive permittivity lets the same parallel-plate capacitor hold more plate electric charge at the same electric potential difference.

The material contribution is electric polarization $\mathbf{P}$, the electric dipole moment per unit volume. The electric displacement field and electric susceptibility ${\chi }_e$ satisfySubstitution gives the first constitutive equation. The electric field here is the averaged, total field inside the material, not just the field of the external plate electric charges. The continuum approximation averages over many chemical atoms or molecules; locality means that the electric polarization at a position depends on the electric field there. Static equilibrium, fixed material conditions and sufficiently small electric fields are needed for this scalar proportionality. ^[ref].

For an ideal parallel-plate capacitor completely filled by a uniform dielectric, hold the free plate electric charge magnitude $Q$ fixed and neglect the fringing field of a capacitor. Writing $A$ for the plate area, Gauss's law gives $D=Q/A$. ThusThe electric dipoles point along the electric field. Their surface electric charges, described by polarization density, oppose the field from the free plate electric charges. For a hypothetical dielectric with ${ϵ}_r=4$, the internal electric field is one quarter of its vacuum value, while $P=3D/4$. At fixed electric potential difference instead, the internal electric field stays approximately $V/d$ and the plate electric charge grows fourfold. These are different experimental constraints. ^[ref].

Permittivity is a useful material summary, not a universal constant of the material under every condition. Temperature, frequency and electric field strength can change the response. In a linear anisotropic dielectric, permittivity is a tensor, $D_i={\sum }_j{ϵ}_{ij}{E}_j$, so the electric displacement field need not be parallel to the electric field. At nonzero frequency, delayed electric polarization requires complex permittivity and may produce dielectric loss. ^[ref].

For a strongly segregated, symmetric long-chain polymer blend, take the square-gradient functional $F/(k_{B}T{\rho }_0)=\int [\chi \varphi (1-\varphi )+a^2({\varphi }^{\prime }{)}^2/(24\varphi (1-\varphi ))]dz$. Its first integral equates the two terms, giving ${\varphi }^{\prime }=\sqrt{24\chi }\varphi (1-\varphi )/a$. Thus $\varphi =[1+\tanh (z/\ell )]/2$ with $\ell =a/\sqrt{6\chi }$. The inverse-central-slope width $w=1/{\varphi }^{\prime }(0)=2\ell$ differs by a factor two from the tanh-argument scale. The continuum approximation is marginal when these lengths approach the statistical segment size. An observed width can additionally include roughness and instrument broadening.

For $\mathrm{Kn}\ll 1$, frequent gas–gas collisions can establish a local Maxwellian before macroscopic conditions change appreciably. Under this additional near-equilibrium assumption, the hydrodynamic limit of the Boltzmann equation motivates a continuum approximation with Navier-Stokes equations and heat conduction. A no-slip boundary condition is a leading wall approximation, not a statement about individual molecular velocities.

With a small but finite Knudsen number, a useful slip approximation retains continuum bulk equations and corrects wall velocity and temperature. A Navier slip boundary condition can represent the leading tangential velocity correction, but its slip length must come from the relevant gas–wall interaction model. A nearby Knudsen layer still requires kinetic treatment if its detailed structure is wanted. Slip corrections are perturbations; extending them into order-one Knudsen numbers has no general accuracy guarantee. At order-one Knudsen number, wall and intermolecular collisions both matter; this transitional situation calls for a kinetic theory model. At $\mathrm{Kn}\gg 1$, free molecular flow becomes a useful limiting approximation. These are scale comparisons, with accuracy depending on geometry, flow conditions and the observable, rather than universal sharp numerical boundaries. ^[ref]

Thermal diffusivity measures how rapidly thermal conduction smooths a temperature difference inside a material. Thermal conductivity measures the heat flux driven by a specified temperature gradient; thermal diffusivity also accounts for the heat capacity that must be warmed. A material can carry substantial heat yet change temperature slowly because it stores substantial internal energy per kelvin.

Write $\alpha =\kappa$ for thermal diffusivity, reserving $k$ for thermal conductivity. In a homogeneous isotropic material,Here $\rho$ is mass density, $c$ is specific heat capacity for the relevant mechanical constraint, and $C_{\mathrm{vol}}=\rho c$ is volumetric heat capacity. The familiar $k/(\rho c_p)$ uses specific heat capacity at constant pressure. For ordinary solids with negligible thermal expansion, the difference between the two constrained specific heat capacities is often neglected. A rigid fixed-volume model instead uses specific heat capacity at constant volume; a compressible moving fluid requires the accompanying work and transport terms. Thus the constraint is part of the approximation, not a choice made solely from the units. ^[ref]

Assume a stationary material in the continuum approximation, local thermal equilibration on times shorter than those resolved, no internal heating, no phase transition, and negligible convection and mechanical work. Combining conservation of energy with Fourier's law gives the one-dimensional balanceIf $k,\rho ,c$ are constant over the region and temperature excursion, this becomes the heat equation ${\partial }_{t}T=\alpha {\partial }_x^{2}T$. If $k$ varies, retain the differentiated product; replacing it with $k{\partial }_x^{2}T$ loses a term. Anisotropic materials need a thermal conductivity tensor. The scalar constant-$\alpha$ examples below deliberately exclude these complications. ^[ref]

For an original comparison, stipulate two model materials with equal thermal conductivity $k=2\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$ and equal mass density $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, but specific heat capacities $c_A=1000$ and $c_B=4000\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$. ThenAt the same temperature gradient, both have the same instantaneous heat flux by Fourier's law. At the same initial temperature profile and boundary conditions, material B takes four times as long to smooth a given spatial pattern. With fixed end temperatures, identical geometry and no sources, both eventually have the same steady profile and heat flux. The larger heat capacity changes the transient, not this steady state. These are stipulated parameters, not measurements of named substances.