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

被引:12
作者
Ahn, Chang-uk [1 ,2 ]
Park, Chanhun [2 ]
Park, Dong Il [2 ]
Kim, Jin-Gyun [1 ]
机构
[1] Kyung Hee Univ, Dept Mech Engn Integrated Engn, 1732 Deogyeong Daero, Yongin 17104, Gyeonggi Do, South Korea
[2] Korea Inst Machinery & Mat, Dept Robot & Mechatron Res, Daejeon 34103, South Korea
基金
新加坡国家研究基金会;
关键词
Inverse heat conduction problem; Tikhonov regularization; Hybrid parameter selection; Ridge estimator; Finite element method; Euler time integrator; Morozov discrepancy principle; SENSITIVITY COEFFICIENT METHOD; GENERALIZED CROSS-VALIDATION; PARTICLE SWARM OPTIMIZATION; OPTIMAL REGULARIZATION; STABILITY ANALYSIS; REDUCTION;
D O I
10.1016/j.ijheatmasstransfer.2021.122076
中图分类号
O414.1 [热力学];
学科分类号
摘要
To deal with the ill-posed nature of the inverse heat conduction problem (IHCP), the regularization pa-rameter alpha can be incorporated into a minimization problem, which is known as Tikhonov regularization method, popular technique to obtain stable sequential solutions. Because alpha is a penalty term, its exces-sive use may cause large bias errors. A ridge regression was developed as an estimator of the optimal alpha 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 alpha 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 consid-ered beta, a parameter to define Forward to Backward Euler time integrators, as a hybrid parameter with alpha. The error amplified by the inverse algorithm can be controlled by alpha first by assuming beta = 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, alpha 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 beta, which is defined to minimize the magnitude of gain coefficient matrix when spectral radius of the amplification matrix is less than one. Consequently, alpha could be updated with new beta 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. (c) 2021 Elsevier Ltd. All rights reserved.
引用
收藏
页数:14
相关论文
共 60 条
[1]  
Ahn J.G., 2021, ARXIV210501304
[2]  
Alifanov O. M., 1989, Journal of Engineering Physics, V56, P308, DOI 10.1007/BF00871169
[3]  
Alifanov O. M., 1994, International series in heat and mass transfer, V1st
[4]  
[Anonymous], 1983, Improperly Posed Problems and their Numerical Treatment
[5]   Exact solution of heat conduction in composite materials and application to inverse problems [J].
Aviles-Ramos, G ;
Haji-Sheikh, A ;
Beck, JV .
JOURNAL OF HEAT TRANSFER-TRANSACTIONS OF THE ASME, 1998, 120 (03) :592-599
[6]  
Bathe K., 1995, Finite Element Procedures
[7]  
Bathe KJ., 1982, FINITE ELEMENT PROCE, DOI DOI 10.1115/1.3264375
[8]   Comparing parameter choice methods for regularization of ill-posed problems [J].
Bauer, Frank ;
Lukas, Mark A. .
MATHEMATICS AND COMPUTERS IN SIMULATION, 2011, 81 (09) :1795-1841
[9]  
Beck J., 1985, Inverse Heat Conduction: Ill-posed Problems
[10]   Inverse problems and parameter estimation: integration of measurements and analysis [J].
Beck, JV ;
Woodbury, KA .
MEASUREMENT SCIENCE AND TECHNOLOGY, 1998, 9 (06) :839-847