FX-PML (Gas Discharge) vs PHAST 9 — orifice leak, choke (non-equilibrium orifice flow)
## Configuration
Pressurized-gas (vapor-phase) leak through a 2" orifice, non-equilibrium orifice choke followed by
atmospheric gas expansion.
```
ChokeModel = ChokeFlow (NonEquilibriumOrificeFlow)
GasExpansionMethod = GasExpansionMethod.ClosestToInitialConditions
ApplySonicCap = false
AirTemperature = 68 F (293.15 K)
WindSpeed = 3.36 MPH (1.502 m/s)
ReleaseHeight = 0.9144 m (3 ft)
InventoryPressure = gauge psi + 14.6959 psi atmospheric
```
- **PHAST** — Phast 9.0 results from `gas-discharge.docx` (OUTPUT *Mass flow rate*; *Final Data after
atmospheric expansion*: Temperature, Liquid mass fraction, Droplet diameter, Expanded diameter, Velocity).
- **FX-PML** — Infer `GasDischarge` with `GasExpansion` (ClosestToInitialConditions, sonic cap off),
CoolProp thermo. Regenerated by `Infer.Testing/ValidationTableHarness.cs`.
Pressurized-gas releases are single-phase vapor: liquid fraction and rainout are ~0 (Ethylene/600psi
and Propane/105psi condense a trace mist after expansion, reported by PHAST as a sub-micron droplet).
## TOC
- [[#Ethylene, gas, 68°F, 300psi, 3ft elevation, 2" hole]]
- [[#Ethylene, gas, 68°F, 600psi, 3ft elevation, 2" hole]]
- [[#Hydrogen, gas, 68°F, 300psi, 3ft elevation, 2" hole]]
- [[#Hydrogen, gas, 68°F, 600psi, 3ft elevation, 2" hole]]
- [[#Methane, gas, 68°F, 100psi, 3ft elevation, 2" hole]]
- [[#Methane, gas, 68°F, 300psi, 3ft elevation, 2" hole]]
- [[#Methane, gas, 68°F, 600psi, 3ft elevation, 2" hole]]
- [[#Propane, gas, 68°F, 105psi, 3ft elevation, 2" hole]]
### Ethylene, gas, 68°F, 300psi, 3ft elevation, 2" hole
| Variable | PHAST | FX-PML |
| --- | --- | --- |
| Total Mdot (lb/s) | 20.31710 | 20.03070 |
| Liquid Fraction | 0.000000 | 0.000000 |
| Drop Diameter (µm) | 0.000 | 0.000 |
| Texp (F) | \-148.3280 | \-146.3121 |
| Vexp (ft/s) | 1640.780 | 1644.378 |
| Expanded Diameter (ft) | 0.35256 | 0.35024 |
### Ethylene, gas, 68°F, 600psi, 3ft elevation, 2" hole
| Variable | PHAST | FX-PML |
| --- | --- | --- |
| Total Mdot (lb/s) | 44.77230 | 42.41999 |
| Liquid Fraction | 0.053688 | 0.053745 |
| Drop Diameter (µm) | 0.835 | 0.000 |
| Texp (F) | \-154.9170 | \-154.7884 |
| Vexp (ft/s) | 1487.400 | 1497.846 |
| Expanded Diameter (ft) | 0.52866 | 0.51172 |
### Hydrogen, gas, 68°F, 300psi, 3ft elevation, 2" hole
| Variable | PHAST | FX-PML |
| --- | --- | --- |
| Total Mdot (lb/s) | 5.19878 | 5.22693 |
| Liquid Fraction | 0.000000 | 0.000000 |
| Drop Diameter (µm) | 0.000 | 0.000 |
| Texp (F) | \-178.8930 | \-184.5620 |
| Vexp (ft/s) | 6456.330 | 6448.738 |
| Expanded Diameter (ft) | 0.32308 | 0.32079 |
### Hydrogen, gas, 68°F, 600psi, 3ft elevation, 2" hole
| Variable | PHAST | FX-PML |
| --- | --- | --- |
| Total Mdot (lb/s) | 10.14960 | 10.22171 |
| Liquid Fraction | 0.000000 | 0.000000 |
| Drop Diameter (µm) | 0.000 | 0.000 |
| Texp (F) | \-188.9340 | \-195.7585 |
| Vexp (ft/s) | 6607.580 | 6593.118 |
| Expanded Diameter (ft) | 0.43817 | 0.43453 |
### Methane, gas, 68°F, 100psi, 3ft elevation, 2" hole
| Variable | PHAST | FX-PML |
| --- | --- | --- |
| Total Mdot (lb/s) | 5.24138 | 5.25065 |
| Liquid Fraction | 0.000000 | 0.000000 |
| Drop Diameter (µm) | 0.000 | 0.000 |
| Texp (F) | \-116.7660 | \-116.3308 |
| Vexp (ft/s) | 2138.990 | 2139.530 |
| Expanded Diameter (ft) | 0.21984 | 0.22017 |
### Methane, gas, 68°F, 300psi, 3ft elevation, 2" hole
| Variable | PHAST | FX-PML |
| --- | --- | --- |
| Total Mdot (lb/s) | 14.79190 | 14.76825 |
| Liquid Fraction | 0.000000 | 0.000000 |
| Drop Diameter (µm) | 0.000 | 0.000 |
| Texp (F) | \-150.9910 | \-149.9998 |
| Vexp (ft/s) | 2261.440 | 2262.725 |
| Expanded Diameter (ft) | 0.34033 | 0.34054 |
### Methane, gas, 68°F, 600psi, 3ft elevation, 2" hole
| Variable | PHAST | FX-PML |
| --- | --- | --- |
| Total Mdot (lb/s) | 29.65070 | 29.46889 |
| Liquid Fraction | 0.000000 | 0.000000 |
| Drop Diameter (µm) | 0.000 | 0.000 |
| Texp (F) | \-169.1580 | \-167.4391 |
| Vexp (ft/s) | 2254.610 | 2256.371 |
| Expanded Diameter (ft) | 0.46772 | 0.46753 |
### Propane, gas, 68°F, 105psi, 3ft elevation, 2" hole
| Variable | PHAST | FX-PML |
| --- | --- | --- |
| Total Mdot (lb/s) | 9.19385 | 9.10842 |
| Liquid Fraction | 0.002659 | 0.007301 |
| Drop Diameter (µm) | 1.081 | 0.000 |
| Texp (F) | \-43.7265 | \-43.8048 |
| Vexp (ft/s) | 1265.870 | 1258.673 |
| Expanded Diameter (ft) | 0.24797 | 0.24626 |
### Parity Plots
All points are PHAST (x) vs FX-PML (y) across the eight gas cases; the red line is perfect parity.
**Total Mdot (lb/s)**
```plotly
data:
- x: [0, 45]
y: [0, 45]
mode: lines
type: scatter
name: parity
line: { color: '#e6194B' }
- x: [0, 45]
y: [0, 49.5]
mode: lines
type: scatter
name: '+10%'
line: { color: 'black', width: 1, dash: dot }
- x: [0, 45]
y: [0, 40.5]
mode: lines
type: scatter
name: '-10%'
line: { color: 'black', width: 1, dash: dot }
- x: [20.3171, 44.7723, 5.19878, 10.1496, 5.24138, 14.7919, 29.6507, 9.19385]
y: [20.0307, 42.42, 5.22693, 10.2217, 5.25065, 14.7683, 29.4689, 9.10842]
mode: markers
type: scatter
name: data
marker: { color: '#e0a800', size: 8, symbol: circle }
layout:
title: Total Mdot (lb/s)
xaxis: { title: 'PHAST 9', showline: false, mirror: false, zeroline: true, ticks: outside }
yaxis: { title: 'FX-PML', showline: false, mirror: false, zeroline: true, ticks: outside }
legend: { orientation: 'h', x: 0.5, xanchor: 'center', y: -0.2 }
margin: { t: 30 }
```
**Texp (F)**
```plotly
data:
- x: [-210, -40]
y: [-210, -40]
mode: lines
type: scatter
name: parity
line: { color: '#e6194B' }
- x: [-210, -40]
y: [-231, -44]
mode: lines
type: scatter
name: '+10%'
line: { color: 'black', width: 1, dash: dot }
- x: [-210, -40]
y: [-189, -36]
mode: lines
type: scatter
name: '-10%'
line: { color: 'black', width: 1, dash: dot }
- x: [-148.328, -154.917, -178.893, -188.934, -116.766, -150.991, -169.158, -43.7265]
y: [-146.312, -154.788, -184.562, -195.758, -116.331, -150, -167.439, -43.8048]
mode: markers
type: scatter
name: data
marker: { color: '#e0a800', size: 8, symbol: circle }
layout:
title: Texp (F)
xaxis: { title: 'PHAST 9', showline: false, mirror: false, zeroline: true, ticks: outside }
yaxis: { title: 'FX-PML', showline: false, mirror: false, zeroline: true, ticks: outside }
legend: { orientation: 'h', x: 0.5, xanchor: 'center', y: -0.2 }
margin: { t: 30 }
```
**Vexp (ft/s)**
```plotly
data:
- x: [0, 7000]
y: [0, 7000]
mode: lines
type: scatter
name: parity
line: { color: '#e6194B' }
- x: [0, 7000]
y: [0, 7700]
mode: lines
type: scatter
name: '+10%'
line: { color: 'black', width: 1, dash: dot }
- x: [0, 7000]
y: [0, 6300]
mode: lines
type: scatter
name: '-10%'
line: { color: 'black', width: 1, dash: dot }
- x: [1640.78, 1487.4, 6456.33, 6607.58, 2138.99, 2261.44, 2254.61, 1265.87]
y: [1644.38, 1497.85, 6448.74, 6593.12, 2139.53, 2262.72, 2256.37, 1258.67]
mode: markers
type: scatter
name: data
marker: { color: '#e0a800', size: 8, symbol: circle }
layout:
title: Vexp (ft/s)
xaxis: { title: 'PHAST 9', showline: false, mirror: false, zeroline: true, ticks: outside }
yaxis: { title: 'FX-PML', showline: false, mirror: false, zeroline: true, ticks: outside }
legend: { orientation: 'h', x: 0.5, xanchor: 'center', y: -0.2 }
margin: { t: 30 }
```
**Expanded Diameter (ft)**
```plotly
data:
- x: [0, 0.7]
y: [0, 0.7]
mode: lines
type: scatter
name: parity
line: { color: '#e6194B' }
- x: [0, 0.7]
y: [0, 0.77]
mode: lines
type: scatter
name: '+10%'
line: { color: 'black', width: 1, dash: dot }
- x: [0, 0.7]
y: [0, 0.63]
mode: lines
type: scatter
name: '-10%'
line: { color: 'black', width: 1, dash: dot }
- x: [0.352556, 0.528662, 0.323079, 0.438174, 0.219842, 0.340331, 0.467718, 0.247967]
y: [0.35024, 0.51172, 0.32079, 0.43453, 0.22017, 0.34054, 0.46753, 0.24626]
mode: markers
type: scatter
name: data
marker: { color: '#e0a800', size: 8, symbol: circle }
layout:
title: Expanded Diameter (ft)
xaxis: { title: 'PHAST 9', showline: false, mirror: false, zeroline: true, ticks: outside }
yaxis: { title: 'FX-PML', showline: false, mirror: false, zeroline: true, ticks: outside }
legend: { orientation: 'h', x: 0.5, xanchor: 'center', y: -0.2 }
margin: { t: 30 }
```