lh2poolx 0.1.0__py3-none-any.whl

lh2poolx/__init__.py

@@ -0,0 +1,19 @@
+"""Evidence-qualified LH2 pool source terms for dispersion-model inputs."""
+
+from .pool import (
+    CONCRETE_CRYOGENIC,
+    CRITICAL_HEAT_FLUX_W_M2,
+    LH2Release,
+    PoolSource,
+    Substrate,
+    evaluate_pool_source,
+    flash_vapour_fraction,
+)
+
+__version__ = "0.1.0"
+
+__all__ = [
+    "__version__", "CONCRETE_CRYOGENIC", "CRITICAL_HEAT_FLUX_W_M2",
+    "LH2Release", "PoolSource", "Substrate", "evaluate_pool_source",
+    "flash_vapour_fraction",
+]

lh2poolx/pool.py

@@ -0,0 +1,148 @@
+"""Quasi-steady LH2 pool source terms with explicit evidence boundaries."""
+
+from __future__ import annotations
+
+from dataclasses import dataclass
+import math
+
+import CoolProp.CoolProp as CP
+
+
+CRITICAL_HEAT_FLUX_W_M2 = 120e3
+
+
+@dataclass(frozen=True)
+class Substrate:
+    """Effective cryogenic substrate properties for the stated evidence case."""
+
+    name: str
+    conductivity_W_mK: float
+    diffusivity_m2_s: float
+
+
+# PRESLHY E3.4 Concrete02 effective cryogenic properties.  Do not use this as
+# a generic concrete correlation without an explicit applicability decision.
+CONCRETE_CRYOGENIC = Substrate("concrete_cryogenic_e34", 2.0, 2.5e-7)
+
+
+@dataclass(frozen=True)
+class LH2Release:
+    """Declared inputs for an LH2 release reaching a possible ground pool."""
+
+    rate_kg_s: float
+    storage_pressure_barg: float
+    deposition_fraction: float = 1.0
+    max_radius_m: float | None = None
+
+    def __post_init__(self) -> None:
+        if self.rate_kg_s <= 0.0:
+            raise ValueError("rate_kg_s must be > 0")
+        if self.storage_pressure_barg < 0.0:
+            raise ValueError("storage_pressure_barg must be >= 0")
+        if not 0.0 <= self.deposition_fraction <= 1.0:
+            raise ValueError("deposition_fraction must be in [0, 1]")
+        if self.max_radius_m is not None and self.max_radius_m <= 0.0:
+            raise ValueError("max_radius_m must be > 0 when supplied")
+
+
+@dataclass(frozen=True)
+class PoolSource:
+    """A source term for one declared time window, not a pool trajectory."""
+
+    elapsed_s: float
+    flash_vapour_fraction: float
+    deposition_fraction: float
+    liquid_to_ground_kg_s: float
+    pool_temperature_K: float
+    latent_heat_J_kg: float
+    evaporative_flux_kg_m2_s: float
+    unconfined_radius_m: float
+    radius_m: float
+    area_m2: float
+    evaporation_rate_kg_s: float
+    liquid_accumulation_rate_kg_s: float
+    confined: bool
+    source_status: str
+
+
+def flash_vapour_fraction(storage_pressure_barg: float, *,
+                          ambient_pressure_Pa: float = 101325.0) -> float:
+    """Isenthalpic saturated-liquid LH2 flash fraction at ambient pressure."""
+    if storage_pressure_barg < 0.0:
+        raise ValueError("storage_pressure_barg must be >= 0")
+    p0 = storage_pressure_barg * 1e5 + ambient_pressure_Pa
+    if p0 >= CP.PropsSI("Pcrit", "Hydrogen"):
+        raise ValueError("storage pressure is at or above H2 critical pressure")
+    h0 = CP.PropsSI("H", "P", p0, "Q", 0, "Hydrogen")
+    x = CP.PropsSI("Q", "P", ambient_pressure_Pa, "H", h0, "Hydrogen")
+    return min(max(float(x), 0.0), 1.0)
+
+
+def _ground_flux(*, substrate: Substrate, latent_heat_J_kg: float,
+                 ambient_temperature_K: float, pool_temperature_K: float,
+                 elapsed_s: float,
+                 critical_heat_flux_W_m2: float | None) -> float:
+    if elapsed_s <= 0.0:
+        raise ValueError("elapsed_s must be > 0")
+    if latent_heat_J_kg <= 0.0:
+        raise ValueError("latent_heat_J_kg must be > 0")
+    if critical_heat_flux_W_m2 is not None and critical_heat_flux_W_m2 <= 0.0:
+        raise ValueError("critical_heat_flux_W_m2 must be > 0 or None")
+    delta_T = ambient_temperature_K - pool_temperature_K
+    if delta_T <= 0.0:
+        return 0.0
+    q = substrate.conductivity_W_mK * delta_T / math.sqrt(
+        math.pi * substrate.diffusivity_m2_s * elapsed_s
+    )
+    if critical_heat_flux_W_m2 is not None:
+        q = min(q, critical_heat_flux_W_m2)
+    return q / latent_heat_J_kg
+
+
+def evaluate_pool_source(release: LH2Release, *, elapsed_s: float,
+                         substrate: Substrate = CONCRETE_CRYOGENIC,
+                         ambient_temperature_K: float = 282.0,
+                         ambient_pressure_Pa: float = 101325.0,
+                         critical_heat_flux_W_m2: float | None =
+                         CRITICAL_HEAT_FLUX_W_M2) -> PoolSource:
+    """Calculate a quasi-steady pool source for one elapsed-time window.
+
+    The pool area instantaneously balances the declared liquid-to-ground rate
+    against heat-limited evaporation.  This is appropriate only after a pool
+    exists and where a quasi-steady source approximation is justified.
+    """
+    x = flash_vapour_fraction(release.storage_pressure_barg,
+                              ambient_pressure_Pa=ambient_pressure_Pa)
+    T_pool = float(CP.PropsSI("T", "P", ambient_pressure_Pa, "Q", 0,
+                              "Hydrogen"))
+    latent = float(CP.PropsSI("H", "P", ambient_pressure_Pa, "Q", 1,
+                              "Hydrogen") - CP.PropsSI(
+                                  "H", "P", ambient_pressure_Pa, "Q", 0,
+                                  "Hydrogen"))
+    flux = _ground_flux(
+        substrate=substrate, latent_heat_J_kg=latent,
+        ambient_temperature_K=ambient_temperature_K,
+        pool_temperature_K=T_pool, elapsed_s=elapsed_s,
+        critical_heat_flux_W_m2=critical_heat_flux_W_m2,
+    )
+    liquid_to_ground = release.rate_kg_s * (1.0 - x) * release.deposition_fraction
+    radius_unconfined = (0.0 if liquid_to_ground == 0.0 else math.sqrt(
+        liquid_to_ground / (math.pi * flux)))
+    radius = radius_unconfined if release.max_radius_m is None else min(
+        radius_unconfined, release.max_radius_m)
+    area = math.pi * radius * radius
+    evaporation = area * flux
+    confined = (release.max_radius_m is not None and
+                 radius_unconfined > release.max_radius_m)
+    return PoolSource(
+        elapsed_s=elapsed_s, flash_vapour_fraction=x,
+        deposition_fraction=release.deposition_fraction,
+        liquid_to_ground_kg_s=liquid_to_ground, pool_temperature_K=T_pool,
+        latent_heat_J_kg=latent, evaporative_flux_kg_m2_s=flux,
+        unconfined_radius_m=radius_unconfined, radius_m=radius, area_m2=area,
+        evaporation_rate_kg_s=evaporation,
+        liquid_accumulation_rate_kg_s=max(liquid_to_ground - evaporation, 0.0),
+        confined=confined,
+        source_status=("quasi_steady_equilibrium" if not confined else
+                       "confined_pool_requires_inventory_model"),
+    )

lh2poolx-0.1.0.dist-info/METADATA

@@ -0,0 +1,92 @@
+Metadata-Version: 2.5
+Name: lh2poolx
+Version: 0.1.0
+Summary: Transparent liquid-hydrogen pool source terms for dispersion models
+Project-URL: Source, https://github.com/lyullee/lh2poolx
+Project-URL: Documentation, https://github.com/lyullee/lh2poolx#readme
+Author: Ugwiyeon Lee
+License: MIT License
+        
+        Copyright (c) 2026 Ugwiyeon Lee
+        
+        Permission is hereby granted, free of charge, to any person obtaining a copy
+        of this software and associated documentation files (the "Software"), to deal
+        in the Software without restriction, including without limitation the rights
+        to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+        copies of the Software, and to permit persons to whom the Software is
+        furnished to do so, subject to the following conditions:
+        
+        The above copyright notice and this permission notice shall be included in all
+        copies or substantial portions of the Software.
+        
+        THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+        IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+        FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+        AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+        LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+        OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+        SOFTWARE.
+License-File: LICENSE
+Keywords: LH2,liquid hydrogen,pool evaporation,process safety,source term
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Programming Language :: Python :: 3
+Classifier: Topic :: Scientific/Engineering
+Requires-Python: >=3.10
+Requires-Dist: coolprop>=6.5
+Provides-Extra: test
+Requires-Dist: pytest>=7; extra == 'test'
+Description-Content-Type: text/markdown
+
+# LH2PoolX
+
+LH2PoolX calculates an **evidence-qualified, quasi-steady liquid-hydrogen
+pool source term** for a declared time window.  It is an upstream input module
+for dispersion models, not a dispersion model and not a general jet-impact or
+pool-spreading solver.
+
+## What it calculates
+
+For a specified release rate, storage pressure, elapsed time, substrate and
+ground-deposition fraction, LH2PoolX calculates:
+
+1. isenthalpic flash fraction using CoolProp;
+2. liquid mass reaching the ground;
+3. heat-limited evaporation flux using a semi-infinite ground-conduction
+   relation with an LH2 critical-heat-flux cap; and
+4. the unconfined equilibrium pool area and evaporation rate.
+
+`deposition_fraction` is deliberately explicit.  The package does **not**
+predict impingement, splashing, droplet transport, drainage, barriers, or
+transient pool spreading.  A confined result reports liquid accumulation and
+must not be used as a steady source without a separate inventory model.
+
+## Install
+
+```bash
+pip install lh2poolx
+```
+
+## Minimal use
+
+```python
+from lh2poolx import LH2Release, evaluate_pool_source
+
+release = LH2Release(rate_kg_s=0.1055, storage_pressure_barg=1.0)
+source = evaluate_pool_source(release, elapsed_s=300.0)
+print(source.area_m2, source.evaporation_rate_kg_s)
+```
+
+## Evidence boundary
+
+The concrete heat-transfer implementation is checked against the four
+published PRESLHY E3.4 Concrete02 mass-loss windows.  That comparison is a
+restricted substrate-and-condition check, not validation for other substrates
+or for jet-to-ground deposition.  Literature supports the physical form of
+the source chain; it does not turn unknown deposition into a fitted value.
+
+## License
+
+MIT.  Cite the accompanying software record when one is released, plus the
+underlying experimental and heat-transfer sources used in an analysis.

lh2poolx-0.1.0.dist-info/RECORD

@@ -0,0 +1,6 @@
+lh2poolx/__init__.py,sha256=UwGw6FCZbN7Wd1EIx5pSktT8StTnPllPgwewDCEriEc,456
+lh2poolx/pool.py,sha256=AexRFuWbCQmJMyfIy1Hzc07WY2Vm9D5C7YEwus9e3VA,6078
+lh2poolx-0.1.0.dist-info/METADATA,sha256=ngOmnjJYaevbSXBlIqAjLb5GiYLKOKIXLUw-S6OEsCo,3850
+lh2poolx-0.1.0.dist-info/WHEEL,sha256=zOwg4jB6zX2kU910N-cMawjivD6tO8NEWvE12je1bVk,87
+lh2poolx-0.1.0.dist-info/licenses/LICENSE,sha256=1NqdQIjRYLmcIHzB2vXntniTEiavLqsAoyiAabNdWtQ,1069
+lh2poolx-0.1.0.dist-info/RECORD,,

lh2poolx-0.1.0.dist-info/WHEEL

@@ -0,0 +1,4 @@
+Wheel-Version: 1.0
+Generator: hatchling 1.32.0
+Root-Is-Purelib: true
+Tag: py3-none-any

lh2poolx-0.1.0.dist-info/licenses/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Ugwiyeon Lee
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.