Cryogenic LNG Tank Boil-Off Gas (BOG) & Stratification Rollover Management Engineering Guide

Interactive Engineering Calculation Tool: Model tank heat leak, boil-off gas generation rate (BOR), dormant time to rollover, and relief valve venting dynamics using the companion LNG Tank Boil-Off Gas & Rollover Predictor.

1. Executive Summary & Cryogenic Storage Physics

Liquefied Natural Gas (LNG) is stored at atmospheric pressure and cryogenic temperatures ($\approx -161.5^\circ\text{C}$ to $-160^\circ\text{C}$ / $111.7\ \text{K}$ to $113.2\ \text{K}$) inside massive Full-Containment Storage Tanks ($160{,}000\ \text{m}^3$ to $260{,}000\ \text{m}^3$). Because the ambient environment is roughly $180\ \text{K}$ hotter than the bulk liquid, heat continuously penetrates the outer concrete containment wall, bottom perlite foam glass insulation, and suspended insulated deck.

This continuous thermal ingress induces two major operational phenomena: 1. Continuous Boil-Off Gas (BOG) Generation: Natural heat leak vaporizes the most volatile components (nitrogen and methane) at a design rate of $0.04\%$ to $0.08\%$ of total liquid volume per day ($\text{BOR \%/day}$). This BOG must be continuously evacuated, compressed, and either recondensed into export send-out headers or routed to a fuel gas network. 2. Double-Diffusive Stratification & Rollover: When a newly delivered LNG cargo with differing density and temperature is unloaded into a partially filled tank without sufficient mixing, two stable, non-mixing liquid layers form. Thermal convection traps superheated liquid in the bottom layer until the densities equalize, triggering an explosive rapid vaporization event known as LNG Tank Rollover. During rollover, BOG vaporization rates can surge by $10\times$ to $30\times$ the normal rate within hours, challenging tank pressure safety valves and risking atmospheric venting.

                  LNG TANK ROLLOVER MECHANISM (STRATIFICATION)
      STAGE 1: STRATIFIED STORAGE             STAGE 2: RAPID ROLLOVER EVENT
     +---------------------------+           +---------------------------+
     |        Vapor Space        |           |  MASSIVE BOG SURGE ^^^^   |
     +===========================+           +===========================+
     | Top Layer (Lighter LNG)   |           |                           |
     | Evaporates to Vapor Space |           |     Violent Convective    |
     | - T_top ~ -161.5°C        |  ====>    |     Overturn & Mixing     |
     | - Density ~ 430 kg/m3     |           |     of Superheated Bottom |
     |---------------------------|           |     Layer                 |
     | Bottom Layer (Trapped)    |           |                           |
     | Superheats from Base Heat |           |     Mass Transfer Spike   |
     | - T_bot ~ -158°C          |           |                           |
     | - Density ~ 435 -> 430    |           |                           |
     +---------------------------+           +---------------------------+

2. Thermodynamics of Heat Ingress & Boil-Off Rate (BOR)

2.1 Thermal Leakage Modeling

The total steady-state heat ingress into a full-containment LNG tank ($Q_{total}$, in $\text{kW}$ or $\text{Btu/hr}$) consists of four primary conductive and convective pathways:

$$Q_{total} = Q_{wall} + Q_{bottom} + Q_{roof} + Q_{penetrations}$$

Where: * Perlite Insulated Annulus Wall ($Q_{wall}$): $$Q_{wall} = \frac{2 \pi k_{perlite} H_{liquid} (T_{amb} - T_{LNG})}{\ln(R_{outer} / R_{inner})}$$ ($k_{perlite} \approx 0.035\text{--}0.045\ \text{W/(m}\cdot\text{K)}$) * Foam Glass Bottom Insulation ($Q_{bottom}$): $$Q_{bottom} = \frac{k_{foam} \pi R_{inner}^2 (T_{heating} - T_{LNG})}{t_{base}}$$ Where $T_{heating} \approx +5^\circ\text{C}$ maintained by electric or glycol base heating cables to prevent soil frost heaving. * Suspended Deck Insulated Roof ($Q_{roof}$): $$Q_{roof} = \frac{k_{deck} \pi R_{inner}^2 (T_{dome} - T_{vapor})}{t_{deck}}$$

2.2 Boil-Off Rate (BOR) Formulation

The normal daily boil-off rate ($\text{BOR}$, in percentage of total tank liquid inventory per 24 hours) is defined by NFPA 59A and BS EN 1473 as:

$$\text{BOR} = \frac{\dot{m}_{BOG} \times 24}{\rho_{LNG} \cdot V_{tank}} \times 100\% = \frac{Q_{total} \times 86{,}400}{\Delta H_{vap} \cdot \rho_{LNG} \cdot V_{tank}} \times 100\%$$

Where: * $\dot{m}{BOG}$ = Boil-off mass vaporization rate ($\text{kg/s}$) * $\Delta H}$ = Latent heat of vaporization of the LNG mixture ($\approx 510\ \text{kJ/kg}$ for pure $\text{CH4$; calculated via Peng-Robinson or GERG-2008 EOS) * $\rho^3$) * $V_{tank}$ = Net stored liquid volume ($\text{m}^3$)}$ = Saturated liquid density at tank operating pressure ($\approx 430\text{--}460\ \text{kg/m

3. Stratification & Rollover Physics (Double-Diffusive Convection)

3.1 The Turner Number and Convective Layer Trapping

When a lighter LNG cargo is added on top of a heavier cargo (or vice versa), the liquid separates into distinct cells due to double-diffusive convection. The stability of the density interface is governed by the Turner Stability Factor ($R_\rho$):

$$R_\rho = \frac{\beta \Delta S}{\alpha \Delta T}$$

Where $\alpha = -\frac{1}{\rho}\frac{\partial \rho}{\partial T}$ is the isobaric thermal expansion coefficient, $\beta = \frac{1}{\rho}\frac{\partial \rho}{\partial S}$ is the compositional contraction coefficient, and $\Delta S$ represents the difference in heavy hydrocarbon/nitrogen mole fractions.

Because heat ingress from the tank floor enters the bottom layer while vaporization at the surface only cools the top layer: 1. The bottom layer cannot vaporize; hydrostatic pressure ($P_{hydro} = \rho g H$) increases its saturation boiling point. 2. The bottom layer absorbs heat, increasing in temperature ($\Delta T \approx 1\text{--}4^\circ\text{C}$). 3. As the bottom layer warms, its thermal expansion decreases its density: $$\rho_{bot}(t) = \rho_{bot,0} \left[ 1 - \alpha (T_{bot}(t) - T_{bot,0}) \right]$$ 4. Simultaneously, the top layer loses volatile methane via surface evaporation, increasing its heavy hydrocarbon fraction and raising its density: $$\frac{d\rho_{top}}{dt} > 0$$

3.2 Time to Rollover ($t_{roll}$) Estimation

When $\rho_{bot}(t) = \rho_{top}(t)$, the stabilizing density barrier vanishes instantaneously. Gravitational buoyant inversion triggers explosive turnover. The dormant time to rollover can be approximated by Germain's thermodynamic model:

$$t_{roll} \approx \frac{\rho_{LNG} \cdot C_{p, LNG} \cdot H_{bot} \cdot (\rho_{bot,0} - \rho_{top,0})}{\alpha \cdot q_{bottom}''}$$

Where $q_{bottom}''$ is bottom heat flux ($\text{W/m}^2$) and $C_{p, LNG} \approx 3.5\ \text{kJ/(kg}\cdot\text{K)}$.

4. Boil-Off Gas (BOG) Handling System Architecture

                 CRYOGENIC BOG RECOVERY & RECONDENSATION TRAIN
+--------------------------------------------------------------------------+
| LNG Storage Tank (Full Containment)                                     |
+-------------------+------------------------------------------------------+
                    | (Low-Pressure BOG, -140°C, 1.05 - 1.15 bar(a))
                    v
          +-------------------+
          | BOG Suction Drum  |
          +---------+---------+
                    |
                    v
          +-------------------+
          | Low-Duty Cryo     |
          | BOG Compressor    | (Pressure: 1.1 bar -> 8-12 bar)
          +---------+---------+
                    |
                    +-----------------------+
                    | (To Recondenser)      | (To Fuel Gas / Gas Turbine)
                    v                       v
          +-------------------+   +-------------------+
          | LNG Sendout Stream|   | Plant Fuel Gas    |
          | Recondenser Drum  |   | System            |
          +---------+---------+   +-------------------+
                    |
                    v
          +-------------------+
          | HP Sendout Pumps  |
          | & Open Rack Vap.  |
          +---------+---------+
                    |
                    v
          +-------------------+
          | Sales Gas Grid    |
          +-------------------+

1. BOG Compressors (Low-Duty vs. High-Duty):
2. Low-Duty (LD) Compressors: Continuously compress normal holding BOG ($10\text{--}25\ \text{t/h}$) from storage tanks.
3. High-Duty (HD) Compressors: Sized for marine jetty unloading operations ($40\text{--}100\ \text{t/h}$) to handle vapor displaced by ship pump rates ($10{,}000\text{--}14{,}000\ \text{m}^3/\text{h}$) plus piston displacement.
4. Direct Contact Recondenser:
5. BOG at $8\text{--}12\ \text{bar(a)}$ is injected into a packed column where it contacts subcooled high-pressure LNG pumped from in-tank send-out pumps.
6. The subcooled LNG completely absorbs and re-condenses the BOG vapor without external refrigeration, eliminating flare losses.

5. Quantitative Technology Comparison Matrix

6. Step-by-Step BOG Sizing & Rollover Mitigation Workflow

   STEP 1: THERMAL INGRESS & BOR CALCULATION
   * Compile tank geometry (D_inner, H_liquid), insulation thicknesses, and k-values.
   * Calculate steady-state Q_total and verify design BOR <= 0.05% / day.
                         |
                         v
   STEP 2: CARGO DENSITY & COMPATIBILITY SCREENING
   * Measure incoming cargo density, temperature, and nitrogen/methane fractions.
   * Compare with current tank heel: delta_rho = |rho_cargo - rho_heel|.
                         |
                         v
   STEP 3: FILLING NOZZLE SELECTION
   * If incoming cargo is HEAVIER than tank heel: Select TOP FILL NOZZLE.
     (Heavy cargo cascades down through light heel, inducing natural turbulent mixing).
   * If incoming cargo is LIGHTER than tank heel: Select BOTTOM FILL NOZZLE.
     (Light cargo ascends through heavy heel, ensuring vigorous convective dispersion).
                         |
                         v
   STEP 4: CONTINUOUS LTD PROFILING
   * Operate motorized Level, Temperature, and Density (LTD) profilometer gauges.
   * Trigger control room alarms if delta_T > 1.0°C or delta_rho > 2.0 kg/m^3.
                         |
                         v
   STEP 5: ROLLOVER ACTIVE JET MIXING / TANK TRANSFER
   * If stratification persists > 48 hours: Activate bottom jet mixing eductors
     or recirculate liquid using in-tank pumps to transfer across tanks.

7. Fully Worked Industrial Example

Problem Statement:

An LNG receiving terminal operates a Full-Containment Storage Tank:

• Tank Inner Diameter ($D_i$): $80.0\ \text{m}$ ($R_i = 40.0\ \text{m}$)
• Maximum Liquid Height ($H_{liq}$): $36.0\ \text{m}$
• Gross Liquid Volume ($V_{tank}$): $180{,}955\ \text{m}^3$ (Operating fill: $160{,}000\ \text{m}^3$)
• Perlite Insulation Annulus Thickness: $1.0\ \text{m}$ ($k_{ins} = 0.038\ \text{W/(m}\cdot\text{K)}$)
• Foam Glass Base Thickness: $0.80\ \text{m}$ ($k_{base} = 0.045\ \text{W/(m}\cdot\text{K)}$)
• Suspended Deck Roof Area: $5{,}026.5\ \text{m}^2$ ($k_{deck} = 0.040\ \text{W/(m}\cdot\text{K)}$, thickness: $0.60\ \text{m}$)
• LNG Bulk Temperature: $-161.0^\circ\text{C}$ ($112.15\ \text{K}$)
• Ambient Design Temperature: $+35.0^\circ\text{C}$ ($308.15\ \text{K}$)
• Heated Base Slab Temperature: $+5.0^\circ\text{C}$ ($278.15\ \text{K}$)
• LNG Properties: Density $\rho = 445\ \text{kg/m}^3$; Latent Heat $\Delta H_{vap} = 510\ \text{kJ/kg}$

Step 1: Heat Ingress Calculations

1. Insulated Wall Heat Leak ($Q_{wall}$): $$R_{inner} = 40.0\ \text{m}, \quad R_{outer} = 41.0\ \text{m}$$ $$\ln(R_o / R_i) = \ln(41.0 / 40.0) = \ln(1.025) = 0.02469$$ $$\Delta T_{wall} = 35 - (-161) = 196\ \text{K}$$ $$Q_{wall} = \frac{2 \pi \times 0.038 \times 36.0 \times 196}{0.02469} = \frac{1{,}686.0}{0.02469} = 68{,}287\ \text{W} = \mathbf{68.29\ \text{kW}}$$

2. Insulated Bottom Heat Leak ($Q_{bottom}$): $$A_{base} = \pi \times (40.0)^2 = 5{,}026.5\ \text{m}^2$$ $$\Delta T_{base} = 5 - (-161) = 166\ \text{K}$$ $$Q_{bottom} = \frac{0.045 \times 5{,}026.5 \times 166}{0.80} = \frac{37{,}548}{0.80} = 46{,}935\ \text{W} = \mathbf{46.94\ \text{kW}}$$

3. Suspended Deck Roof Heat Leak ($Q_{roof}$): $$\Delta T_{roof} \approx 30 - (-150) = 180\ \text{K}$$ $$Q_{roof} = \frac{0.040 \times 5{,}026.5 \times 180}{0.60} = 60{,}318\ \text{W} = \mathbf{60.32\ \text{kW}}$$

4. Piping Penetrations & Structural Attachments ($\approx 10\%$ allowance): $$Q_{misc} = 0.10 \times (68.29 + 46.94 + 60.32) = \mathbf{17.55\ \text{kW}}$$

5. Total Tank Heat Ingress ($Q_{total}$): $$Q_{total} = 68.29 + 46.94 + 60.32 + 17.55 = \mathbf{193.10\ \text{kW}}$$

Step 2: Normal Boil-Off Mass Rate & BOR

• Daily Heat Ingress: $$E_{day} = 193.10\ \text{kW} \times 86{,}400\ \text{s} = 16{,}683{,}840\ \text{kJ/day}$$
• Mass of Boil-Off Gas: $$\dot{m}_{BOG, day} = \frac{16{,}683{,}840\ \text{kJ}}{510\ \text{kJ/kg}} = 32{,}713\ \text{kg/day} = \mathbf{32.71\ \text{tons/day}} \quad (\mathbf{1{,}363\ \text{kg/h}})$$
• Total Mass of Stored LNG: $$M_{LNG} = 160{,}000\ \text{m}^3 \times 445\ \text{kg/m}^3 = 71{,}200{,}000\ \text{kg}$$
• Boil-Off Rate ($\text{BOR}$): $$\text{BOR} = \frac{32{,}713}{71{,}200{,}000} \times 100\% = \mathbf{0.0459\%\ \text{per day}}$$
• Verification: The calculated BOR of $0.046\%/\text{day}$ comfortably complies with the international contractual guarantee threshold of $\le 0.050\%/\text{day}$.

Step 3: Low-Duty BOG Compressor Sizing

• Volumetric flow at compressor suction ($T = -140^\circ\text{C} = 133.15\ \text{K}$, $P = 1.05\ \text{bar(a)}$, $Z \approx 0.98$, $MW = 16.5\ \text{kg/kmol}$): $$\rho_{vapor} = \frac{105 \times 16.5}{8.314 \times 133.15 \times 0.98} = 1.597\ \text{kg/m}^3$$ $$Q_{BOG} = \frac{1{,}363\ \text{kg/h}}{1.597\ \text{kg/m}^3} = \mathbf{853.5\ \text{m}^3/\text{h}} \quad (\mathbf{502\ \text{ACFM}})$$

8. Diagnostic Troubleshooting & Failure Mode Matrix

Engineering Design Verification: To simulate thermal boil-off curves, model dormant stratification rollover times, and evaluate recondensation capacities, run the LNG Tank Boil-Off Gas & Rollover Predictor.