ESTIMATION OF HEAT SOURCES WITHIN TWO DIMENSIONAL SHAPED BODIES
|
|
- Shanon Gillian Summers
- 5 years ago
- Views:
Transcription
1 Proceedigs of 3ICIPE 99 3rd It. Coferece o Iverse Problems i Egieerig Jue 13-18, 1999, Port Ludlow, Washito, USA HT11 ESTIMATION OF HEAT SOURCES WITHIN TWO DIMENSIONAL SHAPED BODIES Abou Khachfe R., Jary Y. Laboratoire de Thermociétique de l ISITEM, UMR CNRS 667 Rue Christia Pauc BP 964, 4436 Nates Cedex 3 jary@isitem.uiv-ates.fr ABSTRACT A iterative algorithm based o the cojugate gradiet method is developed i order to estimate simultaeously the spatial locatio ad the stregth of heat sources withi two dimesioal shaped bodies. The case of puctual sources, together with temperature sesors located o the boudary of the body, is studied. The iterative regularizatio pricilple is used to avoid the amplificatio of the measuremet errors o the computed solutio. NOMENCLATURE A streght of the source c heat capacity J(S) fuctioal equatio (8) J S gradiet of J(S) h heat trasfer coefficiet i ormal uit vector o Γ i T temperature (x y ) locatio of the source Γ i boudary surface i λ thermal coductivity δt sesitivity fuctio ψ adjoit fuctio ρ desity INTRODUCTION The umerical resolutio of the Direct Heat Coductio Problem (DHCP) requires the kowledge of a set of data which ivolves : the geometrical data (shape of the body), the medium data (heat capacity, desity, thermal coductivity, heat trasfer coefficiet), the source data (heat flux desity, fixed temperature o the boudary, heat sources withi the body), the iitial coditio (spatial distributio of the temperature field at the iitial time). Whe some part of this set of data is ukow, the DHCP solutio caot be computed, the umerical resolutio of the Iverse Heat Coductio Problems (IHCP) aims to determie the ukow part of these data from additioal iformatio. These ew data are usually give by temperature sesors which are located o the boudary or iside the body. Durig the last decade, due to the tremedous advacemet i scietific computatio, a icreasig attetio has bee devoted to the techiques for solvig multidimesioal IHCP. Numerical resolutios of IHCP are well kow to be ill - coditioed, therefore some regularizig process has to be cosidered to avoid umerical istabilities o the computed solutios. Several approaches are available (Alifaov, 1994), (Jary et al, 1991), (Beck, 1977). A geeral approach cosists i solvig IHCP with umerical optimizatio tools like the cojugate gradiet method. Combied with the use of a Fiite Elemet Library, It is a powerful method which offers a wide field of practical applicatios i iverse thermal aalysis. I this work, some results are preseted o the estimatio of heat sources withi two dimesioal shaped bodies. This problem has bee cosidered by several authors (Silva Neto et al, 1992), (Le Niliot, 1998),(Park H.M. et al, 1999) especially i the case of puctual sources ad assumig the spatial locatio is kow. I this paper it is show how the simultaeous determiatio of the spatial locatio ad the stregth of poctual heat sources ca be achieved from temperature histories delivered by sesors located o the boudary of the body. Numerical experimets show that the ifluece of the error measuremets 1 Copyright 1999 by ASME
2 o the computed solutio ca be strogly reduced accordig to the iterative regularizatio priciple. DIRECT HEAT CONDUCTION PROBLEM (DHCP) The heat coductio process ivolvig heat source is cosidered withi a arbitrary shaped domai with the boudary Γ 1 [Γ 2. The set of equatios of the DHCP are writte i the form of a parabolic problem where P(t) is the itesity (W =m) ad ω is a parameter used to simulate a puctual heat source, it is chose equal to the mesh size. We cosider ad (x y ) is the source locatio. A(t) = P(t) πω 2 (7) ρc T ; λ T = S(x y t) i t > (1) λ T 1 + ht = ht e o Γ 1 t > (2) λ T 2 = oγ 2 t > (3) T = T i t = (4) Solvig the DHCP cosists i the determiatio of the temperature field T (x y t) o the domai, at each time t, whe the medium data ρc, λ, the heat trasfer coefficiet h, the exteral temperature T e, the heat source S(x y t) ad the iitial field T, are assumed to be kow. I this study, the heat source is cosidered with the followig form S(x y t)=a(t):f(x y;x y ) Fd = 1 where A ad F are respectively the stregth ad the spatial distributio of the source. The goal is to determie the stregth A(t), ad the locatio (x y ) for a prescribed spatial distributio F. (5) FORMULATION OF THE IHCP Additioal iformatio are provided by sesors located withi the spatial domai or o its boudary, o the time iterval [ t f ]. Let Y m (t) (m = 1 ::) be the temperature histories delivered at the locatio (x m y m ). Solvig the IHCP cosists i the determiatio of the source S which matches the solutio of the DHCP T (x m y m t;s) with the data Y m (t). Because of illposedess, the IHCP is solved i the least square sese, by miimizig the residual fuctioal J(S) = 1 2 t f jt (x m y m t;s) ;Y m (t)j 2 dt (8) m=1 with T (x y t;s) solutio of equatios (1)-(4) MINIMIATION OF THE FUNCTIONAL J(S) Miimizatio of J(S) is achieved by usig the cojugate gradiet method (CGM). Let S be the ukow vector to be determied. S T =[x y A k (k= ::Nt) ] (9) where Nt is the umber of time step. The (CGM) is iterative, at each iteratio, the previous estimate S is corrected by S +1 = S ; γ d (1) where the search directio d ad the step legth γ are determied i order to have J(S +1 ) < J(S ) (11) PONCTUAL HEAT SOURCE The spatial distributio of the heat source to be determied is described by a gaussia fuctio S(x y t) = P(t) πω 2 exp (x ; ; x ) 2 ; (y ; y ) 2 ω 2 (6) The vector [d x d y d Ak (k= ::Nt)], is chose accordig to the CGM rules d = JS + β ;1 J S β = β = < J S J S ; J S;1 > kjs k 2 (12) 2 Copyright 1999 by ASME
3 where J S is the gradiet of the J(S) at the iteratio J S =[J x J y J A k (k= ::Nt) ] (13) k:k is the vector orm associated to the scalar product <: :> The step legth γ is give by γ = m=1 t f δt (x m y m t;s )[T (x m y m t;s ) ;Y m (t)]dt [δt (x m y m t;s )] 2 (14) where δt is the sesitivity fuctio which results of the variatio δs of the source. SENSITIVITY EQUATIONS Solvig the sesitivity equatios cosists i the determiatio of the temperature variatio δt which results of a variatio δs of the source. The liearity of the equatios (1)-(4) leads to I order to get J S (x y t;s), equatio (19) has to be put i the form of equatio (2). This is doe by makig statioary the Lagragia associated to the optimizatio problem, equatio (8) LAGRANGIAN AND ADJOINT EQUATIONS Let us itroduce the Lagragia L (T S ψ) t f 1 2 m=1 t f + L((T S ψ)) = (T (x y t;s) ;Y m (t)) 2 δ(x ; x m )δ(y ; y m )dtd (ρc T ; λ T ; S)ψdtd (21) where T S ψ ca be idepedet fuctios. The Lagragia multiplier ψ is, as the fuctio T, a fuctio of x y, ad t. ψ beig fixed, the differetial of L satisfies δl = T δt + δs (22) S δs = FδA + A ρc δt ; λ δt = δs! F δx F + δy i t > x y (15) Let us take ψ solutio of the followig equatio (so called adjoit equatio) δt = 8 δt (23) T the by takig T solutio of the equatios (1)-(4), it comes λ δt 1 + hδt = oγ 1 t > (16) δj = δl = δs (24) S λ δt 2 = oγ 2 t > (17) δt = i t = (18) ad the associated variatio of the fuctioal J(S) is t f m δj(s) = (T (x y t;s) ;Y m (t))δ(x ; x m )δ(y ; y m )δtdtd (19) where δ(:) is the Dirac impulse. By defiitio, J S (x y t;s) satisfies t f δj(s) = J S (x y t;s)δsdtd (2) where S δs = ; t f t f = ; ψδsdtd ψfδadtd ; t f ; t f ψa F y δy dtd ψa F x δx dtd (25) ad the gradiet compoets J S (x y t;s) are easily extracted from equatio (25), J A (t) =; ψf d (26) t f J x = ; ψa F dtd (27) x 3 Copyright 1999 by ASME
4 t f J y = ; ψa F dtd (28) y λ ψ 1 + hψ = oγ 1 t > (34) It ca be oted that J A is a fuctio of time, like the fuctio A(t), J x ad J y are real umbers like x ad y. Let us develop equatio (23). From the defiitio of the Lagragia, (21), we have t f m=1 T δt = (T (x y t;s) ;Y m (t))δ(x ; x m )δ(y ; y m )δtdtd t f + itegratio by parts give t f λ δt ψd = (ρc δt ; λ δt )ψdtd δt ψdt =[ψδt ] t=t t f f t= ; (29) ψ δtdt (3) (λ Γ1 ψ + hψ)δtdγ 1 + λ 1 Γ2 ψ δtdγ λ ψ δtd (31) the with the boudary ad iitial coditios (16), (17), (18), equatio (29) ca be writte i the ew form t f + T δt = m=1 t f (;ρc ψ ; λ ψ)δt dtd (T (x y t;s) ;Y m (t))δ(x ; x m )δ(y ; y m )δt dtd + Γ2 t f Γ1 t f (λ ψ 1 + hψ)δt dtdγ 1 + λ ψ δtdtdγ 2 + ψ(t = t f )δt d 2 (32) Therefore equatio (23) is satisfied by takig the lagrage multiplier ψ solutio of the followig set of equatios m=1 ;ρc ψ ; λ ψ = (T (x y t;s) ;Y m (t))δ(x ; x m )δ(y ; y m ) i t > (33) λ ψ 2 = oγ 2 t > (35) ψ = i t = t f (36) Equatios (33)-(36) are liear, the solutio of which ψ is computed backward i time o the iterval [ t f ] NUMERICAL ALGORITHM Each iteratio of the CGM icludes the umerical resolutio of three distict parabolic set of equatios. The iterative algorithm has the followig structure a) Choose a iitial guess S b) Solve the direct problem, equatios (1)-(4), to compute T (x m y m t;s ), ad J(S ), equatio (8) c) Solve the adjoit problem,equatios (33-36) to compute the vector J S, equatios (26-28) d) Solve the sesitivity problem, equatios (15-18), to compute γ, equatio (14), ad the ew vector S +1 equatio (1) e) Stoppig coditio if J(S +1 ) > J(S ) or J(S ) < ε else set = + 1gotob) f) Ed go to f) ε is defied accordig to the iterative regularizatio priciple which is illustrated i the followig sectio. NUMERICAL EXPERIMENTS The spatial domai cosidered is a square, see figure (1), L = :5m. The values of the thermophysical parameters are : λ = :17W =K:m, a = 1:21 ;7 m 2 =s, h = 1W=m 2 :K, the iitial temperature T = 293K, ad the exteral temperature T e = 293K. A regular mesh with = 121 odes is cosidered. ω = 6:1 ;3 m. The parabolic solver of the Fiite Elemet Library MODULEF is used to compute the solutio of the DHCP. Two cases have bee studied : case#1: x = :1m y = :4m case#2: x = :3m y = :1m The temperature resposes delivered by four sesors located at the middle of each side of the domai,to a time varyig stregth A(t) of the source, are show o figures (2) ad (3). 4 Copyright 1999 by ASME
5 y Γ 1 Source S Sesor 1 Γ Case#1 Iitial 1 locatio Γ 2 S 2 Γ 1 S 3 h Case#2 h S4 x h Temperature ( o C) 3 25 S 1 S 2 S 3 S 4 Figure 1. DOMAIN S 1 S 2 S 3 S 4 Figure 3. TEMPERATURE RESPONSES : case#2 Temperature ( o C) x or y (m) σ =.5 y x Figure 2. TEMPERATURE RESPONSES : case#1 These resposes are corrupted by addig a zero mea gaussia oise, to build the simulated additioal data Y m (t) required to solve the IHCP Iteratio umber Figure 4. IDENTIFICATION OF THE LOCATION : case#1, σ = :5 Y k m = T k m exact + σek m m = 1 ::: k = :::Nt (37) where σ is the stadard deviatio of the oise (assumed to be idetical for each sesor) ad e k m is the ormally distributed radum umber. NUMERICAL RESULTS Case#1 The iitial guess of the source locatio is take at the ceter of the square x = :25m y= :25m. The iitial guess of the stregth of the source is take equal to zero. A alterate iterative research process of the cojugate gradiet has bee performed : istead of modifyig the locatio ad the stregth at 5 Copyright 1999 by ASME
6 A(t) W/m Exact solutio Estimated solutio at = 242 σ =.5 A(t) W/m Exact solutio Estimated solutio at = 155 Estimated solutio at = 186 σ = Figure 5. IDENTIFICATION OF THE STRENGTH : case#1, σ = :5 Figure 7. IDENTIFICATION OF THE STRENGTH : case#1, σ = :2.5.4 σ = σ=.2 σ=.5 x or y (m).3.2 y x J(S ) Iteratio umber Iteratio umber Figure 6. IDENTIFICATION OF THE LOCATION : case#1, σ = :2 Figure 8. MINIMIATION OF J(S ): case#1 each iteratio, oe iteratio is used to correct the locatio keepig the stregth umodified, the followig iteratio is used to correct the stregth while the locatio is hold, ad so o. O figure (4) ad (6) are show the computed locatios resultig of this miimizatio process, with a oise level σ = :5 ad σ = :2 respectively. Figure (8) illustrates the covergece of the sequece J(S ) : asymptotic values J as are observed, they are directly re- lated to the oise level σ as idicate i Table 1. The fial umber of iteratios is chose accordig to the followig equatio J as = J(S )=Nt::σ 2 = ε (38) Whe the iterative process is stopped at the iteratio umber, 6 Copyright 1999 by ASME
7 Table 1. RESULTS OF IDENTIFICATION OF (x, y ), case#1 σ ε J as x y Table 2. RESULTS OF IDENTIFICATION OF (x, y ), case#2 σ ε J as x y A(t) W/m Exact solutio Estimated solutio at = 63 σ = σ =.5 Figure 1. IDENTIFICATION OF THE STRENGTH : case#2, σ = :5 x or y (m) x y Iteratio umber Figure 9. IDENTIFICATION OF THE LOCATION : case#2, σ = :5 x or y (m) x σ =.2 y o istability occur o the computed solutio as it is show o figure (5). Figure (7) illustrates the iterative regularizatio priciple : - with σ = :2 ad = 155, the computed stregth A(t) is stable, close to the exact values; if the iterative process is cotiued, for example = 186, oscillatios appear. - with a lower oise level σ = :5, the fial umber = 242 is greater ad the computed stregth is more accurate figure (5) Case#2 The same iitial guess is cosidered. Figures (9) ad (11) illustrates the covergece of the computed source locatios to the exact values. Figures (1) ad (12) show the ifluece of the oise level o the computed stregth A(t). Results are similar to Iteratio umber Figure 11. IDENTIFICATION OF THE LOCATION : case#2, σ = :2 case#1. It ca be observed that the fial umbers of iteratios depeds o the level oise, like i the previous case. The ifluece of the sesor locatios must be also oted. The distace betwee the earest sesor ad the source is lower tha i case#1, the for σ = :5, istead of = 242 oly = 63 iteratios are required, ad for σ=.2, istead of = 155, = 59 are sufficiet table 2. 7 Copyright 1999 by ASME
8 described by a stiff gaussia fuctio i order to approximate a puctual distributio. The iterative regularizatio priciple together with the use of the cojugate gradiet method has bee illustrated. To avoid istabilities o the computed stregth, the iterative process of miimizatio has to be stopped whe the residual fuctioal reaches a critical value which depeds o the oise level. A(t) W/m Exact solutio Estimated solutio at = 59 Estimated solutio at = 78 σ = Figure 12. IDENTIFICATION OF THE STRENGTH : case#2, σ = : σ=.2 σ=.5 REFERENCES Alifaov O.M.,Iverse Heat trasfert problems Spriger- Verlag,1994 Jary Y., Ozisik M.N., Bardo J.P., A geeral optimizatio method usig adjoit equatio for solvig multidimesioal iverse heat coductio, It. J. Heat Mass Tras.,34 (1991) Beck J. V., Arold K.J., Parameter Estimatio i Egieerig Sciece, Joh Whiley ad Sos, New York, 1977 Silva Neto A.J., Ozisik M.N., Two-dimesioal iverse heat coductio problem of estimatig the time-varyig stregth of a lie heat source, Joural of Applied physics 71 (11) (1992) Le Niliot C., Gallet P., Ifrared thermography applied to the resolutio of iverse heat coductio problems: recovery of heat lie sources ad boudary coditios, Revue Geerale de la Thermique,37 (1998), Park H.M., Chug O.Y., Lee J.H., O the solutio of iverse heat trasfer problem usig the Karhue-Loève Galerki method, It. J. Heat Mass Tras.,42 (1999) J(S ) Iteratio umber Figure 13. MINIMIATION OF J(S ): case#2 CONCLUSION The iverse heat coductio problem which cosists i the determiatio of the locatio ad simultaeously the time varyig stregth of a heat poctual source withi a two dimesioal arbitrary shaped body, has bee computed by the cojugate gradiet method, usig a parabolic solver of the Fiite Elemet Library MODULEF. The spatial distributio of the source has bee 8 Copyright 1999 by ASME
Monte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem
Australia Joural of Basic Applied Scieces, 5(): 097-05, 0 ISSN 99-878 Mote Carlo Optimizatio to Solve a Two-Dimesioal Iverse Heat Coductio Problem M Ebrahimi Departmet of Mathematics, Karaj Brach, Islamic
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationThermal characterization of an insulating material through a tri-layer transient method
Advaced Sprig School «Thermal Measuremets & Iverse techiques», Domaie de Fraço, Biarritz, March -6 05 T6 Thermal characterizatio of a isulatig material through a tri-layer trasiet method V. Félix, Y. Jaot,
More informationChapter 2 The Monte Carlo Method
Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful
More informationPC5215 Numerical Recipes with Applications - Review Problems
PC55 Numerical Recipes with Applicatios - Review Problems Give the IEEE 754 sigle precisio bit patter (biary or he format) of the followig umbers: 0 0 05 00 0 00 Note that it has 8 bits for the epoet,
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationChapter 4 : Laplace Transform
4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic
More informationH05 A 3-D INVERSE PROBLEM IN ESTIMATING THE TIME-DEPENDENT HEAT TRANSFER COEFFICIENTS FOR PLATE FINS
Proceedigs of te 5 t Iteratioal Coferece o Iverse Problems i Egieerig: Teory ad Practice, Cambridge, UK, 11-15 t July 005 A 3-D INVERSE PROBLEM IN ESTIMATING THE TIME-DEPENDENT HEAT TRANSFER COEFFICIENTS
More informationCO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS
CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,
More informationIntroduction to Optimization Techniques. How to Solve Equations
Itroductio to Optimizatio Techiques How to Solve Equatios Iterative Methods of Optimizatio Iterative methods of optimizatio Solutio of the oliear equatios resultig form a optimizatio problem is usually
More informationECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015
ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5
CS434a/54a: Patter Recogitio Prof. Olga Veksler Lecture 5 Today Itroductio to parameter estimatio Two methods for parameter estimatio Maimum Likelihood Estimatio Bayesia Estimatio Itroducto Bayesia Decisio
More informationLinear regression. Daniel Hsu (COMS 4771) (y i x T i β)2 2πσ. 2 2σ 2. 1 n. (x T i β y i ) 2. 1 ˆβ arg min. β R n d
Liear regressio Daiel Hsu (COMS 477) Maximum likelihood estimatio Oe of the simplest liear regressio models is the followig: (X, Y ),..., (X, Y ), (X, Y ) are iid radom pairs takig values i R d R, ad Y
More informationOutline. Linear regression. Regularization functions. Polynomial curve fitting. Stochastic gradient descent for regression. MLE for regression
REGRESSION 1 Outlie Liear regressio Regularizatio fuctios Polyomial curve fittig Stochastic gradiet descet for regressio MLE for regressio Step-wise forward regressio Regressio methods Statistical techiques
More informationa b c d e f g h Supplementary Information
Supplemetary Iformatio a b c d e f g h Supplemetary Figure S STM images show that Dark patters are frequetly preset ad ted to accumulate. (a) mv, pa, m ; (b) mv, pa, m ; (c) mv, pa, m ; (d) mv, pa, m ;
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationPrinciple Of Superposition
ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Trasform The -trasform is a very importat tool i describig ad aalyig digital systems. It offers the techiques for digital filter desig ad frequecy aalysis of digital sigals. Defiitio of -trasform:
More informationMCT242: Electronic Instrumentation Lecture 2: Instrumentation Definitions
Faculty of Egieerig MCT242: Electroic Istrumetatio Lecture 2: Istrumetatio Defiitios Overview Measuremet Error Accuracy Precisio ad Mea Resolutio Mea Variace ad Stadard deviatio Fiesse Sesitivity Rage
More informationParameter Inversion Model for Two Dimensional Parabolic Equation Using Levenberg-Marquardt Method
Iteratioal Cogress o Evirometal Modellig ad Software Brigham Youg Uiversity BYU ScholarsArchive 5th Iteratioal Cogress o Evirometal Modellig ad Software - Ottawa, Otario, Caada - July Jul st, : AM Parameter
More informationOrthogonal Gaussian Filters for Signal Processing
Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationDimension-free PAC-Bayesian bounds for the estimation of the mean of a random vector
Dimesio-free PAC-Bayesia bouds for the estimatio of the mea of a radom vector Olivier Catoi CREST CNRS UMR 9194 Uiversité Paris Saclay olivier.catoi@esae.fr Ilaria Giulii Laboratoire de Probabilités et
More informationMETHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS
Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,
More informationA Genetic Algorithm for Solving General System of Equations
A Geetic Algorithm for Solvig Geeral System of Equatios Győző Molárka, Edit Miletics Departmet of Mathematics, Szécheyi Istvá Uiversity, Győr, Hugary molarka@sze.hu, miletics@sze.hu Abstract: For solvig
More informationData Assimilation. Alan O Neill University of Reading, UK
Data Assimilatio Ala O Neill Uiversity of Readig, UK he Kalma Filter Kalma Filter (expesive Use model equatios to propagate B forward i time. B B(t Aalysis step as i OI Evolutio of Covariace Matrices (
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationwavelet collocation method for solving integro-differential equation.
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More information1 Inferential Methods for Correlation and Regression Analysis
1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet
More informationFinite Difference Method for the Estimation of a Heat Source Dependent on Time Variable ABSTRACT
Malaysia Joural of Matematical Scieces 6(S): 39-5 () Special Editio of Iteratioal Worsop o Matematical Aalysis (IWOMA) Fiite Differece Metod for te Estimatio of a Heat Source Depedet o Time Variable, Allabere
More informationProblem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =
Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,
More informationVariable selection in principal components analysis of qualitative data using the accelerated ALS algorithm
Variable selectio i pricipal compoets aalysis of qualitative data usig the accelerated ALS algorithm Masahiro Kuroda Yuichi Mori Masaya Iizuka Michio Sakakihara (Okayama Uiversity of Sciece) (Okayama Uiversity
More informationTHEORETICAL RESEARCH REGARDING ANY STABILITY THEOREMS WITH APPLICATIONS. Marcel Migdalovici 1 and Daniela Baran 2
ICSV4 Cairs Australia 9- July, 007 THEORETICAL RESEARCH REGARDING ANY STABILITY THEOREMS WITH APPLICATIONS Marcel Migdalovici ad Daiela Bara Istitute of Solid Mechaics, INCAS Elie Carafoli, 5 C-ti Mille
More informationSection 11.8: Power Series
Sectio 11.8: Power Series 1. Power Series I this sectio, we cosider geeralizig the cocept of a series. Recall that a series is a ifiite sum of umbers a. We ca talk about whether or ot it coverges ad i
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationProperties and Hypothesis Testing
Chapter 3 Properties ad Hypothesis Testig 3.1 Types of data The regressio techiques developed i previous chapters ca be applied to three differet kids of data. 1. Cross-sectioal data. 2. Time series data.
More informationCEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering
CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationJacob Hays Amit Pillay James DeFelice 4.1, 4.2, 4.3
No-Parametric Techiques Jacob Hays Amit Pillay James DeFelice 4.1, 4.2, 4.3 Parametric vs. No-Parametric Parametric Based o Fuctios (e.g Normal Distributio) Uimodal Oly oe peak Ulikely real data cofies
More informationMeasurement uncertainty of the sound absorption
Measuremet ucertaity of the soud absorptio coefficiet Aa Izewska Buildig Research Istitute, Filtrowa Str., 00-6 Warsaw, Polad a.izewska@itb.pl 6887 The stadard ISO/IEC 705:005 o the competece of testig
More informationCorrelation Regression
Correlatio Regressio While correlatio methods measure the stregth of a liear relatioship betwee two variables, we might wish to go a little further: How much does oe variable chage for a give chage i aother
More informationMOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND.
XI-1 (1074) MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND. R. E. D. WOOLSEY AND H. S. SWANSON XI-2 (1075) STATISTICAL DECISION MAKING Advaced
More informationFinite Difference Approximation for First- Order Hyperbolic Partial Differential Equation Arising in Neuronal Variability with Shifts
Iteratioal Joural of Scietific Egieerig ad Research (IJSER) wwwiseri ISSN (Olie): 347-3878, Impact Factor (4): 35 Fiite Differece Approimatio for First- Order Hyperbolic Partial Differetial Equatio Arisig
More informationApproximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square
WSEAS TRANSACTONS o BUSNESS ad ECONOMCS S. Khotama, S. Boothiem, W. Klogdee Approimatig the rui probability of fiite-time surplus process with Adaptive Movig Total Epoetial Least Square S. KHOTAMA, S.
More informationTEACHER CERTIFICATION STUDY GUIDE
COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra
More informationPRELIM PROBLEM SOLUTIONS
PRELIM PROBLEM SOLUTIONS THE GRAD STUDENTS + KEN Cotets. Complex Aalysis Practice Problems 2. 2. Real Aalysis Practice Problems 2. 4 3. Algebra Practice Problems 2. 8. Complex Aalysis Practice Problems
More informationThe Phi Power Series
The Phi Power Series I did this work i about 0 years while poderig the relatioship betwee the golde mea ad the Madelbrot set. I have fially decided to make it available from my blog at http://semresearch.wordpress.com/.
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More informationNumerical Solution of the First-Order Hyperbolic Partial Differential Equation with Point-Wise Advance
Iteratioal oural of Sciece ad Research (ISR) ISSN (Olie): 39-74 Ide Copericus Value (3): 4 Impact Factor (3): 4438 Numerical Solutio of the First-Order Hyperbolic Partial Differetial Equatio with Poit-Wise
More informationQuantile regression with multilayer perceptrons.
Quatile regressio with multilayer perceptros. S.-F. Dimby ad J. Rykiewicz Uiversite Paris 1 - SAMM 90 Rue de Tolbiac, 75013 Paris - Frace Abstract. We cosider oliear quatile regressio ivolvig multilayer
More informationDiscrete Orthogonal Moment Features Using Chebyshev Polynomials
Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical
More informationThe Expectation-Maximization (EM) Algorithm
The Expectatio-Maximizatio (EM) Algorithm Readig Assigmets T. Mitchell, Machie Learig, McGraw-Hill, 997 (sectio 6.2, hard copy). S. Gog et al. Dyamic Visio: From Images to Face Recogitio, Imperial College
More informationNumerical Method for Blasius Equation on an infinite Interval
Numerical Method for Blasius Equatio o a ifiite Iterval Alexader I. Zadori Omsk departmet of Sobolev Mathematics Istitute of Siberia Brach of Russia Academy of Scieces, Russia zadori@iitam.omsk.et.ru 1
More informationA collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation
Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios
More informationStochastic Simulation
Stochastic Simulatio 1 Itroductio Readig Assigmet: Read Chapter 1 of text. We shall itroduce may of the key issues to be discussed i this course via a couple of model problems. Model Problem 1 (Jackso
More informationIntroduction to Machine Learning DIS10
CS 189 Fall 017 Itroductio to Machie Learig DIS10 1 Fu with Lagrage Multipliers (a) Miimize the fuctio such that f (x,y) = x + y x + y = 3. Solutio: The Lagragia is: L(x,y,λ) = x + y + λ(x + y 3) Takig
More information(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is
Calculus BC Fial Review Name: Revised 7 EXAM Date: Tuesday, May 9 Remiders:. Put ew batteries i your calculator. Make sure your calculator is i RADIAN mode.. Get a good ight s sleep. Eat breakfast. Brig:
More informationGeneralized Semi- Markov Processes (GSMP)
Geeralized Semi- Markov Processes (GSMP) Summary Some Defiitios Markov ad Semi-Markov Processes The Poisso Process Properties of the Poisso Process Iterarrival times Memoryless property ad the residual
More informationGoal. Adaptive Finite Element Methods for Non-Stationary Convection-Diffusion Problems. Outline. Differential Equation
Goal Adaptive Fiite Elemet Methods for No-Statioary Covectio-Diffusio Problems R. Verfürth Ruhr-Uiversität Bochum www.ruhr-ui-bochum.de/um1 Tübige / July 0th, 017 Preset space-time adaptive fiite elemet
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationTHE KALMAN FILTER RAUL ROJAS
THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationProbability, Expectation Value and Uncertainty
Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such
More informationA statistical method to determine sample size to estimate characteristic value of soil parameters
A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More information11 Correlation and Regression
11 Correlatio Regressio 11.1 Multivariate Data Ofte we look at data where several variables are recorded for the same idividuals or samplig uits. For example, at a coastal weather statio, we might record
More informationAlgebra of Least Squares
October 19, 2018 Algebra of Least Squares Geometry of Least Squares Recall that out data is like a table [Y X] where Y collects observatios o the depedet variable Y ad X collects observatios o the k-dimesioal
More informationComputational Fluid Dynamics. Lecture 3
Computatioal Fluid Dyamics Lecture 3 Discretizatio Cotiued. A fourth order approximatio to f x ca be foud usig Taylor Series. ( + ) + ( + ) + + ( ) + ( ) = a f x x b f x x c f x d f x x e f x x f x 0 0
More information577. Estimation of surface roughness using high frequency vibrations
577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia
More informationCS / MCS 401 Homework 3 grader solutions
CS / MCS 401 Homework 3 grader solutios assigmet due July 6, 016 writte by Jāis Lazovskis maximum poits: 33 Some questios from CLRS. Questios marked with a asterisk were ot graded. 1 Use the defiitio of
More informationNumerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets
Australia Joural of Basic ad Applied Scieces, 5(): 98-5, ISSN 99-878 Numerical Solutio of the Two Poit Boudary Value Problems By Usig Wavelet Bases of Hermite Cubic Splie Wavelets Mehdi Yousefi, Hesam-Aldie
More informationExpectation-Maximization Algorithm.
Expectatio-Maximizatio Algorithm. Petr Pošík Czech Techical Uiversity i Prague Faculty of Electrical Egieerig Dept. of Cyberetics MLE 2 Likelihood.........................................................................................................
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationFourier and Periodogram Transform Temperature Analysis in Soil
ISSN (e): 2250 3005 Volume, 05 Issue, 09 September 2015 Iteratioal Joural of Computatioal Egieerig Research (IJCER) Fourier ad Periodogram Trasform Aalysis i Soil Afolabi O.M. Adekule Ajasi Uiversity,
More informationTopic 5 [434 marks] (i) Find the range of values of n for which. (ii) Write down the value of x dx in terms of n, when it does exist.
Topic 5 [44 marks] 1a (i) Fid the rage of values of for which eists 1 Write dow the value of i terms of 1, whe it does eist Fid the solutio to the differetial equatio 1b give that y = 1 whe = π (cos si
More informationResearch Article A New Second-Order Iteration Method for Solving Nonlinear Equations
Abstract ad Applied Aalysis Volume 2013, Article ID 487062, 4 pages http://dx.doi.org/10.1155/2013/487062 Research Article A New Secod-Order Iteratio Method for Solvig Noliear Equatios Shi Mi Kag, 1 Arif
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationSeptember 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1
September 0 s (Edecel) Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright www.pgmaths.co.uk - For AS, A otes ad IGCSE / GCSE worksheets September 0 Copyright
More informationFinite Difference Approximation for Transport Equation with Shifts Arising in Neuronal Variability
Iteratioal Joural of Sciece ad Research (IJSR) ISSN (Olie): 39-764 Ide Copericus Value (3): 64 Impact Factor (3): 4438 Fiite Differece Approimatio for Trasport Equatio with Shifts Arisig i Neuroal Variability
More informationMathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution
America Joural of Theoretical ad Applied Statistics 05; 4(: 6-69 Published olie May 8, 05 (http://www.sciecepublishiggroup.com/j/ajtas doi: 0.648/j.ajtas.05040. ISSN: 6-8999 (Prit; ISSN: 6-9006 (Olie Mathematical
More informationIntroduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT
Itroductio to Extreme Value Theory Laures de Haa, ISM Japa, 202 Itroductio to Extreme Value Theory Laures de Haa Erasmus Uiversity Rotterdam, NL Uiversity of Lisbo, PT Itroductio to Extreme Value Theory
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More informationStatistical Pattern Recognition
Statistical Patter Recogitio Classificatio: No-Parametric Modelig Hamid R. Rabiee Jafar Muhammadi Sprig 2014 http://ce.sharif.edu/courses/92-93/2/ce725-2/ Ageda Parametric Modelig No-Parametric Modelig
More informationIntroduction to Optimization Techniques
Itroductio to Optimizatio Techiques Basic Cocepts of Aalysis - Real Aalysis, Fuctioal Aalysis 1 Basic Cocepts of Aalysis Liear Vector Spaces Defiitio: A vector space X is a set of elemets called vectors
More information(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?
MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle
More informationSNAP Centre Workshop. Basic Algebraic Manipulation
SNAP Cetre Workshop Basic Algebraic Maipulatio 8 Simplifyig Algebraic Expressios Whe a expressio is writte i the most compact maer possible, it is cosidered to be simplified. Not Simplified: x(x + 4x)
More informationTHE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
More informationExponential Families and Bayesian Inference
Computer Visio Expoetial Families ad Bayesia Iferece Lecture Expoetial Families A expoetial family of distributios is a d-parameter family f(x; havig the followig form: f(x; = h(xe g(t T (x B(, (. where
More informationTRACEABILITY SYSTEM OF ROCKWELL HARDNESS C SCALE IN JAPAN
HARDMEKO 004 Hardess Measuremets Theory ad Applicatio i Laboratories ad Idustries - November, 004, Washigto, D.C., USA TRACEABILITY SYSTEM OF ROCKWELL HARDNESS C SCALE IN JAPAN Koichiro HATTORI, Satoshi
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationADVANCED SOFTWARE ENGINEERING
ADVANCED SOFTWARE ENGINEERING COMP 3705 Exercise Usage-based Testig ad Reliability Versio 1.0-040406 Departmet of Computer Ssciece Sada Narayaappa, Aeliese Adrews Versio 1.1-050405 Departmet of Commuicatio
More information