Mass Conservation Equations for Kraft Digester Organic Black Liquor Balance

Cook

Industrial pulping converts wood chips into cellulose fiber through high-temperature alkaline delignification. Closing an organic mass balance across this extraction requires tracking every material stream entering the continuous digester control volume. The solid feed arrives as wood biomass containing cellulose, hemicellulose, lignin, and extractives, carrying native moisture that dilutes aqueous reagents introduced into the vessel.

White liquor serves as the primary liquid reactant, supplying sodium hydroxide and sodium sulfide in active alkali charges. Weak wash from washer filtrate tanks and direct live steam injection make up the remaining mass flows into the reaction zone.

Establishing the overall mass boundary requires measuring chip moisture content and wood density accurately. Moisture varies by season, species blend, and storage pile retention time. Missing the true chip dry-weight baseline causes systematic errors in every subsequent calculation of dissolved organic yield.

When wet chips pass through the chip bin or atmospheric pre-steamer, direct steam condensation adds liquid straight into the chip pores. Tracking organic mass flow rates across the extraction screens quantifies digester yield. This condensation alters the liquor-to-wood ratio before chips enter the impregnation phase.

Inputs to the Digester Control Volume

Total liquid mass entering the digester includes bound moisture within the wood pores, steam condensate from pre-steaming vessels, primary white liquor, and utility liquor charges. The mass conservation equation for the liquid phase entering the top of the vessel is:

M_liquid_in = M_wood_green (MC / 100) + M_steam_cond + M_WL + M_rec

Where M_wood_green is the green mass flow rate of wood chips, MC is the wet-basis moisture content as a percentage, M_steam_cond is the mass flow rate of pre-steaming condensate, M_WL is the white liquor mass flow rate, and M_rec is the mass flow rate of recycled black or weak wash liquor. The dry biomass mass flow rate entering the vessel is defined by:

M_wood_OD = M_wood_green (1 – MC / 100)

Chemical charges depend directly on M_wood_OD. The white liquor mass flow rate is set by target Effective Alkali and Active Alkali charges. Effective Alkali ~ calculated as sodium hydroxide plus half the sodium sulfide ~ supplies the hydroxyl ions needed to cleave lignin ether bonds and dissolve carbohydrate fragments.

White Liquor Stoichiometric Ratios

White liquor density shifts predictably with temperature and total titratable alkali concentration. Expressing chemical charges as equivalent sodium oxide mass simplifies converting volumetric dosing rates into mass flow rates. Active alkali concentration typically runs between 90 and 120 grams as sodium oxide per liter of liquor.

Sulfidity, the ratio of sodium sulfide to active alkali, drives the kinetics of lignin fragmentation while limiting carbohydrate chain degradation.

Wood moisture directly dictates the allowable white liquor charge. Higher moisture leaves less room for liquid volume at a fixed liquor-to-wood target, forcing reliance on higher chemical concentrations in the white liquor stream. If the target liquid-to-wood mass ratio is 3.5 to 1, total liquid mass inside the digester must equal 3.5 times M_wood_OD.

The white liquor mass flow rate required to deliver an Effective Alkali charge of EA_percent on dry wood is given by:

M_WL = (M_wood_OD (EA_percent / 100)) / (C_EA_mass)

Where C_EA_mass is the mass fraction of effective alkali in white liquor. Balancing this mass ensures that incoming hydroxyl and hydrosulfide ion concentrations meet the chemical demand to break down native lignin without unnecessarily degrading cellulose chains.

In continuous digestion, holding a precise liquor-to-wood ratio prevents channeling and keeps chemical distribution uniform through the chip column. If chip moisture spikes unexpectedly, white liquor volumetric flows or alkali concentrations must be adjusted immediately to prevent localized liquor depletion. Conversely, excess unreacted alkali leaving the impregnation zone triggers carbohydrate peeling, cutting pulp yield while increasing the dissolved organic load in the weak black liquor.

Degradation

Alkaline pulping breaks ether linkages in the macromolecular lignin matrix while attacking wood polysaccharides. The organic mass dissolved into the liquor falls into four main fractions: fragmented kraft lignin, degraded hemicellulose, solubilized cellulose monomer units, and saponified extractives. Lignin breaks down mostly through the rupture of beta-O-4 aryl ether bonds, turning insoluble protolignin into alkali-soluble fragments bearing phenolate and carboxylate groups.

These polar functional groups keep dissolved lignin soluble in the strongly alkaline solution.

Carbohydrate solubilization proceeds through two primary mechanisms: end-wise alkaline peeling and alkaline hydrolysis. Peeling reactions strip sugar units from the reducing ends of polysaccharide chains, forming aliphatic carboxylic acids ~ predominantly glucoisosaccharinic, formic, acetic, and lactic acids. These organic acids rapidly neutralize hydroxyl ions in the liquor, consuming effective alkali and binding sodium as organic sodium salts.

Cellulose resists alkaline degradation much better than hemicellulose, but high temperatures still cause chain cleavage, adding to the total organic carbon load in the liquid phase.

Dissolved organic solids in kraft black liquor carry over sixty percent of the chemical energy entering the recovery boiler.

Lignin Solubilization Kinetics and Mass Yield

The rate of lignin dissolution depends on temperature, hydroxyl and hydrosulfide ion concentrations, and the amount of unreacted lignin remaining in the wood. Kinetic models divide delignification into initial, bulk, and residual phases. Bulk delignification follows a pseudo-first-order rate equation relative to remaining lignin mass:

dM_lignin / dt = – k_bulk M_lignin ^a ^b

Where M_lignin is the unreacted lignin mass remaining in the wood fiber, k_bulk is the temperature-dependent Arrhenius rate constant, is the hydroxyl ion concentration, is the hydrosulfide ion concentration, and exponent coefficients a and b reflect empirical reaction orders. Integrating this kinetic rate over total retention time gives the mass of lignin entering the black liquor:

M_lignin_dissolved = M_lignin_initial – M_lignin_pulp

Where M_lignin_initial comes from the Klason and acid-soluble lignin content of the raw wood, and M_lignin_pulp is derived from the unbleached pulp Kappa number. One Kappa unit represents roughly 0.147 percent residual lignin on oven-dry pulp.

Polysaccharide Dissolution and Hydroxy Acid Formation

Hemicelluloses ~ xylan in hardwoods and galactoglucomannan in softwoods ~ dissolve rapidly early in the cook. Part of this dissolved hemicellulose stays suspended in the black liquor as polymeric carbohydrate chains, while the rest breaks down into monomeric hydroxy acids. As alkali neutralizes these carboxylic acids, the reaction dictates how sodium partitions between active inorganic chemicals and organic sodium complexes.

Carbohydrate mass balances across the digester quantify the dissolved organic fraction coming from carbohydrate degradation. The total mass of organic dissolved solids generated per metric ton of oven-dry wood is defined as:

M_org_solids = M_wood_OD (1 – Y_pulp) + M_ext_sol

Where Y_pulp is total unbleached pulp mass yield as a decimal fraction, and M_ext_sol is the mass of dissolved extractives, including resin acids, fatty acids, and unsaponifiable material. Saponification turns fatty and resin acids into soluble sodium soaps, later skimmed off as tall oil soap in downstream recovery.

Tracking sodium consumption across these dissolved organic fractions is essential to close the coupled inorganic-organic mass balance. Every gram of degraded carbohydrate produces carboxylic acid groups that consume stoichiometric amounts of sodium hydroxide. Underestimating carbohydrate yield loss leads directly to under-dosing white liquor, resulting in raw knots, high screen rejects, and erratic pulp quality.

Vessel

Continuous hydraulic digesters separate pulping into distinct thermodynamic and kinetic zones. As chips travel down the vessel driven by gravity and hydraulic forces, they contact liquor flowing co-currently or counter-currently. Mass conservation equations are set up separately for each zone ~ impregnation, heating, cooking, extraction, and washing ~ with boundary conditions capturing screen extraction around the shell and liquor injection through central pipes.

Extraction screens serve as the main mass splitter in a continuous digester. Spent liquor carrying high concentrations of dissolved organic solids and residual alkali is drawn through the screen plates toward the evaporators. Meanwhile, wash liquor injected at the vessel bottom flows upward counter-currently through the pulp column, displacing dissolved organics from the bed.

This wash zone functions as a differential mass transfer column where organic solids diffuse out of porous fiber walls into the bulk liquor.

At an effective alkali charge of eighteen percent active alkali on oven-dry wood at one hundred and sixty-five degrees Celsius, dissolved lignin accounts for seventy-two percent of total dissolved organic mass.

Zonal Mass Balances across Impregnation and Digestion

In the impregnation zone, chips absorb liquor and heat at temperatures below active delignification. The mass balance matches incoming raw chips and liquor against the impregnated chips and free liquor moving into the cooking zone. Spatial conservation of organic species along digester height z is modeled with the steady-state convection-dispersion-reaction equation:

u_L (dC_i / dz) = D_axial (d^2 C_i / dz^2) + R_i

Where u_L is superficial liquid velocity through the bed, C_i is the concentration of dissolved organic component i (such as dissolved lignin or specific organic acid salts), D_axial is the axial dispersion coefficient for back-mixing, and R_i is the volumetric reaction rate of component i. In co-current zones, liquid and solids move together downward, simplifying convective transport terms.

At the extraction zone, a large fraction of free liquor is withdrawn. The mass conservation equation around the extraction screen envelope is:

M_liquor_in_zone = M_extracted_BL + M_liquor_down

Where M_liquor_in_zone is total liquid mass flow entering the extraction level, M_extracted_BL is extracted weak black liquor sent to the evaporators, and M_liquor_down is residual liquor flowing into the wash zone. Hydraulic stability requires balancing M_extracted_BL against wash liquor upflow to prevent chip column hanging or severe liquor channeling.

Which Extraction Screen Configuration Minimizes Organic Carryover?

Designing extraction screens with stepped slots and controlled face velocities prevents chip mat compaction while keeping solids extraction steady. Splitting extraction across dual screen sets increases liquor removal without exceeding the fluid drag threshold that pulls fiber into recirculation pumps. Proper screen operation balances black liquor leaving the system against wash liquor coming up from the diffusion zone.

When extraction screens plug or scale up, dissolved organics bypass extraction and carry down into the wash zone and unbleached pulp line. This carryover increases the chemical oxygen demand of pulp slurry feeding the bleach plant or paper machine. Unmeasured hemicellulose degradation introduces up to a 3.2 percent error into dissolved organic calculations.

• Screen Blindness Compaction occurs when local screen face velocity exceeds critical hydrodynamics, forcing flexible wood chips against screen slots and stopping black liquor extraction.
• Channeling Short Circuiting develops when non-uniform chip packing allows wash liquor to bypass the central core of the chip column, reducing organic displacement efficiency.
• Temperature Extraction Loss arises when recirculation heat exchangers foul, dropping zonal temperatures and slowing the diffusion rate of macromolecular lignin out of fiber cell walls.
• Dilution Factor Inversion happens when excessive wash water injection forces high volumes of weak liquor back up into the cooking zone, cooling the bulk delignification reaction premature to completion.

High-capacity extraction screens are engineered to maintain uniform liquor draw regardless of chip size distribution, but in practice, chip thickness variations alter bed permeability and hydraulic resistance. Oversized chips slow lignin diffusion out of the core fiber structure, trapping dissolved organics inside the chip matrix until past the extraction zone. That organic material then bleeds into the wash zone, driving up washing losses and diluting black liquor routed to the evaporators.

Concentration

Determining total dissolved solids in weak black liquor requires accounting for organic carbon compounds alongside residual inorganic salts. Organic black liquor solids consist of carbon, hydrogen, oxygen, nitrogen, and sulfur species complexed with sodium cations. Because black liquor density changes with solids concentration and temperature, closing recovery mass balances requires converting volumetric flow meter readings to mass flow rates using non-linear fluid density correlations.

Total Black Liquor Solids content is defined as the dry mass percentage remaining after water is removed under standard drying conditions. The total dissolved solids mass flow rate is calculated as:

M_TBLS = M_BL_total (w_solids / 100)

Where M_BL_total is the total black liquor mass flow rate and w_solids is the dry solids weight fraction. Total solids split into an organic mass fraction w_org and an inorganic mass fraction w_inorg. The organic-to-inorganic mass ratio (O/I ratio) dictates the heating value and combustion behavior of concentrated black liquor in the recovery boiler.

Gravimetric and Refractometric Solids Measurement

Standard laboratory measurement of black liquor solids relies on gravimetric oven drying at 105 degrees Celsius. Volatiles ~ including low-boiling alcohols, organic sulfur compounds, and light fatty acids ~ evaporate with water during drying, slightly understating total organic mass. Online process control relies instead on optical refractometers that measure refractive index and convert the reading to total dissolved solids using empirical calibration curves.

Refractometric readings respond primarily to dissolved organics and sodium salts. Refractometer calibrations must be updated whenever wood species or effective alkali charges shift significantly, since changes in the ratio of dissolved lignin to inorganic salts alter the solution’s optical refractive index.

Organic to Inorganic Mass Ratios in Weak Black Liquor

The organic fraction w_org carries the thermal energy of the liquor, while the inorganic fraction w_inorg supplies the chemical charge needed for smelting and chemical recovery. The balance between these fractions is determined as follows:

1. Sample weak black liquor at the digester extraction line, capturing liquor in a sealed, zero-headspace pressure vessel to prevent flash evaporation of light organics.
2. Measure total dissolved solids content according to TAPPI T 650 by drying a known mass sample on inert sand or filter paper matrices at 105 degrees Celsius for sixteen hours.
3. Perform high-temperature muffle furnace ashing at 525 degrees Celsius to burn off all organic carbon, leaving residual inorganic ash as sodium carbonate, sodium sulfate, and residual inorganic salts.
4. Calculate the organic mass fraction by subtracting the inorganic ash weight from the initial total dry solids weight, correcting for the thermal oxidation of sodium sulfide to sodium sulfate during high-temperature ashing.
5. Determine elemental carbon, hydrogen, nitrogen, and sulfur contents on dry solids using instrumental CHNS elemental analyzers to close the carbon-specific organic mass balance.

Standard test method TAPPI T 650 specifies gravimetric oven drying at one hundred and five degrees Celsius for sixteen hours to determine total black liquor solids content.

When writing liquor supply contracts or comparing mill performance, purchase specifications usually require test reports to state whether total solids reflect as-received refractometer values or gravimetric oven-dry masses corrected for volatile loss. Specifying TAPPI T 650 compliance obligates both parties to base yield claims and thermal guarantees on gravimetric dry weight, avoiding disputes over refractometer calibration drift.

Accounting

Reconciling the organic mass balance across a continuous digester requires a multi-component matrix approach. This framework integrates chip feed rates, active liquor charges, dissolved biomass yield, and wash carryover into a solvable system of simultaneous linear equations, ensuring mass conservation holds across all chemical conversions and fluid separations in the plant.

Modeling the digester control volume with mass conservation matrices treats streams as multi-component vectors where each entry represents a specific species ~ oven-dry cellulose, hemicellulose, protolignin, extractives, sodium hydroxide, sodium sulfide, dissolved kraft lignin, organic acid salts, and water. At steady state, the sum of component mass flow rates across the system boundary equals zero.

Incomplete organic carbon balances in continuous digesters distort downstream steam generation estimates in recovery plant operations.

Matrix Formulation of Multi Component Mass Conservation

The total system mass balance is defined by the matrix vector equation:

A x = b

Where A is an m-by-n matrix of stoichiometric yields and phase separation coefficients, x is an n-dimensional vector of unknown component mass flow rates, and b is an m-dimensional vector of known input mass flow rates. Solving for unknown internal and exit flows requires inverting or decomposing matrix A. The individual biomass component balances are written as:

M_wood_OD = M_pulp_cellulose + M_pulp_hemicellulose + M_pulp_lignin + M_BL_lignin + M_BL_carbohydrate + M_BL_ext

Where M_pulp_i represents the residual mass of component i in unbleached pulp, and M_BL_i is the mass of component i dissolved into black liquor. The complete carbon balance closes when elemental carbon entering in wood and liquor additives equals the carbon leaving in pulp, black liquor, relief gases, and knotter rejects.

Worked Mass Balance Calculation for Industrial Digester

Consider a continuous digester processing 1000 metric tons per day of oven-dry wood chips containing 42.0 percent cellulose, 27.0 percent hemicellulose, 28.0 percent lignin, and 3.0 percent extractives. Target unbleached pulp yield is 47.0 percent at a Kappa number of 30.0 (4.41 percent residual lignin on pulp). White liquor EA charge is 18.0 percent as sodium oxide on OD wood at 30.0 percent sulfidity, with an effective alkali concentration of 105.0 grams per liter as sodium oxide.

Green wood moisture is 50.0 percent on a wet basis.

Step 1: Calculate incoming dry wood mass components per hour:

M_wood_OD = 1000.0 MT / 24.0 h = 41.67 MT/h

M_cell_in = 41.67 0.420 = 17.50 MT/h

M_hemi_in = 41.67 0.270 = 11.25 MT/h

M_lignin_in = 41.67 0.280 = 11.67 MT/h

M_ext_in = 41.67 0.030 = 1.25 MT/h

Step 2: Calculate unbleached pulp mass and component retention:

M_pulp_total = 41.67 0.470 = 19.58 MT/h

M_lignin_pulp = 19.58 (30.0 0.00147) = 0.86 MT/h

Assuming extractives are completely dissolved (M_ext_pulp = 0.0 MT/h) and cellulose retention is 95.0 percent of original cellulose input:

M_cell_pulp = 17.50 0.950 = 16.63 MT/h

M_hemi_pulp = M_pulp_total – M_lignin_pulp – M_cell_pulp = 19.58 – 0.86 – 16.63 = 2.09 MT/h

Step 3: Calculate dissolved organic solids mass generation rates entering the black liquor phase:

M_lignin_BL = M_lignin_in – M_lignin_pulp = 11.67 – 0.86 = 10.81 MT/h

M_cell_diss = M_cell_in – M_cell_pulp = 17.50 – 16.63 = 0.87 MT/h

M_hemi_diss = M_hemi_in – M_hemi_pulp = 11.25 – 2.09 = 9.16 MT/h

M_ext_diss = M_ext_in = 1.25 MT/h

M_carbo_diss_total = M_cell_diss + M_hemi_diss = 0.87 + 9.16 = 10.03 MT/h

Step 4: Calculate total dissolved organic mass flow rate in black liquor:

M_org_BL_total = M_lignin_BL + M_carbo_diss_total + M_ext_diss = 10.81 + 10.03 + 1.25 = 22.09 MT/h

For every 41.67 metric tons of dry wood processed per hour, 22.09 metric tons of dissolved organic solids pass into the black liquor stream. This dissolved organic mass comprises 48.9 percent kraft lignin, 45.4 percent carbohydrate degradation products, and 5.7 percent extractives.

• Mass Reconciliation Audit verifies that incoming biomass organic mass balances against the sum of pulp yield and weak black liquor dry organic solids within a two percent error margin.
• Alkali Stoichiometric Check calculates the ratio of consumed effective alkali to dissolved carbohydrate mass to verify that organic acid neutralization profiles match white liquor dosing rates.
• Organic Carbon Closure validates elemental carbon continuity from wood ultimate analysis down to weak black liquor elemental carbon analysis.
• Evaporator Mass Crosscheck verifies that the calculated total dissolved organic mass entering the evaporation train matches the measured dry solids mass rate leaving the heavy liquor storage tanks.

How does unmeasured lignin condensation alter the calculated ratio of dissolved organic carbon to inorganic sodium in the extraction black liquor line?

Recovery

Thermal conversion in the recovery furnace turns dissolved organic carbon into process steam while recovering cooking chemicals. The organic balance established at the digester dictates the combustion heat balance for the whole mill. Dissolved kraft lignin has a higher heating value of roughly 23.5 megajoules per kilogram of dry organic mass, compared to just 13.5 megajoules per kilogram for degraded carbohydrate hydroxy acids with their higher oxygen content and bound carboxyl groups.

Shifting delignification conditions in the digester alters the lignin-to-carbohydrate ratio in weak black liquor. Pushing effective alkali higher to hit lower Kappa numbers increases carbohydrate peeling, generating low-energy organic acid sodium salts at the expense of energy-dense lignin. This drops the heating value of the black liquor solids, requiring higher firing solids at the recovery boiler to maintain hearth temperatures and achieve proper sodium sulfate reduction.

Calorific Value Dependences on Dissolved Lignin Fractions

The gross calorific value of dry black liquor solids is modeled as a weighted linear combination of its organic and inorganic fractions. The higher heating value equation is:

HHV_BL = w_lignin HHV_lignin + w_acid HHV_acid + w_ext HHV_ext + w_inorg HHV_inorg

Where w_i is the dry mass fraction of component i in total black liquor solids, and HHV_i is the intrinsic higher heating value of that component. Since HHV_inorg is essentially zero, overall heating value depends entirely on organic mass fraction w_org and the proportion of lignin within that organic mass.

Higher lignin fractions raise liquor viscosity at elevated solids concentrations, affecting spray atomization at the recovery boiler liquor guns. Droplet size during spraying determines how liquor dries, pyrolyzes, and burns on the char bed. Liquor rich in carbohydrate degradation products flows with lower viscosity at identical solids levels, but produces lower hearth temperatures and less steam per ton of solids fired.

Chain of Custody Tracing for Organic Dissolved Biomass

Chain of custody standards under FSC and PEFC rules require precise tracking of organic carbon when processing mixtures of certified and uncertified wood. When certified chips enter the digester, yield calculations for pulp and black liquor set the upper limit on certified claims the mill can assign to market pulp, paper, or bio-based byproduct streams.

Under percentage or credit accounting methods, dissolved organic carbon burned in the recovery boiler counts as certified renewable bioenergy proportional to the mass fraction of certified wood input. Proving compliance requires a technical dossier connecting digester mass balance calculations directly to lab reports, flow meter calibrations, and yield audits.

• Wood Input Mass Manifest details the green mass, moisture content, and certified status of all wood chip species mixes fed into the digester impregnation vessel.
• Chemical Analysis Reports provide gravimetric total solids, organic to inorganic ratios, and elemental CHNS analyses for weak black liquor samples gathered according to TAPPI T 650 protocols.
• Digester Yield Reconciliation Logs record daily Kappa numbers, unbleached pulp production rates, reject masses, and calculated chip-to-pulp conversion factors.
• Chain of Custody Credit Ledger tracks the allocation of input certified biomass mass across pulp products, black liquor bioenergy production, and tall oil soap side-streams.

Claims covering bio-based byproducts like crude tall oil or lignin extracted from black liquor rely on these mass conservation balances to substantiate origin. Customs officials and environmental regulators audit these accounts to verify that certified output claims do not exceed physical mass limits. Closing the organic balance ensures every ton of dissolved biomass leaving the continuous digester is tracked accurately across chemical recovery, byproduct sales, and regulatory reporting.