api581-engine 0.1.4__py3-none-any.whl

api581/hydrofluoric_acid.py

@@ -0,0 +1,236 @@
+"""2016 HF tables with explicit concentration, residual and source-conflict policies."""
+from bisect import bisect_left
+from typing import Literal
+from pydantic import Field
+from .core import Model, NonNegative, Positive, Fraction, calculation, ensure
+from .thinning import Posterior, WallLoss
+from .clad_thinning import PressureStrength, CodeStrength
+from .thinning_conditions import ThinningConditions, QualifiedUnclad, qualified_unclad
+
+Units = Literal["printed_si", "printed_us_converted"]
+Lookup = Literal["exact_nodes_only", "bilinear", "maximum_bracketing_value"]
+_C = (.01,.02,.05,.06,.63,.64,.80,.81)
+# Paired low/high velocity rows; last two columns are81wt% low/high residual.
+_CS_US = (
+    (2,150,150,800,800,5,5,2,6),(20,999,999,999,999,50,50,20,60),
+    (10,500,500,999,999,30,30,5,15),(200,999,999,999,999,300,300,50,150),
+    (10,500,500,999,999,30,30,10,30),(100,999,999,999,999,300,300,100,300),
+    (100,999,999,999,999,500,500,20,60),(999,999,999,999,999,999,999,200,600),
+    (100,999,999,999,999,500,500,50,150),(999,999,999,999,999,999,999,500,999),
+    (100,999,999,999,999,500,500,70,210),(999,999,999,999,999,999,999,700,999),
+    (500,999,999,999,999,999,999,100,300),(999,999,999,999,999,999,999,999,999))
+_CS_SI = (
+    (.05,3.81,3.81,20.32,20.32,.13,.13,.05,.15),(.51,25.37,25.37,25.37,25.37,1.27,1.27,.51,1.52),
+    (.25,12.7,12.7,25.37,25.37,.76,.76,.13,.38),(5.08,25.37,25.37,25.37,25.37,7.62,7.62,1.27,3.81),
+    (.25,12.7,12.7,25.37,25.37,.76,.76,.25,.76),(2.54,25.37,25.37,25.37,25.37,7.62,7.62,2.54,7.62),
+    (2.54,25.37,25.37,25.37,25.37,12.7,12.7,.51,1.52),(25.37,25.37,25.37,25.37,25.37,25.37,25.37,5.08,15.24),
+    (2.54,25.37,25.37,25.37,25.37,12.7,12.7,1.27,3.81),(25.37,25.37,25.37,25.37,25.37,25.37,25.37,12.7,25.37),
+    (2.54,25.37,25.37,25.37,25.37,12.7,12.7,1.78,5.33),(25.37,25.37,25.37,25.37,25.37,25.37,25.37,17.78,25.37),
+    (12.7,25.37,25.37,25.37,25.37,25.37,25.37,2.54,7.62),(25.37,25.37,25.37,25.37,25.37,25.37,25.37,25.37,25.37))
+# Paired no/yes aeration rows.
+_A400_US = ((1,1,1,10,10,1,1,2),(10,10,10,25,25,10,10,15),
+    (1,1,1,15,15,5,5,3),(10,10,10,30,30,20,20,15),
+    (5,5,5,20,20,10,10,5),(20,20,20,100,100,50,50,20),
+    (10,10,10,20,20,20,20,10),(100,100,100,200,200,200,200,100))
+_A400_SI = ((.03,.03,.03,.25,.25,.03,.03,.05),(.25,.25,.25,.64,.64,.25,.25,.38),
+    (.03,.03,.03,.38,.38,.13,.13,.08),(.25,.25,.25,.76,.76,.51,.51,.38),
+    (.13,.13,.13,.51,.51,.25,.25,.13),(.51,.51,.51,2.54,2.54,1.27,1.27,.51),
+    (.25,.25,.25,.51,.51,.51,.51,.25),(2.54,2.54,2.54,5.08,5.08,5.08,5.08,2.54))
+
+
+def _qualified(value,message):
+    ensure(bool(value.strip()),message)
+
+
+class HFConcentration(Model):
+    hf_mass_fraction: Fraction
+    concentrated_extension_policy: Literal["reject_above_81wt_percent", "qualified_81wt_percent_node"]
+    concentrated_extension_qualification_basis: str | None = None
+
+
+@calculation("damage/hf-concentration-selection", "HF table concentration with explicit above81wt% extension policy", HFConcentration,
+    ["Part2 Annex2.B §6.1,Tables6.2-6.3M; PDF pp374,376-377"],limitations=[
+        "HF-in-water mass fraction, not water mass fraction. Source body says>=80wt% equivalent to anhydrous but tables distinguish80/81wt% and stop81. Values>81 default reject; qualified81-node extension requires owner basis and flag, never silent extrapolation or81>= label. Below1wt% rejected.80..81 handled by explicit table interpolation policy. No-free-water is not HF immunity: anhydrous acid remains in scope."])
+def concentration_selection(x):
+    qualified = x.concentrated_extension_policy == "qualified_81wt_percent_node"
+    ensure((x.concentrated_extension_qualification_basis is not None and bool(x.concentrated_extension_qualification_basis.strip())) if qualified else x.concentrated_extension_qualification_basis is None,
+        "Concentrated HF extension policy and qualification basis must agree")
+    ensure(x.hf_mass_fraction >= .01,"HF concentration below first1wt% table node")
+    extended = x.hf_mass_fraction > .81
+    ensure(not extended or qualified,"Above81wt% HF requires qualified source-body/table extension","hf_concentrated_extension_required")
+    return {"hf_mass_fraction_for_lookup":min(x.hf_mass_fraction,.81),"concentrated_81_node_extension_used":extended}
+
+
+class HFVelocity(Model):
+    velocity_m_s: NonNegative
+    velocity_basis: Units
+
+
+@calculation("damage/hf-carbon-velocity-band", "HF carbon-steel low/high velocity band", HFVelocity,
+    ["Part2 Tables2.B.6.2/6.2M; PDF p376"],limitations=[
+        "PrintedSI threshold3.05m/s versus10ft/s=3.048m/s explicitly selected. Low strictly below; equality high. No continuous velocity interpolation or generic erosion/film-stripping rate. This rate-table threshold differs from other mechanisms and is not a universal process limit."])
+def velocity_band(x):
+    threshold = 3.05 if x.velocity_basis == "printed_si" else 10*.3048
+    return {"velocity_band": "low" if x.velocity_m_s < threshold else "high","velocity_threshold_m_s":threshold}
+
+
+def _bracket(nodes,value):
+    value = next((n for n in nodes if abs(n-value) <= 1e-10),value)
+    ensure(nodes[0] <= value <= nodes[-1],"Outside printed HF table; no extrapolation")
+    i = bisect_left(nodes,value)
+    return (i,i,0.) if nodes[i] == value else (i-1,i,(value-nodes[i-1])/(nodes[i]-nodes[i-1]))
+
+
+def _lookup(table,us,temps,c,temp,branch,units,policy,residual_high=False):
+    a,b,w = _bracket(_C,c)
+    i,j,v = _bracket(temps,temp)
+    ensure(policy != "exact_nodes_only" or (a == b and i == j),"Exact-node HF policy does not interpolate")
+    def column(k):
+        return 8 if residual_high and k == 7 else k
+    corners = [table[2*t+branch][column(k)] for t in (i,j) for k in (a,b)]
+    if policy == "maximum_bracketing_value":
+        native = max(corners)
+    else:
+        lower = (1-w)*table[2*i+branch][column(a)]+w*table[2*i+branch][column(b)]
+        upper = (1-w)*table[2*j+branch][column(a)]+w*table[2*j+branch][column(b)]
+        native = (1-v)*lower+v*upper
+    return {"corrosion_rate_m_year":native*(.001 if units == "printed_si" else .0000254),
+        "source_999_corner_used":any(us[2*t+branch][column(k)] == 999 for t in (i,j) for k in (a,b))}
+
+
+class HFCarbonRate(HFConcentration,HFVelocity):
+    maximum_temperature_k: Positive
+    residual_category: Literal["low", "high"] | None = None
+    residual_qualification_basis: str | None = None
+    unit_basis: Units
+    lookup_policy: Lookup
+
+
+@calculation("damage/hf-carbon-rate", "HF carbon-steel concentration temperature velocity residual rate table", HFCarbonRate,
+    ["Part2 Tables2.B.6.2/6.2M; PDF p376"],limitations=[
+        "252SI/US cells at7temperatures x2velocity x9columns. Concentrations1,2,5,6,63,64,80,81wt%;81haslow/high residual. Any lookup using81corner needs qualified residual classification; source supplies no numeric Cu/Ni/Cr cutoff, none inferred. Residual fields outside81-corner region rejected as unused. All999mpy/25.37mm/y literal with flag. Source SI/US rounding explicit, no temperature extrapolation. Initial high corrosion before protective fluoride scale, galvanic coupling and disturbed scale not quantified by table. Figure6.1 swapped table numbers resolved by actual table titles, not treated as official erratum."])
+def carbon_rate(x):
+    concentration = concentration_selection(x)
+    c = concentration["hf_mass_fraction_for_lookup"]
+    a,b,_ = _bracket(_C,c)
+    residual_used = b == 7
+    ensure((x.residual_category is not None and x.residual_qualification_basis is not None and bool(x.residual_qualification_basis.strip())) if residual_used else x.residual_category is None and x.residual_qualification_basis is None,
+        "Qualified residual category needed only when lookup touches81wt% column")
+    velocity = velocity_band(x)
+    temps = tuple(t+273.15 for t in (27,41,60,66,75,86,93)) if x.unit_basis == "printed_si" else tuple((t-32)*5/9+273.15 for t in (80,105,140,150,167.5,187.5,200))
+    table = _CS_SI if x.unit_basis == "printed_si" else _CS_US
+    return {**concentration,**velocity,**_lookup(table,_CS_US,temps,c,x.maximum_temperature_k,int(velocity["velocity_band"] == "high"),x.unit_basis,x.lookup_policy,x.residual_category == "high"),
+        "residual_column_used":x.residual_category if residual_used else None,"source_table_id":"2.B.6.2"+("M" if x.unit_basis == "printed_si" else "")}
+
+
+class HFAlloy400Rate(HFConcentration):
+    maximum_temperature_k: Positive
+    aerated_or_oxidants_present: bool
+    unit_basis: Units
+    lookup_policy: Lookup
+
+
+@calculation("damage/hf-alloy400-rate", "HF Alloy400 concentration temperature aeration rate table", HFAlloy400Rate,
+    ["Part2 Tables2.B.6.3/6.3M,Table6.1; PDF pp375,377"],limitations=[
+        "128SI/US cells at4temperatures x2aeration x8HF concentrations. Explicit aeration/oxidants bool, no dissolved-oxygen cutoff supplied or inferred. No blanket aeration multiplier or velocity correction. Actual material and applicable solution conditions external. No extrapolation except flagged qualified above81 extension. Figure6.1 swapped table references resolved using actual Alloy400 table title."])
+def alloy400_rate(x):
+    concentration = concentration_selection(x)
+    temps = tuple(t+273.15 for t in (27,52,79,93)) if x.unit_basis == "printed_si" else tuple((t-32)*5/9+273.15 for t in (80,125,175,200))
+    table = _A400_SI if x.unit_basis == "printed_si" else _A400_US
+    return {**concentration,**_lookup(table,_A400_US,temps,concentration["hf_mass_fraction_for_lookup"],x.maximum_temperature_k,
+        int(x.aerated_or_oxidants_present),x.unit_basis,x.lookup_policy),"source_table_id":"2.B.6.3"+("M" if x.unit_basis == "printed_si" else "")}
+
+
+class HFCarbonService(HFCarbonRate):
+    method: Literal["carbon_steel_table"]
+    established_scale_and_stable_operation_qualified: bool
+    material_residual_content_qualified: bool
+    qualification_basis: str = Field(min_length=1)
+
+
+class HFAlloy400Service(HFAlloy400Rate):
+    method: Literal["alloy400_table"]
+    material_and_aeration_qualified: bool
+    qualification_basis: str = Field(min_length=1)
+
+
+class HFPublishedRate(Model):
+    method: Literal["qualified_published_rate"]
+    material_description: str = Field(min_length=1)
+    hf_mass_fraction: Fraction
+    maximum_temperature_k: Positive
+    velocity_m_s: NonNegative
+    corrosion_rate_m_year: NonNegative
+    published_source_reference: str = Field(min_length=1)
+    qualification_basis: str = Field(min_length=1)
+
+
+class HFAssessment(Model):
+    hf_exposure_present: bool
+    applicability_basis: str = Field(min_length=1)
+    rate_model: HFCarbonService | HFAlloy400Service | HFPublishedRate | None = Field(default=None,discriminator="method")
+    galvanic_coupling_present: bool | None = None
+    thinning_classification_policy: Literal["overview_local", "detailed_general_without_galvanic"] | None = None
+    classification_qualification_basis: str | None = None
+
+
+@calculation("damage/hf-assessment", "Qualified HF exposure material and source-classification corrosion assessment", HFAssessment,
+    ["Part2 Annex2.B §6,Tables6.1-6.3M,Figure6.1,Table1.2; PDF pp323,374-378"],origin="api581_composition_with_explicit_source_policies",limitations=[
+        "HF exposure, not free-water presence controls branch: no free water cannot imply zero because anhydrous HF corrodes. No-HF0 source-specific only. CS table requires stable established fluoride scale/material-residual qualification; table paths reject galvanic coupling whose acceleration is unquantified. Other materials/conditions may use explicit published rate with provenance and applicability qualification per source figure, not an API581 table prediction. Source overviewLocal versus detailedGeneral absent galvanic conflict requires explicit policy/basis; general policy disallowed with galvanic coupling. No HSC/HIC/blistering, SCC, film growth or transient dilution model."])
+def assessment(x):
+    _qualified(x.applicability_basis,"HF exposure and service applicability qualification required")
+    if not x.hf_exposure_present:
+        ensure(x.rate_model is None and x.galvanic_coupling_present is None and x.thinning_classification_policy is None and x.classification_qualification_basis is None,
+            "No-HF branch cannot carry unused service model or classification")
+        return {"corrosion_rate_m_year":0.,"source_table_id":None,"rate_source":None,"source_999_corner_used":False,
+            "hf_mass_fraction_for_lookup":None,"concentrated_81_node_extension_used":False,"velocity_band":None,
+            "velocity_threshold_m_s":None,"residual_column_used":None,"source_thinning_type":None,"calculation_branch":"no_hf_exposure"}
+    ensure(x.rate_model is not None and x.galvanic_coupling_present is not None and x.thinning_classification_policy is not None,
+        "HF service requires model, galvanic review and explicit source classification")
+    ensure(x.classification_qualification_basis is not None and bool(x.classification_qualification_basis.strip()),"Source classification conflict requires qualified policy basis")
+    general = x.thinning_classification_policy == "detailed_general_without_galvanic"
+    ensure(not general or not x.galvanic_coupling_present,"Detailed general-thinning branch excludes potential galvanic attack")
+    m = x.rate_model
+    _qualified(m.qualification_basis,"Material/conditions and source applicability qualification required")
+    if m.method == "carbon_steel_table":
+        ensure(m.established_scale_and_stable_operation_qualified and m.material_residual_content_qualified,
+            "Carbon table requires qualified established scale, stable service and residual metallurgy")
+        ensure(not x.galvanic_coupling_present,"HF table does not quantify galvanic acceleration")
+        result = carbon_rate(m)
+    elif m.method == "alloy400_table":
+        ensure(m.material_and_aeration_qualified,"Alloy400 material and aeration qualification required")
+        ensure(not x.galvanic_coupling_present,"HF table does not quantify galvanic acceleration")
+        result = {**alloy400_rate(m),"velocity_band":None,"velocity_threshold_m_s":None,"residual_column_used":None}
+    else:
+        _qualified(m.published_source_reference,"Published HF rate reference required")
+        _qualified(m.material_description,"Published rate material identification required")
+        ensure(m.hf_mass_fraction > 0,"Exposed HF service cannot use zero HF concentration")
+        result = {"corrosion_rate_m_year":m.corrosion_rate_m_year,"source_table_id":None,"source_999_corner_used":False,
+            "hf_mass_fraction_for_lookup":None,"concentrated_81_node_extension_used":False,"velocity_band":None,
+            "velocity_threshold_m_s":None,"residual_column_used":None}
+    return {**result,"rate_source":m.method,"source_thinning_type":"general" if general else "local","calculation_branch":"qualified_hf_service"}
+
+
+class HFThinning(Model):
+    corrosion: HFAssessment
+    assessment_distribution: Literal["general","local"]
+    distribution_qualification_basis: str = Field(min_length=1)
+    starting_thickness_m: Positive
+    elapsed_years: NonNegative
+    strength: PressureStrength | CodeStrength = Field(discriminator="method")
+    inspections: Posterior
+    conditions: ThinningConditions
+    qualification_basis: str = Field(min_length=1)
+
+
+@calculation("damage/hf-thinning-assessment", "Qualified HF rate through non-AST unclad thinning DF", HFThinning,
+    ["Part2 Annex2.B §6,Table1.2; §4.5.7,Table4.8"],origin="api581_composition_with_explicit_source_policies",limitations=[
+        "Same selected HF rate drives wall loss and qualified HF monitoring row. Source class conflict policy and owner-qualified actual distribution retained separately. No-HF0 not DF0, uncertainty/floor remain. Non-AST unclad only; no added residual/aeration multipliers, no automatic cracking/totalDF/FFS acceptance. Published rate retains supplied provenance in input snapshot, not relabeled API581 table data."])
+def thinning_assessment(x):
+    _qualified(x.distribution_qualification_basis,"Actual general/local inspection distribution qualification required")
+    corrosion = assessment(x.corrosion)
+    ensure(x.conditions.monitoring.row.mechanism == "hydrofluoric_acid","HF calculation requires hydrofluoric-acid monitoring row")
+    result = qualified_unclad(QualifiedUnclad(assessment_scope="non_ast_unclad",qualification_basis=x.qualification_basis,
+        wall=WallLoss(starting_thickness_m=x.starting_thickness_m,corrosion_rate_m_year=corrosion["corrosion_rate_m_year"],elapsed_years=x.elapsed_years),
+        strength=x.strength,inspections=x.inspections,conditions=x.conditions))
+    return {**corrosion,**result,"assessment_distribution":x.assessment_distribution}

api581/ignition.py

@@ -0,0 +1,204 @@
+"""Level2 ignition correlations, explicit temperature reference, figure trees."""
+import math
+from typing import Literal
+from pydantic import Field
+from .core import Model, NonNegative, Positive, Fraction, calculation, ensure
+
+Phase=Literal["liquid","vapor"]
+ReleaseType=Literal["continuous","instantaneous"]
+_COEFFICIENTS={
+    "liquid":(-.1368924,-.007598764,8.282163e-6,-6.124231e-9,.06876128,1.193736e-4,2.081034e-7,-4.057289e-11),
+    "vapor":(-.06053124,-.009958413,1.518603e-5,-1.386705e-8,.04564953,7.912392e-4,-6.489157e-7,7.159409e-10),
+}
+_TABLE={
+    ("continuous","liquid"):(.2,.25), ("continuous","vapor"):(.5,.5),
+    ("instantaneous","liquid"):(.2,.125), ("instantaneous","vapor"):(.1,.25),
+}
+
+
+class FlammableRates(Model):
+    adjusted_mass_rate_kg_s: NonNegative
+    flammable_mass_fraction: Fraction
+    flash_mass_fraction: Fraction
+
+
+@calculation("release/flammable-phase-rates", "Flammable total/liquid/vapor rates for ignition", FlammableRates,
+             ["Part3 eqs(3.106)-(3.108); PDF p546"],
+             limitations=["Flash fraction supplied by qualified thermodynamics. Phase liquid rate is not pool rainout rate; aerosol/rainout partition is separate. Flammable mass fraction, not mole fraction, multiplies total adjusted rate."])
+def flammable_rates(x):
+    total=x.adjusted_mass_rate_kg_s*x.flammable_mass_fraction
+    vapor=total*x.flash_mass_fraction
+    liquid=total*(1-x.flash_mass_fraction)
+    ensure(total > 0 or x.adjusted_mass_rate_kg_s == 0 or x.flammable_mass_fraction == 0,
+           "Positive flammable rate underflow","numerical_failure")
+    ensure(vapor > 0 or total == 0 or x.flash_mass_fraction == 0,
+           "Positive flammable vapor rate underflow","numerical_failure")
+    ensure(liquid > 0 or total == 0 or x.flash_mass_fraction == 1,
+           "Positive flammable liquid rate underflow","numerical_failure")
+    return {"flammable_rate_kg_s":total,"flammable_vapor_rate_kg_s":vapor,"flammable_liquid_rate_kg_s":liquid}
+
+
+class AmbientIgnition(Model):
+    phase: Phase
+    flammable_rate_kg_s: NonNegative
+    flash_point_temperature_k: Positive
+    correlation_basis_source: str = Field(min_length=1,max_length=1000)
+
+
+@calculation("events/ambient-ignition", "Flash-point/rate ambient ignition correlation", AmbientIgnition,
+             ["Part3 eqs(3.109)-(3.110),C4B/C5/C12; PDF pp546-547,647"],
+             limitations=["Printed C12=1.8 converts Kelvin flash point to Rankine; C4B=2.205 converts capped kg/s to lb/s inside natural log. Uses C5=25.2kg/s even though adjacent text also says25.22. Fit range is not certified by accepting source text: requires positive rate exponent and probability<=1, rejects nonphysical extrapolation. Zero rate uses limiting probability0. No LH2 ignition validation."])
+def ambient_ignition(x):
+    ensure(x.correlation_basis_source.strip(),"Ignition correlation basis must be nonblank")
+    t=1.8*x.flash_point_temperature_k
+    c=_COEFFICIENTS[x.phase]
+    a=((c[3]*t+c[2])*t+c[1])*t+c[0]
+    b=((c[7]*t+c[6])*t+c[5])*t+c[4]
+    ensure(math.isfinite(a) and math.isfinite(b) and b > 0,"Ignition polynomial outside finite increasing-rate domain")
+    rate=min(25.2,x.flammable_rate_kg_s)
+    if rate == 0:
+        probability=0
+    else:
+        logp=a+b*(math.log(2.205)+math.log(rate))
+        ensure(logp <= 0,"Ambient ignition correlation exceeds probability1; no implicit clamp")
+        probability=math.exp(logp)
+        ensure(probability > 0,"Ambient ignition probability underflow","numerical_failure")
+    return {"ambient_ignition_probability":probability,"evaluated_rate_kg_s":rate,"rate_capped":x.flammable_rate_kg_s > 25.2}
+
+
+class MaximumIgnition(Model):
+    phase: Phase
+    molecular_weight_kg_kmol: Positive | None = None
+
+
+@calculation("events/ait-ignition-maximum", "Ignition probability at autoignition temperature", MaximumIgnition,
+             ["Part3 eqs(3.111)-(3.112); PDF pp547-548"],
+             limitations=["Liquid maximum1; vapor follows printed molecular-weight interpolation with lower floor.7. MW2 gives.9 and MW>=170 gives.7. This is a release ignition probability correlation, not deterministic autoignition physics or a hydrogen certification."])
+def ignition_maximum(x):
+    if x.phase == "liquid":
+        return {"maximum_ignition_probability":1}
+    ensure(x.molecular_weight_kg_kmol is not None,"Vapor maximum requires molecular weight")
+    return {"maximum_ignition_probability":max(.7,.7+.2*(170-x.molecular_weight_kg_kmol)/168)}
+
+
+class TemperatureReference(Model):
+    reference_basis: Literal["api581_printed_c16","qualified_reference"]
+    qualified_reference_temperature_k: Positive | None = None
+    reference_source: str = Field(min_length=1,max_length=1000)
+
+
+def _reference(x):
+    ensure(x.reference_source.strip(),"Temperature reference source must be nonblank")
+    ensure((x.reference_basis == "qualified_reference") == (x.qualified_reference_temperature_k is not None),
+           "Provide temperature only for qualified_reference basis")
+    return 30.89 if x.reference_basis == "api581_printed_c16" else x.qualified_reference_temperature_k
+
+
+class TemperatureProbability(TemperatureReference):
+    probability_at_reference: Fraction
+    probability_at_ait: Fraction
+    storage_temperature_k: Positive
+    autoignition_temperature_k: Positive
+
+
+@calculation("events/probability-temperature", "Reference-to-AIT probability interpolation", TemperatureProbability,
+             ["Part3 eqs(3.113)-(3.114),(3.116)-(3.117),C16; PDF pp548,648"],
+             origin="api581_interpolation_with_explicit_temperature_reference",
+             limitations=["Printed C16=30.89K and US70R conflict with ambient-temperature discussion; no silent conversion to70F. Required basis chooses literal C16 or a caller-qualified reference. Ts below reference and AIT<=reference rejected; Ts>=AIT returns AIT limit. AIT probability must be at least reference probability. Qualified reference is not an authorized erratum."])
+def temperature_probability(x):
+    reference=_reference(x)
+    ensure(x.autoignition_temperature_k > reference and x.storage_temperature_k >= reference,
+           "Temperature interpolation requires AIT>reference and Ts>=reference")
+    ensure(x.probability_at_ait >= x.probability_at_reference,"AIT maximum is below reference probability")
+    weight=min(1,(x.storage_temperature_k-reference)/(x.autoignition_temperature_k-reference))
+    probability=x.probability_at_ait if weight == 1 else x.probability_at_reference+(x.probability_at_ait-x.probability_at_reference)*weight
+    return {"probability":probability,"temperature_weight":weight,"reference_temperature_k":reference}
+
+
+class ConditionalTable(Model):
+    phase: Phase
+    release_type: ReleaseType
+
+
+@calculation("events/conditional-table", "Table5.3 conditional event probabilities", ConditionalTable,
+             ["Part3 Table5.3; PDF p583"],limitations=["Historical expert-opinion defaults by phase/release type. Immediate probability at reference, not automatically at current storage temperature. VCE and flash columns are conditional on delayed ignition, not unconditional outcomes."])
+def conditional_table(x):
+    immediate,vce=_TABLE[x.release_type,x.phase]
+    return {"immediate_at_reference":immediate,"immediate_at_ait":1,"vce_given_delayed":vce,"flash_given_delayed":1-vce}
+
+
+class TwoPhaseProbability(Model):
+    liquid_probability: Fraction
+    vapor_probability: Fraction
+    flash_mass_fraction: Fraction
+    mixing_basis: Literal["api581_printed_liquid_weight","qualified_vapor_weight"]
+    weighting_source: str = Field(min_length=1,max_length=1000)
+
+
+@calculation("events/two-phase-probability", "Explicit-basis liquid/vapor conditional probability mixing", TwoPhaseProbability,
+             ["Part3 eqs(3.115),(3.118)-(3.119),(3.122); PDF pp548-549"],
+             origin="api581_with_explicit_two_phase_probability_weight",
+             limitations=["Printed equations weight LIQUID by flash fraction and VAPOR by its complement, despite flash-as-vapor wording. Literal and separately qualified vapor weighting are explicit alternatives; neither is silently selected or asserted as an erratum. Inputs must be probabilities of the same conditional quantity. No complete two-phase terminal tree or correlated max(jet,pool) consequence representation here."])
+def two_phase_probability(x):
+    ensure(x.weighting_source.strip(),"Two-phase probability weighting source must be nonblank")
+    liquid_weight=x.flash_mass_fraction if x.mixing_basis == "api581_printed_liquid_weight" else 1-x.flash_mass_fraction
+    return {"probability":liquid_weight*x.liquid_probability+(1-liquid_weight)*x.vapor_probability,"mixing_basis":x.mixing_basis}
+
+
+class SinglePhaseTree(Model):
+    phase: Phase
+    release_type: ReleaseType
+    case_type: Literal["leakage","rupture"]
+    ignition_probability: Fraction
+    immediate_given_ignition: Fraction
+    vce_given_delayed: Fraction
+    branch_probability_source: str = Field(min_length=1,max_length=1000)
+
+
+@calculation("events/single-phase-tree", "Single-phase terminal paths from Figures5.3/5.4", SinglePhaseTree,
+             ["Part3 Figures5.3/5.4; PDF pp586-587; conflicting example eq(3.125) PDF549"],
+             origin="api581_printed_figure_event_tree",
+             limitations=["Follows figure VCE path factor pvcedi, not conflicting printed example1-pvcedi. Figure5.3 labels instantaneous vapor leakage as fireball, although text restricts fireballs to rupture; output explicitly reports figure53_instantaneous applicability when needed. Rupture tree supported only for instantaneous release. Non-ignition is not automatically safe: toxic/NFNT must be evaluated separately."])
+def single_phase_tree(x):
+    ensure(x.branch_probability_source.strip(),"Branch probability source must be nonblank")
+    ensure(x.case_type != "rupture" or x.release_type == "instantaneous","Continuous rupture tree not specified by this implementation")
+    pi,immediate,vce=x.ignition_probability,x.immediate_given_ignition,x.vce_given_delayed
+    probabilities={event:0.0 for event in ["pool_fire","jet_fire","fireball","vce","flash_fire"]}
+    fire="pool_fire" if x.phase == "liquid" else "jet_fire" if x.release_type == "continuous" else "fireball"
+    probabilities[fire]=pi*immediate
+    probabilities["vce"]=pi*(1-immediate)*vce
+    probabilities["flash_fire"]=pi*(1-immediate)*(1-vce)
+    nonignition=1-pi
+    return {"terminal_probabilities":probabilities,"non_ignition_probability":nonignition,
+            "probability_sum":math.fsum([*probabilities.values(),nonignition]),
+            "fireball_applicability_basis":"figure53_instantaneous" if x.case_type == "leakage" and fire == "fireball" else "rupture_only_text"}
+
+
+class SinglePhaseIgnition(AmbientIgnition,TemperatureReference):
+    release_type: ReleaseType
+    case_type: Literal["leakage","rupture"]
+    molecular_weight_kg_kmol: Positive | None = None
+    storage_temperature_k: Positive
+    autoignition_temperature_k: Positive
+
+
+@calculation("events/single-phase-ignition", "Rate/flash-point/temperature to single-phase terminal probabilities", SinglePhaseIgnition,
+             ["Part3 eqs(3.109)-(3.117),Table5.3,Figures5.3/5.4; PDF pp547-549,583,586-587"],
+             origin="api581_single_phase_ignition_composition",
+             limitations=["Composes supplied phase flammable rate and qualified flash point, printed or qualified temperature reference, AIT maximum and Table5.3. Does not infer release phase, flash properties, two-phase mixture or probability validation from process data. Zero source has zero ignition even above AIT. Figure/text fireball conflict is reported via applicability basis. Continuous rupture and below-reference temperature unsupported."])
+def single_phase_ignition(x):
+    ambient=ambient_ignition(AmbientIgnition(**x.model_dump(include=set(AmbientIgnition.model_fields))))
+    maximum=ignition_maximum(MaximumIgnition(phase=x.phase,molecular_weight_kg_kmol=x.molecular_weight_kg_kmol))
+    reference=x.model_dump(include=set(TemperatureReference.model_fields))
+    setup=dict(**reference,storage_temperature_k=x.storage_temperature_k,autoignition_temperature_k=x.autoignition_temperature_k)
+    corrected=temperature_probability(TemperatureProbability(**setup,probability_at_reference=ambient["ambient_ignition_probability"],
+        probability_at_ait=maximum["maximum_ignition_probability"]))
+    table=conditional_table(ConditionalTable(phase=x.phase,release_type=x.release_type))
+    immediate=temperature_probability(TemperatureProbability(**setup,probability_at_reference=table["immediate_at_reference"],probability_at_ait=1))["probability"]
+    pi=corrected["probability"] if x.flammable_rate_kg_s > 0 else 0
+    tree=single_phase_tree(SinglePhaseTree(phase=x.phase,release_type=x.release_type,case_type=x.case_type,
+        ignition_probability=pi,immediate_given_ignition=immediate,vce_given_delayed=table["vce_given_delayed"],
+        branch_probability_source="Composed source correlations/Table5.3 with explicit reference basis"))
+    return {**ambient,**maximum,**tree,"ignition_probability":pi,"immediate_given_ignition":immediate,
+            "vce_given_delayed":table["vce_given_delayed"],"flash_given_delayed":table["flash_given_delayed"],
+            "reference_temperature_k":corrected["reference_temperature_k"],"temperature_weight":corrected["temperature_weight"]}