VPMechanics 1

The VPM Algorithm

Copyright 1999-2001 Eric Maiken

This notebook discusses the mechanics of the VPM decompression algorithm developed by Yount and Hoffman (Y&H). With just a little mirth, I refer to this as the Increased Gradient Bubble Model.

A non-staged-decompression dive to 100 ft for 16 min on air serves to illustrate how the algorithm handles the ascent calculation by sequentially updating the allowable supersaturation gradients.

Program Settings

The dive extends to the 16-min NDL predicted by the VPM, with a linear descent rate of 100 ft/min to 100 feet, followed by a linear ascent at 33 ft/min to the first decompression stop.
Plots of the supersaturation gradients Pss are displayed along with decompression tables for the first and final iteration of the VPM program. The gradients are tracked for 1hr after surfacing.
The constants used in the calculations differ from the Yount and Hoffman values (1986 paper).
I used:
Symmetric ZHL-16 half-time set: {4.,8.,12.5,18.5,27.,38.3,64.3,77.,109.,146.,187.,239.,305.,390.,498.,635.} [min]
[Graphics:Images/VPMech1_gr_1.gif]
pp[Graphics:Images/VPMech1_gr_2.gif] + pp[Graphics:Images/VPMech1_gr_3.gif] = [Graphics:Images/VPMech1_gr_4.gif] [atm]
[Graphics:Images/VPMech1_gr_5.gif] = 0.79
[Graphics:Images/VPMech1_gr_6.gif] = 0.21
pp[Graphics:Images/VPMech1_gr_7.gif] = [Graphics:Images/VPMech1_gr_8.gif] ()

pDeepest = 1 + [Graphics:Images/VPMech1_gr_9.gif] [ata]
pSurface = 1 [ata]

First iteration

On the first iteration, the following supersaturation was used to calculate a conservative ascent schedule
[Graphics:Images/VPMech1_gr_10.gif]
Pss = 0.54 [atm] for all compartments and all stages.

FIRST ITERATION ASCENT SCHEDULE
    (25 min Total Decompression)

          Depth(ft)   Run-time(min)    [Graphics:Images/VPMech1_gr_11.gif]

0.` 0.`
[Graphics:Images/VPMech1_gr_12.gif]
[Graphics:Images/VPMech1_gr_13.gif]
[Graphics:Images/VPMech1_gr_14.gif] [Graphics:Images/VPMech1_gr_15.gif]
[Graphics:Images/VPMech1_gr_16.gif]
[Graphics:Images/VPMech1_gr_17.gif]
[Graphics:Images/VPMech1_gr_18.gif] [Graphics:Images/VPMech1_gr_19.gif]
[Graphics:Images/VPMech1_gr_20.gif]
[Graphics:Images/VPMech1_gr_21.gif]
[Graphics:Images/VPMech1_gr_22.gif] [Graphics:Images/VPMech1_gr_23.gif]
[Graphics:Images/VPMech1_gr_24.gif]
[Graphics:Images/VPMech1_gr_25.gif]
[Graphics:Images/VPMech1_gr_26.gif] [Graphics:Images/VPMech1_gr_27.gif]
[Graphics:Images/VPMech1_gr_28.gif]
[Graphics:Images/VPMech1_gr_29.gif]
[Graphics:Images/VPMech1_gr_30.gif] [Graphics:Images/VPMech1_gr_31.gif]
[Graphics:Images/VPMech1_gr_32.gif]
[Graphics:Images/VPMech1_gr_33.gif]
[Graphics:Images/VPMech1_gr_34.gif] [Graphics:Images/VPMech1_gr_35.gif]
[Graphics:Images/VPMech1_gr_36.gif]
[Graphics:Images/VPMech1_gr_37.gif]
[Graphics:Images/VPMech1_gr_38.gif] [Graphics:Images/VPMech1_gr_39.gif]
[Graphics:Images/VPMech1_gr_40.gif]
[Graphics:Images/VPMech1_gr_41.gif]
0.` [Graphics:Images/VPMech1_gr_42.gif]
[Graphics:Images/VPMech1_gr_43.gif]
[Graphics:Images/VPMech1_gr_44.gif]

The supersaturation gradients for all 16 compartments are shown below, with red corresponding to 4-min half time, and purple corresponding to  the 635 min half time. After surfacing at a run-time of 41 min, the gradients decay to the inherent-unsaturated value of ~ (0.68 -1) Ata, determined by the choice of alvealor water and carbon dioxide partial pressures.

[Graphics:Images/VPMech1_gr_45.gif]

Final Iteration

After iterations, the set of Pss show dispersion by compartment half-time as determined by the Y&H equations  4a-c.  For each compartment, these are:
[Graphics:Images/VPMech1_gr_46.gif]
The resulting PssNew are tabulated below against ZHL16 compartments:

[Graphics:Images/VPMech1_gr_47.gif]

The ascent schedule converges after four iterations, which subsequently reduce the required decompression (tDeco in eqs. 4) as 24,10,3, and finaly 0 min.
CONVERGED ASCENT SCHEDULE
      (0 minutes decompression)

          Depth(ft)   Run-time(min)  [Graphics:Images/VPMech1_gr_48.gif]

0.` 0.`
[Graphics:Images/VPMech1_gr_49.gif]
[Graphics:Images/VPMech1_gr_50.gif]
[Graphics:Images/VPMech1_gr_51.gif] [Graphics:Images/VPMech1_gr_52.gif]
[Graphics:Images/VPMech1_gr_53.gif]
[Graphics:Images/VPMech1_gr_54.gif]
[Graphics:Images/VPMech1_gr_55.gif] [Graphics:Images/VPMech1_gr_56.gif]
[Graphics:Images/VPMech1_gr_57.gif]
[Graphics:Images/VPMech1_gr_58.gif]
[Graphics:Images/VPMech1_gr_59.gif] [Graphics:Images/VPMech1_gr_60.gif]
[Graphics:Images/VPMech1_gr_61.gif]
[Graphics:Images/VPMech1_gr_62.gif]
[Graphics:Images/VPMech1_gr_63.gif] [Graphics:Images/VPMech1_gr_64.gif]
[Graphics:Images/VPMech1_gr_65.gif]
[Graphics:Images/VPMech1_gr_66.gif]
0.` [Graphics:Images/VPMech1_gr_67.gif]
[Graphics:Images/VPMech1_gr_68.gif]
[Graphics:Images/VPMech1_gr_69.gif]

Plot of the Pss gradients corresponding to the converged ascent schedule

[Graphics:Images/VPMech1_gr_70.gif]


Converted by Mathematica      May 27, 2001