AIDT 온디바이스 STEM 고엔트로피 합금(HEA) 레이저 분말 적층(SLM) 멜트풀 서보 엔진

1. 개요: 5개 이상의 주원소가 고용체를 이루는 신개념 금속의 초고속 레이저 3D 적층

전통적인 금속 합금(철강, 알루미늄, 티타늄)은 한 가지 주요 기주 원소에 소량의 부원소를 첨가하는 방식이었습니다. 고엔트로피 합금(High-Entropy Alloy, HEA, 예: FeCoNiCrMn 칸토어 합금)은 5개 이상의 원소가 등원자비(Equiatomic)로 혼합되어 높은 배열 혼합 엔트로피(\(\Delta S_{\text{mix}} > 1.5 R\)) 효과로 인해 복잡한 금속간 화합물 대신 단순 단상 고용체(FCC/BCC)를 형성하여, 극저온 충격인성과 고온 강도를 동시에 극대화합니다.

선택적 레이저 용융(Selective Laser Melting, SLM) 공정으로 HEA를 3D 적층할 때, \(10^6\,\text{K/s}\)에 달하는 급냉 응고(Rapid solidification)가 발생합니다. 레이저 에너지 밀도가 과도하면 금속 증기압에 의해 키홀(Keyhole) 기공이 형성되고, 부족하면 미용융(Lack of Fusion) 결함이 생깁니다. 본 백서에서는 로젠탈(Rosenthal) 3차원 이동 열원 모델을 기반으로 멜트풀(Melt Pool) 형상을 실시간 예측하고 레이저 출력을 변조하는 10,000Hz 0-Alloc C# 서보를 구현합니다.

2. 깁스 혼합 자유에너지 및 로젠탈 멜트풀 수식 모델

다원소 HEA 계의 깁스 혼합 자유에너지 \(\Delta G_{\text{mix}}\)는 혼합 엔탈피 \(\Delta H_{\text{mix}}\)와 이상 혼합 엔트로피 \(\Delta S_{\text{mix}}\)의 경쟁으로 결정됩니다:

\\[ \Delta G_{\text{mix}} = \Delta H_{\text{mix}} - T \Delta S_{\text{mix}}, \quad \Delta S_{\text{mix}} = -R \sum_{i=1}^n c_i \ln c_i \\]

단상 고용체 안정성을 보장하는 \(\Omega\) 파라미터는 다음과 같습니다:

\\[ \Omega = \frac{T_m \Delta S_{\text{mix}}}{|\Delta H_{\text{mix}}|} \ge 1.1, \quad \delta = \sqrt{\sum_{i=1}^n c_i \left( 1 - \frac{r_i}{\bar{r}} \right)^2} \le 6.6\% \\]

이동 속도 \(v_s\), 레이저 흡수율 \(\eta\), 출력 \(P\)를 가진 반무한 고체 표면의 준정상 로젠탈(Rosenthal) 온도장 방정식은 다음과 같습니다:

\\[ T(x, y, z) - T_0 = \frac{\eta P}{2 \pi k R} \exp\left( -\frac{v_s (R + x)}{2 \alpha_{\text{th}}} \right), \quad R = \sqrt{x^2 + y^2 + z^2} \\]

여기서 \(k\)는 열전도도, \(\alpha_{\text{th}} = \frac{k}{\rho C_p}\)는 열확산율입니다. 용융 온도 \(T_m\) 등온면이 정의하는 멜트풀 깊이 \(d_{\text{mp}}\)와 폭 \(w_{\text{mp}}\)의 종횡비 \(\text{AR} = d_{\text{mp}} / w_{\text{mp}}\)가 \(0.8\)을 초과하면 키홀 불안정 붕괴(Keyhole collapse)로 인한 가스 기공이 갇히게 되므로 능동 파워 제어가 필수적입니다.

3. 실시간 SLM 멜트풀 동축 광학 서보 아키텍처 다이어그램

graph TD A[Yb 파이버 레이저 1070nm 파워 P_laser] --> B[검류계 갈바노 미러 스캐너 속도 v_s] B --> C[HEA 분말 베드 멜트풀 형성 T > T_liquidus] C --> D[고속 광다이오드 동축 멜트풀 열복사 모니터링] D --> E[2파장 고온계 멜트풀 온도 T_mp & 면적 A_mp 계측] E --> F[로젠탈 역산 멜트풀 종횡비 AR = d / w] F --> G{AR > 0.75 키홀 붕괴 위험?} G -- Yes --> H[레이저 파워 P_laser 급속 감발] G -- No --> I[표준 전도 용융 출력 유지] H & I --> J[10,000Hz 0-Alloc C# 에너지 밀도 보정기] J --> A

4. 고성능 0-Alloc C# HEA 멜트풀 형상 제어 서보 구현

스캔 속도와 동축 열복사 센서값으로부터 멜트풀 치수를 역산하고 레이저 파워를 보정하는 0-Alloc C# 클래스입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.MetalAdditiveManufacturing
{
    public readonly record struct HeaSlmParameters(
        double LiquidusTempK,     // Tm (예: FeCoNiCrMn ~ 1600 K)
        double DensityRho,        // kg/m^3 (8000)
        double SpecificHeatCp,    // J/(kg K) (500)
        double ThermalCondK,      // W/(m K) (15)
        double LaserAbsorptionEta // 0.35
    );

    public ref struct MeltPoolServoState
    {
        public double CommandLaserPowerW;
        public double MeltPoolWidthUm;
        public double MeltPoolDepthUm;
        public double AspectRatio;
        public double LinearEnergyDensityJmm;
        public bool IsKeyholeRisk;
    }

    public static class HighEntropyAlloySlmServo
    {
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void StepServo(
            in HeaSlmParameters p,
            double nominalPowerW,
            double scanSpeedMps,    // m/s (예: 1.0 m/s = 1000 mm/s)
            double preheatTempK,    // 기판 예열 온도 (예: 473 K)
            ref MeltPoolServoState state)
        {
            state.CommandLaserPowerW = nominalPowerW;

            // 1. 선형 에너지 밀도 LED = P / v (J/mm)
            state.LinearEnergyDensityJmm = state.CommandLaserPowerW / (scanSpeedMps * 1000.0);

            // 2. 열확산율 alpha = k / (rho * Cp)
            double alphaTh = p.ThermalCondK / (p.DensityRho * p.SpecificHeatCp);

            // 3. 로젠탈 기반 무차원 용융 치수 스케일링
            // R_crit ~ eta * P / (2 * pi * k * (Tm - T0))
            double deltaT = Math.Max(100.0, p.LiquidusTempK - preheatTempK);
            double rCritM = (p.LaserAbsorptionEta * state.CommandLaserPowerW) / (2.0 * Math.PI * p.ThermalCondK * deltaT);

            // 페클레 수 Pe = v * r_crit / (2 * alpha)
            double pe = (scanSpeedMps * rCritM) / (2.0 * alphaTh);
            double elongationFactor = Math.Sqrt(1.0 + pe);

            // 멜트풀 폭 및 깊이 근사 (um 단위)
            state.MeltPoolWidthUm = (rCritM / elongationFactor) * 2.0 * 1e6;
            state.MeltPoolDepthUm = (rCritM / Math.Pow(elongationFactor, 1.2)) * 1e6;

            // 4. 종횡비 계산
            state.AspectRatio = state.MeltPoolDepthUm / Math.Max(10.0, state.MeltPoolWidthUm);

            // 5. 키홀 모드 위험성 판정 (AR > 0.75)
            state.IsKeyholeRisk = state.AspectRatio > 0.75;

            // 6. 능동 파워 피드백 보정: 키홀 진입 시 파워 감발
            if (state.IsKeyholeRisk)
            {
                double excessRatio = state.AspectRatio / 0.75;
                state.CommandLaserPowerW = nominalPowerW / excessRatio;
            }
        }
    }
}

5. 교육 연수 및 극초음속 비행체/우주 터빈 국방 적용 지침

방통대 직무연수 및 국방 초고온 내열 부품 제조

물리/재료 교과 지도: 깁스 자유에너지와 배열 엔트로피 계산, 금속의 3차원 열전도 및 레이저-물질 상호작용(전도 용융 대 키홀 용융)을 가상 3D 프린팅으로 교육.
국방 극초음속 순항미사일 노즐 및 로켓 연소실: 마하 5 이상 비행 시 발생하는 1,500도 이상의 공력 가열을 견디는 고엔트로피 합금 터빈 블레이드의 내부 냉각 유로 일체형 3D 프린팅.