Darcy-Forchheimerův zákon
Darcy-Forchheimer koriguje zákon Darcyho pro průtoky tekutiny v porézních médiích s přihlédnutím setrvačný účinek. Vychází z tvorby Julesa Dupuita (1863) a Philippa Forchheimera (1901) a ve své současné podobě jej dal JC Ward (1964).
Darcy-Forchheimerův zákon
Když je průtok vysoký, lze vzít v úvahu setrvačné účinky pomocí korekce zahrnující Reynoldsovo číslo na základě charakteristické délkyK.{\ displaystyle {\ sqrt {K}}}
REk=ρ||PROTI||K.μ{\ displaystyle Re_ {k} = {\ frac {\ rho \, || \ mathbf {V} || {\ sqrt {K}}} {\ mu}}}s
-
ρ{\ displaystyle \ rho} hustota kapaliny,
-
PROTI{\ displaystyle \ mathbf {V}}průměrná rychlost v porézním médiu, definovaná kde je hmotnostní tok,q=ρPROTI{\ displaystyle \ mathbf {q} = \ rho \ mathbf {V}}q{\ displaystyle \ mathbf {q}}
-
K.{\ displaystyle K} propustnost , předpokládaná jako skalární,
-
μ{\ displaystyle \ mu} dynamická viskozita .
Tato oprava je vyjádřena v Darcyho-Forchheimerově zákoně následovně
K.(∇p-ρG)=-μPROTI(1+αREk){\ displaystyle K (\ nabla p- \ rho \ mathbf {g} \,) = - \ mu \ mathbf {V} \, (1+ \ alpha Re_ {k})}nebo
-
∇p{\ displaystyle \ nabla p} je tlakový gradient,
-
G{\ displaystyle \ mathbf {g}} je pole zrychlení,
-
α{\ displaystyle \ alpha}je číslo Ergün (nebo číslo Ward). Je to řádově 0,5.
Měření charakteristik prostředí
Lze získat pro plyn současně a ze série laboratorních experimentů, kde je zanedbána jednorozměrná situace .
K.{\ displaystyle K}α{\ displaystyle \ alpha}G{\ displaystyle \ mathbf {g}}
Darcy-Forchheimerova rovnice je napsána:
-dpdX=μK.qρ+αK.q2ρ{\ displaystyle - {\ frac {\ mathrm {d} p} {\ mathrm {d} x}} = {\ frac {\ mu} {K}} {\ frac {q} {\ rho}} + {\ frac {\ alpha} {\ sqrt {K}}} {\ frac {q ^ {2}} {\ rho}}}Stavovou rovnici údajně dokonalého plynu používáme ve tvaru:
p=ρRTM{\ displaystyle p = {\ frac {\ rho RT} {M}}}nebo
Vynásobením výše uvedeného vztahu vlevo a jeho hodnoty vpravo to přijde:
p{\ displaystyle p}
-pdpdX=RTMq(μK.+αK.q){\ displaystyle -p {\ frac {\ mathrm {d} p} {\ mathrm {d} x}} = {\ frac {RT} {M}} q \ vlevo ({\ frac {\ mu} {K} } + {\ frac {\ alpha} {\ sqrt {K}}} {q} \ vpravo)}což dává odhad
-M2RTqμΔ(p2)ΔX=1K.+αK.qμ{\ displaystyle - {\ frac {M} {2RTq \ mu}} {\ frac {\ Delta (p ^ {2})} {\ Delta x}} = {\ frac {1} {K}} + {\ frac {\ alpha} {\ sqrt {K}}} {\ frac {q} {\ mu}}}Změnou tlaku a měřením průtoku v experimentu získáme řadu bodů, ze kterých extrahujeme průsečík y a sklon (viz křivka).
1K.{\ displaystyle {\ frac {1} {K}}}αK.{\ displaystyle {\ frac {\ alpha} {\ sqrt {K}}}}
Reference
-
(in) Ward, JC, Turbulent Flow in Porous Media , Journal Hydraulics Division Proceedings of ASCE , sv. 5, s. 1-12 (1964)
-
(in) Donald A. Nield, Adrian Bejan , Konvekce v porézních médiích , Springer , 2006 ( ISBN 0-387-29096-6 )
-
(in) Cornell, D. a Katz, LD, Flow of Gases Through Porous Media Consolidated , Industrial & Engineering Chemistry Research , sv. 45, 1953, str. 2145–2152
<img src="https://fr.wikipedia.org/wiki/Special:CentralAutoLogin/start?type=1x1" alt="" title="" width="1" height="1" style="border: none; position: absolute;">