Heat generation in laser irradiated tissue

Many medical applications involving lasers rely upon the generation of heat within the tissue for the desired therapeutic effect. Determination of the absorbed light energy in tissue is difficult in many cases. Although UV wavelengths of the excimer laser and 10.6 microns wavelength of the CO2 laser are absorbed within the first 20 microns of soft tissue, visible and near infrared wavelengths are scattered as well as absorbed. Typically, multiple scattering is a significant factor in the distribution of light in tissue and the resulting heat source term. An improved model is presented for estimating heat generation due to the absorption of a collimated (axisymmetric) laser beam and scattered light at each point r and z in tissue. Heat generated within tissue is a function of the laser power, the shape and size of the incident beam and the optical properties of the tissue at the irradiation wavelength. Key to the calculation of heat source strength is accurate estimation of the light distribution. Methods for experimentally determining the optical parameters of tissue are discussed in the context of the improved model.     

Estimation of tumor characteristics in a breast tissue with known skin surface temperature

Presence of a tumor and its characteristics like location, size and properties are estimated. Estimation is based on the measurement of the skin surface temperature of the breast. Consideration is given to a 2-D breast tissue infected with a tumor. Heat transfer in the breast tissue modeled using the Pennes bioheat equation is solved by the finite volume method. Skin surface temperature profile of the breast is characteristic of the tumor location, its size and grade. In the inverse analysis, the objective function is minimized using the genetic algorithm. Exact estimation is obtained if one parameter is estimated at a time. However, the accuracies are acceptable even when blood perfusion rate, location and size are estimated simultaneously.     

FDTD analysis of a noninvasive hyperthermia system for brain tumors

ABSTRACT: BACKGROUND: Hyperthermia is considered one of the new therapeutic modalities for cancer treatment and is based on the difference in thermal sensitivity between healthy tissues and tumors. During hyperthermia treatment, the temperature of the tumor is raised to 40--4[DEGREE SIGN]C for a definite period resulting in the destruction of cancer cells. This paper investigates design, modeling and simulation of a new non-invasive hyperthermia applicator system capable of effectively heating deep seated as well as superficial brain tumors using inexpensive, simple, and easy to fabricate components without harming surrounding healthy brain tissues. METHODS: The proposed hyperthermia applicator system is composed of an air filled partial half ellipsoidal chamber, a patch antenna, and a head model with an embedded tumor at an arbitrary location. The irradiating antenna is placed at one of the foci of the hyperthermia chamber while the center of the brain tumor is placed at the other focus. The finite difference time domain (FDTD) method is used to compute both the SAR patterns and the temperature distribution in three different head models due to two different patch antennas at a frequency of 915 MHz. RESULTS: The obtained results suggest that by using the proposed noninvasive hyperthermia system it is feasible to achieve sufficient and focused energy deposition and temperature rise to therapeutic values in deep seated as well as superficial brain tumors without harming surrounding healthy tissue. CONCLUSIONS: The proposed noninvasive hyperthermia system proved suitable for raising the temperature in tumors embedded in the brain to therapeutic values by carefully selecting the systems components. The operator of the system only needs to place the center of the brain tumor at a pre-specified location and excite the antenna at a single frequency of 915 MHz. Our study may provide a basis for a clinical applicator prototype capable of heating brain tumors.     

Procedure to estimate thermophysical and geometrical parameters of embedded cancerous lesions using thermography

Based on the fact that malignant cancerous lesions (neoplasms) develop high metabolism and use more blood supply than normal tissue, infrared thermography (IR) has become a reliable clinical technique used to indicate noninvasively the presence of cancerous diseases, e.g., skin and breast cancer. However, to diagnose cancerous diseases by IR, the technique requires procedures that explore the relationship between the neoplasm characteristics (size, blood perfusion rate and heat generated) and the resulting temperature distribution on the skin surface. In this research work the dual reciprocity boundary element method (DRBEM) has been coupled with the simulated annealing technique (SA) in a new inverse procedure, which coupled to the IR technique, is capable of estimating simultaneously geometrical and thermophysical parameters of the neoplasm. The method is of an evolutionary type, requiring random initial values for the unknown parameters and no calculations of sensitivities or search directions. In addition, the DRBEM does not require any re-meshing at each proposed solution to solve the bioheat model. The inverse procedure has been tested considering input data for simulated neoplasms of different sizes and positions in relation to the skin surface. The successful estimation of unknown neoplasm parameters validates the idea of using the SA technique and the DRBEM in the estimation of parameters. Other estimation techniques, based on genetic algorithms or sensitivity coefficients, have not been capable of obtaining a solution because the skin surface temperature difference is very small.     

Simulation and investigation of quantum dot effects as internal heat-generator source in breast tumor site

The Pennes bio-heat transfer equation, which introduces the exchange magnitude of heat transfer between tissue and blood, is often used to solve the temperature distribution for thermal imaging and sensing. Near-infrared light has the ability to be used as a non-invasive means of diagnostic imaging within the woman's breast. Due to the diffusive nature of light in different tissue, computational model-based methods are required for functional imaging within the breast. In this article, the time-dependent bio-heat transfer is solved by a numerical method. In our model, the heat generation source (intrinsic and extrinsic) involves laser, metabolism, and quantum dot that the metabolism and heat generated by QDs are considered as intrinsic. We supposed the injected quantum dots would target the tumor site by a passive targeting process and then by interaction of laser radiation and quantum dot, the photoluminescence of quantum dot is converted to heat in the tumor site. The extra generated heat can impact on the extracted heat profile. One of the important applications of this research has led to a sensitivity improvement of the imaging system, which is potentially useful in the diagnosis and detection of breast cancer.     

内热源影响下的均质组织表面温度分布研究

现有的医用红外热像技术几乎都采自体表,为了提高诊断准确性,必须进一步确定体表温度分布与体内热源间的关系。基于人体内部传热机制和人体与环境的传热机制,推导出稳态下均质（各向同性）组织的表面温度分布与内部球状热源之间的关系,通过模拟计算可验证结论的准确性。     

Locating multiple tumors by moving shape analysis

A methodology to determine the unknown shape of an embedded tumor is proposed. A functional that represents the mismatch between a measured experimental temperature profile, which may be obtained by infrared thermography at skin surface, and the solution of an appropriate boundary problem is defined. Using the Pennes’s bioheat transfer equations, the temperature in a section of healthy tissue with a tumor region is modeled by a boundary problem. The functional is related to the shape of the tumor through the solution of the boundary problem, in such a way that finding the minimum of the functional form also means finding the unknown shape of the embedded tumor. The shape derivative of the functional is computed in each node of an approximation of the solution by the method of Finite Elements using similar methods considered by Pironneau [7]. The algorithm presented include an adaptive strategy to improve the error of the objective function. Numerical results with multiple connected tumors are considered to illustrate the potential of the proposed methodology.     

Numerical analysis for determination of the presence of a tumor and estimation of its size and location in a tissue

This article deals with the numerical analysis to ascertain the presence of a tumor and to estimate its size and location in a tissue. Heat transfer in the tissue is modeled using the Pennes bioheat transfer equation, and is solved using the finite volume method. Consideration is given to 1-D brain and breast tissues. Temperature distributions in the tissues are specific to the tumor grades, its locations and sizes, and these are different than that of a normal tissue. With temperature distribution known a priori, estimations of the position and the size of a tumor are done using the inverse analysis. The proposed approach gives a correct estimation of the presence of a tumor and its location and size.     

Numerical study of the effect of ultrasound frequency on temperature distribution in layered tissue

The high intensity focused ultrasound (HIFU) technology can produce therapeutic benefits in deep-seated tissues of interest, selectively and noninvasively. In order to control the treatment process, it is important to recognize the heat generation in biological tissue and the parameters that have an effect on temperature rising. This study investigates the influence of frequency and source intensity on temperature distribution during high-intensity focused ultrasound (HIFU). A nonlinear full wave equation model is simulated to compute the pressure field. Additionally, the absorbed coefficient of tissue is added to the nonlinear equations to simulate accurately the wave propagation in tissue with high absorbed coefficient. In addition, temperature distribution was solved by the Pennes bio-heat equation. Conclusively, frequencies in the range of 1–1.5 MHz are prescribed to have maximum heat absorption in the focal region.     

Radio-frequency ablation in a realistic reconstructed hepatic tissue

This study uses a reconstructed vascular geometry to evaluate the thermal response of tissue during a three-dimensional radiofrequency (rf) tumor ablation. MRI images of a sectioned liver tissue containing arterial vessels are processed and converted into a finite-element mesh. A rf heat source in the form of a spherically symmetric Gaussian distribution, fit from a previously computed profile, is employed. Convective cooling within large blood vessels is treated using direct physical modeling of the heat and momentum transfer within the vessel. Calculations of temperature rise and thermal dose are performed for transient rf procedures in cases where the tumor is located at three different locations near the bifurcation point of a reconstructed artery. Results demonstrate a significant dependence of tissue temperature profile on the reconstructed vasculature and the tumor location. Heat convection through the arteries reduced the steady-state temperature rise, relative to the no-flow case, by up to 70% in the targeted volume. Blood flow also reduced the thermal dose value, which quantifies the extent of cell damage, from approximately 3600 min, for the no-flow condition, to 10 min for basal flow (13.8 cms). Reduction of thermal dose below the threshold value of 240 min indicates ablation procedures that may inadequately elevate the temperature in some regions, thereby permitting possible tumor recursion. These variations are caused by vasculature tortuosity that are patient specific and can be captured only by the reconstruction of the realistic geometry.     

Heat transfer analysis of skin during thermal therapy using thermal wave equation

Specifying exact geometry of vessel network and its effect on temperature distribution in living tissues is one of the most complicated problems of the bioheat field. In this paper, the effects of blood vessels on temperature distribution in a skin tissue subjected to various thermal therapy conditions are investigated. Present model consists of counter-current multilevel vessel network embedded in a three-dimensional triple-layered skin structure. Branching angles of vessels are calculated using the physiological principle of minimum work. Length and diameter ratios are specified using length doubling rule and Cube law, respectively. By solving continuity, momentum and energy equations for blood flow and Pennes and modified Pennes bioheat equations for the tissue, temperature distributions in the tissue are measured. Effects of considering modified Pennes bioheat equation are investigated, comprehensively. It is also observed that blood has an impressive role in temperature distribution of the tissue, especially at high temperatures. The effects of different parameters such as boundary conditions, relaxation time, thermal properties of skin, metabolism and pulse heat flux on temperature distribution are investigated. Tremendous effect of boundary condition type at the lower boundary is noted. It seems that neither insulation nor constant temperature at this boundary can completely describe the real physical phenomena. It is expected that real temperature at the lower levels is somewhat between two predicted values. The effect of temperature on the thermal properties of skin tissue is considered. It is shown that considering temperature dependent values for thermal conductivity is important in the temperature distribution estimation of skin tissue; however, the effect of temperature dependent values for specific heat capacity is negligible. It is seen that considering modified Pennes equation in processes with high heat flux during low times is significant.     

The heat transfer analysis of nanoparticle heat source in alanine tissue by molecular dynamics

The purpose of this study is to simulate the heat transfer problem when the 3-D Alanine tissue is heated by the gold nanoparticle in the field of molecular dynamics. In this paper, the Alanine molecule is adopted and its parameters are available in the GROMACS protein data bank. A computing algorithm is developed to evaluate the heat transfer phenomena in the nano-scale biological system based on the molecular dynamics and the protein data bank. The value of the thermal conductivity of Alanine is calculated from the autocorrelation function of the Green-Kubo formalism and this result has a roughly approximation with the bulk thermal conductivity reported by experimental data . Two kinds of problems are investigated in the paper. One is the Alanine tissue heated by the constant heat source and the other is by the time-varying heat source. The numerical results show that a temperature jump exists around the source and the temperature profiles drop to the environmental temperature within a very short distance. It concludes that only a small region around the nano-scale heat source is affected by the heated process. Therefore, the results of the nanoparticle-heated method could be applied to the clinical therapy of tumor, and the normal cells are destroyed only within a smaller region than those of chemotherapy or surgery.     

3D in silico study of magnetic fluid hyperthermia of breast tumor using Fe3O4 magnetic nanoparticles

Modeling and simulation of the temperature distribution, the mass concentration, and the heat transfer in the breast tissue are hot issues in magnetic fluid hyperthermia treatment of cancer. The breast tissue can be visualized as a porous matrix with saturated blood. In this paper, 3D in silico study of breast cancer hyperthermia using magnetic nanoparticles (MNPs) is conducted. The 3D FEM models are incorporated to investigate the infusion and backflow of nanofluid in the breast tumor, the diffusion of nanofluid, temperature distribution during the treatment, and prediction of the fraction of tumor necrosis while dealing with the thermal therapy. All the hyperthermia procedures are simulated and analyzed on COMSOL Multiphysics. The sensitivity of frequency and amplitude of the applied magnetic field (AMF) is investigated on the heating effect of the tumor. The mesh dependent solution of Penne's bioheat model is also analyzed. The simulated results demonstrate successful breast cancer treatment using MNPs with minimum side effects. Validation of current simulations results with experimental studies existing in literature advocates the success of our therapy. The increase in the amplitude and frequency of the AMF increases of the temperature in the tumor. The variation of mesh from coarser to finer increased the temperature through small fractions. We have also simulated the magnetic induction problem where the magnetic field is generated by current-carrying coil conductors induce heat in nearby breast tumors due to excitation of MNPs by magnetic flux. This research will aid treatment protocols and real-time clinical breast cancer treatments.     

A self-heated thermistor technique to measure effective thermal properties from the tissue surface

A microcomputer based instrument to measure effective thermal conductivity and diffusivity at the surface of a tissue has been developed. Self-heated spherical thermistors, partially embedded in an insulator, are used to simultaneously heat tissue and measure the resulting temperature rise. The temperature increase of the thermistor for a given applied power is a function of the combined thermal properties of the insulator, the thermistor, and the tissue. Once the probe is calibrated, the instrument accurately measures the thermal properties of tissue. Conductivity measurements are accurate to 2 percent and diffusivity measurements are accurate to 4 percent. A simplified bioheat equation is used which assumes the effective tissue thermal conductivity is a linear function of perfusion. Since tissue blood flow strongly affects heat transfer, the surface thermistor probe is quite sensitive to perfusion.     

Heat transfer normal to paired arterioles and venules embedded in perfused tissue during hyperthermia

A numerical model of the heat transer normal to an arteriole-venule pair embedded in muscle tissue has been constructed. Anatomical data describing the blood vessel size, spacing, and density have been incorporated into the model. This model computes temperatures along the vessel walls as well as the temperature throughout the tissue which comprises an infinitely long Krogh cylinder around the vessel pair. Tissue temperatures were computed in the steady-state under resting conditions, while transient calculations were made under hyperthermic conditions. Results show that for both large- (1st generation) and medium-sized (5th generation) vessel pairs, the mean tissue temperature within the tissue cylinder is not equal to the mean of the arteriole and venule blood temperatures under both steady-state and transient conditions. The numerical data were reduced so that a comparison could be made with the predictions of a simple two-dimensional superposition of line sources and sinks presented by Baish et al. This comparison reveals that the superposition model accurately describes the heat transfer effects during hyperthermia, permitting subsequent incorporation of this theory into a realistic three-dimensional model of heat transfer in a whole limb during hyperthermia.     

A review on numerical modeling for magnetic nanoparticle hyperthermia: Progress and challenges

Recent progress in nanotechnology has advanced the development of magnetic nanoparticle (MNP) hyperthermia as a potential therapeutic platform for treating diseases. Due to the challenges in reliably predicting the spatiotemporal distribution of temperature in the living tissue during the therapy of MNP hyperthermia, critical for ensuring the safety as well as efficacy of the therapy, the development of effective and reliable numerical models is warranted. This article provides a comprehensive review on the various mathematical methods for determining specific loss power (SLP), a parameter used to quantify the heat generation capability of MNPs, as well as bio-heat models for predicting heat transfer phenomena and temperature distribution in living tissue upon the application of MNP hyperthermia. This article also discusses potential applications of the bio-heat models of MNP hyperthermia for therapeutic purposes, particularly for cancer treatment, along with their limitations that could be overcome.     

Estimation of surface heat flux and temperature distributions in a multilayer tissue based on the hyperbolic model of heat conduction

In this study, an inverse algorithm based on the conjugate gradient method and the discrepancy principle is applied to solve the inverse hyperbolic heat conduction problem in estimating the unknown time-dependent surface heat flux in a skin tissue, which is stratified into epidermis, dermis, and subcutaneous layers, from the temperature measurements taken within the medium. Subsequently, the temperature distributions in the tissue can be calculated as well. The concept of finite heat propagation velocity is applied to the modeling of the bioheat transfer problem. The inverse solutions will be justified based on the numerical experiments in which two different heat flux distributions are to be determined. The temperature data obtained from the direct problem are used to simulate the temperature measurements. The influence of measurement errors on the precision of the estimated results is also investigated. Results show that an excellent estimation on the time-dependent surface heat flux can be obtained for the test cases considered in this study.     

A three-dimensional variable geometry countercurrent model for whole limb heat transfer.

A new formulation of the combined macro and microvascular model for heat transfer in a human arm developed in Song et al. [1] is proposed using a recently developed approximate theory for the heat exchange between countercurrent vessels embedded in a tissue cylinder with surface convection [2]. The latter theory is generalized herein to treat an arm with an arbitrary variation in cross-sectional area and continuous bleed off from the axial vessels to the muscle and cutaneous tissue. The local microvascular temperature field is described by a "hybrid" model which applies the Weinbaum-Jiji [3] and Pennes [4] equations in the peripheral and deeper tissue layers, respectively. To obtain reliable end conditions at the wrist and other model input parameters, a plethysmograph-calorimeter has been used to measure the blood flow distribution between the arm and hand circulations, and hand heat loss. The predictions of the model show good agreement with measurements for the axial surface temperature distribution in the arm and confirm the minimum in the axial temperature variation first observed by Pennes [4] for an arm in a warm environment.     

Visualization of body thermoregulation by infrared imaging

This paper analyzes two mechanisms applied on the human body in order to study the thermoregulatory system according to heat generation and heat loss. Two approaches are presented. The first approach is based on plethysmography, where an armband is placed on the forearm in order to modulate the blood flow. The second approach uses a cold stimulation. The visualization is achieved using infrared imaging devices. The resulting images reveal a temperature balance between the stimulated and the non-stimulated hands. The thermal behavior and typical thermographic recordings on each subject are discussed and analyzed in response to different stimulations.     

利用δ声波场和近红外光实现乳腺肿瘤的精确定位

本文提出一种利用δ声波场和近红外光漫射理论实现球形乳腺肿瘤精确定位的新颖思路。通过构建一个δ波形的声波场,作用到乳腺组织中从而改变组织内某一点的光学特性参数,这种改变对组织表面光分布的影响可以视为微扰,通过控制其作用点在深度上扫描,测得一系列乳腺组织表面一级微扰光分布,从中提取肿瘤与正常组织的差异特性,实现乳腺肿瘤的精确定位和大小测定。该方法具有广泛的临床医学应用前景,为乳腺癌的早期检测提供一种全新思路。