IMPLICIT FORMULATION FOR 1D AND 2D ST VENANT EQUATIONS PRESENTATION OF THE METHOD, VALIDATION AND APPLICATIONS
|
|
- Anis Dixon
- 5 years ago
- Views:
Transcription
1 IMPLICIT FORMULATION FOR 1D AND 2D ST VENANT EQUATIONS PRESENTATION OF THE METHOD, VALIDATION AND APPLICATIONS Tierry Lepelletier 1, Quentin Araud 2 1 Hydratec Immeuble Central Seine - 42 quai de la Râpée Paris - lepelletier@ ydra.setec.fr 2 Hydratec 1 rue de la Course Strasbourg araud@ydra.setec.fr KEY WORDS St Venant equations, Implicit formulation, solving procedure, sock capture, bencmarking. ABSTRACT St Venant equations form te teoretical basis to ydraulic modeling in te field of river flow and urban ydraulics. Numerical scemes may prove difficult to implement due to te yperbolic nature of tese equations and to te variety of configurations met in practice. Existing ydro-informatic software are usually specialized in a given domain. Optimized solving scemes are selected to best fit te practical performance wic is required in tis domain of application. In order to meet te needs to a large variety of domains an original formulation as been developed to solve te full St Venant equations in one and two dimensions, wit optional simplifications to optimize some computations, wile keeping te same code arcitecture. Te underlying formulation is based upon an original metod for solving full St Venant equations troug finite volume space discretization and implicit time marcing solving procedure. Tis metod is unconditionally stable, time step is variable and is adjusted automatically witin calculation in order to preserve numerical accuracy in case of discontinuities, suc as sock formation. Tis numerical sceme is implemented in Hydra software, wic is developed by HYDRATEC. Various 1D and 2D tests, dam break problems and recirculation in 2D domains are presented and are compared wit oter referenced software in teir respective domain of application, suc as Telemac (developed by EDF). 1 DISCRETIZATION OF TWO DIMENSIONAL ST VENANT EQUATIONS AND SOLVING PROCEDURE St Venant equations are first expressed in classical integral form witin a cell : UU Ωii ddωii + (FF Γ ii nn xx + GGnn yy )ddγ ii = SS ff Ωii ddωii + ggrad(z b ) dωi wit : (1) Ωi UU = qq xx uu xx qq yy uu yy qq xx qq xx ² FF = + gg ² 2 GG = qq xx qq yy qq yy ² qq yy qq xx qq yy + gg ² 2 n Notations are defined in appended glossary. Pressure and gravity terms are combined and are approximated as follows : gg ( Γ ²)nnddΓ ii 2 ii + ggrad(z b )dωi gg ıııı Ω jj zz ii 2 jj zz ii nnddγ Γ iiii were ıııı = ii+ jj iiii 2 (2)
2 Tis expression is valid as long as differences between water depts between two adjacent cells remain small. It turns out tat it remains valid in presence of socks as will be sown in section 2 below. zi i zj j Mes i Mes j Spatial discretization of equation (1) consists in decomposing te wole domain into quadrangular or triangular cells and to average components of te unknown U vector witin eac cell. Connectivity condition between cells must be ensured. Equation (1) is rewritten as follows in discretized form, using te approximation of equation (2) : AA ii UU ii FF = + AA iiss ffff + jj FFnn xx + GGnn yy ll iiii = 0 (3) were : qq xx ² qq xx + 0.5gg zz ıııı jj zz ii GG = qq xx qq yy qq yy ² qq yy qq xx qq yy + 0.5gg zz ıııı jj zz ii Te terms jj FFnn xx + GGnn yy ll iiii Flux ij represent te flux terms across two adjacent cells. Tis alternative formulation yields following advantages : - gravity terms are source terms but tey are included in te line integral terms F and G wic are expressed along eac cell boundary : expression (3) is best fitted to link togeter equations between different domains, suc as 2D and 1D domain. Linkage is provided by te flux terms Flux ij wic ave different expressions according to te pysical link type. It is also easy to combine simplified and complete formulations witin a model. Simplified treatment of a sub domain just requires to ignore te convective terms in vectors F and G. - It provides numerical robustness, - It is easy to introduce singularities witin adjacent cells. Witin a 2D domain te flux components are treated by : - centered formulation for te volume fluxes, - upwind formulation for te momentum fluxes. 0.5(qq ii + qq jj ) 0 FF iiii. nn xx + GG iiii. nn yy = αα ii uu xxxx qq ii + αα jj uu xxxx qq jj + 0.5gg zz ıııı jj zz ii nn xx wit : αα ii uu yyyy qq ii + αα jj uu yyyy q jj 0.5gg zz ıııı jj zz ii nn yy uu xxxx uu yyyy : average flow velocity witin one cell i qq ii q xx,ii. nn xx + q yy,ii. nn yy qq jj q xx,jj. nn xx + q yyyy. nn yy αα ii = 1 si qq ii > 0 αα ii = 0 si qq ii < 0 αα jj = 1 si qq jj > 0 αα jj = 0 si qq jj < 0
3 Time integration of equation (3) is treated in a fully implicit way : Flux terms : FF nn+1 iiii = FF nn iiii + FF iiii UU UU ii + FF iiii UU ii UU jj jj GG ii nn+1 = GG ii nn + GG iiii UU ii UU ii + GG iiii UU jj UU jj UU ii nn+11 = UU ii nn + UU ii UU jj nn+11 = UU jj nn + UU jj Te Jacobian matrices : FF iiii UU ii, FF iiii UU jj, GG iiii UU ii, GG iiii UU jj are ranked 3x3 et contain differentiated terms wit respect to natives variables defined in vector U Inertia term: UU AA ii ii = AA ii UU dddd ii Friction term : SS nn+1 ffi = SS nn ffff + SS ffff UU UU ii ii Te final set of equations to be solved at eac time step is put in matrix form : MM iiii UU ii + KK iiii jj UU ii = FF ii (4) In practice tis matrix system is formed by calculating contribution of eac boundary to flux terms, ten contributions of volume terms witin eac cell. Solution of matrix system (4) is calculated using Pardiso Matrix Solver. Tis formulation encompasses te 1D case, it can terefore be specialized to 1D problems or extended to mixed problems including 1D and 2D sub domains witin a model witout difficulty. 2. SHOCK CAPTURE Te above formulation is based on native conservation laws of yperbolic St Venant equations. It contains all te terms wic are supposed to produce discontinuities depending on pysical configurations. In order to demonstrate te ability of te above algoritm to describe socks te restricted 1D case wit 3 adjacent cells is considered : M1 M2 M3 Te discretized equation (3) along te x axis is written as follows in steady state flow regime for te cell M2. Te following equation is obtained by neglecting friction terms : QQuu 1 + QQuu 2 + gg 2 (zz 2 zz 1 ) 12 + gg 2 (zz 3 zz 2 ) 23 = 0 (5) 12 = 0.5( ) 23 = 0.5( )
4 Let us assume tat a ydraulic jump occurs between cell 1 and 2 : One as : zz 2 zz 3 z 1 z z 3 3 Q ( flow rate) Equation (5) becomes : QQuu 1 + QQuu 2 + gg ( ) 12 = 0 QQuu 2 + gg = QQuu 1 + gg (6) Equation (6) describes a jump equation wic is linearized wit respect to term 12. Tis sows tat our formulation allows to capture socks, albeit introduction of some distortion for strong socks. Te practical validity or tis simplification is discussed wit te tests case presented in capter 4. Te above analysis can be extended easily to transient sock waves, suc as in te dam break problem, by taking in account te inertia terms of equation (3) : it can be cecked tat te rigt expression for te sock velocity of te moving sock wave is obtained. Tis sock analysis can extended to te 2D case. A furter simplification is introduced by assuming tat te local direction of a 2D sock wave is perpendicular to te boundary between two adjacent cells. Starting from equation (3) one obtains te equations for te oblique ydraulic jump, assuming tat te boundary line between two cells is a sock line. Tis additional assumption yields wrong results locally but it as been found tat te distortion levels out wen one considers te resulting wave pattern on a scale involving a few cells togeter. 3. EXTENDED FORMULATION WITH HYDRAULIC SINGULARITIES So far te formulation above applies to an omogenous 2D domain : two adjacent cells are connected by relations involving flux terms. It is fairly straigtforward to generalize tis concept by adapting te flux terms to ydraulic excange laws involving singular ead losses: Excange between te two cells is controlled by te flow rate Q et te momentum flux vector F. In case of a singularity linking te cells, Q and F are expressed as : Q(Hi, Hj) and F(Hi, Hj) were H stands for te averaged total ead witin one cell. Te expressions for Q and F depend on te nature of te singularity.
5 Assuming a flow from cell i to cell j te flux leaving cell i is written as : Q qq FFFFFFFF iiii = ii uu ii qq ii vv ii And te flux entering flux j is written as : FFFFFFFF jjjj = Q FF xx FF yy As an example te gate let us consider te gate type singularity wic is governed by a variety of laws according to te gate geometry and ydraulic upstream and downstream conditions. Figure 1 : definition sketc for te gate/weir singularity Next diagrams illustrates all te laws wic govern tis singularity and te range of validity of eac. Te flow equations are defined suc tat flow continuity is ensured at every transition between two flow regimes. Equations are given below for 4 basic regimes: Weir flow regime: HH 1 pp 1 < 3 2 vv : QQQQ DD = 2 3 CC ccbb 2gg(HH 1 pp 1 ) 3 2 QQQQ NN = bb 2gg[(HH 1 pp 1 ) (HH 2 pp 2 )] 1 2(HH 2 pp 2 ) Gate flow regime: HH 1 pp 1 > 3 2 vv QQQQ DD = bbcc cc vv 2gg(HH 1 pp 1 ) QQQQ NN = KKKK vv 2gg((HH 1 pp 1 ) (HH 2 pp 2 )) 1/2 Figure 2: delimitation of flow regimes for te gate/weir singularity Special consideration must be given to te momentum flux just downstream of te gate. Iterative calculations are necessary to determine te proper flow regime and to apply te corresponding flow equations. In case of an incident flow oblique to te gate te upstream momentum flux vector is decomposed into two components, one normal to te gate, te oter one parallel to te gate. Te flux parallel to te gate is supposed uncanged, wile te oter flux component is treated by te 1D gate equation system.
6 3. BENCHMARKING Te above formulation as been implemented into te Hydra computational software and as been evaluated against analytical solutions, field data and existing referenced codes suc as Telemac ( EDF) and Rubard (IRSTEA). Some selected examples are presented below to demonstrate te capabilities of te code: 3.1 Hydraulic jump in rectangular cannel Cannel widt is discretized into 3 meses 5m eac. Cannel lengt is discretized into 60 meses : - 20 cells 10x5m in te upstream part of te canal wic as a 2% slope, - 40 cells 5x5m in te downstream part wic is flat. Flow rate is 60 m 3 /s. Figure below sows te computed water profile along te centerline. It is compared wit te analytical solution: Analytical solution Canal invert 2D simulation Figure 3: water profile elevation along a rectangular cannel : computed and analytical solutions Te two curves are in good agreement, given te spatial discretization used for te computations. Difference on te upstream end is due to te fact tat no upstream limit condition is imposed : te code selects critical flow regime. Te water profile ten connects gradually towards uniform flow. 3.2 flow past a gate Cannel discretization and flow conditions are te same as above but a gate is introduced at abscissa 0.4km. Gate opening dept is 50cm. te grap below sows te computed water profiles using complete formulation and simplified formulation ( witout convective terms) :
7 Ca na l in ve rt Canal invert Complete numerical solution Simplifiednumerical solution Figure 4: computed water profile elevation along a rectangular cannel wit a gate singularity Water profile computed wit te full equations exibits a ydraulic jump in te upper section of te cannel and anoter one about 25m downstream of te gate. Te flow past te gate is not influenced by te downstream conditions. In te simplified calculations convective terms are ignored : super critical flow is not modelled, wic results in a significant over estimates of te water profiles upstream of te gate D flow in a river past a navigation weir. A section of te Aisne river as been modelled to analyze velocity field past a navigation weir and to investigate resulting geo-morpological consequences. Te model consists of about meses arranged as follows:
8 Navigation Weir Flow input ndition z(q) Figure 5 : 2D mesing for te river modelling pas a navigation weir Te weir is closed except along a 25m section next to te rigt bank. It is oblique to te upstream flow. River flow rate is 100 m 3 /s. Te weir elevation is adjusted so tat te flow past it is critical. A close up view of te flow velocity distribution is sown below : ydra code Figure 6 : flow distribution computed by ydra numerical code One can see clearly te two recirculation zones induced by te jet flux past te weir. Tis pattern was actually observed on te site. Te same calculations was performed using Telemac code. Te computed flow distribution exibit a very similar pattern as sown below. In te Telemac model te
9 mes size around te weir is muc more refined tan in te Hydra model because te weir is modelled as a geometric singularity : it is totally immerged witin te 2D mesing. Telemac code Figure 7: flow distribution computed by telemac numerical code 3.4 Transient dam break problem in a rectangular cannel. Tis problem is classic and as an analytical solution. Te case investigated is defined by te following definition sketc : Mesing Dam position expansion wave Sock wave Figure 8 : definition sketc for te dam break problem
10 Initially water dept is 5m upstream of te dam and 1m downstream. Te dam breaks at time=0. Figures below sow te longitudinal water and velocity profiles along te canal at time t = 288seconds. Numerical solution Analytical solution Figure 9 : water elevation and velocity profiles at time 288seconds after sudden dam break. Examination of tese graps sow near perfect agreement between computation and analytical solution. Tis good agreement suggests tat te approximation introduced in equation (3) is valid to solve practical problems involving socks.
11 4 DISCUSSION AND CONCLUSION Te implicit formulation presented in tis paper to solve full equations of free surface flow is believed to mark a progress wic departs from existing metods. It combines several advantages: unconditional numerical stability, computational robustness, accuracy, flexibility for solving problems in a wide variety of fields: estuary and river ydraulics, flood propagation and networks ydraulics. Finite volume formulation is based on te native conservative form of te St Venant equations witout any assumptions on te solution structure across a ydraulic sock. It includes any kind of ydraulics singularities in transient as well as steady state flow regimes and to combine simplified and full resolution of equations in te same model. Formulation is also quite suited to model 1D-2D interactions witin a single model. Calculations are fast: inclusion of singularities and constraints lines into a model enables to coose fairly loose mesing as compared wit full 2D codes for te same accuracy. Tis results in considerable time savings in computations. Tis upgraded formulation is implemented into Hydra software and as sown particular efficiency in all modelling problems involving floods and inundations in urban areas. Hydra is now coupled to a user s interface immerged in te Open Source Qgis GIS ( ref[4]). Coupling between tis interface and te computational code, sould make Hydra attractive among te ydro informatic software available on te market. BIBLIOGRAPHY 1. J.M. Hervouet - Hydrodynamics of free surface flow. Modelling wit Finite Element metod - Wiley ed A Paquier - Rupture de barrage : validation des modèles numériques du CEMAGREF dans le cadre de CADAM Ingenierie IRSTEA edition E. Toro Sock capturing metods for free surface flows Wiley ed Development of a ydraulic modelling platform witin an Open Source environment C Duran, T. Lepelletier, V. Mora Symydro symposium. GLOSSARY ( ux uy ) : velocity components in te x and y direction. g : acceleration of gravity : water dept i : mean water dept witin cell i lij : segment lengt between two adjacent cells Ai : cell area zb : bottom elevation U : vector composed of unknown variables ( qx qy ) ( qx qy ) : = ( ux uy ) H : pressure ead = + uu xx 2 +uu2 yy 2gg SS ff : bottom friction vector wit components : SS ff,xx = gg (nn² xx uu xx uu 2 xx + uu2 yy 4/3 ) et SS ff,yy = gg(nn² yy uu 2 xx + uu2 yy 4/3 ) ( nx ny ) : Manning coefficients in Ox et Oy direction.
Quantization of electrical conductance
1 Introduction Quantization of electrical conductance Te resistance of a wire in te classical Drude model of metal conduction is given by RR = ρρρρ AA, were ρρ, AA and ll are te conductivity of te material,
More informationManipulator Dynamics (1) Read Chapter 6
Manipulator Dynamics (1) Read Capter 6 Wat is dynamics? Study te force (torque) required to cause te motion of robots just like engine power required to drive a automobile Most familiar formula: f = ma
More informationPart 2: Introduction to Open-Channel Flow SPRING 2005
Part : Introduction to Open-Cannel Flow SPRING 005. Te Froude number. Total ead and specific energy 3. Hydraulic jump. Te Froude Number Te main caracteristics of flows in open cannels are tat: tere is
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationChapter 1 Functions and Graphs. Section 1.5 = = = 4. Check Point Exercises The slope of the line y = 3x+ 1 is 3.
Capter Functions and Graps Section. Ceck Point Exercises. Te slope of te line y x+ is. y y m( x x y ( x ( y ( x+ point-slope y x+ 6 y x+ slope-intercept. a. Write te equation in slope-intercept form: x+
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More information1watt=1W=1kg m 2 /s 3
Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,
More information6. Non-uniform bending
. Non-uniform bending Introduction Definition A non-uniform bending is te case were te cross-section is not only bent but also seared. It is known from te statics tat in suc a case, te bending moment in
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More information3. Gradually-Varied Flow
5/6/18 3. Gradually-aried Flow Normal Flow vs Gradually-aried Flow Normal Flow /g EGL (energy grade line) iction slope Geometric slope S Normal flow: Downslope component of weigt balances bed friction
More informationKeywords: quadratic programming, sparse matrix computation, symmetric indefinite factorization
A New Pivot Selection Algoritm for Symmetric Indefinite Factorization Arising in Quadratic Programming wit Block Constraint Matrices Duangpen Jetpipattanapong* and Gun Srijuntongsiri Scool of Information,
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationCFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E.
CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E. Publised in: IEA Annex 41 working meeting, Kyoto, Japan Publised:
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationCFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow. Annex 41 Kyoto, April 3 rd to 5 th, 2006
CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Annex 41 Kyoto, April 3 rd to 5 t, 2006 Adam Neale 1, Dominique Derome 1, Bert Blocken 2 and Jan Carmeliet 2,3
More informationChapter 4 Derivatives [ ] = ( ) ( )= + ( ) + + = ()= + ()+ Exercise 4.1. Review of Prerequisite Skills. 1. f. 6. d. 4. b. lim. x x. = lim = c.
Capter Derivatives Review of Prerequisite Skills. f. p p p 7 9 p p p Eercise.. i. ( a ) a ( b) a [ ] b a b ab b a. d. f. 9. c. + + ( ) ( + ) + ( + ) ( + ) ( + ) + + + + ( ) ( + ) + + ( ) ( ) ( + ) + 7
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationlim 1 lim 4 Precalculus Notes: Unit 10 Concepts of Calculus
Syllabus Objectives: 1.1 Te student will understand and apply te concept of te limit of a function at given values of te domain. 1. Te student will find te limit of a function at given values of te domain.
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationA PHYSICAL MODEL STUDY OF SCOURING EFFECTS ON UPSTREAM/DOWNSTREAM OF THE BRIDGE
A PHYSICA MODE STUDY OF SCOURING EFFECTS ON UPSTREAM/DOWNSTREAM OF THE BRIDGE JIHN-SUNG AI Hydrotec Researc Institute, National Taiwan University Taipei, 1617, Taiwan HO-CHENG IEN National Center for Hig-Performance
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationExercise 19 - OLD EXAM, FDTD
Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationAN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES
Ice in te Environment: Proceedings of te 1t IAHR International Symposium on Ice Dunedin, New Zealand, nd t December International Association of Hydraulic Engineering and Researc AN ANALYSIS OF AMPLITUDE
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationarxiv: v2 [cs.na] 22 Dec 2016
MONOTONICITY-PRESERVING FINITE ELEMENT SCHEMES BASED ON DIFFERENTIABLE NONLINEAR STABILIZATION SANTIAGO BADIA AND JESÚS BONILLA arxiv:66.8743v2 [cs.na] 22 Dec 26 Abstract. In tis work, we propose a nonlinear
More informationHydraulic validation of the LHC cold mass heat exchanger tube.
Hydraulic validation o te LHC cold mass eat excanger tube. LHC Project Note 155 1998-07-22 (pilippe.provenaz@cern.c) Pilippe PROVENAZ / LHC-ACR Division Summary Te knowledge o te elium mass low vs. te
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationHydraulic Evaluation of Discharge Over Rock Closing Dams on the Upper Mississippi River
ydraulic Evaluation of Discarge Over Rock Closing Dams on te Upper Mississippi River Jon endrickson, P.E. Senior ydraulic Engineer, St Paul District Introduction Prototype data was used for calibrating
More informationLines, Conics, Tangents, Limits and the Derivative
Lines, Conics, Tangents, Limits and te Derivative Te Straigt Line An two points on te (,) plane wen joined form a line segment. If te line segment is etended beond te two points ten it is called a straigt
More informationGRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES
GRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES Cristoper J. Roy Sandia National Laboratories* P. O. Box 5800, MS 085 Albuquerque, NM 8785-085 AIAA Paper 00-606 Abstract New developments
More informationRobotic manipulation project
Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationCritical control in transcritical shallow-water flow over two obstacles
Lougboroug University Institutional Repository Critical control in transcritical sallow-water flow over two obstacles Tis item was submitted to Lougboroug University's Institutional Repository by te/an
More informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More information2019 Chapter Competition Solutions
019 Capter Competition Solutions Are you wondering ow we could ave possibly tougt tat a Matlete would be able to answer a particular Sprint Round problem witout a calculator? Are you wondering ow we could
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More informationSimulations of the turbulent channel flow at Re τ = 180 with projection-based finite element variational multiscale methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Met. Fluids 7; 55:47 49 Publised online 4 Marc 7 in Wiley InterScience (www.interscience.wiley.com). DOI:./fld.46 Simulations of te
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationSeepage Analysis through Earth Dam Based on Finite Difference Method
J. Basic. Appl. Sci. Res., (11)111-1, 1 1, TetRoad Publication ISSN -44 Journal of Basic and Applied Scientific Researc www.tetroad.com Seepage Analysis troug Eart Dam Based on Finite Difference Metod
More informationCHAPTER 2 MODELING OF THREE-TANK SYSTEM
7 CHAPTER MODELING OF THREE-TANK SYSTEM. INTRODUCTION Te interacting tree-tank system is a typical example of a nonlinear MIMO system. Heiming and Lunze (999) ave regarded treetank system as a bencmark
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationNOTES ON OPEN CHANNEL FLOW
NOTES ON OPEN CANNEL FLOW Prof. Marco Pilotti Facoltà di Ingegneria, Università degli Studi di Brescia Profili di moto permanente in un canale e in una serie di due canali - Boudine, 86 OPEN CANNEL FLOW:
More informationINTRODUCTION DEFINITION OF FLUID. U p F FLUID IS A SUBSTANCE THAT CAN NOT SUPPORT SHEAR FORCES OF ANY MAGNITUDE WITHOUT CONTINUOUS DEFORMATION
INTRODUCTION DEFINITION OF FLUID plate solid F at t = 0 t > 0 = F/A plate U p F fluid t 0 t 1 t 2 t 3 FLUID IS A SUBSTANCE THAT CAN NOT SUPPORT SHEAR FORCES OF ANY MAGNITUDE WITHOUT CONTINUOUS DEFORMATION
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More informationUNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics
UNIVERSITY O SASKATCHEWAN Department of Pysics and Engineering Pysics Pysics 117.3 MIDTERM EXAM Regular Sitting NAME: (Last) Please Print (Given) Time: 90 minutes STUDENT NO.: LECTURE SECTION (please ceck):
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More informationWeek #15 - Word Problems & Differential Equations Section 8.2
Week #1 - Word Problems & Differential Equations Section 8. From Calculus, Single Variable by Huges-Hallett, Gleason, McCallum et. al. Copyrigt 00 by Jon Wiley & Sons, Inc. Tis material is used by permission
More information, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim
Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More information11.6 DIRECTIONAL DERIVATIVES AND THE GRADIENT VECTOR
SECTION 11.6 DIRECTIONAL DERIVATIVES AND THE GRADIENT VECTOR 633 wit speed v o along te same line from te opposite direction toward te source, ten te frequenc of te sound eard b te observer is were c is
More informationADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS
ADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS Yasuo NIHEI (1) and Takeiro SAKAI (2) (1) Department of Civil Engineering, Tokyo University of Science, 2641 Yamazaki,
More informationMODIFIED DIFFERENTIAL EQUATIONS. Dedicated to Prof. Michel Crouzeix
ESAIM: PROCEEDINGS, September 2007, Vol.21, 16-20 Gabriel Caloz & Monique Dauge, Editors MODIFIED DIFFERENTIAL EQUATIONS Pilippe Cartier 1, Ernst Hairer 2 and Gilles Vilmart 1,2 Dedicated to Prof. Micel
More informationNotes: Most of the material in this chapter is taken from Young and Freedman, Chap. 12.
Capter 6. Fluid Mecanics Notes: Most of te material in tis capter is taken from Young and Freedman, Cap. 12. 6.1 Fluid Statics Fluids, i.e., substances tat can flow, are te subjects of tis capter. But
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationComment on Experimental observations of saltwater up-coning
1 Comment on Experimental observations of saltwater up-coning H. Zang 1,, D.A. Barry 2 and G.C. Hocking 3 1 Griffit Scool of Engineering, Griffit University, Gold Coast Campus, QLD 4222, Australia. Tel.:
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationDistribution of reynolds shear stress in steady and unsteady flows
University of Wollongong Researc Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 13 Distribution of reynolds sear stress in steady
More informationTheoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump
TECHNICAL PAPE Teoretical Analysis of Flow Caracteristics and Bearing Load for Mass-produced External Gear Pump N. YOSHIDA Tis paper presents teoretical equations for calculating pump flow rate and bearing
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationTeaching Differentiation: A Rare Case for the Problem of the Slope of the Tangent Line
Teacing Differentiation: A Rare Case for te Problem of te Slope of te Tangent Line arxiv:1805.00343v1 [mat.ho] 29 Apr 2018 Roman Kvasov Department of Matematics University of Puerto Rico at Aguadilla Aguadilla,
More informationarxiv: v1 [math.na] 9 Mar 2018
A simple embedded discrete fracture-matrix model for a coupled flow and transport problem in porous media Lars H. Odsæter a,, Trond Kvamsdal a, Mats G. Larson b a Department of Matematical Sciences, NTNU
More information7.1 Using Antiderivatives to find Area
7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationRecent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm
Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationChapters 19 & 20 Heat and the First Law of Thermodynamics
Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,
More informationEntropy and the numerical integration of conservation laws
Pysics Procedia Pysics Procedia 00 2011) 1 28 Entropy and te numerical integration of conservation laws Gabriella Puppo Dipartimento di Matematica, Politecnico di Torino Italy) Matteo Semplice Dipartimento
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More information