1import math
2
3def calculate_combustion(fuel_type, excess_air_percentage, oxidizer="air"):
4    # 定义燃料成分
5    fuel_compositions = {
6        "methane": {"C": 1, "H": 4, "O": 0, "N": 0},
7        "ethane": {"C": 2, "H": 6, "O": 0, "N": 0},
8        "propane": {"C": 3, "H": 8, "O": 0, "N": 0},
9        "butane": {"C": 4, "H": 10, "O": 0, "N": 0},
10        "octane": {"C": 8, "H": 18, "O": 0, "N": 0}
11    }
12    
13    # 获取燃料成分
14    fuel = fuel_compositions.get(fuel_type)
15    if not fuel:
16        raise ValueError(f"未知的燃料类型: {fuel_type}")
17    
18    # 计算化学计量氧气需求
19    stoich_oxygen = fuel["C"] + fuel["H"]/4 - fuel["O"]/2
20    
21    # 计算实际氧气与过量空气
22    actual_oxygen = stoich_oxygen * (1 + excess_air_percentage/100)
23    
24    # 计算氮气（如果使用空气）
25    nitrogen = 0
26    if oxidizer == "air":
27        nitrogen = actual_oxygen * 3.76
28    
29    # 计算产物
30    co2 = fuel["C"]
31    h2o = fuel["H"]/2
32    excess_o2 = actual_oxygen - stoich_oxygen
33    n2 = nitrogen + fuel["N"]
34    
35    # 计算空气-燃料比
36    molar_masses = {"C": 12.01, "H": 1.008, "O": 16.00, "N": 14.01, "O2": 32.00, "N2": 28.01}
37    fuel_mass = (fuel["C"] * molar_masses["C"] + 
38                fuel["H"] * molar_masses["H"] + 
39                fuel["O"] * molar_masses["O"] + 
40                fuel["N"] * molar_masses["N"])
41    
42    air_mass = actual_oxygen * molar_masses["O2"] + nitrogen * molar_masses["N2"]
43    air_fuel_ratio = air_mass / fuel_mass
44    
45    return {
46        "air_fuel_ratio": air_fuel_ratio,
47        "products": {
48            "CO2": co2,
49            "H2O": h2o,
50            "O2": excess_o2,
51            "N2": n2
52        }
53    }
54
55# 示例用法
56result = calculate_combustion("methane", 20, "air")
57print(f"空气-燃料比: {result['air_fuel_ratio']:.2f}")
58print("产物（摩尔）:")
59for product, amount in result["products"].items():
60    print(f"  {product}: {amount:.2f}")
61

1function calculateCombustion(fuelType, excessAirPercentage, oxidizer = "air") {
2  // 摩尔质量
3  const molarMass = {
4    C: 12.01,
5    H: 1.008,
6    O: 16.00,
7    N: 14.01,
8    O2: 32.00,
9    N2: 28.01
10  };
11  
12  // 预定义燃料
13  const fuels = {
14    methane: { carbon: 1, hydrogen: 4, oxygen: 0, nitrogen: 0 },
15    ethane: { carbon: 2, hydrogen: 6, oxygen: 0, nitrogen: 0 },
16    propane: { carbon: 3, hydrogen: 8, oxygen: 0, nitrogen: 0 },
17    butane: { carbon: 4, hydrogen: 10, oxygen: 0, nitrogen: 0 },
18    octane: { carbon: 8, hydrogen: 18, oxygen: 0, nitrogen: 0 }
19  };
20  
21  const fuel = fuels[fuelType];
22  if (!fuel) {
23    throw new Error(`未知燃料类型: ${fuelType}`);
24  }
25  
26  // 计算化学计量氧气需求
27  const stoichOxygen = fuel.carbon + fuel.hydrogen/4.0 - fuel.oxygen/2.0;
28  
29  // 计算实际氧气与过量空气
30  const actualOxygen = stoichOxygen * (1 + excessAirPercentage/100.0);
31  
32  // 根据氧化剂计算氮气
33  let nitrogen = 0;
34  if (oxidizer === "air") {
35    nitrogen = actualOxygen * 3.76;
36  }
37  
38  // 计算产物
39  const products = {
40    CO2: fuel.carbon,
41    H2O: fuel.hydrogen/2,
42    O2: actualOxygen - stoichOxygen,
43    N2: nitrogen + fuel.nitrogen
44  };
45  
46  // 计算空气-燃料比
47  const fuelMass = (
48    fuel.carbon * molarMass.C +
49    fuel.hydrogen * molarMass.H +
50    fuel.oxygen * molarMass.O +
51    fuel.nitrogen * molarMass.N
52  );
53  
54  const airMass = actualOxygen * molarMass.O2 + nitrogen * molarMass.N2;
55  const airFuelRatio = airMass / fuelMass;
56  
57  // 生成平衡方程
58  const formatCoefficient = (num) => {
59    return num === 1 ? "" : num.toFixed(2);
60  };
61  
62  const equation = `${formatCoefficient(1)}C${fuel.carbon}H${fuel.hydrogen} + ${formatCoefficient(actualOxygen)}O₂ + ${formatCoefficient(nitrogen)}N₂ → ${formatCoefficient(products.CO2)}CO₂ + ${formatCoefficient(products.H2O)}H₂O + ${formatCoefficient(products.O2)}O₂ + ${formatCoefficient(products.N2)}N₂`;
63  
64  return {
65    balancedEquation: equation,
66    airFuelRatio,
67    products
68  };
69}
70
71// 示例用法
72const result = calculateCombustion("methane", 20, "air");
73console.log(`平衡方程: ${result.balancedEquation}`);
74console.log(`空气-燃料比: ${result.airFuelRatio.toFixed(2)}`);
75

1public class CombustionCalculator {
2    // 摩尔质量
3    private static final double C_MOLAR_MASS = 12.01;
4    private static final double H_MOLAR_MASS = 1.008;
5    private static final double O_MOLAR_MASS = 16.00;
6    private static final double N_MOLAR_MASS = 14.01;
7    private static final double O2_MOLAR_MASS = 32.00;
8    private static final double N2_MOLAR_MASS = 28.01;
9    
10    public static class FuelComposition {
11        public int carbon;
12        public int hydrogen;
13        public int oxygen;
14        public int nitrogen;
15        
16        public FuelComposition(int carbon, int hydrogen, int oxygen, int nitrogen) {
17            this.carbon = carbon;
18            this.hydrogen = hydrogen;
19            this.oxygen = oxygen;
20            this.nitrogen = nitrogen;
21        }
22    }
23    
24    public static class CombustionResult {
25        public String balancedEquation;
26        public double airFuelRatio;
27        public double co2;
28        public double h2o;
29        public double o2;
30        public double n2;
31    }
32    
33    public static CombustionResult calculateCombustion(FuelComposition fuel, 
34                                                      double excessAirPercentage, 
35                                                      boolean useAir) {
36        // 计算化学计量氧气
37        double stoichOxygen = fuel.carbon + fuel.hydrogen/4.0 - fuel.oxygen/2.0;
38        
39        // 计算实际氧气与过量空气
40        double actualOxygen = stoichOxygen * (1 + excessAirPercentage/100.0);
41        
42        // 根据氧化剂计算氮气
43        double nitrogen = 0;
44        if (useAir) {
45            nitrogen = actualOxygen * 3.76;
46        }
47        
48        // 计算产物
49        double co2 = fuel.carbon;
50        double h2o = fuel.hydrogen/2.0;
51        double excessO2 = actualOxygen - stoichOxygen;
52        double n2 = nitrogen + fuel.nitrogen;
53        
54        // 计算空气-燃料比
55        double fuelMass = (
56            fuel.carbon * C_MOLAR_MASS +
57            fuel.hydrogen * H_MOLAR_MASS +
58            fuel.oxygen * O_MOLAR_MASS +
59            fuel.nitrogen * N_MOLAR_MASS
60        );
61        
62        double airMass = actualOxygen * O2_MOLAR_MASS + nitrogen * N2_MOLAR_MASS;
63        double airFuelRatio = airMass / fuelMass;
64        
65        // 生成平衡方程
66        String equation = String.format(
67            "C%dH%d + %.2fO₂ + %.2fN₂ → %.2fCO₂ + %.2fH₂O + %.2fO₂ + %.2fN₂",
68            fuel.carbon, fuel.hydrogen, actualOxygen, nitrogen, 
69            co2, h2o, excessO2, n2
70        );
71        
72        // 创建并返回结果
73        CombustionResult result = new CombustionResult();
74        result.balancedEquation = equation;
75        result.airFuelRatio = airFuelRatio;
76        result.co2 = co2;
77        result.h2o = h2o;
78        result.o2 = excessO2;
79        result.n2 = n2;
80        
81        return result;
82    }
83    
84    public static void main(String[] args) {
85        // 示例：计算甲烷的燃烧，过量空气为20%
86        FuelComposition methane = new FuelComposition(1, 4, 0, 0);
87        CombustionResult result = calculateCombustion(methane, 20, true);
88        
89        System.out.println("平衡方程: " + result.balancedEquation);
90        System.out.println("空气-燃料比: " + String.format("%.2f", result.airFuelRatio));
91    }
92}
93