Krzysztof Grysa and Artur Maciag
The purpose of this paper is to derive the Trefftz functions (T-functions) for the Pennes’ equation and for the single-phase-lag (SPL) model (hyperbolic equation) with perfusion…
Abstract
Purpose
The purpose of this paper is to derive the Trefftz functions (T-functions) for the Pennes’ equation and for the single-phase-lag (SPL) model (hyperbolic equation) with perfusion and then comparing field of temperature in a flat slab made of skin in the case when perfusion is taken into account, with the situation when a Fourier model is considered. When considering the process of heat conduction in the skin, one needs to take into account the average values of its thermal properties. When in biological bodies relaxation time is of the order of 20 s, the thermal wave propagation appears. The initial-boundary problems for Pennes’ model and SPL with perfusion model are considered to investigate the effect of the finite velocity of heat in the skin, perfusion and thickness of the slab on the rate of the thermal wave attenuation. As a reference model, the solution of the classic Fourier heat transfer equation for the considered problems is calculated. A heat flux has direction perpendicular to the surface of skin, considered as a flat slab. Therefore, the equations depend only on time and one spatial variable.
Design/methodology/approach
First of all the T-functions for the Pennes’ equation and for the SPL model with perfusion are derived. Then, an approximate solutions of the problems are expressed in the form of a linear combination of the T-functions. The T-functions satisfy the equation modeling the problem under consideration. Therefore, approximating a solution of a problem with a linear combination of n T-functions one obtains a function that satisfies the equation. The unknown coefficients of the linear combination are obtained as a result of minimization of the functional that describes an inaccuracy of satisfying the initial and boundary conditions in a mean-square sense.
Findings
The sets of T-functions for the Pennes’ equation and for the SPL model with perfusion are derived. An infinite set of these functions is a complete set of functions and stands for a base functions layout for the space of solutions for the equation used to generate them. Then, an approximate solutions of the initial-boundary problem have been found and compared to find out the effect of finite velocity of heat in the skin, perfusion and thickness of the slab on the rate of the thermal wave attenuation.
Research limitations/implications
The methods used in the literature to find an approximate solution of any bioheat transfer problems are more complicated than the one used in the presented paper. However, it should be pointed out that there is some limitation concerning the T-function method, namely, the greater number of T-function is used, the greater condition number becomes. This limitation usually can be overcome using symbolic calculations or conducting calculations with a large number of significant digits.
Originality/value
The T-functions for the Pennes’ equation and for the SPL equation with perfusion have been reported in this paper for the first time. In the literature, the T-functions are known for other linear partial differential equations (e.g. harmonic functions for Laplace equation), but for the first time they have been derived for the two aforementioned equations. The results are discussed with respect to practical applications.
Details
Keywords
Michał Ciałkowski, Aleksander Olejnik, Magda Joachimiak, Krzysztof Grysa and Andrzej Frąckowiak
To reduce the heat load of a gas turbine blade, its surface is covered with an outer layer of ceramics with high thermal resistance. The purpose of this paper is the selection of…
Abstract
Purpose
To reduce the heat load of a gas turbine blade, its surface is covered with an outer layer of ceramics with high thermal resistance. The purpose of this paper is the selection of ceramics with such a low heat conduction coefficient and thickness, so that the permissible metal temperature is not exceeded on the metal-ceramics interface due to the loss ofmechanical properties.
Design/methodology/approach
Therefore, for given temperature changes over time on the metal-ceramics interface, temperature changes over time on the inner side of the blade and the assumed initial temperature, the temperature change over time on the outer surface of the ceramics should be determined. The problem presented in this way is a Cauchy type problem. When analyzing the problem, it is taken into account that thermophysical properties of metal and ceramics may depend on temperature. Due to the thin layer of ceramics in relation to the wall thickness, the problem is considered in the area in the flat layer. Thus, a one-dimensional non-stationary heat flow is considered.
Findings
The range of stability of the Cauchy problem as a function of time step, thickness of ceramics and thermophysical properties of metal and ceramics are examined. The numerical computations also involved the influence of disturbances in the temperature on metal-ceramics interface on the solution to the inverse problem.
Practical implications
The computational model can be used to analyze the heat flow in gas turbine blades with thermal barrier.
Originality/value
A number of inverse problems of the type considered in the paper are presented in the literature. Inverse problems, especially those Cauchy-type, are ill-conditioned numerically, which means that a small change in the inputs may result in significant errors of the solution. In such a case, regularization of the inverse problem is needed. However, the Cauchy problem presented in the paper does not require regularization.