A rigorous derivation of the bioheat equation for local tissue heat transfer based on a volume averaging theory

A general three-dimensional bioheat equation for local tissue heat transfer has been derived with less assumptions, exploiting a volume averaging theory commonly used in fluid-saturated porous media. The volume averaged energy equations obtained for the arterial blood, venous blood and tissue were combined together to form a single energy equation in terms of the tissue temperature alone. The resulting energy equation turns out to be remarkably simple as we define the effective thermal conductivity tensor, which accounts not only for the countercurrent heat exchange mechanism but also for the thermal dispersion mechanism. The present equation for local tissue heat transfer naturally reduces to the Weinbaum-Jiji equation for the unidirectional case.     

Vasculature based model for characterizing the oxygen transport in skin tissues – analogy to the Weinbaum-Jiji bioheat equation

Based on the conceptual three-layer microvascular structure of skin tissues proposed by Weinbaum et al. [20–25] and in analogy to the well known Weinbaum-Jiji (W-J) bioheat equation, a new oxygen transport model was established in this paper, which collectively included the contributions of the vascular geometry and the blood flow condition. The new one-dimensional three-layer oxygen transport model was then applied to predict the average oxygen concentration distribution in skin tissues and numerical solutions for the boundary value problem coupling the three layers were obtained. A simple expression for the tensor diffusivity (D_eff) of oxygen transport over the deep tissue layer was presented, which was orders of magnitude higher than the intrinsic diffusivity (D_t) in tissue without blood flow. Effects of blood flow velocity and vascular geometry to the oxygen transport were investigated. Calculations indicated that the vascular geometry had significant effects on oxygen transport. The oxygen exchange between the arteries and veins was relatively small for the deep tissue layer. Further, the average oxygen concentration gradient appears low in intermediate layer due to large capillary perfusion. The theoretical results were implemented to interpret some previous experimental results and a better understanding on the oxygen transport across the vascularized living tissues was obtained. The strategy proposed in this paper may provide a feasible way to comprehensively characterize the oxygen transport behaviors in living tissues with real and complex vasculature.     

Thermodynamical consistency of the dual-phase-lag heat conduction equation

Dual-phase-lag equation for heat conduction is analyzed from the point of view of non-equilibrium thermodynamics. Its first-order Taylor series expansion is consistent with the second law as long as the two relaxation times are not negative.     

Analysis of heating a three-dimensional porous bed utilizing the two energy equation model

In this paper heating a three-dimensional porous packed bed by a non-thermal equilibrium flow of incompressible fluid is analytically investigated. A two energy equation model is employed to simulate the temperature difference between the fluid and solid phases. Using the perturbation technique, an analytical solution for the problem is obtained. It is shown that the temperature difference between the fluid and solid phases forms a wave localized in space and propagating from the fluid inlet boundary in the direction of the flow. The amplitude of the wave decreases while the wave propagates downstream.     

Analytical solution of the poiseuille flow problem using the ellipsoidal-statistical model of the Boltzmann kinetic equation

An analytical solution (in the form of a Neumann series) of the problem of rarefied gas flow in a plane channel with infinite walls in the presence of a pressure gradient (Poiseuille flow) parallel to them is constructed within the framework of the kinetic approach in an isothermal approximation. The ellipsoidal-statistical model of the Boltzmann kinetic equation and the diffuse reflection model are used as the basic equation and the boundary condition, respectively. Using the resulting distribution function, the mass and heat flux densities in the direction of the pressure gradient per unit channel length in the y′ direction are calculated, and profiles of the gas mass velocity and heat flux in the channel are constructed. The results obtained for the continuum and free-molecular flow models are analyzed and compared with similar results obtained by numerical methods.     

Heat transfer mechanisms during short-pulse laser heating of two-layer composite thin films

 Using a simple perturbation technique, an analytical investigation is presented for the heat transfer mechanisms during ultrafast laser heating of two-layer composite thin slabs from a microscopic point of view. The composite slab consists of two thin metal films which are in perfect thermal contact. The microscopic parabolic two-step model is adopted to describe the behavior of the composite slab. In the microscopic two-step model, the heating process is modeled by the deposition of radiation energy on electrons, the transport of energy by electrons, and the heating of the material lattice through electron–phonon interactions. The proposed perturbation technique is used when the normalized temperature difference between the solid lattice and the electron gas is relatively a small perturbed quantity. Received on 20 September 2000 / Published online: 29 November 2001     

Thermal stability of composite superconducting tape under the effect of a two-dimensional dual-phase-lag heat conduction model

Thermal stability of composite superconducting tape subjected to a thermal disturbance is numerically investigated under the effect of a two-dimensional dual-phase-lag heat conduction model. It is found that the dual-phase-lag model predicts a wider stable region as compared to the predictions of the parabolic and the hyperbolic heat conduction models. The effects of different design, geometrical and operating conditions on superconducting tape thermal stability were also studied.a conductor width, (m) - A conductor cross sectional area of, (m²) - As conductor aspect ratio, (a/b) - b conductor thickness, (m) - Bi Biot number - B dimensionless disturbance Intensity - C heat capacity, (J m^–3 K^–1) - D disturbance energy density, (W m^–3) - f volume fraction of the stabilizer in the conductor - g(T) steady capacity of the Ohmic heat source, (W m^–3) - g_max Ohmic heat generation with the whole current in the stabilizer, (W m^–3) - G_max dimensionless maximum Joule heating - h convective heat transfer coefficient, (W m^–2 K^–1) - J current density, (A m^–2) - k thermal conductivity of conductor, (W m^–1 K^–1) - q conduction heat flux vector, (W m^–2) - Q dimensionless Joule heating - R relaxation times ratio (_T/2_q) - t rime, (s) - T temperature, (K) - T_c critical temperature, (K) - T_c1 current sharing temperature, (K) - T_i initial temperature, (K) - T_o ambient temperature, (K) - x, y co-ordinate defined in Fig. 1, (m) - thermal diffusivity (m² s^–1) - dimensionless time - _i dimensionless duration time - dimensionless y-variable - _o superconductor dimensionless thickness - dimensionless temperature - _c1 dimensionless current sharing temperature - 1 dimensionless maximum temperature - dimensionless disturbance energy - numerical tolerance - x width of conductor subjected to heat disturbances, (m) - y thickness of conductor subjected to heat disturbances, (m) - dimensionless x-variable - _o superconductor dimensionless width - stabilizer electrical resistivity, () - _q relaxation time of heat flux, (s) - _T relaxation time of temperature gradient, (s) - i initial - sc current sharing - max maximum - o ambient     

Finite element analysis of conductive and radiative heating of a thin skin calorimeter

The main purpose of the present work is to investigate the influence of normal and lateral conduction on the temperature distribution and heat transfer coefficient on the surface of a typical sounding rocket. A two-dimensional heat conduction equation with a time dependent aerodynamic heating condition at one surface and a radiation boundary condition at the other end is solved using finite element method.     

Heat transfer during transient cooling of high temperature surface with an impinging jet

 An experimental study of transient boiling heat transfer during a cooling of a hot cylindrical block with an impinging water jet has been made at atmospheric pressure. The experimental data were taken for the following conditions: a degree of subcooling of ΔT _sub = 20–80 K, a jet velocity of u _j  = 5–15 m/s, a nozzle diameter of d _j  = 2 mm and three materials of copper, brass and carbon steel. The block was initially and uniformly heated to about 250 °C and the transient temperatures in the block were measured at eight locations in r-direction at two different depths from the surface during the cooling of hot block. The surface heat flux distribution with time was evaluated using a numerical analysis of 2-D heat conduction. Behavior of the wetting front, which is extending the nucleate boiling region outward, is observed with a high-speed video camera. A position of wetting region is measured and it is correlated well with a power function of time. The changes in estimated heat flux and temperature were compared with the position of wetting region to clarify the effects of subcooling, jet velocity and thermal properties of block on the transient cooling. Received on 17 March 2000     

Analytical study of the transition curves in the bi-linear Mathieu equation

Nonlinear Dynamics - The current work is primarily devoted to the asymptotic analysis of the instability zones existing in the bi-linear Mathieu equation. In this study, we invoke the common...     

Two-fluid model for transient analysis of slug flow in oil wells

In this work it is presented a transient, one-dimensional, adiabatic model for slug flow simulation, which appears when liquid (mixture of oil and water) and gas flow simultaneously through pipes. The model is formed by space and time averaged conservation equations for mass, momentum and energy for each phase, the numerical solution is based on the finite difference technique in the implicit scheme. Velocity, pressure, volumetric fraction and temperature profiles for both phases were predicted for inclination angles from the horizontal to the vertical position (unified model) and ascendant flow. Predictions from the model were validated using field data and ten correlations commonly used in the oil industry. The effects of gas heating or cooling, due to compression and expansion processes, on the predictions and numerical stability, were studied. It was found that when these effects are taken into account, a good behavior of temperature predictions and numerical stability are obtained. The model presents deviations lower than 14% regarding field data and it presents better predictions than most of the correlations.     

Thermal analysis of thin multi-layer metal films during femtosecond laser heating

Multi-layer metals films are widely used in modern engineering applications such as gold-coated metal mirrors used in high power laser systems. A transient heat flux model is derived to analyze multi-layer metal films under laser heating. The two separate system composed of electrons and the lattice is considered to take into account the electron–lattice interaction. The present model predicted the effects of underlying chromium’s thermal properties on temperature rise of the top gold layer. The effects of two adjacent and different metals with different electron–lattice coupling factors are analyzed for the heating mechanism of different lattices. The derived transient model combined with the two different conservation equations for the lattice and electrons are applied for the ultra short-pulse laser heating of a multi-layer film composed of gold and chromium.     

Moisture transfer analysis during drying of slab woods

This article presents an analytical technique for determining the moisture diffusivities and moisture transfer coefficients for slab shaped woods subjects to drying process. The analysis of transient moisture diffusion is carried out on the basis of two important practical criteria: 0.1<Bi<100 and Bi>100. The drying coefficients and lag factors are defined for wood-drying applications and incorporated into the models. In order to verify the present models, the model results are compared with experimental measurements taken from the literature and good agreement was found. Results show that the technique presented here is capable of determining the moisture diffusivities and moisture transfer coefficients for slab woods in a simple and accurate manner for practical applications and will be beneficial to the relevant wood␣drying industries. This approach can be extended to␣different wood products of regular and irregular shapes. Received on 23 January 1998     

Numerical analysis of the transient conjugated heat transfer in a circular duct with a power-law fluid

We treat numerically in this paper, the transient analysis of a conjugated heat transfer process in the thermal entrance region of a circular tube with a fully developed laminar power-law fluid flow. We apply the quasi-steady approximation for the power-law fluid, identifying the suitable time scales of the process. Thus, the energy equation in the fluids is solved analytically using the well-known integral boundary layer technique. This solution is coupled to the transient energy equation for the solid where the transverse and longitudinal heat conduction effects are taken into account. The numerical results for the temporal evolution of the average temperature of the tube wall, _av, is plotted for different nondimensional parameters such as conduction parameter, , the aspect ratios of the tube, and ₀ and the index of power-law fluid, n.     

Use of pseudo-concentrations to follow creeping viscous flows during transient analysis

Creeping viscous flows are followed through finite element meshes by use of pseudo-concentrations which define material position. The concentrations, assumed to be transported only by convection, serve as material markers. Illustrations are presented related to industrial forming processes and the slow deformation of geological structures.     

Two dimensional analysis of transient flow and heat transfer in a natural circulation loop

The transient heat transfer, fluid flow and pressure in a natural circulation loop have been studied under laminar flow conditions. Most studies of these systems have utilized a onedimensional approach which requires a priori specifications of the friction and the heat-transfer coefficients. In the present work the variation of the friction and heat-transfer coefficients are determined. Detailed pressure, temperature and velocity distributions are presented.     

The generalized Hamiltonian model for the shafting transient analysis of the hydro turbine generating sets

Traditional rotor dynamics mainly focuses on the steady- state behavior of the rotor and shafting. However, for systems such as hydro turbine generating sets (HTGS) where the control and regulation is frequently applied, the shafting safety and stabilization in transient state is then a key factor. The shafting transient state inevitably involves multiparameter domain, multifield coupling, and coupling dynamics. In this paper, the relative value form of the Lagrange function and its equations have been established by defining the base value system of the shafting. Taking the rotation angle and the angular speed of the shafting as a link, the shafting lateral vibration and generator equations are integrated into the framework of generalized Hamiltonian system. The generalized Hamiltonian control model is thus established. To make the model more general, additional forces of the shafting are taken as the input excitation in proposed model. The control system of the HTGS can be easily connected with the shafting model to form the whole simulation system of the HTGS. It is expected that this study will build a foundation for the coupling dynamics theory using the generalized Hamiltonian theory to investigate coupling dynamic mechanism among the shafting vibration, transient of hydro turbine generating sets, and additional forces of the shafting.