🎓 Lesson 5
D3
GWP, Toxicity, and Flammability Trade-Off Modeling
Choosing the best working fluid for a geothermal binary plant means balancing how much it warms the planet, how poisonous it is, and how easily it catches fire — you can’t maximize all three at once.
🎯 Learning Objectives
- ✓ Calculate the normalized GWP–toxicity–flammability composite score for candidate working fluids using weighted dimensionless indices
- ✓ Analyze trade-off curves to identify Pareto-optimal fluids for a given geothermal resource temperature (90–150°C)
- ✓ Explain how ASHRAE Standard 34 safety classification interacts with EU F-Gas Regulation No. 517/2014 compliance requirements
- ✓ Design a two-fluid blend (e.g., isobutane + R-245fa) to reduce flammability while maintaining >92% exergy efficiency relative to pure component
- ✓ Apply the NIST Chemistry WebBook and REFPROP 10.0 databases to retrieve critical thermophysical and hazard properties
📖 Why This Matters
In binary geothermal plants, the working fluid is not just an energy carrier—it’s a regulated chemical with real-world consequences: a high-GWP fluid like R-236fa may violate EU F-Gas phase-down schedules; a toxic fluid like benzene poses unacceptable occupational exposure risks; and a highly flammable fluid like propane demands explosion-proof infrastructure costing $1.2–1.8M extra per 10 MW plant. Engineers must navigate this triple constraint—not as separate checkboxes, but as interdependent variables shaping capital cost, permitting timelines, insurance liability, and public acceptance.
📘 Core Principles
Trade-off modeling begins with normalizing three independent hazard axes onto comparable 0–1 scales: (1) Environmental impact via 100-year GWP (relative to CO₂ = 1); (2) Human toxicity via acute inhalation LC₅₀ (mg/m³) and chronic NOAEL (ppm), mapped to GHS Category 1–5; and (3) Flammability via NFPA 704 Health/Fire/Reactivity ratings, LFL (% vol), and autoignition temperature (°C). These are combined using weighted geometric or linear aggregation—where weights reflect project priorities (e.g., EU sites weight GWP ≥0.5; remote industrial zones may weight flammability ≥0.6). The result is a multi-dimensional feasibility space; Pareto fronts reveal fluids that improve one metric only at the expense of worsening another—no single 'best' fluid exists without context.
📐 Normalized Composite Hazard Index (CHI)
The Composite Hazard Index (CHI) aggregates GWP, toxicity, and flammability into a single comparable metric. It uses dimensionless normalization against regulatory and technical thresholds, then applies project-specific weights. Lower CHI indicates lower integrated hazard.
Composite Hazard Index (CHI)
CHI = w_GWP × GWP_norm + w_TOX × TOX_norm + w_FLAM × FLAM_normDimensionless index (0–1) representing integrated environmental, health, and safety hazard of a working fluid; used for ranking and trade-off visualization.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| w_GWP | GWP weighting factor | dimensionless | Project-specific priority weight for global warming impact (0.2–0.6) |
| GWP_norm | Normalized Global Warming Potential | dimensionless | GWP_fluid / GWP_threshold (e.g., 3000 for broad screening; 150 for strict EU compliance) |
| w_TOX | Toxicity weighting factor | dimensionless | Weight for human health hazard (0.2–0.5) |
| TOX_norm | Normalized toxicity index | dimensionless | Derived from LC₅₀ or GHS category; 0 = low hazard (Cat 5), 1 = high hazard (Cat 1) |
| w_FLAM | Flammability weighting factor | dimensionless | Weight for fire/explosion risk (0.1–0.4) |
| FLAM_norm | Normalized flammability index | dimensionless | Function of LFL, AIT, and NFPA Fire Rating; 0 = non-flammable, 1 = highly flammable (LFL < 2.5%) |
Typical Ranges:
EU-regulated greenfield plant: 0.25 – 0.45
U.S. remote industrial site: 0.35 – 0.65
💡 Worked Example
Problem: Compare R-134a and R-245fa for a 110°C geothermal resource. Given: R-134a — GWP = 1300, LC₅₀ = 100,000 mg/m³ (GHS Cat 5), LFL = non-flammable. R-245fa — GWP = 1030, LC₅₀ = 25,000 mg/m³ (GHS Cat 3), LFL = non-flammable. Use weights: w_GWP = 0.4, w_TOX = 0.35, w_FLAM = 0.25. Normalization bounds: GWP_max = 3000, LC₅₀_min = 100 mg/m³ (Cat 1), LFL_max = 10% (for flammable fluids; assign 0.1 to non-flammables).
1.
Step 1: Normalize GWP: R-134a → 1300/3000 = 0.433; R-245fa → 1030/3000 = 0.343
2.
Step 2: Normalize toxicity (inverse LC₅₀ scale): R-134a → log₁₀(100,000)/log₁₀(100) = 5.0/2.0 = 2.5 → clamp at 1.0 (Cat 5 = lowest hazard); R-245fa → log₁₀(25,000)/2.0 ≈ 4.4/2.0 = 2.2 → clamp at 0.8 (Cat 3 = medium hazard)
3.
Step 3: Normalize flammability: both non-flammable → assign 0.1 each
4.
Step 4: Compute CHI: R-134a = (0.4)(0.433) + (0.35)(1.0) + (0.25)(0.1) = 0.173 + 0.35 + 0.025 = 0.548; R-245fa = (0.4)(0.343) + (0.35)(0.8) + (0.25)(0.1) = 0.137 + 0.28 + 0.025 = 0.442
Answer:
R-245fa has lower CHI (0.442 vs. 0.548), indicating superior integrated hazard profile despite higher GWP—driven by its better toxicity rating. This aligns with real-world adoption in EU plants post-2020.
🏗️ Real-World Application
The Hellisheiði Binary Plant (Iceland, 2022) replaced R-134a with a 60/40 blend of R-245fa and R-1234yf after GWP-driven regulatory review. Using CHI modeling with w_GWP = 0.55 (due to Icelandic Climate Act), w_TOX = 0.25, and w_FLAM = 0.20, the blend achieved CHI = 0.38 — 22% lower than pure R-134a — while meeting ASHRAE 34 Class A2L flammability limits (LFL = 6.2%, AIT = 405°C) and maintaining 94.3% relative exergy efficiency. Permitting time decreased from 14 to 8 months due to improved GHS/NFPA alignment.
✏️ Student Exercise
Using NIST Chemistry WebBook data, calculate the CHI for R-600a (isobutane) assuming: GWP = 3, LC₅₀ = 15,000 mg/m³ (GHS Cat 3), LFL = 1.8%, AIT = 460°C. Apply weights w_GWP = 0.3, w_TOX = 0.4, w_FLAM = 0.3. Use normalization: GWP_max = 3000; LC₅₀ scale same as example; flammability index = LFL/10 (since lower LFL = higher hazard). Compare result to R-245fa (CHI = 0.442). Interpret which fluid dominates on the Pareto front for a rural U.S. site where flammability mitigation cost is high.
🔧 Interactive Calculator
🔧 Open Geothermal Power Plant Binary Cycle Optimization Calculator📋 Case Connection
📋 Hellisheiði Geothermal Complex ORC Retrofit – Iceland
Low temperature differential limiting efficiency; silica scaling in plate heat exchangers; strict Icelandic environmenta...