TY - JOUR

T1 - Optimal hybrid parameter selection for stable sequential solution of inverse heat conduction problem

AU - Ahn, Chang uk

AU - Park, Chanhun

AU - Park, Dong Il

AU - Kim, Jin Gyun

N1 - Funding Information:
This research was funded by National Research Foundation of Korea ( NRF-2021R1A2C4087079 ), grant from Kyung Hee University in 2020 (KHU-20210488), and MOTIE, Korea Government (grant No. 20001228).
Publisher Copyright:
© 2021 Elsevier Ltd

PY - 2022/2

Y1 - 2022/2

N2 - To deal with the ill-posed nature of the inverse heat conduction problem (IHCP), the regularization parameter α can be incorporated into a minimization problem, which is known as Tikhonov regularization method, popular technique to obtain stable sequential solutions. Because α is a penalty term, its excessive use may cause large bias errors. A ridge regression was developed as an estimator of the optimal α to minimize the magnitude of a gain coefficient matrix appropriately. However, the sensitivity coefficient matrix included in the gain coefficient matrix depends on the time integrator; thus, certain parameters of the time integrators should be carefully considered with α to handle instability. Based on this motivation, we propose an effective iterative hybrid parameter selection algorithm to obtain stable inverse solutions. We considered the Euler time integrator to solve IHCP using the finite element method. We then considered β, a parameter to define Forward to Backward Euler time integrators, as a hybrid parameter with α. The error amplified by the inverse algorithm can be controlled by α first by assuming β=1. The total error is then classified into bias and variance errors. The bias error can be computed using the maximum heat flux change, and the variance error can be calculated using the measurement noise error generated by prior information. Therefore, α can initially be efficiently defined by the summation of the bias and variance errors computed in a time-independent manner. Reducing the total error for better stability of the inverse solutions is also available by adjusting β, which is defined to minimize the magnitude of gain coefficient matrix when spectral radius of the amplification matrix is less than one. Consequently, α could be updated with new β in the iteration process. The proposed efficient ridge estimator is essential to implement the iterative hybrid parameter selection algorithm in engineering practice. The possibility and performance of the hybrid parameter selection algorithm were evaluated by well-constructed 1D and 2D numerical examples.

AB - To deal with the ill-posed nature of the inverse heat conduction problem (IHCP), the regularization parameter α can be incorporated into a minimization problem, which is known as Tikhonov regularization method, popular technique to obtain stable sequential solutions. Because α is a penalty term, its excessive use may cause large bias errors. A ridge regression was developed as an estimator of the optimal α to minimize the magnitude of a gain coefficient matrix appropriately. However, the sensitivity coefficient matrix included in the gain coefficient matrix depends on the time integrator; thus, certain parameters of the time integrators should be carefully considered with α to handle instability. Based on this motivation, we propose an effective iterative hybrid parameter selection algorithm to obtain stable inverse solutions. We considered the Euler time integrator to solve IHCP using the finite element method. We then considered β, a parameter to define Forward to Backward Euler time integrators, as a hybrid parameter with α. The error amplified by the inverse algorithm can be controlled by α first by assuming β=1. The total error is then classified into bias and variance errors. The bias error can be computed using the maximum heat flux change, and the variance error can be calculated using the measurement noise error generated by prior information. Therefore, α can initially be efficiently defined by the summation of the bias and variance errors computed in a time-independent manner. Reducing the total error for better stability of the inverse solutions is also available by adjusting β, which is defined to minimize the magnitude of gain coefficient matrix when spectral radius of the amplification matrix is less than one. Consequently, α could be updated with new β in the iteration process. The proposed efficient ridge estimator is essential to implement the iterative hybrid parameter selection algorithm in engineering practice. The possibility and performance of the hybrid parameter selection algorithm were evaluated by well-constructed 1D and 2D numerical examples.

KW - Euler time integrator

KW - Finite element method

KW - Hybrid parameter selection

KW - Inverse heat conduction problem

KW - Morozov discrepancy principle

KW - Ridge estimator

KW - Tikhonov regularization

UR - http://www.scopus.com/inward/record.url?scp=85122514543&partnerID=8YFLogxK

U2 - 10.1016/j.ijheatmasstransfer.2021.122076

DO - 10.1016/j.ijheatmasstransfer.2021.122076

M3 - Article

AN - SCOPUS:85122514543

SN - 0017-9310

VL - 183

JO - International Journal of Heat and Mass Transfer

JF - International Journal of Heat and Mass Transfer

M1 - 122076

ER -