열탄성의 상호 에너지 교환

온도장을 입력 하나 더 넣는 것으로 충분한가

고체를 가열하면 팽창하려는 경향이 있다. 그래서 먼저 열 방정식을 풀고 그 온도로 변형을 계산하는 방법을 생각할 수 있다. 어떤 조건에서는 유용한 근사다. 그러나 변형도 온도를 바꿀 수 있다. 결합 모델이라면 어느 방향의 영향을 포함하는지, 두 장 사이에서 에너지가 어떻게 교환되는지 밝혀야 한다.

여기서는 선형 열탄성체의 공간 모드 하나를 다룬다. 변위·속도·온도 편차의 세 진폭만으로 상호 교환을 확인하고 일방향 근사와 비교할 수 있다. 실제 공간·시간 적분까지 대조하므로 과학 AI 대리모델을 평가하는 작은 실험이 된다. 두 곡선을 잘 맞추는 것만으로 두 장의 교환 항이 일치한다고 결론 내릴 수는 없다.

물성·기준 상태·영역을 선언한다

축방향 변위를 u,그 공간 미분을 변형률 ε,양의 기준 절대온도 T0에서 벗어난 양을 θ라 하자. 탄성계수 E,밀도 ρ,체적 열용량 C는 양수다. 열전도도 κ는0 이상이고 열응력 계수 β는 상수다. C는 일정 변형률에서의 열용량이며 일정 응력의 열용량과 구별한다. 결합된 물질에서는 두 열용량이 다를 수 있다.

$$\varepsilon=u_x,\qquad T=T_0+\theta,\qquad [E]=\mathrm{Pa},\quad[C]=\mathrm{J/(m^3K)},\quad [\beta]=\mathrm{Pa/K},\quad[\kappa]=\mathrm{W/(mK)}.$$

작은 변형률과 작은 온도 편차,상수 물성을 가정한다. 소성·큰 변형·기하 변화나 실제3차원 시편의 보정된 물성은 다루지 않는다. 양의 파수 k와 길이2π/k인 주기 셀,단면적 A를 사용한다. k는 파수이고 κ는 열전도도다. k=0으로 바꾸면 유한한 주기 셀이 사라지므로 이 가족의 입력에 포함하지 않는다.

하나의 열역학 선택에서 양방향을 유도한다

이차 Helmholtz 자유에너지 전개를 선택하면 응력과 엔트로피 편차를 일관되게 정할 수 있다. 여기의 에너지·엔트로피는 기준 체적당 양이다. 기준 상태에서 일정하거나 일차인 항은 이 섭동 표현에서 생략했다.

$$\psi=\frac E2\varepsilon^2-\beta\varepsilon\theta -\frac C{2T_0}\theta^2,\qquad \sigma=\psi_{\varepsilon}=E\varepsilon-\beta\theta,\qquad s=-\psi_{\theta}=\beta\varepsilon+\frac C{T_0}\theta.$$

자유에너지의 온도 곡률이 음수인 것은 양의 열용량과 일치한다. 저장된 열의 총량이 음수라는 뜻은 아니다. 자유에너지와 내부 에너지는 사용하는 독립 변수와 역할이 다르다. 특정 에너지 함수의 부호만으로 모든 물리적 저장량을 판단하지 않도록 주의하자.

운동량 보존에서 역학식을 얻는다. 엔트로피 관계를 기준 온도 근처에서 선형화하고 Fourier 열전도를 적용하면 열 방정식을 얻는다. 응력에 나타난 것과 같은 β가 역방향 열 결합에도 들어간다.

$$\rho u_{tt}=E u_{xx}-\beta\theta_x,\qquad C\theta_t=\kappa\theta_{xx}-T_0\beta\varepsilon_t.$$

β가 양수인 단열 과정에서는 늘어나는 변형이 온도 편차를 낮출 수 있다. 그래프의 에너지를 감소시키려고 부호를 임의로 고른 것이 아니다. 선언한 열역학 관계에서 나온 부호다. 관례를 바꾼다면 연결된 모든 양도 일관되게 바꾸어야 한다.

PDE를 실제 상태 변수 세 개로 줄인다

변위에 sine 모드,온도에 cosine 모드를 선택한다. 공간 미분이 서로 결합하는 데 필요한 위상을 만들기 때문이다. 원래 방정식에 직접 대입하면 다음 닫힌 선형 진폭 계를 얻는다.

$$u=a(t)\sin(kx),\quad\theta=b(t)\cos(kx),\quad v=\dot a, \qquad \dot a=v,\quad \rho\dot v=-Ek^2a+\beta kb,\quad C\dot b=-\kappa k^2b-T_0\beta kv.$$

초기 변위 진폭·속도 진폭·온도 진폭은 서로 다른 초기 데이터다. 변위의 한 장면만으로 속도와 온도까지 결정할 수 없다. b가 음수라는 것은 기준보다 낮은 온도 편차이며 음의 절대온도가 아니다. 교육용 입력은 작은 변형·온도 편차로 제한한다. 상호 결합 모델의 에너지 상한은 전체 이력의 작은 변형 조건을 추가로 점검하는 데 사용할 수 있다.

실제 위치의 변위와 온도는 진폭 자체와도 다르다. 예를 들어 kx=0에서는 변위 sine 값은0이지만 온도 편차는 b다. 한 점의 변위가0이라고 전체 모드나 운동이 없다고 판단하면 안 된다.

의미가 분명한 에너지를 선택한다

이 선형화 문제에서는 운동에너지와 ψ+θs의 이차 부분을 합친다. 혼합 항은 상쇄된다. 실제 주기 셀에서 sine²와 cosine²를 적분하면 각각 평균1/2을 얻는다. 셀 체적 V=2πA/k를 사용한 결과가 다음 가용 에너지다.

$$\mathcal A_m=\frac V4(\rho v^2+Ek^2a^2),\qquad \mathcal A_T=\frac{VC}{4T_0}b^2,\qquad \mathcal A=\mathcal A_m+\mathcal A_T.$$

가용 에너지는 기준 상태에 대한 교란의 에너지를 측정한다. 전체 내부 에너지,누적 열량 또는 C에 온도 진폭을 곱한 값과 같지 않다. b의 부호와 무관하게 열 가용 에너지는 비음수다. 열전도는 교란을 평탄하게 만들며 가용 에너지를 줄일 수 있지만 물리 에너지가 사라진다는 뜻은 아니다. 이 선형 섭동 모델은 평균 가열과 엔트로피의 모든 이차 기여를 추적하지 않는다.

단면적만 바꾸면 진폭 운동은 그대로이고 두 적분 에너지가 함께 바뀐다. 파수를 바꾸면 진동·확산 속도뿐 아니라 실제 주기 길이와 체적도 바뀐다. 체적 인자를 잊으면 서로 다른 실험의 에너지를 잘못 비교하게 된다. 밀도와 체적을 이미 포함한 값에 다시 곱하는 중복 계산도 피해야 한다.

교환 항은 상쇄되고 전도 소산은 남는다

진폭 방정식으로 두 에너지를 미분하자. 역학 교환은 양수나 음수일 수 있어 어느 방향으로도 가용 에너지를 전달할 수 있다. 상호 열 교환은 그 반대 부호이고 전도에 의한 가용 에너지 소산은 비음수다.

$$\dot{\mathcal A}_m=\frac{V\beta k}{2}vb,\qquad \dot{\mathcal A}_T=-\frac{V\beta k}{2}vb -\frac{V\kappa k^2}{2T_0}b^2,\qquad \dot{\mathcal A}=-\frac{V\kappa k^2}{2T_0}b^2.$$

열전도도가0이면 총 가용 에너지는 보존되어도 역학·열 부분은 각각 변할 수 있다. β가0이면 탄성 진동과 확산하는 온도 모드가 분리된다. 두 극한은 다르다. 전도를 제거한다고 상호 교환도 없어지는 것은 아니다. 역학 부분만 일시적으로 증가했다는 이유로 전체 계가 능동적으로 에너지를 만든다고 해석하지 않는다.

상호 결합에서는 각 에너지의 양수성으로 변형률 진폭의 상한도 얻는다. 이 상한은 이 모델에 적용되며,다른 구동을 받거나 일방향으로 바꾼 계에 자동으로 적용되지는 않는다.

$$|ka(t)|\le\sqrt{\frac{4\mathcal A(0)}{VE}}.$$

작은 변화를 지우지 않고 계산한다

에너지로 정규화한 좌표를 사용하면 결합 행렬의 의미가 선명해진다. 역학·열 사이의 상호 항은 반대 부호를 갖는다.

$$y=\begin{pmatrix}k\sqrt E\,a\\\sqrt\rho\,v\\\sqrt{C/T_0}\,b\end{pmatrix},\qquad \dot y=\begin{pmatrix}0&\omega&0\\-\omega&0&g\\0&-g&-d\end{pmatrix}y, \quad\omega=k\sqrt{E/\rho},\quad g=\beta k\sqrt{T_0/(\rho C)},\quad d=\kappa k^2/C.$$

구현은 작은 구간으로 조절한 Taylor 전개와 반복 제곱으로 행렬 지수를 수치 계산한다. 지수 행렬에서 항등행렬을 뺀 변화를 직접 전파하므로,반올림된 두 상태를 뺀 결과로 작은 변화를 정의하지 않는다. 전체 전파 행렬도 별도로 유지한다. 그래야−1에 가까운 변화를1에 다시 더할 때 사라질 수 있는 작은 감쇠 꼬리를 보존한다. 부호 있는 로그는 일반 부동소수점 범위 밖의 비영 진폭과 이차 양을 보존한다. 이 방식이 모든 실수 입력의 정확한 반올림을 증명하는 것은 아니다.

유한 시간 구간에서는 실제 이력의 교환·전도 소산을 별도로 수치 적분한다. 수지로부터 얻은 가용 에너지 변화와 끝점 뺄셈을 구분해서 반환한다. 같은 항등식에서 구성한 두 양이 맞는다는 사실은 진단이며 독립 정확도 인증이 아니다. 참조 계산에서는 물리 변수의 행렬 지수,별도 고유계와 실제 공간 적분을 사용한다.

초기 변화율을 손으로 계산한다

E=2Pa,ρ=1kg/m³,C=1J/(m³ K),T0=300K,β=0.02Pa/K,κ=0.1W/(m K),k=1/m,A=1m²를 택하자. 초기 a=0.001m,v=0.001m/s,b=0.02K다. 실제 재료의 보정된 물성이 아니라 교육용 계수다.

셀 체적은2πm³다. 역학 가용 에너지는 약4.71239μJ,열 가용 에너지는 약2.09440μJ다. 위상0의 응력은0.0016Pa다. 초기 역학 교환률은 약+1.25664μW,전도에 의한 가용 에너지 소산은 약0.418879μW다. 열 부분의 변화율은 약−1.67552μW이므로 합은 전도 소산의 음수와 같다. 역학 부분이 처음 증가해도 전체 수지는 일치한다.

직접 숫자를 확인할 때에는 μJ와 μW를 구별하자. 한쪽은 저장량이고 다른 쪽은 시간당 변화다. 현재 시점의 소산율을 시간 구간 전체의 누적 소산과 바꾸어 읽으면 안 된다.

일방향 근사를 정직하게 비교한다

온도는 확산만으로 진화시키면서 역학 응력은 계속 온도의 영향을 받도록 바꾸어보자. 열 진폭 식은 b′=−db가 된다. 역학 교환은 남으므로 닫힌 상호 결합의 이전 수지도 달라진다.

$$\text{one-way:}\quad \dot b=-db,\qquad \dot{\mathcal A}=\frac{V\beta k}{2}vb -\frac{V\kappa k^2}{2T_0}b^2.$$

이는 다른 모델이고 실제 궤적도 달라진다. 역작용이 약한 근사로 쓰거나 외부 제어기가 온도를 유지하는 상황으로 해석할 수 있다. 후자의 해석이라면 적절한 외부 열 입력을 선언해야 한다. 같은 닫힌 상호 결합 계라고 부르면서 설명되지 않은 에너지 증가를 남기면 안 된다. 기존 궤적에서 누락 항의 잔차만 계산하는 것과 변경한 계를 실제로 푸는 것도 구분해야 한다.

시작 직후의3차 소산을 이해한다

변위와 온도 진폭은0이고 초기 속도 진폭만 비영인 경우를 생각하자. 처음 생성되는 온도는 시간에 일차다. 전도 소산율은 이차이고 누적 소산은 삼차가 된다.

$$b(t)\sim-\frac{T_0\beta k v_0}{C}t,\qquad \int_0^t\!D(s)\,ds \sim\frac{V\kappa T_0\beta^2k^4v_0^2}{6C^2}t^3.$$

끝점 에너지가 초기 값과 같게 반올림되어도 양의3차 소산은 수학적으로 남을 수 있다. 안정 적분은 일반 숫자 표시가 불가할 때도 그 로그를 보존한다. 정확한 영 결합,영 전도도,영 초기 상태는 별도의 물리적 경우이며 해당하는 양의 정확한0을 유지해야 한다.

풀이가 있는 여섯 문제

문제 1 — 부족한 초기 데이터

초기 변위만 알면 이 진화를 유일하게 정할 수 있는가?

풀이: 아니다. 변위·속도·온도 진폭을 각각 지정해야 한다. 세 값은 다른 상태 좌표에 속하며 서로를 자동으로 결정하지 않는다.

문제 2 — 음의 온도 진폭

b<0이면 열 가용 에너지나 절대온도가 음수인가?

풀이: 둘 다 아니다. 열 가용 에너지는 b²에 비례한다. 절대온도는 양의 T0에 작은 공간 편차를 더한 값이다.

문제 3 — 단열 교환

전도도를0으로 하고 β를 유지하면 역학 가용 에너지가 일정해야 하는가?

풀이: 아니다. 상호 교환으로 두 부분이 각각 변할 수 있다. 합은 일정하고 전도 소산은 정확히0이다.

문제 4 — 면적 배율

물성·파수·초기 진폭을 유지하고 면적만2배 하면 무엇이 바뀌는가?

풀이: 진폭 궤적은 그대로다. 셀 체적,적분 가용 에너지,교환 전력과 소산 전력은2배다.

문제 5 — 일방향 온도

확산만 따르는 온도가 역학을 구동하는 모델에 닫힌 상호 가용 에너지 수지를 그대로 써도 되는가?

풀이: 안 된다. 역방향 교환을 제거했다. 근사 또는 외부 열 제어기를 선언하고 남은 교환 항을 수지에 포함해야 한다.

문제 6 — 같게 반올림된 에너지

초기 온도가0이고 속도가 비영일 때 끝점 에너지가 같게 표시되면 초기 전도 소산이0이라는 증거인가?

풀이: 아니다. 생성 온도는 처음 시간에 일차,누적 소산은 삼차다. 끝점 뺄셈과 별도로 안정 적분과 로그를 확인해야 한다.

자료와 연구 적용 범위

P. Kelly의 Solid Mechanics Part IV §2.2는 열탄성 자유에너지와 열전도식의 변형 기여를 전개한다. MOOSE의 열·역학 튜토리얼은 두 물리 구성요소의 연결 예다. 이 단원의 일차원 주기 가족과 가용 에너지 수지는 선언한 조건에서 직접 유도했다. AI 연구에서는 상호 교환·초기 데이터·매개변수 척도·실제 공간/시간 적분을 검사한 뒤 누락 모드와 비선형 효과로 범위를 넓혀야 한다. 단일 모드의 적합 성공을 범용 열탄성 해석기로 확대하지 않는다.