An efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy
|
|
- Abel Craig
- 5 years ago
- Views:
Transcription
1 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 and more suitable numerical scheme to solve efficiently the Allievi's equations governing the dynamics of fluids inside a pipe. The method arises from a proper combination of the Method Of Characteristics with a spline "shape preserving" interpolation and a nonlinear filtering procedure designed in order to eliminate the residual oscillations. 1 Introduction Problems of fluid-structure interaction or safety monitoring of penstocks need accurate and efficient numerical schemes for hydraulic transient computations. For the solution of these hyperbolic problems, the most common algorithm used is the well known "Method Of Characteristics" (MOC) because of its easy implementation as well as a low time consuming capability. Unfortunately with high gradients and discontinuities the accuracy is affected by significant overshootings. Among the enhancements proposed in literature, we recall "high resolution schemes" based on concepts like "upwinding" and "field by field" decomposition, [1]. These methods, in general, suffer from the disadvantage of a low computational efficiency and notable complexity, especially for practical applications. On the other hand, interpolation of the solution is required using MOC when grid and computational points does not match. Preliminary results, only based on splines "shape preserving" interpolation technique (SSP) have shown a lack of accuracy due to residual overshootings near discontinuities. Such spurious oscillations can be suppressed by means of a spatial nonlinear filtering procedure nested with MOC and SSP
2 388 Hydraulic Engineering Software 2 Numerical methods for the solution of the AllievFs equations The dynamics of a fluid in a straight circular pipe is described in time (t) and space (jc) by its velocity (V) and hydraulic head (H) fields, governed by the Allievi's equations: OH + v + = v dt dx g dx doc (i) g dt g dx dx 2Dg where c, D, z, /, g are the wave celerity, pipe diameter, vertical coordinate, friction coefficient, gravity acceleration respectively. In its simplest version MOC subdivides (usually uniformly) the space domain into intervals of amplitude Ac; the time step At is related to grid spacing Ac through the well known Courant condition: Af < A*Cf n\ W c + V When C,.= l, the positive characteristic lines intersect the time-space grid exactly in computational nodes, where V and H are already known from the previous step. If C,< 1 the need of fields V and H interpolation arises. Indeed this is the case of fluid-structure interaction problems (steel-water) where, for semplicity, spatial grid sizes are the same for both fluid and structural domains; consequently, time step amplitude is forced by the structural dynamics, resulting about 1/4 of thefluidone. For this purpose several schemes have been proposed in literature, leading to some accuracy problems. For example a linear interpolation scheme is a very dissipative one while quadratic interpolation schemes create strong overshootings. A very interesting approach has been proposed by Sibetheros et al. [2]. They proposed a spatial interpolation of V and H by means of spline functions with "shape preserving" capability; these special functions have the peculiarity of preserving the concavity of the data. Such an effect is obtained computing suitable additional grid points whenever very steep slopes in the unknown fields are encountered. A particular feature of this approach is that the Courant condition is no more a stability limit. In fact numerical tests have pointed out that solutions are stable and accurate also for At> Ax/(c+V). In order to check the validity of SSP method several benchmarks have been designed and faced. A comparison of the results has been done with a TVD method showing that the accuracy of SSP method is usually good and in particular that: the amplitude accuracy of SSP depends on the value of the ratio A/7 Ac; in particular is high for low values of the ratio and tends to have a local minimum
3 Hydraulic Engineering Software 389 for ratio values near unity; in this range residual overshooting can be observed, while TVD method is always accurate, the phase accuracy is very good over the whole range of At/Ax (from 0.2 to 1.2), while for TVD method is not, showing a considerable phase shift for low values (0.2). Consequently, some final considerations can be drawn: SSP is an interesting method, with a very good phase behaviour, allows the use of large time steps, shows high frequency residual oscillations. These drawbacks can be overcome by mean of a suitable numerical filtering algorithm. 3 Nonlinear Filtering Approach for the Suppression of Numerical Oscillations The idea of using nonlinear filtering procedures to eliminate or reduce numerically generated oscillations was early investigated by B. Engquist et al. [3] in More recently, Shyy et al. [4] (1992), evaluated the effectiveness of the filter carrying out a detailed study of the development of numerical oscillations generated by different dispersive and diffusive numerical methods, and their interactions with the filter action. Following these previous experiences, we decided to employ a similar nonlinear filtering approach to solve efficiently the Allievi's equations, (1). To check the general effectiveness of the filter, we analysed the application of several extended versions of the original scheme, such as TVD limiters or relaxation procedures, [4] Fundamentals The purpose of our approach is the coupling of a generic numerical integration scheme G( -) with a given filter operator P( ), [3]: where: v"+* represents the first approximated solution and u^^ the filtered one. For a filter P to be effective it should possess some basic properties to provide a minimal effect on an already smooth solution and to achieve an efficient action, with a minimal computational effort, to guarantee no spurious oscillations near discontinuities, [3]. Let us consider briefly the fundamental steps in the nonlinear filtering procedure, starting from the simpler version, with no TVD constraint. Given the vector solution: \ = (v^\ v^,..., v,/"^)\ from the first of eqns (3), we start with the following procedure: 1) search for the spatial components of the vector v, v/^\ j = l,...,n, with corrections of local extrema: (A+Vj"~^)(A-Vj"+*)<0, with: A+, A., the usual forward and backward differences respectively: A±VJ=±(VJ+I-VJ),
4 390 Hydraulic Engineering Software 2) the value corrections in step 1) are made by decreasing maxima and increasing minima, 3) when a correction is added to a grid pointy*, the same correction must be subtracted from a neighboring point in such a way to preserve a conservation criterion. More precisely, the corrected neighbor is the grid point j + 7*, or j- 7*, if the absolute maximum of differences: A+ and A coincides with: A+, or A., respectively, 4) no value may be corrected to overcome its neighbors. The filtering scheme at a given time step may be summarized as follow: 6 = min(6_, 6+/2) v. = v. + s 5 where the corrections are performed inside an iterative cycle in the spatial domain, and: 6* / I A W + l I I * W + l I _ = Engquist et al. showed that the above no TVD filtering scheme does not affects the convergence property of the original numerical scheme. The nonlinear filtering procedure has proved to be very effective in the suppression of spurious, numerically generated, oscillations of 2 Ax and 4Ax wavelengths, [4]. As the wavelengths become longer the filter effectiveness become smaller. The inability of the filter to treat long wavelength oscillations may results a pleasing feauture, because generally we aim to eliminate short wavelengths oscillations produced by the numerical scheme while only slightly to affect physically expected longer wavelength components. Regarding the procedure of practical filter implementation, since the initially generated oscillations tend to develop with time, the filtering action is most effective when applied at each time step of the computation. 3.2 Filtering improvements Many extensions of the above filter are in principle possible. In fact, the filtering procedure does not yield unique solutions with different degrees of dissipation and dispersion characteristics of a given numerical scheme. A relaxation procedure, for example, has been proposed by Shyy et al. [4] in order to enhance the effectiveness of the filtering. More precisely that improvement consists in the introduction, at the correction level, of a multiplicative parameter, w (relaxation parameter), mostly adjusted empirically. In [4] the authors showed that a w> 1, further improves the filtered solution. Indeed, the relaxation procedure introduces extra adjustments by removing the plateaus, i.e. the multipoints extrema:
5 Hydraulic Engineering Software = min(6_, 5J2) = v. + a) s 6 with: w>0. Another possible extension of the filtering consists in the generalization of the reduction coefficient for the quantity 5+, with numerical values different than 0.5. That version further enhances the effectiveness of the procedure, especially when residual plateaus are present: 6=min(<5_,a<5+) with w>0 and In order to eliminate some unwelcome properties of the basic filter, in the present work we consider also an important modification, first suggested by Engquist [3]. The aim is to modify the algorithm to obtain a TVD-enforcing filtering. Now the extremum condition at grid point JQ is replaced by the following one: min(v,\, va v,\) < v."""\ if. v"** is a local minimum, JQ * JQ JQ JQ JQ 6 max(vy"_i, Vy", v"+i) > v"*, if. v"* is a local maximum, for all extrema: 1^""^ of v"+\ Moreover for a TVD constraint a condition of no consecutive extrema is also necessary; that implies the replacement of the expression: 6=min(6_,a<5+) with the following one: where: 6 = min(6_, && i M+l, n n n \ d^ = Vj - max(v,-_i, Vj, v^ n n ^ n+l _i, v,, ^+J - v^ The final result is an accurate nonlinear filtering procedure with TVD condition. As already pointed out in [3,4], the above more sophisticated versions of the filtering procedure appear less robust than the first proposed one. Indeed at present does not exist general theorems proving that the TVD version is a stable finite algorithm, i.e. the filter terminates in a finite number of iterations on each time level, even if there are very small changes in the algorithm structure. Therefore it seems very opportune to carry out extensive theoretical and numerical analyses to clarify and to enlighten some important properties of the above appealing extensions of nonlinear filters.
6 392 Hydraulic Engineering Software 4 Applications Our test problem consists in a straight pipe (D=0.4 [m], lenght L=100 [m]) conveying water flux at V^ =.5 [m/s] from a reservoir to a valve suddenly closed at the beginning of the transient. Other parameters are: c=1450 [m/s], /=0.5, Ax=L/20 [m], C,=0.9. A parametric analysis has been carried out using different combinations of TVD and relaxation parameter, as well as different grid spacings. On the other hand only a representative result will be shown in the following. In Fig. 1 typical pressure time behaviours, obtained with SSP and SSPF (with relaxation, o> = 1.2, and TVD constraint), are reported. As anticipated SSP exhibits a strong residual overshootings near discontinuities, expecially at the beginning of the transient. Instead SSPF almost eliminates all the spurious oscillations without affecting the phase properties of the solution. Fig. 2 shows an enlargement of the transient start of Fig. 1: it is evident how residual oscillations are almost completely suppressed. Moreover the capability to make very sharp bends near all discontinuities is clearly pointed out. Filter action can be better enlighted through time evolution of the wave spatial profiles. In fact, Fig. 3 shows how the spurious oscillations arise just the transient starts and how they tend to amplify. On the other hand, in Fig. 4 spurious oscillations are completely absent. It is worth noting that this notable result has been reached applying the filtering technique each time step. 5 Conclusions The principal features of the proposed method make it suitable for all the wave propagation-like problems where a high level of accuracy and efficiency is required. The new scheme has been checked facing some significant test problems, mainly hydraulic shock wave propagation ones, obtaining efficiently, very accurate results. Some applications of the numerical technique to real systems, related to fluid-structure interaction and penstock monitoring cases, are in progress. References 1. Rider, W.J. Methods for extending high-resolution schemes to nonlinear systems of hyperbolic conservation laws, Int. Jour. Num. Meth. Fluids, 1993, Sibetheros, I.A., Hooley, E.R. and Branski, J.M. Spline interpolations for water hammer analysis Journal of Hydraulic Engineering, Vol. 117, No 10, Engquist, B.,L6tstedt, P. and Sjogreen, B. Nonlinear filters for efficient shock computation, Mathematics of Computation, 1989, 52, Shyy, W., Chen, M.H., Mittal, R. and Udaykumar, H.S. On the suppression of numerical oscillations using a nonlinear filter, Jour. Comp. Phys., 1992, 102.
7 Hydraulic Engineering Software 393 SRSF MetH->o<d ~WD o r-ne<ga = 1.2 Figure 1: Pressure time evolution at closed end computed by SSP (dashed line) and SSPF (continuous line). Pressure is expressed in Pa (N/nf) and time in sec. PSF Metnoa TVD ome<ga=1.2 Figure 2: Particular of Fig. 1: SSP (dashed line) vs. SSPF (continuous line).
8 394 Hydraulic Engineering Software SRS Method Space Profile Figure 3: Pressure wave spatial profiles at different time steps with SSP. x 1 O SRSF Method (TVO) Space Spaco Figure 4: Pressure wave spatial profiles at different time steps with SSPF (TVD)
Comparison of MOC and Lax FDE for simulating transients in Pipe Flows
International Research Journal of Engineering and Technology (IRJET) e-issn: 395-0056 Volume: 04 Issue: 03 Mar -07 www.irjet.net p-issn: 395-007 Comparison of MOC and Lax FDE for simulating transients
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 informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
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 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 informationSalmon: Lectures on partial differential equations
4 Burger s equation In Lecture 2 we remarked that if the coefficients in u x, y,! "! "x + v x,y,! "! "y = 0 depend not only on x,y but also on!, then the characteristics may cross and the solutions become
More informationComparison of (Some) Algorithms for Edge Gyrokinetics
Comparison of (Some) Algorithms for Edge Gyrokinetics Greg (G.W.) Hammett & Luc (J. L.) Peterson (PPPL) Gyrokinetic Turbulence Workshop, Wolfgang Pauli Institute, 15-19 Sep. 2008 w3.pppl.gov/~hammett Acknowledgments:
More informationEmpirical Models Interpolation Polynomial Models
Mathematical Modeling Lia Vas Empirical Models Interpolation Polynomial Models Lagrange Polynomial. Recall that two points (x 1, y 1 ) and (x 2, y 2 ) determine a unique line y = ax + b passing them (obtained
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
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 informationIMPORTANCE OF NUMERICAL EFFICIENCY FOR REAL TIME CONTROL OF TRANSIENT GRAVITY FLOWS IN SEWERS
1106 September 11~16, 2005, Seoul, Korea IMPORTANCE OF NUMERICAL EFFICIENCY FOR REAL TIME CONTROL OF TRANSIENT GRAVITY FLOWS IN SEWERS ARTURO S. LEON 1, MOHAMED S. GHIDAOUI 2, ARTHUR R. SCHMIDT 3 and MARCELO
More informationFlood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum
Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (29), No. 3 4, pp. 121 137 IBW PAN, ISSN 1231 3726 Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Dariusz
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 informationSolution Methods. Steady convection-diffusion equation. Lecture 05
Solution Methods Steady convection-diffusion equation Lecture 05 1 Navier-Stokes equation Suggested reading: Gauss divergence theorem Integral form The key step of the finite volume method is to integrate
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. #12 Fundamentals of Discretization: Finite Volume Method
More informationLecture 7: Edge Detection
#1 Lecture 7: Edge Detection Saad J Bedros sbedros@umn.edu Review From Last Lecture Definition of an Edge First Order Derivative Approximation as Edge Detector #2 This Lecture Examples of Edge Detection
More informationGENG2140, S2, 2012 Week 7: Curve fitting
GENG2140, S2, 2012 Week 7: Curve fitting Curve fitting is the process of constructing a curve, or mathematical function, f(x) that has the best fit to a series of data points Involves fitting lines and
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More informationIterative Controller Tuning Using Bode s Integrals
Iterative Controller Tuning Using Bode s Integrals A. Karimi, D. Garcia and R. Longchamp Laboratoire d automatique, École Polytechnique Fédérale de Lausanne (EPFL), 05 Lausanne, Switzerland. email: alireza.karimi@epfl.ch
More informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
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 informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationChapter 8: Flow in Pipes
8-1 Introduction 8-2 Laminar and Turbulent Flows 8-3 The Entrance Region 8-4 Laminar Flow in Pipes 8-5 Turbulent Flow in Pipes 8-6 Fully Developed Pipe Flow 8-7 Minor Losses 8-8 Piping Networks and Pump
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 informationLimits and Continuity. 2 lim. x x x 3. lim x. lim. sinq. 5. Find the horizontal asymptote (s) of. Summer Packet AP Calculus BC Page 4
Limits and Continuity t+ 1. lim t - t + 4. lim x x x x + - 9-18 x-. lim x 0 4-x- x 4. sinq lim - q q 5. Find the horizontal asymptote (s) of 7x-18 f ( x) = x+ 8 Summer Packet AP Calculus BC Page 4 6. x
More informationTutorial for the heated pipe with constant fluid properties in STAR-CCM+
Tutorial for the heated pipe with constant fluid properties in STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities
More informationFraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, Darmstadt, Germany
Scale Space and PDE methods in image analysis and processing Arjan Kuijper Fraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, 64283
More informationStudy of nonlinear phenomena in a tokamak plasma using a novel Hilbert transform technique
Study of nonlinear phenomena in a tokamak plasma using a novel Hilbert transform technique Daniel Raju, R. Jha and A. Sen Institute for Plasma Research, Bhat, Gandhinagar-382428, INDIA Abstract. A new
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 informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationSolution of Unsteady Flow Equations in High Pressure Pipe
Solution of Unsteady Flow Equations in High Pressure Pipe Bharati Medhi Das 1, Madan Mohan Das 2, Bibhash Sarma 3 Assistant Professor, Department of Civil Engineering, Assam Engineering College, Guwahati,
More informationME 305 Fluid Mechanics I. Part 8 Viscous Flow in Pipes and Ducts. Flow in Pipes and Ducts. Flow in Pipes and Ducts (cont d)
ME 305 Fluid Mechanics I Flow in Pipes and Ducts Flow in closed conduits (circular pipes and non-circular ducts) are very common. Part 8 Viscous Flow in Pipes and Ducts These presentations are prepared
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationNonlinear seismic imaging via reduced order model backprojection
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, Vladimir Druskin 2 and Mikhail Zaslavsky 2 University of Houston, 2 Schlumberger-Doll Research Center Mamonov, Druskin,
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationExtremum-Preserving Limiters for MUSCL and PPM
arxiv:0903.400v [physics.comp-ph] 7 Mar 009 Extremum-Preserving Limiters for MUSCL and PPM Michael Sekora Program in Applied and Computational Mathematics, Princeton University Princeton, NJ 08540, USA
More informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
More informationa( i), where N i N ) is accurately d( i), where N i N. The derivative at ( ) ( ) ( )
More on implementing the derivative filter John C. Bancroft ABSTRACT More on implementing the derivative filter An update is presented on finding faster and or more accurate implementations of differentiator
More informationME 305 Fluid Mechanics I. Chapter 8 Viscous Flow in Pipes and Ducts
ME 305 Fluid Mechanics I Chapter 8 Viscous Flow in Pipes and Ducts These presentations are prepared by Dr. Cüneyt Sert Department of Mechanical Engineering Middle East Technical University Ankara, Turkey
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 information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationLOSSES DUE TO PIPE FITTINGS
LOSSES DUE TO PIPE FITTINGS Aim: To determine the losses across the fittings in a pipe network Theory: The resistance to flow in a pipe network causes loss in the pressure head along the flow. The overall
More informationGradient and smoothing stencils with isotropic discretization error. Abstract
Gradient and smoothing stencils with isotropic discretization error Tommy Anderberg (Dated: April 30, 2012) Abstract Gradient stencils with isotropic O(h 2 ) and O(h 4 ) discretization error are constructed
More informationPIPING SYSTEMS. Pipe and Tubing Standards Sizes for pipes and tubes are standardized. Pipes are specified by a nominal diameter and a schedule number.
PIPING SYSTEMS In this chapter we will review some of the basic concepts associated with piping systems. Topics that will be considered in this chapter are - Pipe and tubing standards - Effective and hydraulic
More informationFormulation and performance of the Variable-Cubic Atmospheric Model
Formulation and performance of the Variable-Cubic Atmospheric Model John McGregor CSIRO Marine and Atmospheric Research Aspendale, Melbourne Southern Hemisphere PDEs on the Sphere NCAR 11 April 2014 CSIRO
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
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 informationCurve Fitting. Objectives
Curve Fitting Objectives Understanding the difference between regression and interpolation. Knowing how to fit curve of discrete with least-squares regression. Knowing how to compute and understand the
More informationTutorial for the supercritical pressure pipe with STAR-CCM+
Tutorial for the supercritical pressure pipe with STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with
More informationChapter 8: Flow in Pipes
Objectives 1. Have a deeper understanding of laminar and turbulent flow in pipes and the analysis of fully developed flow 2. Calculate the major and minor losses associated with pipe flow in piping networks
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 informationLesson 6 Review of fundamentals: Fluid flow
Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass
More informationChapter 9. Non-Parametric Density Function Estimation
9-1 Density Estimation Version 1.1 Chapter 9 Non-Parametric Density Function Estimation 9.1. Introduction We have discussed several estimation techniques: method of moments, maximum likelihood, and least
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
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 informationStructural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake
Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Seung Keun Park and Hae Sung Lee ABSTRACT This paper presents a system identification (SI) scheme
More informationFluid Mechanics II 3 credit hour. Fluid flow through pipes-minor losses
COURSE NUMBER: ME 323 Fluid Mechanics II 3 credit hour Fluid flow through pipes-minor losses Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Losses in Noncircular
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationOn the diminishing of spurious oscillations in explicit finite element analysis of linear and non-linear wave propagation and contact problems
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16315 On the diminishing of spurious oscillations
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 informationCurve Fitting. 1 Interpolation. 2 Composite Fitting. 1.1 Fitting f(x) 1.2 Hermite interpolation. 2.1 Parabolic and Cubic Splines
Curve Fitting Why do we want to curve fit? In general, we fit data points to produce a smooth representation of the system whose response generated the data points We do this for a variety of reasons 1
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 informationThe Principle of Least Action
The Principle of Least Action In their never-ending search for general principles, from which various laws of Physics could be derived, physicists, and most notably theoretical physicists, have often made
More information4. Multilayer Perceptrons
4. Multilayer Perceptrons This is a supervised error-correction learning algorithm. 1 4.1 Introduction A multilayer feedforward network consists of an input layer, one or more hidden layers, and an output
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
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 informationINTERPOLATION. and y i = cos x i, i = 0, 1, 2 This gives us the three points. Now find a quadratic polynomial. p(x) = a 0 + a 1 x + a 2 x 2.
INTERPOLATION Interpolation is a process of finding a formula (often a polynomial) whose graph will pass through a given set of points (x, y). As an example, consider defining and x 0 = 0, x 1 = π/4, x
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 informationApplication of the relaxat ion met hod to model hydraulic jumps
Application of the relaxat ion met hod to model hydraulic jumps P. J. Montgomery Mathematics and Computer Science Program, University of Northern British Columbia, Prince George, Canada. Abstract A finite
More informationProcessing Well Log Data Using Empirical Mode Decomposition (EMD) and Diffusion Equation*
Processing Well Log Data Using Empirical Mode Decomposition (EMD) and Diffusion Equation* Sangmyung D. Kim 1, Seok Jung Kwon 1, and Kwang Min Jeon 1 Search and Discovery Article #120188 (2015) Posted April
More informationDiscretization of Convection Diffusion type equation
Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
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 informationFrancis X. Giraldo,
1 Time-Integrators Francis X. Giraldo, giraldo@nrlmry.navy.mil, www.nrlmry.navy.mil/~giraldo/projects/nseam.html 1.1 Introduction Roughly speaking, there are 2 classes of TIs: 1. EulerianMethods(fixed-frame-e.g.,arockatthebottomofaflowing
More informationPipe Flow. Lecture 17
Pipe Flow Lecture 7 Pipe Flow and the Energy Equation For pipe flow, the Bernoulli equation alone is not sufficient. Friction loss along the pipe, and momentum loss through diameter changes and corners
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationComparison of Averaging Methods for Interface Conductivities in One-dimensional Unsaturated Flow in Layered Soils
Comparison of Averaging Methods for Interface Conductivities in One-dimensional Unsaturated Flow in Layered Soils Ruowen Liu, Bruno Welfert and Sandra Houston School of Mathematical & Statistical Sciences,
More informationLecture Hilbert-Huang Transform. An examination of Fourier Analysis. Existing non-stationary data handling method
Lecture 12-13 Hilbert-Huang Transform Background: An examination of Fourier Analysis Existing non-stationary data handling method Instantaneous frequency Intrinsic mode functions(imf) Empirical mode decomposition(emd)
More informationHilbert-Huang and Morlet wavelet transformation
Hilbert-Huang and Morlet wavelet transformation Sonny Lion (sonny.lion@obspm.fr) LESIA, Observatoire de Paris The Hilbert-Huang Transform The main objective of this talk is to serve as a guide for understanding,
More informationToward Binary Black Hole Simulations in Numerical Relativity
Toward Binary Black Hole Simulations in Numerical Relativity Frans Pretorius California Institute of Technology BIRS Workshop on Numerical Relativity Banff, April 19 2005 Outline generalized harmonic coordinates
More informationEFFECT 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 informationLinear Diffusion. E9 242 STIP- R. Venkatesh Babu IISc
Linear Diffusion Derivation of Heat equation Consider a 2D hot plate with Initial temperature profile I 0 (x, y) Uniform (isotropic) conduction coefficient c Unit thickness (along z) Problem: What is temperature
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationLecture 12: Transcritical flow over an obstacle
Lecture 12: Transcritical flow over an obstacle Lecturer: Roger Grimshaw. Write-up: Erinna Chen June 22, 2009 1 Introduction The flow of a fluid over an obstacle is a classical and fundamental problem
More informationNonlinear Solitary Wave Basis Functions for Long Time Solution of Wave Equation
Nonlinear Solitary Wave Basis Functions for Long Time Solution of Wave Equation John Steinhoff Subhashini Chitta Nonlinear Waves --Theory and Applications Beijing - June, 2008 Objective Simulate short
More informationVarberg 8e-9e-ET Version Table of Contents Comparisons
Varberg 8e-9e-ET Version Table of Contents Comparisons 8th Edition 9th Edition Early Transcendentals 9 Ch Sec Title Ch Sec Title Ch Sec Title 1 PRELIMINARIES 0 PRELIMINARIES 0 PRELIMINARIES 1.1 The Real
More informationSeptember Math Course: First Order Derivative
September Math Course: First Order Derivative Arina Nikandrova Functions Function y = f (x), where x is either be a scalar or a vector of several variables (x,..., x n ), can be thought of as a rule which
More information