입사각과 측면 전도 효과를 고려한 전자기 피부 가열의 고차 비대칭 해법

초록

본 논문은 전자기 빔이 피부에 입사할 때 발생하는 온도 분포를, 입사각과 측면 열전도 효과를 동시에 포함한 고차 비대칭 해석으로 다룬다. 깊이‑횡방향 비율 ε가 작을 때 온도장을 ε의 거듭제곱으로 전개하고, ε⁰, ε¹, ε² 항을 각각 구해 닫힌 형태의 해를 제시한다. ε≈0.1 수준에서도 두 번째 차수(ε²) 항이 온도 예측에 크게 기여함을 수치 검증한다.

상세 분석

이 연구는 밀리미터파(MMW) 대역의 전자기 빔이 인체 피부에 미치는 열효과를 정량화하기 위해, 비대칭 다중 스케일 접근법을 채택하였다. 핵심 가정은 전자기 파가 침투하는 깊이(≈1 mm)와 빔 횡단면 크기(수 cm) 사이의 스케일 차이가 매우 크다는 점이며, 이를 비율 ε = zₛ/rₛ(깊이/횡방향 길이)로 정의한다. ε가 작을수록 깊이 방향 열전도는 지배적이며, 횡방향 전도는 ε² 차수에만 나타난다.

입사각 θ₁은 두 가지 물리적 메커니즘을 통해 온도장에 영향을 미친다. 첫째, 빔이 표면에 투사되는 면적이 1/cosθ₁ 만큼 늘어나면서 전력 밀도가 cosθ₁ 배로 감소한다. 둘째, 굴절각 θ₂(스넬 법칙에 의해 결정) 때문에 파가 깊이 z를 진행할 때 실제 전파 거리는 z cosθ₂이며, 이는 흡수율을 1/ cosθ₂ 로 늘린다. 이 두 효과는 ε⁰ 및 ε¹ 항에 모두 포함된다.

또한, 굴절에 따라 빔이 깊이와 동시에 x‑방향으로 이동한다(이동량 = z tanθ₂). 이 “깊이‑의존 횡방향 이동”은 ε 차수와 결합되어 ε¹ 항에 나타나지만, 실제 온도분포에 미치는 영향은 ε² 항에서 측면 열전도와 결합될 때 비로소 정량화된다.

수학적으로 저자들은 온도 T(x,y,z,t)를
T = T⁽⁰⁾ + ε T⁽¹⁾ + ε² T⁽²⁾ + …
와 같이 전개하고, 각 차수마다 초기·경계값 문제(IBVP)를 도출한다. T⁽⁰⁾와 T⁽¹⁾는 1‑차원 열전도 방정식(횡방향 라플라시안이 사라짐)으로 축소되며, 해는 각각 W⁽⁰⁾(z,t), W⁽¹⁾(z,t) 형태의 함수로 표현된다. ε² 차수는 두 개의 서브항 A, B 로 나뉘어, A는 횡방향 라플라시안(ε²∇⊥²T⁽⁰⁾) 효과, B는 입사각에 의한 비선형 이동(θ₂·∂ₓT⁽¹⁾) 효과를 포함한다. 저자들은 변수를 분리하고 라플라시안 연산자를 적용해 각각의 ODE/PDE를 풀어, 폐쇄형 해를 얻는다.

수치 검증에서는 중앙 차분법으로 3‑차원 열 방정식을 고해상도 격자에서 직접 풀어 “정확한” 해를 만든 뒤, ε = 0.1, θ₁ = 30° 등 네 가지 시나리오(횡방향 전도 유무·입사각 유무)를 비교한다. 결과는 ε⁰ 해만 사용할 경우 오차가 20 % 이상 발생하지만, ε¹·ε² 항을 포함하면 오차가 2 % 이하로 감소한다. 특히, 횡방향 전도가 존재할 때 ε² 항이 온도 피크와 전파 속도에 미치는 영향이 두드러진다.

이 논문의 주요 공헌은 (1) 입사각에 따른 전력 밀도 감소와 굴절에 의한 깊이‑횡방향 결합을 정확히 모델링한 점, (2) ε² 차수까지 전개해 횡방향 전도 효과를 정량화한 점, (3) 폐쇄형 해를 제공해 실시간 온도 예측이나 역문제(표면 온도 기반 내부 온도 추정)에 바로 적용 가능하도록 한 점이다. 한계로는 피부 비균질성(다중 층 구조)과 표면 열손실을 무시했으며, ε가 0.2 이상일 경우 고차 항이 더 필요할 수 있다는 점을 들 수 있다. 향후 연구에서는 다층 모델, 비선형 열전도, 그리고 실제 인체 실험 데이터를 통한 검증이 요구된다.