Governing equations on differential form: Difference between revisions

From Flowpedia
← Older edit
Nian (talk | contribs)
Bureaucrats, Interface administrators, Administrators
580 edits
VisualWikitext
No edit summary
No edit summary
 
(7 intermediate revisions by the same user not shown)
Line 12: Line 12:
--><noinclude><!--
--><noinclude><!--
-->{{#vardefine:secno|2}}<!--
-->{{#vardefine:secno|2}}<!--
-->{{#vardefine:eqno|31}}<!--
-->{{#vardefine:eqno|28}}<!--
--></noinclude><!--
--></noinclude><!--


Line 22: Line 22:
The continuity equation on integral form reads
The continuity equation on integral form reads


{{NumEqn|<math>
{{InfoBox|<math>
\frac{d}{dt}\iiint_{\Omega} \rho dV+\iint_{\partial \Omega}\rho \mathbf{v}\cdot \mathbf{n} dS=0
\frac{d}{dt}\iiint_{\Omega} \rho dV+\iint_{\partial \Omega}\rho \mathbf{v}\cdot \mathbf{n} dS=0
</math>|nonumber=1}}
</math>}}


Apply Gauss's divergence theorem on the surface integral gives
Apply Gauss's divergence theorem on the surface integral gives
Line 30: Line 30:
{{NumEqn|<math>
{{NumEqn|<math>
\iint_{\partial \Omega}\rho \mathbf{v}\cdot \mathbf{n} dS=\iiint_{\Omega}\nabla\cdot(\rho\mathbf{v})dV
\iint_{\partial \Omega}\rho \mathbf{v}\cdot \mathbf{n} dS=\iiint_{\Omega}\nabla\cdot(\rho\mathbf{v})dV
</math>|number=5.34|spar=1|prefix=ekv.|border=1|padding=2em}}
</math>}}
 


Also, if <math>\Omega</math> is a fixed control volume
Also, if <math>\Omega</math> is a fixed control volume
Line 49: Line 48:
{{NumEqn|<math>
{{NumEqn|<math>
\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho\mathbf{v})=0
\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho\mathbf{v})=0
</math>}}
</math>|label=eq-cont-pde}}


which is the continuity equation on partial differential form.
which is the continuity equation on partial differential form.
Line 57: Line 56:
The momentum equation on integral form reads
The momentum equation on integral form reads


{{NumEqn|<math>
{{InfoBox|<math>
\frac{d}{dt}\iiint_{\Omega} \rho \mathbf{v} dV+\iint_{\partial \Omega} \left[(\rho\mathbf{v}\cdot\mathbf{n})\mathbf{v}+p\mathbf{n}\right]dS=\iiint_{\Omega}\rho \mathbf{f}dV
\frac{d}{dt}\iiint_{\Omega} \rho \mathbf{v} dV+\iint_{\partial \Omega} \left[(\rho\mathbf{v}\cdot\mathbf{n})\mathbf{v}+p\mathbf{n}\right]dS=\iiint_{\Omega}\rho \mathbf{f}dV
</math>|nonumber=1}}
</math>}}


As for the continuity equation, the surface integral terms are rewritten as volume integrals using Gauss's divergence theorem.
As for the continuity equation, the surface integral terms are rewritten as volume integrals using Gauss's divergence theorem.
Line 87: Line 86:
{{NumEqn|<math>
{{NumEqn|<math>
\frac{\partial}{\partial t}(\rho \mathbf{v}) + \nabla\cdot(\rho \mathbf{v}\mathbf{v}) + \nabla p = \rho \mathbf{f}  
\frac{\partial}{\partial t}(\rho \mathbf{v}) + \nabla\cdot(\rho \mathbf{v}\mathbf{v}) + \nabla p = \rho \mathbf{f}  
</math>}}
</math>|label=eq-mom-pde}}


which is the momentum equation on partial differential form
which is the momentum equation on partial differential form
Line 95: Line 94:
The energy equation on integral form reads
The energy equation on integral form reads


{{NumEqn|<math>
{{InfoBox|<math>
\frac{d}{dt}\iiint_{\Omega}\rho e_o dV+\iint_{\partial \Omega}\rho h_o(\mathbf{v}\cdot\mathbf{n})dS=</math><br><br><math>\iiint_{\Omega}\rho\mathbf{f}\cdot\mathbf{v}dV+\iiint_{\Omega} \dot{q}\rho dV
\frac{d}{dt}\iiint_{\Omega}\rho e_o dV+\iint_{\partial \Omega}\rho h_o(\mathbf{v}\cdot\mathbf{n})dS=</math><br><br><math>\iiint_{\Omega}\rho\mathbf{f}\cdot\mathbf{v}dV+\iiint_{\Omega} \dot{q}\rho dV
</math>|nonumber=1}}
</math>}}


Gauss's divergence theorem applied to the surface integral term in the energy equation gives
Gauss's divergence theorem applied to the surface integral term in the energy equation gives
Line 121: Line 120:
{{NumEqn|<math>
{{NumEqn|<math>
\frac{\partial}{\partial t}(\rho e_o) + \nabla\cdot(\rho h_o\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
\frac{\partial}{\partial t}(\rho e_o) + \nabla\cdot(\rho h_o\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
</math>}}
</math>|label=eq-energy-pde}}


which is the energy equation on partial differential form
which is the energy equation on partial differential form
Line 129: Line 128:
The governing equations for compressible inviscid flow on partial differential form:
The governing equations for compressible inviscid flow on partial differential form:


{{NumEqn|<math>
<div style="border: solid 1px;">
{{OpenInfoBox|<math>
\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho\mathbf{v})=0
\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho\mathbf{v})=0
</math>}}
</math>|description=Continuity:}}


{{NumEqn|<math>
{{OpenInfoBox|<math>
\frac{\partial}{\partial t}(\rho \mathbf{v}) + \nabla\cdot(\rho \mathbf{v}\mathbf{v}) + \nabla p = \rho \mathbf{f}
\frac{\partial}{\partial t}(\rho \mathbf{v}) + \nabla\cdot(\rho \mathbf{v}\mathbf{v}) + \nabla p = \rho \mathbf{f}
</math>}}
</math>|description=Momentum:}}


{{NumEqn|<math>
{{OpenInfoBox|<math>
\frac{\partial}{\partial t}(\rho e_o) + \nabla\cdot(\rho h_o\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
\frac{\partial}{\partial t}(\rho e_o) + \nabla\cdot(\rho h_o\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
</math>}}
</math>|description=Energy:}}
</div>


=== The Differential Equations on Non-Conservation Form ===
=== The Differential Equations on Non-Conservation Form ===
Line 149: Line 150:
{{NumEqn|<math>
{{NumEqn|<math>
\frac{D}{Dt}=\frac{\partial}{\partial t}+\mathbf{v}\cdot\nabla
\frac{D}{Dt}=\frac{\partial}{\partial t}+\mathbf{v}\cdot\nabla
</math>}}
</math>|label=eq-cont-pde-non-cons}}


where the first term of the right hand side is the local derivative and the second term is the convective derivative.
where the first term of the right hand side is the local derivative and the second term is the convective derivative.
Line 205: Line 206:
{{NumEqn|<math>
{{NumEqn|<math>
\frac{D\mathbf{v}}{Dt}+\frac{1}{\rho}\nabla p = \mathbf{f}
\frac{D\mathbf{v}}{Dt}+\frac{1}{\rho}\nabla p = \mathbf{f}
</math>}}
</math>|label=eq-mom-pde-non-cons}}


==== Conservation of Energy ====
==== Conservation of Energy ====


The last equation on non-conservation differential form is the energy equation. We start by rewriting the energy equation on differential form (Eqn. \ref{eq:governing:energy:pde}), repeated here for convenience
The last equation on non-conservation differential form is the energy equation. We start by rewriting the energy equation on differential form {{EquationNote|label=eq-energy-pde}}, repeated here for convenience


{{NumEqn|<math>
{{NumEqn|<math>
\frac{\partial}{\partial t}(\rho e_o) + \nabla\cdot(\rho h_o\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
\frac{\partial}{\partial t}(\rho e_o) + \nabla\cdot(\rho h_o\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
</math>}}
</math>|nonumber=1}}


Total enthalpy, <math>h_o</math>, is replaced with total energy, <math>e_o</math>
Total enthalpy, <math>h_o</math>, is replaced with total energy, <math>e_o</math>
Line 230: Line 231:


{{NumEqn|<math>
{{NumEqn|<math>
\rho\frac{\partial e_o}{\partial t} + e_o\frac{\partial \rho}{\partial t} + \rho\mathbf{v}\cdot\nabla e_o + e_o\nabla\cdot(\rho \mathbf{v}) + \nabla\cdot(p\mathbf{v})= \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
\rho\frac{\partial e_o}{\partial t} + e_o\frac{\partial \rho}{\partial t} + \rho\mathbf{v}\cdot\nabla e_o + e_o\nabla\cdot(\rho \mathbf{v}) + \nabla\cdot(p\mathbf{v})=</math><br><br><math>= \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
</math>}}
</math>}}


Line 243: Line 244:
{{NumEqn|<math>
{{NumEqn|<math>
\rho\frac{De_o}{Dt} + \nabla\cdot(p\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
\rho\frac{De_o}{Dt} + \nabla\cdot(p\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
</math>}}
</math>|label=eq-energy-pde-non-cons}}


==== Summary ====
==== Summary ====


Continuity:
<div style="border: solid 1px;">
 
{{OpenInfoBox|<math>
{{NumEqn|<math>
\frac{D\rho}{Dt}+\rho(\nabla\cdot\mathbf{v})=0
\frac{D\rho}{Dt}+\rho(\nabla\cdot\mathbf{v})=0
</math>}}
</math>|description=Continuity:}}


Momentum:
{{OpenInfoBox|<math>
 
{{NumEqn|<math>
\frac{D\mathbf{v}}{Dt}+\frac{1}{\rho}\nabla p = \mathbf{f}
\frac{D\mathbf{v}}{Dt}+\frac{1}{\rho}\nabla p = \mathbf{f}
</math>}}
</math>|description=Momentum:}}


Energy:
{{OpenInfoBox|<math>
 
{{NumEqn|<math>
\rho\frac{De_o}{Dt} + \nabla\cdot(p\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
\rho\frac{De_o}{Dt} + \nabla\cdot(p\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho
</math>}}
</math>|description=Energy:}}
</div>


=== Alternative Forms of the Energy Equation ===
=== Alternative Forms of the Energy Equation ===
Line 275: Line 272:
</math>}}
</math>}}


Inserted in Eqn. \ref{eq:governing:energy:non}, this gives
Inserted in {{EquationNote|label=eq-energy-pde-non-cons|nopar=1}}, this gives


{{NumEqn|<math>
{{NumEqn|<math>
Line 281: Line 278:
</math>}}
</math>}}


Now, let's replace the substantial derivative <math>D\mathbf{v}/Dt</math> using the momentum equation on non-conservation form (Eqn. \ref{eq:governing:mom:non}).
Now, let's replace the substantial derivative <math>D\mathbf{v}/Dt</math> using the momentum equation on non-conservation form {{EquationNote|label=eq-mom-pde-non-cons}}.


{{NumEqn|<math>
{{NumEqn|<math>
Line 290: Line 287:


{{NumEqn|<math>
{{NumEqn|<math>
\rho\frac{De}{Dt} \cancel{-\mathbf{v}\cdot\nabla p} + \cancel{\mathbf{v}\cdot\nabla p} +  p(\nabla\cdot\mathbf{v}) = \dot{q}\rho\Rightarrow \rho\frac{De}{Dt} + p(\nabla\cdot\mathbf{v}) = \dot{q}\rho
\rho\frac{De}{Dt} \cancel{-\mathbf{v}\cdot\nabla p} + \cancel{\mathbf{v}\cdot\nabla p} +  p(\nabla\cdot\mathbf{v}) = \dot{q}\rho\Rightarrow</math><br><br><math>\Rightarrow\rho\frac{De}{Dt} + p(\nabla\cdot\mathbf{v}) = \dot{q}\rho
</math>}}
</math>}}


Line 297: Line 294:
{{NumEqn|<math>
{{NumEqn|<math>
\frac{De}{Dt} + \frac{p}{\rho}(\nabla\cdot\mathbf{v}) = \dot{q}
\frac{De}{Dt} + \frac{p}{\rho}(\nabla\cdot\mathbf{v}) = \dot{q}
</math>}}
</math>|label=eq-energy-pde-non-cons-b}}


Conservation of mass gives
Conservation of mass gives
Line 305: Line 302:
</math>}}
</math>}}


Insert in Eqn. \ref{eq:governing:energy:non:b}
Insert in {{EquationNote|label=eq-energy-pde-non-cons-b|nopar=1}}


{{NumEqn|<math>
{{NumEqn|<math>
Line 323: Line 320:
</math>}}
</math>}}


with <math>De/Dt</math> from Eqn. \ref{eq:governing:energy:non:b}
with <math>De/Dt</math> from {{EquationNote|label=eq-energy-pde-non-cons-b|nopar=1}}


{{NumEqn|<math>
{{NumEqn|<math>
Line 331: Line 328:
{{NumEqn|<math>
{{NumEqn|<math>
\frac{Dh}{Dt}=\dot{q} + \frac{1}{\rho}\frac{Dp}{Dt}
\frac{Dh}{Dt}=\dot{q} + \frac{1}{\rho}\frac{Dp}{Dt}
</math>}}
</math>|label=eq-energy-pde-non-cons-c}}


==== Total Enthalpy Formulation ====
==== Total Enthalpy Formulation ====
Line 339: Line 336:
</math>}}
</math>}}


From the momentum equation (Eqn. \ref{eq:governing:mom:non})
From the momentum equation {{EquationNote|label=eq-mom-pde-non-cons}}


{{NumEqn|<math>
{{NumEqn|<math>
Line 351: Line 348:
</math>}}
</math>}}


Inserting <math>Dh/Dt</math> from Eqn. \ref{eq:governing:energy:non:c} gives
Inserting <math>Dh/Dt</math> from {{EquationNote|label=eq-energy-pde-non-cons-c|nopar=1}} gives


{{NumEqn|<math>
{{NumEqn|<math>
\frac{Dh_o}{Dt}=\dot{q} + \frac{1}{\rho}\frac{Dp}{Dt}+\mathbf{v}\cdot\mathbf{f} -\frac{1}{\rho}\mathbf{v}\cdot\nabla p = \frac{1}{\rho}\left[\frac{Dp}{Dt}-\mathbf{v}\cdot\nabla p\right] + \dot{q} + \mathbf{v}\cdot\mathbf{f}
\frac{Dh_o}{Dt}=\dot{q} + \frac{1}{\rho}\frac{Dp}{Dt}+\mathbf{v}\cdot\mathbf{f} -\frac{1}{\rho}\mathbf{v}\cdot\nabla p =</math><br><br><math>=\frac{1}{\rho}\left[\frac{Dp}{Dt}-\mathbf{v}\cdot\nabla p\right] + \dot{q} + \mathbf{v}\cdot\mathbf{f}
</math>}}
</math>}}


Retrieved from "https://wiki.g3dflow.com/index.php?title=Governing_equations_on_differential_form"