40 C.F.R. § 1065.643

Current through September 30, 2024

Section 1065.643 - Carbon balance error verification calculations

This section describes how to calculate quantities used in the carbon balance error verification described in § 1065.543 . Paragraphs (a) through (c) of this section describe how to calculate the mass of carbon for a test interval from carbon-carrying fluid streams, intake air into the system, and exhaust emissions, respectively. Paragraph (d) of this section describes how to use these carbon masses to calculate four different quantities for evaluating carbon balance error. Use rectangular or trapezoidal integration methods to calculate masses and amounts over a test interval from continuously measured or calculated mass and molar flow rates.

(a)Fuel and other fluids. Determine the mass of fuel, DEF, and other carbon-carrying fluid streams, other than intake air, flowing into the system, m_fluidj, for each test interval. Note that § 1065.543 allows you to omit all flows other than fuel. You may determine the mass of DEF based on ECM signals for DEF flow rate. You may determine fuel mass during field testing based on ECM signals for fuel flow rate. Calculate the mass of carbon from the combined carbon-carrying fluid streams flowing into the system as follows:

View Image

Where:

j = an indexing variable that represents one carbon-carrying fluid stream.

N = total number of carbon-carrying fluid streams into the system over the test interval.

w_C = carbon mass fraction of the carbon-carrying fluid stream as determined in § 1065.655(d) .

m_fluid = the mass of the carbon-carrying fluid stream determined over the test interval.

Example:

N = 2

w_Cfuel = 0.869

w_CDEF = 0.065

m_fuel = 1119.6 g

m_DEF = 36.8 g

m_Cfluid = 0.869·1119.6 + 0.065·36.8 = 975.3 g

(b)Intake air. Calculate the mass of carbon in the intake air, m_Cair, for each test interval using one of the methods in this paragraph (b). The methods are listed in order of preference. Use the first method where all the inputs are available for your test configuration. For methods that calculate m_Cair based on the amount of CO₂ per mole of intake air, we recommend measuring intake air concentration, but you may calculate x_CO2int using Eq. 1065.655-10 and letting x_CO2intdry = 375 [MICRO]mol/mol.

(1) Calculate m_Cair, using the following equation if you measure intake air flow:

View Image

Where:

M_C = molar mass of carbon.

n_int = measured amount of intake air over the test interval.

x_CO2int = amount of intake air CO₂ per mole of intake air.

Example:

M_C = 12.0107 g/mol

n_int = 62862 mol

x_CO2int = 369 [MICRO]mol/mol = 0.000369 mol/mol

m_Cair = 12.0107·62862·0.000369 = 278.6 g

(2) Calculate m_Cair, using the following equation if you measure or calculate raw exhaust flow and you calculate chemical balance terms:

View Image

Where:

M_C = molar mass of carbon.

n_exh = calculated or measured amount of raw exhaust over the test interval.

x_H2Oexh = amount of H₂O in exhaust per mole of exhaust.

x_CO2int = amount of intake air CO₂ per mole of intake air.

x_dil/exhdry = amount of excess air per mole of dry exhaust. Note that excess air and intake air have the same composition, so x_CO2dil = x_CO2int and x_H2Odil = x_H2Oint for the chemical balance calculation for raw exhaust.

x_int/exhdry = amount of intake air required to produce actual combustion products per mole of dry exhaust.

Example:

M_C = 12.0107 g/mol

n_exh = 62862 mol

x_H2Oexh = 0.034 mol/mol

x_CO2int = 369 [MICRO]mol/mol = 0.000369 mol/mol

x_dil/exhdry = 0.570 mol/mol

x_int/exhdry = 0.465 mol/mol

m_Cair = 12.0107·62862·(1 - 0.034)·0.000369·(0.570 + 0.465) = 278.6 g

(3) Calculate m_Cair, using the following equation if you measure raw exhaust flow:

View Image

Where:

M_C = molar mass of carbon.

n_exh = measured amount of raw exhaust over the test interval.

x_CO2int = amount of intake air CO₂ per mole of intake air.

Example:

M_C = 12.0107 g/mol

n_exh = 62862 mol

x_CO2int = 369 [MICRO]mol/mol = 0.000369 mol/mol

m_Cair = 12.0107·62862·0.000369 = 278.6 g

(4) Calculate m_Cair, using the following equation if you measure diluted exhaust flow and dilution air flow:

View Image

Where:

M_C = molar mass of carbon.

n_dexh = measured amount of diluted exhaust over the test interval as determined in § 1065.642 .

n_dil = measured amount of dilution air over the test interval as determined in § 1065.667(b) .

x_CO2int = amount of intake air CO₂ per mole of intake air.

Example:

M_C = 12.0107 g/mol

n_dexh = 942930 mol

n_dil = 880068 mol

x_CO2int = 369 [MICRO]mol/mol = 0.000369 mol/mol

m_Cair = 12.0107·(942930 - 880068)·0.000369 = 278.6 g

(5) Determined m_Cair based on ECM signals for intake air flow as described in paragraph (b)(1) of this section.

(6) If you measure diluted exhaust, determine m_Cair as described in paragraph (b)(4) of this section using a calculated amount of dilution air over the test interval as determined in § 1065.667(d) instead of the measured amount of dilution air.

(c)Exhaust emissions. Calculate the mass of carbon in exhaust emissions, m_Cexh, for each test interval as follows:

View Image

Where:

M_C = molar mass of carbon.

m_CO2 = mass of CO₂ over the test interval as determined in § 1065.650(c) .

M_CO2 = molar mass of carbon dioxide.

m_CO = mass of CO over the test interval as determined in § 1065.650(c) .

M_CO = molar mass of carbon monoxide.

m_THC = mass of THC over the test interval as determined in § 1065.650(c) .

M_THC = effective C₁ molar mass of total hydrocarbon as defined in § 1065.1005(f)(2) .

Example:

M_C = 12.0107 g/mol

m_CO2 = 4567 g

M_CO2 = 44.0095 g/mol

m_CO = 0.803 g

M_CO = 28.0101 g/mol

m_THC = 0.537 g

M_THC = 13.875389 g/mol

View Image

(d)Carbon balance error quantities. Calculate carbon balance error quantities as follows:

(1) Calculate carbon mass absolute error, [EPSILON]_aC, for a test interval as follows:

View Image

Where:

m_Cexh = mass of carbon in exhaust emissions over the test interval as determined in paragraph (d) of this section.

m_Cfluid = mass of carbon in all the carbon-carrying fluid streams flowing into the system over the test interval as determined in paragraph (a) of this section.

m_Cair = mass of carbon in the intake air flowing into the system over the test interval as determined in paragraph (b) of this section.

Example:

= 1247.2 g

= 975.3 g

= 278.6 g

= 1247.2 - 975.3 - 278.6

= -6.7 g

(2) Calculate carbon mass rate absolute error, [EPSILON]_aCrate, for a test interval as follows:

View Image

Where:

t = duration of the test interval.

Example:

= -6.7 g

t = 1202.2 s = 0.3339 hr

View Image

= -20.065 g/hr

(3) Calculate carbon mass relative error, [EPSILON]_rC, for a test interval as follows:

View Image

Example:

= -6.7 g

= 975.3 g

= 278.6 g

View Image

= -0.0053

(4) Calculate composite carbon mass relative error, [EPSILON]_rCcomp, for a duty cycle with multiple test intervals as follows:

(i) Calculate [EPSILON]_rCcomp using the following equation:

View Image

Where:

i = an indexing variable that represents one test interval.

N = number of test intervals.

WF = weighting factor for the test interval as defined in the standard-setting part.

m_Cexh = mass of carbon in exhaust emissions over the test interval as determined in paragraph (c) of this section.

m_Cfluid = mass of carbon in all the carbon-carrying fluid streams that flowed into the system over the test interval as determined in paragraph (a) of this section.

m_Cair = mass of carbon in the intake air that flowed into the system over the test interval as determined in paragraph (b) of this section.

t = duration of the test interval. For duty cycles with multiple test intervals of a prescribed duration, such as cold-start and hot-start transient cycles, set t = 1 for all test intervals. For discrete-mode steady-state duty cycles with multiple test intervals of varying duration, set t equal to the actual duration of each test interval.

(ii) The following example illustrates calculation of [EPSILON]_rCcomp, for cold-start and hot-start transient cycles:

N = 2

= ¹/₇

= ⁶/₇

= 1255.3 g

= 1247.2 g

= 977.8 g

= 975.3 g

= 280.2 g

= 278.6 g

View Image

= -0.0049

(iii) The following example illustrates calculation of [EPSILON]_rCcomp for multiple test intervals with varying duration, such as discrete-mode steady-state duty cycles:

N = 2

= 0.85

= 0.15

= 2.873 g

= 0.125 g

= 2.864 g

= 0.095 g

= 0.023 g

= 0.024 g

= 123 s

= 306 s

View Image

= -0.0047

40 C.F.R. §1065.643

86 FR 34557, June 29, 2021, as amended at 87 FR 64865, Oct. 26, 2022; 88 FR 4680, Jan. 24, 2023

86 FR 34557, 6/29/2021; as amended at 87 FR 64865, 10/26/2022; 88 FR 4680, 3/27/2023