Soustava rovnic popisující přenos tepla v Zemi a její různé aproximace

Z ωικι.matfyz.cz
Přejít na: navigace, hledání
  • ZZHmoty:
$ \frac{D}{Dt}\int_V \rho dV = 0 $
$ \frac{\partial \rho}{\partial t}+\Delta\cdot(\rho v) = 0 $

a

$ [\rho(v-\nu)]_-^+\cdot n =0 $
  • ZZHybnosti
$ \frac{D}{Dt}\int_V\rho v dV = \int_V \rho g dV - 2\int_V \rho \Omega \times v dV - \int_V \rho \Omega \times(\Omega \times r)dV+ \int_S \tau \cdot n dS $
$ \nabla \cdot \tau^T+\rho g-2\rho\Omega\times v-\rho\Omega\times(\Omega\times r)=\rho\frac{Dv}{Dt}=\rho(\frac{\partial v}{\partial t}+v\cdot \nabla v) $

a $ [\rho v(v-\nu)-\tau]_-^+\cdot n = 0 $

  • ZZMomentu hybnosti
$ \frac{D}{Dt}\int_V r\times\rho v dV = \int_V r\times \rho g dV - 2\int_V r\times(\rho \Omega \times v) dV - \int_V r\times( \rho \Omega \times(\Omega \times r))dV+ \int_S r\times \tau \cdot n dS $
$ I\times \tau = 0 $ <=> $ \tau = \tau^T $
  • ZZEnergie
$ \frac{D}{Dt}\int_V (\rho\epsilon+1/2 \rho v \cdot v)dV = \int_S v\cdot \tau \cdot n dS+\int_V \rho g \cdot v dV-2\int_V \rho(\Omega\times v)\cdot v dV - \int_V \rho(\Omega\times(\Omega\times r))\cdot v dV - \int_S \q\cdot n dS + \int_V HdV $

Gibbsova relace

$ \rho T\frac{Ds}{Dt}=\rho\frac{D\epsilon}{Dt}+p\nabla\cdot v $

a z toho

$ \rho T\frac{Ds}{Dt}=\nabla\cdot(k\cdot \nabla T)+\sigma : \nabla v+H $

a

$ [k\cdot\nabla T]_-^+\cdot n = -[v\cdot \tau]_-^+\cdot n+[(\rho\epsilon+1/2\rho v\cdot v)(v-\nu)]_-^+\cdot n $

Dále tepelná rovnice ve stavových proměnných

$ \rho T \frac{Ds}{Dt} = \rho T\frac{\partial s}{\partial T}_V \frac{DT}{Dt}+\rho T\frac{\partial s}{\partial V}_T \frac{DV}{Dt}=\rho c_V \frac{DT}{Dt}+\mu\nabla\cdot v $

kde $ \mu = T(\frac{\partial s}{\partial V})_T = T (\dfrac{\partial p}{\partial T})_V = T [\frac{1}{V}(\frac{\partial V}{\partial T})_p][-V(\frac{\partial p}{\partial V})_T]=\alpha K_T $ kde termální expanze $ \alpha = -\frac{1}{p}(\frac{\partial p}{\partial T})_p $ a izotermální objemový modul $ K_T = -V(\frac{\partial p}{\partial V})_T $

Pak získáme $ \mu=T\alpha K_T = \rho c_V T \gamma $ pro Grüneisenův parametr $ \gamma = \alpha K_T / \rho c_V $

a z toho plyne termální rovnice

$ \rho c_V \frac{\partial T}{\partial t} = \nabla \cdot (k\cdot \nabla T)-\rho c_V v\cdot \nabla T - \rho c_V T\gamma \nabla\cdot v+\sigma: \nabla v+H $
  • Termální rovnice v proměnných p,T
$ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k\cdot \nabla T)-\rho c_p v\cdot \nabla T - \alpha T(\frac{\partial p}{\partial t}+v \cdot \nabla p)+\sigma: \nabla v+H $

kde

$ c_p = c_V(1+\gamma \alpha T) $

a rovnice kontinuity

$ \nabla \cdot v = -\frac{1}{\rho}\frac{D\rho}{Dt}=-K_T^{-1}\frac{Dp}{Dt}+\alpha \frac{DT}{Dt} $
  • Termální rovnice v kontinuu s dominantním hydrostatickým tlakem

platí

$ \nabla p_0 = \rho_0 g_0 - \rho_0\Omega \times(\Omega\times r) $

nechť tedy

$ \frac{\partial p}{\partial t}+v cdot \nabla p = -v_r\rho g $

v radiální složce, kde g je

$ g=|\vec g| = \vec g_0-\Omega\times(\Omega\times r) $

pak dostaneme obvyklou termální rovnici

$ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k\cdot \nabla T)-\rho c_p v\cdot \nabla T - \rho v_r\alpha Tg+\sigma: \nabla v+H $

kde figuruje adiabatický gradient / ohřívání

$ (\frac{\partial T}{\partial r})_s = -\frac{\alpha Tg}{c_p} $
Citováno z „http://wiki.matfyz.cz/index.php?title=Soustava_rovnic_popisující_přenos_tepla_v_Zemi_a_její_různé_aproximace&oldid=21246“