온실가스 $CO_2$를 대기 중에서 포집하여 신재생 전력으로 석유화학 산업의 쌀인 에틸렌($C_2H_4$)으로 직접 합성하는 탄소 포집 및 활용(CCU) 기술의 핵심은 다전자-다양성자 전이 반응을 유도하는 구리(Cu) 촉매 표면의 선택성입니다. 단일 탄소 화합물(일산화탄소, 메탄)로의 분기를 막고 두 개의 일산화탄소 중간체($*CO$)가 결합하는 C-C 커플링 활성화 에너지를 낮추는 표면 흡착 화학과 기체 확산 전극의 국소 pH 효과를 1,000Hz 주기로 0-Alloc C# 수치해석하는 온디바이스 탄소중립 반응기 시뮬레이터를 공개합니다.

1. 2022 개정 교육과정 교과 연계 분석

  • 화학 (탄소 화합물의 산화-환원 반응, 촉매 표면 흡착 메커니즘, 유기 분자 합성): $CO_2$의 탄소 원자 산화수(+4)가 에틸렌($C_2H_4$)의 탄소(-2)로 환원되는 12전자 전이 반응 $2\text{CO}_2 + 8\text{H}_2\text{O} + 12e^- \to \text{C}_2\text{H}_4 + 12\text{OH}^-$, 구리 나노 입자 표면의 랭뮤어-힌셜우드 반응 경로.
  • 물리학 (전기 화학 전위차, 전기 이중층 및 물질 확산 플럭스): 전극 전위(Overpotential vs. RHE)에 따른 전류 밀도 응답, 다공성 탄소 가스 확산층(GDL)을 통한 $CO_2$ 기체 확산 및 수산화물 이온 확산-영동 수송.
  • 수학 (비선형 표면 커버리지 평형 및 다성분 패러데이 효율 연립방정식 해석): 표면 흡착률 $\theta_{CO}$와 전압에 따른 반응 속도식 $r_{C2H4} = k_{cc} \theta_{CO}^2 \exp(-\alpha F \eta / RT)$, 생성물별 패러데이 효율 $\text{FE}_i = \frac{n_i F \dot{n}_i}{I_{total}}$ 산출.
  • 정보 (무할당 실시간 1,000Hz CCU 전해조 에틸렌 선택도 및 전력 수율 연산 코어): 1ms 실시간 루프 내에서 가비지 컬렉터 없이 스택 메모리 구조체 `ref struct`로 최적 전위점을 0-Alloc 제어하는 C# 온디바이스 엔진.

2. CO2 전기화학 환원 및 C-C 커플링 수식 모델

구리 촉매 표면에서 에틸렌 합성의 속도 결정 단계(RDS)인 두 흡착종 간의 화학적 C-C 커플링 반응 속도 $r_{cc}$는 표면 커버리지 $\theta_{CO}$의 2차식에 비례합니다.

$$r_{cc} = k_{cc} \theta_{CO}^2 \exp\left( -\frac{\Delta G_{act}^{\ddagger}}{R T} \right)$$

인가 전위 $E\,(\text{V vs. RHE})$에 따른 에틸렌 부분 전류 밀도 $j_{C2H4}$와 경쟁 수소 발생 반응(HER) 전류 밀도 $j_{HER}$는 다음과 같습니다.

$$j_{C2H4} = 12 F \cdot r_{cc}, \quad j_{HER} = 2 F k_{HER} (1 - \theta_{CO}) C_{H^+} \exp\left( -\frac{\alpha_{HER} F \eta_{HER}}{R T} \right)$$

국소 알칼리화(Local Alkalinization, $\text{pH} \approx 13$)를 통해 수소 이온 농도 $C_{H^+}$를 억제하고 $\theta_{CO}$를 극대화함으로써 에틸렌 패러데이 선택도 $\text{FE}_{C2H4} = \frac{j_{C2H4}}{j_{total}} \ge 70\%$를 달성하는 최적 운전 윈도우를 실시간 도출합니다.

3. 1,000Hz 0-Alloc C# CCU 에틸렌 합성 반응기 코어

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Education.StemChemistry.CcuCo2Reduction
{
    public readonly record struct ReactorOperatingState(
        float CathodePotentialVsRhe, // 음극 전위 (V vs. RHE, 예: -0.95 V)
        float LocalPhValue,          // 촉매 계면 국소 pH (10.0 ~ 14.0)
        float Co2PartialPressureAtm, // CO2 분압 (atm)
        float TemperatureCelsius     // 운전 온도 (°C)
    );

    public ref struct CcuReactorTelemetry
    {
        public float TotalCurrentDensityMaCm2;
        public float EthyleneCurrentDensityMaCm2;
        public float EthyleneFaradaicEfficiencyPercent;
        public float HydrogenFaradaicEfficiencyPercent;
        public float ProductionRateGramsPerHourPerM2;
        public bool OptimalSelectivityWindow;
    }

    public static class Co2ReductionEngineCore
    {
        private const float FaradayConstantF = 96485.33f;
        private const float EthyleneMolarMass = 28.054f; // g/mol

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void EvaluateReactionStep(
            in ReactorOperatingState state,
            ref CcuReactorTelemetry telemetry)
        {
            float overpotential = MathF.Abs(state.CathodePotentialVsRhe) - 0.17f; // 열역학 평형 전위(0.17V) 대비 과전압
            if (overpotential < 0.05f) overpotential = 0.05f;

            // 1. 랭뮤어 흡착 등온식 표면 *CO 커버리지 추정 (0-Alloc)
            float kAdsorption = 4.5f * state.Co2PartialPressureAtm;
            float thetaCo = Math.Clamp(kAdsorption / (1f + kAdsorption), 0.1f, 0.95f);

            // 2. 국소 pH에 따른 HER 억제율 산출 (알칼리성일수록 H+ 감소로 HER 급감)
            float hConcentrationMolar = MathF.Pow(10f, -state.LocalPhValue);
            float herFactor = Math.Clamp(hConcentrationMolar * 1e12f, 0.05f, 1.0f);

            // 3. 에틸렌 및 수소 부분 전류 밀도 계산 (mA/cm^2)
            float jC2H4 = 85.0f * (thetaCo * thetaCo) * MathF.Exp(0.85f * overpotential * 38.9f * 0.45f);
            float jHer = 30.0f * (1f - thetaCo) * herFactor * MathF.Exp(0.5f * overpotential * 38.9f * 0.4f);
            float jTotal = jC2H4 + jHer + 5.0f; // 기타 C1 부산물 5 mA

            telemetry.TotalCurrentDensityMaCm2 = jTotal;
            telemetry.EthyleneCurrentDensityMaCm2 = jC2H4;

            // 4. 패러데이 효율 및 시간당 에틸렌 생산 속도 (g / h / m^2)
            float feC2H4 = Math.Clamp((jC2H4 / jTotal) * 100f, 0f, 95f);
            telemetry.EthyleneFaradaicEfficiencyPercent = feC2H4;
            telemetry.HydrogenFaradaicEfficiencyPercent = Math.Clamp((jHer / jTotal) * 100f, 0f, 95f);

            // 몰 유량 = (I * FE) / (12 * F), 면적 1 m^2 = 10,000 cm^2
            float currentAmpsPerM2 = jC2H4 * 10.0f;
            float molPerSecM2 = currentAmpsPerM2 / (12f * FaradayConstantF);
            telemetry.ProductionRateGramsPerHourPerM2 = molPerSecM2 * EthyleneMolarMass * 3600f;

            telemetry.OptimalSelectivityWindow = feC2H4 >= 65.0f && jTotal >= 150.0f;
        }
    }
}

4. 교육적 가치 및 국가 탄소중립 화학산업 연계

본 교구는 기후 변화의 원인인 $CO_2$를 단순한 폐기물이 아니라 전 세계 플라스틱 및 정밀 화학 제품의 핵심 원료로 탈바꿈시키는 최첨단 전기화학 촉매 공학을 학생들에게 교육합니다. 학생들은 전위와 계면 pH를 조절하며 귀중한 에틸렌 대신 수소 기체가 낭비되는 원인을 탐구하고, 왜 전 세계 과학자들이 나노 구리 촉매 표면을 원자 단위로 제어하려 하는지를 소프트웨어 실시간 그래프를 통해 명쾌하게 이해하게 됩니다.