Stack/system-level flowsheet modeling and techno-economic analysis (TEA) of green hydrogen produced by high-temperature solid oxide electrolysis cells (HT-SOEC) at industrial scale, with the produced hydrogen routed to Underground Hydrogen Storage (UHS). The study couples a process-simulator flowsheet with a Monte Carlo financial model to quantify the levelized cost of hydrogen and the economic risk of the system.
This repository is part of my PhD research. A companion repository, SOEC-Cell-First-Principles-Modeling, covers the cell-level, first-principles electrochemical model at laboratory scale.
- Process model & validation — an HT-SOEC stack/system model whose performance is validated against published experimental data (e.g. operation at 950 °C and Sorensen et al.).
- Techno-economic model — CAPEX/OPEX build-up and discounted cash-flow analysis producing the Net Present Value (NPV) and levelized cost of hydrogen.
- Monte Carlo simulation (uncertainty & risk) — propagates uncertainty in key techno-economic inputs (electricity price, capital cost, efficiency, etc.) using truncated normal distributions to produce probability distributions of NPV and identify the main cost drivers.
- Process flowsheet & simulation files — the underlying SOEC process simulation
(
.simx) and flowsheet (.drawio) used to size equipment and generate mass/energy balances.
TEA-UHS-Hydrogen/
├── notebooks/
│ ├── monte_carlo_v2.ipynb # Main Monte Carlo NPV / risk analysis
│ ├── MCS_HTSOEC.ipynb # Monte Carlo simulation of the HT-SOEC TEA
│ └── Model_validation.ipynb # SOEC performance model vs. literature data
├── data/
│ ├── HT-SOEC_TEA_optimized.xlsx # Optimized techno-economic model
│ ├── TEA_HTSOEC_H2_NEW.xlsx # Techno-economic input/results workbook
│ └── CashFlow_Modeling.xlsx # Discounted cash-flow model
├── simulation/
│ ├── *.simx # Process-simulator files (SOEC electrolysis rig & TEA)
│ └── SOEC_FS.drawio # Process flowsheet diagram
├── results/ # Generated figures (NPV distributions, CAPEX/OPEX, validation)
├── requirements.txt
├── LICENSE
└── README.md
Note on
.simxfiles: these are binaries for a proprietary process simulator and can only be opened in that software. They are versioned here for completeness; the file ending in…_finalreflects the most complete model.
Figures produced by the notebooks live in results/:
| File | Shows |
|---|---|
npv_probability_distribution_updated.png |
Monte Carlo probability distribution of project NPV |
output_capex_opex.png |
CAPEX / OPEX breakdown |
output_component_comparison.png |
Cost contribution by component |
Model_validation_950 C.png, Model_validation_soren_etal.png |
Model vs. experimental validation |
# (recommended) create a virtual environment
python -m venv .venv
source .venv/bin/activate # Windows: .venv\Scripts\activate
pip install -r requirements.txt
jupyter labOpen notebooks/monte_carlo_v2.ipynb to reproduce the techno-economic and risk analysis.
The notebooks read the workbooks in data/ — keep that folder alongside the
notebooks so the relative paths resolve.
| Area | Tools |
|---|---|
| Data & numerics | numpy, pandas, scipy |
| Financial modeling | numpy-financial (NPV / IRR) |
| Statistics & uncertainty | scipy.stats (truncated normal), Monte Carlo |
| Model validation | scikit-learn (R²), scipy.optimize (curve fitting) |
| Visualization | matplotlib, seaborn |
| Spreadsheets | openpyxl |
Krishnamurthy Busam — PhD research, techno-economic analysis of hydrogen systems.
Released under the MIT License. Excel workbooks, simulation files, and figures are shared for transparency and reuse within the terms of that license.