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.

Advanced thermal analysis of microelectronics using spreading resistance models

Thermal analysis of electronic devices is one of the most important steps for designing of modern devices. Precise thermal analysis is essential for designing an effective thermal management system of modern electronic devices such as batteries, LEDs, microelectronics, ICs, circuit boards, semiconductors and heat spreaders. For having a precise thermal analysis, the temperature profile and thermal spreading resistance of the device should be calculated by considering the geometry, property and boundary conditions. Thermal spreading resistance occurs when heat enters through a portion of a surface and flows by conduction. It is the primary source of thermal resistance when heat flows from a tiny heat source to a thin and wide heat spreader. In this thesis, analytical models for modeling the temperature behavior and thermal resistance in some common geometries of microelectronic devices such as heat channels and heat tubes are investigated. Different boundary conditions for the system are considered. Along the source plane, a combination of discretely specified heat flux, specified temperatures and adiabatic condition are studied. Along the walls of the system, adiabatic or convective cooling boundary conditions are assumed. Along the sink plane, convective cooling with constant or variable heat transfer coefficient are considered. Also, the effect of orthotropic properties is discussed. This thesis contains nine chapters. Chapter one is the introduction and shows the concepts of thermal spreading resistance besides the originality and importance of the work. Chapter two reviews the literatures on the thermal spreading resistance in the past fifty years with a focus on the recent advances. In chapters three and four, thermal resistance of a twodimensional flux channel with non-uniform convection coefficient in the heat sink plane is studied. The non-uniform convection is modeled by using two functions than can simulate a wide variety of different heat sink configurations. In chapter five, a non-symmetrical flux channel with different heat transfer coefficient along the right and left edges and sink plane is analytically modeled. Due to the edge cooling and non-symmetry, the eigenvalues of the system are defined using the heat transfer coefficient on both edges and for satisfying the orthogonality condition, a normalized function is calculated. In chapter six, thermal behavior of two-dimensional rectangular flux channel with arbitrary boundary conditions on the source plane is presented. The boundary condition along the source plane can be a combination of the first kind boundary condition (Dirichlet or prescribed temperature) and the second kind boundary condition (Neumann or prescribed heat flux). The proposed solution can be used for modeling the flux channels with numerous different source plane boundary conditions without any limitations in the number and position of heat sources. In chapter seven, temperature profile of a circular flux tube with discretely specified boundary conditions along the source plane is presented. Also, the effect of orthotropic properties are discussed. In chapter 8, a three-dimensional rectangular flux channel with a non-uniform heat convection along the heat sink plane is analytically modeled. In chapter nine, a summary of the achievements is presented and some systems are proposed for the future studies. It is worth mentioning that all the models and case studies in the thesis are compared with the Finite Element Method (FEM).

Analytical thermal study on nonlinear fundamental heat transfer cases using a novel computational technique

In this paper a novel computational technique called Parameterized Perturbation Method (PPM) is used to obtain the solutions of nonlinear fundamental heat conduction equations. Three well known problems in the area of heat transfer are addressed to be solved. An analytical investigation is carried out for: (a) the temperature distribution in a fin with a temperature-dependent thermal conductivity, (b) the cooling of the lumped system with variable specific heat, and (c) the temperature distribution of a convective-radiative fin. The validity of the results of PPM solution was verified via comparison with numerical results obtained using a fourth order Runge-Kutta method. These comparisons revealed that PPM is a powerful approach for solving these problems. Also, the results showed that the main attributions of this method are very straightforward calculations and low computational burden compared to previous analytical and numerical approaches.

Defining the flammability of cylindrical metal rods through characterization of the thermal effects of the ignition promoter

All relevant international standards for determining if a metallic rod is flammable in oxygen utilize some form of “promoted ignition” test. In this test, for a given pressure, an overwhelming ignition source is coupled to the end of the test sample and the designation flammable or nonflammable is based upon the amount burned, that is, a burn criteria. It is documented that (1) the initial temperature of the test sample affects the burning of the test sample both (a) in regards to the pressure at which the sample will support burning (threshold pressure) and (b) the rate at which the sample is melted (regression rate of the melting interface); and, (2) the igniter used affects the test sample by heating it adjacent to the igniter as ignition occurs. Together, these facts make it necessary to ensure, if a metallic material is to be considered flammable at the conditions tested, that the burn criteria will exclude any region of the test sample that may have undergone preheating during the ignition process. A two-dimensional theoretical model was developed to describe the transient heat transfer occurring and resultant temperatures produced within this system. Several metals (copper, aluminum, iron, and stainless steel) and ignition promoters (magnesium, aluminum, and Pyrofuze®) were evaluated for a range of oxygen pressures between 0.69 MPa (100 psia) and 34.5 MPa (5,000 psia). A MATLAB® program was utilized to solve the developed model that was validated against (1) a published solution for a similar system and (2) against experimental data obtained during actual tests at the National Aeronautics and Space Administration White Sands Test Facility. The validated model successfully predicts temperatures within the test samples with agreement between model and experiment increasing as test pressure increases and/or distance from the promoter increases. Oxygen pressure and test sample thermal diffusivity were shown to have the largest effect on the results. In all cases evaluated, there is no significant preheating (above about 38°C/100°F) occurring at distances greater than 30 mm (1.18 in.) during the time the ignition source is attached to the test sample. This validates a distance of 30 mm (1.18 in.) above the ignition promoter as a burn length upon which a definition of flammable can be based for inclusion in relevant international standards (that is, burning past this length will always be independent of the ignition event for the ignition promoters considered here. KEYWORDS: promoted ignition, metal combustion, heat conduction, thin fin, promoted combustion, burn length, burn criteria, flammability, igniter effects, heat affected zone.

Crystal growth from undercooling melts under strong gradients of temperatures generating in short-duration microgravity

Crystal growth of bulk CdTe in short-duration microgravity is performed by the unidirectional cooling method. The largest growth grains in microgravity samples are 4X2mm. The cooling profiles indicate undercooling melts in microgravity. Cooling melt samples in microgravity generate strong gradient of temperature due to stop thermal convections. Temperature distribution in the melt is calculated by the one-dimensional equation of heat conduction, and about 100 K-undercooling is considered to occur at the cooling surface.

Natural convection flow in a porous enclosure with localized heating from below

Natural convection flow in a two-dimensional fluid saturated porous enclosure with localized heating from below, symmetrical cooling from the sides and the top and rest of the bottom walls are insulated, has been investigated numerically. Darcy’s law for porous media along with the energy equation based on the 1st law of thermodynamics has been considered. Implicit finite volume method with TDMA solver is used to solve the governing equations. Localized heating is simulated by a centrally located isothermal heat source on the bottom wall, and four different values of the dimensionless heat source length, 1/5, 2/5, 3/5 and 4/5 are considered. The effect of heat source length and the Rayleigh number on streamlines and isotherms are presented, as well as the variation of the local rate of heat transfer in terms of the local Nusselt number from the heated wall. Finally, the average Nusselt number at the heated part of the bottom wall has been shown against Rayleigh number for the non-dimensional heat source length.

Heating effects in nanofocusing metal wedges

This paper investigates theoretically and numerically local heating effects in plasmon nanofocusing structures with a particular focus on the sharp free-standing metal wedges. The developed model separates plasmon propagation in the wedge from the resultant heating effects. Therefore, this model is only applicable where the temperature increments in a nanofocusing structure are sufficiently small not to result in significant variations of the metal permittivity in the wedge. The problem is reduced to a one-dimensional heating model with a distributed heat source resulting from plasmon dissipation in the metal wedge. A simple heat conduction equation governing the local heating effects in a nanofocusing structure is derived and solved numerically for plasmonic pulses of different lengths and reasonable energies. Both the possibility of achieving substantial local temperature increments in the wedge (with a significant self-influence of the heating plasmonic pulses), and the possibility of relatively weak heating (to ensure the validity of the previously developed nanofocusing theory) are demonstrated and discussed, including the future applications of the obtained results. Applicability conditions for the developed model are also derived and discussed.

The hyperthermic patient

Concepts used in this chapter include: Thermoregulation:- Thermoregulation refers to the body’s sophisticated, multi-system regulation of core body temperature. This hierarchical system extends from highly thermo-sensitive neurons in the preoptic region of the brain proximate to the rostral hypothalamus, down to the brain stem and spinal cord. Coupled with receptors in the skin and spine, both central and peripheral information on body temperature is integrated to inform and activate the homeostatic mechanisms which maintain our core temperature at 37oC1. Hyperthermia:- An imbalance between the metabolic and external heat accumulated in the body and the loss of heat from the body2. Exertional heat stroke:- A disorder of excessive heat production coupled with insufficient heat dissipation which occurs in un-acclimated individuals who are engaging in over-exertion in hot and humid conditions. This phenomenon includes central nervous system dysfunction and critical dysfunction to all organ systems including renal, cardiovascular, musculoskeletal and hepatic functions. Non-exertional heat stroke:- In contrast to exertional heatstroke as a consequence of high heat production during strenuous exercise, non-exertional heatstroke results from prolonged exposure to high ambient temperature. The elderly, those with chronic health conditions and children are particularly susceptible.3 Rhabdomylosis:- An acute, sometimes fatal disease characterised by destruction of skeletal muscle. In exertional heat stroke, rhabdomylosis occurs in the context of strenuous exercise when mechanical and/or metabolic stress damages the skeletal muscle, causing elevated serum creatine kinease. Associated with this is the potential development of hyperkalemia, myoglobinuria and renal failure. Malignant hyperthermia:- Malignant hyperthermia is “an inherited subclinical myopathy characterised by a hypermetabolic reaction during anaesthesia. The reaction is related to skeletal muscle calcium dysregulation triggered by volatile inhaled anaesthetics and/or succinylcholine.”4 Presentation includes skeletal muscle rigidity, mixed metabolic and respiratory acidosis, tachycardia, hyperpyrexia, rhabdomylosis, hyperkalaemia, elevated serum creatine kinease, multi-organ failure, disseminated intravascular coagulation and death.5

Catalytic probes with nanostructured surface for gas/discharge diagnostics : a study of a probe signal behaviour

Catalytic probes are used for plasma diagnostics in order to quantify the density of neutral atoms. The probe response primarily depends on the probe material and its surface morphology. Here we report on the design, operation and modelling of the response of niobium pentoxide sensors with a flat and nanowire (NW) surfaces. These sensors were used to detect neutral oxygen atoms in the afterglow region of an inductively coupled rf discharge in oxygen. A very different response of the flat-surface and NW probes to the varying densities of oxygen atoms was explained by modelling heat conduction and taking into account the associated temperature gradients. It was found that the nanostructure probe can measure in a broader range than the flat oxide probe due to an increase in the surface to volume ratio, and the presence of nanostructures which act as a thermal barrier against sensor overheating. These results can be used for the development of the new generation of catalytic probes for gas/discharge diagnostics in a range of industrial and environmental applications.

Observation of a thermally reversing window in bulk Ge15Te85-xInx glasses

Alternating differential scanning calorimetry measurements have been undertaken on the Ge15Te85-xInx (1 <= x <= 11) series of glasses. It is found that there is a marginal decrease in the glass transition temperature (T-g) in the composition range 1 <= x <= 3. Above x = 3, a monotonic increase is seen in T-g which indicates a continuous build-up in network connectivity and absence of any nanophase separation. The non-reversing heat flow (Delta H-NR) has been found to exhibit a broad trough between the compositions x = 3 and 7, which clearly indicates the presence of a thermally reversing window in Ge15Te85-xInx glasses in the composition range 3 <= x <= 7.

Steady incompressible flow of cohesionless granular materials through a wedge-shaped hopper: Frictional 14kinetic solution to the smooth wall, radial gravity problem

Hybrid frictional-kinetic equations are used to predict the velocity, grain temperature, and stress fields in hoppers. A suitable choice of dimensionless variables permits the pseudo-thermal energy balance to be decoupled from the momentum balance. These balances contain a small parameter, which is analogous to a reciprocal Reynolds number. Hence an approximate semi-analytical solution is constructed using perturbation methods. The energy balance is solved using the method of matched asymptotic expansions. The effect of heat conduction is confined to a very thin boundary layer near the exit, where it causes a marginal change in the temperature. Outside this layer, the temperature T increases rapidly as the radial coordinate r decreases. In particular, the conduction-free energy balance yields an asymptotic solution, valid for small values of r, of the form T proportional r-4. There is a corresponding increase in the kinetic stresses, which attain their maximum values at the hopper exit. The momentum balance is solved by a regular perturbation method. The contribution of the kinetic stresses is important only in a small region near the exit, where the frictional stresses tend to zero. Therefore, the discharge rate is only about 2.3% lower than the frictional value, for typical parameter values. As in the frictional case, the discharge rate for deep hoppers is found to be independent of the head of material.

Decoupled three-dimensional finite element computation of thermoelastic damping using Zener's approximation

We consider three dimensional finite element computations of thermoelastic damping ratios of arbitrary bodies using Zener's approach. In our small-damping formulation, unlike existing fully coupled formulations, the calculation is split into three smaller parts. Of these, the first sub-calculation involves routine undamped modal analysis using ANSYS. The second sub-calculation takes the mode shape, and solves on the same mesh a periodic heat conduction problem. Finally, the damping coefficient is a volume integral, evaluated elementwise. In the only other decoupled three dimensional computation of thermoelastic damping reported in the literature, the heat conduction problem is solved much less efficiently, using a modal expansion. We provide numerical examples using some beam-like geometries, for which Zener's and similar formulas are valid. Among these we examine tapered beams, including the limiting case of a sharp tip. The latter's higher-mode damping ratios dramatically exceed those of a comparable uniform beam.

Estimation of surface temperature from MONTBLEX data

It is observed that the daily mean temperature of the soil is linear with depth and the variation of the temperature is sinusoidal with a period of a day. Based on these observations the one-dimensional heat conduction equation for the soil can be solved which gives the amplitude and phase variation of the temperature wave with depth. Given the temperature data at three levels below the surface, the amplitude and phase variation and hence the surface temperature variation over the day are estimated. The daily mean temperature of the surface is estimated from linear extrapolation of the daily means at the three levels below the surface. Estimated values of soil thermal diffusivity show a subtantial change after sudden and heavy rains.

Oscillatory flow and temperature fields in an open tube with temperature difference across the ends

The oscillating flow and temperature field in an open tube subjected to cryogenic temperature at the cold end and ambient temperature at the hot end is studied numerically. The flow is driven by a time-wise sinusoidally varying pressure at the cold end. The conjugate problem takes into account the interaction of oscillatory flow with the heat conduction in the tube wall. The full set of compressible flow equations with axisymmetry assumption are solved with a pressure correction algorithm. Parametric studies are conducted with frequencies of 5-15 Hz, with one end maintained at 100 K and other end at 300 K. The flow and temperature distributions and the cooldown characteristics are obtained. The frequency and pressure amplitude have negligible effect on the time averaged Nusselt number. Pressure amplitude is an important factor determining the enthalpy flow through the solid wall. The frequency of operation has considerable effect on penetration of temperature into the tube. The density variation has strong influence on property profiles during cooldown. The present study is expected to be of interest in applications such as pulse tube refrigerators and other cryocoolers, where oscillatory flows occur in open tubes. (C) 2011 Elsevier Ltd. All rights reserved.