Modelling the Zero-Inertia, Horizontal Viscous Dam-Break Problem
|
|
- Shanna Lamb
- 6 years ago
- Views:
Transcription
1 r WSEAS International Conference on APPLIED an TEORETICAL MECANICS, Spain, December 4-6, 7 8 Moelling the Zero-Inertia, orizontal Viscous Dam-Break Problem BLAISE NSOM, WILFRIED NDONG AND BLAISE RAVELO LIME/LBMS/PMB - Université e Bretagne Occientale IUT e Brest BP 969. Rue e Kergoat.. 9 BREST Ceex FRANCE blaise. nsom@univ-brest.fr Abstract: - Debris flows such as avalanches an lahars iffer from the classical am-break problem of hyraulics ue to the relative importance of viscous versus inertial forces in the momentum balance. An equation of motion escribing ebris flow in the limit of zero inertia is evelope an solve analytically in two limits: short time an long time. These solutions are then combine into a single, universal moel. Limitations of the moel are examine by comparison to a converge finite ifference numerical solution of the flow equation. Key-Wors: - Dam failure; Finite ifference metho; Flui height, Numerical moels; One imensional flow; Shallow water; Viscous flow, Wavefront. Introuction Consier a am obstructing a horizontal smooth channel, ry ownstream an with a given quantity of flui upstream (with heighth ), containe between a fixe plate an a am. At initial time, the am collapses an the flui is release ownstream (positive wave), while a negative wave propagates upstream (negative wave). From the amcollapse ate to the time when the negative wave reaches the fixe plate, Ritter [] gives the so-calle inertial solution, stating that the wavefront avances with a constant spee of gh, while the negative wave moves back with a constant spee of gh. Between these two extremities, the average spee U an the hyrograph are respectively given by ( x ) t ( gh x) U + gh = () gh () = t where the am is assume to be locate at x =. This configuration generally represents a flow generate by a am failure cause by an exceptional rainfall (e.g. Malpasset, France in 99) or by an act of war (e.g. Dnieproghes, Ukraine in 94). The flui is water an the flow is escribe by the Navier Stokes an continuity equations, together with the non slip conition. Assuming the shallow water approximation, this system of equations leas to the Barré e Saint-Venant equations [], a one-imensional hyperbolic system. The complete hyroynamic equations escribing this unsteay flow in an open channel were solve by Faure an Nahas [], using the metho of characteristics. unt [4], comparing a one-imensional turbulent flow moel own a slope with its viscous counterpart, conclue that the viscous flow moel gives the best escription for ebris flows. In this work, a -D moel is presente, aiming to provie practical laws useful to engineers. A priori knowlege of the spee of the floo wave is inee important because this will etermine the available time in which forecast an rescue measures nee to be effecte. Assuming the shallowwater approximation, equations of motion governing viscous am-break flow are built an put in non-imensional form, an the initial an bounary conitions are state. Then, an analytical solution is presente both for short time an long time behaviour. Finally, an explicit proceure was use here, which oes not take into account turbulence generate by ambreak wave, as the flow evelops over a ry smooth be []. Numerical results are shown an compare with the analytical ones in each regime.
2 r WSEAS International Conference on APPLIED an TEORETICAL MECANICS, Spain, December 4-6, Problem statement. - Equations of motion Let h enote the height of flui at negative time in a smooth horizontal rectangular channel, g the gravity, ρ an µ the flui ensity an viscosity, respectively. Using a cartesian system of coorinates with the origin at the am site, x-axis lying on the channellength an the z-axis in the increasing vertical irection (fig. ). Wall Flui h Dam -L X Fig : Configuration of horizontal am-break flow at negative time The flui is assume to flow mainly in the irection of x-axis with height h at the given control section of the abscissa x, at time t. So, the vertical velocities are negligibly small, an therefore the pressure is hyrostatic, the pressure in the flow is given by p= p + ρ g( h ) () z where p enotes the (constant) pressure at the free surface. The balance between the pressure graient an the viscous forces is thus expresse by p h u = g = ν ρ x x z (4) where the horizontal erivatives have been neglecte in comparison with the vertical erivatives on the right-han sie of (.4) because the length of the current is very much greater than its thickness. At the base of the flui layer the non slip conition writes u ( x,, t) = () Consiering that the shear stress at the top of the current is very much less than its value within the current, it can be approximate as u (6) z ( x, h, t) = the solution of (4), () an (6) is g u( x, z, t) = h z( h z) (7) ν x A complete etermination of the unknowns u an h requires the equation of continuity which can be written here as h + = h uz (8) t x Substituting (7) into (8) we obtain ( 4 h ρg h ) = t µ x (9) If l enotes the reservoir length, we can assume the following set of non imensional variables: x gh ( ', ', ', ') (,,, t) h h h x f ρ h x x f t h µ l = () where subscript f enotes the wave-front, the equation of motion (9) then becomes, in the non imensional form: ( h' 4) h ' ' = () x t Eq. is similar to the equation of motion obtaine by Schwarz [6] an Barthes-Biesel [7], escribing the evolution of a thin liqui layer flowing own a horizontal plane when surface tension effects can be neglecte.. - Initial an bounary conitions Using (), the flui height at initial time is given by: h ' ', ' x t= = for x ' () h ' ', ' x t= = otherwise () Furthermore, a complementary bounary conition shoul be impose upstream, assuming that a short time or an asymptotic solution is sought. These bounary conitions are suggeste by experimental observation. For the short time case, it is written as: '( x' = l', t' ) = h with l= ' h l (4) which means that only a given flui quantity in the upper part of the reservoir is release ownstream the very few moments following the am collapse. While for the long time case, it is written as: h () ( x' = l', t' ) ' ' = x
3 r WSEAS International Conference on APPLIED an TEORETICAL MECANICS, Spain, December 4-6, 7 8 which means that there is no flow at the fixe wall; so at that site, the free surface is horizontal. For convenience, in the rest of the paper, non imensional quantities will be enote by corresponing capital letters with no primes. - Analytical solution Dam-break flow belongs to the general class of gravity currents; so the solution epens on the time scale [8]. First of all, the inertial regime, characterize by a fixe height at the am-site hols immeiately after the am collapse []. Then, a solution ominate by viscous effects appears an tens to an asymptotic form. The solution sought here will give the analytical expression for a short time ( T << ) an a long time ( T >> ) viscous solutions, as well as the critical height (a fixe point at the am-site in the inertial regime).. - Short time solution Assume a solution of the form ( X, f'η [ ] = = X η (6) ( where the variable η is efine in the interval [,]. If such solution exists, the flui height at the am-site ( X=, is constant. Substituting (6) in () gives = ct (7) where the constant c is etermine from the solution of the following orinary ifferential equation: f 4( η) f( η) + cη = (8) η η f ( η ) = f ( η + ) = (9) A polynomial approximation of (7) in the form of a Taylor series at η = in the form ( ) [( ) ( ) ( ) ] η η +... f( η) = c 4 8 η + 4 () an substitution of (7) in () prouce c =.48 () which etermines the solution from which (constant) flui height at the am site is obtaine: η= =. () ( ) 664 For the long time solution, a conition expressing the mass conservation may be written as: ( X, X= () which means that the flow only retains the initial (non imensional) volume of the reservoir V =. L an not the etails of its initial geometry. We seek a similar solution of the form ( X, = (, ( η) η= X + ( T+ ) (4) where variable η is efine in the interval [,], while the functions (, (, an ( η) are etermine by irect substitution in the equation of motion () an in the equation of mass conservation (). By setting = A h ' ( η) η we obtain: X ( ) ( ) f T = B T+ D ( ), T = A[ ( + ] 4 B = A () (6) To etermine constants A an C, the following ifferential equation has to be solve: {[ ( )]} ( ) 8 4 η η + η + ( η) = η η with the bounary conitions ( η) = Finally, we get ( η) = ( ) η ( η) = η, =. 84 (7) η = (8) A, B=. 86 (9) To etermine the transition point from the short time regime to the long time regime, both solutions previously obtaine are equate an the equation of motion of the front wave is finally obtaine as X f =.86 T (). - Long time solution
4 r WSEAS International Conference on APPLIED an TEORETICAL MECANICS, Spain, December 4-6, Numerical solution To buil a numerical proceure, it is necessary to efine the total channel length l t. The non imensional extreme (ecreasing) abscissa is L e lt h l = () This point is locate so far from the am site that the flow is assume to never reach it uring a given experiment (D assumption) for a total uration τ. Then, the problem to solve numerically is escribe by the following equation of motion = T X 4 ( ) () associate with the following initial conitions: ( X,) = if X [ L,] () ( X,) = otherwise (4) an the bounary conitions: ( L, = T X () ( X e, = T (6) This problem is solve by a finite ifference metho. To o so, the function ( X, T ) is compute in the set [ L, Le] [,τ] an in turn iscretize into a finite number of ientical small rectangles with sies T an X. The equation will be approximate at the gri points locate at the following coorinates in set: the [ L, Le] [,τ] ( X i, T j ) = ( L+ i X, j [ L ] X L e, j [, τ ] i, + (7) T Notice that eq.() can be put in the form: = 4 + (8) T X X An heuristic approach [9-] consiers the prouct ( 4 ) in the right-han sie of eq.(8) as a coefficient of iffusion. Inee, the following equations are consiere : V V = 4 (9) T X an V ( X,) = ( X,) (4) The numerical scheme of the equations (9-4) is teste using the von Neumann metho to provie a stability criterion which is necessary to ensure the convergence of our non-linear problem. The erivative of the unknown function an the non linear term in equ. () were estimate using Taylor s formula. Then, the finite ifference equation to solve, which uses a first orer time scheme an a centre secon orer spatial scheme, is written as T 4 [( ) ( ) ] 4 i, j+ = i, j + i+, j + i, j ( ) T (46) i j ( ) (, )4 In orer to have a stability criterion, the equation (4) is iscretize following the same numerical scheme, i.e. a first orer time scheme an secon orer centre spatial scheme. The numerical problem is written as: V,, 4 T i j+ = Vi j+ V i+, j + V i, j V i, j ( ) ( ) (49) with the same bounary conitions as, i.e. : V = V an V = max, +., j+, j+ i By applying i max ˆ i kp /( i max ) V = V e π + () p, q k, q k= where V ˆk, q is the Fourier component corresponing to the wave number k at time T = q T, efine by : i max ˆ iπ kp /( i max+ ) Vk, q = Vp, qe () i max+ p= equation () is rewritten as: ˆ ˆ 6 T Vk, q Vk, q sin ( π k + = i max ) + () The stability criterion consists of consiering that: 6 sin ( kπ ) < XT i max+, k () Since sin ( kπ ) <, we obtain i max+ finally the following stability criterion j
5 r WSEAS International Conference on APPLIED an TEORETICAL MECANICS, Spain, December 4-6, 7 8 ( ) 8 X T (4) We can see that the numerical scheme escribe by the equation (46) gives + + if,. The front wave I, j I j X f velocity efine as Vf = must then t satisfy: Vf () T A McCormack finite ifference scheme can improve the accuracy of the solution when Ve. In our case, the time step must be T chosen small enough to satisfy the conition (). 4 Numerical results Notice that a complete escription of the flow shoul inclue surface tension, introucing a complementary term in the equation of motion, such as: = 4 T X X B X (6) where B enotes the Bon number, efine as ρ gl σ B= (7) an σ the flui surface tension. Computation of eq.(4) was carrie out using the proceure escribe in the previous section. The flui s physical properties (ensity, viscosity an surface tension) were taken in Weast et al. []. For a similar flow configuration, results were quite ientical to those obtaine from eq.(), i.e. when the surface tension is neglecte. In fact, the surface tension woul affect the viscous am-break flow only uner film lubrication conitions (e.g.: []). The results obtaine analytically an numerically are concorant. Moreover, they also agree with those obtaine experimentally by Nsom et al. [] an Nsom [4]. Notably, fig. shows that for a very short time after the am collapse (the short time solution), the flow height remains constant at the am-site with the value =. 664 an with a ecrease of the flui height in the long time regime. Upstream the flui height ecreases both in the short time an long time regimes, while ownstream at a given station the flui height abruptly increases, reaches a maximum an then ecreases. This maximum height an the corresponing time were accurately etermine. They are governe by the following equations: max ( ) ( ) ( ) X = α( X+ ) / ( X) ( ) ( X+ ) max = + T c ( αγ) (8) T (9) γ where α =. 84, γ =. 86 an T c =. 8. The time evolution of the wavefront was characterize. Three flow regimes coul be pointe out with concoring equations both analytically an numerically. An inertial regime vali immeiately after the am collapse corresponing to. 664 at the am-site an governe by T = 4 for T.4L ( ) T (6) meets Ritter s invisci solution []. Then, a viscous regime results with the following equation of motion: ( = T where.4l T.8L (6) This regime can be observe for a long perio for large values of L, which means that the asymptotic solution correspons to an infinite reservoir length. The short time viscous regime evolves into a long time viscous regime (fig.) T= T=. T=.4 T=.7 T=. T= T= X Fig.: Time variation of the free surface profile
6 r WSEAS International Conference on APPLIED an TEORETICAL MECANICS, Spain, December 4-6, 7 86 with the asymptotic equation ( T >> ) ( =.86T. 9 (6) The short time viscous regime evolves into a long time viscous regime (fig.) with the T >> asymptotic equation ( ) ( =.86T. 9 (6) Conclusion The horizontal viscous am-break flow was consiere. The channel was smooth an ry at the initial time. By applying the conservation of mass an momentum with the shallow water approximation, an equation of motion was erive an mae non imensional when the viscous forces were assume to be the ominant ones. Then, an analytical solution of the equation of motion was built. Immeiately after the am collapse, an inertial regime results. In this regime, the flow height at the am site has a (fix) characteristic value an the equation of motion of the front wave is a Ritter type. As the flow height insie the reservoir reaches this characteristic value, a short time viscous regime is establishe without any influence from the reservoir extening onto the front wave s equation of motion. Finally, a long time viscous results epening on the reservoir length. Also, a convex free surface profile with a self-similar form an a time variation of the flow height at both the upstream an ownstream stations were furnishe. Furthermore, the problem was consiere numerically. The previous equation of motion was of porous meium type an was approximate using an explicit finite ifference metho. The stability an convergence of the computations were insure using a criteria base on a heuristic approach. The very goo agreement between numerical an asymptotic solutions shows the consistency of the numerical scheme for both the short time an the asymptotic form of the long time solution which correspon respectively to the long an short reservoir length L.. References [] Ritter, A., Die fortpflanzung er wesser wellen. Ver Deutsch Ingenieure Zeitschrift, Vol.6, N, 89, pp (in German) [] e Saint-Venant, B., Théorie u mouvement non permanent es eaux, C. R. Aca. Sci., Vol.7, N, 87, pp.47-4 (in French) [] Faure, J. an Nahas, N., Etue numérique et expérimentale intumescences à forte courbure u front, La ouille Blanche, Vol.6, N, 96, pp (in French) [4] unt, B., Newtonian flui mechanics treatment of ebris flows an avalanches, J. yr. Div. (ASCE), Vol., N, 994, pp.-6 [] Shigematsu, T., Liu, P.L.F. an Oa, K., Numerical moeling of the initial stages of am-break waves, J. yr. Res. (IAR), Vol.4, N, 4, pp.8-9 [6] Schwarz, L.W., Viscous flows own an incline plane: Instability an finger formation, Phys. Fluis A, Vol., N, 989, pp [7] Barthes-Biesel, D., Rectification un film liquie sous l effet e la gravité et e la tension superficielle (in French) < ofs/informatique/georges.gonthier/pi98/film.h tml> (June, 998) [8] Simpson, J.E., Gravity Currents in the Environment an the Laboratory, John Wiley & Sons, New York, 987. [9] Smith, G.D., Numerical Solution of Partial Differential Equations, Oxfor University Press, New York, 969. [] Forsythe, G.E. an Wasow, W.R., Finite- Difference Methos for Partial Differential Equations, John Wiley an Sons, New York, 967. [] Richtmyer, R.D. an Morton, K.W., Difference Methos for Initial Value Problems, John Wiley an Sons, New York, 967. [] Weast, R.C., Astle,M.J. an Beyer, W.., es. anbook of chemistry an physics, CRC Boca Raton, Fla., 987. [] Nsom, B., Debiane, K. an Piau, J.M., Be slope effect in the am-break problem, J. yr. Res. (IAR), Vol.8, N 6,, pp [4] Nsom, B., orizontal Viscous Dam-Break Floo: Experiments an Theory, J. yr. Eng. (ASCE), Vol.8, N,, pp.4-46
Equipment Safety in Renewable Energies Exploitation
European Association or the Development o Renewable Energies, Environment and Power Quality (EAEPQ International Conerence on Renewable Energies and Power Quality (ICREPQ 11 Las Palmas de Gran Canaria
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
More informationMAE 210A FINAL EXAM SOLUTIONS
1 MAE 21A FINAL EXAM OLUTION PROBLEM 1: Dimensional analysis of the foling of paper (2 points) (a) We wish to simplify the relation between the fol length l f an the other variables: The imensional matrix
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationA Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion
Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationSection 2.7 Derivatives of powers of functions
Section 2.7 Derivatives of powers of functions (3/19/08) Overview: In this section we iscuss the Chain Rule formula for the erivatives of composite functions that are forme by taking powers of other functions.
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationTo understand how scrubbers work, we must first define some terms.
SRUBBERS FOR PARTIE OETION Backgroun To unerstan how scrubbers work, we must first efine some terms. Single roplet efficiency, η, is similar to single fiber efficiency. It is the fraction of particles
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More informationCompletely passive natural convection
Early View publication on wileyonlinelibrary.com (issue an page numbers not yet assigne; citable using Digital Object Ientifier DOI) ZAMM Z. Angew. Math. Mech., 1 6 (2011) / DOI 10.1002/zamm.201000030
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationAbstract A nonlinear partial differential equation of the following form is considered:
M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationComparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope
Proceeings of the Worl Congress on Engineering Vol WCE, July 6-8,, Lonon, U.K. Comparative Approaches of Calculation of the Back Water Curves in a Trapezoial Channel with Weak Slope Fourar Ali, Chiremsel
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More information6 Wave equation in spherical polar coordinates
6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More information2.20 Marine Hydrodynamics Lecture 3
2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationThermal Modulation of Rayleigh-Benard Convection
Thermal Moulation of Rayleigh-Benar Convection B. S. Bhaauria Department of Mathematics an Statistics, Jai Narain Vyas University, Johpur, Inia-3400 Reprint requests to Dr. B. S.; E-mail: bsbhaauria@reiffmail.com
More informationUnit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule
Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin
More informationMULTISCALE FRICTION MODELING FOR SHEET METAL FORMING
MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of
More informationAsymptotics of a Small Liquid Drop on a Cone and Plate Rheometer
Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More informationMany problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx
Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More informationObjective: To introduce the equations of motion and describe the forces that act upon the Atmosphere
Objective: To introuce the equations of motion an escribe the forces that act upon the Atmosphere Reaing: Rea pp 18 6 in Chapter 1 of Houghton & Hakim Problems: Work 1.1, 1.8, an 1.9 on pp. 6 & 7 at the
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationSPE Copyright 1999, Society of Petroleum Engineers Inc.
SPE 664 Effect of Flow Through a Choke Valve on Emulsion Stability M.J. van er Zane, SPE, K.R. van Heuven, J.H. Muntinga, SPE, an W.M.G.T. van en Broek, SPE, Delft University of Technology Copyright 1999,
More informationIntroduction to Markov Processes
Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav
More informationA 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS
A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS Peifeng Ma 1 an Ole Secher Masen Accurate preiction of current velocity an bottom shear stress, which both can be significantly influence by
More informationA SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS
International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical
More informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More informationAn Analytical Expression of the Probability of Error for Relaying with Decode-and-forward
An Analytical Expression of the Probability of Error for Relaying with Decoe-an-forwar Alexanre Graell i Amat an Ingmar Lan Department of Electronics, Institut TELECOM-TELECOM Bretagne, Brest, France Email:
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More information(2012) , ISBN
Cregan, V. an O'Brien, Stephen B.G. an McKee, Sean (2012) Asymptotics of a small liqui rop on a cone an plate rheometer. In: Progress in Inustrial Mathematics at ECMI 2010. Mathematics in Inustry: Progress
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationProblem 1 (20 points)
ME 309 Fall 01 Exam 1 Name: C Problem 1 0 points Short answer questions. Each question is worth 5 points. Don t spen too long writing lengthy answers to these questions. Don t use more space than is given.
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationFLUID MECHANICS UNIVERSITY OF LEEDS. May/June Examination for the degree of. BEng/ MEng Civil Engineering. Time allowed: 2 hours
This question paper consists of printe pages, each of which is ientifie by the Coe Number CIVE 4 UNIVERSITY OF LEEDS May/June Examination for the egree of BEng/ MEng Civil Engineering FLUID MECANICS Time
More informationA note on the Mooney-Rivlin material model
A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel
More informationFinal Exam Study Guide and Practice Problems Solutions
Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making
More informationfv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n
Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie
More informationNON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS
ANZIAM J. 46(2004), 1 16 NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS A. C. MCINTOSH 1,R.O.WEBER 2 ang.n.mercer 2 (Receive 14 August, 2002) Abstract This
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationSolutions to Math 41 Second Exam November 4, 2010
Solutions to Math 41 Secon Exam November 4, 2010 1. (13 points) Differentiate, using the metho of your choice. (a) p(t) = ln(sec t + tan t) + log 2 (2 + t) (4 points) Using the rule for the erivative of
More information