Understanding DPMFoam/MPPICFoam
|
|
- Carmel Golden
- 5 years ago
- Views:
Transcription
1 Understanding DPMFoam/MPPICFoam Jeroen Hofman March 18, 2015 In this document I intend to clarify the flow solver and at a later stage, the article-fluid and article-article interaction forces as imlemented in DPMFoam in OenFoam 2.3.x. 1 The Flow Solver In this section we first set u the incomressible Navier-Stokes equations for the gas flow solver, then we describe the derivation of the ressure equation in semi-discretized form, which is used in the PIMPLE algorithm emloyed by OenFOAM in the last section. 1.1 Incomressible Navier-Stokes We start with the incomressible Navier-Stokes equation for laminar flow with gas volume fraction ɛ, constant viscosity µ, gas density ρ, gravity g, stress tensor τ and momentum transfer term F: ρ t + ρ ɛuu) ρ ɛτ = + ρɛg F 1) where the stress tensor with unit tensor δ and kinematic viscosity ν = µ ρ is given by: τ = ν ) u + u) T + 2 ν u) δ 2) 3 The interhase momentum transfer term F contains both momentum transfer due to drag and buoyancy [2], details will be given in section 1.2. Therefore it is not fully a source term since the drag is roortional to u. When dividing by the density and defining the kinematic ressure as P = ρ we obtain: t + ɛuu) ɛτ = P + ɛg F ρ 3) 1
2 We also define the continuity equation as: ɛ + ɛu) = 0 4) t 1.2 Momentum transfer In this section we look at different ways of defining momentum transfer to either the article or the fluid hase. By switching from the discrete article descrition to a continuum descrition we try to link the different ways of defining the momentum transfer term since regardless of definition of F equation 3 the resulting PDEs should be equal. We will first look at the way the momentum transfer term is usually defined as a combination of buoyancy and drag for articles in a comutational cell: F = f = 1 V f drag V i ) = F drag 1 V V i 5) where V is the volume of the comutational cell, V the volume of the article and i the locally averaged ressure gradient. In a continuum descrition we can write 1 V V i = 1 ɛ) and rewrite equation 3: t + ɛuu) ɛτ = ɛ + ɛg F drag ρ 6) Another way of incororating buoyancy in stead of via the ressure gradient is by exlicitly adding a local acceleration term lus an adjustment to the drag [6], section 2.5.1): F = f = 1 V fdrag where [ ] Dt gravity field is uniform and we use 1 V ɛ ρv g )) = F drag + 1 Dt ɛ V ρv g is the locally averaged fluid velocity at the location of article. term and we can write equation 3 as: [ ] ) Dt Using that the ρv g = 1 ɛ)ρg we take this term out of the interaction 7) t + ɛuu) ɛτ = P + g F drag ɛρ 1 V [ ] V Dt 8) Furthermore realizing that, by using continuity in equation 4 and using roduct rules, we have: t + ɛuu) = ɛ Dt = ɛ u + ɛu u 9) t 2
3 and using that in the continuum descrition we have 1 V second form of the momentum equation: V Dt ɛτ = + g F drag ɛρ [ ] Dt = 1 ɛ) Dt we obtain a 10) Exressions 6 and 10 should be the same when the latter is multilied by ɛ, however there are clearly not because of the treatment of the stress term. I have sent quite some time digging into this roblem and it turns out this is a very subtle roblem arising from the recise definitions of both fluid and solid equations in a continuum descrition. See for instance [7], where the authors show that many well-cited) aers even mess this u. Let us first consider how to fix this roblem: 1. Omitting the viscuous stress is the aroach that is taken in the MP-PIC aers [8], this will make equations 6 and 10 equal. Arguably this is valid as the driving factors are often drag and article ressure, and esecially for gases the viscosity is very small. 2. The aroach alied in both [7] and [6] is changing the ansatz of the derivation, namely the momentum equation; equation 3. Similar to accounting for ressure in both the gas hase momentum equation and the article momentum transfer term, we can instead also slit the viscuous stress in this way. We then obtain: t + ɛuu) τ = P + ɛg F 11) ρ F = F drag 1 V i + τ i ) 12) V Combining the above equations and using the continuum descrition we end u with: ɛ Dt ɛ τ = ɛ + ɛg F drag ρ 13) Similarly if we use equation 11 we can follow the same stes when incororating buoyancy via the local acceleration as we did before and then we obtain the following relacement of equation 10: Dt τ = + g F drag ɛρ 14) Now we see that the aroaches do line u. Some comments can be made here: first of all in this aroach we end u with the hase fraction ɛ outside the divergence term of the viscuous stress. Arguably this is valid when ɛ 0 since ɛτ = ɛ τ + τ ɛ, however this is clearly not the case in for instance fluidized beds. Furthermore a simle derivation using control volumes and flux evaluation will not give you an answer to this question. Instead, one can theoretically derive the momentum equations in a continuum framework) in which there are two aroaches [9], chater 3): the Jackson aroach, in which the viscuous stress 3
4 has ɛ outside the divergence term, but also includes 1 ɛ) τ in the article momentum equations; and secondly the Ishii aroach, in which ɛ is inside the divergence term of the vicuous stress but not included in the article equations. These aroaches also known as model A and model B in literature) cause a high amount of confusion and often are wrongly imlemented in for instance TFM, however the difference can be subtle as again viscuous stresses are deemed not so imortant. If we omit viscuous stress the models are the same. The confusion is further enhanced when the article hase is discrete, since we introduce an extra layer of comlexity. Returning to the key-oint of this reort, namely understanding OenFOAM, it seems to use a mixed aroach. The momentum equation is as equation 3 and it uses the locally averaged fluid velocity as substitute for the ressure gradient, so eventually equation 8 is used, where the drag F drag is defined as: 1 1 ɛ)v V β u u ) 15) where β is defined according to the Ergun-Wen-Yu drag model [9], u is the article velocity and u is the fluid velocity interolated to the article osition. 1.3 Discretization There are many choices regarding discretization of the terms in equation 3, which we will not discuss in deth here. We only highlight that when looking at the convection term ɛuu) we see that the discretized system would be non-linear. A ossible solution would be to solve this system using non-linear solvers but this is in general quite exensive and imractical. A second solution would be to linearize the convection term. Using the discretized version of Gauss theorem [1] we can write for a control volume V P with faces f centered around a oint P, where all variables are defined at cell centres: V P ɛuu)dv = f = f = f S ɛuu) f S ɛu) f u f F u f 16) where we have defined F = S ɛu) f as the hase) mass face flux and values are assumed to be interolated from cell centers 1. F has to satisfy continuity, which is equation 4. 1 We make the assumtion that ab) f = a f b f, which is not necessarily true for interolated values, but often assumed to make interolations easier, see for instance the discussion in [3] 4
5 By linearizing the convection term it follows that we define the fluxes F as a time-lagged quantity, i.e. we assume we have an existing velocity flux field F satisfying equation 4. We can now write the discretized version of the convection term as a linear combination of the velocity: V P ɛuu)dv = f F u f u P + N a N u N 17) where the coefficients a N and contain the known fluxes F and the velocity at the face is assumed to be interolated from the control volume and neighboring control volumes. This imlicitly also assumes the hase fraction ɛ to be known when discretizing u. 1.4 Derivation of the ressure equation Since equation 3 can be discretized in a linear fashion following the treatment of the convection term above, we can semi-discretize the momentum equation. F can be slit as a combination of a contribution to the velocity and a source term, i.e. F = f P u P + f. Details of the slit were discussed in section 1.2. For now it is sufficient that the slit can be made. The semi-discretized equation has the following form, combining the diagonal coefficients in one coefficient : a P f ) P u P = ρ N a N u N + a 0 u 0 P + g + f ρ P u P = Hu) P u P = Hu) P 18) where the subscrit 0 indicates the revious value. Hu) contains all off-diagonal contributions of the linear system, including source terms. The semi-discretized version of the continuity equation equation 4) can be written as: ɛ t + f S ɛu) f = 0 19) where we assume some exlicit discretization of the temoral term, not shown here. We first exress the velocity of the cell face as: ) Hu) u f = f ) P 20) where P, Hu) and are interolated to the cell faces. Interolation of P is erformed by simle differencing of both control volumes adjacent to the face, i.e. P ) f P P P E x, where control volumes with centres P and E are adjacent to face f and x is the distance between P and E. By 5
6 directly calculating the gradient of ressure on the cell faces by the ressure in adjacent cells and not interolating the gradient itself, checkerboard roblems for the ressure are avoided. By combining equation 19 and 20 we obtain the ressure equation: f S ɛ f P ) f = ɛ t + f ) Hu) S ɛ f 21) The flux F is defined as: [ Hu) ) F = S ɛu) f = S ɛ f f ) ] P 22) This flux is guaranteed to be conservative by construction. 1.5 Rhie-Chow interolation It is well-known that solving the incomressible Navier-Stokes equations on a collocated grid gives rise to the so-called checkerboard roblem when using central differencing, where the solutions for velocity and ressure field are no longer unique u to a constant in the case of ressure), see [4], chater 7.5. A well-known cure to this roblem is the so called Rhie-Chow interolation rocedure, where the face velocity is adjusted as: ) 1 [ ] u f = u f A Pf P f 23) The overline denotes interolated values. In other words, the interolated velocity is corrected by the difference in the gradient of the ressure at the cell face and the interolated gradient. It can be shown that this correction is roortional to the third derivative of the ressure, corresonding to a fourth order correction in the continuity equation. We will now show something similar haens when deriving the ressure equation as we did in the revious section. Consider again equation 18 where we exlicitly state that the gradient of ressure is evaluated at oint P : u P = Hu) P ) P 24) By defining the cell-faced velocity as in equation 20 and combining with the above we have: u f = u P + P ) ) ) P P 25) Writing out the interolation oeration with a f = a E ) we get an exression that is very similar to equation 23: 6
7 u f = 1 2 u P + u E ) + 1 [ P ) P + P ) E ] 1 P P P E 2a f a f x = 1 2 u P + u E ) + 1 [ PEE P P + P ] E P W 1 P P P E 2a f 2 x 2 x a f x = 1 2 u P + u E ) + 1 4a f x [P EE 3P E + 3P P P W ] 26) the last term on the right hand side of the above equation is indeed the discretization of the third derivative of the ressure, the same as the Rhie Chow interolation. 1.6 PIMPLE in OenFOAM We are now ready to describe the combined PISO + SIMPLE algorithm that is used in OenFOAM. The following stes can be identified: 1. Set u the discretized equation for u on the grid without any source terms, i.e. we set u a coefficient matrix according to the first line of equation 18, excluding the source terms. 2. Relax this equation. 3. Solve equation 18 with an initial ressure field estimate P guess to obtain the momentum redictor u at the cell centres. P guess is usually chosen equal to the ressure field in the revious time-ste. 4. Using the redictor u we can assemble the oerator Hu ) and interolate both and Hu ) to the cell faces. 5. Solve the ressure equation equation 21) on the cell faces. 6. Udate the flux field at the cell faces using the obtained ressure equation 22). 7. Relax the ressure. 8. Udate the velocity at the cell centres using the obtained ressure equation 18). 9. Udate the boundary conditions for consistency. Stes 1-9 can be reeated any number of times. This, together with the relaxation, is the essence of the SIMPLE algorithm [5]. The momentum correction stes 4-9) can also be reeated any number of times, this is the essence of the PISO algorithm [5]. OenFOAM defines the PIMPLE SIMPLE + PISO) algorithm alying both loos, including relaxation. Combining the methods in general has roven to give faster convergence. 7
8 2 Particle contributions in MPPICFoam 2.1 Particle-Particle forces 3 MPPICFoam/DPMFoam code 3.1 Parcel Evolution 3.2 Flow Solving References [1] htt://owerlab.fsb.hr/ed/kturbo/oenfoam/docs/hrvojejasakphd.df [2] Andrews, M.J., O Rourke, P.J., 1996, The multihase article-in-cell MP-PIC) method for dense articulate flows, Int. J. Multihase Flow 22, [3] htt:// [4] Ferziger, J. H. and Peric, M., Comutational Methods for Fluid Dynamics, 2nd ed., Sringer- Verlag 2001) [5] Versteeg H.K., Malalasekera W., An introduction to comutational fluid dynamics, the finite volume method 2nd ed.), [6] htt://gradworks.umi.com/33/03/ html [7] Z. Y. ZHOU, S. B. KUANG, K. W. CHU and A. B. YU 2010). Discrete article simulation of article-fluid flow: model formulations and their alicability. Journal of Fluid Mechanics, 661, doi: /s x [8] O Rourke, P.J., Snider, D.M., 2010, An imroved collision daming time for MP-PIC calculations of dense article flows with alications to olydiserse sedimenting beds and colliding jet articles, Chemical Engineering Science 65, [9] Wachem, B., Derivation, imlementation, and validation of comuter simulation models for gas-solid fluidized beds, PhD thesis of our grou,
5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:
5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum
More informationModeling Volumetric Coupling of the Dispersed Phase using the Eulerian-Lagrangian Approach
ILASS Americas 21st Annual Conference on Liquid Atomization and Sray Systems, Orlando, FL, May 28 Modeling Volumetric Couling of the Disersed Phase using the Eulerian-Lagrangian Aroach Ehsan Shams and
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationJohn Weatherwax. Analysis of Parallel Depth First Search Algorithms
Sulementary Discussions and Solutions to Selected Problems in: Introduction to Parallel Comuting by Viin Kumar, Ananth Grama, Anshul Guta, & George Karyis John Weatherwax Chater 8 Analysis of Parallel
More informationCFD Modelling of Mass Transfer and Interfacial Phenomena on Single Droplets
Euroean Symosium on Comuter Arded Aided Process Engineering 15 L. Puigjaner and A. Esuña (Editors) 2005 Elsevier Science B.V. All rights reserved. CFD Modelling of Mass Transfer and Interfacial Phenomena
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationAntony Jameson. Stanford University Aerospace Computing Laboratory Report ACL
Formulation of Kinetic Energy Preserving Conservative Schemes for Gas Dynamics and Direct Numerical Simulation of One-dimensional Viscous Comressible Flow in a Shock Tube Using Entroy and Kinetic Energy
More informationNUMERICAL ANALYSIS OF THE IMPACT OF THE INLET AND OUTLET JETS FOR THE THERMAL STRATIFICATION INSIDE A STORAGE TANK
NUMERICAL ANALYSIS OF HE IMAC OF HE INLE AND OULE JES FOR HE HERMAL SRAIFICAION INSIDE A SORAGE ANK A. Zachár I. Farkas F. Szlivka Deartment of Comuter Science Szent IstvÆn University Æter K. u.. G d llı
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationCaveats on the Implementation of the Generalized Material Point Method
Coyright c 2008 Tech Science Press CMES, vol.31, no.2,.85-106, 2008 Caveats on the Imlementation of the Generalized Material Point Method O. Buzzi 1, D. M. Pedroso 2 and A. Giacomini 1 Abstract: The material
More informationNumerical Simulation of Particle Concentration in a Gas Cyclone Separator *
2007 Petroleum Science Vol.4 No.3 Numerical Simulation of Particle Concentration in a Gas Cyclone Searator * Xue Xiaohu, Sun Guogang **, Wan Gujun and Shi Mingxian (School of Chemical Science and Engineering,
More informationSolutions to Problem Set 5
Solutions to Problem Set Problem 4.6. f () ( )( 4) For this simle rational function, we see immediately that the function has simle oles at the two indicated oints, and so we use equation (6.) to nd the
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationImplementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis
CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationWhen solving problems involving changing momentum in a system, we shall employ our general problem solving strategy involving four basic steps:
10.9 Worked Examles 10.9.1 Problem Solving Strategies When solving roblems involving changing momentum in a system, we shall emloy our general roblem solving strategy involving four basic stes: 1. Understand
More informationParticipation Factors. However, it does not give the influence of each state on the mode.
Particiation Factors he mode shae, as indicated by the right eigenvector, gives the relative hase of each state in a articular mode. However, it does not give the influence of each state on the mode. We
More informationNumerical Methods for Particle Tracing in Vector Fields
On-Line Visualization Notes Numerical Methods for Particle Tracing in Vector Fields Kenneth I. Joy Visualization and Grahics Research Laboratory Deartment of Comuter Science University of California, Davis
More informationCFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS
CFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS M. BORGHI - M. MILANI Diartimento di Scienze dell Ingegneria Università degli Studi di Modena Via Cami, 213/b 41100 Modena E-mail: borghi@omero.dsi.unimo.it
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationNotes on pressure coordinates Robert Lindsay Korty October 1, 2002
Notes on ressure coordinates Robert Lindsay Korty October 1, 2002 Obviously, it makes no difference whether the quasi-geostrohic equations are hrased in height coordinates (where x, y,, t are the indeendent
More informationThe Role of Momentum Interpolation Mechanism of the Roe. Scheme in the Shock Instability
The Role of Momentum Interolation Mechanism of the Roe Scheme in the Shock Instabilit Xue-song Li Ke Laborator for Thermal Science and Power Engineering of Ministr of Education, Deartment of Thermal Engineering,
More informationAn Analysis of Reliable Classifiers through ROC Isometrics
An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit
More informationWave Drift Force in a Two-Layer Fluid of Finite Depth
Wave Drift Force in a Two-Layer Fluid of Finite Deth Masashi Kashiwagi Research Institute for Alied Mechanics, Kyushu University, Jaan Abstract Based on the momentum and energy conservation rinciles, a
More informationUniversal Finite Memory Coding of Binary Sequences
Deartment of Electrical Engineering Systems Universal Finite Memory Coding of Binary Sequences Thesis submitted towards the degree of Master of Science in Electrical and Electronic Engineering in Tel-Aviv
More informationModel checking, verification of CTL. One must verify or expel... doubts, and convert them into the certainty of YES [Thomas Carlyle]
Chater 5 Model checking, verification of CTL One must verify or exel... doubts, and convert them into the certainty of YES or NO. [Thomas Carlyle] 5. The verification setting Page 66 We introduce linear
More informationON POLYNOMIAL SELECTION FOR THE GENERAL NUMBER FIELD SIEVE
MATHEMATICS OF COMPUTATIO Volume 75, umber 256, October 26, Pages 237 247 S 25-5718(6)187-9 Article electronically ublished on June 28, 26 O POLYOMIAL SELECTIO FOR THE GEERAL UMBER FIELD SIEVE THORSTE
More informationMODULAR LINEAR TRANSVERSE FLUX RELUCTANCE MOTORS
MODULAR LINEAR TRANSVERSE FLUX RELUCTANCE MOTORS Dan-Cristian POPA, Vasile IANCU, Loránd SZABÓ, Deartment of Electrical Machines, Technical University of Cluj-Naoca RO-400020 Cluj-Naoca, Romania; e-mail:
More informationComputer arithmetic. Intensive Computation. Annalisa Massini 2017/2018
Comuter arithmetic Intensive Comutation Annalisa Massini 7/8 Intensive Comutation - 7/8 References Comuter Architecture - A Quantitative Aroach Hennessy Patterson Aendix J Intensive Comutation - 7/8 3
More informationEffect of geometry on flow structure and pressure drop in pneumatic conveying of solids along horizontal ducts
Journal of Scientific LAÍN & Industrial SOMMERFELD Research: PNEUMATIC CONVEYING OF SOLIDS ALONG HORIZONTAL DUCTS Vol. 70, February 011,. 19-134 19 Effect of geometry on flow structure and ressure dro
More informationOptimization of Gear Design and Manufacture. Vilmos SIMON *
7 International Conference on Mechanical and Mechatronics Engineering (ICMME 7) ISBN: 978--6595-44- timization of Gear Design and Manufacture Vilmos SIMN * Budaest Universit of Technolog and Economics,
More informationHomogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers
Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas
More informationRecent Developments in Multilayer Perceptron Neural Networks
Recent Develoments in Multilayer Percetron eural etworks Walter H. Delashmit Lockheed Martin Missiles and Fire Control Dallas, Texas 75265 walter.delashmit@lmco.com walter.delashmit@verizon.net Michael
More informationRobust Predictive Control of Input Constraints and Interference Suppression for Semi-Trailer System
Vol.7, No.7 (4),.37-38 htt://dx.doi.org/.457/ica.4.7.7.3 Robust Predictive Control of Inut Constraints and Interference Suression for Semi-Trailer System Zhao, Yang Electronic and Information Technology
More informationStatistical Models approach for Solar Radiation Prediction
Statistical Models aroach for Solar Radiation Prediction S. Ferrari, M. Lazzaroni, and V. Piuri Universita degli Studi di Milano Milan, Italy Email: {stefano.ferrari, massimo.lazzaroni, vincenzo.iuri}@unimi.it
More informationAn-Najah National University Civil Engineering Departemnt. Fluid Mechanics. Chapter [2] Fluid Statics
An-Najah National University Civil Engineering Deartemnt Fluid Mechanics Chater [2] Fluid Statics 1 Fluid Statics Problems Fluid statics refers to the study of fluids at rest or moving in such a manner
More informationGOOD MODELS FOR CUBIC SURFACES. 1. Introduction
GOOD MODELS FOR CUBIC SURFACES ANDREAS-STEPHAN ELSENHANS Abstract. This article describes an algorithm for finding a model of a hyersurface with small coefficients. It is shown that the aroach works in
More informationFault Tolerant Quantum Computing Robert Rogers, Thomas Sylwester, Abe Pauls
CIS 410/510, Introduction to Quantum Information Theory Due: June 8th, 2016 Sring 2016, University of Oregon Date: June 7, 2016 Fault Tolerant Quantum Comuting Robert Rogers, Thomas Sylwester, Abe Pauls
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationYixi Shi. Jose Blanchet. IEOR Department Columbia University New York, NY 10027, USA. IEOR Department Columbia University New York, NY 10027, USA
Proceedings of the 2011 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelsach, K. P. White, and M. Fu, eds. EFFICIENT RARE EVENT SIMULATION FOR HEAVY-TAILED SYSTEMS VIA CROSS ENTROPY Jose
More informationCode_Aster. Connection Harlequin 3D Beam
Titre : Raccord Arlequin 3D Poutre Date : 24/07/2014 Page : 1/9 Connection Harlequin 3D Beam Summary: This document exlains the method Harlequin develoed in Code_Aster to connect a modeling continuous
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More information2-D Analysis for Iterative Learning Controller for Discrete-Time Systems With Variable Initial Conditions Yong FANG 1, and Tommy W. S.
-D Analysis for Iterative Learning Controller for Discrete-ime Systems With Variable Initial Conditions Yong FANG, and ommy W. S. Chow Abstract In this aer, an iterative learning controller alying to linear
More informationChapter 10: Flow Flow in in Conduits Conduits Dr Ali Jawarneh
Chater 10: Flow in Conduits By Dr Ali Jawarneh Hashemite University 1 Outline In this chater we will: Analyse the shear stress distribution across a ie section. Discuss and analyse the case of laminar
More informationReal Analysis 1 Fall Homework 3. a n.
eal Analysis Fall 06 Homework 3. Let and consider the measure sace N, P, µ, where µ is counting measure. That is, if N, then µ equals the number of elements in if is finite; µ = otherwise. One usually
More information0.6 Factoring 73. As always, the reader is encouraged to multiply out (3
0.6 Factoring 7 5. The G.C.F. of the terms in 81 16t is just 1 so there is nothing of substance to factor out from both terms. With just a difference of two terms, we are limited to fitting this olynomial
More informationON THE NORM OF AN IDEMPOTENT SCHUR MULTIPLIER ON THE SCHATTEN CLASS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 000-9939XX)0000-0 ON THE NORM OF AN IDEMPOTENT SCHUR MULTIPLIER ON THE SCHATTEN CLASS WILLIAM D. BANKS AND ASMA HARCHARRAS
More informationRole of Momentum Interpolation Mechanism of the Roe. Scheme in Shock Instability
Role of Momentum Interolation Mechanism of the Roe Scheme in Shock Instability Xiao-dong Ren 1,2 Chun-wei Gu 1 Xue-song Li 1,* 1. Key Laboratory for Thermal Science and Power Engineering of Ministry of
More informationMeshless Methods for Scientific Computing Final Project
Meshless Methods for Scientific Comuting Final Project D0051008 洪啟耀 Introduction Floating island becomes an imortant study in recent years, because the lands we can use are limit, so eole start thinking
More informationComparison Of 2D Numerical Schemes For Modelling Supercritical And Transcritical Flows Along Urban Floodplains
City University of New York (CUNY) CUNY Academic Works International Conference on Hydroinformatics 8-1-014 Comarison Of D Numerical Schemes For Modelling Suercritical And Transcritical Flows Along Urban
More informationSystematizing Implicit Regularization for Multi-Loop Feynman Diagrams
Systematizing Imlicit Regularization for Multi-Loo Feynman Diagrams Universidade Federal de Minas Gerais E-mail: alcherchiglia.fis@gmail.com Marcos Samaio Universidade Federal de Minas Gerais E-mail: msamaio@fisica.ufmg.br
More informationProbability Estimates for Multi-class Classification by Pairwise Coupling
Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics
More informationu = 1 (B 2 + E2 E B (16.2) + N = j E (16.3) One might be tempted to put u and N into a 4-vector N and write the equation in the form
Chater 6 Energy-momentum tensor (Version., 3 November 7 Earlier, we obtained for the energy density and flux u = (B + E µ c (6. We also had a continuity equation N = µ E B (6. u t + N = j E (6.3 One might
More informationNotes on Instrumental Variables Methods
Notes on Instrumental Variables Methods Michele Pellizzari IGIER-Bocconi, IZA and frdb 1 The Instrumental Variable Estimator Instrumental variable estimation is the classical solution to the roblem of
More informationAsymptotically Optimal Simulation Allocation under Dependent Sampling
Asymtotically Otimal Simulation Allocation under Deendent Samling Xiaoing Xiong The Robert H. Smith School of Business, University of Maryland, College Park, MD 20742-1815, USA, xiaoingx@yahoo.com Sandee
More informationLORENZO BRANDOLESE AND MARIA E. SCHONBEK
LARGE TIME DECAY AND GROWTH FOR SOLUTIONS OF A VISCOUS BOUSSINESQ SYSTEM LORENZO BRANDOLESE AND MARIA E. SCHONBEK Abstract. In this aer we analyze the decay and the growth for large time of weak and strong
More informationCSE 599d - Quantum Computing When Quantum Computers Fall Apart
CSE 599d - Quantum Comuting When Quantum Comuters Fall Aart Dave Bacon Deartment of Comuter Science & Engineering, University of Washington In this lecture we are going to begin discussing what haens to
More informationSTABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS
STABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS Massimo Vaccarini Sauro Longhi M. Reza Katebi D.I.I.G.A., Università Politecnica delle Marche, Ancona, Italy
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationChapter 8 Internal Forced Convection
Chater 8 Internal Forced Convection 8.1 Hydrodynamic Considerations 8.1.1 Flow Conditions may be determined exerimentally, as shown in Figs. 7.1-7.2. Re D ρumd μ where u m is the mean fluid velocity over
More informationWe E Multiparameter Full-waveform Inversion for Acoustic VTI Medium with Surface Seismic Data
We E16 4 Multiarameter Full-waveform Inversion for Acoustic VI Medium with Surface Seismic Data X. Cheng* (Schlumberger) K. Jiao (Schlumberger) D. Sun (Schlumberger) & D. Vigh (Schlumberger) SUMMARY In
More informationTransport Phenomena Coupled by Chemical Reactions in Methane Reforming Ducts
I. J. Trans. Phenomena, Vol. 11,. 39 50 Rerints available directly from the ublisher Photocoying ermitted by license only 009 Old City Publishing, Inc. Published by license under the OCP Science imrint,
More informationTurbulent Flow Simulations through Tarbela Dam Tunnel-2
Engineering, 2010, 2, 507-515 doi:10.4236/eng.2010.27067 Published Online July 2010 (htt://www.scirp.org/journal/eng) 507 Turbulent Flow Simulations through Tarbela Dam Tunnel-2 Abstract Muhammad Abid,
More information+++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic +++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic
Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic Daniel Leitner, Johannes Krof, Günther Zauner, TU Vienna, Austria, dleitner@osiris.tuwien.ac.at Yuri Karov, Yuri Senichenkov, Yuri
More informationNotes on duality in second order and -order cone otimization E. D. Andersen Λ, C. Roos y, and T. Terlaky z Aril 6, 000 Abstract Recently, the so-calle
McMaster University Advanced Otimization Laboratory Title: Notes on duality in second order and -order cone otimization Author: Erling D. Andersen, Cornelis Roos and Tamás Terlaky AdvOl-Reort No. 000/8
More informationParabolic Flow in Parallel Plate Channel ME 412 Project 4
Parabolic Flow in Parallel Plate Channel ME 412 Project 4 Jingwei Zhu April 12, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of this project is to develop and apply a computer
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More information2.6 Primitive equations and vertical coordinates
Chater 2. The continuous equations 2.6 Primitive equations and vertical coordinates As Charney (1951) foresaw, most NWP modelers went back to using the rimitive equations, with the hydrostatic aroximation,
More informationAlgorithm for Solving Colocated Electromagnetic Fields Including Sources
Algorithm for Solving Colocated Electromagnetic Fields Including Sources Chris Aberle A thesis submitted in artial fulfillment of the requirements for the degree of Master of Science in Aeronautics & Astronautics
More informationHeat Transfer Analysis in the Cylinder of Reciprocating Compressor
Purdue University Purdue e-pubs International Comressor Engineering Conference School of Mechanical Engineering 2016 Heat Transfer Analysis in the Cylinder of Recirocating Comressor Ján Tuhovcák Brno University
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell October 25, 2009 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationDeriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.
Deriving ndicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deutsch Centre for Comutational Geostatistics Deartment of Civil &
More informationAn Ant Colony Optimization Approach to the Probabilistic Traveling Salesman Problem
An Ant Colony Otimization Aroach to the Probabilistic Traveling Salesman Problem Leonora Bianchi 1, Luca Maria Gambardella 1, and Marco Dorigo 2 1 IDSIA, Strada Cantonale Galleria 2, CH-6928 Manno, Switzerland
More informationFinding Shortest Hamiltonian Path is in P. Abstract
Finding Shortest Hamiltonian Path is in P Dhananay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune, India bstract The roblem of finding shortest Hamiltonian ath in a eighted comlete grah belongs
More informationSets of Real Numbers
Chater 4 Sets of Real Numbers 4. The Integers Z and their Proerties In our revious discussions about sets and functions the set of integers Z served as a key examle. Its ubiquitousness comes from the fact
More informationNumerical Study of Heat Transfer and Flow Bifurcation of CuO Nanofluid in Sudden Expansion Microchannel Using Two-Phase Model
Modern Mechanical Engineering, 2017, 7, 57-72 htt://www.scir.org/journal/mme ISSN Online: 2164-0181 ISSN Print: 2164-0165 Numerical Study of Heat Transfer and Flow Bifurcation of CuO Nanofluid in Sudden
More informationintegral invariant relations is not limited to one or two such
The Astronomical Journal, 126:3138 3142, 2003 December # 2003. The American Astronomical Society. All rights reserved. Printed in U.S.A. EFFICIENT ORBIT INTEGRATION BY SCALING AND ROTATION FOR CONSISTENCY
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationTHE APPLICATION OF DISCRETE TRANSFER METHOD IN THE SIMULATION OF RADIATION HEATING OR COOLING
, Volume, Number 1,.1-5, 2003 THE APPLICATION OF DISCRETE TRANSFER METHOD IN THE SIMULATION OF RADIATION HEATING OR COOLING Guangcai Gong, Kongqing Li and Gengxin Xie Civil Engineering College, Hunan University,
More informationChapter 6. Thermodynamics and the Equations of Motion
Chater 6 hermodynamics and the Equations of Motion 6.1 he first law of thermodynamics for a fluid and the equation of state. We noted in chater 4 that the full formulation of the equations of motion required
More informationA numerical implementation of a predictor-corrector algorithm for sufcient linear complementarity problem
A numerical imlementation of a redictor-corrector algorithm for sufcient linear comlementarity roblem BENTERKI DJAMEL University Ferhat Abbas of Setif-1 Faculty of science Laboratory of fundamental and
More informationCEE 452/652. Week 11, Lecture 2 Forces acting on particles, Cyclones. Dr. Dave DuBois Division of Atmospheric Sciences, Desert Research Institute
CEE 452/652 Week 11, Lecture 2 Forces acting on articles, Cyclones Dr. Dave DuBois Division of Atmosheric Sciences, Desert Research Institute Today s toics Today s toic: Forces acting on articles control
More informationA numerical assessment of the random walk particle tracking method for heterogeneous aquifers
288 Calibration and Reliability in Groundwater Modelling: From Uncertainty to Decision Making (Proceedings of ModelCARE 2005, The Hague, The Netherlands, June 2005). IAHS Publ. 304, 2006. A numerical assessment
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationRhie-Chow interpolation in OpenFOAM 1
Rhie-Chow interpolation in OpenFOAM 1 Fabian Peng Kärrholm Department of Applied Mechanics Chalmers Univesity of Technology Göteborg, Sweden, 2006 1 Appendix from Numerical Modelling of Diesel Spray Injection
More informationAnalysis of execution time for parallel algorithm to dertmine if it is worth the effort to code and debug in parallel
Performance Analysis Introduction Analysis of execution time for arallel algorithm to dertmine if it is worth the effort to code and debug in arallel Understanding barriers to high erformance and redict
More informationNumerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model
Numerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model Q. Sun, Y.J. Liu, R.H, Jin School of Aeronautics, Northwestern Polytechnical University, Xi an 710072,
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationImproving AOR Method for a Class of Two-by-Two Linear Systems
Alied Mathematics 2 2 236-24 doi:4236/am22226 Published Online February 2 (htt://scirporg/journal/am) Imroving AOR Method for a Class of To-by-To Linear Systems Abstract Cuixia Li Shiliang Wu 2 College
More informationMath 4400/6400 Homework #8 solutions. 1. Let P be an odd integer (not necessarily prime). Show that modulo 2,
MATH 4400 roblems. Math 4400/6400 Homework # solutions 1. Let P be an odd integer not necessarily rime. Show that modulo, { P 1 0 if P 1, 7 mod, 1 if P 3, mod. Proof. Suose that P 1 mod. Then we can write
More informationRigorous bounds on scaling laws in fluid dynamics
Rigorous bounds on scaling laws in fluid dynamics Felix Otto Notes written by Steffen Pottel Contents Main Result on Thermal Convection Derivation of the Model. A Stability Criterion........................
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More information