Fourier 열전도와 열 수지

막대의 온도가 균일해져도 열은 사라지지 않는다

따뜻한 부분과 차가운 부분이 번갈아 있는 원형 막대를 생각해봅시다. 주변과는 열을 주고받지 않습니다. 열은 따뜻한 부분에서 차가운 이웃으로 이동하고 시간이 지나면 온도가 거의 균일해집니다. 열이 사라진 걸까요? 닫힌 주기 모델에서는 전체 현열 에너지가 일정한 채 온도 차이가 줄어듭니다. 단지 ‘에너지’라고 붙인 곡선이 내려간다고 해서 실제 열에너지가 없어졌다고 해석하면 안 됩니다. 온도 불균일성의 수학적 측도와 구분해야 합니다.

이 단원은 공간 Fourier 모드 하나를 따라가며 이 차이를 직접 계산합니다. 열유속과 국소 수지를 유도하고 임의의 관측 창과 한 주기 전체를 비교합니다. 매끄러운 애니메이션이나 작은 PDE 잔차 하나만으로 모델이 어떤 물리량을 보존하는지는 알 수 없습니다.

선수 지식, 가정과 물리 단위

기울기, 적분, 보존 수지와 지수 감쇠를 복습합니다. 정지한 일차원 물질에서 체적 열용량 C는 양의 상수이고 열전도도 kappa는 음이 아닌 상수입니다. 내부 열원, 상전이와 복사는 없습니다. 일정한 실제 단면적을 지정합니다. 온도는K의 절대 온도입니다.1K의 온도 차이는 섭씨1도의 차이와 크기가 같지만 절대 온도와 섭씨 눈금의 수치는 다릅니다.

물리량	단위	의미
T	K	절대 온도
C	J/(m³ K)	체적당 단위 온도 변화에 필요한 에너지
kappa	W/(m K)	열전도도
q	W/m²	단위 면적을 통과하는 방향 있는 열유속
단면적	m²	막대의 실제 단면적

$$C=\rho c_p,\qquad [C]=\mathrm{J/(m^3K)},\quad[\kappa]=\mathrm{W/(mK)},\quad[q]=\mathrm{W/m^2}.$$

손계산의 물성 수치는 교육용 척도이며 특정 금속의 측정값이 아닙니다. 일정한 열용량은 명시한 온도 범위에서 사용하는 구성 근사입니다.

방향 있는 열유속에서 열방정식 유도하기

Fourier 구성법칙은 온도 기울기가 내려가는 쪽으로 열이 흐른다고 정합니다. 이 법칙을 국소 에너지 보존과 결합합니다. q가 양수라는 것은 증가하는 x 방향으로 흐른다는 뜻이지 모든 관측 창으로 열이 들어온다는 뜻이 아닙니다. 고정 창의 왼쪽 끝 바깥 법선은−1, 오른쪽 끝은+1입니다. 국소 열유속을 순수출 전력으로 바꿀 때 이 방향을 놓치면 부호가 뒤집힙니다.

$$q=-\kappa T_x,\qquad C T_t+q_x=0,\qquad T_t=\alpha T_{xx},\qquad\alpha=\frac{\kappa}{C}.$$

열확산율의 단위는m²/s입니다. 시간 척도를 정하지만 열전도도와 같은 물리량은 아닙니다. 두 물질의 열용량과 열전도도를 같은 배율로 바꾸면 열확산율과 온도 진행은 같아도 주어진 기울기에서 운반하는 열의 전력은 달라집니다.

양의 온도를 갖는 Fourier 모드 구성하기

양의 균일 평균 위에 부호 있는 코사인 섭동을 놓습니다. 양의 파수와 위상은 따뜻한 곳과 차가운 곳의 간격과 위치를 정합니다. 대비의 부호를 바꾸면 두 영역이 교환됩니다. 상대 대비의 절댓값을0.2 이하로 제한하므로 온도는 양수입니다. 주기 막대 모델이며 양끝을 특정 온도에 고정한 유한 막대의 해가 아닙니다.

$$T(x,t)=\overline T+B(t)\cos(kx+\phi),\quad B(t)=\overline T a e^{-\alpha k^2t},\qquad |a|\le0.2.$$

한 주기의 길이는2pi/k입니다. 코사인의 한 주기 적분이 영이므로 평균 온도는 일정합니다. 임의의 더 짧은 구간에는 온도 섭동의 적분이 남을 수 있고 주변 구간과 열을 주고받습니다. 주기 전체의 보존식을 작은 창에 그대로 적용하지 않습니다.

감쇠를 추측하지 말고 검증하기

코사인을 공간으로 두 번 미분하면 음의 파수 제곱이 나옵니다. 지수 진폭을 시간 미분하면 열확산율을 곱한 같은 인자가 생깁니다. 따라서 직접 대입으로 모든 점에서 열방정식이 성립함을 확인합니다. 초기 대비가 같더라도 파수가 클수록 더 미세한 공간 구조가 빨리 줄어듭니다.

$$T_t=-\alpha k^2B\cos(kx+\phi),\qquad T_{xx}=-k^2B\cos(kx+\phi),\qquad\tau=\frac1{\alpha k^2}.$$

감쇠 시간 tau는 열확산율이 양수일 때 유한합니다. 전도도가 영이면 초기 온도가 불균일해도 이 가족에서는 시간 감쇠가 없습니다. 초기 대비를 영으로 놓아 모든 위치의 온도를 같게 만드는 대조와는 다른 설정입니다.

열유속과 온도를 다른 관측으로 읽기

열유속은 국소 온도 자체가 아니라 그 기울기에 의존합니다. 온도가 최대인 점에서는 기울기가 영이므로 열유속도 영입니다. 반대로 평균 온도를 지나는 점에서 큰 열유속이 생길 수 있습니다. 열유속은 사인 무늬이므로 코사인 온도와 공간 위상이 다릅니다. 두 곡선의 최대점이 같은 위치에 있어야 한다고 요구하면 안 됩니다.

$$T_x=-kB\sin(kx+\phi),\qquad q=\kappa kB\sin(kx+\phi),\qquad qT_x=-\kappa T_x^2\le0.$$

마지막 부등식은 열전도의 방향을 확인합니다. 모든 유한 창이 같은 부호로 열을 얻거나 잃어야 한다는 뜻은 아닙니다. 각 창의 두 경계 열유속이 순교환을 결정합니다.

W로 계산하는 손계산 예

평균 온도300K, 대비0.1, 체적 열용량2J/(m³ K), 열전도도1W/(m K), 파수1rad/m, 단면적2m²를 선택합니다. 시간 영, x=pi/2 m에서 온도 섭동은 영이지만 기울기는−30K/m입니다. 따라서 양의 x 방향 열유속은30W/m²입니다. 실제 단면적을 곱하면 그 단면을 통과하는 전력은60W입니다.

$$B(0)=30\mathrm{\ K},\quad\alpha=0.5\mathrm{\ m^2/s},\quad x=\pi/2:\quad T=300\mathrm{\ K},\quad T_x=-30\mathrm{\ K/m},\quad q=30\mathrm{\ W/m^2}.$$

온도 최대점 x=0의 국소 열유속은 영입니다. 그래도 그 점을 중심으로 한 대칭 창은 양쪽 경계로 열을 내보낼 수 있습니다. 한 점의 열유속 측정으로 적분 수지를 대신할 수는 없습니다.

실제 유한 창 적분하기

중심과 음이 아닌 반폭을 실제m로 지정합니다. 창의 현열 에너지는 기준 온도와의 차이를 적분하고 체적 열용량과 단면적을 곱한 값입니다. 에너지 변화율은 오른쪽에서 왼쪽을 뺀 열 전력과 맞습니다. 양의 순수출 전력은 창 내부 에너지의 감소에 대응합니다. 이 수지는W를 사용하며W/m²와 구분합니다.

$$E_W=CA\int_{x_c-h}^{x_c+h}(T-T_{\rm ref})\,dx,\qquad\frac{dE_W}{dt}+A(q_R-q_L)=0.$$

섭동 적분은 거의 같은 사인 값을 빼는 대신 대칭 항등식으로 계산할 수 있습니다. sinc는 영에서1로 정의합니다. 부동소수점 끝점 좌표가 같은 수로 반올림되더라도 처음 지정한 물리적 반폭은 따로 유지해야 합니다.

$$E'_W=CA\,2hB\,\operatorname{sinc}(kh)\cos(kx_c+\phi),\qquad\operatorname{sinc}(z)=\frac{\sin z}{z},\quad\operatorname{sinc}(0)=1.$$

여기서 에너지의 프라임은 미분이 아니라 섭동 기여입니다. 손계산 설정에서 창을−pi/2부터+pi/2 m로 잡으면 이 기여는240J입니다. 초기 변화율은−120W이고 순수출 전력+120W와 맞습니다. 경계 방향을 빠뜨리면 양쪽에서 나가는 열을 잘못 상쇄하게 됩니다.

한 주기 전체의 현열 에너지 보존

주기 셀의 실제 체적은 단면적과 길이의 곱입니다. 기준 온도를273.15K로 정하고 그 기준에 대한 현열 에너지를 보고합니다. 에너지의 영점은 관례입니다. 이 기준보다 차가운 물질의 기준 상대 에너지가 음수라고 해서 절대 온도가 음수라는 뜻은 아닙니다. 일정한 열용량을 절대 영도까지 외삽할 필요도 없습니다.

$$L=\frac{2\pi}{k},\quad V=AL,\quad E_{\rm cell}=CV(\overline T-T_{\rm ref}),\quad\frac{dE_{\rm cell}}{dt}=0,\quad T_{\rm ref}=273.15\mathrm{\ K}.$$

이 정확 해 가족에서 섭동의 한 주기 순열 기여는 영입니다. 그래도 내부에는 열이 흐릅니다. 전체 에너지 보존과 비영의 내부 열유속은 서로 모순되지 않습니다.

다른 단위를 가진 이차 노름은 감소한다

온도 불균일성은 평균에서 벗어난 온도 제곱의 체적 적분의 절반으로 측정할 수 있습니다. 코사인 제곱의 평균은 절반입니다. 이 노름을 시간 미분하거나 기울기 제곱을 적분하면 같은 감쇠 관계를 얻습니다. 수학적 안정성을 보는 데 유용하지만 단위는J가 아니라K²m³입니다.

$$N=\frac12\int_V(T-\overline T)^2\,dV=\frac{VB^2}{4},\qquad\dot N=-\alpha\int_VT_x^2\,dV=-2\alpha k^2N.$$

손계산 설정의 초기 노름은 이 단위로900pi입니다. 그 감소는 온도가 매끄러워지는 것을 나타냅니다. 이 감소를 열에너지의 파괴라고 부르면 앞의 현열 보존식과 충돌합니다. 엔트로피 생성률 역시 별도의 물리량입니다. 이 노름의 이름만 바꾸어 엔트로피를 계산했다고 하지 않습니다.

대조 실험과 잔차가 말할 수 있는 범위

불균일한 초기 온도를 유지하고 전도도만 끄면 선택한 구성법칙의 열유속과 진화가 영입니다. 대신 대비를 영으로 놓으면 양의 전도도가 있어도 온도가 균일하므로 열유속과 진화가 사라집니다. 마지막으로 실제 비균일 전도장을 그대로 두고 창의 경계 전력만 평가에서 빼봅니다. 남는 에너지 변화율은 누락된 기여를 드러내며 원래 장이 보존을 위반했다는 뜻은 아닙니다.

$$R_{\rm omit\ boundary}=\frac{dE_W}{dt},\qquad R_{\rm full}=\frac{dE_W}{dt}+A(q_R-q_L).$$

같은 해석 식에서 유도한 두 표현이 일치하는 것은 일관성 진단입니다. 독립 미분, 공간 적분과 실제 온도의 유한 차분이 추가 근거가 됩니다. 어느 하나도 일반적인 불균질 열 솔버나 임의의 학습 구성법칙 전체를 인증하지는 않습니다.

작은 변화와 산술 관측

양의 감쇠 지수가 매우 작으면 표시된 지수 인자가1로 반올림될 수 있습니다. expm1로 변화를 직접 계산하면 반올림 인자에서1을 빼는 방법이 잃는 정보를 유지할 수 있습니다. 상대 변화와 절대 온도 진폭 변화의 표현 범위도 다릅니다. 작은 상대 변화에 큰 온도 진폭이 곱해지면 표현 가능한 절대 변화가 되기도 합니다.

$$z=\alpha k^2t,\qquad\Delta B=B(0)\operatorname{expm1}(-z),\qquad\Delta N=N(0)\operatorname{expm1}(-2z).$$

반대로 지수 인자만 먼저 계산하면 언더플로하지만 큰 진폭과 결합한 값은 표현 가능할 수 있습니다. 로그 공간에서 배율을 함께 처리합니다. 비영 관측을 표현하지 못하면 계산 불가로 남기고, 평균과 같을 수 있는 반올림된 전체 온도를 따로 둡니다. 폭이 정확히 영인 창의 섭동과 변화율은 영입니다. 반면 양의 폭의 두 끝점이 같은 수로 반올림되는 경우는 별개이며 지정한 실제 폭을 보존합니다.

표현 가능 여부는 관측마다 따로 판단합니다. 평균이300K이고 양의 온도 섭동이 표현 가능한 최소 수보다 작아졌다고 합시다. 전체 온도는 여전히300K로 반올림되어 표시할 수 있지만 섭동은 계산 불가입니다. 이 결과가 정확히 균일한 장이라는 뜻은 아닙니다. 유한 창에서도 모드의 에너지 기여를 표시할 수 없더라도 기준 상대 현열 에너지는 유한한 값으로 남을 수 있습니다. 기준 기여와 모드 기여를 먼저 합친 다음 총에너지의 표현 가능 여부를 판단합니다. 평균이 기준 온도273.15K와 정확히 같으면 기준 기여가 사라지므로, 남은 비영 모드가 표현 범위 아래일 때 전체 창 에너지도 계산 불가입니다. 대비나 폭을 정확히 영으로 둔 대조와 구분해야 합니다.

실험 순서와 AI 연구 활용

시간 영의 온도와 열유속 곡선부터 확인합니다. 온도 최대점과 평균을 지나는 점을 찾고 열이 흐르는 방향을 예측합니다. 따뜻한 영역 주위에 유한 창을 두고 에너지 변화율과 방향 있는 경계 전력을 비교합니다. 관측을 한 주기 전체로 넓혀 보존되는 현열과 감소하는 이차 노름을 대조합니다. 열확산율을 유지한 채 파수를 두 배로 늘리면 감쇠율이 네 배가 되는지 확인합니다. 전도도 영, 균일 초기 온도, 경계 전력 누락과 작은 비영 감쇠를 각각 비교합니다.

학습 모델에서는 온도 재구성과 실제 열유속을 따로 평가합니다. 그럴듯하게 매끄러워지는 온도를 예측해도 열전도도, 단면적이나 경계 방향이 틀릴 수 있습니다. 보지 않은 열확산율과 파수를 평가용으로 남기고 보존과 불균일성 감소를 따로 보고합니다. 거시 Fourier 모델의 수학적 순간 전파가 미시적·상대론적 보편 법칙을 뜻하지는 않습니다. 변하는 물성, 벽, 열원과 유체 수송 결합에는 뒤의 별도 가정과 모델이 필요합니다.

해설이 있는 연습

연습 1: 같은 열확산율

초기 온도를 유지하면서 열전도도와 체적 열용량을 둘 다 두 배로 바꾸면 무엇이 달라지나요?

해설: 열확산율과 온도 진행은 같습니다. 열유속, 저장 현열 에너지와 경계 전력은 두 배가 됩니다. 온도 영상만 맞추는 것으로 두 물성 계수를 독립적으로 식별할 수는 없습니다.

연습 2: 더 미세한 공간 구조

열확산율과 대비를 고정하고 파수를 두 배로 바꾸면 지수 감쇠율은 어떻게 되나요?

해설: 파수 제곱에 비례하므로 네 배입니다. 실제 주기는 절반이 되므로 셀 에너지를 비교할 때 그 기하 변화도 명시해야 합니다.

연습 3: 국소 열유속 영

따뜻한 창 중심의 온도 기울기가 영이면 창의 에너지 변화율도 반드시 영인가요?

해설: 아닙니다. 양쪽 경계의 열유속이 수지를 정합니다. 손계산의 대칭 창은 중심 열유속이 영이어도 두 끝으로 열을 내보냅니다.

연습 4: 감소하는 노름

이차 온도 노름은 줄고 전체 현열 에너지는 일정합니다. 보존 위반인가요?

해설: 아닙니다. 정의와 단위가 다른 양입니다. 열이 재분배되어 온도가 균일해집니다. K²m³의 노름을J의 열에너지로 보고하면 안 됩니다.

연습 5: 경계 전력 누락

고정 창의 에너지가 감소하는데 평가에서 경계 수송을 빠뜨렸습니다. 어떤 잔차가 남나요?

해설: 비영의 에너지 변화율이 남습니다. 바깥으로 흐르는 경계 전력을 복원하면 맞습니다. 누락 대조는 평가를 바꾸며 실제 온도장을 다른 것으로 바꾸지 않습니다.

연습 6: 변하지 않은 듯한 인자

양의 열전도도에서 지수 인자가1로 반올림되었습니다. 정확한 전도도 영 대조로 바꾸어도 되나요?

해설: 안 됩니다. 안정적 expm1 변화와 물성 입력, 산술적으로 계산 불가한 관측을 남깁니다. 반올림된 같음, 표현 범위 아래의 비영 변화, 정확히 영인 물리적 감쇠는 서로 다릅니다.

자료와 적용 범위

Gilbert Strang의 MIT OCW 열방정식 노트는 Fourier 모드 감쇠와 적분 보존·제곱 노름 감소의 구별을 다룹니다. 여기서는 차원을 가진 구성 모델, 기준 상대 현열 에너지와 유한 창 대조를 직접 유도했습니다. 다음 Gaussian 확산 단원은 관련된 수학 방정식을 다른 물리 단위의 농도에 적용합니다.