Computational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
|
|
- Cathleen Jennifer Carr
- 5 years ago
- Views:
Transcription
1 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 equatons whch govern moton of common fluds such as water, honey, and so on. If you just want to see runnng code, t s on GtHub. Smulatons Due to the complex nature of the Naver-Stokes equatons, analytcal solutons range from dffcult to practcally mpossble, so the most common way to make use of the Naver-Stokes equatons s through smulaton and approxmaton. A number of smulaton methods exst, and n the next secton, we wll examne one of the algorthms often used n computer graphcs and nteractve applcatons. Smoothed Partcle Hydrodynamcs One of the most researched and most used methods for flud smulatons, especally n nteractve applcatons, s a partcle-based approach known as Smoothed Partcle Hydrodynamcs (SPH). SPH s, at the most basc level, an attempt to dscretze the Naver-Stokes equatons. In addton to the flud partcles themselves, a smoothng functon (called the smoothng kernel) s used whch spreads the propertes of each partcle over a contnuum. In fact, SPH s often called an nterpolaton method for partcle systems. We wll now derve and demonstrate the use of SPH for modelng a smple ncompressble Newtonan flud. Smoothng Kernels and Bass of SPH Snce SPH s essentally an nterpolaton method for partcle systems, one of the key deas behnd t s the ablty to evaluate a scalar quantty n any pont n space, even though the quantty only truly has a value at the dscrete partcles. In order to evaluate a quantty A at a pont r, we use the followng equaton: A( r) = j W ( r r j, h). Ths s a weghted summaton of all the partcles around the pont. s the value of A at the partcle, s the mass of the partcle, s the densty of the partcle, and W ( r r j, h) s the weghtng factor, and h s the maxmum dstance. The functon W s also called the smoothng kernel, because t smoothes the partcle values throughout contnuous space. The smoothng functon W has a few notable propertes. Frst of all, the functon should decrease as the dstance ncreases; ths s only logcal, snce the ponts nearby should have more nfluence than those far away. Second of all, the kernel should evaluate to 0 at dstances equal to or greater than h. Fnally, the smoothng kernel should be normalzed. A normalzed smoothng kernel satsfes the equaton W ( r, h)dr = 1. Dervatves of Quanttes The Naver-Stokes equatons often requre calculatng the frst or second dervatves of varous quanttes. Ths s partcularly easy n SPH, because the dervatve of a quantty s r A( r) = r j W ( r r j, h). Equvalently, Thus, r A( r) = j A( r) = j r W ( r r j, h). W ( r r j, h) 1
2 and 2 A( r) = j 2 W ( r r j, h). Calculatng Densty As you may have notced, the equaton for the SPH value of a scalar quantty ncludes both a mass and a densty term. Ths s because each partcle represents a small volume V = m ρ. The Naver-Stokes equatons for an ncompressble flud (such as water) requre that the densty term s constant; however, due to the nature of SPH, mantanng a constant densty s very dffcult. Snce n many applcatons accuracy s less mportant than realsm and speed, we gnore the requrement that densty s constant. We can later compensate for that by makng the densty restorng force (pressure) very strong. Snce densty ρ s non-constant, we may need to calculate densty at a gven poston. Snce densty s yet another ntensve property, we can apply the SPH equaton above to get = j W ( r r j, h) = j W ( r r j, h). The problem of densty s one of the bggest reasons not to use SPH for models whch requre precson or accuracy. SPH s meant to provde a realstc-lookng flud smulaton, at the cost of gnorng certan requrements, such as symmetry, conservaton of momentum, or constant densty. Smulaton The Naver-Stokes equatons can now be appled to the SPH partcles. We have two man equatons. The frst s the conservaton of mass equaton, v = 0. The second equaton stems from the conservaton of momentum, and states that ρ Dv Dt = p + µ 2 v + f. The frst equaton s guaranteed for partcle systems, because we have a constant amount of partcles, each wth a constant mass. Therefore, the total mass remans constant, and we can gnore the frst equaton. The second equaton dctates the moton of the partcles. The left hand term s the change n momentum. Therefore, we can express the acceleraton of each partcle as a = f ρ where f s the sum of all the force terms on the rght hand sde. In order to be able to smulate the partcle moton, we must express each of the force terms through a dscrete summaton. Force Terms: Pressure If we apply our SPH equaton to pressure, we arrve at the equaton f pressure = p( r) = p j W ( r r j, h). However, ths sn t symmetrc. Pressure should be symmetrc, snce every acton should have an equal and opposte reacton; however, f we examne a 2-partcle system, the pressure of partcle depends on the pressure of partcle j and vce versa, so symmetry s not guaranteed n the general case. In order to rectfy ths, we can use a dfferent pressure equaton, whch uses an average of partcle pressures to ensure symmetry: f pressure = p( r) = p + p j W ( r r j, h). 2 We must also remember that the only state a partcle keeps around s ts mass, densty, and velocty. Thus, we cannot just use pressure, we must calculate t. Usng the deal gas state equaton, we can say that p = kρ. 2
3 In order to make the smulaton more numercally stable, we can make pressure smaller by usng the dfference. Snce the equaton uses the gradent of pressure, the math s unchanged f we use the equaton p = k(ρ ρ o ) where ρ o s the rest densty. By dong ths substtuton, we reduce the value of p, and thus decrease the effects of floatng pont roundng errors ntroduced by the computatons. Force Terms: Vscosty Applyng SPH prncples to vscosty, we arrve at the equaton f vscosty = µ 2 v( r ) = µ v j 2 W ( r r j, h). However, we encounter the same problem as before, namely that the vscosty force s non-symmetrc. We can resolve t n a smlar manner, usng the dfference n veloctes nstead. The resultng equaton takes the form of f vscosty = µ v j v 2 W ( r r j, h). In physcal terms, ths equaton states that a partcle s accelerated n the drecton that ts envronment s movng. Thnkng about a floatng partcle n a stream of lqud, ths makes sense. Force Terms: Surface Tenson In a flud, molecules experence ntramolecular attracton forces. For a partcle embedded nsde the flud, the force s evened out; however, for partcles on the edge of the flud, the force does cause an acceleraton, and thus s non-neglgble. Snce each partcle s the cause of the ntramolecular attracton, we defne another quantty (called color for convenence) such that c s equal to 1 at every partcle and 0 otherwse. Therefore, the equaton that descrbes the SPH quantty for c would be c( r) = c j W ( r r j, h). Snce c s 1 at all partcle locatons, we can smplfy ths to say that c( r) = 1 W ( r r j, h). For physcal reasons that are beyond the scope of ths analyss, the surface tenson force s equal to f surface = σκ n = σ 2 c n n. We try to avod calculatng the force when n s small, because that causes problems because we are dvdng by a small floatng pont number. We only calculate the surface tenson force when n s greater than a gven preset. The fnal equaton for the surface tenson force (when we calculate t) would be Force Terms: External Forces f surface = σ 2 c c c. External forces are smply appled to every partcle as they would be to a normal mass. If a partcle touches (or goes nto) a sold object, t s reflected (pushed out wth t s velocty vector reflected over the object normal vector). Smoothng Kernels The smoothng kernel functons are requred to do the actual computatons for smulaton. The qualty of smulaton, as well as the speed and accuracy, hghly depend on the choce of smoothng kernels. In the orgnal SPH publcaton that nvestgated smulatng fluds n real-tme [2], the authors suggested the use of three dfferent kernels for dfferent thngs. Note that the smoothng kernels are pece-wse functons; the functons presented subsequently defne the values of W for 0 r h. In other cases, W = 0. Note that the followng smoothng kernels are not normalzed, and should be pror to use. 3
4 General Smoothng Kernel The general smoothng kernel, used n all but two cases, s a sxth-degree polynomal wth the equaton W g (r, h) = πh 9 ( h 2 r 2) 3. One of the benefts of ths emprcally desgned kernel s that t uses r 2, whch means that we never need to calculate r. Calculatng r would be an expensve operaton, because t nvolves a square root; calculatng r 2, however, nvolves only multplcaton, and s thus much faster. Pressure Kernel For calculatons nvolvng pressure, the general kernel becomes zero when partcles are very close together. Ths s not what fluds actually act lke - n fact, as partcles come closer together, they should repel more due to pressure nstead of less. Thus, Desbrun [6] suggests a dfferent smoothng kernel of the form Vscosty Kernel W p (r, h) = 15 πh 6 (h r)3. The fnal smoothng kernel we are usng s specfc to the vscosty term. If we were to use the general smoothng term for vscosty, two partcles whch are very close to each other could possbly ncrease n speed. Ths s ncorrect, snce vscosty s a frcton force, whch should not be ncreasng the energy n the system. Thus, the authors of the orgnal flud smulaton SPH paper suggest usng the kernel of Ongong Research W v (r, h) = r3 2h 3 + r2 h 2 + h 2r 1. Although the Naver-Stokes equatons were orgnally formulated nearly 200 years ago, flud mechancs s a rapdly evolvng feld. The advent of computer technology has allowed physcsts to accurately model the behavour of fluds through teratve approxmatons of ther behavour, makng the study of the Naver-Stokes equatons stll more mportant. Current research can be subdvded nto two categores: one form of research s theoretcal, seekng mathematcal proof of the exstence and smoothness of the soluton to Naver-Stokes; the other branch of ongong research s more focused on practcal applcaton of Naver-Stokes to real-world problems. Though n some cases the dstncton s purely artfcal, t s nonetheless useful n dscussng the varous topcs currently beng nvestgated. Theoretcal Research The fundamental theoretcal problem of Naver-Stokes s one of the seven Mllenum Prze Problems. The Mllenum Prze Problems are seven mathematcal problems put forth by the Clay Mathematcs Insttute n the year 2000; the frst person to offer a vald soluton to any of the problems s promsed a mllon-dollar reward. The problem of Naver-Stokes, as stated by the Clay Mathematcs Insttute, s as follows: Prove or gve a counter-example of the followng statement: In three space dmensons and tme, gven an ntal velocty feld, there exsts a vector velocty and a scalar pressure feld, whch are both smooth and globally defned, that solve the Naver-Stokes equatons. Though ths statement tself has not been proven or ndcated false, attempts at proof have come close. It s known that for Naver-Stokes n two dmensons, the solutons exst and are smooth. In addton, t s known that for low ntal veloctes, smooth solutons to 3-dmensonal Naver-Stokes exst; however, for hgher ntal veloctes, the exstence of smooth solutons has not been proven. Practcal Research Though the mathematcal problem posed by Naver-Stokes has not yet been resolved, computer-aded approxmatons are sutable for many purposes. From the desgn of aerodynamc arcraft to smulaton of water for Pxar flms, flud smulatons are hghly used n ndustry. Much current research n the feld of computatonal 4
5 flud dynamcs (CFD) attempts to fnd faster or more-accurate approxmaton algorthms n order to make these smulatons more useful stll. The methods used n CFD can be broken down nto two categores: dscretzaton models and turbulence models. Dscretzaton models attempt to break the flud nto a dscrete set of elements, and teratvely smulate the moton of the flud through tme usng the Naver-Stokes equatons. Turbulence models deal wth turbulence, a chaotc type of flud flow, whch s nearly mpossble to smulate as wth the dscretzaton models; nstead, turbulence models attempt to track or obtan the value of specfc propertes of a flud over tme. Though the basc equatons used are the same, approxmaton methods often take dfferent approaches, dependng on ther objectve. Some methods, especally ones amed towards very fast or nteractve flud smulaton, use partcle based modelng. A volume of flud s subdvded nto a fnte number of partcles, and each partcle s subject to the velocty felds dctated by the Naver-Stokes equatons as well as external forces. Ths model typcally smplfes computaton and mplementaton, and s commonly used n computer graphcs and computer games. The partcle nature s usually hdden by nterpolatng partcle postons through a process called pont splattng. Smoothed Partcle Hydrodynamcs, the method nvestgated n the prevous secton, falls under that category. Other methods, however, use a three-dmensonal grd (a mesh ) to model the flud. The nternals and boundares are tracked through teratve Naver-Stokes approxmatons; the beneft of ths method s a hgher precson and hgher flexblty n terms of nteracton wth other bodes, at the cost of computatonal speed. 5
Numerical 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 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 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 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 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 informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationCONTROLLED FLOW SIMULATION USING SPH METHOD
HERI COADA AIR FORCE ACADEMY ROMAIA ITERATIOAL COFERECE of SCIETIFIC PAPER AFASES 01 Brasov, 4-6 May 01 GEERAL M.R. STEFAIK ARMED FORCES ACADEMY SLOVAK REPUBLIC COTROLLED FLOW SIMULATIO USIG SPH METHOD
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 informationBindel, Spring 2014 Applications of Parallel Computers (CS 5220)
Bndel, Sprng 2014 Applcatons of Parallel Computers (CS 5220) 1 Dervng SPH The Naver-Stokes equatons wth gravty are ρa = p + µ 2 v + ρg. The acceleraton s the materal dervatve of velocty, and we usually
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 informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space 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 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 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 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 information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
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 informationBuckingham s pi-theorem
TMA495 Mathematcal modellng 2004 Buckngham s p-theorem Harald Hanche-Olsen hanche@math.ntnu.no Theory Ths note s about physcal quanttes R,...,R n. We lke to measure them n a consstent system of unts, such
More informationTemperature. Chapter Heat Engine
Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the
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 informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
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 informationWeek 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product
The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the
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 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 informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
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 informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
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 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 informationTurbulent Flow. Turbulent Flow
http://www.youtube.com/watch?v=xoll2kedog&feature=related http://br.youtube.com/watch?v=7kkftgx2any http://br.youtube.com/watch?v=vqhxihpvcvu 1. Caothc fluctuatons wth a wde range of frequences and
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
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 informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
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 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 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 informationN-Body Simulation. Typical uncertainty: π = 4 Acircle/Asquare! 4 ncircle/n. πest π = O(n
N-Body Smulaton Solvng the CBE wth a 6-D grd takes too many cells. Instead, we use a Monte-Carlo method. Example: Monte-Carlo calculaton of π. Scatter n ponts n square; count number ncrcle fallng wthn
More informationTHE SUMMATION NOTATION Ʃ
Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the
More informationFoundations of Arithmetic
Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an
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 informationNewton s Laws of Motion
Chapter 4 Newton s Laws of Moton 4.1 Forces and Interactons Fundamental forces. There are four types of fundamental forces: electromagnetc, weak, strong and gravtatonal. The frst two had been successfully
More informationLecture 20: Noether s Theorem
Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external
More informationSMOOTHED PARTICLE HYDRODYNAMICS METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM
The 5th Internatonal Symposum on Computatonal Scences (ISCS) 28 May 2012, Yogyakarta, Indonesa SMOOTHED PARTICLE HYDRODYNAMICS METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM Dede Tarwd 1,2 1 Graduate School
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 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 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 informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationMA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials
MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have
More informationGravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)
Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng
More informationPhysics 207 Lecture 6
Physcs 207 Lecture 6 Agenda: Physcs 207, Lecture 6, Sept. 25 Chapter 4 Frames of reference Chapter 5 ewton s Law Mass Inerta s (contact and non-contact) Frcton (a external force that opposes moton) Free
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
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 informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
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 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 informationThermodynamics General
Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
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 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 informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
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 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 informationLinear Momentum. Center of Mass.
Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html
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 informationTurbulence and its Modelling
School of Mechancal Aerospace and Cvl Engneerng 3rd Year Flud Mechancs Introducton In earler lectures we have consdered how flow nstabltes develop, and noted that above some crtcal Reynolds number flows
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 information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6
REVIEW of Lecture 5 2.29 Numercal Flud Mechancs Fall 2011 Lecture 6 Contnuum Hypothess and conservaton laws Macroscopc Propertes Materal covered n class: Dfferental forms of conservaton laws Materal Dervatve
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
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 informationFirst Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.
Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act
More informationRotor Noise Modeling Kenneth S. Brentner Penn State University
Rotor Nose Modelng Kenneth S. Brentner Penn State Unversty Joby Avaton S4 www.jobyavaton.com 2018 Kenneth S. Brentner. All rghts reserved. 5 th Transformatve Vertcal Flght Workshop, January 18-19, 2018
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 informationwhere v means the change in velocity, and t is the
1 PHYS:100 LECTURE 4 MECHANICS (3) Ths lecture covers the eneral case of moton wth constant acceleraton and free fall (whch s one of the more mportant examples of moton wth constant acceleraton) n a more
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationThermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)
hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In
More informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
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 informationComputational Astrophysics
Computatonal Astrophyscs Solvng for Gravty Alexander Knebe, Unversdad Autonoma de Madrd Computatonal Astrophyscs Solvng for Gravty the equatons full set of equatons collsonless matter (e.g. dark matter
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 informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
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 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. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
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 informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationOn the correction of the h-index for career length
1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat
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 information