Skip to main content
Have a personal or library account? Click to login
Uncertainty Quantification Implementation for Quantitative Risk Assessments of Hydrogen Systems Cover

Uncertainty Quantification Implementation for Quantitative Risk Assessments of Hydrogen Systems

Open Access
|Oct 2026

Figures & Tables

Table 1

Input uncertainty distributions.

NAMEDISTRIBUTIONDESCRIPTION
valve-ii = 0.01; Log–normal (–5.9, 0.2)
i = 0.1; Log–normal (–7.4, 0.4)
i = 1.0; Log–normal (–9.8, 1.1)
i = 10; Log–normal (–10.6, 0.6)
i = 100; Log–normal (–12.2, 1.4)
Default probability distribution for valves for each leak size i ∈ {0.01,0.1,1,10,100}
instrument-ii = 0.01; Log–normal (–7.4, 0.7)
i = 0.1; Log–normal (–8.5, 0.8)
i = 1.0; Log–normal (–9.1, 0.9)
i = 10; Log–normal (–9.2, 1.1)
i = 100; Log–normal (–10.2, 1.5)
Default probability distribution for instruments for each leak size i ∈ {0.01,0.1,1,10,100}
joint-ii = 0.01; Log–normal (–10.3, 0.2)
i = 0.1; Log–normal (–12.3, 0.9)
i = 1.0; Log–normal (–11.8, 0.5)
i = 10; Log–normal (–11.8, 0.6)
i = 100; Log–normal (–12, 0.7)
Default probability distribution for joints for each leak size i ∈ {0.01,0.1,1,10,100}
hose-ii = 0.01; Log–normal (–7.5, 0.4)
i = 0.1; Log–normal (–8.5, 0.6)
i = 1.0; Log–normal (–8.7, 0.6)
i = 10; Log–normal (–8.8, 0.6)
i = 100; Log–normal (–9.7, 1.0)
Default probability distribution for hoses for each leak size i ∈ {0.01,0.1,1,10,100}
pipe-ii = 0.01; Log–normal (–11.7, 0.7)
i = 0.1; Log–normal (–12.5, 0.7)
i = 1.0; Log–normal (–13.9, 1.3)
i = 10; Log–normal (–14.6, 1.2)
i = 100; Log–normal (–15.7, 1.8)
Default probability distribution for pipes for each leak size i ∈ {0.01,0.1,1,10,100}. Every 1 meter of pipe is considered a separate pipe component.
noz_po_distBeta (0.5, 610415.5)Default probability distribution for nozzle pop-off (disconnection)
overp_distBeta (3.5, 310289.5)Default probability distribution for an overpressure during fueling
pvalve_fto_distLog–normal (–11.74, 0.67)Default probability distribution for pressure relief valve failure to open
driveoff_distBeta (31.5, 610384.5)Default probability distribution for accidents due to drive-off during refueling
coupling_ftc_distBeta (0.5, 5031)Default probability distribution for breakaway coupling failure to close
x_i_locationUniform (1, 20) [m]Independently sampled for each occupant i ∈ {1,2,…N}
z_i_locationUniform (1, 12) [m]Independently sampled for each occupant I ∈ {1,2,…N}
Figure 1

Illustration of the Gaussian Process model concept: for data points near each other in the input space, there is less uncertainty, which results in less variation between functions. For interpolations between distant points, there is more variation.

Table 2

Descriptions of the different cases included in the sensitivity analysis.

SAMPLE STRUCTURECASE 1: BASE CASECASE 2: OCCUPANTS MOVED CLOSER TO LEAKSCASE 3: DOUBLING OF SYSTEM PRESSURECASE 4: OCCUPANTS MOVED CLOSER TO LEAK AND DOUBLING OF SYSTEM PRESSURE
SINGLE LOOPDUAL LOOPSINGLE LOOPDUAL LOOPSINGLE LOOPDUAL LOOPSINGLE LOOPDUAL LOOP
Number of Occupants19191919
Number of Samples2000N = 125
M = 40
2000N = 125
M = 40
2000N = 125
M = 40
2000N = 125
M = 40
X Loc. ~ uniform [m]1–200.01–11–200.01–1
Z-Loc. ~ uniform [m]1–120.01–0.021–120.01–0.02
Pressure [MPa]35357070

[i] Note: N is the sample size for the location sampling loop; M is the sample size for the non-location parameter loop.

Figure 2

Uncertainty quantification results generated by propagating 1,000 LHS samples through HyRAM+ QRA calculations. The left plot shows a histogram of the PLL, while the right plot shows the empirical cumulative distribution frequency. Median, 5th percentile, and deterministic HyRAM+ results are highlighted on both plots.

Figure 3

Comparison between the empirical CDFs of the simulation results and the regression predictions for the Case 1 single loop analysis.

Figure 4

Main (top) and total (bottom) effect indices for the Case 1 single-loop analysis identify the instrument rupture frequency and location of the occupant relative to the leak as the most important factors affecting PLL.

Figure 5

Comparison of the simulation results and regression results from the Case 1 single-loop analysis for the instrument rupture frequency (left) and z-location (right) versus the standardized PLL.

Figure 6

Comparison between the empirical cumulative distribution functions of the simulation results and the surrogate predictions for the Case 1 dual-loop analysis indicates a good fit in the regression analysis for PLL.

Figure 7

Comparison of the simulation results and regression results from the Case 1 dual-loop analysis for the standardized instrument rupture frequency (left) and standardized joint rupture frequency (right) versus the standardized PLL.

Figure 8

Total effect sensitivity indices for the Case 1 dual-loop analysis.

Figure 9

Total effect sensitivity indices for the Case 2 dual-loop analysis with different colors denoting different occupant location realizations.

Language: English
Page range: 317 - 330
Submitted on: May 6, 2026
Accepted on: Oct 15, 2026
Published on: Oct 30, 2026
In partnership with: Paradigm Publishing Services

© 2026 Dusty Brooks, Benjamin Schroeder, Brian Ehrhart, published by KIT Scientific Publishing
This work is licensed under the Creative Commons Attribution 4.0 License.