diff options
| author |  Jan Darowski <jan.darowski@gmail.com> | 2015-07-03 22:45:29 +0200 |
|---|---|---|
| committer | Jan Darowski <jan.darowski@gmail.com> | 2015-07-03 22:45:29 +0200 |
| commit | 2435d79c0e6e4e30ace3b23765005c90181ddae5 (patch) | |
| tree | 54ca29170d5cf372975729e85ca164a8b593be9a /deco.c | |
| parent | 0180d2eb1e8101eb362efb00edaae9eb9bcb24aa (diff) | |
| download | subsurface-2435d79c0e6e4e30ace3b23765005c90181ddae5.tar.gz | |
VPM-B: add nuclei size calculation and nuclear regeneration.
This function calculates the size of nuclei at the end of deco,
then simulates their regeneration, to the moment before the deco.
This is redundant as nuclear regeneration is a very slow process.
Function should be called with time in seconds, just before the ascent.
Signed-off-by: Jan Darowski <jan.darowski@gmail.com>
Diffstat (limited to 'deco.c')
| -rw-r--r-- | deco.c | 15 |
1 files changed, 15 insertions, 0 deletions
@@ -206,6 +206,21 @@ double he_factor(int period_in_seconds, int ci) return cache[ci].last_factor; } +void nuclear_regeneration(double time) +{ + time /= 60.0; + int ci; + double crushing_radius_N2, crushing_radius_He; + for (ci = 0; ci < 16; ++ci) { + //rm + crushing_radius_N2 = 1.0 / (max_n2_crushing_pressure[ci] / (2.0 * (vpmb_config.skin_compression_gammaC - vpmb_config.surface_tension_gamma)) + 1.0 / vpmb_config.crit_radius_N2); + crushing_radius_He = 1.0 / (max_he_crushing_pressure[ci] / (2.0 * (vpmb_config.skin_compression_gammaC - vpmb_config.surface_tension_gamma)) + 1.0 / vpmb_config.crit_radius_He); + //rs + n2_regen_radius[ci] = crushing_radius_N2 + (vpmb_config.crit_radius_N2 - crushing_radius_N2) * (1.0 - exp (-time / vpmb_config.regeneration_time)); + he_regen_radius[ci] = crushing_radius_He + (vpmb_config.crit_radius_He - crushing_radius_He) * (1.0 - exp (-time / vpmb_config.regeneration_time)); + } +} + // Calculates the nucleons inner pressure during the impermeable period double calc_inner_pressure(double crit_radius, double onset_tension, double current_ambient_pressure) { |