• Journal of Korean Medicine Rehabilitation
• /
• v.24 no.3
• /
• pp.111-119
• /
• 2014

Objectives This study aims to analyze a thermal distribution in biological living tissue during warm needling therapy by using a finite element method. The analysis provides an understanding of warm needling's efficacy and safety. Methods A model which consisted of four-layered tissue and stainless steel needle was adopted to analyze the thermal distribution in living tissue with a bioheat transfer analysis. The governing equation for the analysis was a Pennes' bioheat equation. A heat source characteristic of warm needling therapy was obtained by previous experimental measurements. The first analysis of the time-dependent temperature distribution was conducted through points on a boundary between the needle and the tissue. The second analysis was conducted to visualize the horizontal temperature distribution. Results When heat source's peak temperatures was above $500^{\circ}C$ and temperature rising rates were relatively slow, the peak temperature at skin surface exceeded a threshold of pain and tissue damage ($45^{\circ}C$), whereas when the peak temperature was around $400^{\circ}C$, the peak temperature at the skin surface was within a safe limit. In addition, the conduction of combustion energy from the moxa was limited to the skin layer around the needle. Conclusions The results suggest that the skin layer around the needle can be heated effectively by warm needling therapy, but it appears to have little effect at the deeper tissue. These findings enhance our understanding of the efficacy and the safety of the warm needling therapy.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2007.11a
• /
• pp.509-510
• /
• 2007

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.8
• /
• pp.757-762
• /
• 2011

We investigate the transient temperature response in biological tissue whose surface is exposed to alternately varying sinusoidal oscillation. Based on the Pennes bio-heat equation, we apply numerical analysis using a finite element method to find the effects of the physical properties of the skin layers. Three layers of tissue-epidermis, dermis, and subcutaneous-are considered as the solution region. We investigate the effects of different properties of the skin layers on the temperature profile. We also investigate the effects of the perfusion rate for the dermis, which is the most sensitive layer. The results show that the temperature profile of tissue depth has a discontinuous point when different physical properties are used.

• Journal of applied mathematics & informatics
• /
• v.9 no.2
• /
• pp.731-744
• /
• 2002

The temperature distribution in human skin and subdermal tissue layer is presented using bioheat transfer equation. The body temperature is determined by the balance between heat produced and heat lost by our body. The time-dependent solutions have been found to be affected by the metabolic heat generation rate, blood mass flow, the rate of evaporation of perspiration and also by the atmospheric temperature. The analytic solutions for different layers have been calculated numerically and are also shown graphically.

• Journal of the Korean Physical Society
• /
• v.73 no.9
• /
• pp.1289-1294
• /
• 2018

The present study aims to investigate experimentally and theoretically thermal ablation in soft tissues by using high intensity focused ultrasound (HIFU) to assess tissue damage during HIFU thermotherapy. The HIFU field was calculated by solving the axisymmetric Khokhlov-Zabolotskaya-Kuznetsov equation from the frequency-domain perspective. The temperature field was calculated by solving Pennes' bioheat transfer equation, and the thermal dose required to create a thermal lesion was calculated by using the thermal dose formula based on the thermal dose of a 240-min exposure at $43^{\circ}C$. In order to validate the simulation results, we performed thermal ablation experiments in a tissue-mimicking phantom and ex-vivo porcine liver for two different HIFU source conditions by using a 1.1-MHz, single-element, spherically focused HIFU transducer. The small difference between the measured and the predicted lesion sizes suggests that the implementation of the numerical model used here should be modified to iteratively allow for temperature-dependent changes in the physical properties of tissues.

• Journal of Biomedical Engineering Research
• /
• v.23 no.6
• /
• pp.467-475
• /
• 2002

Hyperthermia using transrectal thermal probes has been used for a noninvasive treatment of prostate diseases. However it is known that heating the rectal wall at excessively high temperature can lead to destruction of the rectal mucous membrane. and it is difficult to maintain an optimum temperature over the entire prostate. Thus, a more accurate understanding of the heat transfer mechanism between prostate and hyperthermia system is needed Numerical analysis was performed to investigate how the cold/warm stimulations on the prostate surface affect the temperature distribution in the prostate model. The general purpose software "FLUENT" was used for obtaining a finite volume solution to the unsteady conduction equation and to calculate the time-varying temperature in the prostate. Effects of the warm/cold stimulations and the stimulation frequency on the temperature distribution were simulated. and we visualized how hyperthermia affected the inside of the prostate. It was found that the effect of hyperthermia by using a typical heating method is limited due to the low thermal conductivity of the prostate. Consecutive repetitions of warm and cold stimulations were considered to provide the thermal irritations inside a prostate. The effects of temperature difference and duration of warm/cold stimulations were investigated, and basic data for the optimum period and effective patterns of stimulations were obtained. A simplified bioheat equation was also solved to describe effects of the blood flow on the blood-tissue heat transfer. The effect of blood flow was not dominant compared to that of warm/cold stimulations. These results might be used as data for design of prostate treating probe, prostatic therapy and thermal stimulation effects on the prostate.

• Journal of Navigation and Port Research
• /
• v.39 no.4
• /
• pp.285-291
• /
• 2015

Exposure to cold sea water can be life-threatening to the drowned individual. Although appropriate life jacket can be usually be provided for the buoyance at the drowning accident, heat loss can make the drowned individual experience the hypothermia. Inflatable life jackets filled with inflatable air pocket can increase the thermal protection as well as the buoyancy force. Because it is important to know how the human body behaves unde the different life jacket, present study compares the thermal insulation capacity of solid type life jacket with that of inflatable life jacket. In order to represent the insulation capacity of life jacket, thermal resistance is estimated based on the assumption of steady-state. Also, a transient three-dimensional thermal distribution of the thigh is analyzed by using finite element method implementing the Pennes bioheat equation. The finite element model is a segmental, multi-layered representation of the body section which considers the heat conduction within tissue, bone, fat and local blood flow rate.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.22 no.5
• /
• pp.644-653
• /
• 2019

With the development of themal observatory device(TOD), thermal camouflage system has been applied not only to the weapon system but also to the combat suit for soldiers. In this paper, the characteristic of thermal radiation of human body depending on the clothing material properties was analyzed through numerical simulations. The bioheat equation with thermoregulatory model was solved to obtain the realistic surface temperature of human body and these results are combined with the emissivity of human skin and clothing in order to calculate the thermal signature from the human body. According to each thermal resistance of clothing, the optimal background radiance which makes contrast radiance intensity(CRI) be lowest is different. Also, the average CRI variation per thermal resistance change is about twice as much as the case of evaporative resistance change.

• The Journal of Korean Physical Therapy
• /
• v.22 no.6
• /
• pp.77-83
• /
• 2010

Purpose: The aim of this study was to do numerical analysis of the wavelength dependence in low level laser therapy (LLLT) using a finite element method (FEM). Methods: Numerical analysis of heat transfer based on a Pennes' bioheat equation was performed to assess the wavelength dependence of effects of LLLT in a single layer and in multilayered tissue that consists of skin, fat and muscle. The three different wavelengths selected, 660 nm, 830 nm and 980 nm, were ones that are frequently used in clinic settings for the therapy of musculoskeletal disorders. Laser parameters were set to the power density of 35.7 W/$cm^2$, a spot diameter of 0.06 cm, and a laser exposure time of 50 seconds for all wavelengths. Results: Temperature changes in tissue based on a heat transfer equation using a finite element method were simulated and were dominantly dependent upon the absorption coefficient of each tissue layer. In the analysis of a single tissue layer, heat generation by fixed laser exposure at each wavelength had a similar pattern for increasing temperature in both skin and fat (980 nm > 660 nm > 830 nm), but in the muscle layer 660nm generated the most heat (660 nm ${\gg}$ 980 nm > 830 nm). The heat generation in multilayered tissue versus penetration depth was shown that the temperature of 660 nm wavelength was higher than those of 830 nm and 980 nm Conclusion: Numerical analysis of heat transfer versus penetration depth using a finite element method showed that the greatest amount of heat generation is seen in multilayered tissue at = 660 nm. Numerical analysis of heat transfer may help lend insight into thermal events occurring inside tissue layers during low level laser therapy.

• Journal of Biomedical Engineering Research
• /
• v.44 no.1
• /
• pp.85-91
• /
• 2023

Stress urinary incontinence (SUI) occurs when abdominal pressure increases, such as sneezing, exercising, and laughing. Surgical and non-surgical treatments are the common methods of SUI treatment; however, the conventional treatments still require continuous and invasive treatment. Laser have been used to treat SUI, but excessive temperature increase often causes thermal burn on urethra tissue. Therefore, the optimal conditions must be considered to minimize the thermal damage for the laser treatment. The current study investigated the feasibility of the laser irradiation condition for SUI treatment using non-ablative 980 nm laser from a safety perspective through numerical simulations. COMSOL Multiphysics was used to analyze the numerical simulation model. The Pennes bioheat equation with the Beer's law was used to confirm spatio-temporal temperature distributions, and Arrhenius equation defined the thermal damage caused by the laser-induced heat. Ex vivo porcine urethral tissue was tested to validate the extent of both temperature distribution and thermal damage. The temperature distribution was symmetrical and uniformly observed in the urethra tissue. A muscle layer had a higher temperature (28.3 ℃) than mucosal (23.4 ℃) and submucosal layers (25.5 ℃). MT staining revealed no heat-induced collagen and muscle damage. Both control and treated groups showed the equivalent thickness and area of the urethral mucosal layer. Therefore, the proposed numerical simulation can predict the appropriate irradiation condition (20 W for 15 s) for the SUI treatment with minimal temperature-induced tissue.