mifr/subsurface
1
0
Fork 0
mirror of https://github.com/subsurface/subsurface.git synced 2025-02-19 22:16:15 +00:00
Code Issues Projects Releases Packages Wiki Activity

VPM-B: Use correct gamma values

Browse source 
In other implementations of VPM-B, surface_tension_gamma and
skin_compression_gammaC are taken as 0.0179 N/msw and 0.257 N/msw respectively.
We do pressure calculations in bar, not msw, so our gamma values need to reflect
that.  Previously we had used 0.179 and 2.57, which are close but not quite
correct (10 msw == 1.01325 bar).

Signed-off-by: Rick Walsh <rickmwalsh@gmail.com>
Signed-off-by: Dirk Hohndel <dirk@hohndel.org>
This commit is contained in:
Rick Walsh 2015-08-29 21:34:19 +10:00 • committed by Dirk Hohndel
parent 1e7dd49234
commit becc1e6dc4
1 changed files with 4 additions and 4 deletions

8
deco.c
View file

@ -52,8 +52,8 @@ struct vpmb_config {
double crit_radius_He; //! Critical radius of He nucleon (microns).
double crit_volume_lambda; //! Constant corresponding to critical gas volume (bar * min).
double gradient_of_imperm; //! Gradient after which bubbles become impermeable (bar).
double surface_tension_gamma; //! Nucleons surface tension constant (N / 10m).
double skin_compression_gammaC; //! Skin compression gammaC (N / 10m).
double surface_tension_gamma; //! Nucleons surface tension constant (N / bar = m2).
double skin_compression_gammaC; //! Skin compression gammaC (N / bar = m2).
double regeneration_time; //! Time needed for the bubble to regenerate to the start radius (min).
double other_gases_pressure; //! Always present pressure of other gasses in tissues (bar).
};
@ -63,8 +63,8 @@ struct vpmb_config vpmb_config = {
.crit_radius_He = 0.45,
.crit_volume_lambda = 199.58,
.gradient_of_imperm = 8.2,
.surface_tension_gamma = 0.179,
.skin_compression_gammaC = 2.57,
.surface_tension_gamma = 0.18137175, // = 0.0179 N/msw
.skin_compression_gammaC = 2.6040525, // = 0.257 N/msw
.regeneration_time = 20160.0,
.other_gases_pressure = 0.1359888
};