40 CFR § 1065.655 - Chemical balances of fuel, DEF, intake air, and exhaust.

§ 1065.655 Chemical balances of fuel, DEF, intake air, and exhaust.

(a) General. Chemical balances of fuel, intake air, and exhaust may be used to calculate flows, the amount of water in their flows, and the wet concentration of constituents in their flows. With one flow rate of either fuel, intake air, or exhaust, you may use chemical balances to determine the flows of the other two. For example, you may use chemical balances along with either intake air or fuel flow to determine raw exhaust flow. Note that chemical balance calculations allow measured values for the flow rate of diesel exhaust fluid for engines with urea-based selective catalytic reduction.

(b) Procedures that require chemical balances. We require chemical balances when you determine the following:

(1) A value proportional to total work, W

, when you choose to determine brake-specific emissions as described in § 1065.650(f).

(2) Raw exhaust molar flow rate either from measured intake air molar flow rate or from fuel mass flow rate as described in paragraph (f) of this section.

(3) Raw exhaust molar flow rate from measured intake air molar flow rate and dilute exhaust molar flow rate, as described in paragraph (g) of this section.

(4) The amount of water in a raw or diluted exhaust flow, χH2Oexh, when you do not measure the amount of water to correct for the amount of water removed by a sampling system. Correct for removed water according to § 1065.659.

(5) The calculated total dilution air flow when you do not measure dilution air flow to correct for background emissions as described in § 1065.667(c) and (d).

(c) Chemical balance procedure. The calculations for a chemical balance involve a system of equations that require iteration. We recommend using a computer to solve this system of equations. You must guess the initial values of up to three quantities: the amount of water in the measured flow, xH2Oexh, fraction of dilution air in diluted exhaust, xdil/exh, and the amount of products on a C1 basis per dry mole of dry measured flow, xCcombdry. You may use time-weighted mean values of intake air humidity and dilution air humidity in the chemical balance; as long as your intake air and dilution air humidities remain within tolerances of ±0.0025 mol/mol of their respective mean values over the test interval. For each emission concentration, x, and amount of water, xH2Oexh, you must determine their completely dry concentrations, xdry and xH2Oexhdry. You must also use your fuel mixture's atomic hydrogen-to-carbon ratio, α, oxygen-to-carbon ratio, β, sulfur-to-carbon ratio, γ, and nitrogen-to-carbon ratio, δ; you may optionally account for diesel exhaust fluid (or other fluids injected into the exhaust), if applicable. You may calculate α, β,γ, and δ based on measured fuel composition or based on measured fuel and diesel exhaust fluid (or other fluids injected into the exhaust) composition together, as described in paragraph (e) of this section. You may alternatively use any combination of default values and measured values as described in paragraph (e) of this section. Use the following steps to complete a chemical balance:

(1) Convert your measured concentrations such as, xCO2meas, xNOmeas, and xH2Oint, to dry concentrations by dividing them by one minus the amount of water present during their respective measurements; for example: xH2OxCO2meas, xH2OxNOmeas, and xH2Oint. If the amount of water present during a “wet” measurement is the same as the unknown amount of water in the exhaust flow, xH2Oexh, iteratively solve for that value in the system of equations. If you measure only total NOX and not NO and NO2 separately, use good engineering judgment to estimate a split in your total NOX concentration between NO and NO2 for the chemical balances. For example, if you measure emissions from a stoichiometric spark-ignition engine, you may assume all NOX is NO. For a compression-ignition engine, you may assume that your molar concentration of NOX, xNOx, is 75% NO and 25% NO2. For NO2 storage aftertreatment systems, you may assume xNOx is 25% NO and 75% NO2. Note that for calculating the mass of NOX emissions, you must use the molar mass of NO2 for the effective molar mass of all NOX species, regardless of the actual NO2 fraction of NOX.

(2) Enter the equations in paragraph (c)(4) of this section into a computer program to iteratively solve for xH2Oexh, xCcombdry, and xdil/exh. Use good engineering judgment to guess initial values for xH2Oexh, xCcombdry, and xdil/exh. We recommend guessing an initial amount of water that is about twice the amount of water in your intake or dilution air. We recommend guessing an initial value of xCcombdry as the sum of your measured CO2, CO, and THC values. We also recommend guessing an initial xdil/exh between 0.75 and 0.95, such as 0.8. Iterate values in the system of equations until the most recently updated guesses are all within ±1% of their respective most recently calculated values.

(3) Use the following symbols and subscripts in the equations for performing the chemical balance calculations in this paragraph (c):

Table 1 of § 1065.655—Symbols and Subscripts for Chemical Balance Equations

x dil/exh Amount of dilution gas or excess air per mole of exhaust
X dil/exh amount of dilution gas or excess air per mole of exhaust.
x Ccombdry amount of carbon from fuel and any injected fluids in the exhaust per mole of dry exhaust
x H2dry amount of H2 in exhaust per amount of dry exhaust
K H2Ogas water-gas reaction equilibrium coefficient; you may use 3.5 or calculate your own value using good engineering judgment
x H2Oexhdry amount of H2O in exhaust per dry mole of dry exhaust
x prod/intdry amount of dry stoichiometric products per dry mole of intake air
x dil/exhdry amount of dilution gas and/or excess air per mole of dry exhaust
x int/exhdry amount of intake air required to produce actual combustion products per mole of dry (raw or diluted) exhaust
x raw/exhdry amount of undiluted exhaust, without excess air, per mole of dry (raw or diluted) exhaust
x O2int amount of intake air O2 per mole of intake air
x CO2intdry amount of intake air CO2 per mole of dry intake air; you may use x CO2intdry = 375 µmol/mol, but we recommend measuring the actual concentration in the intake air
x H2Ointdry amount of intake air H2O per mole of dry intake air
x CO2int amount of intake air CO2 per mole of intake air
x CO2dil amount of dilution gas CO2 per mole of dilution gas
x CO2dildry amount of dilution gas CO2 per mole of dry dilution gas; if you use air as diluent, you may use x CO2dildry = 375 µmol/mol, but we recommend measuring the actual concentration in the intake air
x H2Odildry amount of dilution gas H2O per mole of dry dilution gas
x H2Odil amount of dilution gas H2O per mole of dilution gas
x [emission]meas amount of measured emission in the sample at the respective gas analyzer
x [emission]dry amount of emission per dry mole of dry sample
x H2O[emission]meas amount of H2O in sample at emission-detection location; measure or estimate these values according to § 1065.145(e)(2)
x H2Oint amount of H2O in the intake air, based on a humidity measurement of intake air
α atomic hydrogen-to-carbon ratio of the fuel (or mixture of test fuels) and any injected fluids
β atomic oxygen-to-carbon ratio of the fuel (or mixture of test fuels) and any injected fluids
γ atomic sulfur-to-carbon ratio of the fuel (or mixture of test fuels) and any injected fluids
δ atomic nitrogen-to-carbon ratio of the fuel (or mixture of test fuels) and any injected fluids

(4) Use the following equations to iteratively solve for xdil/exh, xH2Oexh, and xCcombdry:

x dil / exh = 1 x raw / exhdry 1 + x H 2 Oexhdry Eq. 1065.655-1

x H 2 Oexh = x H 2 Oexhdry 1 + x H 2 Oexhdry Eq. 1065.655-2

x Ccombdry = x CO 2 dry + x COdry + x THCdry x CO 2 dil · x dil / exhdry x CO 2 int · x int / exhdry Eq. 1065.655-3

x H 2 dry = x COdry · ( x H 2 Oexhdry x H 2 Odil · x dil / exhdry ) K H2O-gas · ( x CO 2 dry x CO 2 dil · x dil / exhdry ) Eq. 1065.655-4

x H 2 Oexhdry = α 2 ( x Coombdry x THCdry ) + x H 2 Odil · x dil / exhdry + x H 2 Oint · x int / exhdry x H 2 dry Eq. 1065.655-5

x dil / exhdry = x dil / exh 1 x H 2 Oexh Eq. 1065.655-6

x int / exhdry = 1 2 · x O 2 int ( ( α 2 β + 2 + 2 γ ) ( x Ccombdry x THCdry ) ( x COdry x NOdry 2 x NO 2 dry + x H 2 dry ) ) Eq. 1065.655-7

x raw / exhdry = 1 2 ( ( α 2 + β + δ ) ( x Ccombdry x THCdry ) + ( 2 x THCdry + x COdry x NO 2 dry + x H 2 dry ) ) + x int / exhdry Eq. 1065.655-8

x O 2 int = 0.209820 x CO 2 intdry 1 + x H 2 O intdry Eq. 1065.655-9

x CO 2 int = x CO 2 intdry 1 + x H 2 O intdry Eq. 1065.655-10

x H 2 O intdry = x H 2 O int 1 x H 2 O int Eq. 1065.655-11

x CO 2 dil = x CO 2 dildry 1 + x H 2 O dildry Eq. 1065.655-12

x H 2 O dildry = x H 2 O dil 1 x H 2 O dil Eq. 1065.655-13

x CO dry = x CO meas 1 x H 2 O CO meas Eq. 1065.655-14

x CO 2 dry = x CO 2 meas 1 x H 2 O CO 2 meas Eq. 1065.655-15

x NO dry = x NO meas 1 x H 2 O NO meas Eq. 1065.655-16

x NO 2 dry = x NO 2 meas 1 x H 2 O NO 2 meas Eq. 1065.655-17

x THC dry = x THC meas 1 x H 2 O THC meas Eq. 1065.655-18

(5) The following example is a solution for xdil/exh,x, xH2Oexh, and xCcombdry using the equations in paragraph (c)(4) of this section:

x dil / exh = 1 0.184 1 + 35.37 1000 = 0.822 mol / mol

x H 2 Oexh = 35.37 1 + 35.37 1000 = 34.16 mmol / mol x Ccombdry = 0.0252 + 29.3 1000000 + 47.6 1000000 0.371 1000 · 0.851 0.369 1000 · 0.172 = 0.0249 mol / mol x H 2 dry = 29.3 · ( 0.034 0.012 · 0.851 ) 3.5 · ( 25.2 1000 0.371 1000 · 0.851 ) = 8.5 μ mol / mol x H 2 Oexhdry = 1.8 2 ( 0.0249 47.6 1000000 ) + 0.012 · 0.851 + 0.017 · 0.172 8.5 1000000 = 0.0353 mol / mol x dil / exhdry = 0.822 1 0.034 = 0.851 mol / mol x int / exhdry = 1 2 · 0.206 ( 1.8 2 0.050 + 2 + 2 · 0.0003 ) ( 0.0249 47.6 1000000 ) 29.3 1000000 50.4 1000000 2 · 12.1 1000000 + 8.5 1000000 = 0.172 mol / mol x raw / exhdry = 1 2 ( ( 1.8 2 + 0.050 + 0.0001 ) ( 0.0249 47.6 1000000 ) + 2 · 47.6 1000000 + 29.3 1000000 12.1 1000000 + 8.5 1000000 ) + 0.172 = 0.184 mol / mol x O 2 int = 0.209820 0.000375 1 + 17.22 1000 = 0.206 mol / mol x CO 2 int = 0.000375 · 1000 1 + 17.22 1000 = 0.369 mmol / mol x H 2 Ointdry = 16.93 1 16.93 1000 = 17.22 mmol / mol

x CO 2 dil = 0.375 1 + 12.01 1000 = 0.371 mmol / mol x H 2 Odildry = 11.87 1 11.87 1000 = 12.01 mmol / mol x COdry = 29.0 1 8.601 1000 = 29.3 mmol / mol x CO 2 dryl = 24.98 1 8.601 1000 = 25.2 mmol / mol x NOdry = 50.0 1 8.601 1000 = 50.4 mmol / mol x NO 2 dry = 12.0 1 8.601 1000 = 12.1 mmol / mol x THCdry = 46 1 34.16 1000 = 47.6 mmol / mol

a = 1.8
b = 0.05
g = 0.0003
d = 0.0001

(d) Carbon mass fraction of fuel. Determine carbon mass fraction of fuel, wC, based on the fuel properties as determined in paragraph (e) of this section, optionally accounting for diesel exhaust fluid's contribution to α, β, γ, and δ, or other fluids injected into the exhaust, if applicable (for example, the engine is equipped with an emission control system that utilizes DEF). Calculate wC using the following equation:

w C = 1 · M C 1 · M C + α · M H + β · M O + γ · M S + Δ · M N Eq. 1065.655-19

Where:
wC = carbon mass fraction of the fuel (or mixture of test fuels) and any injected fluids.
MC = molar mass of carbon.
a = atomic hydrogen-to-carbon ratio of the fuel (or mixture of test fuels) and any injected fluids.
MH = molar mass of hydrogen.
b = atomic oxygen-to-carbon ratio of the fuel (or mixture of test fuels) and any injected fluids.
MO = molar mass of oxygen.
g = atomic sulfur-to-carbon ratio of the fuel (or mixture of test fuels) and any injected fluids.
MS = molar mass of sulfur.
d = atomic nitrogen-to-carbon ratio of the fuel (or mixture of test fuels) and any injected fluids.
MN = molar mass of nitrogen.
Example:
a = 1.8
b = 0.05
g = 0.0003
d = 0.0001
MC = 12.0107
MH = 1.00794
MO = 15.9994
MS = 32.065
MN = 14.0067

w c = 1 · 12.0107 1 · 12.0107 + 1.8 · 1.00794 + 0.05 · 15.9994 + 0.0003 · 32.065 + 0.0001 · 14.0067

wC = 0.8206

(e) Fuel and diesel exhaust fluid composition. Determine fuel and diesel exhaust fluid composition represented by α, β, γ, and δ as described in this paragraph (e). When using measured fuel or diesel exhaust fluid properties, you must determine values for α and β in all cases. If you determine compositions based on measured values and the default value listed in Table 2 of this section is zero, you may set γ and δ to zero; otherwise determine γ and δ (along with α and β) based on measured values. Determine elemental mass fractions and values for α, β, γ, and δ as follows:

(1) For liquid fuels, use the default values for α, β, γ, and δ in Table 2 of this section or determine mass fractions of liquid fuels for calculation of α, β, γ, and δ as follows:

(i) Determine the carbon and hydrogen mass fractions according to ASTM D5291 (incorporated by reference in § 1065.1010). When using ASTM D5291 to determine carbon and hydrogen mass fractions of gasoline (with or without blended ethanol), use good engineering judgment to adapt the method as appropriate. This may include consulting with the instrument manufacturer on how to test high-volatility fuels. Allow the weight of volatile fuel samples to stabilize for 20 minutes before starting the analysis; if the weight still drifts after 20 minutes, prepare a new sample). Retest the sample if the carbon, hydrogen, oxygen, sulfur, and nitrogen mass fractions do not add up to a total mass of 100 ±0.5%; you may assume oxygen has a zero mass contribution for this specification for diesel fuel and neat (E0) gasoline. You may also assume that sulfur and nitrogen have a zero mass contribution for this specification for all fuels except residual fuel blends.

(ii) Determine oxygen mass fraction of gasoline (with or without blended ethanol) according to ASTM D5599 (incorporated by reference in § 1065.1010). For all other liquid fuels, determine the oxygen mass fraction using good engineering judgment.

(iii) Determine the nitrogen mass fraction according to ASTM D4629 or ASTM D5762 (incorporated by reference in § 1065.1010) for all liquid fuels. Select the correct method based on the expected nitrogen content.

(iv) Determine the sulfur mass fraction according to subpart H of this part.

(2) For gaseous fuels and diesel exhaust fluid, use the default values for α, β, γ, and δ in Table 2 of this section, or use good engineering judgment to determine those values based on measurement.

(3) For nonconstant fuel mixtures, you must account for the varying proportions of the different fuels. This paragraph (e)(3) generally applies for dual-fuel and flexible-fuel engines, but it also applies if diesel exhaust fluid is injected in a way that is not strictly proportional to fuel flow. Account for these varying concentrations either with a batch measurement that provides averaged values to represent the test interval, or by analyzing data from continuous mass rate measurements. Application of average values from a batch measurement generally applies to situations where one fluid is a minor component of the total fuel mixture, for example dual-fuel and flexible-fuel engines with diesel pilot injection, where the diesel pilot fuel mass is less than 5% of the total fuel mass and diesel exhaust fluid injection; consistent with good engineering judgment.

(4) Calculate α, β,γ, and δ using the following equations:

Where:
N = total number of fuels and injected fluids over the duty cycle.
j = an indexing variable that represents one fuel or injected fluid, starting with j = 1.
m
j = the mass flow rate of the fuel or any injected fluid j. For applications using a single fuel and no DEF fluid, set this value to 1. For batch measurements, divide the total mass of fuel over the test interval duration to determine a mass rate.
wHj = hydrogen mass fraction of fuel or any injected fluid j.
wCj = carbon mass fraction of fuel or any injected fluid j.
wOj = oxygen mass fraction of fuel or any injected fluid j.
wSj = sulfur mass fraction of fuel or any injected fluid j.
wNj = nitrogen mass fraction of fuel or any injected fluid j.
Example:
N = 1
j = 1
m
1 = 1
wH1 = 0.1239
wC1 = 0.8206

wO1 = 0.0547

wS1 = 0.00066

wN1 = 0.000095

MC = 12.0107 g/mol

MH = 1.00794 g/mol

MO = 15.9994 g/mol

MS = 32.065 g/mol

MN = 14.0067

α = 1.799
β = 0.05004
γ = 0.0003012
δ = 0.0001003

(5) Table 2 follows:

(f) Calculated raw exhaust molar flow rate from measured intake air molar flow rate or fuel mass flow rate. You may calculate the raw exhaust molar flow rate from which you sampled emissions, n

exh, based on the measured intake air molar flow rate, n
int, or the measured fuel mass flow rate, m
fuel, and the values calculated using the chemical balance in paragraph (c) of this section. The chemical balance must be based on raw exhaust gas concentrations. Solve for the chemical balance in paragraph (c) of this section at the same frequency that you update and record or n
int or m
fuel. For laboratory tests, calculating raw exhaust molar flow rate using measured fuel mass flow rate is valid only for steady-state testing. See § 1065.915(d)(5)(iv) for application to field testing.

(1) Crankcase flow rate. If engines are not subject to crankcase controls under the standard-setting part, you may calculate raw exhaust flow based on n

int or m
fuel using one of the following:

(i) You may measure flow rate through the crankcase vent and subtract it from the calculated exhaust flow.

(ii) You may estimate flow rate through the crankcase vent by engineering analysis as long as the uncertainty in your calculation does not adversely affect your ability to show that your engines comply with applicable emission standards.

(iii) You may assume your crankcase vent flow rate is zero.

(2) Intake air molar flow rate calculation. Calculate n

exh based on n
int using the following equation:

n · exh = n · int ( 1 + ( x int / exhdry x raw / exhdry ) ( 1 + x H 2 Oexhdry ) ) Eq. 1065.655-24

Where:
n
exh = raw exhaust molar flow rate from which you measured emissions.
n
int = intake air molar flow rate including humidity in intake air.
Example:
n
int = 3.780 mol/s
xint/exhdry = 0.69021 mol/mol
xraw/exhdry = 1.10764 mol/mol
xH20exhdry = 107.64 mmol/mol = 0.10764 mol/mol

n · exh = 3.780 ( 1 + ( 0.69021 1.10764 ) ( 1 + 0.10764 ) ) n · exh = 6.066 mol / s (3) Fuel mass flow rate calculation may be used only for steady-state laboratory testing. See § 1065.915(d)(5)(iv) for application to field testing. Calculate n · exh based on m · fuel using the following equation: n · exh = m · fuel · w c · ( 1 + x H 2 Oexhdry ) M c · x Ccombdry Eq. 1065.655-25

(3) Fluid mass flow rate calculation. This calculation may be used only for steady-state laboratory testing. You may not use this calculation if the standard-setting part requires carbon balance error verification as described in § 1065.543. See § 1065.915(d)(5)(iv) for application to field testing. Calculate n

exh based on m
j using the following equation:

Where:
n
exh = raw exhaust molar flow rate from which you measured emissions.
j = an indexing variable that represents one fuel or injected fluid, starting with j = 1.
N = total number of fuels and injected fluids over the duty cycle.
m
j = the mass flow rate of the fuel or any injected fluid j.
wCj = carbon mass fraction of the fuel and any injected fluid j.
Example:
N = 1
j = 1
m
1 = 7.559 g/s
wC1 = 0.869 g/g
MC = 12.0107 g/mol
xCcombdry1 = 99.87 mmol/mol = 0.09987 mol/mol
xH20exhdry1 = 107.64 mmol/mol = 0.10764 mol/mol
nexh = 6.066 mol/s

(g) Calculated raw exhaust molar flow rate from measured intake air molar flow rate, dilute exhaust molar flow rate, and dilute chemical balance. You may calculate the raw exhaust molar flow rate, n

exh, based on the measured intake air molar flow rate, n
int, the measured dilute exhaust molar flow rate, n
dexh, and the values calculated using the chemical balance in paragraph (c) of this section. Note that the chemical balance must be based on dilute exhaust gas concentrations. For continuous-flow calculations, solve for the chemical balance in paragraph (c) of this section at the same frequency that you update and record n
int and n
dexh. This calculated n
exh may be used for the PM dilution ratio verification in § 1065.546; the calculation of dilution air molar flow rate in the background correction in § 1065.667; and the calculation of mass of emissions in § 1065.650(c) for species that are measured in the raw exhaust.

(1) Crankcase flow rate. If engines are not subject to crankcase controls under the standard-setting part, calculate raw exhaust flow as described in paragraph (e)(1) of this section.

(2) Dilute exhaust and intake air molar flow rate calculation. Calculate n

exh as follows:

n · exh = ( x raw / exhdry x int / exhdry ) · ( 1 x H 2 Oexh ) · n · dexh + n · int Eq. 1065.655-26

Example:
n
int = 7.930 mol/s
xraw/exhdry = 0.1544 mol/mol
xint/exhdry = 0.1451 mol/mol
xH20/exh = 32.46 mmol/mol = 0.03246 mol/mol
n
dexh = 49.02 mol/s
n
exh = (0.1544 − 0.1451) · (1 − 0.03246) · 49.02 + 7.930 = 0.4411 + 7.930 = 8.371 mol/s
[73 FR 37331, June 30, 2008, as amended at 73 FR 59334, Oct. 8, 2008; 75 FR 23051, Apr. 30, 2010; 76 FR 57458, Sept. 15, 2011; 79 FR 23799, Apr. 28, 2014; 81 FR 74182, Oct. 25, 2016; 86 FR 34563, June 29, 2021; 87 FR 64866, Oct. 26, 2022; 88 FR 4684, Jan. 24, 2023]