Supplementary Information

Similar documents
Outline. Review Homework Problem. Review Homework Problem II. Review Dimensionless Problem. Review Convection Problem

Comparing Different Estimators for Parameters of Kumaraswamy Distribution

Analysis of Stress in PD Front End Solenoids I. Terechkine

Transistor configurations: There are three main ways to place a FET/BJT in an architecture:

ECSE Partial fraction expansion (m<n) 3 types of poles Simple Real poles Real Equal poles

On a Z-Transformation Approach to a Continuous-Time Markov Process with Nonfixed Transition Rates

On imploding cylindrical and spherical shock waves in a perfect gas

INF 5460 Electronic noise Estimates and countermeasures. Lecture 13 (Mot 10) Amplifier Architectures

The Central Limit Theorems for Sums of Powers of Function of Independent Random Variables

Capítulo. of Particles: Energy and Momentum Methods

Spectrum of The Direct Sum of Operators. 1. Introduction

Available online at J. Math. Comput. Sci. 2 (2012), No. 4, ISSN:

Analytical solution of tank drainage for electrically conducting power law fluid

GENERALIZED FRACTIONAL INTEGRAL OPERATORS AND THEIR MODIFIED VERSIONS

BINOMIAL THEOREM OBJECTIVE PROBLEMS in the expansion of ( 3 +kx ) are equal. Then k =

Existence and Smoothness of Solution of Navier-Stokes Equation on R 3

A Weighted Moving Average Process for Forcasting

Comparison between Fourier and Corrected Fourier Series Methods

COST OPTIMIZATION OF SLAB MILLING OPERATION USING GENETIC ALGORITHMS

Statistical Optics and Free Electron Lasers

Exact Solution of Unsteady Tank Drainage for Ellis Fluid

( ) 1 Comparison Functions. α is strictly increasing since ( r) ( r ) α = for any positive real number c. = 0. It is said to belong to

Consider the time-varying system, (14.1)

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

ELECTROOSMOTIC FLOW IN A MICROCHANNEL PACKED WITH MICROSPHERES

ABSOLUTE INDEXED SUMMABILITY FACTOR OF AN INFINITE SERIES USING QUASI-F-POWER INCREASING SEQUENCES

Numerical Solution of Sine-Gordon Equation by Reduced Differential Transform Method

ON POINTWISE APPROXIMATION OF FUNCTIONS BY SOME MATRIX MEANS OF FOURIER SERIES

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION

PRESSURE AND PRESSURE DERIVATIVE ANALYSIS FOR PSEUDOPLASTIC FLUIDS IN VERTICAL FRACTURED WELLS

THE SOIL STRUCTURE INTERACTION ANALYSIS BASED ON SUBSTRUCTURE METHOD IN TIME DOMAIN

1 Notes on Little s Law (l = λw)

PRICING AMERICAN PUT OPTION WITH DIVIDENDS ON VARIATIONAL INEQUALITY

Effect of Material Gradient on Stresses of Thick FGM Spherical Pressure Vessels with Exponentially-Varying Properties

Latticed pentamode acoustic cloak (supplementary Info)

S, we call the base curve and the director curve. The straight lines

Available online at J. Math. Comput. Sci. 4 (2014), No. 4, ISSN:

Lecture 17: Kinetics of Phase Growth in a Two-component System:

Cameras and World Geometry

Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation

Generalized Fibonacci-Type Sequence and its Properties

CAPACITY ANALYSIS OF ASYMPTOTICALLY LARGE MIMO CHANNELS. Georgy Levin

Calculation of maximum ground movement and deformation caused by mining

The Non-Truncated Bulk Arrival Queue M x /M/1 with Reneging, Balking, State-Dependent and an Additional Server for Longer Queues

Optimal Operating Policies of EPQ Model for Deteriorating Items with Time Dependent Production having Production Quantity Dependent Demand

Pulsatile Flow of Hershel- Bulkley Fluid in Tapered Blood Vessels

Multivector Functions

arxiv: v4 [math.pr] 20 Jul 2016

The Alpha-Logarithmic Series Distribution of Order k and Some of Its Applications

هقارنت طرائق تقذير هعلواث توزيع كاها ري املعلوتني

Fujii, Takao; Hayashi, Fumiaki; Iri Author(s) Oguro, Kazumasa.

SUSTAINABLE HEAT FARMING OF GEOTHERMAL SYSTEMS: A CASE STUDY OF HEAT EXTRACTION AND THERMAL RECOVERY IN A MODEL EGS FRACTURED RESERVOIR

EXACT ANALYSIS OF UNSTEADY CONVECTIVE DIFFUSION IN A HERSCHEL BULKLEY FLUID IN AN ANNULAR PIPE

Neutron Slowing Down Distances and Times in Hydrogenous Materials. Erin Boyd May 10, 2005

Engineering Accreditation. Heat Transfer Basics. Assessment Results II. Assessment Results. Review Definitions. Outline

Orthotropic Materials

One of the common descriptions of curvilinear motion uses path variables, which are measurements made along the tangent t and normal n to the path of

Solutions Manual 4.1. nonlinear. 4.2 The Fourier Series is: and the fundamental frequency is ω 2π

Lecture 15: Three-tank Mixing and Lead Poisoning

CHAPTER 5 : SERIES. 5.2 The Sum of a Series Sum of Power of n Positive Integers Sum of Series of Partial Fraction Difference Method

A Weighted Moving Average Process for Forecasting

F.Y. Diploma : Sem. II [AE/CH/FG/ME/PT/PG] Applied Mathematics

6.2 Improving Our 3-D Graphics Pipeline

L-functions and Class Numbers

MODERN CONTROL SYSTEMS

, on the power of the transmitter P t fed to it, and on the distance R between the antenna and the observation point as. r r t

Class 36. Thin-film interference. Thin Film Interference. Thin Film Interference. Thin-film interference

ON GENERALIZED FRACTIONAL INTEGRAL OPERATORS. ( ρ( x y ) T ρ f(x) := f(y) R x y n dy, R x y n ρ( y )(1 χ )

Indoor Navigation without the use of GPS utilizing Intelligent Data Algorithms

FBD of SDOF Base Excitation. 2.4 Base Excitation. Particular Solution (sine term) SDOF Base Excitation (cont) F=-(-)-(-)= 2ζω ωf

Highly connected coloured subgraphs via the Regularity Lemma

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

F D D D D F. smoothed value of the data including Y t the most recent data.

Inverse Heat Conduction Problem in a Semi-Infinite Circular Plate and its Thermal Deflection by Quasi-Static Approach

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM

Mathematical Models and the Soil Hydraulic Properties

Relations on the Apostol Type (p, q)-frobenius-euler Polynomials and Generalizations of the Srivastava-Pintér Addition Theorems

Stochastic Design of Enhanced Network Management Architecture and Algorithmic Implementations

The Hyperbolic Model with a Small Parameter for. Studying the Process of Impact of a Thermoelastic. Rod against a Heated Rigid Barrier

L8b - Laplacians in a circle

Computer Propagation Analysis Tools

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

SUMMATION OF INFINITE SERIES REVISITED

Research Article On Pointwise Approximation of Conjugate Functions by Some Hump Matrix Means of Conjugate Fourier Series

Boundary Integral Equation Method for Unsteady Two Dimensional Flow in Unconfined Aquifer

Let s express the absorption of radiation by dipoles as a dipole correlation function.

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

ECE 570 Session 7 IC 752-E Computer Aided Engineering for Integrated Circuits. Transient analysis. Discuss time marching methods used in SPICE

A Complex Neural Network Algorithm for Computing the Largest Real Part Eigenvalue and the corresponding Eigenvector of a Real Matrix

Appendix A. Expressions for equation of state parameters of TNT explosion products

Real-time TDDFT simulations within SIESTA. Daniel Sánchez-Portal, Rafi Ullah, Fabiano Corsetti, Miguel Pruneda and Emilio Artacho

UNSTEADY HELICAL FLOWS OF A MAXWELL FLUID

S n. = n. Sum of first n terms of an A. P is

On the Precision of Time-of-Flight Shear Wave Speed Estimation in Homogeneous Soft Solids: Initial Results Using a Matrix Array Transducer

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

Problems and Solutions for Section 3.2 (3.15 through 3.25)

KEY. Math 334 Midterm II Fall 2007 section 004 Instructor: Scott Glasgow

EFFECT OF GRADIENTS IN INHOMOGENEOUS MEDIA ON PLANE WAVE REFLECTION COEFFICIENT RESONANCE

Transcription:

Supplemeay Ifomaio No-ivasive, asie deemiaio of he coe empeaue of a hea-geeaig solid body Dea Ahoy, Daipaya Saka, Aku Jai * Mechaical ad Aeospace Egieeig Depame Uivesiy of Texas a Aligo, Aligo, TX, USA. * Coespodig Auho: email: jaia@ua.edu; 5 W Fis S, Rm, Aligo, TX, USA 769 Ph: + (87) 7-9338; Fax: + (87) 7 95 Theoeical Deivaio of Equaio () This secio deives a elaioship bewee he asie coe empeaue, T coe () a = ad he asie suface empeaue disibuio, T (θ,) of a hea-geeaig ifiie cylide. While simila soluios have bee deived i he pas usig he mehod of Gee s fucios (fo example, Özişik, M.N., Hea coducio, d Ed., Joh Wiley & Sos, 993), his deivaio akes io accou cicumfeeial vaiaio i T, as well as aisoopic hemal coducio wihi he cylide. This deivaio leads o a expessio fo T coe () ha foms he basis fo he echique o o-iusively deemie he coe empeaue. Cosa hea geeaio ae This sub-secio cosides a ifiie cylide of adius R geeaig hea a a uifom ad cosa ae Q, show schemaically i Supplemeay Figue S. The cylide is assumed o

have ohoopic hemal coduciviies, wih values of k ad k θ i he adial ad cicumfeeial diecios especively. The cicumfeeially-vayig, asie empeaue disibuio alog he oue suface a =R, give by T (θ,) is assumed o be kow, fo example, hough a ifaed measueme. The cylide is assumed o be a uifom empeaue iiially. The asie empeaue disibuio i he cylide ca be deemied by solvig he goveig eegy cosevaio equaio subjec o appopiae bouday codiios. Usig he sepaaio of vaiables appoach (Özişik, M.N., Hea Coducio, d Ed., Joh Wiley & Sos, 98), he empeaue field may be spli io wo pas ad assumed o be of he followig fom: Q R T(,, ) T (, ) T (,, ) A J e Bm ( ) J m k cos( m (S.) ) 4k m m m k whee T, compisig he fis wo ems, accous fo ieal hea geeaio, ad T, compisig he hid em accous fo ime- ad θ-depede empeaue a =R. Hee, J efes o he Bessel fucio of he fis kid, ad he adial om N, is give by N, R J( R) (S.) ad he eigevalues λ m ae obaied fom he oos of J m k k. Specifically, λ ae obaied fom he oos of J. Usig he ohogoaliy piciple (Özişik, M.N., Hea Coducio, d Ed., Joh Wiley & Sos, 98), A ae give by A R R J Q 4k N, d (S.3)

The coefficie fucios B m () mus be deemied such ha he soluio saisfies he goveig equaio ad bouday codiios. Noe ha he coe empeaue T coe () may be expessed i ems of he ukow coefficie fucios B m () by puig = i equaio (S.) T coe ( ) T QR, ( ), ( ) coe T coe Ae B ( ) 4k (S.4) Equaio (S.4) shows ha eve hough all coefficie fucios B m () ae eeded o deemie he geeal empeaue disibuio T(,θ,), oly he coefficie fucios B () ae eeded fo measuig he coe empeaue. To do so, he expessio fo T (,θ,) fom equaio (S.) is iegaed wih espec o θ o elimiae all B m () fo m>, followed by use of he mehod of udeemied paamees (Myes, G.E., Aalyical Mehods i Coducio Hea Tasfe, d Ed., AMCHT Publicaios, 998) o deive he followig odiay diffeeial equaio fo B () db d B RJ N R, T, d (S.5) Fuhe, based o he iiial codiio fo he empeaue field, he followig iiial codiio applies fo B () B () (S.6) A soluio fo equaio (S.5) subjec o (S.6) is give by RJ R B ( ) w I ( ) exp d N, (S.7) 3

whee w I (τ) is he cicumfeeially aveaged value of he measued empeaue a he ouside suface, give by w I ( ) T, d (S.8) This complees he deemiaio of he coe empeaue T coe (), which is give by equaio (S.4), whee he coefficie fucios ae give by equaios (S.) ad (S.7). Effec of ime-vayig hea geeaio ae The pevious sub-secio assumed a ime-ivaia ieal hea geeaio ae i he cylide. I his sub-secio, a moe geeal case is cosideed whee he hea geeaio vaies wih ime, give by Q(). Time-vayig hea geeaio is ecoueed i seveal egieeig applicaios, whee he measueme of coe empeaue i espose o Q() is of iees. To develop a echique fo ieal empeaue measueme fo his case, he goveig eegy equaio eeds o be evisied. As show i equaio (S.), he compoes T ad T of he oveall empeaue disibuio accou fo empeaue ise due o ieal hea geeaio ad due o he ime-depede suface empeaue especively. As a esul, vaiaio i Q wih ime affecs oly he T compoe of he empeaue field. A Laplace asfom appoach is used o deive a soluio fo T (,) fo whe Q vaies wih ime. Takig Laplace asfom of he goveig equaio ad bouday codiio fo T esuls i he followig: d T dt s T d d Q( s) k (S.9) 4

5 Whee ) (s Q is he Laplace asfom of he hea geeaio ae Q(). T also saisfies d dt a = (S.) ad T a =R (S.) The soluio o his odiay diffeeial equaio ca be show o be give by R s I s I sk Q T (S.) whee I is he modified Bessel fucio of he fis kid. As a esul, he coe empeaue T coe () fo a ifiie cylide ca be wie as,, ) ( ) ( ) ( ) ( coe coe coe B R s I sk Q L T T T (S.3) whee L - efes o he ivese Laplace asfom. Noe ha B () is give by equaio (S.7). Fo a give Q(), he ivese Laplace asfom may be deemied aalyically whe possible, o umeically usig ivese Laplace asfom umeical algoihms. This povides a echique o pedic he coe empeaue of he cylide i pesece of a ime-vayig hea geeaio wihi he body.

Supplemeay Figues 6

Supplemeay Figue S: Schemaic of he geomey of a hea-geeaig ifiie cylide. The goal of he measueme is o deemie he coe empeaue T(=,) as a fucio of ime fom measued suface empeaue T (θ,) as a fucio of ime. 7

Supplemeay Figue S: Compaiso of asie coe empeaue pediced by he heoeical model wih fiie-eleme simulaio esuls fo wo suface empeaue codiios, (a) cosa T = C, (b) liealy iceasig T ()= C. Boh model ad ( s) fiie-eleme simulaios assume a 3 mm adius cylide wih 57,954 W/m 3 heaig, ad he same hemal popeies as he expeimeal hemal es cell. 8

Supplemeay Figue S3: Measued suface empeaue disibuio aoud he cylide T (θ,τ) a a umbe of imes followig he sa of heaig fo.86 W heaig powe i he hemal es cell i fee covecio coolig codiios. 9

Supplemeay Figue S4: Plo of he maximum esidual ove he eie expeime duaio as a fucio of heaig powe fo Figue 3.

Supplemeay Figue S5: Pediced vaiaio of coe empeaue as a fucio of ime fo diffee umbe of eigevalues cosideed fo he B () em i equaio (). Fo compaiso, he embedded hemocouple measueme is also show as a boke lie.

Supplemeay Figue S6: Covegece of he ifiie seies i equaio (7) wih espec o he umbe of eigevalues cosideed, showig slow seies covegece.

Supplemeay Figue S7: Compaiso of IR camea measued empeaue of a suface wih a kow suface empeaue as a fucio of ime, duig a empeaue amp fom oom empeaue o 6 C a a ae of 5 C/mi. IR daa ae show boh wih ad wihou calibaio i ed lies. The kow suface empeaue is also ploed i black fo compaiso. 3