Solving Transient Nonlinear Heat Conduction Problems by Proper Orthogonal Decomposition and the Finite-Element Method

ABSTRACT

A method of reducing the number of degrees of freedom and the overall computing time by combining proper orthogonal decomposition (POD) with the finite-element method (FEM) has been devised. The POD-FEM technique can be applied both to linear and nonlinear problems. At the first stage of the method a standard FEM time-stepping procedure is invoked. The temperature fields obtained for the first few time steps undergo statistical analysis, yielding an optimal set of globally defined trial and weighting functions for the Galerkin solution of the problem at hand. The resulting set of ordinary differential equations (ODEs) is of greatly reduced dimensionality when compared with the original FEM formulation. For linear problems, the set can be solved either analytically, resorting to the modal analysis technique, or by time stepping. In the case of nonlinear problems, only time stepping can be applied. The focus of this article is on the time-stepping approach, in which the generation of the FEM-POD matrices, requiring some additional matrix manipulations, can be embedded in the assembly of standard FEM matrices. The gain in execution times comes from the significantly shorter time of solution of the set of algebraic equations at each time step. Numerical results are presented for both linear and nonlinear problems. In the case of linear problems, the derived time-stepping technique is compared with the standard FEM and the modal analysis. For nonlinear problems the proposed POD-FEM approach is compared with the standard FEM. Good accuracy of the POD-FEM solver has been observed. Controlling the error introduced by the reduction of the degrees of freedom in POD is also discussed.

This work was supported by the National Committee for Fundamental Research, Poland, under grant 135/E-367/SPB/5.PR UE/DZ 348.

A preliminary version of the article was presented at CHT-04: An ICHMT International Symposium on Advances in Computational Heat Transfer, April 2004, G. de Vahl Davis and E. Leonardi (eds.), CD-ROM Proceedings, ISBN 1-5670-174-2, Begell House, New York, 2004.

Notes

^a Capacitance matrix was calculated only once.