Optimal Control of Temperature in Fluid Flow
|
|
- Grace Warren
- 5 years ago
- Views:
Transcription
1 Kawahara Lab. 5 March. 27 Optmal Control of Temperature n Flud Flow Dasuke YAMAZAKI Department of Cvl Engneerng, Chuo Unversty Kasuga -3-27, Bunkyou-ku, Tokyo 2-855, Japan E-mal : d33422@educ.kc.chuo-u.ac.jp Abstract Ths paper presents an optmal control problem of temperature usng the adjont equaton method of an optmal control theory and the fnte element method. The three dmensonal model s used to analyze the thermal flud flow, whch s descrbed by the ncompressble Naver-Stokes equaton and the energy equaton. The formulaton s based on the optmal control theory n whch the performance functon s expressed by the computed and target temperatures. The optmal value can be obtaned by mnmzng the performance functon. The fractonal step method s appled to solve the ncompressble Naver-Stokes equaton. The weghted gradent method s employed as a mnmzaton technque. Keyword : T hermal F lud F low, Optmal Control T heory, F ractonal Step M ethod, W eghted Gradent M ethod. Introducton The thermal flud flow s a phenomenon generated by non-homogenety of temperature. The thermal flud flow s dvded nto natural flud flow and forced flud flow. The Natural flud flow s caused by the dfference of the specfc gravty between warmed flud and cool flud. Forced flud flow s compulsorly caused by fan or other tems. In recent years, the optmal control of thermal flud flow s used n varous places n our lfe. For nstance, the optmal control of the thermal flud flow s appled to ar-condtonng, heatng, refrgerator and bath and so on. Therefore, the optmal control of thermal flud flow s ndspensable and very mportant problem. Especally, the optmal control n the natural flud flow s one of most-attracted recent problem. Ths s because the control brngs low cost and good for the ecology. Floor heatng system and temperature management of grass feld n the football stadum are good examples. Therefore, ths research covers the problem n the natural flud flow. The purpose of ths study s the control of temperature n the natural flud flow. The temperature s controlled usng the three dmensonal feld as a computatonal doman. The ncompressble Naver-Stokes and the energy equatons are employed for the state equaton. In the optmal control theory, the control varable that makes the optmal state can be obtaned by mnmzng the performance functon. The performance functon s composed of the square sum of dfference between computed and target temperatures. The control varable can be converged to the target temperature n case that the performance functon s mnmzed. The weghted gradent method s appled as the mnmzaton technque. The fractonal step method s appled to solve the Naver-Stokes equaton. The fractonal step method n ths research uses the ntermedate velocty that does not completely satsfy the contnuty equaton. The Crank-Ncolson method and the Galerkn method expanded by the mxed nterpolaton are appled to temporal and spatal dscretzaton, respectvely. The stablzed bubble functon s appled to velocty and temperature felds. The lnear nterpolaton s utlzed to pressure feld. Ths study presents an analyss and a control problem of thermal flud flow, n whch two numercal studes are carred out. In numercal study, thermal flud flow s analyzed as the three dmensonal model. In numercal study 2, the optmal control of thermal flud flow s performed out n the three dmensonal model.
2 2 State Equaton The thermal flud flow s descrbed by the ncompressble Naver-Stokes equaton and the energy equaton. Therefore, these equatons are employed for the state equaton. The ncompressble Naver-Stokes equaton s appled to the flud feld. The energy equaton s appled to the temperature feld. In the ncompressble flow, densty of flud s assumed constant. The Naver-Stokes equaton and energy equaton are wrtten as follows: where ν, κ and f are wrtten as follows: u + u j u,j + p, ν(u,jj + u j,j ) = f θ n Ω, () u, = n Ω, (2) θ + u j θ,j κθ,jj = n Ω, (3) ν = Pr, f = PrRa, (4) where u, p, θ and f are the velocty, the pressure, the temperature and the gravtatonal acceleraton, ν and κ are the knematc vscosty coeffcent and the dffuson coeffcent, respectvely. The Raylegh number and the Prandtl number are denoted by Ra and Pr, respectvely. The ntal condtons are gven as follows: u (t ) = û, (5) θ(t ) = ˆθ, (6) where û and ˆθ are the ntal known value for velocty and temperature, respectvely. The boundary condtons are gven as follows: u = û on Γ u D, (7) t = pδ j + ν(u,j + u j, )n j = ˆt on Γ u N, (8) θ = ˆθ on Γ θ D, (9) b = κθ, n = ˆb on Γ θ N, () θ = ˆθ cont on Γ θ C, () where subscrpts D, N and C of Γ are the Drchlet boundary, the Neumann boundary and the control boundary, respectvely. The control varable on the control boundary Γ θ C s denoted by ˆθ cont. 3 Fractonal Step Method To solve the state equaton and the adjont equaton, the Crank-Ncolson method s used for the temporal dscretzaton. Ths method s capable of takng long tme ncrement and superor n stablty. Therefore, a lot of tme cycles can be taken n the computaton. The momentum equaton, contnuty equaton and energy equaton are expressed as follows: u n+ u n t + u j u n+ 2,j + p n+, ν(u n+ 2,jj + u n+ 2 j,j ) = f θ n Ω, (2) u n+, = n Ω, (3) θ n+ θ n t + u j θ n+ 2,j κθ n+ 2,jj = n Ω, (4) u n+ 2 = 2 (un + un+ ), (5) θ n+ 2 = 2 (θn + θn+ ). (6) 2
3 The fractonal step method s appled to solve the ncompressble Naver-Stokes equaton. Ths method s a technque that s dvded nto the velocty and the pressure by dervng the pressure Posson equaton. The pressure Posson equaton s derved after obtanng the ntermedated velocty ũ. The pressure at the prevous step p n s assumed to be approxmate pressure. The pressure of the momentum equaton s replaced wth p n. Unknown velocty s replaced wth the ntermedate velocty. Then, the momentum equaton n Eq.(2) s as follows: ũ n+ u n t + u j ũ n+ 2,j + p n, ν(ũ n+ 2,jj + ũ n+ 2 j,j ) = f θ n Ω. (7) The dfference between Eq.(2) and Eq.(7) s taken, the followng equaton s obtaned. u n+ ũ n t + 2 u j(u n+,j ũ n+,j ) + (p n+, p n, ) 2 ν{(u,jj + u j,j ) n+ (ũ,jj + ũ j,j ) n+ } = f θ n Ω. (8) The pressure Posson equaton s obtaned from the transpraton of Eq.(8) consderng the contnuty equaton as follows: t(p n+, p n,) = ũ n+, n Ω. (9) 4 Spatal Dscretzaton 4. Bubble Functon Interpolaton The Galerkn method s appled to the spatal dscretzaton. The mxed nterpolaton based on the bubble functon and lnear functon s used for the spatal dscretzaton. The bubble functon nterpolaton s appled to the velocty and the temperature felds and the lnear nterpolaton s appled to the pressure feld as follows: For bubble functon nterpolaton: u = Φ u + Φ 2 u 2 + Φ 3 u 3 + Φ 4 u 4 + Φ 5 ũ 5, (2) ũ 5 = u 5 4 (u + u 2 + u 3 + u 4 ), (2) θ = Φ θ + Φ 2 θ 2 + Φ 3 θ 3 + Φ 4 θ 4 + Φ 5 θ5, (22) θ 5 = θ 5 4 (θ + θ 2 + θ 3 + θ 4 ), (23) and for lnear nterpolaton: Φ = L, Φ 2 = L 2, Φ 3 = L 3, Φ 4 = L 4, Φ 5 = 256L L 2 L 3 L 4, (24) p = Ψ p + Ψ 2 p 2 + Ψ 3 p 3 + Ψ 4 p 4, (25) Ψ = L, Ψ 2 = L 2, Ψ 3 = L 3, Ψ 4 = L 4, (26) 3
4 Fg. : Bubble functon element 2 Fg.2 : Lnear element 2 where u u 4 and θ θ 4 are the velocty and the temperature at nodes 4 of each fnte element. Φ α (α = 5) s the bubble functon for the velocty and the temperature, n whch p p 4 are the pressure at nodes 4 of each fnte element. The lnear nterpolaton functon s denoted by Ψ α (α = 4) for pressure. Uncoordnated s expressed by L. 4.2 Stablzed Form In the bubble functon, the numercal stablzaton s not enough. Therefore, the stablzed parameter whch determnes the magntude of streamlne s used for stablzaton. In equaton, the stablzed parameter τ eb can be wrtten as follows: τ eb = φ e, 2 Ω e ν φ e,j 2 Ω e A e, (27) where Ω e s element doman and < u, v >= uvdω, Ω e u 2 Ω e = uudω, Ω e Ae = dω. Ω e (28) The ntegral of bubble functon s expressed as follows: < φ e, > Ωe = Ae 6, u e,j 2 Ω e = 2A e g, g = 2 Ψ α, 2, (29) = where ν s control parameter for the stablzng acton. Ths value s determned to become equvalent to τ es usng the stablzed fnte element method. ( τ es = ) 2 2 u + h e ( ) 2 4ν h 2 e 2, u = u + u 2 + u 3, (3) where h e s an element sze. In generally, the stablzed parameter n Eq.(3) s not equal to the optmal parameter n Eq.(3). Thus, the bubble functon whch gves the optmal vscosty satsfes the followng equaton. φ e, 2 (ν + ν ) φ e,j 2 Ω e A e = τ es. (3) Eq.(32) adds stablzed operator control term only of the barycenter pont to the equaton of moton as follows: N e e= ν φ e,j 2 Ω e be, (32) where N e and be are the total number of element and the barycenter pont. In the energy equaton, ν s replaced wth κ. 4
5 5 Optmal Control Theory 5. Performance Functon In the optmal control theory, control varable can be obtaned by mnmzng performance functon. The performance functon J s composed of the square sum of dfference between computed and target temperatures. The control varable are computed to mnmze the performance functon under the state equaton and boundary condtons. Mnmzng the performance functon means that computed temperature should be as close as possble to the objectve temperature. The performance functon s wrtten as follows: J = 2 tf t Ω (θ θ obj ) T Q(θ θ obj )dωdt, (33) where θ, θ obj and Q are the computed temperature, the target temperature at the object pont and the weghtng dagonal matrx. Superscrpt T means transpose. 5.2 Adjont Equaton The Lagrange multpler method s appled to mnmze the performance functon. The Lagrange multpler method s used for the mnmzaton problem wth constrants. The extended performance functon J s expressed as follows: tf J = J u T ( u + u j u,j + p, ν(u,jj + u j,j ) f θ)dωdt t Ω tf p T (u, )dωdt t tf t Ω Ω θ T ( θ + u j θ,j κθ,jj )dωdt, (34) where u, p and θ denote the Lagrange multpler for velocty, pressure and temperature, respectvely. The extended performance functon J s dvded nto the Hamltonan H and the tme dfferentaton term as follows: J = tf where the Hamltonan H s defned as follows: H = 2 (θ θ obj) T Q(θ θ obj ) t Ω {H + (u T u + θ T θ)}dωdt, (35) u T {u j u,j + p, ν(u,jj + u j,j ) f θ} p T u, θ T (u j θ,j κθ,jj ). (36) The statonary condton s needed to mnmze performance functon. The frst varaton of J should be zero as follows, δj =. (37) Therefore, the adjont equaton and termnal condtons are obtaned. The adjont equaton and the termnal condtons are wrtten as follows: u T = H u, H p =, θ T = H θ, (38) u (t f ) =, p (t f ) =, θ (t f ) =. (39) 5
6 6 Mnmzaton Technque 6. Weghted Gradent Method The weghted gradent method s appled to the mnmzaton technque. Ths method s a technque that control varable are renewed by the statonary condton of a modfed performance functon. The modfed performance functon K (l) to whch the penalty term s added s defned as follows: K (l) = J + 2 ( θ (l+) (t) θ (l) (t)) T W (l) ( θ (l+) (t) θ (l) (t)), (4) where θ(t), W (l) and l are control varable, weghtng dagonal matrx and teraton number. To mnmze the modfed performance functon, the penalty term should be zero. To mnmze the extended performance functon s equal to mnmze the modfed performance functon. Then, the followng equaton s obtaned. K θ The control varable s renewed by the followng equaton: =. (4) θ (l+) (t) = θ (l) (t) W (l) grad(j)(l). (42) 7 Numercal Study 7. Numercal Study The numercal analyss of natural flud flow s carred out usng the three dmensonal model. The fnte element mesh and computatonal model are shown n Fgs.3 and 4. The total number of nodes and elements are 46 and 73728, respectvely. On the bottom surface, temperature at ponts s specfed to be equal to.. On the top surface, temperature s set to be equal to.. The boundary condton for velocty s non-slp condton on all the surface. The Raylegh number and the Prandtl number are set to. and.7, respectvely. The tme ncrement t s.3. The natural flud flow usng the fnte element method s carred out. The temperature of non-dmensonal tme t=.45 s shown n Fg.7. The velocty of non-dmensonal tme t=.45 s llustrated n Fg.8. It s seen that phenomenon of thermal flud flow s smulated. 3 T=. T=. Z X Y Fg.3 : Fnte Element Mesh Fg.4 : Computatonal Model 6
7 7.2 Numercal Study 2 The optmal control of temperature n natural flud flow s performed. The fnte element mesh and computaton model are shown n Fgs.5 and 6. The total number of nodes and elements are 2975 and 356, respectvely. The control and the specfed boundares are shown n Fg.6. The temperature on the specfed boundary s shown n Fg.9. The boundary condton for velocty s gven as non-slp condton on the all surface. In ths study, the Raylegh number and the Prandtl number are set to. and.7, respectvely. The tme ncrement t s.. The object pont s set as shown n Fg.6. Usng the specfed temperature shown n Fg.9, the natural flud flow s occurred. The temperature at the object pont shown n Fg. can be obtaned. To reduce the temperature at the object pont from ntal to target temperature, the optmal control s carred out. The target temperature s set to.. Ths optmal control problem s to fnd the control temperature on the control boundares so as to mnmze the performance functon. The weghted gradent method s appled as the mnmzaton technque. The tme hstory of the performance functon s shown n Fg.. The tme hstory of the control temperature at the control boundary s shown n Fg.2. The tme hstory of temperature at the object pont s shown n Fg.3. At the target pont, the temperature could be controlled to the objectve temperature. The temperatures on Y=.5 and Z=.25 are shown n Fg.4. The temperature and the velocty wthout control at non-dmensonal tme t=2. are shown n Fgs.5 and 7. The temperature and the velocty wth control at non-dmensonal tme t=2. are shown n Fgs.6 and 8. The temperature at the same alttude can be reduced by settng only one object pont as shown n shown Fg.4. The optmal control of temperature n natural flud flow can be carred out..5 Z X Y Control Boundary Z Specfed Boundary Object Pont X Y Fg.5 : Fnte Element Mesh Fg.6 : Computatonal Model 8 Concluson The optmal control of temperature n the natural flud flow s presented. The control temperature can be found so as to mnmze the performance functon usng the weghted gradent method. It s shown that the temperature at object pont can be reduced usng the control temperature. The all temperature of the same alttude s reduced by settng only one object pont. The temperature at the object pont s well controlled usng the present method. References. M.Kawahara, K.Sasak and Y.Sano, Parameter Identfcaton and Optmal Control of Ground Temperature, Int.J.N umer.m eth.f l., vol.2,pp.789-8, S.Suzuk, A.Anju and M.Kawahara, Management of Ground Temperature by Bang-Bang Control Based on Fnte Element Applcaton, Int.J.N umer.m esh.eng., vol.39, pp.885-9, J.Matsumoto, A.A.Khan, S.S. Wang, and M.Kawahara, Shallow Water Flow Analyss wth Movng Boundary Technque Usng Least-Squares Bubble Functon, Int.J.Comp.FludDyn., vol.6(2), pp.29-34, J.Matsumoto, T.Umetsu and M.Kawahara, Stablzed Bubble Functon Method for Shallow Water Long Wave Equaton,Int.J.Comp.F luddyn., vol.7(4), pp , T. E. Tezduyar, Stablzed Fnte Element Formulatons for Incompressble Flow Computatons, AdvancesnAppl.M ech., vol.28, pp.-42,992 7
8 Fg.7 : Temperature ( t=.45 ) Fg.8 : Velocty ( t=.45 )..9 Specfed Temperature.2. Wthout Control Temperature Temperature Non-Dmensonal Tme Fg.9 : Specfed Temperature Non-Dmensonal Tme Fg. : Temperature at Object Pont wthout control Performance Functon. -. Control Temperature Performance Functon 2.5 Temperature Itaraton Number Non-Dmensonal Tme Fg.: Performance Functon Fg.2 : Control Temperature at Control boundary 8
9 .2. Wth Control Wthout Control.5.4 Wth Control Wthout Control.8.3 Temperature.6.4 Temperature Non-Dmensonal Tme Fg.3 : Temperature at Object Pont X-axs Fg.4 : Temperature Dstrbuton (t=2.) Fg.5 : Temperature Wthout Control (t=.6) Fg.6 : Temperature Wth Control (t=.6) 9
10 Fg.7 : Velocty Wthout Control (t=.6) Fg.8 : Velocty Wth Control (t=.6)
Optimal Control of Dissolved Oxygen in Shallow Water Flow
Kawahara Lab. 21 Mar. 2009 Optmal Control of Dssolved Oxygen n Shallow Water Flow Toshk SEKINE Department of Cvl Engneerng, Chuo Unversty, Kasuga 1-13-27, Bunkyo-ku, Tokyo 112-8551, Japan E-mal : d34115@educ.kc.chuo-u.ac.jp
More informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationSolution of the Navier-Stokes Equations
Numercal Flud Mechancs Fall 2011 Lecture 25 REVIEW Lecture 24: Soluton of the Naver-Stokes Equatons Dscretzaton of the convectve and vscous terms Dscretzaton of the pressure term Conservaton prncples Momentum
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationHigher Order Wall Boundary Conditions for Incompressible Flow Simulations
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN KOREA OCTOBER 7-30 003 Hgher Order Wall Boundary Condtons for Incompressble Flow Smulatons Hdetosh Nshda. Department of Mechancal and System Engneerng
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationThe interface control domain decomposition (ICDD) method for the Stokes problem. (Received: 15 July Accepted: 13 September 2013)
Journal of Coupled Systems Multscale Dynamcs Copyrght 2013 by Amercan Scentfc Publshers All rghts reserved. Prnted n the Unted States of Amerca do:10.1166/jcsmd.2013.1026 J. Coupled Syst. Multscale Dyn.
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationThe Finite Element Method: A Short Introduction
Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationSolutions to exam in SF1811 Optimization, Jan 14, 2015
Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationSIO 224. m(r) =(ρ(r),k s (r),µ(r))
SIO 224 1. A bref look at resoluton analyss Here s some background for the Masters and Gubbns resoluton paper. Global Earth models are usually found teratvely by assumng a startng model and fndng small
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationArmy Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More informationVisco-Rubber Elastic Model for Pressure Sensitive Adhesive
Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA
More informationAdjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R.
th World Congress on Structural and Multdscplnary Optmsaton 7 th - th, June 5, Sydney Australa Adjont Methods of Senstvty Analyss for Lyapunov Equaton Bopng Wang, Kun Yan Department of Mechancal and Aerospace
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationRelaxation Methods for Iterative Solution to Linear Systems of Equations
Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary
More informationMATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)
1/16 MATH 829: Introducton to Data Mnng and Analyss The EM algorthm (part 2) Domnque Gullot Departments of Mathematcal Scences Unversty of Delaware Aprl 20, 2016 Recall 2/16 We are gven ndependent observatons
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationNumerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline
IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP 6-30 www.osrournals.org Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal
More informationThe Minimum Universal Cost Flow in an Infeasible Flow Network
Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran
More informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationNovember 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project
SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω
More informationUniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity
Int. Journal of Math. Analyss, Vol. 6, 212, no. 22, 195-114 Unqueness of Weak Solutons to the 3D Gnzburg- Landau Model for Superconductvty Jshan Fan Department of Appled Mathematcs Nanjng Forestry Unversty
More informationA Cartesian-grid integrated-rbf method for viscoelastic flows
Home Search Collectons Journals About Contact us My IOPscence A Cartesan-grd ntegrated-rbf method for vscoelastc flows Ths artcle has been downloaded from IOPscence. Please scroll down to see the full
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationA new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems
Home Search Collectons Journals About Contact us My IOPscence A new ntegrated-rbf-based doman-embeddng scheme for solvng flud-flow problems Ths artcle has been downloaded from IOPscence. Please scroll
More informationarxiv: v1 [physics.flu-dyn] 16 Sep 2013
Three-Dmensonal Smoothed Partcle Hydrodynamcs Method for Smulatng Free Surface Flows Rzal Dw Prayogo a,b, Chrstan Fredy Naa a a Faculty of Mathematcs and Natural Scences, Insttut Teknolog Bandung, Jl.
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationExtension of Smoothed Particle Hydrodynamics (SPH), Mathematical Background of Vortex Blob Method (VBM) and Moving Particle Semi-Implicit (MPS)
Amercan Journal of Computatonal athematcs, 04, 5, 44-445 Publshed Onlne December 04 n ScRes. http://www.scrp.org/ournal/acm http://dx.do.org/0.436/acm.04.45036 Extenson of Smoothed Partcle Hydrodynamcs
More informationNumerical Solution of two dimensional coupled viscous Burgers Equation using the Modified Cubic B-Spline Differential Quadrature Method
umercal Soluton of two dmensonal coupled vscous Burgers Equaton usng the odfed Cubc B-Splne Dfferental Quadrature ethod H. S. Shukla 1, ohammad Tamsr 1*, Vneet K. Srvastava, Ja Kumar 3 1 Department of
More informationSolutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution.
Solutons HW #2 Dual of general LP. Fnd the dual functon of the LP mnmze subject to c T x Gx h Ax = b. Gve the dual problem, and make the mplct equalty constrants explct. Soluton. 1. The Lagrangan s L(x,
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationNUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD
THERMAL SCIENCE: Year 2018, Vol. 22, No. 5, pp. 1955-1962 1955 NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD Introducton by Tomok IZUMI a* and
More informationMultigrid Methods and Applications in CFD
Multgrd Metods and Applcatons n CFD Mcael Wurst 0 May 009 Contents Introducton Typcal desgn of CFD solvers 3 Basc metods and ter propertes for solvng lnear systems of equatons 4 Geometrc Multgrd 3 5 Algebrac
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationThree-dimensional eddy current analysis by the boundary element method using vector potential
Physcs Electrcty & Magnetsm felds Okayama Unversty Year 1990 Three-dmensonal eddy current analyss by the boundary element method usng vector potental H. Tsubo M. Tanaka Okayama Unversty Okayama Unversty
More informationElectrostatic Potential from Transmembrane Currents
Electrostatc Potental from Transmembrane Currents Let s assume that the current densty j(r, t) s ohmc;.e., lnearly proportonal to the electrc feld E(r, t): j = σ c (r)e (1) wth conductvty σ c = σ c (r).
More informationVARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES
VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES BÂRZĂ, Slvu Faculty of Mathematcs-Informatcs Spru Haret Unversty barza_slvu@yahoo.com Abstract Ths paper wants to contnue
More informationFlow Induced Vibration
Flow Induced Vbraton Project Progress Report Date: 16 th November, 2005 Submtted by Subhrajt Bhattacharya Roll no.: 02ME101 Done under the gudance of Prof. Anrvan Dasgupta Department of Mechancal Engneerng,
More informationA large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method
A large scale tsunam run-up smulaton and numercal evaluaton of flud force durng tsunam by usng a partcle method *Mtsuteru Asa 1), Shoch Tanabe 2) and Masaharu Isshk 3) 1), 2) Department of Cvl Engneerng,
More information1 Derivation of Point-to-Plane Minimization
1 Dervaton of Pont-to-Plane Mnmzaton Consder the Chen-Medon (pont-to-plane) framework for ICP. Assume we have a collecton of ponts (p, q ) wth normals n. We want to determne the optmal rotaton and translaton
More informationPublication 2006/01. Transport Equations in Incompressible. Lars Davidson
Publcaton 2006/01 Transport Equatons n Incompressble URANS and LES Lars Davdson Dvson of Flud Dynamcs Department of Appled Mechancs Chalmers Unversty of Technology Göteborg, Sweden, May 2006 Transport
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationHybridizable discontinuous Galerkin (HDG) method for Oseen flow
Hybrdzable dscontnuous Galerkn (HDG) method for Oseen flow Matteo Gacomn Laborator de Càlcul Numèrc, E.T.S. de Ingeneros de Camnos, Unverstat Poltècnca de Catalunya BarcelonaTech, Jord Grona 1, E-08034
More informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationIntroduction to Computational Fluid Dynamics
Introducton to Computatonal Flud Dynamcs M. Zanub 1, T. Mahalakshm 2 1 (PG MATHS), Department of Mathematcs, St. Josephs College of Arts and Scence for Women-Hosur, Peryar Unversty 2 Assstance professor,
More information6.3.7 Example with Runga Kutta 4 th order method
6.3.7 Example wth Runga Kutta 4 th order method Agan, as an example, 3 machne, 9 bus system shown n Fg. 6.4 s agan consdered. Intally, the dampng of the generators are neglected (.e. d = 0 for = 1, 2,
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationSimulation of Turbulent Flow Using FEM
Internatonal Journal of Engneerng and Technology Volume 2 No. 8, August, 2012 Smulaton of Turbulent Flow Usng FEM Sabah Tamm College of Computng, AlGhurar Unversty, Duba, Unted Arab Emrates. ABSTRACT An
More informationMEM Chapter 4b. LMI Lab commands
1 MEM8-7 Chapter 4b LMI Lab commands setlms lmvar lmterm getlms lmedt lmnbr matnbr lmnfo feasp dec2mat evallm showlm setmvar mat2dec mncx dellm delmvar gevp 2 Intalzng the LMI System he descrpton of an
More informationVector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.
Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm
More informationSolutions for Euler and Navier-Stokes Equations in Powers of Time
Solutons for Euler and Naver-Stokes Equatons n Powers of Tme Valdr Montero dos Santos Godo valdr.msgodo@gmal.com Abstract We present a soluton for the Euler and Naver-Stokes equatons for ncompressble case
More informationAcceleration of Gas Bubble-Free Surface Interaction Computation Using Basis Preconditioners
Acceleraton of Gas Bubble-Free Surface Interacton Computaton Usng Bass Precondtoners K. L. Tan, B. C. Khoo and J. K. Whte ABSTRACT The computaton of gas bubble-free surface nteracton entals a tme-steppng
More informationJournal of Fluid Science and Technology
Journal of Flud Scence and Technology Numercal Smulaton of Incompressble Flows around a Fsh Model at Low Reynolds Number Usng Seamless Vrtual Boundary Method * Hdetosh NISHIDA ** and Kyohe TAJIRI ** **Department
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationAnalytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures
IMS 2 Workshop Analytcal Gradent Evaluaton of Cost Functons n General Feld Solvers: A Novel Approach for Optmzaton of Mcrowave Structures P. Harscher, S. Amar* and R. Vahldeck and J. Bornemann* Swss Federal
More informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More informationA numerical simulation on seepage-induced collapse of caisson breakwater using a stabilized ISPH method
A numercal smulaton on seepage-nduced collapse of casson breakater usng a stablzed ISPH method *Tetsuro Goda 1) and Mtsuteru Asa ) 1), ) Department of Cvl and Structural Engneerng, Kyushu Unversty, 744
More informationHomogeneous model: Horizontal pipe and horizontal well. Flow loops can't duplicate field conditions. Daniel D. Joseph. April 2001
Homogeneous model of producton of heavy ol through horzontal ppelnes and wells based on the Naver-Stokes equatons n the ppelne or the well and Darcy's law n the reservor Homogeneous model: Danel D. Joseph
More informationLagrangian Field Theory
Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationLecture 12. Modeling of Turbulent Combustion
Lecture 12. Modelng of Turbulent Combuston X.S. Ba Modelng of TC Content drect numercal smulaton (DNS) Statstcal approach (RANS) Modelng of turbulent non-premxed flames Modelng of turbulent premxed flames
More informationComputational Electromagnetics in Antenna Analysis and Design
Computatonal Electromagnetcs n Antenna Analyss and Desgn Introducton It s rare for real-lfe EM problems to fall neatly nto a class that can be solved by the analytcal methods presented n the precedng lectures.
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationUncertainty as the Overlap of Alternate Conditional Distributions
Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant
More information