テイト方程式のソースを表示

'''テイト方程式'''（テイトほうていしき、{{lang-en-short|Tait equation}}）は、[[流体力学]]において、液体[[密度]]と静水圧[[圧力]]を関連付ける[[状態方程式 (熱力学)|状態方程式]]である。この方程式は1888年に[[ピーター・ガスリー・テイト]]によって発表され、以下の形で表される<ref>{{cite book |last1=Tait |first1=P. G. |title=Physics and Chemistry of the Voyage of H.M.S. Challenger |date=1888 |volume=II, part IV |chapter=Report on some of the physical properties of fresh water and of sea water}}</ref>。

:<math>
   \frac{V_0 - V}{P V_0} = \frac{A}{\Pi + P}
 </math>

ここで、<math>P</math>は大気圧に加わる静水圧、<math>V_0</math>は大気圧での体積、<math>V</math>は追加圧力<math>P</math>下での体積を表し、<math>A</math>および<math>\Pi</math>は実験により決定されるパラメータである。<math>A</math>および<math>\Pi</math>の2つのパラメータの物理的解釈を含むテイト方程式の詳細な歴史的研究が参照されている<ref name="ISTE - WILEY">{{cite book |last1=Aitken |first1=Frederic |last2=Foulc |first2=Jean-Numa |title=From Deep Sea to Laboratory 3:From Tait's Work on the Compressibility of Seawater to Equations-of-State for Liquids |date=2019 |url=http://www.iste.co.uk/book.php?id=1534 |publisher=ISTE - WILEY |location=London, UK |isbn=9781786303769}}</ref>。

== テイト・タンマン状態方程式 ==
1895年には<ref>{{cite journal |last1=Tammann |first1=G. |title=Über die Abhängigkeit der volumina von Lösungen vom druck |date=1895 |journal=Zeitschrift für Physikalische Chemie |volume=17 |pages=620–636}}</ref><ref name=Hayward/>、元の等温テイト方程式がタンマンにより以下のように置き換えられた。

:<math>
  {K} = -{V_0} \left ( \frac{\partial P}{\partial V} \right )_T = {V_0}\frac{(B + P)}{C}
 </math>

ここで、<math>K</math>は等温の混合体積弾性率を示す。この上記の方程式は一般にテイト方程式として知られている。また、積分形は次のように表される。

:<math> V = V_0 - C \log \left(\frac{B+P}{B+P_0}\right)</math>

ここで、<math>V</math>は物質の[[比容積]]（単位は[[リットル|ml]]/[[グラム|g]]またはm<sup>3</sup>/kg）、<math>V_0</math>は<math>P = P_0</math>のときの比容積、<math>B</math>（<math>P_0</math>と同じ単位）および<math>C</math>（<math>V_0</math>と同じ単位）は[[温度]]の関数である。

=== 圧力式 ===
比容積から見た圧力の式は次のようになる。

:<math>
   P = (B + P_0)\exp\left(-\frac{V - V_0}{C}\right) - B \,
 </math>

テイト・タンマン状態方程式に関する研究では、経験的パラメータ<math>C</math>と<math>B</math>の物理的解釈が、参考文献の第3章に記述されている<ref name="ISTE - WILEY"/>。まず、水、海水、ヘリウム-4、ヘリウム-3について、液相全体における臨界温度<math>T_c</math>までの<math>C</math>および<math>B</math>の温度依存式についてである。次に、水の過冷却相については、参考文献の付録Dで論じられている<ref>{{cite journal |last1=Aitken |first1=F. |last2=Volino |first2=F. |title=A new single equation of state to describe the dynamic viscosity and self-diffusion coefficient for all fluid phases of water from 200 to 1800 K based on a new original microscopic model |journal=Physics of Fluids |date=November 2021 |volume=33 |issue=11 |pages=117112 |doi=10.1063/5.0069488|arxiv=2108.10666 |bibcode=2021PhFl...33k7112A |s2cid=237278734 }}</ref>。さらに、液体アルゴンについては、[[三重点]]温度から148Kまでの範囲で詳細に扱われており、これは参考文献の第6節に記載されている<ref>{{cite journal |last1=Aitken |first1=Frédéric |last2=Denat |first2=André |last3=Volino |first3=Ferdinand |title=A New Non-Extensive Equation of State for the Fluid Phases of Argon, Including the Metastable States, from the Melting Line to 2300 K and 50 GPa |journal=Fluids |date=24 April 2024 |volume=9 |issue=5 |pages=102 |doi=10.3390/fluids9050102 |doi-access=free |arxiv=1504.00633 }}</ref>。

== テイト・マーナハン状態方程式 ==
[[File:WaterTaitMacDonaldEOSSpVol.pdf|thumb|テイト・マーナハン状態方程式により予測される、圧力の関数としての比容積。]]
もう一つのテイト方程式として広く知られる<ref name=Thompson>Thompson, P. A., & Beavers, G. S. (1972). Compressible-fluid dynamics. Journal of Applied Mechanics, 39, 366.</ref><ref>Kedrinskiy, V. K. (2006). Hydrodynamics of Explosion: experiments and models. Springer Science & Business Media.</ref>、等温状態方程式は、マーナハンモデルである<ref name=MacDonald>Macdonald, J. R. (1966). Some simple isothermal equations of state. Reviews of Modern Physics, 38(4), 669.</ref>。 これは、以下の形で表される。

:<math>
   \frac{V}{V_0} = \left[1 + \frac{n}{K_0}\,(P - P_0)\right]^{-1/n}
 </math>

ここで、<math>V</math>は圧力<math>P</math>における比容積、<math>V_0</math>は圧力<math>P_0</math>における比容積、<math>K_0</math>は圧力<math>P_0</math>における体積弾性率、<math>n</math>は材料の特性を示すパラメータである。

=== 圧力式 ===
この式は、圧力式で次のように置き換えられる。

:<math>
   P = \frac{K_0}{n} \left[\left(\frac{V_0}{V}\right)^n - 1\right] + P_0 = \frac{K_0}{n} \left[\left(\frac{\rho}{\rho_0}\right)^n - 1\right] + P_0
 </math>

ここで、<math>\rho, \rho_0</math>はそれぞれ圧力<math>P, P_0</math>における質量密度を示す。純粋な水に対して、典型的なパラメータは以下の通りである{{Citation needed|reason=Compressibility of water seems rather around 6?|date=November 2016}}。

* <math>P_0</math> = 101,325 Pa
* <math>\rho_0</math> = 1000&nbsp;kg/cu.m
* <math>K_0</math> = 2.15 GPa
* <math>n</math> = 7.15

この形のテイト状態方程式は、[[マーナハンの状態方程式]]の形と同一であることに注意する必要がある。

=== 体積弾性率の式 ===
マクドナルド・テイトモデルによって予測される接線体積弾性率は以下の通りである。

:<math>
   K = K_0\left(\frac{V_0}{V}\right)^n \,
 </math>

== トゥムリルツ・テイト・タンマン状態方程式 ==
[[File:WaterTaitEOS.svg|250px|thumbnail|right|純水の実験データへのフィットに基づくトゥムリルツ・テイト・タンマン状態方程式。]]
液体のモデル化に使用できる関連する状態方程式として、トゥムリルツ方程式（時にタンマン方程式とも呼ばれる。）がある。この方程式は、最初にトゥムリルツが1909年に提案し、タンマンが1911年に純水のために再提案した<ref name=Hayward/><ref name=Fisher>Fisher, F. H., and O. E. Dial Jr. Equation of state of pure water and sea water. No. MPL-U-99/67. SCRIPPS INSTITUTION OF OCEANOGRAPHY LA JOLLA CA MARINE PHYSICAL LAB, 1975. http://www.dtic.mil/dtic/tr/fulltext/u2/a017775.pdf</ref>。この方程式の形は次の通りである。

:<math>
   V(P,S,T) = V_\infty - K_1 S + \frac{\lambda}{P_0 + K_2 S + P}
 </math>

<math>V(P,S,T)</math>は比容積、<math>P</math>は圧力、<math>S</math>は塩分濃度、<math>T</math>は温度、<math>V_\infty</math>は圧力が無限大（<math>P=\infty</math>）のときの比容積を表す。また、<math>K_1</math>、<math>K_2</math>、および<math>P_0</math>は、実験データに基づいて調整可能なパラメーターである。

淡水（<math>S=0</math>）の場合のトゥムリルツ・タンマン形式のテイト方程式は以下の通りである。

:<math>
   V = V_\infty + \frac{\lambda}{P_0 + P} \,
 </math>

純水の場合、<math>V_\infty</math>、<math>\lambda</math>、<math>P_0</math>の値は、以下の関係式で表される<ref name=Fisher/>。

:<math>
   \begin{align}
     \lambda &= 1788.316 + 21.55053\, T - 0.4695911\, T^2 + 3.096363 \times 10^{-3}\, T^3 - 0.7341182 \times 10^{-5}\, T^4 \\
     P_0 &= 5918.499 + 58.05267\, T - 1.1253317\, T^2 + 6.6123869 \times 10^{-3}\, T^3 - 1.4661625 \times 10^{-5}\, T^4 \\
     V_\infty &= 0.6980547 - 0.7435626 \times 10^{-3}\, T + 0.3704258 \times 10^{-4}\, T^2 - 0.6315724 \times 10^{-6}\, T^3 \\
               & + 0.9829576 \times 10^{-8}\, T^4 - 0.1197269 \times 10^{-9} \,T^5 + 0.1005461 \times 10^{-11}\,T^6 \\
               & - 0.5437898 \times 10^{-14} \,T^7 +  0.169946 \times 10^{-16}\, T^8 
               - 0.2295063 \times 10^{-19}\, T^9
   \end{align}
 </math>

温度<math>T</math>は摂氏、<math>P_0</math>はバー（bars）、<math>V_\infty</math>は立方センチメートル毎グラム（cc/gm）、<math>\lambda</math>はバー-立方センチメートル毎グラム（bars-cc/gm）単位である。

=== 圧力式 ===
圧力の比容積に対する逆トゥムリルツ・テイト・タンマン関係は次のようになる。

:<math>
   P =  \frac{\lambda}{V - V_\infty} - P_0 \,
 </math>

=== 体積弾性率の式 ===

純水の即時接線[[体積弾性率]]のトゥムリルツ・テイト・タンマン公式は、圧力<math>P</math>の二次関数として表現される<ref name=Hayward>Hayward, A. T. J. (1967). Compressibility equations for liquids: a comparative study. British Journal of Applied Physics, 18(7), 965. http://mitran-lab.amath.unc.edu:8081/subversion/Lithotripsy/MultiphysicsFocusing/biblio/TaitEquationOfState/Hayward_CompressEqnsLiquidsComparative1967.pdf</ref>。

:<math>
   K = -V\,\frac{\partial P}{\partial V} = \frac{V\,\lambda}{(V - V_\infty)^2} = (P_0 + P) + \frac{V_\infty}{\lambda}(P_0 + P)^2\,
 </math>

== 修正したテイト状態方程式 ==
水中爆発の研究、特に発生する衝撃波に関連して、[[ジョン・G・カークウッド]]は1965年に<ref>{{cite book |last1=Cole |first1=R. H. |title=Underwater Explosions |date=1965 |publisher=Dover Publications |location=New York}}</ref>、1 kbarを超える高圧を記述するために、より適切な状態方程式の形式を提案した。この形式では、等エントロピー圧縮率係数を以下のように表現する。

:<math>
   -\frac{1}{V} \left ( \frac{\partial V}{\partial P} \right )_S = \frac{1}{n(B + P)}
 </math>

ここで、<math>S</math>はエントロピーを表す。さらに、2つの経験的パラメータ<math>n</math>および<math>B</math>はエントロピーの関数であり、それぞれ以下の性質を持つ。

* <math>n</math>：無次元

* <math>B</math>：<math>P</math>（圧力）と同じ単位

積分を行うことで、等エントロピー条件（<math>S</math>）における体積<math>V(P, S)</math>に関する次の式が得られる。

:<math>
   \frac{V}{V_0} = \left(1 + \frac{P_0}{B}\right)^{1/n} \left(1 + \frac{P}{B}\right)^{-1/n}
 </math>

このとき、

:<math>V_0 = V(P_0, S)</math>

である。

=== 圧力式 ===
等エントロピー条件（<math>S</math>）における比容積に基づく圧力<math>P(V, S)</math>の式は以下の通りである。

:<math>
   P = (B + P_0)\,\left(\cfrac{V_0}{V}\right)^n - B \,
 </math>

修正されたテイト状態方程式に関する詳細な研究では、2つの経験的パラメータ<math>n</math>と<math>B</math>の物理的解釈が第4章で説明されている<ref name="ISTE - WILEY"/>。

また、水、ヘリウム3、ヘリウム4に対して、これらの経験的パラメータ<math>n</math>と<math>B</math>のエントロピーに基づく式が示されている。

== 関連項目 ==
* [[状態方程式 (熱力学)]]

== 脚注 ==
{{脚注ヘルプ}}
{{Reflist}}

{{DEFAULTSORT:ていとほうていしき}}
[[Category:状態方程式 (熱力学)]]
[[Category:流体力学の方程式]]
[[Category:数学に関する記事]]
[[Category:人名を冠した数式]]
[[Category:流体力学]]
[[Category:化学工学]]