[논문 리뷰] Implicit Euler time discretization and FDM with Newton method in nonlinear heat transfer modeling

Stefan M. Filipov, István Faragó|arXiv (Cornell University)|2018. 11. 15.

Advanced Numerical Analysis Techniques참고 문헌 3인용 수 11

한 줄 요약

이 논문은 시간에 따라 변화하는 비선형 열전달 문제를 해결하기 위한 새로운 수치적 접근법을 제안한다. 먼저 시간에 대해 음성 오일러 방법을 적용하여 편미분방정식(PDE)을 비선형 두점 경계값 문제의 연속으로 변환한다. 이후 공간적 해석에는 유한차분법(FDM)을 사용하고, 유도된 비선형 시스템은 뉴턴-라프슨 방법으로 해결하며, 비선형 항의 해석적 도함수를 유도하여 수렴의 안정성을 확보한다. MATLAB을 사용하여 반도체 열적 데이터를 기반으로 검증하였으며, 복잡한 열적 거동을 모델링하는 데 있어 안정성과 정확성을 입증하였다.

ABSTRACT

This paper considers one-dimensional heat transfer in a media with temperature-dependent thermal conductivity. To model the transient behavior of the system, we solve numerically the one-dimensional unsteady heat conduction equation with certain initial and boundary conditions. Contrary to the traditional approach, when the equation is first discretized in space and then in time, we first discretize the equation in time, whereby a sequence of nonlinear two-point boundary value problems is obtained. To carry out the time-discretization, we use the implicit Euler scheme. The second spatial derivative of the temperature is a nonlinear function of the temperature and the temperature gradient. We derive expressions for the partial derivatives of this nonlinear function. They are needed for the implementation of the Newton method. Then, we apply the finite difference method and solve the obtained nonlinear systems by Newton method. The approach is tested on real physical data for the dependence of the thermal conductivity on temperature in semiconductors. A MATLAB code is presented.

연구 동기 및 목표

• 온도에 따라 변화하는 열전도율을 갖는 재료에서 일시적 열전도를 모델링하는 데 도전하는 것.
• 시간과 공간 이산화 단계의 순서를 재정렬함으로써 비선형 PDE에서 발생하는 수치적 안정성과 정확도를 향상시키는 것.
• 시간 이산화된 PDE에서 유도된 비선형 시스템을 뉴턴 방법과 해석적 자코비안을 사용하여 안정적인 해법을 개발하는 것.
• 실제 물리적 데이터를 기반으로 한 반도체 열전도율 데이터를 사용하여 방법의 실용적 관련성을 확보하는 것.
• 재현 가능성과 공 ingeneering 응용을 위해 완전한 MATLAB 구현을 제공하는 것.

제안 방법

• 시간 도함수를 이산화하기 위해 음성 오일러 스킴을 적용하여 정 steady-state 비선형 두점 경계값 문제의 연속으로 비정상 상태 PDE를 변환한다.
• 공간적 제2도미분을 유한차분법으로 이산화하여 열전도율이 온도와 온도 기울기에 따라 비선형적으로 의존하는 특성을 유지한다.
• 열전도율을 나타내는 비선형 함수의 해석적 편도함수를 유도하여 뉴턴-라프슨 업데이트의 정확성을 보장한다.
• 시간 및 공간 이산화된 PDE에서 유도된 비선형 방정식 시스템을 수립하고, 유도된 자코비안을 사용하여 반복적으로 뉴턴의 방법으로 해결한다.
• 각 시간 단계에서 비선형 시스템을 뉴턴 방법으로 해결하여 이차 수렴성과 안정성을 확보한다.
• 시간 스텝, 공간 이산화, 선형 검색을 포함한 선형화된 뉴턴 반복을 포함한 전체 알고리즘을 MATLAB에 구현한다.

실험 결과

연구 질문

• RQ1 음성 오일러 시간 이산화 후 FDM과 뉴턴의 방법을 적용한 방법이 전통적인 연산자 분리 접근법에 비해 비선형 열전달 문제 해결에 있어 어떻게 비교되는가?
• RQ2 이 맥락에서 비선형 열전도율 함수의 해석적 도함수가 뉴턴의 방법의 수렴성과 안정성에 개선을 가져오는가?
• RQ3 실제 반도체 열전도율 데이터에 적용했을 때 제안된 방법의 정확도와 안정성은 어떠한가?
• RQ4 선형화나 단순화 없이 강한 온도 의존성에 의한 열전도율 변화를 어떻게 다루는가?
• RQ5 표준 MATLAB 프레임워크를 사용하여 제안된 방법을 효율적으로 구현하고 재현할 수 있는가?

주요 결과

• 제안된 방법은 시간-공간 이산화 전략을 사용하여 온도에 따라 변화하는 열전도율을 갖는 일시적 열전도를 성공적으로 모델링한다.
• 뉴턴-라프슨 방법에 해석적 도함수를 사용함으로써 각 시간 단계에서 비선형 해법의 빠른 이차 수렴성을 보장한다.
• 실제 반도체 열전도율 데이터에 적용했을 때 수치적 안정성과 정확성이 검증되었으며 결과로 입증되었다.
• MATLAB에서의 구현은 연구자와 엔지니어들이 사용할 수 있고 재현 가능한 프레임워크를 제공한다.
• 이 방법은 비선형 문제에서 자주 발생하는 연산자 분리 방법의 정확성 부족과 불안정성 문제를 피한다.
• 수치적 결과는 강한 열전도율 비선형성 특성을 갖는 재료에서 열전달의 물리적 거동을 정확히 포착하고 있음을 확인한다.