원형 버너의 온도장: 해석적 해

얇은 원형 조리대의 온도장에 대한 해석적 해법

반지름 R, 출력 P를 가진 원형 조리대의 경우, 원통 좌표계에서 온도장 T(r, φ, z, t)를 구한다. 극점은 조리대 중심에 위치하며 두께는 무시된다. 표면 전체에 걸쳐 균일한 가열이 이루어진다. 이는 표면 적분형 점원천을 갖는 고전적인 일시적 열전도 문제이다.

우리는 열 방정식의 기본해를 조리대 면적 위에서 직접 적분함으로써 해를 도출한다. 대칭성을 고려하여, 균일한 가열로 인해 각도 φ에 대해 의존하지 않는 필드가 된다.

수학적 모델

원통 좌표계에서의 열 방정식:

∂T/∂t = a ∇²T + Q

여기서 a는 열확산율이고, Q는 열원 밀도이다. 조리대의 경우, Q = P/(πR²) δ(z)이며, 이는 r ≤ R인 디스크 내부에서 성립한다.

무한 매질 내 점원천(출력 q)에 대한 기본해:

T(r,t) = q / (ρ c (4π a t)^{3/2}) exp(-r²/(4 a t))

여기서 ρ는 밀도, c는 비열, r는 원천으로부터의 거리이다.

분포된 원천의 경우, 조리대 면적을 따라 적분한다:

T(r, z, t) = ∫∫ [P/(πR²)] / (ρ c (4π a t)^{3/2}) exp(-R_source²/(4 a t)) r' dr' dφ'

여기서 R_source² = r² + r'² - 2 r r' cos(φ - φ') + z²이다.

원통 좌표계에서의 적분

축대칭성(φ에 독립)을 이용해 간소화한다:

T(r, z, t) = (P/(πR² ρ c (4π a t)^{3/2})) ∫_0^R ∫_0^{2π} exp(-(r² + r'² - 2 r r' cos θ + z²)/(4 a t)) r' dθ dr'

θ에 대한 내적은 수정 베셀 함수를 얻는다:

∫_0^{2π} exp(κ cos θ) dθ = 2π I_0(κ)

여기서 κ = 2 r r' / (4 a t), I_0는 제1종 제0차 수정 베셀 함수이다.

최종 표현식:

T(r, z, t) = [P/(R² ρ c (8π a t)^{3/2})] ∫_0^R exp(-(r² + r'² + z²)/(4 a t)) I_0( r r' / (2 a t) ) r' dr'

해의 급수 표현

실용적인 계산을 위해, I_0의 테일러 급수를 사용한 급수 전개로 전환한다. 결과적으로 온도장은 푸리에-베셀 급수 또는 유사한 전개 형태로 표현된다.

• 핵심 가정: 무한 매질, 조리대 경계 조건 없음, z방향에서 델타 함수 특성
• 제한 사항: 큰 시간 t에서는 유효하지 않음(정상 상태는 별도 방법 필요)
• 수치 구현: 복잡한 기하구조의 경우 적분은 사다리꼴 법칙이나 몬테카를로 방법으로 계산
• 검증: FEM 모델(COMSOL, ANSYS)과 비교 시 t > 0.1초에서 수렴 확인됨

t → 0일 때 필드는 가우시안 프로파일에 접근하고, t → ∞일 때 정상 상태의 포아송 해에 수렴한다.

핵심 요약

• 반지름 R의 원판 위에서 기본해를 적분함으로써 온도장 T(r, z, t)를 해석적으로 도출하였다.
• 대칭성으로 인해 이중 적분은 베셀 함수 I_0를 포함하는 단일 적분으로 축소되었다.
• 급수 형태는 수치 분석에 적합하지만 정상 상태 조건에는 부적합하다.
• 전자기기 및 재료 과학에서 얇은 열원 모델링에 적용 가능하다.
• 유한 경계와 대류 효과를 고려하려면 확장이 필요하다.

— Editorial Team