A method of fundamental solutions for two-dimensional heat conduction

We investigate an application of the method of fundamental solutions (MFS) to heat conduction in two-dimensional bodies, where the thermal diffusivity is piecewise constant. We extend the MFS proposed in Johansson and Lesnic [A method of fundamental solutions for transient heat conduction, Eng. Anal. Bound. Elem. 32 (2008), pp. 697–703] for one-dimensional heat conduction with the sources placed outside the space domain of interest, to the two-dimensional setting. Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate results can be obtained efficiently with small computational cost.

A method of fundamental solutions for transient heat conduction

In this paper we investigate an application of the method of fundamental solutions (MFS) to transient heat conduction. In almost all of the previously proposed MFS for time-dependent heat conduction the fictitious sources are located outside the time-interval of interest. In our case, however, these sources are instead placed outside the space domain of interest in the same manner as is done for stationary heat conduction. A denseness result for this method is discussed and the method is numerically tested showing that accurate numerical results can be obtained. Furthermore, a test example with boundary singularities shows that it is advisable to remove such singularities before applying the MFS.

A comparative study on applying the method of fundamental solutions to the backward heat conduction problem

We investigate an application of the method of fundamental solutions (MFS) to the backward heat conduction problem (BHCP). We extend the MFS in Johansson and Lesnic (2008) [5] and Johansson et al. (in press) [6] proposed for one and two-dimensional direct heat conduction problems, respectively, with the sources placed outside the space domain of interest. Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate and stable results can be obtained efficiently with small computational cost.

A method of fundamental solutions for transient heat conduction in layered materials

In this paper we investigate an application of the method of fundamental solutions (MFS) to transient heat conduction in layered materials, where the thermal diffusivity is piecewise constant. Recently, in Johansson and Lesnic [A method of fundamental solutions for transient heat conduction. Eng Anal Boundary Elem 2008;32:697–703], a MFS was proposed with the sources placed outside the space domain of interest, and we extend that technique to numerically approximate the heat flow in layered materials. Theoretical properties of the method, as well as numerical investigations are included.

Recovering boundary data in planar heat conduction using boundary integral equation method

We consider a Cauchy problem for the heat equation, where the temperature field is to be reconstructed from the temperature and heat flux given on a part of the boundary of the solution domain. We employ a Landweber type method proposed in [2], where a sequence of mixed well-posed problems are solved at each iteration step to obtain a stable approximation to the original Cauchy problem. We develop an efficient boundary integral equation method for the numerical solution of these mixed problems, based on the method of Rothe. Numerical examples are presented both with exact and noisy data, showing the efficiency and stability of the proposed procedure and approximations.

Computer methods for the heat conduction equation

The merits of various numerical methods for the solution of the one and two dimensional heat conduction equation with a radiation boundary condition have been examined from a practical standpoint in order to determine accuracies and efficiencies. It is found that the use of five increments to approximate the space derivatives gives sufficiently accurate results provided the time step is not too large; further, the implicit backward difference method of Liebmann (27) is found to be the most accurate method. On this basis, a new implicit method is proposed for the solution of the three-dimensional heat conduction equation with radiation boundary conditions. The accuracies of the integral and analogue computer methods are also investigated.

A method of fundamental solutions for radially symmetric and axisymmetric backward heat conduction problems

We propose and investigate an application of the method of fundamental solutions (MFS) to the radially symmetric and axisymmetric backward heat conduction problem (BHCP) in a solid or hollow cylinder. In the BHCP, the initial temperature is to be determined from the temperature measurements at a later time. This is an inverse and ill-posed problem, and we employ and generalize the MFS regularization approach [B.T. Johansson and D. Lesnic, A method of fundamental solutions for transient heat conduction, Eng. Anal. Boundary Elements 32 (2008), pp. 697–703] for the time-dependent heat equation to obtain a stable and accurate numerical approximation with small computational cost.

Modelling of transient heat transfer in annealing furnaces

This is a study of heat transfer in a lift-off furnace which is employed in the batch annealing of a stack of coils of steel strip. The objective of the project is to investigate the various factors which govern the furnace design and the heat transfer resistances, so as to reduce the time of the annealing cycle, and hence minimize the operating costs. The work involved mathematical modelling of patterns of gas flow and modes of heat transfer. These models are: Heat conduction and its conjectures in the steel coils;Convective heat transfer in the plates separating the coils in the stack and in other parts of the furnace; and Radiative and convective heat transfer in the furnace by using the long furnace model. An important part of the project is the development of numerical methods and computations to solve the transient models. A limited number of temperature measurements was available from experiments on a test coil in an industrial furnace. The mathematical model agreed well with these data. The model has been used to show the following characteristics of annealing furnaces, and to suggest further developments which would lead to significant savings: - The location of the limiting temperature in a coil is nearer to the hollow core than to the outer periphery. - Thermal expansion of the steel tends to open the coils, reduces their thermal conductivity in the radial direction, and hence prolongs the annealing cycle. Increasing the tension in the coils and/or heating from the core would overcome this heat transfer resistance. - The shape and dimensions of the convective channels in the plates have significant effect on heat convection in the stack. An optimal design of a channel is shown to be of a width-to-height ratio equal to 9. - Increasing the cooling rate, by using a fluidized bed instead of the normal shell and tube exchanger, would shorten the cooling time by about 15%, but increase the temperature differential in the stack. - For a specific charge weight, a stack of different-sized coils will have a shorter annealing cycle than one of equally-sized coils, provided that production constraints allow the stacking order to be optimal. - Recycle of hot flue gases to the firing zone of the furnace would produce a. decrease in the thermal efficiency up to 30% but decreases the heating time by about 26%.

Mechanisms of drying of particulate slurries

The literature on the evaporation of pure liquid drops and the drying of drops of solutions and slurries has been reviewed with particular reference to spray drying. A 0.1-0.2 mm glass filament-thermocouple was constructed and used to study simultaneously, heat and mass transfer from a single suspended drop placed in a humidity and temperature controlled, 28 mm OD vertical wind tunnel. Heat conduction through the filament was minimised eg at 100¦C it accounted for only 9.3% of the total heat transferred to a drop. Evaporation of single water drops was also studied in a 101 mm OD vertical wind tunnel. The Nusselt number was found to be a function of the Reynolds, Prandtl and Transfer number over an air temperature range of 17¦C to 107¦C. The proposed correlation is: Nu = 2+(-12.96B+0.76)Re¦-5Pr0-33 Experimental drying studies were carried out on single suspended 1 to 2.5 mm diameter drops of aqueous sodium sulphate decahydrate, sodium chloride, potassium sulphate, copper sulphate and sodium acetate solutions and slurries at temperatures of 20¦C to 124¦C. Dried crusts were examined by Scanning Electron Microscopy. The drying history of any material depended upon the nature of the crust formed. Sodium acetate formed a non-rigid skin prior to the formation of a rigid crust. A modified receding evaporation interface model was proposed for the drying of solutions and slurries. This covered both the constant rate period prior to crust formation and the subsequent falling rate period. The model was solved numerically for the variation in core temperature, drop weight and crust thickness. Good agreement was obtained between model predictions and experimental results for materials forming rigid crusts i.e. sodium sulphate decahydrate, sodium chloride, potassium sulphate and copper sulphate. However, the drying histories of drops of 10-20% weight initial concentration sodium acetate were unpredictable since formation of a non-rigid skin deviated from the model assumption of a rigid outer surface. At higher initial concentrations (40% weight) where a rigid crust was formed for sodium acetate, good agreement was obtained between experimental results and model predictions. Single suspended drop studies are concluded to provide a valuable insight into the drying mechanisms of specific solutions and slurries.

Identification of a time-dependent bio-heat blood perfusion coefficient

In the paper the identification of the time-dependent blood perfusion coefficient is formulated as an inverse problem. The bio-heat conduction problem is transformed into the classical heat conduction problem. Then the transformed inverse problem is solved using the method of fundamental solutions together with the Tikhonov regularization. Some numerical results are presented in order to demonstrate the accuracy and the stability of the proposed meshless numerical algorithm.

A method of fundamental solutions for theaxisymmetric backward heat equation

We propose and investigate an application of the method of fundamental solutions (MFS) to the radially symmetric and axisymmetric backward heat conduction problem (BHCP) in a solid or hollow cylinder. In the BHCP, the initial temperature is to be determined from the temperature measurements at a later time. This is an inverse and ill-posed problem, and we employ and generalize the MFS regularization approach [B.T. Johansson and D. Lesnic, A method of fundamental solutions for transient heat conduction, Eng. Anal. Boundary Elements 32 (2008), pp. 697–703] for the time-dependent heat equation to obtain a stable and accurate numerical approximation with small computational cost.

The drying of small drops of particulate slurries

The literature on the evaporation of drops of pure liquids, drops containing solids and droplet sprays has been critically reviewed. An experimental study was undertaken on the drying of suspended drops of pure water and aqueous sodium sulphate decahydrate with concentrations varying from 5 to 54. 1 wt. %. Individual drops were suspended from a glass filament balance in a 26 mm I.D. vertical wind tunnel, designed and constructed to supply hot de-humidified air, to simulate conditions encountered in commercial spray driers. A novel thin film thermocouple was developed to facilitate the simultaneous measurement of drop weight and core temperature. The heat conduction through the thermocouple was reduced because of its unique design; using essentially a single 50μ diameter nickel wire. For pure water drops, the Nusselt number was found to be a function of the Reynolds, Prandtl and Transfer numbers for a temperature range between 19 to 79°C.                  Nu = 2 + 0.19 (1/B)0.24 Re0.5 Pr0.33 Two distinct periods were observed during the drying of aqueous sodium sulphate decahydrate. The first period was characterised by the evaporation from a free liquid surface, whilst drying in the second period was controlled by the crust resistance. Fracturing of the crust occurred randomly but was more frequent at higher concentrations and temperatures. A model was proposed for the drying of slurry drops, based on a receding evaporation interface. The model was solved numerically for the variation of core temperature, drop weight and crust thickness as a function of time. Experimental results were in excellent agreement with the model predictions although at higher temperatures modifications to the model had to be made to accommodate the unusual behaviour of sodium sulphate slurries, i.e. the formation of hydrates.

Mechanisms of drying of skin forming materials

The literature relating to evaporation from single droplets of pure liquids, and to the drying of droplets containing solids and of droplet sprays has been reviewed. The heat and mass transfer rates for a single droplet suspended from a nozzle were studied within a 42mm I.D. horizontal wind tunnel designed to supply hot dry air, to simulate conditions encountered in a practical spray dryer. A novel rotating glass nozzle was developed to facilitate direct measurements of droplet weight and core temperature. This design minimised heat conduction through the nozzle. Revised correlations were obtained for heat and mass transfer coefficients, for evaporation from pure water droplets suspended from a rotating nozzle. Nu = 2.0 + 0.27 (l/B)°-18Re°-5Pr°-83 Sh = 2.0 + 0.575 ((T0-T.)/Tomfc) -o.o4Reo.5 ^0.33 Experimental drying studies were carried out on single droplets of different types of skin-forming materials, namely, custard, gelatin, skim milk and fructose at air temperatures ranging from 19°C to 198°C. Dried crusts were recovered and examined by Scanning Electron Microscopy. Skin-forming materials were classified into three types according to the mechanisms of skin formation. In the first type (typified by droplets of custard and starch) skin formed due to gelatinisation at high temperatures. Increasing the drying temperature resulted in increased crust resistance to mass transfer due to increased granule swelling and the crust resistance was completely transferred to a skin resistance at drying temperatures > 150°C. In the second type e.g. gelatin droplets the skin formed immediately drying had taken place at any drying temperature. At drying temperature > 60° C a more resistant skin was formed. In the third type (typified by droplets of skim milk and fructose) the skin appeared on the droplet surface at a certain stage of the drying process under any drying conditions. As the drying temperature was increased the resistance of the skin to mass transfer increased. The drying rate history of any material depended upon the nature of the skin formed which, in turn, depended upon the drying conditions. A mathematical model was proposed for the drying of the first type of skin-forming material. This was based on the assumption that, once all the granules gelatinised at the gelatinisation temperature, a skin appeared instantaneously on the droplet surface. The experimentally-observed times at which the skin appeared on the droplets surfaces were in excellent agreement with those predicted from the model. The work should assist in understanding the fundamentals of paniculate drying processes, particularly when skin-formation occurs and may be a crucial factor in volatiles retention.

Population lensing vs thermal lensing in Yb:YAG rods and disks

The induced lenses in the Yb:YAG rods and disks end-pumped by a Gaussian beam were analyzed both analytically and numerically. The thermally assisted mechanisms of the lens formation were considered to include: the conventional volume thermal index changes ("dn/dT"), the bulging of end faces, the photoelastic effect, and the bending (for a disk). The heat conduction equations (with an axial heat flux for a disk and a radial heat flux for a rod), and quasi-static thermoelastic equations (in the plane-stress approximation with free boundary conditions) were solved to find the thermal lens power. The population rate equation with saturation (by amplified spontaneous emission or an external wave) was examined to find the electronic lens power in the active elements.

Laminar heat transfer to non-Newtonian fluids in pipes with abrupt changes in diameter

In this thesis the results of experimental work performed to determine local heat transfer coefficients for non-Newtonian fluids in laminar flow through pipes with abrupt discontinuities are reported. The fluids investigated were water-based polymeric solutiorrs of time-indpendent, pseudoplastic materials, with flow indices "n" ranging from 0.39 to 0.9.The tube configurations were a 3.3 :1 sudden convergence, and a 1: 3.3 sudden divergence.The condition of a prescribed uniform wall heat flux was considered, with both upstream and downstream tube sections heated. Radial temperature traverses were also under­ taken primarily to justify the procedures used in estimating the tube wall and bulk fluid temperatures and secondly to give further insight into the mechanism of heat transfer beyond a sudden tube expansion. A theoretical assessment of the influence of viscous dissipation on a non-Newtonian pseudoplastic fluid of' arbitrary index "n" was carried out. The effects of other secondary factors such as free convection and temperature-dependent consistency were evaluated empirically. In the present investigations, the test conditions were chosen to minimise the effects of natural convection and the estimates of viscous heat generation showed the effect to be insignificant with the polymeric concentrations tested here. The final results have been presented as the relationships between local heat transfer coef'ficient and axial distance downstream of the discontinuities and relationships between dimensionless wall temperature and reduced radius. The influence of Reynolds number, Prandtl number, non-Newtonian index and heat flux have been indicated.