我正在尝试计算化学反应产物的体积分数。我已经弄清楚了空气体积分数的计算。
我知道化学计量系数的公式是从平衡基于 C、H、O、N 原子的化学反应中提取的。我不明白的部分是 CASE(3) 的最后一部分（包括图片）。
公式如下：VOLUME_FRACTION( 3 ) = RN % NU_H2O + SPECIES_MIXTURE( 1 ) % VOLUME_FRACTION(H2O_INDEX)*RN % NU_O2 / &
SPECIES_MIXTURE( 1 ) % VOLUME_FRACTION(O2_INDEX)。
公式中突出显示的部分是否代表空气转化为水蒸气的部分？
有人可以解释一下，为什么突出显示的部分包含在公式中？
代码:
CASE(3)
IF (N_SIMPLE_CHEMISTRY_REACTIONS==1) RN=>REACTION(1)
IF (N_SIMPLE_CHEMISTRY_REACTIONS==2) RN=>REACTION(2)
ID = 'PRODUCTS'
FORMULA = 'Z2'
SPEC_ID(1) = 'CARBON MONOXIDE'
SPEC_ID(2) = 'SOOT'
SPEC_ID(3) = 'WATER VAPOR'
SPEC_ID(4) = 'CARBON DIOXIDE'
SPEC_ID(5) = 'HYDROGEN CYANIDE'
SPEC_ID(6) = 'NITROGEN'
RN%NU_CO = (SPECIES_MIXTURE(FUEL_SMIX_INDEX)%MW/MW_CO) *REACTION(1)%CO_YIELD
RN%NU_HCN = (SPECIES_MIXTURE(FUEL_SMIX_INDEX)%MW/MW_HCN) *REACTION(1)%HCN_YIELD
RN%NU_SOOT = (SPECIES_MIXTURE(FUEL_SMIX_INDEX)%MW/RN%MW_SOOT)*REACTION(1)%SOOT_YIELD
RN%NU_H2O = 0.5_EB*REACTION(1)%H - 0.5_EB*(RN%NU_SOOT*REACTION(1)%SOOT_H_FRACTION + RN%NU_HCN)
IF (ABS(RN%NU_H2O) < TWO_EPSILON_EB) RN%NU_H2O = 0._EB
RN%NU_CO2 = REACTION(1)%C - RN%NU_CO - RN%NU_HCN - RN%NU_SOOT*(1._EB-REACTION(1)%SOOT_H_FRACTION)
IF (ABS(RN%NU_CO2) < TWO_EPSILON_EB) RN%NU_CO2 = 0._EB
RN%NU_O2 = RN%NU_CO2 + 0.5_EB*(RN%NU_CO+RN%NU_H2O-REACTION(1)%O)
RN%NU_N2 = 0.5_EB*(REACTION(1)%N - RN%NU_HCN)
VOLUME_FRACTION(1) = RN%NU_CO
VOLUME_FRACTION(2) = RN%NU_SOOT
VOLUME_FRACTION(3) = RN%NU_H2O + SPECIES_MIXTURE(1)%VOLUME_FRACTION(H2O_INDEX)*RN%NU_O2 / &
SPECIES_MIXTURE(1)%VOLUME_FRACTION(O2_INDEX)
VOLUME_FRACTION(4) = RN%NU_CO2 + SPECIES_MIXTURE(1)%VOLUME_FRACTION(CO2_INDEX)*RN%NU_O2 / &
SPECIES_MIXTURE(1)%VOLUME_FRACTION(O2_INDEX)
VOLUME_FRACTION(5) = RN%NU_HCN
VOLUME_FRACTION(6) = RN%NU_N2 + SPECIES_MIXTURE(1)%VOLUME_FRACTION(N2_INDEX)*RN%NU_O2 / &
SPECIES_MIXTURE(1)%VOLUME_FRACTION(O2_INDEX)
VOLUME_FRACTION = VOLUME_FRACTION/SUM(VOLUME_FRACTION)
SPECIES(SOOT_INDEX)%ATOMS=0._EB
SPECIES(SOOT_INDEX)%ATOMS(1)=REACTION(1)%SOOT_H_FRACTION
SPECIES(SOOT_INDEX)%ATOMS(6)=1._EB-REACTION(1)%SOOT_H_FRACTION