TIME DEPENDENT TEMPERATURE DISTRIBUTION MODEL IN LAYERED HUMAN DERMAL PART
|
|
- Nora Moore
- 6 years ago
- Views:
Transcription
1 VOL. 8, No. II, DECEMBER, 0, pp TIME DEPENDENT TEMPERATURE DISTRIBUTION MODEL IN LAYERED HUMAN DERMAL PART Saraswati Acharya*, D. B. Gurung, V. P. Saxena Department of Natura Sciences (Mathematics), Schoo of Science, Kathmandu University Sagar Institute of Research and Technoogy, Bhopa, India *Corresponding address: Received 7 January, 0; Revised August, 0 ABSTRACT The paper deveoped appication of finite eement method with inear function in the study of temperature distribution in the ayers of derma partstratum corneum, stratum germinativum, papiary region, reticuar region and subcutaneous tissues as eements. The method is appied to obtain the numerica soution of governing differentia equation for one dimensiona unsteady state bio-heat transfer using suitabe vaues of parameters that effect the heat transfer in human body. The numerica resuts obtained are exhibited graphicay for various atmospheric temperatures for comparative study of temperature distribution profies. The oss of heat from the outer surface of the body to the environment is taken due to convection, radiation and sweat evaporation. Keywords: finite eement method, human derma part, bio-heat equation INTRODUCTION Human skin is a very compex tissue consisting of severa distinct ayers and components. Hence it exhibits compex materia behavior. It is the argest organ of the human body and it has severa functions. The most obvious function is protection of the body against externa infuences and it heps to reguate the body temperature. Skin temperature distribution of the human body is a compex interaction of physica heat exchange processes and the potentia for physioogica adjustment. Temperature infuences the functioning of bioogica systems. For humans, the core temperature varies within narrow bounds around 37 C. The human body acts ike an unsteady state system with a continuous input and not a reguar output. This continuous input provokes a mass storage in the body, eading for a reaction in order to burn it. This process is caed metaboism rate. The norma range of human body temperature varies due to metaboic rate of an individua. Higher metaboic rate causes higher norma body temperature and ower metaboic rate causes ower the norma body temperature. Other factors that might affect the body temperature of an individua may be the time of day or the part of body where the temperature is measured. The body temperature is ower in the morning due to the resting condition and higher at night after a day of muscuar activity and after food intake as we as the thermoreguatory processes of the body are carried on by different methods: convection, conduction, radiation and evaporation. Most of the heat of the body is ost due to convection because it is in contact with fuid such as air and water. Convective transfer by bood pays a key roe. Densities of bood vesses are very high and hence it is hard to compute a detaied temperature distribution in even a sma part of the body whie accounting for a of the bood vesses individuay. Fortunatey, the effects of the bood vesses can be described coectivey with some success by the hep of Pennes bio heat equation. So externa (cimate) and interna (metaboic) heat sources infuences body temperature. Heavy exercise, iness and not hot and humid but aso cod and windy environments ater body 66
2 VOL. 8, No. II, DECEMBER, 0, pp temperature outside the norma range. Ambient temperature, humidity, air movement and radiant heat from the sun as we as warm and cod surfaces contribute to cimate heat stress. Body temperature refects the carefu baance between heat production and heat oss. There is a continuous heat exchange between the body and the environment. Bio- directiona routes for heat exchange are: convection (Cv), conduction (Cd), and radiation (R). There are aso two unidirectiona routes: Metaboic heat (M) increases the therma oad, evaporation (E) decreases the oad. The net heat storage (S) formua is: S= M+/- Cv+/-Cd+/-R-E When the net heat storages (S) is positive, body temperature wi rise and when S is negative, it wi go down. The temperature mode on human body deas the study of temperature distributions on the ayers of human skin exposed to sources of temperature, e.g. surrounding temperature or any sources. So any irreguarities in the temperature distribution in derma ayers due to abnorma environment cause the disturbances in thermoreguation. Hence the study of temperature distribution has the cinica and theoretica importance. Temperature pays an important roe in the functioning of bioogica system. To predict tissue temperatures as we as the temperature of the skin ayers infuence of bood fow must be accounted. The temperature dependence of bioogica processes can be used to cinica effect. For exampe, hyperthermia treatment against cancer, cooing of the head to prevent hair oss as a side effect of chemotherapy and cooing of patients during major surgery to protect the brain are such exampes. Mathematica mode for temperature distributions have a roe both in treatment and diagnosis. They can aid in predicting the time in the course of treatment or giving information on the temperature where thermometry is acking. So numerica modeing and anaysis of the skin has numerous appications in medicine but aso in aesthetics, in computer, graphic, etc. MATHEMATICAL MODEL AND MODEL ASSUMPTIONS The mathematica mode used for bioheat transfer is based on the Pennes equation [8] which incorporates the effect of metaboism and bood perfusion in to the standard therma diffusion equation. This equation is written in simpified form as ρc T t = K. T + M(T a T) + S () where, ρ = Tissue density c = Tissue specific heat K = Tissue therma conductivity m b = Bood mass fow rate c b = Bood specific heat ω b = The bood perfusion rate T a = Arteria bood temperature M = ω b m b c b 67
3 VOL. 8, No. II, DECEMBER, 0, pp S = Metaboic heat generation rate The oss of heat from the skin surface due to convection, radiation and evaporation is considered. So the mixed boundary condition is K T x Skin surface = h (T T ) + LE () where, h = Combined heat transfer coefficient due to convection and radiation T L E = Surrounding temperature = Latent heat of evaporation = Rate of sweat evaporation The inner body core temperature T b is assumed to be 37 0 C. The thickness of stratum corneum, stratum germinativum, papiary region, reticuar region and subcutaneous tissue have been considered as,, 3, 3, 5 respectivey and T 0,, T, and T 5 = T b are the noda temperatures at a distances x = 0, x =, x =, x = 3, x = and x = 5. T (i), i =,, 3,, 5 be the temperature functions in the ayers stratum corneum, stratum germinativum, papiary region, reticuar region and subcutaneous tissue respectivey. T o 3 T T () T () Stratum Corneum Stratum Germinativum T (3) Papiary Region 5 T () Reticuar Region T 5 T (5) Subcutaneous Tissue Body core temperature = T b Figure : Schematic diagram of five ayers of derma part with noda points The anatomica structure of human derma part makes it reasonabe to consider M and S zero in stratum corneum. In the mode, the therma conductivity in the ayers of derma part is 68
4 VOL. 8, No. II, DECEMBER, 0, pp considered as constant. A the assumptions for parameters in the ayers of derma part can be summed up as: Quantity T T a K M S Outer Boundary Stratum Corneum (0 x ) Stratum Germinativum ( x ) Papiary Region ( x 3 ) Reticuar Region ( 3 x ) Subcutaneous Tissue ( x 5 ) T 0 T () = T 0 + T0 x T a () T = T + T x = 3T 3 + T 3 x = x = 0 K () () T a = 0 (3) T a = T b T a () = T b T 5 = 5 T 5 + T 5 (5) x T a 5 5 = T b M () = 0 S () = 0 K () M () = 0 S () K (3) M 3 = x K () M = x m m = x S 3 = x S = x K (5) M 5 = m S 5 = s s s s SOLUTION OF THE PROBLEM The variationa integra form of () in one dimensiona unsteady state case together with outer skin boundary condition () is given by I = 0 K dt dx + M(Ta T) ST + ρc T t dx + h T T + LET We write I separatey for the Five ayers; I for stratum corneum, I for stratum germinativum, I 3 for papiary region, I for reticuar region and I 5 for subcutaneous tissue, so () () () dt () () () () T I K M T T S T c dx dx A t 0 ht 0 Ta LET0 I () () dt () () () () T K M T A T S T c dx dx t () I 3 (3) (3) dt (3) (3) (3) (3) T K M T A T S T c dx dx t 3 (3) 69
5 VOL. 8, No. II, DECEMBER, 0, pp I () () dt T K M TA T S T c dx dx t () () () () () 3 I 5 (5) (5) dt T K M TA T S T c dx dx t (5) (5) (5) (5) (5) 5 Evauating the integra I i, i =,,,5 with the hep of ayers wise assumptions, we get the foowing system of equations given beow I A BT 0 DT 0 ET FT 0T T0 T T0T t I A BT CT DT ET FTT T T TT t I3 A3 B3T C3T3 D3T E3T3 F3T T3 T T3 TT 3 t I A BT3 CT DT3 ET FT T 3 T T3T t I5 A5 B5T C5T5 D5T E5T5 F5T5 T T5 TT 5 t where A i, B i, D i, E i, F i, i 5 and C j, j 5 are a constants depending upon the vaue of physica and physioogica parameters of derma part, and A = ht ; B = LE ht ; D = K () A = 0; B = P ; C = P ; D = R ; ; E = R A 3 = T b m( 3 ) ( ) ; + h ; E = K() ; F = R ; F = K() B 3 = P ( 3 ) N ( ); C 3 = P ( 3 3 ) + N ( ); D 3 = R + Q ( ); E 3 = R + 8 Q ( 3 ) 3 ; F 3 = R + Q ( 3 ) 3 ; A = T b m ( ) ( ); B = P ( ) N ( ); C = P ( ) + N ( ); 70
6 VOL. 8, No. II, DECEMBER, 0, pp D = R 3 + Q ( ); E = R 3 + Q ( ); F = R 3 Q ( ); A 5 = mt b ( 5 ); B 5 = P 3 ( 5 ); C 5 = P 3 ( 5 ); D 5 = R + Q 3 ; E 5 = R + Q 3 ; F 5 = R + Q 3. Where, P = s( ) ; 6( ) R = Q = K () ; ( ) P = m ; ( 3 )( ) Q 3 = m( 5 ); 6 R = K(3) ; ( 3 ) N = T b m ; P 6( ) 3 = (T b m+s); Q = R 3 = K() ; ( 3 ) m ; ( 3 )( ) s. 6( ) R = K(5) ; ( 5 ) α = ρc ; β = ρc ( ) ; γ = ρc ( 3 ) ; μ = ρc ( 3 ) ; η = ρc ( 5 ) We differentiate I with regard to the noda temperatures T o,, T,, and T 5 and set di dt j = 0, j = 0,,, 3,. Since T 5 = T b (the body core temperature), we get the system of equations in matrix form: CT + PT = W (3) Where, C = 7
7 VOL. 8, No. II, DECEMBER, 0, pp T = T o T D F O O O F (E +D ) F O O, P = O F (E + D 3 ) F 3 O O O F 3 (E 3 + D ) F O O O F (E +D 5 ), T = T 0 T t T t t t NUMERICAL RESULTS AND DISCUSSION Parameter vaues used in the mode are as considered in Tabe beow. Tabe : Parameter vaues considered [3, ]. Parameter Vaue Unit K () w/m 0 C K () w/m 0 C K (3) w/m 0 C K () w/m 0 C K (5) w/m 0 C L.0 6 J/kg h 6.80 w/m 0 C E.90 - w/m ρ 050 kg/m C J/kg M w 0 C S 56.0 w/m 3 Tabe : Two sets of derma ayers Sets (m) (m) 3 (m) (m) 5 (m) I II And, for noda temperatures at t = 0, we are taking the equation T(x, t) = T(x, 0) + px considering initia temperature.87 0 C at skin surface because at norma atmospheric temperatures skin surface temperature is around.87 0 C. To sove the system of ordinary 7
8 Noda temperatures ( o C) Noda temperatures ( o C) KATHMANDU UNIVERSITY JOURNAL OF SCIENCE, ENGINEERING AND TECHNOLOGY VOL. 8, No. II, DECEMBER, 0, pp differentia equation (3), we use Crank-Nicoson method. According to the method, the system of equation (3) can be written as C + t P T(i+) = C t P T(i) + tw () where t is the time interva. Using the above numerica vaues in equation () the profies for temperature distribution in the ayers of derma part so obtained is as shown in figures beow taking various atmospheric temperatures T T T 0 0 T Fig.: Noda Temperature for Set-I at T a = 5 0 C Fig. : Noda Temperature for Set-II at T a = 5 0 C The temperature distribution profies for noda temperatures T i (i=0,,, 3, ) for the two set of skin thicknesses at different atmospheric temperatures with time dependence are shown in figures-6. From these figures it has been observed that the curve for T i rise faster in set I and reach steady state case earier than set II. It is aso observed that steady state temperatures at these depths are higher for set I than to set II. These are due to ower thickness of set I than to set II. These figures revea that the temperatures Ti at high atmospheric temperatures are higher than ower atmospheric temperatures. This is due to the more effect on body temperature at higher atmospheric temperatures. So the thicknesses of derma ayers and different atmospheric temperatures have the significant effects in temperature profies. 73
9 Noda temperatures ( o C) Noda temperatures ( o C) Noda temperatures ( o C) Noda temperatures ( o C) KATHMANDU UNIVERSITY JOURNAL OF SCIENCE, ENGINEERING AND TECHNOLOGY VOL. 8, No. II, DECEMBER, 0, pp T 0 8 T T 0 T Fig. 3: Noda Temperature for Set-I at T a = 37 0 C Fig. : Noda Temperature for Set-II at T a = 37 0 C T T 0 T 0 5 T Fig. 5: Noda Temperature for Set-I at T a = 5 0 C Fig. 6: Noda Temperature for Set-II at T a = 5 0 C 7
10 VOL. 8, No. II, DECEMBER, 0, pp CONCLUSION The resuts are comparabe with those obtained in [] who considered ony three ayers of skin epidermis, dermis and subcutaneous tissue. The differences in the resuts whatsoever may be, due to the extension up to five ayers of derma part Stratum Corneum, Stratum Germinativum, Papiary Region, Reticuar Region and Subcutaneous Tissue. Our mode gives better profies for temperature distribution in derma ayers. This is because the mode has incorporated more feasibe ayers and has taken significant biophysica parameters. The epidermis ayer contains five ayers and dermis contains two ayers. So bioogicay derma region has eight ayers. We may further extended the mode considering the eight ayers as eight eements in one dimensiona case. It may be extended for two dimensions and three dimensions. The study may be carried out for predication of therma response for humans exposed to hot or cod environments. The mode may be extended by considering the fied variabe as quadratic or of higher order poynomia in FEM technique to get better temperature profies. The different temperature profie can be viewed by considering this mode at cinica situation, indoor cimate and risk under stressfu condition. It is hepfu for the medica scientist to suggest and evauate for new treatment strategies. It is aso hepfu to maintain indoor temperature. The presented mode can be extended to expore severa probems mentioned above. REFERENCES [] Cooper T E & Trezek G J, A probe technique for determining the therma conductivity of tissue, J. Heat transfer, ASME, 9(97)33. [] Guyton C & Ha E, Text book of medica physioogy, Esevier, India, 009. [3] Gurung D B, Mathematica study of abnorma thermo-reguation in human derma parts, Ph. D. Thesis, Kathmandu University, 007. [] Gurung D B, Saxena V P & Adhikary P R, FEM approach to one dimensiona unsteady state temperature distribution in the derma parts with quadratic shape function, J. App. Math and Informatics, 7(-)(009). [5] Gurung D B, Thermoreguation through skin at ow atmospheric temperatures, KUSET, 5() (009). [6] Gurung D B & Saxena V P, Periphera temperature distribution in human subjects exposed to wind fow, Internationa ejournas of Mathematics and Engineering 73(00) 78. [7] Myers G E, Anaytic methods in conduction heat transfer, MC. Graw Hi Book Co., New York, 97. [8] Pennes H H, Anaysis of tissue and arteria bood temperature in the resting human forearm, J. App. Physio., (98)93. 75
11 VOL. 8, No. II, DECEMBER, 0, pp [9] Crank J, The mathematics of diffusion, Oxford Science Pubications, second edition, 975. [0] Rao S S, The finite eement method in engineering, Esevier, a division of Reed Esevier India Private Limited, (009). [] Per W, Heat and matter distribution in body tissues and the determination of tissue bood fow by oca cearance method, J. Theor. Bio, (96)0. [] Saxena V P, Adhikary P R & Gurung D B, Mathematica estimation of unsteady state burn damage due to hot temperature, The Nepai Math Sc. Report, 7(-)(007). [3] Rao Y V C, Theory and probems of thermodynamics, New Age Internationa (P) Limited, second edition, (000)3. 76
4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationNumerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet
Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport
More informationEffect of Oxygen Injection into Argon Induction Plasmas on Chemically Non-Equilibrium Conditions
Proceedings of 17th Internationa Symposium on Pasma Chemistry, Toronto, Canada, August 7-12, 25 Effect of Oxygen Injection into Argon Induction Pasmas on Chemicay Non-Equiibrium Conditions Nobuhiko Atsuchi
More informationWilson Ogoh and Dominic Groulx Department of Mechanical Engineering Dalhousie University, Halifax, NS, Canada
Presented at the COMSOL Conference 2010 Boston Wison Ogoh and Dominic Groux Department of Mechanica Engineering Dahousie University, Haifax, NS, Canada Sensibe Heat Storage: A heat storage system that
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More informationNumerical simulation of javelin best throwing angle based on biomechanical model
ISSN : 0974-7435 Voume 8 Issue 8 Numerica simuation of javein best throwing ange based on biomechanica mode Xia Zeng*, Xiongwei Zuo Department of Physica Education, Changsha Medica University, Changsha
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.
More informationMass Transport 2: Fluids Outline
ass Transport : Fuids Outine Diffusivity in soids, iquids, gases Fick s 1st aw in fuid systems Diffusion through a stagnant gas fim Fick s nd aw Diffusion in porous media Knudsen diffusion ass Transfer
More informationMA 201: Partial Differential Equations Lecture - 10
MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary
More informationEffect of transport ratio on source term in determination of surface emission coefficient
Internationa Journa of heoretica & Appied Sciences, (): 74-78(9) ISSN : 975-78 Effect of transport ratio on source term in determination of surface emission coefficient Sanjeev Kumar and Apna Mishra epartment
More informationMONTE CARLO SIMULATIONS
MONTE CARLO SIMULATIONS Current physics research 1) Theoretica 2) Experimenta 3) Computationa Monte Caro (MC) Method (1953) used to study 1) Discrete spin systems 2) Fuids 3) Poymers, membranes, soft matter
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.
Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationDYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE
3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses
More information14 Separation of Variables Method
14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt
More informationImproving the Reliability of a Series-Parallel System Using Modified Weibull Distribution
Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationHYDROGEN ATOM SELECTION RULES TRANSITION RATES
DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationSupplementary Material: An energy-speed-accuracy relation in complex networks for biological discrimination
Suppementary Materia: An energy-speed-accuracy reation in compex networks for bioogica discrimination Feix Wong, 1,2 Arie Amir, 1 and Jeremy Gunawardena 2, 1 Schoo of Engineering and Appied Sciences, Harvard
More informationApproximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation
Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische
More informationCFD MODELLING OF DIRECT CONTACT STEAM INJECTION
Fifth Internationa Conference on CFD in the Process Industries CSIRO, Meourne, Austraia 13-15 Decemer 006 CFD MODELLING OF DIRECT CONTACT STEAM INJECTION Curtis MARSH 1 and Denis WITHERS 1 CFD Design &
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More informationRelated Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage
Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is
More informationInternational Conference on Mechanical, Industrial and Energy Engineering December, 2014, Khulna, BANGLADESH
International Conference on Mechanical, Industrial and Energy Engineering 214 25-2 December, 214, Khulna, BANGLADESH ICMIEE-PI-13581 Analysis of Bio-heat Transfer Problem Using Finite Element Approach
More informationComponentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems
Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,
More informationEstablishment of Weak Conditions for Darboux- Goursat-Beudon Theorem
Georgia Southern University Digita Commons@Georgia Southern Mathematica Sciences Facuty Pubications Department of Mathematica Sciences 2009 Estabishment of Weak Conditions for Darboux- Goursat-Beudon Theorem
More informationCrystallisation of a supercooled spherical nodule in a flow
EUROTHERM 69 Heat and Mass Transfer in Soid-Liquid Phase hange Processes June 25-27, 2003, Bistra caste, Ljubjana, Sovenia Eds.: B. Sarer, D. Gobin rystaisation of a supercooed spherica nodue in a fow
More informationHomotopy Perturbation Method for Solving Partial Differential Equations of Fractional Order
Int Journa of Math Anaysis, Vo 6, 2012, no 49, 2431-2448 Homotopy Perturbation Method for Soving Partia Differentia Equations of Fractiona Order A A Hemeda Department of Mathematics, Facuty of Science
More informationarxiv:nlin/ v2 [nlin.cd] 30 Jan 2006
expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationA Two-Dimensional Mathematical Model to Analyze Thermal Variations in Skin and Subcutaneous Tissue Region of Human Limb during Surgical Wound Healing
Applied Mathematics, 016, 7, 145-158 Published Online February 016 in SciRes. http://www.scirp.org/journal/am http://dx.doi.org/10.436/am.016.7014 A Two-Dimensional Mathematical Model to Analyze Thermal
More informationIn-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017
In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More informationAn approximate method for solving the inverse scattering problem with fixed-energy data
J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationUI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE
UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,
More informationTwo Kinds of Parabolic Equation algorithms in the Computational Electromagnetics
Avaiabe onine at www.sciencedirect.com Procedia Engineering 9 (0) 45 49 0 Internationa Workshop on Information and Eectronics Engineering (IWIEE) Two Kinds of Paraboic Equation agorithms in the Computationa
More informationHigh Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings
Purdue University Purdue e-pubs Internationa Compressor Engineering Conference Schoo of Mechanica Engineering 2010 High Efficiency Deveopment of a Reciprocating Compressor by Carification of Loss Generation
More informationTorsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg
Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................
More informationVTU-NPTEL-NMEICT Project
MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid
More informationOn Integrals Involving Universal Associated Legendre Polynomials and Powers of the Factor (1 x 2 ) and Their Byproducts
Commun. Theor. Phys. 66 (216) 369 373 Vo. 66, No. 4, October 1, 216 On Integras Invoving Universa Associated Legendre Poynomias and Powers of the Factor (1 x 2 ) and Their Byproducts Dong-Sheng Sun ( 孙东升
More informationElements of Kinetic Theory
Eements of Kinetic Theory Statistica mechanics Genera description computation of macroscopic quantities Equiibrium: Detaied Baance/ Equipartition Fuctuations Diffusion Mean free path Brownian motion Diffusion
More informationMassive complex scalar field in the Kerr Sen geometry: Exact solution of wave equation and Hawking radiation
JOURNAL OF MATHEMATICAL PHYSICS VOLUME 44, NUMBER 3 MARCH 2003 Massive compex scaar fied in the Kerr Sen geometry: Exact soution of wave equation and Hawking radiation S. Q. Wu a) Interdiscipinary Center
More informationNonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationExperimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure
Experimenta Investigation and Numerica Anaysis of New Muti-Ribbed Sab Structure Jie TIAN Xi an University of Technoogy, China Wei HUANG Xi an University of Architecture & Technoogy, China Junong LU Xi
More informationHEAT EXCHANGER NETWORK SYNTHESIS CONSIDERING CHANGING PHASE STREAMS
HEAT EXCHANGER NETWORK SYNTHESIS CONSIDERING CHANGING PHASE STREAMS F. S. Liporace a, F. L. P. Pessoa b, and E. M. Queiroz b, a PETROBRAS/CENPES/EB/SAP Cidade Universitária - Iha do Fundão 21949-900, Rio
More informationDISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE
DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More information12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes
Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.
More informationResearch Article Building Infinitely Many Solutions for Some Model of Sublinear Multipoint Boundary Value Problems
Abstract and Appied Anaysis Voume 2015, Artice ID 732761, 4 pages http://dx.doi.org/10.1155/2015/732761 Research Artice Buiding Infinitey Many Soutions for Some Mode of Subinear Mutipoint Boundary Vaue
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 7
Strauss PDEs 2e: Section 4.3 - Exercise 1 Page 1 of 7 Exercise 1 Find the eigenvaues graphicay for the boundary conditions X(0) = 0, X () + ax() = 0. Assume that a 0. Soution The aim here is to determine
More informationarxiv: v1 [math.ca] 6 Mar 2017
Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear
More informationD. Prémel, J.M. Decitre and G. Pichenot. CEA, LIST, F Gif-sur-Yvette, France
SIMULATION OF EDDY CURRENT INSPECTION INCLUDING MAGNETIC FIELD SENSOR SUCH AS A GIANT MAGNETO-RESISTANCE OVER PLANAR STRATIFIED MEDIA COMPONENTS WITH EMBEDDED FLAWS D. Préme, J.M. Decitre and G. Pichenot
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationOn Some Basic Properties of Geometric Real Sequences
On Some Basic Properties of eometric Rea Sequences Khirod Boruah Research Schoar, Department of Mathematics, Rajiv andhi University Rono His, Doimukh-791112, Arunacha Pradesh, India Abstract Objective
More informationNumerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process
ISIJ Internationa Vo 54 (2014) No 2 pp 254 258 Numerica Simuation for Optimizing Temperature Gradients during Singe Crysta Casting Process Aeksandr Aeksandrovich INOZEMTSEV 1) Aeksandra Sergeevna DUBROVSKAYA
More informationMathematical Scheme Comparing of. the Three-Level Economical Systems
Appied Mathematica Sciences, Vo. 11, 2017, no. 15, 703-709 IKAI td, www.m-hikari.com https://doi.org/10.12988/ams.2017.7252 Mathematica Scheme Comparing of the Three-eve Economica Systems S.M. Brykaov
More information), enthalpy transport (i.e., the heat content that moves with the molecules)
Steady-state conseration statements for a composite of ces and airspace In steady state, conseration of moecues requires that the tota fux into a representatie oume of mesophy is equa to the fux out storage
More informationOn the evaluation of saving-consumption plans
On the evauation of saving-consumption pans Steven Vanduffe Jan Dhaene Marc Goovaerts Juy 13, 2004 Abstract Knowedge of the distribution function of the stochasticay compounded vaue of a series of future
More informationHigher dimensional PDEs and multidimensional eigenvalue problems
Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u
More informationWavelet Galerkin Solution for Boundary Value Problems
Internationa Journa of Engineering Research and Deveopment e-issn: 2278-67X, p-issn: 2278-8X, www.ijerd.com Voume, Issue 5 (May 24), PP.2-3 Waveet Gaerkin Soution for Boundary Vaue Probems D. Pate, M.K.
More informationPolar Snakes: a fast and robust parametric active contour model
Poar Snakes: a fast and robust parametric active contour mode Christophe Coewet To cite this version: Christophe Coewet. Poar Snakes: a fast and robust parametric active contour mode. IEEE Int. Conf. on
More informationAdjustment of automatic control systems of production facilities at coal processing plants using multivariant physico- mathematical models
IO Conference Series: Earth and Environmenta Science AER OEN ACCESS Adjustment of automatic contro systems of production faciities at coa processing pants using mutivariant physico- mathematica modes To
More informationMultiscale Multilevel Approach to Solution of Nanotechnology Problems
Mathematica Modeing and Computationa Physics 2017 Mutiscae Mutieve Approach to Soution of Nanotechnoogy Probems Sergey Poyakov 1,2, and Viktoriia Podryga 1,3, 1 Kedysh Institute of Appied Mathematics of
More informationLecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling
Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest
More informationTheory of Generalized k-difference Operator and Its Application in Number Theory
Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication
More informationСРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS
СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА Милко Стоянов Милошев 1, Константин Савков Казаков 2 Висше Строително Училище Л. Каравелов - София COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS
More informationMA 201: Partial Differential Equations Lecture - 11
MA 201: Partia Differentia Equations Lecture - 11 Heat Equation Heat conduction in a thin rod The IBVP under consideration consists of: The governing equation: u t = αu xx, (1) where α is the therma diffusivity.
More informationMultigrid Method for Elliptic Control Problems
J OHANNES KEPLER UNIVERSITÄT LINZ Netzwerk f ür Forschung, L ehre und Praxis Mutigrid Method for Eiptic Contro Probems MASTERARBEIT zur Erangung des akademischen Grades MASTER OF SCIENCE in der Studienrichtung
More information3.10 Implications of Redundancy
118 IB Structures 2008-9 3.10 Impications of Redundancy An important aspect of redundant structures is that it is possibe to have interna forces within the structure, with no externa oading being appied.
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationAvailable online at ScienceDirect. Energy Procedia 49 (2014 ) SolarPACES R. Yogev a, A.
Avaiabe onine at www.sciencedirect.com ScienceDirect Energy Procedia 49 (2014 ) 1061 1070 SoarPACES 2013 PCM storage system with integrated active heat pipe R. Yogev a, A. Kribus a * a Schoo of Mechanica
More informationA NOVEL WAVELET-GALERKIN METHOD FOR MOD- ELING RADIO WAVE PROPAGATION IN TROPO- SPHERIC DUCTS
Progress In Eectromagnetics Research B, Vo 36, 35 52, 2012 A NOVEL WAVELET-GALERKIN METHOD FOR MOD- ELING RADIO WAVE PROPAGATION IN TROPO- SPHERIC DUCTS A Iqba * and V Jeoti Department of Eectrica and
More informationThermophoretic interaction of heat releasing particles
JOURNAL OF APPLIED PHYSICS VOLUME 9, NUMBER 7 1 APRIL 200 Thermophoretic interaction of heat reeasing partices Yu Doinsky a) and T Eperin b) Department of Mechanica Engineering, The Pearstone Center for
More informationSTOCHASTIC Volterra integral equations (SVIEs) is a
IAENG Internationa Journa of Appied Mathematics, 44:4, IJAM_44_4_2 The Petrov-Gaerkin Method for Numerica oution of tochastic Voterra Integra Equations F. Hosseini hekarabi 1, M. Khodabin 2 and K. Maeknejad
More informationRadiation Fields. Lecture 12
Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence
More informationSEA Subsystem Attribute Correction Based on Stiffness Multipliers
Internationa Conference on Eectrica, Mechanica and Industria Engineering (ICEMIE 016) SEA Subsystem Attribute Correction Based on Stiffness Mutipiers Jie Gao1,* and Yong Yang 1 Coege of hysics and Information
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationA Solution to the 4-bit Parity Problem with a Single Quaternary Neuron
Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria
More informationInternational Journal of Advance Engineering and Research Development
Scientific Journa of Impact Factor (SJIF): 4.4 Internationa Journa of Advance Engineering and Research Deveopment Voume 3, Issue 3, March -206 e-issn (O): 2348-4470 p-issn (P): 2348-6406 Study and comparison
More information6 Wave Equation on an Interval: Separation of Variables
6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.
More informationSimplified analysis of EXAFS data and determination of bond lengths
Indian Journa of Pure & Appied Physics Vo. 49, January 0, pp. 5-9 Simpified anaysis of EXAFS data and determination of bond engths A Mishra, N Parsai & B D Shrivastava * Schoo of Physics, Devi Ahiya University,
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationNumerical methods for PDEs FEM - abstract formulation, the Galerkin method
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Numerica methods for PDEs FEM - abstract formuation, the Gaerkin method Dr. Noemi Friedman Contents of the course Fundamentas of functiona
More informationWave Equation Dirichlet Boundary Conditions
Wave Equation Dirichet Boundary Conditions u tt x, t = c u xx x, t, < x 1 u, t =, u, t = ux, = fx u t x, = gx Look for simpe soutions in the form ux, t = XxT t Substituting into 13 and dividing
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationScienceDirect. Numerical modelling of debris bed water quenching
Avaiabe onine at www.sciencedirect.com ScienceDirect Procedia IUTAM 15 2015 64 71 IUTAM Symposium on Mutiphase fows with phase change: chaenges and opportunities, Hyderabad, India December 08 December
More informationMULTI-PERIOD MODEL FOR PART FAMILY/MACHINE CELL FORMATION. Objectives included in the multi-period formulation
ationa Institute of Technoogy aicut Department of echanica Engineering ULTI-PERIOD ODEL FOR PART FAILY/AHIE ELL FORATIO Given a set of parts, processing requirements, and avaiabe resources The objective
More information2.1. Cantilever The Hooke's law
.1. Cantiever.1.1 The Hooke's aw The cantiever is the most common sensor of the force interaction in atomic force microscopy. The atomic force microscope acquires any information about a surface because
More informationData Search Algorithms based on Quantum Walk
Data Search Agorithms based on Quantum Wak Masataka Fujisaki, Hiromi Miyajima, oritaka Shigei Abstract For searching any item in an unsorted database with items, a cassica computer takes O() steps but
More informationA Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem
Journa o mathematics and computer science 14 (15) 87-96 A Fictitious ime Integration Method or a One-Dimensiona Hyperboic Boundary Vaue Probem Mir Saad Hashemi 1,*, Maryam Sariri 1 1 Department o Mathematics,
More informationChapter 5. Wave equation. 5.1 Physical derivation
Chapter 5 Wave equation In this chapter, we discuss the wave equation u tt a 2 u = f, (5.1) where a > is a constant. We wi discover that soutions of the wave equation behave in a different way comparing
More information