VPM-B Gradients as Gradient factors

In a comment, Rick suggested that given that in the end, VPM-B computes depth dependent maximal gradients, one could compare those to the depth dependent gradients computed from the Buehlmann model and expressed as gradient factors (i.e. as percentage of the gradient of a plain vanilla Buehlmann gradient). Here is what I found when I did this.

 

For concreteness, I computed deco schedules with Subsurface for a dive to 120m leaving the bottom at runtime 20min with TMX18/50 (yes, I know this is far too much pO2, but for deco that does not matter), EAN50 and EAN100. For some reason, my ascent speeds are set to 9m/min up to 75% average depth, then 6m/min up to 6m and 1m/min for the final ascent. With Bühlmann with GF 30/80, I obtain a plan that reaches the surface at a runtime of 188min while for VPM-B+1 it takes only 141min.

Actually, I patched subsurface a little bit (see this branch on github) that, during deco, it prints out the allowed gradient for all tissues at all deco stops. This output I ran thru a mathematica notebook  to produce these plots:

 

On the x-axis, we have depth in meters while on the y-axis, we have the current gradient factor. The different colored dots are the different tissues (red being those with the shortest halft-time) with the thick dot indicating the leading tissue, i.e. the tissue that is causing the current stop.

The first diagram is Bühlmann while the second is VPM-B.

A few remarks are in order: First, why is Bühlmann not a single line from GF 30 to GF 40? This is because of the way I compute things: The exact Bühlmann a and b coefficients depend on the gas mixture currently in the tissue. In addition, what is actually computed is the maximal ambient pressure at the ceiling (which is the current depth only for the leading tissue) for that tissue and the gradient I plot is the difference between that and the current ambient pressure. So this is not exactly the right thing for tissues that are not leading. But still we see, it is a nice band of values and they all nicely move between 30% and 80%.

For the VPM-B plot, on the other hand, we see that the variation within a tissue is much less but the gradient factor varies a lot more between tissues. After all, the actual variation in gradient factor with depth comes about mainly because the leading tissue is changing and much less from the cubic curve that implements the Boyle compensation.

The other thing that is striking is the scale: At the last deco stop, the effective gradient factor is 124% for the leading tissue (while for some slower ones it is as high 180% (which never matters as that tissue never leads). In fact, a Bühlmann plan with GF 20/125 leads to a very similar profile.

Here are the actual profiles: First, Bühlmann GF 30/80

Depth

Duration

Runtime

Gas

120m 7min 7min (18/50)
120m 13min 20min
51m 9min 29min
51m 2min 31min
48m 1min 32min
48m 1min 32min
45m 1min 33min
45m 2min 34min
42m 1min 35min
42m 2min 36min
39m 1min 37min
39m 4min 40min
36m 1min 41min
36m 3min 43min
33m 1min 44min
33m 4min 47min
30m 1min 48min
30m 6min 53min
27m 1min 54min
27m 6min 59min
24m 1min 60min
24m 10min 69min
21m 1min 70min
21m 5min 74min EAN50
18m 1min 75min
18m 7min 81min
15m 1min 82min
15m 8min 89min
12m 1min 90min
12m 13min 102min
9m 1min 103min
9m 20min 122min
6m 1min 123min
6m 20min 142min EAN100
3m 3min 145min
3m 40min 185min
0m 3min 188min

Then VPM-B+1: 

Depth

Duration

Runtime

Gas

120m 7min 7min (18/50)
120m 13min 20min
63m 7min 27min
63m 1min 28min
57m 1min 29min
57m 1min 30min
54m 1min 31min
54m 1min 31min
51m 1min 32min
51m 2min 33min
48m 1min 34min
48m 1min 34min
45m 1min 35min
45m 2min 36min
42m 1min 37min
42m 2min 38min
39m 1min 39min
39m 3min 41min
36m 1min 42min
36m 3min 44min
33m 1min 45min
33m 3min 47min
30m 1min 48min
30m 4min 51min
27m 1min 52min
27m 5min 56min
24m 1min 57min
24m 7min 63min
21m 1min 64min
21m 4min 67min EAN50
18m 1min 68min
18m 5min 72min
15m 1min 73min
15m 6min 78min
12m 1min 79min
12m 8min 86min
9m 1min 87min
9m 12min 98min
6m 1min 99min
6m 15min 113min EAN100
3m 3min 116min
3m 22min 138min
0m 3min 141min

Finally Bühlmann 20/125:

Depth

Duration

Runtime

Gas

120m 7min 7min (18/50)
120m 13min 20min
54m 8min 28min
54m 1min 29min
51m 1min 30min
51m 2min 31min
48m 1min 32min
48m 1min 32min
45m 1min 33min
45m 2min 34min
42m 1min 35min
42m 2min 36min
39m 1min 37min
39m 2min 38min
36m 1min 39min
36m 3min 41min
33m 1min 42min
33m 3min 44min
30m 1min 45min
30m 5min 49min
27m 1min 50min
27m 5min 54min
24m 1min 55min
24m 6min 60min
21m 1min 61min
21m 4min 64min EAN50
18m 1min 65min
18m 4min 68min
15m 1min 69min
15m 7min 75min
12m 1min 76min
12m 8min 83min
9m 1min 84min
9m 13min 96min
6m 1min 97min
6m 12min 108min EAN100
3m 3min 111min
3m 21min 132min
0m 3min 135min

Leave a Reply

Your email address will not be published. Required fields are marked *