일상 속 인체 보행, 파도, 미세 진동과 같은 버려지는 기계적 에너지를 전기에너지로 변환하는 마찰대전 나노발전기(Triboelectric Nanogenerator, TENG)는 배터리 교체가 불가능한 웨어러블 헬스케어 디바이스와 자가발전 스마트 센서 네트워크의 혁신 동력입니다. 두 유전체 표면의 접촉 대전(Contact Electrification)과 정전 유도(Electrostatic Induction) 현상, 그리고 맥스웰 변위 전류(Maxwell's Displacement Current) 메커니즘을 1,000Hz 주기로 0-Alloc C# 수치해석하는 고성능 물리 에너지 하베스팅 시뮬레이터를 공개합니다.

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

  • 물리학 (정전기 유도, 전하 이동, 맥스웰 변위 전류): 이종 물질 간 마찰 시 페르미 준위 차이로 인한 표면 전하 밀도 $\sigma$ 형성, 시간에 따라 변하는 전기장에 의해 유도되는 변위 전류 밀도 $J_D = \varepsilon_0 \frac{\partial E}{\partial t}$.
  • 화학/신소재 (마찰 대전 서열, 불소계 고분자 PTFE 및 표면 나노패터닝): 전자 공여체(금속 전극 Al/Cu)와 강한 전자 수용체(PTFE, PDMS) 계면 분자의 접촉 화학 반응 및 나노 피라미드 구조에 의한 접촉 표면적 극대화.
  • 수학 (시간 가변 커패시터 미분방정식 및 전하-전압 $V-Q-x$ 관계식): 전극 간격 $x(t)$의 주기적 변위에 따른 등가 정전용량 $C(t) = \frac{\varepsilon_0 A}{d_0 + x(t)}$의 1차 비선형 미분방정식 수치해법.
  • 정보 (무할당 실시간 1,000Hz 펄스 피크 검출 및 정류 축전 에너지 적분 코어): 불규칙한 기계적 외력 충격에 의해 생성되는 수백 볼트 임펄스 전압을 1ms 내에 정류하고 커패시터 축전량을 0-Alloc 계산하는 C# 온디바이스 엔진.

2. TENG 맥스웰 변위 전류 및 $V-Q-x$ 수식 모델

두 대전체 판 사이의 간격 $x(t)$, 유전체 두께 $d$, 상대 유전율 $\varepsilon_r$, 마찰 표면 전하 밀도 $\sigma$, 전극 면적 $S$에 대해 TENG 개방 회로 전압(Open-Circuit Voltage) $V_{OC}$와 단락 회로 전하량(Short-Circuit Charge) $Q_{SC}$는 다음과 같이 정의됩니다:

$$V_{OC}(t) = \frac{\sigma \cdot x(t)}{\varepsilon_0}, \quad Q_{SC}(t) = \frac{S \sigma x(t)}{d_0 + x(t)} \quad \left( d_0 = \frac{d}{\varepsilon_r} \right)$$

외부 부하 저항 $R$이 연결되었을 때, 전하 이동량 $Q(t)$에 대한 동적 지배방정식은 키르히호프 법칙에 의해 다음과 같이 도출됩니다:

$$R \frac{dQ}{dt} = - \frac{1}{C(x)} Q + \frac{\sigma x(t)}{\varepsilon_0} = - \left( \frac{d_0 + x(t)}{\varepsilon_0 S} \right) Q + \frac{\sigma x(t)}{\varepsilon_0}$$

전기장 변화에 따른 맥스웰 변위 전류 밀도 $J_D$는 발전기 내부를 가로지르는 비도전성 변위 전류를 정량화하며, 피크 순간 전력 밀도 $P = I^2 R / S$를 최대화하는 최적 정합 임피던스 $R_{opt} \approx \frac{1}{\omega C_{avg}}$가 존재합니다.

3. 1,000Hz 0-Alloc C# TENG 에너지 수확 해석 코어

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Education.StemEnergy.Triboelectric
{
    public readonly record struct TengDeviceConfig(
        float ContactAreaCm2,        // 접촉 면적 (cm^2, 예: 25.0)
        float DielectricThicknessUm, // 유전체 필름 두께 (um, 예: 50.0)
        float RelativePermittivity,  // 상대 유전율 (PTFE: 2.1)
        float SurfaceChargeDensity,  // 표면 전하 밀도 (uC/m^2, 예: 65.0)
        float LoadResistanceMegaOhm  // 외부 부하 저항 (MOhm, 예: 10.0)
    );

    public ref struct TengHarvestingTelemetry
    {
        public float GapDistanceMm;
        public float OpenCircuitVoltageV;
        public float OutputCurrentMicroAmp;
        public float InstantaneousPowerMilliWatt;
        public float AccumulatedEnergyMicroJoule;
        public float DisplacementCurrentDensity;
    }

    public static class TengEnergyEngineCore
    {
        private const float Eps0 = 8.8541878e-12f; // F/m

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void SimulateStep(
            in TengDeviceConfig config,
            float gapDistanceMeters,
            float gapVelocityMps,
            ref float inOutTransferredChargeCoulomb,
            ref float inOutAccumulatedEnergyJ,
            float dtSeconds,
            ref TengHarvestingTelemetry telemetry)
        {
            // 1. 디바이스 파라미터 변환 (0-Alloc)
            float areaM2 = config.ContactAreaCm2 * 1e-4f;
            float sigma = config.SurfaceChargeDensity * 1e-6f; // C/m^2
            float d0 = (config.DielectricThicknessUm * 1e-6f) / MathF.Max(1.0f, config.RelativePermittivity);
            float x = MathF.Max(1e-7f, gapDistanceMeters);

            // 2. 개방 전압 V_oc = sigma * x / eps0
            float voc = (sigma * x) / Eps0;
            telemetry.GapDistanceMm = x * 1000.0f;
            telemetry.OpenCircuitVoltageV = voc;

            // 3. 등가 커패시턴스 C(t) = eps0 * S / (d0 + x)
            float cap = (Eps0 * areaM2) / (d0 + x);

            // 4. 동적 지배방정식 dQ/dt = (Voc - Q / C) / R
            float loadOhm = config.LoadResistanceMegaOhm * 1e6f;
            float vCap = inOutTransferredChargeCoulomb / MathF.Max(1e-15f, cap);
            float current = (voc - vCap) / MathF.Max(1.0f, loadOhm);

            // 오일러 수치 적분 (0-Alloc)
            inOutTransferredChargeCoulomb += current * dtSeconds;
            telemetry.OutputCurrentMicroAmp = current * 1e6f;

            // 5. 순간 전력 P = I^2 * R (mW) 및 누적 에너지 적분 (uJ)
            float powerWatts = current * current * loadOhm;
            telemetry.InstantaneousPowerMilliWatt = powerWatts * 1000.0f;

            inOutAccumulatedEnergyJ += powerWatts * dtSeconds;
            telemetry.AccumulatedEnergyMicroJoule = inOutAccumulatedEnergyJ * 1e6f;

            // 맥스웰 변위 전류 밀도 J_D = sigma * v / (d0 + x) 근사
            float jd = (sigma * d0 * gapVelocityMps) / MathF.Max(1e-12f, (d0 + x) * (d0 + x));
            telemetry.DisplacementCurrentDensity = jd;
        }
    }
}

4. 교육 및 연구 현장 활용 제언

본 시뮬레이터는 한국기계연구원(KIMM) 나노융합장비 연구단 및 방통대 우호성 교수 연구실의 차세대 친환경 자가발전 에너지 교원 연수에서 기계적 운동 에너지를 마찰 전기로 변환하는 과정을 정량적으로 학습하는 실험 교구로 활용될 수 있습니다. 학습자는 전극 간격 변위 속도와 유전체 상대 유전율에 따른 변위 전류 생성 효율을 1,000Hz 실시간으로 시뮬레이션하며 에너지 하베스팅의 물리적 원리를 완벽히 체득할 수 있습니다.