Hydraulic Modelling for Drilling Automation
|
|
- Matthew Shelton
- 6 years ago
- Views:
Transcription
1 Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts
2 Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators: 1. Norway : G. O. Kaasa (Kelda Drilling Controls, NTNU) 2. France : F. de Meglio (Ecole de Mines) PhD Group : M.H. Abbasi (affiliated to CASA and Kelda Drilling Controls) Project Sponsors The project HYMODRA (HYdraulic MOdelling for DRilling Automation) is sponsored by Shell and NWO-I under the aegis of Shell NWO/FOM PhD Programme in Computational Sciences for Energy Research. 2 Supervisors and Collaborators
3 Application Perspective Mathematical Notion of Drift Flux Model Physical and Numerical Boundary Conditions Sound Speed Model Numerical Schemes Numerical results Conclusions and perspectives 3 Outline
4 Main Goal : Develop hydraulic models and supporting model reduction techniques that are 1. accurate enough 2. simple enough to be employed in the context of drilling scenario simulations and real time estimation and control in case of gas influx. Managed Pressure Drilling!! Figure: Drilling Schematic 4 Objective
5 Characteristics/Features: The downhole pressure must be kept within allowable limits Delays in transmission of information No downhole measurements during certain phases of drilling operations Top-side measurements are available Figure: Managed Pressure Drilling Hydraulic Model: Serve as a model for controller/ estimator Aid to design operations before hand 5 Characteristics/Features
6 Complexities from physical perspective Timescales in the drilling process Slow transient corresponding to mass transport Fast transient corresponding to the propagation of the acoustic waves Nonlinearities: Acoustic velocity changes very rapidly in the one-phase to two-phase transition regions and vice versa. Disappearance and Appearance of Phases. Various flow regimes across different sections of the well. Distributed non-linearities due to source terms. 6 Complexities from physical perspective
7 Governing Equations of Drift Flux Model t (ρ l α l ) + x (ρ l α l v l ) = Γ l t (ρ g α g ) + x (ρ g α g v g ) = Γ g t (ρ l α l v l + ρ g α g v g ) + x (ρ l α l v 2 l + ρ g α g v 2 g + P) = Q g + Q v Q g = g(ρ l α l + ρ g α g )sin(θ); Q v = 32µv d 2 α l = Liquid Void Fraction ; α g = Gas Void Fraction ρ l = Liquid Density ; ρ g = Gas Density v l = Liquid Velocity ; v g = Gas Velocity ; v = Mixture Flow Velocity Γ l and Γ g are the phase change terms P = Pressure ; d = hydraulic diameter ; µ = mixture viscosity θ is the well inclination ; g = acceleration due to gravity 7 Governing Equations of Drift Flux Model
8 Closure Laws α g + α l = 1 ρ g = P /a 2 g ρ l = ρ l0 + (P P l0 )/a 2 l v g = (Kv l α l + S)/(1 K α g ) K and S are flow dependent parameters. There is singularity in the slip law when we approach pure gas region. a l is the speed of sound in the liquid phase a g is the speed of sound in the liquid phase P l0 : standard atmospheric pressure ρ l0 : density of liquid at standard atmospheric pressure 8 Closure Laws Modelling
9 The 1-D non-linear conservation law: 9 System of Conservation Laws w t + (f (w)) x = s is hyperbolic if the Jacobian matrix f w is diagonalizable with real eigenvalues for each physically relevant w. ρ l α l ρ l α l v l w = ρ g α g ρ l α l v l + ρ g α g v g, f (w) = ρ g α g v g ρ l α l v 2 l + ρ g α g v 2 g + P Γ l s = Γ g Q g + Q v For further discussion, Γ l = 0, Γ g = 0
10 Eigenvalues of the Jacobian Matrix The corresponding eigenvalues are given by: λ 1 = v l ω, λ 2 = v l + ω, λ 3 = v g where, ω is the speed of sound in two-phase mixture Two eigenvalues are linked to the compressibility effects Third eigenvalue is coincident to the gas velocity One pressure pulse propagates downstream and the other pressure pulse propagates upstream. Gas Volume wave travels downstream. 10 Eigenvalues of the Jacobian Matrix
11 Physical Boundary Conditions: (ρ l α l v l )(0, t) = f (t) (ρ g α g v g )(0, t) = h (t) P(L, t) = r(t) or (α l v l )(0, t) = f (t) (α g v g )(0, t) = h (t) P(L, t) = r(t) 11 Physical Boundary Conditions
12 Compatibility relations for multi phase system, which are: Characteristic 1 and 2: Compatibility relation corresponding to the pressure wave propagating in the upstream direction and downstream direction of the flow: d dt p + ρ l ω(v g v l ) d dt α g ρ l α l (v g v l + ω) d dt v l = q(v g v l + ω) where, d dt = t + (v l + ω) x is the directional derivative Characteristic 3: Compatibility relation corresponding to the gas volume wave: d dt p + p d α g (1 K α g ) dt α g = 0 where, d dt = t + (v g) x is the directional derivative 12 Numerical Boundary Conditions
13 Sound Speed Model Approximate Sound Speed Model is written as: a l if α g < ɛ ω = c(p, α g, ρ l, K) if ɛ α g 1 ɛ if α g > 1 ɛ a g where ɛ is a small parameter P c(p, α g, ρ l, K) = α g ρ l (1 K α g ) a l and a g are the sound speeds in liquid and gas medium respectively 13 Sound Speed Model
14 Assumptions for Sound Speed in the two phase mixture Liquid is incompressible αg ρ g << α l ρ l Why is sound speed in the two phase mixture important? Numerical flux computations are heavily dependent on the mixture sound speed Numerical dissipation depends on the sound speed of the two phase mixture Enable correct determination of locations and speeds of the wave fronts 14 Sound Speed Model
15 Reasons for Model Improvement Drilling fluids are highly compressible Existing models for sound speed in two phase mixture are singular at low and high void fractions Existing models become singular before rendering the Drift Flux Model non-hyperbolic Existing models also fail in modelling the realistic effects at high operating pressures Need of a unified model for single phase flow and two phase flow modelling 15 Sound Speed Model
16 ρ l ρ C = g (ρ g ( Kv l (1 αg )+S 1 K αg ) ρ l v l ) ( (1 αg )v l a 2 l (1 αg )v ρ l v l l a l 2 D = ρ g ( Kv l (1 αg )+S αg ( Kv l (1 αg )+S 1 K ) αg 1 K ) αg ag 2 ( v l 2ρ l + ( Kv l (1 αg )+S 1 K ) 2 ρ αg g ) (1 + ( Kv l (1 αg )+S 2 1 K ) αg αg ag 2 16 Modified Sound Speed Model 1 αg a l 2 0 αg ag αg ( Kv l (1 αg )+S 1 K ) αg ag 2 ) ((1 α g )ρ l + α g ρ g ( K(1 αg ) 1 K αg )) ρ l (1 α g ) α g ρ g ( K(1 αg ) 1 K αg ) + v2 l (1 αg ) a l 2 ) (2ρ l (1 α g )v l + 2ρ g α g v g ( K(1 αg ) 1 K αg )) Eigenvalues of the Jacobian matrix C 1 D can be computed numerically. In particular for, K=1 and S=0 i.e. assuming zero slip between the liquid and gaseous phase. Modified sound speed comes out to be: 1/2 ω new (ρ g ρ = a g a l ( l ) ((ρ l + α g ρ g α g ρ l )(ag 2 ρ g α g ag 2 ρ g + α g al 2ρ l )) )
17 Figure: Comparative plot Figure: Zooming the comparative plot 17
18 W n+1 i = W n i w t + (f (w)) x = s t { } F n (W x i + 1 L, W R ) F n (W 2 i 2 1 L, W R ) + ts n i W L and W R are estimated value of variables at left and right cell interface respectively 18 Full Discretization Figure: Stencil for discretization in space and time
19 Features of Hyperbolic PDE: Information propagates with finite speed and has preferred direction Discontinuities or shock waves develop in a finite time and propagate even if initial and boundary data are smooth. Requirements from the Numerical Method: Sharp Resolution of discontinuities No spurious oscillations Minimal smearing effect Consistent, Stable and Convergent Conservation property in discrete sense 19 Numerical Methods
20 Approximation of Numerical Flux Liquid Gas {}}{{}}{ Fi FVS +1/2 (w L, w R ) = (α l ρ l ) L Ψ + l,l + (α l ρ l ) R Ψ l,r + (α g ρ g ) L Ψ + g,l + (αg ρg ) R Ψ g,r } {{ } Numerical Convective Flux + (F p ) i +1/2 }{{} Numerical Pressure Flux Liquid Contribution Gas Contribution Pressure Contribution Ψ + l,l = Ψ+ (v l l,l, ω i +1/2 ) Ψ l,r = Ψ (v l l,r, ω i +1/2 ) Ψ + (v, ω) = V + (v, ω) 1 0 l v Ψ (v, ω) = V (v, ω) 1 0 l v Ψ + g,l = Ψ+ g (v g,l, ω i +1/2 ) Ψ g,r = Ψ g (v g,r, ω i +1/2 ) Ψ g + (v, ω) = V + (v, ω) 0 1 v Ψg (v, ω) = V (v, ω) 0 1 v (F p ) i +1/2 = ( 0 0 p i +1/2 ) T p i +1/2 = P + (v L, ω i +1/2 )p L + P (v R, ω i +1/2 )p R v = mixture fluid velocity Splitting Functions V ± and P ± are the functions that satisfy the consistency, upwinding, monotonicity, differentiability and positivity property 20 Approximation of Numerical Flux
21 Numerical Test Cases No analytical results exist for Drift Flux Model. We try out numerical benchmark tests for multiphase flow problems. 1. Shock Tube:Shock capturing due to pressure difference 2. Fast Transients: Propagation of pressure pulses. 3. Slow Transients: Propagation of mass transport wave Correct description of fluid transport and pressure waves requires high resolution schemes possessing little numerical diffusion. Both first order and second order schemes were investigated. 21 Numerical Test Cases
22 Figure: Wave Fronts 22 Numerical Benchmarking Figure: Numerical Example
23 23 Shock Tube Figure: Behaviour of Gas Void Fraction using first order FVS
24 24 Shock Tube Figure: Behaviour of Pressure using first order FVS
25 25 Shock Tube Figure: Behaviour of Liquid Velocity using first order FVS
26 Fast Transients : Test Case 1 Figure: Fast Transient Test Case Capturing fast transients allows the modelling of water hammer effects. 26 Fast Transients
27 Figure: Snapshots of fast transients test case using first order FVS at CFL = 0.25; Gas Volume Fraction(left), Liquid Velocity(middle), Pressure(right) 27 Numerical Results
28 Fast Transients : Test Case 2 Figure: Fast Transients Test Case Capturing fast transients allows the modelling of water hammer effects. 28 Fast Transients
29 Figure: Comparison between second order AUSM scheme and second order FVS scheme at CFL = 0.25; Pressure(left) and Liquid Velocity(right) 29 Numerical Results
30 30 Numerical Results Figure: Simulation of Fast Transient using AUSM scheme
31 Slow Transients : Test Case 1 Figure: Slow Transients Test Case Models transient behaviour induced by injecting gas and liquid at the inlet 31 Slow Transients
32 Future Work Non-linear Stability Analysis For hyperbolic conservation laws, the spectrum of the upwind spatial differential operator constitutes eigenvalues that lie in the left half plane near the imaginary axis The absolute stability region of the forward Euler method intersects the imaginary axis only at the origin Forward Euler is typically not a stable choice of time discretization; furthermore it is only first order accurate Nonlinear stability conditions become critical for the convergence in the presence of shocks or sharp gradients Establish order of merit of the numerical scheme 32 Future Work
33 Future Work Model Order Reduction Based on the properties of the fully discretized or semi discretized models, an appropriate model order reduction technique needs to be obtained, which: Handle non-linearities and delays (due to wave propagation) Preserves stability characteristics of the original model Preserves multiple time scales involved in the problem Preserves input-output behaviour of the original system 33 Future Work
34 Future Work Numerical Modelling of Tripping Benchmark Scenario Challenges from Simulation Perspective Cross sectional area changes dynamically as the pipe moves Regridding of the annular region Higher than 1D model would be more accurate Figure: Drilling Process 34 Future Work
35 35 Thank You for your attention!!
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationA WIMF Scheme for the Drift-Flux Two-Phase Flow Model
A WIMF Scheme for the Drift-Flux Two-Phase Flow Model Steinar Eve A,B,D, Tore Flåtten B,A and Svend Tollak Munkeord C A Centre of Mathematics for Applications (CMA), 153 Blindern, NO-316 Oslo, Norway B
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationFlux corrected transport solver for solving 1D multiphase equations for drilling applications
Computational Methods in Multiphase Flow VII 169 Flu corrected transport solver for solving 1D multiphase equations for drilling applications P. A. Varadarajan & P. S. Hammond Schlumberger Gould Research,
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
More informationNodalization. The student should be able to develop, with justification, a node-link diagram given a thermalhydraulic system.
Nodalization 3-1 Chapter 3 Nodalization 3.1 Introduction 3.1.1 Chapter content This chapter focusses on establishing a rationale for, and the setting up of, the geometric representation of thermalhydraulic
More informationTHE numerical simulation of the creation and evolution
Proceedings of the World Congress on Engineering Vol III WCE, July 4-6,, London, U.K. Numerical Simulation of Compressible Two-phase Flows Using an Eulerian Type Reduced Model A. Ballil, Member, IAENG,
More informationRECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS
RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS Dean Wang April 30, 2015 24.505 Nuclear Reactor Physics Outline 2 Introduction and Background Coupled T-H/Neutronics Safety Analysis Numerical schemes
More informationComputational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid
Computational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid Prasanna Welahettige 1, Bernt Lie 1, Knut Vaagsaether 1 1 Department of Process, Energy and Environmental
More informationAn efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy
An efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy Abstract The aim of this paper is to present a new
More informationThe multi-stage centred-scheme approach applied to a drift-flux two-phase flow model
The multi-stage centred-scheme approach applied to a drift-flux two-phase flow model Svend Tollak Munkejord 1 Steinar Evje 2 Tore Flåtten 2 1 Norwegian University of Science and Technology (NTNU), Department
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationMass flow determination in flashing openings
Int. Jnl. of Multiphysics Volume 3 Number 4 009 40 Mass flow determination in flashing openings Geanette Polanco Universidad Simón Bolívar Arne Holdø Narvik University College George Munday Coventry University
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationTCCS Published by Elsevier Ltd. Selection and/or peer-review under responsibility of [name organizer]
Energy Procedia 00 (2011) 000 000 Energy Procedia www.elsevier.com/locate/procedia TCCS-6 Depressurization of CO 2 a numerical benchmark study Sigmund Clausen a,c* and Svend Tollak Munkejord b a Gassco
More informationA Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems
A Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems Dia Zeidan,a) and Eric Goncalves 2 School of Basic Sciences and Humanities, German Jordanian University, Amman, Jordan
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationA General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations
A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations Ronald P. Fedkiw Xu-Dong Liu Stanley Osher September, 2 Abstract Standard
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS 1D-time averaged models
A SHORT INTRODUCTION TO TWO-PHASE FLOWS 1D-time averaged models Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +3304 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationInvestigation of slug flow characteristics in inclined pipelines
Computational Methods in Multiphase Flow IV 185 Investigation of slug flow characteristics in inclined pipelines J. N. E. Carneiro & A. O. Nieckele Department of Mechanical Engineering Pontifícia Universidade
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationNumerical investigation of cavitation-regimes in a converging-diverging nozzle
Numerical investigation of cavitation-regimes in a converging-diverging nozzle 1 Polina Gorkh, 1 Steffen J. Schmidt, and 1 Nikolaus A. Adams 1 Institute of Aerodynamics and Fluid Mechanics, Department
More informationSlug tracking simulation of severe slugging experiments
Slug tracking simulation of severe slugging experiments Tor Kindsbekken Kjeldby, Ruud Henkes and Ole Jørgen Nydal Abstract Experimental data from an atmospheric air/water terrain slugging case has been
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationNumerical Schemes for 1-D Two-Phase Flows
Numerical Schemes for 1-D Two-Phase Flows by Qiulan Zeng A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Doctor of Philosophy in Mechanical and
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationTransient Phenomena in Liquid/Gas Flow in Pipelines
Proceedings of the International Conference on Heat Transfer and Fluid Flow Prague, Czech Republic, August 11-12, 214 Paper No. 71 Transient Phenomena in Liquid/Gas Flow in Pipelines Zohra Ouchiha, S.
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More informationTarget Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu
Muon Collider/Neutrino Factory Collaboration Meeting May 26 28, CERN, Geneva U.S. Department of Energy Target Simulations Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu Center for Data Intensive
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
More informationCapability of CFD-tools for supersonic multiphase flows
OpenFOAM Capability of CFD-tools for supersonic multiphase flows by József Nagy and Michael Harasek Simulations for both multiphase flows and supersonic single phased flows are well known, however the
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
More informationComparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows
Comparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows Dr Argiris I. Delis Dr Ioannis K. Nikolos (TUC) Maria
More informationFluid Dynamics. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/14
Fluid Dynamics p.1/14 Fluid Dynamics Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/14 The equations of fluid dynamics are coupled PDEs that form an IVP (hyperbolic). Use
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationEulerian interface-sharpening methods for hyperbolic problems
Eulerian interface-sharpening methods for hyperbolic problems Application to compressible multiphase flow Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan 11:55-12:25, March 05,
More informationCranfield ^91. College of Aeronautics Report No.9007 March The Dry-Bed Problem in Shallow-Water Flows. E F Toro
Cranfield ^91 College of Aeronautics Report No.9007 March 1990 The Dry-Bed Problem in Shallow-Water Flows E F Toro College of Aeronautics Cranfield Institute of Technology Cranfield. Bedford MK43 OAL.
More informationGodunov methods in GANDALF
Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
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 informationSECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS. Min Tang. (Communicated by the associate editor name)
Manuscript submitted to AIMS Journals Volume X, Number X, XX 2X Website: http://aimsciences.org pp. X XX SECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS Min Tang Department of Mathematics
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More information1 Energy dissipation in astrophysical plasmas
1 1 Energy dissipation in astrophysical plasmas The following presentation should give a summary of possible mechanisms, that can give rise to temperatures in astrophysical plasmas. It will be classified
More informationA phase field model for the coupling between Navier-Stokes and e
A phase field model for the coupling between Navier-Stokes and electrokinetic equations Instituto de Matemáticas, CSIC Collaborators: C. Eck, G. Grün, F. Klingbeil (Erlangen Univertsität), O. Vantzos (Bonn)
More informationTwo-Fluid Model 41. Simple isothermal two-fluid two-phase models for stratified flow:
Two-Fluid Model 41 If I have seen further it is by standing on the shoulders of giants. Isaac Newton, 1675 3 Two-Fluid Model Simple isothermal two-fluid two-phase models for stratified flow: Mass and momentum
More informationNumerically Solving Partial Differential Equations
Numerically Solving Partial Differential Equations Michael Lavell Department of Applied Mathematics and Statistics Abstract The physics describing the fundamental principles of fluid dynamics can be written
More informationSevere slugging: modeling, simulation and stability criteria
Severe slugging: modeling, simulation and stability criteria Jorge Luis Baliño jlbalino@usp.br Departamento de Engenharia Mecânica Escola Politécnica Outline Introduction Air-water model Simulation results
More informationThe Hopf equation. The Hopf equation A toy model of fluid mechanics
The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationHead loss coefficient through sharp-edged orifices
Head loss coefficient through sharp-edged orifices Nicolas J. Adam, Giovanni De Cesare and Anton J. Schleiss Laboratory of Hydraulic Constructions, Ecole Polytechnique fédérale de Lausanne, Lausanne, Switzerland
More informationSound Waves Sound Waves:
3//18 Sound Waves Sound Waves: 1 3//18 Sound Waves Linear Waves compression rarefaction Inference of Sound Wave equation: Sound Waves We look at small disturbances in a compressible medium (note: compressible
More informationAeroacoustic and Aerodynamics of Swirling Flows*
Aeroacoustic and Aerodynamics of Swirling Flows* Hafiz M. Atassi University of Notre Dame * supported by ONR grant and OAIAC OVERVIEW OF PRESENTATION Disturbances in Swirling Flows Normal Mode Analysis
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationExperimental and Numerical Investigation of Two- Phase Flow through Enlarging Singularity
Purdue University Purdue e-pubs International Refrigeration and Air Conditioning Conference School of Mechanical Engineering 212 Experimental and Numerical Investigation of Two- Phase Flow through Enlarging
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationFluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows
Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 17 Laminar and Turbulent flows Welcome back to the video course on fluid mechanics. In
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
More informationNUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00786 NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A
More informationFigure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m
1. For the manometer shown in figure 1, if the absolute pressure at point A is 1.013 10 5 Pa, the absolute pressure at point B is (ρ water =10 3 kg/m 3, ρ Hg =13.56 10 3 kg/m 3, ρ oil = 800kg/m 3 ): (a)
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationAN EFFICIENT NUMERICAL SCHEME FOR MODELING TWO-PHASE BUBBLY HOMOGENEOUS AIR-WATER MIXTURES
AN EFFICIENT NUMERICAL SCHEME FOR MODELING TWO-PHASE BUBBLY HOMOGENEOUS AIR-WATER MIXTURES ARTURO S. LEÓN 1, MOHAMED S. GHIDAOUI 2, ARTHUR R. SCHMIDT 3 and MARCELO H. GARCÍA 4 1 Ph.D. Candidate, Dept.
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationNumerical and mathematical analysis of a five-equation model for two-phase flow
SINTEF Energy Research Master s Thesis PFE-Master Numerical and mathematical analysis of a five-equation model for two-phase flow Author Pedro José Martínez Ferrer Supervisors Svend Tollak Munkejord Tore
More informationA FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS
A FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS Ruben Juanes 1 and Knut-Andreas Lie 2 1 Stanford University, Dept. Petroleum Engineering, USA 2 SINTEF IKT, Dept., Norway ECMOR IX, August 30 September
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 34 Outline 1 Lecture 7: Recall on Thermodynamics
More informationNUMERICAL METHOD FOR THREE DIMENSIONAL STEADY-STATE TWO-PHASE FLOW CALCULATIONS
' ( '- /A NUMERCAL METHOD FOR THREE DMENSONAL STEADY-STATE TWO-PHASE FLOW CALCULATONS P. Raymond,. Toumi (CEA) CE-Saclay DMT/SERMA 91191 Gif sur Yvette, FRANCE Tel: (331) 69 08 26 21 / Fax: 69 08 23 81
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationInternal boundary layers in the ocean circulation
Internal boundary layers in the ocean circulation Lecture 9 by Andrew Wells We have so far considered boundary layers adjacent to physical boundaries. However, it is also possible to find boundary layers
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationModel adaptation in hierarchies of hyperbolic systems
Model adaptation in hierarchies of hyperbolic systems Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France February 15th, 2012 DFG-CNRS Workshop Nicolas Seguin (LJLL, UPMC) 1 / 29 Outline of the
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More informationAdvection, Conservation, Conserved Physical Quantities, Wave Equations
EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University Contents EP711 Supplementary
More informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
More information