Calculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
|
|
- Madeline Griffin
- 5 years ago
- Views:
Transcription
1 Ŕ periodica polytechnica Mechanical Engineering 54/ doi: /pp.me web: me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary layer beide wall with time-dependent heat tranfer coefficient ttila Fodor Received btract Thi paper propoe to invetigate the change in the temperature of external wall boundary layer of building when the heat tranfer coefficient reache it tationary tate in time exponentially. We eek the olution to the one-dimenional parabolic partial differential equation decribing the heat tranfer proce under pecial boundary condition. The earch for the olution originate from the olution of a Volterra integral equation of the econd kind. The kernel of the Volterra integral equation i lightly ingular therefore it olution i calculated numerically by one of the mot efficient collocation metho. Uing the Euler approach an iterative calculation algorithm i obtained, to be implemented through a programme written in the Maple computer algebra ytem. Change in the temperature of the external boundary of brick wall and wall inulated with polytyrene foam are calculated. The concluion i reached that the external temperature of the inulated wall matche the air temperature ooner than that of the brick wall. Keywor heat conduction lightly ingular Volterra integral equation numerical olution heat tranfer coefficient relaxation time ttila Fodor BME, H-1521 Budapet, Hungary fodor@pmmk.pte.hu 1 Introduction The heat tranfer and thermal conductivity of the wall boundary layer of building enure the temperature potential inide and outide the wall boundarie of building. The majority of the relevant model focu on temperature change at internal point in the wall boundary depending on the variou inulation layer. It i aumed that the temperature of the external wall boundary and the average temperature of the outide air urrounding the wall are the ame i.e. the outide temperature become contant and the temperature of the internal wall boundary i alo contant. Thi problem i not dicued by thi paper. We invetigate how quickly the external wall boundary adopt the temperature of the air boundary layer under a contant temperature of the internal wall boundary. Thi only require the examination of the implified model of the wall boundary preented in Fig. 1. In the mathematical model the wall i regarded a one of infinite thickne. It temperature i denoted by ux, t at poition x > and time t >. It i ufficient to deal with the one dimenional cae a it can be aumed that along the length of the wall the temperature i ditributed equally at each vertical croection. In accordance with the denotation u, t denote the temperature of the external boundary wall layer. The external layer of the wall i in contact with the ambient air. The internal layer of the wall mixe with the air of the room, aumed to have a contant temperature. The wall reit the incoming heat flow perpendicular to the wall therefore in the "air -wall" boundary layer the temperature of the air i different than the urface temperature of the wall. The temperature of the external wall boundary take ome time to adopt the average temperature of the outdoor air. denote the average temperature of the air flow. In winter it i typically < µ, t i.e. the outdoor air cool the wall. In ummer, t <, i.e. the outdoor air warm the wall. The model can be applied to both cae but implified criteria u, = i ued. The initial temperature t = ec of the wall layer i thu C and <, i.e. the air flow gradually warm the external wall boundary. Then, if t ten to infinity the temperature of the external wall boundary reache air tem- Time-dependent heat tranfer coefficient
2 Fig. 1. Simplified model of wall boundary Fig. 2. Curve of the α heat tranfer coefficient in the boundary layer with relaxation time T= 2, 5, 1, 15 ec downwar from top perature. Thi change will be calculated numerically. The thermal impact of the air on the wall i characterized by convective heat tranfer coefficient α. It i aumed that the thermal tranmittance coefficient αt of the wall change in time αt = α 1 e t T 1 exponentially [1] where α i the contant convective heat tranfer coefficient of the external ide, α = 23 W/m 2 K and T > i the o-called relaxation time with a dimenion inec. The change of the heat tranfer coefficient in time ipreented in Fig. 2. In Fig. 2 the graph how the change of the convective heat tranfercoefficient for relaxation time T = 2, 5, 1 and 15. T increae the graph go downwar i.e. the function reache a tate of balance later if T i higher. Time t required to reach the ame degree of balance i in direct proportion to the length of T. For example time t required to reach,95 time balance α i attained by olving equation, 95 = 1 e t T where the olution i t = T. Our objective i to determine by continual approache function u, t which how in time the change of temperature of the external wall boundary layer x =, taking into account that the air flow parallel to the wall modifie in time the αt heat tranfer factor of the wall-air boundary layer according to exponential equation preented in formula 1. The ubject of thi paper therefore belong to the topic of boundary layer theory. The partial equation i written up and will be olved uing the Volterra equation. In ection 2 of thi paper the partial differential equation determining temperature ux, t of the wall will be preented together with the neceary boundary condition. To eek the olution the Laplace tranformation of the olution in time t will be ued. To attain the reult we need to olve a Volterra integral equation of the econd kind. Section 3 will preent the formula determining the iterative olution of the obtained Volterra integral equation. Uing thee formula a Maple computer algebra programme will be written calculating the approaching value of the olution in an iterative way. Iue relating to the numerical preciion and tability of the approaching value can be conulted in the bibliography ee [1, 3] and [4]. Section 4 deal with two type of wall: brick wall and wall inulated with polytyrene. The neceary parameter are calculated for thee two cae and the numerical olution are written up. Finally a concluion will be drawn about the behaviour of the boundary layer of the two type of wall. 2 The mathematical model of the problem and deduction of the Volterra integral equation 2.1 Governing equation and olution procedure Temperature ux, t of the wall atifie the regular heat conduction equation 2 u x, t = a u x, t t x2 with the following boundary condition: u x, = u, t = x > ; t > 2.a 2.b 2.c 16 Per. Pol. Mech. Eng. ttila Fodor
3 t λ x x= u x, = α 1 e T t u, t 2.d where: λ i the thermal conductivity of the wall and the inulation material, W/mK, α i the tationary convective heat tranfer coefficient, W/m 2 K, a = ρc λ thermal diffuivity, m2 / if the temperature i not evenly ditributed it define the peed of equalization ρ denity, kg/m 3, c pecific heat capacity, J/kg K. Dynamic equation 2.d play a crucial role in the proce. The peed of the change in temperature perpendicularly impacting the wall i proportional to the difference in the temperature of the wall urface and the air and the proportion factor i heat tranfer coefficient α obtained by formula Step of deducing the Volterra integral equation Temperature u, t of the external wall boundary layer i obtained from Eq. 2.a-d and it i not neceary to calculate temperature ux, t x > of the internal point of the wall. We need to deduce from Eq. 2.a-d a o-called Volterra integral equation only containing function u, t, which ha a lightly ingular kernel of the econd kind. From the theory of Volterra integral equation of the econd kind the method of numerical olution i ued where the Euler approach i applied to calculate the value of the integral. Thi cheme i the implet verion of the collocation method. The integral equation to calculate temperature u, t of the wall-air boundary layer i derived through the following tep. working of thi problem Garbai [2] too. Step 1 Variable t, x and ux, t in Eq. 2.a-d are replaced tranformation X = x ; τ = t u x, t ; U X, t = T a T by new dimenionle variable X, τ and dependent variable UX,τ. Thee variable turn Eq. 2.a-d into the following impler verion: U X, τ = 2 t X U X, τ x= where = α at λ. U X, τ X > ;τ > 3.a X 2 U X, = U, τ = 3.b 3.c = 1 e τ U, τ 1 3.d In the following ection we will deal with the olution to Eq. 3.a-d but variable X, τ and UX, τ are replaced by the original variable x, t and ux, t. in Eq. 3.a-d boundary layer U,τ tarted from the initial τ = C and following the tranformation the mean air temperature went up to 1 C we will invetigate the curve along which the function of temperature U,τ increae from to 1. If we wih to return to the original temperature [ C] the calculated U, τ value nee to be multiplied by becaue u x, t = U X, τ fter the tranformation the variable value of τ how how many time it excee relaxation time T. U, τ = r mean the temperature that the boundary layer adopt at point of time t = r T, i.e. r time of relaxation time T. Step 2 To find olution ux,t to parabolic partial differential equation 3.a-d the Laplace tranformation i applied to variable t. The Laplace tranform of wx, for t of function x,t i defined by formula L t u x, t = e t u x, t dt = w x, The Laplace tranform of partial differential Eq. 3.a turn into L t t u x, t 2 = L t u x, t x2 w x, u x, = 2 w x, x2 w x, = 2 w x, 4 x2 ordinary differential equation where the uniformity for the derived Laplace tranform and initial condition 3.b ux, = are ued. Step 3 With regard to variable x Eq. 4 i an ordinary linear differential equation of the econd kind with a contant coefficient where i one of the parameter. It general olution i w x, = F1 e x + F2 e x where F1 and F2 are the dicretionary function of parameter. in the cae of x ux, t i zero according to 3.c therefore lim w x, = alo applie to the Laplace tranform of x wx, for t. Thi will only be atified if F1= becaue the wall i therefore w x, = F2 e x 5 Step 4 We derive the Laplace tranform of wx, from 5 by x and re-ubtitute wx, from Eq. 5 into the obtained equation x w x, = F2 e x x w x, = w x, 6 Time-dependent heat tranfer coefficient
4 Step 5 Since we are intereted in the temperature of the boundary layer ubtitution x= nee be performed in Eq. 6. x x= w x, = w, t i 1 e g = 1 + π i 1 j= t j+1 1 e g where F = L t f t and and H = L t h t. fter arranging the equation Step 6 If the Laplace tranformation i ued in variable t for boundaryeq. 3.d g t = 1 t 1 e g t 13 π x x= w x, = L t 1 e t u, t 8 i obtained where the Laplace tranformation equality wa ued L t 1 e t = = the left ide of Eq. 7 and 8 are the ame the right ide are alo equal w, = L t 1 e t u, t Both ide of the equation hould be divided by expreion and it hould be taken into account that w, = L t u, t L t u, t = L t 1 e t u, t 1 Step 7 Eq. 1 include the Laplace tranform L t u, t of temperature function u, t of the wall boundary layer and the Laplace tranform of function 1 e t u, t. In function Eq can be brought to zero if u, t i ubtituted by u, t = g t + 1 t 11 Baed on the linearity of Laplace tranformation and correlation 9 L t g t + 1 = 1 + L t 1 e t g t + 1 L t g t + 1 = L t 1 e t g t 12 Step 8 The invere Laplace tranformation hould be ued for both ide of Eq. 12 and the invere Laplace tranformation corre- lation L t 1 = 1 t hould be applied a well a the convolution formula written up with invere g t i = Volterra integral equation of the econd kind i obtained in which kernel K t, = 1 t i not interpreted at the = t end point of the integral interval. 3 Numerical olution to Volterra integral equation uing the Euler approach We have een that function g t = u, t 1 atifie the lightly ingular Volterra integral equation of the econd kind 13 where function ux, t i the olution to group of Eq. 3.a- 3.d. Let ee < t <H time interval and it ditribution to N>1 equal part B N t = < t 1 < t 2 < t 3 < < t i < t i+1 <... < t N = H where t i = iδ i the i-th diviion point and δ = H N i the pacing i =, 1, 2,..., N. We eek the approaching value g t 1, g t 2,..., g t N to olution gt of the Volterra integral Eq.13 at B N diviion point. We know the initial g=-1 value. It i aumed that approaching value g t 1, g t 2,..., g t i 1, are known. To obtain the value of g t i we need to break up the integral from the integral equation into the um of integral from partial interval B N g t i = t i i 1 j= 1 e g = π t j+1 1 e g 14 1 where it may be: i = 1, 2, 3,..., N. In integral from [, +1 ] partial interval multiplying factor f = 1 e g hould be approached by contant value f recorded in the left end point = while ingular kernel hould be integrated by in the given integral. Thi give u the following +1 1 e g = 1 e g +1 ti 1 e g 2 ti t i +1 1 ti = Per. Pol. Mech. Eng. ttila Fodor
5 for all j =, 1, 2,..., i 1. Therefore to obtain the approaching value for gt i the following recurive calculation formula i ued g t i = i 1 j= 1 e g ti t i On the right ide of recurive formula 16 value g t = 1, g t 1, g t 2,..., g t i 1 are preented therefore the value of function g can be calculated tep by tep. To calculate gt 2 i only neceary to know g t = 1. To calculate g t 2 the previouly calculated g t 1 i alo required. To calculate g t 3 both g t 1 and g t 2 previouly calculated are needed. The calculation of the next temperature figure alway ue the previouly obtained value of temperature. 4 Numerical reult for different wall urface There i a ingle parameter in Volterra integral Eq. 13 which can be calculated uing correlation = α at λ with the earlier mentioned data of the urface material of the wall. In the cae of a B3 brick wall λ=.647, ρ=146 and c=.88. If relaxation time T = 1 ec i ued in the calculation then =.7977 i obtained. For a B3 brick wall relaxation time T = 2 ec = i obtained. For inulating material polytyrene foam the ame data are λ=.47 with parameter ρ=15, c=1.46. If T = 1 then = while for T = 2 the parameter i = It can be thu tated that the value of coefficient in the Volterra integral Eq. 13 increae if better inulting material are ued relaxation time T i increaed. The following numerical and graphical reult were attained through the Maple computer algebra ytem. Maple function wa written for recurive correlation 16. Uing element g =-1, g 1, g 2,..., g i 1 in a lit it calculate the next g i element through formula Fig. 3. Surface temperature function u, t in the cae of B3 brick wall 1 C outdoor air temperature of C to 1 C outdoor air temperature =1 C. The horizontal line i the t-time axi meaured in ec and the vertical line i wall urface temperature u,t meaured in C. In Fig. 4 the graph for function u, t ha alo been drawn in a coordinate ytem where logarithm formula x = log 1 2 T t wa ued on the horizontal axi. Thi figure illutrate better the change of the temperature function for mall t time. Fig. 4. Function u,t with parameter =.8 and horizontal cale x = log 1 2 T t The approaching calculation are repeated for parameter correponding to =22 polytyrene inulation. The reult i preented in Fig. 5. iteration : = i 1 i e jδ g j i j δ i j 1 δ j= 17 Chooing parameter =.8 which correpon to brick wall B3 with relaxation time T = 1 [ec] the graph how the approaching value for urface temperature u, t = gt + 1. The Fig. 3 preent the curve along which the external layer of the brick wall increae from the initial temperature of C to ued Fig. 5. Surface temperature function u,t if polytyrene foam =22 i The graph for function u, t with parameter value =22 wa written up according to the Fig. 6 where the horizontal axi Time-dependent heat tranfer coefficient
6 i logarithmically caled. Fig. 6. Function u,t with parameter =22 and horizontal cale x = log 1 2 T t Drawing up the graph in Fig. 3 and 5 in a common coordinate ytem the difference between the two graph i highly viible a hown in Fig. 7. The higher graph how how the external wall boundary layer adapt to air temperature =1 C in the cae of a B3 brick wall with polytyrene inulation and the lower graph how the curve for a wall without inulation. Fig. 8. Surface temperature u, t hown on a logarithmic cale The logarithmic cale ued on the horizontal axi revere the place of the two graph compared to Fig. 7. In Fig. 8 the left graph applie to the B3 brick wall while the right graph applie to the polytyrene inulation. wall Fig. 7. Surface temperature u,t for inulated and uninulated B3 brick Fig. 6 how that the temperature of the external urface of the inulated wall almot immediately adopt the temperature of the outdoor air while thi proce take longer for an uninulated brick wall. The explanation for thi fact i that inulation doe not carry the outdoor temperature inide but form an obtacle to the preading of the heat o there i no equalization between the temperature of the external and internal layer in term of thermal tranmittance. The brick wall, however, tranmit the external volume of heat better and a kind of equalization proce tart between the indoor and outdoor temperature. a reult the temperature of the external urface of the B3 brick wall i impacted by indoor temperature C i.e. the brick wall carrie the volume of heat inide and reache outdoor temperature 1 C later than the external urface of the inulated wall. The graph of the temperature function preented in Fig. 4 and 6 were put into the ame coordinate ytem Fig. 8 where the horizontal axi i logarithmically caled. 5 Concluion Baed on Fig. 7 the following concluion can be drawn about the adaptation capability of the external wall urface: the temperature of the wall urface adapt more quickly in time to the air temperature if the wall i inulated. Thi phenomenon i explained by the fact that brick in an uninulated wall carry the heat inide and thu the temperature of the external wall urface adapt more lowly to the air temperature. In thi cae the indoor temperature influence the outdoor urface temperature. If the wall i inulated, however, the indoor temperature cannot impact the temperature of the external wall urface to the ame extent a the inulation provide a thermal barrier between the external and internal temperature and adopt the outdoor air temperature almot immediately. Reference 1 Becker N M, Bivin R L, Hu Y C, Murphy H D, Jr.White B, Wing G.M, Heat diffuion with time-dependent convective boundary condition, International Journal for numerical metho in energinering , Garbai L, Fénye T, The newet reult of heat conduction theory, Periodica polytechnica, Ser. Mech. Eng. 44 2, no. 2, Diogo T, Lima P, Rebelo M, Computational metho for a nonlinear Volterra integral equation, Proceeding of Hercma 25, Ori P, Product Integration for Volterra integral equation of the econd kind with weakly ingular kernel, Mathematic of Computation , no. 215, Per. Pol. Mech. Eng. ttila Fodor
7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationSolving Differential Equations by the Laplace Transform and by Numerical Methods
36CH_PHCalter_TechMath_95099 3//007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method OBJECTIVES When you have completed thi chapter, you hould be able to: Find the
More informationME 375 EXAM #1 Tuesday February 21, 2006
ME 375 EXAM #1 Tueday February 1, 006 Diviion Adam 11:30 / Savran :30 (circle one) Name Intruction (1) Thi i a cloed book examination, but you are allowed one 8.5x11 crib heet. () You have one hour to
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationLaplace Transformation
Univerity of Technology Electromechanical Department Energy Branch Advance Mathematic Laplace Tranformation nd Cla Lecture 6 Page of 7 Laplace Tranformation Definition Suppoe that f(t) i a piecewie continuou
More informationGiven the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is
EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationME 375 FINAL EXAM SOLUTIONS Friday December 17, 2004
ME 375 FINAL EXAM SOLUTIONS Friday December 7, 004 Diviion Adam 0:30 / Yao :30 (circle one) Name Intruction () Thi i a cloed book eamination, but you are allowed three 8.5 crib heet. () You have two hour
More informationSECTION x2 x > 0, t > 0, (8.19a)
SECTION 8.5 433 8.5 Application of aplace Tranform to Partial Differential Equation In Section 8.2 and 8.3 we illutrated the effective ue of aplace tranform in olving ordinary differential equation. The
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationPractice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions
For Quetion -6, rewrite the piecewie function uing tep function, ketch their graph, and find F () = Lf(t). 0 0 < t < 2. f(t) = (t 2 4) 2 < t In tep-function form, f(t) = u 2 (t 2 4) The graph i the olid
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationEvolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis
Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne
More informationME 375 FINAL EXAM Wednesday, May 6, 2009
ME 375 FINAL EXAM Wedneday, May 6, 9 Diviion Meckl :3 / Adam :3 (circle one) Name_ Intruction () Thi i a cloed book examination, but you are allowed three ingle-ided 8.5 crib heet. A calculator i NOT allowed.
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationCHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS
Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical
More informationJump condition at the boundary between a porous catalyst and a homogeneous fluid
From the SelectedWork of Francico J. Valde-Parada 2005 Jump condition at the boundary between a porou catalyt and a homogeneou fluid Francico J. Valde-Parada J. Alberto Ochoa-Tapia Available at: http://work.bepre.com/francico_j_valde_parada/12/
More informationRecent progress in fire-structure analysis
EJSE Special Iue: Selected Key Note paper from MDCMS 1 1t International Conference on Modern Deign, Contruction and Maintenance of Structure - Hanoi, Vietnam, December 2007 Recent progre in fire-tructure
More informationThe Hassenpflug Matrix Tensor Notation
The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of
More informationReading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions
Chapter 4 Laplace Tranform 4 Introduction Reading aignment: In thi chapter we will cover Section 4 45 4 Definition and the Laplace tranform of imple function Given f, a function of time, with value f(t
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world leading publiher of Open Acce book Built by cientit, for cientit 3,5 8,.7 M Open acce book available International author and editor Download Our author are among the 5 Countrie
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationChapter 13. Root Locus Introduction
Chapter 13 Root Locu 13.1 Introduction In the previou chapter we had a glimpe of controller deign iue through ome imple example. Obviouly when we have higher order ytem, uch imple deign technique will
More informationModeling of Transport and Reaction in a Catalytic Bed Using a Catalyst Particle Model.
Excerpt from the Proceeding of the COMSOL Conference 2010 Boton Modeling of Tranport and Reaction in a Catalytic Bed Uing a Catalyt Particle Model. F. Allain *,1, A.G. Dixon 1 1 Worceter Polytechnic Intitute
More informationA SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho
A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES Sanghyun Cho Abtract. We prove a implified verion of the Nah-Moer implicit function theorem in weighted Banach pace. We relax the
More informationDesign By Emulation (Indirect Method)
Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal
More informationMODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS
MODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS S. VAN BAARS Department of Science, Technology and Communication, Univerity of Luxembourg, Luxembourg ABSTRACT: In oil dynamic, the oil i
More informationSpot-on: Safe Fuel/Air Compression
Spot-on: Safe Fuel/Air Compreion Problem preented by Robert Hart and Kevin Hughe Veeder-Root Participant: Jeffrey Bank Joeph Fehribach Alitair Fitt John Ockendon Colin Pleae Don Schwendeman Burt Tilley
More informationHomework #7 Solution. Solutions: ΔP L Δω. Fig. 1
Homework #7 Solution Aignment:. through.6 Bergen & Vittal. M Solution: Modified Equation.6 becaue gen. peed not fed back * M (.0rad / MW ec)(00mw) rad /ec peed ( ) (60) 9.55r. p. m. 3600 ( 9.55) 3590.45r.
More informationEE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505)
EE 4443/5329 LAB 3: Control of Indutrial Sytem Simulation and Hardware Control (PID Deign) The Inverted Pendulum (ECP Sytem-Model: 505) Compiled by: Nitin Swamy Email: nwamy@lakehore.uta.edu Email: okuljaca@lakehore.uta.edu
More informationtwo equations that govern the motion of the fluid through some medium, like a pipe. These two equations are the
Fluid and Fluid Mechanic Fluid in motion Dynamic Equation of Continuity After having worked on fluid at ret we turn to a moving fluid To decribe a moving fluid we develop two equation that govern the motion
More informationReading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions
Chapter 4 Laplace Tranform 4 Introduction Reading aignment: In thi chapter we will cover Section 4 45 4 Definition and the Laplace tranform of imple function Given f, a function of time, with value f(t
More informationA Single Particle Thermal Model for Lithium Ion Batteries
A Single Particle Thermal Model for Lithium Ion Batterie R. Painter* 1, B. Berryhill 1, L. Sharpe 2 and S. Keith Hargrove 2 1 Civil Engineering, Tenneee State Univerity, Nahville, TN, USA 2 Mechanical
More informationChapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem
Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationNONISOTHERMAL OPERATION OF IDEAL REACTORS Plug Flow Reactor
NONISOTHERMAL OPERATION OF IDEAL REACTORS Plug Flow Reactor T o T T o T F o, Q o F T m,q m T m T m T mo Aumption: 1. Homogeneou Sytem 2. Single Reaction 3. Steady State Two type of problem: 1. Given deired
More informationPRESSURE WORK EFFECTS IN UNSTEADY CONVECTIVELY DRIVEN FLOW ALONG A VERTICAL PLATE
Proceeding of 3ICCHMT 3 rd International Conference on Computational Heat and Ma Tranfer May 6 3, 3, Banff, CANADA Paper Number 87 PRESSURE WORK EFFECTS IN UNSTEADY CONVECTIVELY DRIVEN FLOW ALONG A VERTICAL
More informationSource slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis
Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar
DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationConduction Heat transfer: Unsteady state
Conduction Heat tranfer: Unteady tate Chapter Objective For olving the ituation that Where temperature do not change with poition. In a imple lab geometry where temperature vary alo with poition. Near
More informationHyperbolic Partial Differential Equations
Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a
More informationMath Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK
ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI
More informationNew bounds for Morse clusters
New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität
More informationConstant Force: Projectile Motion
Contant Force: Projectile Motion Abtract In thi lab, you will launch an object with a pecific initial velocity (magnitude and direction) and determine the angle at which the range i a maximum. Other tak,
More informationThermal Stress in a Half-Space with Mixed Boundary Conditions due to Time Dependent Heat Source
IOSR Journal of Mathematic (IOSR-JM) e-issn: 2278-5728, p-issn: 239-765X Volume, Iue 6 Ver V (Nov - Dec 205), PP 9-25 wwwiorjournalorg Thermal Stre in a Half-Space with Mixed Boundary Condition due to
More informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More information5.5 Application of Frequency Response: Signal Filters
44 Dynamic Sytem Second order lowpa filter having tranfer function H()=H ()H () u H () H () y Firt order lowpa filter Figure 5.5: Contruction of a econd order low-pa filter by combining two firt order
More informationDigital Control System
Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)
More informationConvective Heat Transfer
Convective Heat Tranfer Example 1. Melt Spinning of Polymer fiber 2. Heat tranfer in a Condener 3. Temperature control of a Re-entry vehicle Fiber pinning The fiber pinning proce preent a unique engineering
More informationAdvanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment
Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,
More informationNonequilibrium Distribution Function Theory of Many-Particle Effects in the Reversible Reactions of the Type A+B C+B
986 Bull Korean Chem Soc 5, Vol 6, No Jinuk Lee et al Nonuilibrium Ditribution Function Theory of Many-Particle Effect in the Reverible Reaction of the Type AB CB Jinuk Lee, Jeik Uhm, Woojin Lee, Sangyoub
More informationAn estimation approach for autotuning of event-based PI control systems
Acta de la XXXIX Jornada de Automática, Badajoz, 5-7 de Septiembre de 08 An etimation approach for autotuning of event-baed PI control ytem Joé Sánchez Moreno, María Guinaldo Loada, Sebatián Dormido Departamento
More informationRao Transforms: A New Approach to Integral and Differential Equations
Rao Tranform: A New Approach to Integral and Differential Equation Dr. Muralidhara SubbaRao (Rao) SUNY at Stony Brook, murali@ece.unyb.edu,, rao@integralreearch.net Rao Tranform (RT) provide a brand new
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationOnline Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat
Online Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat Thi Online Appendix contain the proof of our reult for the undicounted limit dicued in Section 2 of the paper,
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationSolutions to homework #10
Solution to homework #0 Problem 7..3 Compute 6 e 3 t t t 8. The firt tep i to ue the linearity of the Laplace tranform to ditribute the tranform over the um and pull the contant factor outide the tranform.
More informationEmittance limitations due to collective effects for the TOTEM beams
LHC Project ote 45 June 0, 004 Elia.Metral@cern.ch Andre.Verdier@cern.ch Emittance limitation due to collective effect for the TOTEM beam E. Métral and A. Verdier, AB-ABP, CER Keyword: TOTEM, collective
More informationLinear Motion, Speed & Velocity
Add Important Linear Motion, Speed & Velocity Page: 136 Linear Motion, Speed & Velocity NGSS Standard: N/A MA Curriculum Framework (006): 1.1, 1. AP Phyic 1 Learning Objective: 3.A.1.1, 3.A.1.3 Knowledge/Undertanding
More informationMulti-dimensional Fuzzy Euler Approximation
Mathematica Aeterna, Vol 7, 2017, no 2, 163-176 Multi-dimenional Fuzzy Euler Approximation Yangyang Hao College of Mathematic and Information Science Hebei Univerity, Baoding 071002, China hdhyywa@163com
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4
More informationTo appear in International Journal of Numerical Methods in Fluids in Stability analysis of numerical interface conditions in uid-structure therm
To appear in International Journal of Numerical Method in Fluid in 997. Stability analyi of numerical interface condition in uid-tructure thermal analyi M. B. Gile Oxford Univerity Computing Laboratory
More informationDetermination of Flow Resistance Coefficients Due to Shrubs and Woody Vegetation
ERDC/CL CETN-VIII-3 December 000 Determination of Flow Reitance Coefficient Due to hrub and Woody Vegetation by Ronald R. Copeland PURPOE: The purpoe of thi Technical Note i to tranmit reult of an experimental
More informationSampling and the Discrete Fourier Transform
Sampling and the Dicrete Fourier Tranform Sampling Method Sampling i mot commonly done with two device, the ample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquire a CT ignal at
More informationLinearteam tech paper. The analysis of fourth-order state variable filter and it s application to Linkwitz- Riley filters
Linearteam tech paper The analyi of fourth-order tate variable filter and it application to Linkwitz- iley filter Janne honen 5.. TBLE OF CONTENTS. NTOCTON.... FOTH-OE LNWTZ-LEY (L TNSFE FNCTON.... TNSFE
More informationThermal Resistance Measurements and Thermal Transient Analysis of Power Chip Slug-Up and Slug-Down Mounted on HDI Substrate
Intl Journal of Microcircuit and Electronic Packaging Thermal Reitance Meaurement and Thermal Tranient Analyi of Power Chip Slug-Up and Slug-Down Mounted on HDI Subtrate Claudio Sartori Magneti Marelli
More informationLaplace Adomian Decomposition Method for Solving the Nonlinear Volterra Integral Equation with Weakly Kernels
Studie in Nonlinear Science (4): 9-4, ISSN -9 IDOSI Publication, Laplace Adomian Decompoition Method for Solving the Nonlinear Volterra Integral Equation with Weakly Kernel F.A. Hendi Department of Mathematic
More informationApproximate Analytical Solution for Quadratic Riccati Differential Equation
Iranian J. of Numerical Analyi and Optimization Vol 3, No. 2, 2013), pp 21-31 Approximate Analytical Solution for Quadratic Riccati Differential Equation H. Aminikhah Abtract In thi paper, we introduce
More informationLTV System Modelling
Helinki Univerit of Technolog S-72.333 Potgraduate Coure in Radiocommunication Fall 2000 LTV Stem Modelling Heikki Lorentz Sonera Entrum O heikki.lorentz@onera.fi Januar 23 rd 200 Content. Introduction
More informationFactor Analysis with Poisson Output
Factor Analyi with Poion Output Gopal Santhanam Byron Yu Krihna V. Shenoy, Department of Electrical Engineering, Neurocience Program Stanford Univerity Stanford, CA 94305, USA {gopal,byronyu,henoy}@tanford.edu
More informationNCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011
NCAAPMT Calculu Challenge 011 01 Challenge #3 Due: October 6, 011 A Model of Traffic Flow Everyone ha at ome time been on a multi-lane highway and encountered road contruction that required the traffic
More informationChapter 7. Root Locus Analysis
Chapter 7 Root Locu Analyi jw + KGH ( ) GH ( ) - K 0 z O 4 p 2 p 3 p Root Locu Analyi The root of the cloed-loop characteritic equation define the ytem characteritic repone. Their location in the complex
More information1. Basic introduction to electromagnetic field. wave properties and particulate properties.
Lecture Baic Radiometric Quantitie. The Beer-Bouguer-Lambert law. Concept of extinction cattering plu aborption and emiion. Schwarzchild equation. Objective:. Baic introduction to electromagnetic field:
More information696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of
Vol 12 No 7, July 2003 cfl 2003 Chin. Phy. Soc. 1009-1963/2003/12(07)/0695-05 Chinee Phyic and IOP Publihing Ltd Lie ymmetrie and conerved quantitie of controllable nonholonomic dynamical ytem Fu Jing-Li(ΛΠ±)
More informationLecture 15 - Current. A Puzzle... Advanced Section: Image Charge for Spheres. Image Charge for a Grounded Spherical Shell
Lecture 15 - Current Puzzle... Suppoe an infinite grounded conducting plane lie at z = 0. charge q i located at a height h above the conducting plane. Show in three different way that the potential below
More informationCONTROL SYSTEMS. Chapter 2 : Block Diagram & Signal Flow Graphs GATE Objective & Numerical Type Questions
ONTOL SYSTEMS hapter : Bloc Diagram & Signal Flow Graph GATE Objective & Numerical Type Quetion Quetion 6 [Practice Boo] [GATE E 994 IIT-Kharagpur : 5 Mar] educe the ignal flow graph hown in figure below,
More informationAutomatic Control Systems. Part III: Root Locus Technique
www.pdhcenter.com PDH Coure E40 www.pdhonline.org Automatic Control Sytem Part III: Root Locu Technique By Shih-Min Hu, Ph.D., P.E. Page of 30 www.pdhcenter.com PDH Coure E40 www.pdhonline.org VI. Root
More informationThe continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.
1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker
More informationCONTROL SYSTEMS. Chapter 5 : Root Locus Diagram. GATE Objective & Numerical Type Solutions. The transfer function of a closed loop system is
CONTROL SYSTEMS Chapter 5 : Root Locu Diagram GATE Objective & Numerical Type Solution Quetion 1 [Work Book] [GATE EC 199 IISc-Bangalore : Mark] The tranfer function of a cloed loop ytem i T () where i
More information722 Chen Xiang-wei et al. Vol. 9 r i and _r i are repectively the poition vector and the velocity vector of the i-th particle and R i = dm i dt u i; (
Volume 9, Number 10 October, 2000 1009-1963/2000/09(10)/0721-05 CHINESE PHYSICS cfl 2000 Chin. Phy. Soc. PERTURBATION TO THE SYMMETRIES AND ADIABATIC INVARIANTS OF HOLONOMIC VARIABLE MASS SYSTEMS * Chen
More informationSMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD
SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem
More informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
More informationOptimization model in Input output analysis and computable general. equilibrium by using multiple criteria non-linear programming.
Optimization model in Input output analyi and computable general equilibrium by uing multiple criteria non-linear programming Jing He * Intitute of ytem cience, cademy of Mathematic and ytem cience Chinee
More informationHalliday/Resnick/Walker 7e Chapter 6
HRW 7e Chapter 6 Page of Halliday/Renick/Walker 7e Chapter 6 3. We do not conider the poibility that the bureau might tip, and treat thi a a purely horizontal motion problem (with the peron puh F in the
More informationLecture #9 Continuous time filter
Lecture #9 Continuou time filter Oliver Faut December 5, 2006 Content Review. Motivation......................................... 2 2 Filter pecification 2 2. Low pa..........................................
More informationin a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim
Quantization of electromagnetic eld in a circular cylindrical cavity K. Kakazu Department of Phyic, Univerity of the Ryukyu, Okinawa 903-0, Japan Y. S. Kim Department of Phyic, Univerity of Maryland, College
More informationThe Laplace Transform , Haynes Miller and Jeremy Orloff
The Laplace Tranform 8.3, Hayne Miller and Jeremy Orloff Laplace tranform baic: introduction An operator take a function a input and output another function. A tranform doe the ame thing with the added
More informationImproving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization
1976 MONTHLY WEATHER REVIEW VOLUME 15 Improving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization PETER LYNCH Met Éireann, Dublin, Ireland DOMINIQUE GIARD CNRM/GMAP, Météo-France,
More information