Numerical Heat and Mass Transfer


 Shannon Riley
 4 years ago
 Views:
Transcription
1 Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06FnteDfference Method (Onedmensonal, steady state heat conducton) Fausto Arpno
2 Introducton Why we use models and how Many physcal systems of nterest are extremely complex, wth many nteractng elements governng ther behavor. The way the complete system behaves may be explored by experments. However, f t s a nonlnear system the effect of two changes made together s very dfferent from the sum of the changes made separately so that t may be mpossble to derve enough nformaton from emprcal parameter analyss. Alternatvely, buldng or runnng experments may be expensve or dangerous. Numercal models ncorporate the knowledge ganed by experments allowng the runnng of new tests by computer smulatons.
3 Introducton Man steps of numercal smulaton Physcs approxmaton ü Is the mathematcal model avalable and able to accurately descrbe the physcs of the system under nvestgaton? ü In buldng computatonal models the skll requred s to smulate as well as possble a physcal system n order to nvestgate ts behavor. ü The term as well as possble s subject to a whole range of condtons (on one hand the physcs of the system mght be qute well known but the system may be so large that t s not possble to smulate t perfectly; on the other hand, the physcs of the system s not so well known and we known just general behavor). Numercal approxmaton ü Preprocessng: conssts of the nput of a problem to CFD program by means of an userfrendly nterface and the subsequent transformaton of the nput nto a form sutable for the solver. ü Processng: conssts of the numercal resoluton of governng Partal Dfferental Equatons (fnte dfference, fnte element, fnte volume, spectral methods, etc. ü Postprocessng: conssts of the results analyss (verfcaton, valdaton, ) 3
4 Mathematcal model The quanttes of nterest are descrbed by Partal Dfferental Equatons (PDEs) n the form: ( ) F x,θ,ϕ, ϕ L ϕ,g θ, ϕ x,..., ϕ x,...,g = 0 Smlar PDEs can be used for the resoluton of dfferent problems. The governng PDEs equatons represent the mathematcal statements of the conservaton law of physcs. z EQUATIONS EXAMPLES... ü Transport equaton ü Heat equaton ü Waves equaton ϕ θ + uϕ ( ) = 0 ϕ θ ϕ = G ϕ θ + ϕ = 0 Ω Ω. q x. q y x. q z+dz. q z y. q y+dy. ρc q x+ u. T θ 4
5 Numercal modelng Snce the exact soluton φ(x,y,z,ϑ) can be very dffcult to obtan, thought numercal models t s possble to obtan an approxmated soluton φ N (x,y,z,ϑ ), wth =1,,N, wth the use of a computer L(ϕ,g) = 0 Exact operator ϕ(x, ϑ) Numercal methods L N (ϕ N,g N ) = 0 Approxmated operator ϕ = 1,...N L N ( ϕ,g ) = L( ϕ,g ) + E T Roundoff error A numercal method s convergent f: A numercal method s consstent f: lm N lm L N N ϕ ϕ N = 0 ( ϕ,g ) = 0 A numercal method s stable small roundoff errors cause small (and convergent) soluton oscllatons. 5
6 Numercal modelng Lax s equvalence theorem Consstency + Stablty Convergence Convergence s usually very dffcult to establsh theoretcally and n practce we use the Lax s equvalence theorem whch states that for lnear problems a necessary and suffcent condton for convergence s that the method s both consstent and stable. In CFD methods such theorem s of lmted use snce the governng equatons are nonlnear. In such problems, consstency and stablty are necessary condtons of convergence, but not suffcent. 6
7 Numercal modelng The objectve s then the constructon of an equvalent, approxmated operator, usually as a system of lnear algebrac equatons :... θ = 0 How s t possble to acheve ths goal? a 11 ϕ 1 + a 1 ϕ... + a 1n ϕ n = 0 a n1 ϕ 1 + a n ϕ... + a nn ϕ n = 0 Fnte dfference: Fntedfference methods approxmate the solutons to dfferental equatons by approxmatng the dervatve expressons n the governng equaton; the method s easy to understand, but dffcult to mplement. Fnte element: Based on the weak formulaton and on the nterpolaton, the fnte element method s less ntutve, but powerful, sutable for multphyscs and smple to mplement. Fnte volume: The Fnte Volume method s a refned verson of the fnte dfference method and has became popular n CFD. 7
8 Fntedfference method for solvng heat conducton problems The numercal method of soluton s used n practcal applcatons to determne the temperature dstrbuton and heat flow n solds havng complcated geometres, boundary condtons, and temperaturedependent propertes. A commonly used numercal scheme (especally n the past) s the fntedfference method. In ths method, the partal dfferental equaton of heat conducton s approxmated by a set of algebrac equatons for temperature at a number of nodal ponts over the regon. Snce the method transforms the analyss of heat conducton problem to the soluton of a set of coupled algebrac equatons, t s also mportant to manage the methods of solvng smultaneous algebrac equatons. When a heat conducton problem s solved exactly by an analytcal method, the resultng soluton satsfes the governng dfferental equaton at every doman pont. When the problem s solved by a numercal method, such as fntedfference, the dfferental equaton s satsfed only n correspondence of dscrete number of ponts, called nodes. The fntedfference method can be developed by replacng the partal dervatves n the heat conducton equaton wth ther equvalent fntedfference forms or wrtng an energy balance for a dfferental volume element. 8
9 1D Fntedfference method: Mathematcal formulaton Consder the followng onedmensonal, steady state heat conducton equaton wthout energy generaton dt = 0 The computatonal doman s dscretzed n space usng a unform grd (or mesh). The nodal temperature T s defned n each node of the computatonal grd. T L M x It s assumed that the temperature n the doman s contnuous, dervable wth contnuous and lmted dervatve. Ths assumpton s necessary to fnd the approxmate expresson of the temperatures dervatves, wth the ad of a Taylor seres expanson. 9
10 1D Fntedfference method: Mathematcal formulaton T L M x Wth reference to the above dscretzed doman, the temperature T +1 at the node +1 can be expressed as a functon of the temperature T at the node by usng a Taylor seres expanson. dt d T T = T + Δ x+ + o + 1 ( ) 3 (1) Neglectng the second order terms, t s possble to obtan an approxmated expresson for the dervatve at the node (forward approxmaton): dt T T ( x) + 1 = + o Δ CUT OFF ERROR 10
11 1D Fntedfference method: Mathematcal formulaton T L M x Wth reference to the above dscretzed doman, the temperature T 1 at the node 1 can be expressed as a functon of the temperature T at the node by usng a Taylor seres expanson. dt d T T = T Δ x+ + o 1 ( ) 3 () Neglectng the second order terms, t s possble to obtan an approxmated expresson for the dervatve at the node (backward approxmaton): dt T T ( x) 1 = + o Δ CUT OFF ERROR 11
12 1D Fntedfference method: Mathematcal formulaton T L M x Both the dervatve approxmatons (forward and backward) mply a leadng error of order. To reduce the cut off error t s necessary to reduce the elements dmenson (). It can be observed that subtractng eq. () from eq. (1) T T = T T dt d T + Δ x + dt + o ( ) 3 t s possble to obtan an expresson of the temperature dervatve wth a second order approxmaton level (central approxmaton). dt T+ 1 T 1 = + o ( ) CUT OFF ERROR 1
13 1D Fntedfference method: Mathematcal formulaton Graphc representaton of the dervatve forward approxmaton (frst order) T(x) T +1 dt d T T = T + Δ x+ + o + 1 ( ) 3 T dt α x x +1 dt T T x ( x) + 1 = + o Δ 13
14 1D Fntedfference method: Mathematcal formulaton Graphc representaton of the dervatve central approxmaton (second order) T(x) T +1 T T 1 dt T T ( x) = + o Δ α x 1 x x +1 x 14
15 1D Fntedfference method: Mathematcal formulaton In a smlar way, the second temperature dervatve can be approxmated by addng the members of the followng equatons: 3 3 dt d T d T T = T + Δ x+ + + o dt T+ 1 T + T 1 The fnal equaton represents the approxmated expresson of the second dervatve (central second dervatve). In the case of nternal nodes, the above analyss leads: ( x) 3 ( ) 3 3 dt d T d T T = T Δ x+ + o 6 + = + o Δ ( ) 4 4 T T + T =
16 FnteDfference from energy balance Consder, agan, a onedmensonal, steady state heat conducton wthout energy generaton n a fnte regon 0 x L. We dvde the regon un a number of subregons (volumes). The center of the volumes corresponds to the nodes n the mathematcal fntedfference formulaton; each node presents a temperature value whch s representatve of the whole subregon temperature. To develop the fntedfference equaton, we consder a dfferental volume element about the node. The steady state energy balance equaton for ths volume element can be stated as:! Q 1 Q!! Q Q netta = Q 1 Q +1 Q = 0 ka ka ( T T ) + ( T T) hp( T T ) = T+ 1 T + T+ 1 hp = Δ x T T = m Δ x T T ka ( ) ( ) 16
17 FnteDfference: Boundary Condtons The boundary condtons for heat conducton problem may be a prescrbed temperature, prescrbed heat flux or convecton boundary condton. Prescrbed temperature The temperature T 0 and T M at nodes x=0 and x=l are known and ths provdes the two addtonal relatons needed to make the number of equatons equal the number of unknown nodal temperatures Prescrbed heat flux (adabatc end) Supposng that the heat flux s prescrbed at the boundary x=l, to develop the fntedfference form of ths boundary condton, we need to wrte the energy balance equaton for a dfferental volume Dx/ at node M. ka M M 1 = M ( T T ) hp ( T T ) T M T M 1 Q M 1 M = Q M = m ( T M T ) Q! M 1 M M1 M / x Q! = 0 17
18 Example: steady state heat conducton n a fn wth adabatc tp T s T (surroundng flud temperature) h=h c +h =cost A (Adabatc fn tp) L A=cross secton area P=cross secton permeter d T = m T T T ( ) 0 Ts ka ( ) where m = hp dt = = 0 L 18
The 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 informationOnesided finitedifference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note Onesded fntedfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
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 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 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 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 informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationNONCENTRAL 7POINT 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 NONCENTRAL 7POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
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 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 WomenHosur, Peryar Unversty 2 Assstance professor,
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscretetme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
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 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 informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng  ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng  ABCM, Curtba, Brazl, Dec. 58, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECONDORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECONDORDER HYPERBOLIC EUATION
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 informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
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 informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 19329466 Vol. 10, Issue 1 (June 2015), pp. 223229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationTHE STURMLIOUVILLE EIGENVALUE PROBLEM  A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 736 www.amcm.pcz.pl piss 999965 DOI: 0.75/jamcm.06..4 eiss 3530588 THE STURMLIOUVILLE EIGEVALUE PROBLEM  A UMERICAL SOLUTIO USIG THE
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
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 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 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 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 informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons Â and ˆb avalable. Then the best thng we can do s to solve Âˆx ˆb exactly whch
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 informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
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 information3.1 Expectation of Functions of Several Random Variables. )' be a kdimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
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 informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the RungeKutta method one employs an teratve mplct technque. Ths s because the structure
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTERVI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 008 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
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 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 informationSolving Nonlinear Differential Equations by a Neural Network Method
Solvng Nonlnear Dfferental Equatons by a Neural Network Method Luce P. Aarts and Peter Van der Veer Delft Unversty of Technology, Faculty of Cvlengneerng and Geoscences, Secton of Cvlengneerng Informatcs,
More informationApplication of BSpline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of BSplne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationThermalFluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)
hermalfluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 4106323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In
More informationThe Twoscale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites
7 AsaPacfc Engneerng Technology Conference (APETC 7) ISBN: 9786595443 The Twoscale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fntedmensonal 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 informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationON MECHANICS WITH VARIABLE NONCOMMUTATIVITY
ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY CIPRIAN ACATRINEI Natonal Insttute of Nuclear Physcs and Engneerng P.O. Box MG6, 07725Bucharest, Romana Emal: acatrne@theory.npne.ro. Receved March 6, 2008
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 informationProcedia Computer Science
Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5*  NOTE (Summary) ECON 5*  NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v101570100y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationThe Discretization Process
FMIA F Moukalled L Mangan M Darwsh An Advanced Introducton wth OpenFOAM and Matlab Ths textbook explores both the theoretcal foundaton of the Fnte Volume Method (FVM) and ts applcatons n Computatonal Flud
More informationA Cartesiangrid integratedrbf method for viscoelastic flows
Home Search Collectons Journals About Contact us My IOPscence A Cartesangrd ntegratedrbf method for vscoelastc flows Ths artcle has been downloaded from IOPscence. Please scroll down to see the full
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 informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
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 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 informationIrregular vibrations in multimass discretecontinuous systems torsionally deformed
(2) 4 48 Irregular vbratons n multmass dscretecontnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscretecontnuous systems consstng of an arbtrary number rgd bodes connected
More informationNumerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method
Appled Mathematcs, 6, 7, 54 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationAERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTINGLINE THEORY
LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTINGLINE THEORY The BotSavart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot Savart formula
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 informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
More informationMATH 5630: Discrete TimeSpace Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete TmeSpace Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 50757 HIKARI Ltd, wwwmhkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
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 informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 97457X (Onlne), ISSN 9745718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 1314 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationAsymptotics of the Solution of a Boundary Value. Problem for OneCharacteristic Differential. Equation Degenerating into a Parabolic Equation
Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5  HIKARI Ltd, www.mhar.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for OneCharacterstc Dfferental Equaton
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 ONEDIMENSIONAL, ONEPHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the faketest data; fxed
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationIntegrals and Invariants of EulerLagrange Equations
Lecture 16 Integrals and Invarants of EulerLagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
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 informationME 501A Seminar in Engineering Analysis Page 1
umercal Solutons of oundaryvalue Problems n Os ovember 7, 7 umercal Solutons of oundary Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work1
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 informationOn a direct solver for linear least squares problems
ISSN 20666594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationSuppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl
RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for ConsumptionSaving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for ConsumptonSavng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
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 Durgapur713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationDiscretization. Consistency. Exact Solution Convergence x, t > 0. Figure 5.1: Relation between consistency, stability, and convergence.
Chapter 5 Theory The numercal smulaton of PDE s requres careful consderatons of propertes of the approxmate soluton. A necessary condton for the scheme used to model a physcal problem s the consstency
More informationA large scale tsunami runup simulation and numerical evaluation of fluid force during tsunami by using a particle method
A large scale tsunam runup 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 informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LSDYNA
14 th Internatonal Users Conference Sesson: ALEFSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53535, 0 ISSN: 0407467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More informationComputational Electromagnetics in Antenna Analysis and Design
Computatonal Electromagnetcs n Antenna Analyss and Desgn Introducton It s rare for reallfe EM problems to fall neatly nto a class that can be solved by the analytcal methods presented n the precedng lectures.
More informationarxiv: v1 [math.ho] 18 May 2008
Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc
More informationStat260: Bayesian Modeling and Inference Lecture Date: February 22, Reference Priors
Stat60: Bayesan Modelng and Inference Lecture Date: February, 00 Reference Prors Lecturer: Mchael I. Jordan Scrbe: Steven Troxler and Wayne Lee In ths lecture, we assume that θ R; n hgherdmensons, reference
More informationGroup Analysis of Ordinary Differential Equations of the Order n>2
Symmetry n Nonlnear Mathematcal Physcs 997, V., 64 7. Group Analyss of Ordnary Dfferental Equatons of the Order n> L.M. BERKOVICH and S.Y. POPOV Samara State Unversty, 4430, Samara, Russa Emal: berk@nfo.ssu.samara.ru
More information1 Derivation of Rate Equations from SingleCell Conductance (HodgkinHuxleylike) Equations
Physcs 171/271 Davd Klenfeld  Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from SngleCell Conductance (HodgknHuxleylke) Equatons We consder a network of many neurons, each of whch obeys
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 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 secondorder lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More information9 Derivation of Rate Equations from SingleCell Conductance (HodgkinHuxleylike) Equations
Physcs 171/271  Chapter 9R Davd Klenfeld  Fall 2005 9 Dervaton of Rate Equatons from SngleCell Conductance (HodgknHuxleylke) Equatons We consder a network of many neurons, each of whch obeys a set
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 informationIntroduction.  The Second Lyapunov Method.  The First Lyapunov Method
Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of
More information