W dniu 2022-01-29 o 20:36, WM pisze:
>
> I tu utknąłem, bo mam równanie ciśnienia względnego, a potrzebuję ciśnienia bezwzględnego.
>
> Muszę założyć ciśnienie w środku Ziemi i wtedy mam określony jego przebieg po promieniu dla
> izotermy :)))
Może należy założyć, że na powierzchni powinno być 101325 Pa?
Ponieważ brakuje konkretnych odpowiedzi na tytułowe
pytanie, proponuję taki symulator i jego ulepszanie
(można uruchomić np. w kompilatorze online
https://www.onlinegdb.com/online_c++_compiler
- wkleić kod w okienku z kodem i nacisnąć "Run"):
#include <iostream>
int main()
{
double depth = 0; // głębokość pod powierzchnią Ziemi
double p0 = 101325; // ciśnienie na powierzchni Ziemi
double R = 6371000; // promień Ziemi
double g0 = 9.81; // przyspieszenie na pow. Ziemi
double dh = 10; // przyrost (krok) głębokości
double ro0 = 1.293; // gęstość powietrza na pow. Ziemi
double pkr = 7.4e6; // ciśnienie krytyczne
double ro; // gęstość powietrza na głębokości dh
double m; // masa słupa powietrza o wysokości dh
// na głębokości depth
double g; // przyspieszenie na głębokości depth
double p; // ciśnienie na głębokości depth
p = p0; // zaczynamy od powierzchni z ciśnieniem p0
// schodzimy głębiej o dh metrów
while (depth <= 3500) // max głębokość w metrach, <= R
{
g = g0 * (R - depth) / R; // liniowo do 0 w środku Ziemi
if (p < pkr)
ro = (p / p0) * ro0; // gęstość powietrza na głęb. depth
else
ro = (pkr / p0) * ro0; // powyżej pkr nie zwiększamy gęstości
m = dh * ro; // masa słupa powietrza o wysokości dh i przekroju 1m2
p = p + g * m; // ciśnienie większe o nacisk słupa na m2
depth = depth + dh; // jesteśmy o krok niżej
std::cout << "depth = " << depth << " m: "
<< "p = " << p / p0 << " atm\n";
}
return 0;
}
Zakładam, że po osiągnięciu ciśnienia 7.4 MPa powietrze już się
nie spręża, temperatura jest stała itp. (bardzo uproszczone).
Wyniki:
depth = 0 m: p = 1.00000 atm
depth = 10 m: p = 1.00125 atm
depth = 20 m: p = 1.00251 atm
depth = 30 m: p = 1.00376 atm
...
depth = 3500 m: p = 1.54922 atm
--
depth = 10000 m: p = 3.49063 atm
...
depth = 6371000 m: p = 28883 atm
Może ktoś sprawdzi / poprawi tę symulację?
P.