MaximaLabs
All unit operations

LNG storage tank simulation

1 inlet -> stored liquid + boil-off gas (reports weathering)

Governing equations

The exact equations the solver works for a lng storage tank — the same math shown in the app's "Theory" panel, not a black box.

n˙BOG=Q˙leakΔhvap,Ttank=Tbub(x,P)\dot n_{BOG} = \frac{\dot Q_{leak}}{\Delta h_{vap}},\qquad T_{tank} = T_{bub}(x, P)
yi=Ki(x,Tbub,P)xijKjxj(N2 and CH4 leave first)y_i = \frac{K_i(x, T_{bub}, P)\,x_i}{\sum_j K_j x_j}\quad(\text{N}_2\text{ and CH}_4\text{ leave first})
xistored=n˙inxin˙BOGyin˙inn˙BOG(weathering)x_i^{stored} = \frac{\dot n_{in} x_i - \dot n_{BOG} y_i}{\dot n_{in} - \dot n_{BOG}}\quad(\text{weathering})
n˙BOG\dot n_{BOG}
boil-off gas rate [mol/s]
Q˙leak\dot Q_{leak}
heat leaking into the tank [W]
Δhvap\Delta h_{vap}
enthalpy of the equilibrium vapour minus the liquid [J/mol]
TtankT_{tank}
tank temperature = bubble point at tank pressure [K]
TbubT_{bub}
bubble-point temperature [K]
PP
tank pressure [Pa]
xix_i
liquid mole fraction of component i [-]
yiy_i
boil-off vapour mole fraction of component i [-]
KiK_i
equilibrium K-value of component i [-]
n˙in\dot n_{in}
LNG fed to the tank [mol/s]
xistoredx_i^{stored}
mole fraction left in the stored liquid after weathering [-]

Parameters

storage tank whose boil-off rate comes from the heat leak ('heat_leak' [W], or a guaranteed 'boil_off_rate' [1/day]) and whose boil-off COMPOSITION comes from a real equilibrium flash -- N2 and methane leave first, so the BOG is light and the stored liquid weathers heavier. Outlets: stored liquid + BOG. Well-mixed and steady-state; the differential-vapour composition is capped (with a warning) once the boil-off cut gets large, and stratification is a different model -- see analysis/lng_rollover.py

Example flowsheets that use it

Stop fighting legacy software. Build your first flowsheet in 60 seconds.