2022 개정 교육과정 연계 온디바이스 극지해양·빙하동역학 STEM 시리즈 #72

빙상 착지선(Grounding Line) 동역학 및 빙붕 하부 융해 서보 엔진

해양성 빙상 불안정성(MISI) 착지선 후퇴율 추정 및 0-Alloc C# 자율무인잠수정(AUV) 해류-융해 연동 제어 루프

ZeliDesk 극지해양학·빙하물리팀 남극스웨이츠탐사·해수면변동 분과

2026년 9월 12일

1. 부력 정역학 평형 및 셸프-빙상 착지선 질량 유속 수리 모델

빙상이 바다에 떠오르는 착지선에서의 정역학적 부력 평형 조건은 얼음 두께 $H$와 기반암 수심 $b(x)$ 사이에서 성립합니다:

여기서 $\rho_i \approx 917\text{ kg/m}^3$은 빙하 밀도, $\rho_w \approx 1028\text{ kg/m}^3$은 남극 해수 밀도입니다. 쉔펠트(Schoof) 경계층 이론에 따른 착지선 통과 빙하 질량 플럭스 $q_g$는 빙하 두께의 약 4.5제곱에 비례합니다:

빙붕 하부 공동의 난류 열 교환에 의한 국소 기저 융해율 $\dot{m}_b$는 열 구동력(Thermal Driving $\Delta T = T_{ocean} - T_{freeze}$)과 마찰 속도 $u_*$의 곱으로 결정되는 3방정식 경계층 모델을 따릅니다:

2. 0-Alloc C# 스택 기반 극지 AUV 착지선 융해 추종 서보

빙붕 하부 자율무인잠수정(AUV)의 음향 거리계 및 수온 센서 데이터로부터 실시간 기저 융해율을 산출하고 추진체 안전 심도를 유지하는 ref struct입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.PolarOceanography
{
    public readonly ref struct CavityOceanTelemetry
    {
        public readonly double AmbientWaterTempC;
        public readonly double InSituSalinityPsu;
        public readonly double AcousticIceCeilingDistanceM;
        public readonly double BathymetryDepthM;
        public readonly double CurrentVelocityMPerSec;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public CavityOceanTelemetry(double tempC, double salPsu, double iceDistM, double bathyM, double velMps)
        {
            AmbientWaterTempC = tempC;
            InSituSalinityPsu = salPsu;
            AcousticIceCeilingDistanceM = iceDistM;
            BathymetryDepthM = bathyM;
            CurrentVelocityMPerSec = velMps;
        }
    }

    public readonly ref struct GroundingLineAnalysisResult
    {
        public readonly double BasalMeltRateMPerYear;
        public readonly double DistanceToGroundingLineKm;
        public readonly double TargetAuvDepthM;
        public readonly bool IsMisiWarningActive;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public GroundingLineAnalysisResult(double meltRate, double distGl, double targetDepth, bool misiWarning)
        {
            BasalMeltRateMPerYear = meltRate;
            DistanceToGroundingLineKm = distGl;
            TargetAuvDepthM = targetDepth;
            IsMisiWarningActive = misiWarning;
        }
    }

    public ref struct GroundingLineMeltServo
    {
        private const double RhoW = 1028.0;
        private const double RhoI = 917.0;
        private const double Lf = 334000.0; // 잠열 J/kg
        private const double Cp = 3986.0;   // 해수 비열 J/(kg*K)
        private const double GammaT = 1.1e-4; // 무차원 열전달 계수

        private readonly double _bedSlopeGradient;
        private double _prevMeltRate;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public GroundingLineMeltServo(double bedSlope)
        {
            _bedSlopeGradient = bedSlope;
            _prevMeltRate = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public GroundingLineAnalysisResult Step(in CavityOceanTelemetry telemetry)
        {
            // 1. 현장 수심에서의 어는점 계산 (압력 및 염분 의존: T_f = -0.0575 * S + 0.0901 - 7.61e-4 * Depth)
            double tFreeze = -0.0575 * telemetry.InSituSalinityPsu + 0.0901 - 7.61e-4 * telemetry.BathymetryDepthM;
            double thermalDriving = Math.Max(0.0, telemetry.AmbientWaterTempC - tFreeze);

            // 2. 난류 마찰 속도 u_* = sqrt(Cd) * U (Cd ~ 0.0025)
            double uStar = Math.Sqrt(0.0025) * Math.Abs(telemetry.CurrentVelocityMPerSec);

            // 3. 3방정식 기반 국소 기저 융해율 (m/s -> m/year 변환: * 31,536,000)
            double meltRateMps = (RhoW * Cp * GammaT * uStar * thermalDriving) / (RhoI * Lf);
            double meltRateMpy = meltRateMps * 31536000.0;
            _prevMeltRate = meltRateMpy;

            // 4. 착지선 잔여 거리 추정 (빙상 천장과 해저 바닥 사이 간극이 좁아지는 기하학적 수렴)
            double waterColumnThickness = telemetry.AcousticIceCeilingDistanceM;
            double distToGroundingLineKm = Math.Max(0.1, waterColumnThickness / Math.Max(1e-4, Math.Abs(_bedSlopeGradient)) / 1000.0);

            // 5. MISI (역경사 침식 가속) 경보 조건 (융해율 > 40 m/yr 이고 착지선 접근 2km 이내)
            bool misiWarning = meltRateMpy > 40.0 && distToGroundingLineKm < 2.0;

            // 6. AUV 충돌 회피 목표 수심 (천장과 바닥의 정중앙 유지)
            double targetDepth = telemetry.BathymetryDepthM - (waterColumnThickness * 0.5);

            return new GroundingLineAnalysisResult(meltRateMpy, distToGroundingLineKm, targetDepth, misiWarning);
        }
    }
}

3. 교육적 가치 및 중고교 실험 시뮬레이션 적용