A multiple/scale/direction polynomial Trefftz method for solving the BHCP in high-dimensional arbitrary simply-connected domains. (January 2016)

• Record Type: Journal Article
• Authors: Liu, Chein-Shan
• Is Part Of: International journal of heat and mass transfer. Volume 92(2016:Jan.)
• Page Start: 970
• Page End: 978
• Publication Date: 2016-01
• Subjects: High-dimensional backward heat conduction problem -- Multi-dimensional heat polynomial -- Multiple/scale/direction polynomial Trefftz method -- Post-conditioner; Heat -- Transmission -- Periodicals; Mass transfer -- Periodicals; Chaleur -- Transmission -- Périodiques; Transfert de masse -- Périodiques; Electronic journals; 621.4022
• DOI: 10.1016/j.ijheatmasstransfer.2015.09.057
• Languages: English
• ISSNs: 0017-9310

Abstract

Highlights: Develop a Trefftz polynomial method for heat equation. Derive multiple/scale/direction polynomial Trefftz method for multi-dimensional backward heat conduction problem. The MSDPTM can retrieve the missing initial data very well. Abstract: The Trefftz method, endowing with polynomial type trial solutions of the high-dimensional heat conduction equation in an arbitrary simply-connected domain, is a cheap boundary-type meshless method to solve the high-dimensional backward heat conduction problem (BHCP), which is well-known a highly ill-posed problem. For solving the BHCP we develop a multiple/scale/direction polynomial Trefftz method (MSDPTM), of which the directions are uniformly distributed on a hyper-sphere and the scales are determined in advance by the collocation points on boundary. A post-conditioner for the resultant linear system is introduced by using the derived multiple-scale, and then use the conjugate gradient method (CGM) to solve the post-conditioned linear system, of which the expansion coefficients are determined fast. The MSDPTM can find the missing initial data very well, with several two- and three-dimensional numerical examples of the BHCP to evaluate the performance. Although under a very large relative noise from 30% to 80% with the maximum absolute noisy error being 2.17, the presented method still recovers the unknown initial temperature quite accurately. It is worthy to notice that the CPU time does not exceed two seconds for 2D examples and seven seconds for 3D examples.