otex 0.3.0

OTEX - Ocean Thermal Energy eXchange: OTEC plant design, simulation, and analysis

OTEX - Ocean Thermal Energy eXchange

A Python library for OTEC plant design, simulation, and techno-economic analysis

Features • Installation • Quick Start • Documentation • Citation

---

Overview

OTEX (Ocean Thermal Energy eXchange) is a Python library for designing, simulating, and analyzing Ocean Thermal Energy Conversion (OTEC) power plants. It integrates with global oceanographic databases to enable site-specific techno-economic assessments anywhere in the tropical oceans.

OTEX enables researchers and engineers to:

• Design OTEC plants with multiple thermodynamic cycles and working fluids
• Analyze regional and global potential using CMEMS or HYCOM oceanographic data
• Perform uncertainty analysis with Monte Carlo simulations and sensitivity studies
• Compare scenarios across different locations, plant sizes, and configurations

Features

Thermodynamic Cycles

Cycle	Description	Status
Rankine Closed	Ammonia/organic fluid closed loop	✅ Stable
Rankine Open	Flash evaporation of seawater	✅ Stable
Rankine Hybrid	Combined closed/open cycle	✅ Stable
Kalina	Ammonia-water mixture	✅ Stable
Uehara	Advanced ammonia-water cycle	✅ Stable

Working Fluids

• Ammonia (NH₃) - Default, polynomial or CoolProp
• R134a - Requires CoolProp
• R245fa - Requires CoolProp
• Propane - Requires CoolProp
• Isobutane - Requires CoolProp

Data Sources

Source	Resolution	Depth Levels	Period	Authentication
CMEMS	0.083°	50	1993–present	Required (free account)
HYCOM	0.08°	40	1994–2015, 2019–2024	Not required

Analysis Capabilities

• Regional Analysis: Site-specific LCOE maps and power profiles
• Multi-year Simulations (0.2.0): Run a continuous N-year simulation with NPV-based LCOE, configurable degradation (constant / logistic / step), and OPEX escalation (flat / fixed-rate / indexed). Inter-annual variability of AEP is reported per site.
• Site Screening: Optional exclusion of protected areas (WDPA) and busy shipping lanes (World Bank vessel density), plus seismic and cyclone risk-based cost multipliers (GEM PGA, NOAA IBTrACS)
• Uncertainty Analysis: Monte Carlo with Latin Hypercube Sampling
• Sensitivity Analysis: Sobol indices and Tornado diagrams
• Off-design Performance: Time-resolved power output profiles

Installation

Basic Installation

pip install otex

With Optional Dependencies

# High-accuracy fluid properties
pip install otex[coolprop]

# Uncertainty analysis (Sobol indices)
pip install otex[uncertainty]

# Siting layers (protected areas, shipping lanes, seismic and cyclone hazards)
pip install otex[siting]

# All optional dependencies
pip install otex[all]

Development Installation

git clone https://github.com/msotocalvo/OTEX.git
cd OTEX
pip install -e ".[dev]"

Oceanographic Data Access

HYCOM (no credentials needed):

from otex.regional import run_regional_analysis

otec_plants, sites = run_regional_analysis(
    studied_region='Jamaica',
    data_source='HYCOM',
    year_start=2020,
    year_end=2020,
)

CMEMS (requires free account):

1. Create account at Copernicus Marine
2. Configure credentials:

   copernicusmarine login
   

See Installation Guide for detailed instructions.

Quick Start

Basic Plant Configuration

from otex.config import parameters_and_constants

# Configure a 100 MW OTEC plant
inputs = parameters_and_constants(
    p_gross=-100000,           # 100 MW (negative = power output)
    cost_level='low_cost',
    cycle_type='rankine_closed',
    fluid_type='ammonia',
    year_start=2020,
    year_end=2020,
)

print(f"Cycle: {inputs['cycle_type']}")
print(f"Discount rate: {inputs['discount_rate']:.1%}")
print(f"Plant lifetime: {inputs['lifetime']} years")

Regional Analysis

# Analyze Cuba for 2020 with a 50 MW plant (CMEMS, default)
otex-regional Cuba --year 2020 --power -50000

# Using HYCOM data (no credentials needed)
otex-regional Philippines --year 2020 --data-source HYCOM

# Analyze with Kalina cycle
otex-regional Philippines --cycle kalina --year 2021

Multi-Year Analysis

Since 0.2.0, a single run can span multiple calendar years. The oceanographic data is concatenated along the time axis, LCOE is computed via discounted cashflow NPV (with configurable degradation and OPEX escalation), and the output CSV includes inter-annual AEP statistics:

# Continuous 4-year simulation, NPV LCOE, inter-annual AEP variability
otex-regional Jamaica --year-start 2020 --year-end 2023

from otex.regional import run_regional_analysis

otec_plants, sites = run_regional_analysis(
    studied_region='Jamaica',
    year_start=2020,
    year_end=2023,
)
# sites includes columns: LCOE (NPV), LCOE_legacy (CRF), AEP_min/p50/max/std
# A second CSV `OTEC_sites_yearly_*.csv` reports per-(site, year) energy.

Uncertainty Analysis

from otex.analysis import (
    MonteCarloAnalysis,
    UncertaintyConfig,
    TornadoAnalysis,
    plot_histogram,
    plot_tornado
)

# Monte Carlo analysis
config = UncertaintyConfig(n_samples=1000, seed=42)
mc = MonteCarloAnalysis(T_WW=28.0, T_CW=5.0, config=config)
results = mc.run()

# Get statistics
stats = results.compute_statistics()
print(f"LCOE: {stats['lcoe']['lcoe_mean']:.2f} ± {stats['lcoe']['lcoe_std']:.2f} ct/kWh")
print(f"90% CI: [{stats['lcoe']['lcoe_p5']:.2f}, {stats['lcoe']['lcoe_p95']:.2f}]")

# Tornado diagram
tornado = TornadoAnalysis(T_WW=28.0, T_CW=5.0)
tornado_results = tornado.run()
plot_tornado(tornado_results)

Command Line Interface

# Tornado analysis
python scripts/uncertainty_analysis.py --T_WW 28 --T_CW 5 --method tornado

# Monte Carlo with 500 samples
python scripts/uncertainty_analysis.py --T_WW 28 --T_CW 5 --method monte-carlo --samples 500

# Full analysis with plots
python scripts/uncertainty_analysis.py --T_WW 28 --T_CW 5 --method all --samples 200 --save-plots

Documentation

Document	Description
Installation Guide	Detailed setup instructions
Quick Start Tutorial	Get started in 10 minutes
Regional Analysis	Analyze specific regions
Uncertainty Analysis	Monte Carlo and sensitivity
API Reference	Complete API documentation
01 - Quick Start	Basic plant sizing and cost analysis
02 - Regional Analysis	Analyze OTEC potential for a region
03 - Uncertainty Analysis	Monte Carlo, Tornado, Sobol

Project Structure

OTEX/
├── otex/                    # Main package
│   ├── core/               # Thermodynamic cycles and fluids
│   ├── plant/              # Plant sizing and operation
│   ├── economics/          # Cost models and LCOE
│   ├── analysis/           # Uncertainty and sensitivity
│   ├── data/               # Data loading (CMEMS, HYCOM, NetCDF)
│   └── config.py           # Configuration management
├── scripts/                 # CLI scripts
│   ├── regional_analysis.py
│   ├── global_analysis.py
│   └── uncertainty_analysis.py
├── tests/                   # Test suite
├── docs/                    # Documentation
└── data/                    # Reference data files

Configuration Options

Parameter	Options	Default
cycle_type	rankine_closed, rankine_open, rankine_hybrid, kalina, uehara	rankine_closed
fluid_type	ammonia, r134a, r245fa, propane, isobutane	ammonia
cost_level	'low_cost', 'high_cost', or a CostScheme object	'low_cost'
p_gross	Any negative value (kW)	-136000
data_source	'CMEMS', 'HYCOM'	'CMEMS'
year	1993–present (CMEMS), 1994–2015 / 2019–2024 (HYCOM)	2020

Custom Cost Schemes

Beyond the two built-in scenarios you can define your own cost parameters with CostScheme and Python's standard dataclasses.replace():

from otex.economics import CostScheme, LOW_COST
from dataclasses import replace

# Modify specific parameters of an existing scheme
my_scheme = replace(LOW_COST, turbine_coeff=400, opex_fraction=0.04)

# Use it everywhere cost_level is accepted
inputs = parameters_and_constants(p_gross=-100000, cost_level=my_scheme)
costs, capex, opex, lcoe = capex_opex_lcoe(plant, inputs, my_scheme)

All existing code that uses cost_level='low_cost' or cost_level='high_cost' continues to work unchanged.

Requirements

• Python >= 3.9
• NumPy, Pandas, SciPy, Matplotlib
• xarray, netCDF4 (oceanographic data)
• tqdm (progress bars)

Optional:

• CoolProp (additional working fluids)
• SALib (Sobol sensitivity analysis)

Acknowledgments

OTEX builds upon pyOTEC by Langer et al. For the original methodology, see:

Langer, J., Blok, K. The global techno-economic potential of floating, closed-cycle ocean thermal energy conversion. J. Ocean Eng. Mar. Energy (2023). https://doi.org/10.1007/s40722-023-00301-1

Citation

If you use OTEX in your research, please cite:

• Soto Calvo M, and Lee HS., 2025. Ocean Thermal Energy Conversion (OTEC) Potential in Central American and Caribbean Regions: A Multicriteria Analysis for Optimal Sites. Applied Energy. 394: 126182. https://doi.org/10.1016/j.apenergy.2025.126182

Contributing

We welcome contributions! Please see CONTRIBUTING.md for guidelines.

License

This project is licensed under the MIT License - see the LICENSE file for details.

---