PyThermoCalcDB — Python Thermodynamic Calculation Engine powered by PyThermoDB
PyThermoCalcDB is a companion package to PyThermoDB, designed to provide a robust set of thermodynamic calculations, utilities, and workflows — using thermodynamic properties, constants, and equations built with PyThermoDB.
- Seamless integration with PyThermoDB: use the database of thermodynamic data (pure component properties, correlations, etc.) as the foundation.
- Reusable calculation tools: functions for various calculations.
- Consistent with your thermodynamic ecosystem: fits with existing packages like PyThermoModels, PyPhaseEQ, PyThermoFlash — offering a central “calculation engine” layer.
- Minimal dependencies: lean design, using PyThermoDB as core data source; optionally depend on well-known numerical libraries (e.g. NumPy, SciPy) if needed for complex calculations.
- Flexible use-cases: from quick property lookups to advanced process calculations, or integration in larger workflows or multi-agent systems.
pip install pyThermoCalcDBPyThermoCalcDB provides high-level thermodynamic routines built on PyThermoDB. It is optimized for scripts, notebooks, and pipelines ? the examples/ directory contains runnable workflows.
- Collect ThermoDB pickle paths for your components and load them with
pyThermoLinkDB.load_and_build_model_source.
import os
import pyThermoLinkDB as ptdblink
from pythermodb_settings.models import Component, ComponentRule, ComponentThermoDBSource, Temperature, Pressure
CO2 = Component(name="carbon dioxide", formula="CO2", state="g")
thermodb_dir = "/path/to/thermodb"
CO2_src = ComponentThermoDBSource(
component=CO2,
source=os.path.join(thermodb_dir, "carbon dioxide-g.pkl"),
)
model_source = ptdblink.load_and_build_model_source(
thermodb_sources=[CO2_src],
original_equation_label=False
)- Thermodynamic properties (
pyThermoCalcDB.docs.thermo):calc_En,calc_GiFrEn,calc_En_range,calc_GiFrEn_range,calc_dEn(Delta H between temperatures),calc_dEnt(Delta S with pressures and phase),calc_En_mix(mixture enthalpy with optional departure/excess terms andoutput_unit). - Saturation & phase change (
pyThermoCalcDB.docs.sat):calc_VaPr,calc_VaPr_range,calc_EnVap,calc_EnVap_range,calc_T_sat,calc_VaPr_sensitivity,calc_VaPr_sensitivity_range. - Reactions (
pyThermoCalcDB.reactions.reactions):dEn_rxn_STDanddGiFrEn_rxn_STDfor standard reaction energetics.
Example:
from pyThermoCalcDB.docs import thermo, sat
from pythermodb_settings.models import Temperature, Pressure
T = Temperature(value=300.0, unit="K")
P = Pressure(value=1.5e6, unit="Pa")
enthalpy = thermo.calc_En(component=CO2, model_source=model_source, temperature=T)
gibbs = thermo.calc_GiFrEn(component=CO2, model_source=model_source, temperature=T, phase="IG")
VaPr = sat.calc_VaPr(component=CO2, model_source=model_source, temperature=T)
Tsat = sat.calc_T_sat(
component=CO2,
model_source=model_source,
pressure=P,
temperature_guess=T,
method="least_squares",
)Notes:
- All helpers expect typed
Component,Temperature, andPressureinputs plus amodel_sourcefrom PyThermoLinkDB. - The optional
modekwarg acceptssilent,log, orattachto control logging and timing.
Contributions are highly welcome — bug fixes, new calculation routines, mixture models, extended unit tests, documentation, etc.
This project is distributed under the Apache License, Version 2.0, which grants you broad freedom to use, modify, and integrate the software into your own applications or projects, provided that you comply with the conditions outlined in the license. Although Apache 2.0 does not require users to retain explicit author credit beyond standard copyright and license notices, I kindly request that if you incorporate this work into your own software, you acknowledge Sina Gilassi as the original author. Referencing the original repository or documentation is appreciated, as it helps recognize the effort invested in developing and maintaining this project.
For any question, contact me on LinkedIn