Английская Википедия:Bulk modulus
Шаблон:Short description Шаблон:Confusing
The bulk modulus (<math>K</math> or <math>B</math> or <math>k</math>) of a substance is a measure of the resistance of a substance to bulk compression. It is defined as the ratio of the infinitesimal pressure increase to the resulting relative decrease of the volume.[1]
Other moduli describe the material's response (strain) to other kinds of stress: the shear modulus describes the response to shear stress, and Young's modulus describes the response to normal (lengthwise stretching) stress. For a fluid, only the bulk modulus is meaningful. For a complex anisotropic solid such as wood or paper, these three moduli do not contain enough information to describe its behaviour, and one must use the full generalized Hooke's law. The reciprocal of the bulk modulus at fixed temperature is called the isothermal compressibility.
Definition
The bulk modulus <math>K</math> (which is usually positive) can be formally defined by the equation
- <math>K=-V\frac{dP}{dV} ,</math>
where <math>P</math> is pressure, <math>V</math> is the initial volume of the substance, and <math>dP/dV</math> denotes the derivative of pressure with respect to volume. Since the volume is inversely proportional to the density, it follows that
- <math>K=\rho \frac{dP}{d\rho} ,</math>
where <math>\rho</math> is the initial density and <math>dP/d\rho</math> denotes the derivative of pressure with respect to density. The inverse of the bulk modulus gives a substance's compressibility. Generally the bulk modulus is defined at constant temperature as the isothermal bulk modulus, but can also be defined at constant entropy as the adiabatic bulk modulus.
Thermodynamic relation
Strictly speaking, the bulk modulus is a thermodynamic quantity, and in order to specify a bulk modulus it is necessary to specify how the pressure varies during compression: constant-temperature (isothermal <math>K_T</math>), constant-entropy (isentropic <math>K_S</math>), and other variations are possible. Such distinctions are especially relevant for gases.
For an ideal gas, an isentropic process has:
- <math>PV^\gamma=\text{constant} \Rightarrow P\propto \left(\frac{1}{V}\right)^\gamma\propto \rho ^\gamma,
</math>
where <math>\gamma </math> is the heat capacity ratio. Therefore, the isentropic bulk modulus <math>K_S</math> is given by
- <math>K_S=\gamma P.</math>
Similarly, an isothermal process of an ideal gas has:
- <math>PV=\text{constant} \Rightarrow P\propto \frac{1}{V} \propto \rho,
</math>
Therefore, the isothermal bulk modulus <math>K_T</math> is given by
- <math>K_T = P </math> .
When the gas is not ideal, these equations give only an approximation of the bulk modulus. In a fluid, the bulk modulus <math>K</math> and the density <math>\rho</math> determine the speed of sound <math>c</math> (pressure waves), according to the Newton-Laplace formula
- <math>c=\sqrt{\frac{K}{\rho}}.</math>
In solids, <math>K_S</math> and <math>K_T</math> have very similar values. Solids can also sustain transverse waves: for these materials one additional elastic modulus, for example the shear modulus, is needed to determine wave speeds.
Measurement
It is possible to measure the bulk modulus using powder diffraction under applied pressure. It is a property of a fluid which shows its ability to change its volume under its pressure.
Selected values
Material | Bulk modulus in GPa | Bulk modulus in Mpsi |
---|---|---|
Diamond (at 4K) [2] | Шаблон:Val | Шаблон:Val |
Alumina (γ phase)[3] | Шаблон:Val ± 14 | Шаблон:Val |
Steel | Шаблон:Val | Шаблон:Val |
Limestone | Шаблон:Val | Шаблон:Val |
Granite | Шаблон:Val | Шаблон:Val |
Glass (see also diagram below table) | Шаблон:Val to Шаблон:Val | Шаблон:Val |
Graphite 2H (single crystal)[4] | Шаблон:Val | Шаблон:Val |
Sodium chloride | Шаблон:Val | Шаблон:Val |
Shale | Шаблон:Val | Шаблон:Val |
Chalk | Шаблон:Val | Шаблон:Val |
Rubber[5] | Шаблон:Val to Шаблон:Val | Шаблон:Val to Шаблон:Val |
Sandstone | Шаблон:Val | Шаблон:Val |
A material with a bulk modulus of 35 GPa loses one percent of its volume when subjected to an external pressure of 0.35 GPa (~Шаблон:Val) (assumed constant or weakly pressure dependent bulk modulus).
β-Carbon nitride | Шаблон:Val[7] (predicted) |
Water | Шаблон:Val (Шаблон:Val) (value increases at higher pressures) |
Methanol | Шаблон:Val (at 20 °C and 1 Atm) |
Solid helium | Шаблон:Val (approximate) |
Air | Шаблон:Val (adiabatic bulk modulus [or isentropic bulk modulus]) |
Air | Шаблон:Val (isothermal bulk modulus) |
Universe (space-time) | Шаблон:Val (for typical gravitational wave frequencies of 100Hz) [8] |
Microscopic origin
Interatomic potential and linear elasticity
Since linear elasticity is a direct result of interatomic interaction, it is related to the extension/compression of bonds. It can then be derived from the interatomic potential for crystalline materials.[9] First, let us examine the potential energy of two interacting atoms. Starting from very far points, they will feel an attraction towards each other. As they approach each other, their potential energy will decrease. On the other hand, when two atoms are very close to each other, their total energy will be very high due to repulsive interaction. Together, these potentials guarantee an interatomic distance that achieves a minimal energy state. This occurs at some distance a0, where the total force is zero:
- <math>F=-{\partial U \over \partial r}=0</math>
Where U is interatomic potential and r is the interatomic distance. This means the atoms are in equilibrium.
To extend the two atoms approach into solid, consider a simple model, say, a 1-D array of one element with interatomic distance of a, and the equilibrium distance is a0. Its potential energy-interatomic distance relationship has similar form as the two atoms case, which reaches minimal at a0, The Taylor expansion for this is:
- <math>u(a)=u(a_0)+ \left({\partial u \over \partial r} \right )_{r=a_0}(a-a_0)+{1 \over 2} \left ({\partial^2\over\partial r^2}u \right )_{r=a_0}(a-a_0)^2+O \left ((a-a_0)^3 \right )</math>
At equilibrium, the first derivative is 0, so the dominant term is the quadratic one. When displacement is small, the higher order terms should be omitted. The expression becomes:
- <math>u(a)=u(a_0)+{1 \over 2} \left ({\partial^2\over\partial r^2}u \right )_{r=a_0}(a-a_0)^2</math>
- <math>F(a)=-{\partial u \over \partial r}= \left ({\partial^2\over\partial r^2}u \right )_{r=a_0}(a-a_0)</math>
Which is clearly linear elasticity.
Note that the derivation is done considering two neighboring atoms, so the Hook's coefficient is:
- <math>K=a_0{dF \over dr}=a_0 \left ({\partial^2\over\partial r^2}u \right )_{r=a_0}</math>
This form can be easily extended to 3-D case, with volume per atom(Ω) in place of interatomic distance.
- <math>K=\Omega_0 \left ({\partial^2\over\partial \Omega^2}u \right )_{\Omega=\Omega_0}</math>
Relationship with atomic radius
As derived above, the bulk modulus is directly related to the interatomic potential and the volume per atom. We can further evaluate the interatomic potential to connect K with other properties. Usually, the interatomic pair potential can be expressed as a function of distance that has two terms, one term for attraction and another term for repulsion. For example,
- <math>u=-Ar^{-n}+Br^{-m}</math>
where the term involving A represents the attraction term and the B term represents repulsion. A and B are both chosen to be positive and n and m are usually integers, with m usually larger than n due to the short-ranged nature of repulsion. At the equilibrium position, u is at its minimum and so the first derivative is 0. We have
- <math>\left ({\partial u \over \partial r} \right )_{r_0}=Anr^{-n-1}+-Bmr^{-m-1}=0</math>
- <math>{B \over A}={n \over m}r_0^{m-n}</math>
- <math>u=-Ar^{-n} \left (1-{B \over A}r^{n-m} \right )=-Ar^{-n} \left (1-{n \over m}r_0^{m-n}r^{n-m} \right )</math>
when r is close to, recall that the n (usually 1 to 6) is smaller than m (usually 9 to 12), ignore the second term, evaluate the second derivative
- <math>\left ({\partial^2\over\partial r^2}u \right )_{r=a_0}=-An(n+1)r_0^{-n-2}</math>
Recall the relationship between r and Ω
- <math>\Omega ={4\pi \over 3}r^3 </math>
- <math>\left ({\partial^2\over\partial \Omega^2}u \right )= \left ({\partial^2\over\partial r^2}u \right ) \left ({\partial r \over \partial \Omega } \right )^2=\left ({\partial^2\over\partial r^2}u \right )\Omega^{-\frac{4}{3}}</math>
- <math>K=\Omega_0 \left ({\partial^2 u \over \partial r^2} \right )_{\Omega =\Omega_0} \propto r_0^{-n-3}</math>
In many cases, such as in metal or ionic material, the attraction force is electrostatic, so n = 1, we have
- <math>K\propto r_0^{-4}</math>
This applies to atoms with similar bonding nature. This relationship is verified within alkali metals and many ionic compounds.[10]
See also
References
Further reading
Шаблон:Elastic moduli Шаблон:Authority control
- ↑ Шаблон:Cite web
- ↑ Page 52 of "Introduction to Solid State Physics, 8th edition" by Charles Kittel, 2005, Шаблон:ISBN
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite web
- ↑ Шаблон:Cite web
- ↑ Шаблон:Cite web
- ↑ Liu, A. Y.; Cohen, M. L. (1989). "Prediction of New Low Compressibility Solids". Science. 245 (4920): 841–842.
- ↑ Beau, M. R. (2018). "On the nature of space-time, cosmological inflation, and expansion of the universe". Preprint. DOI:10.13140/RG.2.2.16796.95364
- ↑ Шаблон:Cite book
- ↑ Шаблон:Cite book