AIDT 온디바이스 STEM 전기유체역학(EHD) 테일러 콘(Taylor Cone) 분무 서보 엔진

1. 개요: 노즐 직경의 1/1000 수준 미세 액적을 토출하는 전기유체역학 기술

전통적인 잉크젯 프린팅은 기계적 압력이나 열 팽창을 이용하므로 노즐 오리피스 크기보다 작은 액적을 형성하기 어렵고, 고점도 나노 페이스트 토출 시 노즐 막힘이 빈번합니다. 전기유체역학(Electrohydrodynamics, EHD) 젯 프린팅은 전도성 유체 표면에 강한 정전기장을 인가하여 액적 계면을 원추형인 테일러 콘(Taylor Cone) 형태로 변형시킨 후, 그 정점(Apex)에서 서브마이크론 크기의 초미세 젯을 분사합니다.

반도체 패키징 미세 배선, 바이오 센서 어레이, 차세대 디스플레이 퀀텀닷 패터닝 등 초정밀 적층 제조 분야에서 EHD 젯의 안정적인 제어는 필수적입니다. 본 논문에서는 정전기 맥스웰 응력과 표면장력 사이의 상호작용을 정밀 모델링하고, 전압 펄스와 유량 서보를 10,000Hz 주기로 제어하는 0-Alloc C# 실시간 엔진을 제시합니다.

2. 테일러 콘 형성 조건 및 레일리 전하 한계 수식 모델

유체 표면에 수직으로 작용하는 정전기적 맥스웰 응력(Maxwell stress) \(T_e\)는 다음과 같습니다:

\\[ T_e = \\frac{1}{2} \\epsilon_0 \\epsilon_r E^2 \\]

유체 곡률 \(\\kappa = \\kappa_1 + \\kappa_2\)에 의한 라플라스 모세관 압력 \(\\Delta P_{\\gamma} = \\gamma \\kappa\)와 정전기 응력이 평형을 이룰 때 유체는 반각(half-angle) \(\\theta_0 = 49.3^\\circ\)인 정적 테일러 콘을 형성합니다. 인가 전압이 임계 전압 \(V_c\)를 초과하면 정점에서 하전된 미세 젯이 방출됩니다:

\\[ V_c = \\sqrt{\\frac{\\gamma r_c}{\\epsilon_0}} \\ln\\left(\\frac{4h}{r_c}\\right) \\]

여기서 \(r_c\)는 노즐 반경, \(h\)는 노즐-기판 간 거리, \(\\gamma\)는 표면장력입니다. 액적이 분열하지 않고 안정적으로 전하를 유지할 수 있는 레일리 전하 한계(Rayleigh Limit, \(Q_R\))는 다음과 같습니다:

\\[ Q_R = 8 \\pi \\sqrt{\\epsilon_0 \\gamma R_d^3} \\]

방출 액적의 전하량 \(q\)가 \(Q_R\)에 도달하면 쿨롱 척력이 표면장력을 압도하여 쿨롱 분열(Coulombic Fission)이 일어나 2차 미세 위성 액적(Satellite droplets)이 비산하므로, 서보 시스템은 \(q / Q_R < 0.85\) 상태를 엄격히 유지해야 합니다.

3. 실시간 테일러 콘 EHD 젯 서보 아키텍처 다이어그램

graph TD A[노즐 유량 펌프 서보 Q_in] --> B[유체 펜던트 드롭 계면 형성] C[고전압 펄스 앰프 V_t] --> D[맥스웰 응력 인가 T_e = 0.5 * eps * E^2] B & D --> E[테일러 콘 정점 Apex 왜곡 및 전기장 집중] E --> F{전기장 E > E_crit ?} F -- Yes --> G[서브마이크론 젯 분사 Drop-on-Demand] F -- No --> H[메니스커스 완화 및 전하 축적 대기] G --> I[비행 시간 TOF 및 전하 센서 피드백] I --> J[10,000Hz 0-Alloc C# EHD 젯 서보 제어기] J --> C

4. 고성능 0-Alloc C# 테일러 콘 서보 구현

실시간 EHD 젯의 전압 펄스 폭 변조 및 쿨롱 분열 방지를 위한 0-Alloc C# 엔진입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.EhdJetting
{
    public readonly record struct EhdParameters(
        double SurfaceTensionGamma,    // N/m (예: 0.0728 물 기준, 에탄올 0.022)
        double PermittivityEpsR,       // 비유전율
        double NozzleRadiusRc,         // m (예: 25 um)
        double SubstrateDistanceH,     // m (예: 500 um)
        double FlowRateQ               // m^3/s
    );

    public ref struct EhdServoState
    {
        public double AppliedVoltage;
        public double ConeApexField;
        public double DropletRadius;
        public double DropletCharge;
        public double RayleighLimit;
        public double JetVelocity;
        public bool IsJetting;
    }

    public static class EhdTaylorConeServo
    {
        private const double Eps0 = 8.8541878128e-12; // F/m

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void StepServo(
            in EhdParameters p,
            double commandVoltage,
            double dt,
            ref EhdServoState state)
        {
            // 1. 임계 전압 V_c 계산
            double lnTerm = Math.Log(4.0 * p.SubstrateDistanceH / p.NozzleRadiusRc);
            double vCritical = Math.Sqrt((p.SurfaceTensionGamma * p.NozzleRadiusRc) / Eps0) * lnTerm;

            state.AppliedVoltage = commandVoltage;

            // 2. 테일러 콘 정점 전기장 집중도 모델링
            double effectiveGap = Math.Max(1e-6, p.SubstrateDistanceH);
            double baseField = commandVoltage / effectiveGap;
            double enhancementFactor = 1.0 + (lnTerm * 1.8);
            state.ConeApexField = baseField * enhancementFactor;

            // 3. 젯팅 발생 판정
            if (commandVoltage >= vCritical)
            {
                state.IsJetting = true;
                // 방출 액적 반경: 스케일링 법칙 R_d ~ (rho * eps0 * Q^3 / gamma^2)^(1/6)
                double rEst = Math.Cbrt(p.FlowRateQ / (Math.PI * 400.0)); // 근사 스케일링
                state.DropletRadius = Math.Clamp(rEst, 50e-9, p.NozzleRadiusRc);

                // 레일리 한계 계산 Q_R = 8 * pi * sqrt(eps0 * gamma * R^3)
                state.RayleighLimit = 8.0 * Math.PI * Math.Sqrt(Eps0 * p.SurfaceTensionGamma * Math.Pow(state.DropletRadius, 3));

                // 축적 전하량 추정
                state.DropletCharge = state.RayleighLimit * 0.75; // 레일리 한계 75% 수준 유지

                // 젯 분사 속도: v_jet = (2 * q * V / m_d)^0.5
                double mass = (4.0 / 3.0) * Math.PI * Math.Pow(state.DropletRadius, 3) * 1000.0;
                state.JetVelocity = Math.Sqrt((2.0 * state.DropletCharge * commandVoltage) / mass);
            }
            else
            {
                state.IsJetting = false;
                state.DropletRadius = 0.0;
                state.DropletCharge = 0.0;
                state.RayleighLimit = 0.0;
                state.JetVelocity = 0.0;
            }
        }
    }
}

5. 교육 연수 및 방위 산업 적용 지침

방통대 직무연수 및 마이크로 적층 전자전 응용

물리/화학 교과 지도: 정전기 유도와 표면장력 사이의 계면 불안정성(Taylor instability), 레일리 전하 한계를 시각적 시뮬레이션 교구로 구현.
국방 MUM-T 전자전 패키징: 전파 흡수체(RAM) 마이크로 패턴 인쇄 및 유연 기판 초소형 RF 안테나의 현장 주문형(On-demand) 3D 수리 제작에 적용.