Math 575-Lecture 25 Deep Water Waves, Dispersion Relation

Total Page:16

File Type:pdf, Size:1020Kb

Math 575-Lecture 25 Deep Water Waves, Dispersion Relation Math 575-Lecture 25 Deep water waves, dispersion relation The full 2D equations describing the water are given by the Euler equations 1 @ u + u@ u + v@ u = − @ p; t x y ρ x 1 @ v + u@ v + v@ v = − @ p − g; t x y ρ y @xu + @yv = 0: We aim to find the equations for the surface waves. Assume that the water is deep and the flows are irrotational: u = rφ. Then, the Euler equations are reduced to the Bernoulli equations 1 @ φ + juj2 + (p/ρ + gy) = 0; t 2 and ∆φ = 0: At the surface y = h(x; t), we have v = ht + uhx; at y = h(x; t): and if we impose p = 0, 1 @ φ + (u2 + v2) + gh = 0; at y = h(x; t): t 2 Since the water is deep, h = h0 + η; with = η=h0 1. This is called the small amplitude assumption. By throwing away a totally time dependent constant in φ, we are able to keep the leading order effects: 2 v(x; h0 + η) = ht + O( ); 2 @tφ(x; h0 + η) + gη + O( ) = 0 Further, by Taylor expansion, we find to leading order, @yφ(x; h0; t) = v(x; h0) = ηt;@tφ(x; h0; t) + gη = 0 1 The deep water is then described by @tφ(x; t) + gη = 0; y = h0 ηt = @yφ(x; t); y = h0; φxx + φyy = 0: Inserting the ansatz η = A sin(kx − !t), we find the dispersion relation !2 = gjkj: The energy is tranported with group velocity d! v = g dk If k ∼ 0, the wavelength is big, the approximation we did p = 0 at y = h is reasonable and the equations we have are good. Such water waves are called the gravity waves. If k ∼ 1, p = 0 at y = h is no longer good. The surface tension matters . The first equation should be revised to T @ φ(x; t) + gη − η = 0: t ρ xx These waves are called `capillary waves'. This is the case when a drop falls onto the surface of lake. The dispersion relation is then !2 = gjkj + T k2/ρ. Shallow water equations Consider the water above the ground y = 0. The free surface is given by y = h(x; t). h0 is the typical value of h. We are not assuming that the amplitude of the wave is small. Instead, we assume h0=L 1 where L is the typical length of the domain (or wavelength). This is the shallowness assumption. Consider the Euler equations again: 1 @ u + u@ u + v@ u = − @ p; t x y ρ x 1 @ v + u@ v + v@ v = − @ p − g; t x y ρ y @xu + @yv = 0: 2 0 0 0 0 As before, we scale x; y by x = x=L and y = y=h0. u = u=U and v = v=(Uh0=L), 0 p t = t=(L=U), U = gh0. We assume we are considering the regime where all prime numbers are O(1). By finding the leading order terms, we have 1 @ u + u@ u + v@ u = − @ p; t x y ρ x 1 0 = − @ p − g; ρ y @xu + @yv = 0: By the second equation: p = p0 − ρg(y − h(x; t)) Hence, the first equation is reduced to @tu + u@xu + v@yu = −ghx Consider that u(x; y; 0) = u0(x), i.e. it does not depend on y. Then, according to the above equation, u(x; y; t) = u(x; t) for all time. Hence, @tu + u@xu = −ghx By the incompressibility condition, v = −@xuy; because v = 0 for y = 0. Lastly, consider a fluid particle y = h(x; t). It then implies v = ht + uhx for y = h(x; t), or ht + uhx + @xuh(x; t) = 0. We then have the Shallow water equations ut + uux + ghx = 0; ht + uhx + uxh(x; t) = 0: In conservation law form 2 ut + (u =2 + gh)x = 0; ht + (uh)x = 0: or (ρht) + (ρhu)x = 0; g (ρh2) + (ρhu2) + ( (ρh)2) = 0: t x 2ρ x 3 p Letting c = gh, we can rewrite the equations into (@t + (u + c)@x)(u + 2c) = 0; (@t + (u − c)@x)(u − 2c) = 0: Such kind of nonlinear hyperbolic systems can develop discontinuity, called `shock' even though the initial value is smooth. 4.
Recommended publications
  • Dispersion of Tsunamis: Does It Really Matter? and Physics and Physics Discussions Open Access Open Access S
    EGU Journal Logos (RGB) Open Access Open Access Open Access Advances in Annales Nonlinear Processes Geosciences Geophysicae in Geophysics Open Access Open Access Nat. Hazards Earth Syst. Sci., 13, 1507–1526, 2013 Natural Hazards Natural Hazards www.nat-hazards-earth-syst-sci.net/13/1507/2013/ doi:10.5194/nhess-13-1507-2013 and Earth System and Earth System © Author(s) 2013. CC Attribution 3.0 License. Sciences Sciences Discussions Open Access Open Access Atmospheric Atmospheric Chemistry Chemistry Dispersion of tsunamis: does it really matter? and Physics and Physics Discussions Open Access Open Access S. Glimsdal1,2,3, G. K. Pedersen1,3, C. B. Harbitz1,2,3, and F. Løvholt1,2,3 Atmospheric Atmospheric 1International Centre for Geohazards (ICG), Sognsveien 72, Oslo, Norway Measurement Measurement 2Norwegian Geotechnical Institute, Sognsveien 72, Oslo, Norway 3University of Oslo, Blindern, Oslo, Norway Techniques Techniques Discussions Open Access Correspondence to: S. Glimsdal ([email protected]) Open Access Received: 30 November 2012 – Published in Nat. Hazards Earth Syst. Sci. Discuss.: – Biogeosciences Biogeosciences Revised: 5 April 2013 – Accepted: 24 April 2013 – Published: 18 June 2013 Discussions Open Access Abstract. This article focuses on the effect of dispersion in 1 Introduction Open Access the field of tsunami modeling. Frequency dispersion in the Climate linear long-wave limit is first briefly discussed from a the- Climate Most tsunami modelers rely on the shallow-water equations oretical point of view. A single parameter, denoted as “dis- of the Past for predictions of propagationof and the run-up. Past Some groups, on persion time”, for the integrated effect of frequency dis- Discussions the other hand, insist on applying dispersive wave models, persion is identified.
    [Show full text]
  • Fluid Mechanics
    FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER ANKARA UNIVERSITY FACULTY OF AGRICULTURE DEPARTMENT OF AGRICULTURAL MACHINERY AND TECHNOLOGIES ENGINEERING 1 1. INTRODUCTION Mechanics is the oldest physical science that deals with both stationary and moving bodies under the influence of forces. Mechanics is divided into three groups: a) Mechanics of rigid bodies, b) Mechanics of deformable bodies, c) Fluid mechanics Fluid mechanics deals with the behavior of fluids at rest (fluid statics) or in motion (fluid dynamics), and the interaction of fluids with solids or other fluids at the boundaries (Fig.1.1.). Fluid mechanics is the branch of physics which involves the study of fluids (liquids, gases, and plasmas) and the forces on them. Fluid mechanics can be divided into two. a)Fluid Statics b)Fluid Dynamics Fluid statics or hydrostatics is the branch of fluid mechanics that studies fluids at rest. It embraces the study of the conditions under which fluids are at rest in stable equilibrium Hydrostatics is fundamental to hydraulics, the engineering of equipment for storing, transporting and using fluids. Hydrostatics offers physical explanations for many phenomena of everyday life, such as why atmospheric pressure changes with altitude, why wood and oil float on water, and why the surface of water is always flat and horizontal whatever the shape of its container. Fluid dynamics is a subdiscipline of fluid mechanics that deals with fluid flow— the natural science of fluids (liquids and gases) in motion. It has several subdisciplines itself, including aerodynamics (the study of air and other gases in motion) and hydrodynamics (the study of liquids in motion).
    [Show full text]
  • Fluid Viscosity and the Attenuation of Surface Waves: a Derivation Based on Conservation of Energy
    INSTITUTE OF PHYSICS PUBLISHING EUROPEAN JOURNAL OF PHYSICS Eur. J. Phys. 25 (2004) 115–122 PII: S0143-0807(04)64480-1 Fluid viscosity and the attenuation of surface waves: a derivation based on conservation of energy FBehroozi Department of Physics, University of Northern Iowa, Cedar Falls, IA 50614, USA Received 9 June 2003 Published 25 November 2003 Online at stacks.iop.org/EJP/25/115 (DOI: 10.1088/0143-0807/25/1/014) Abstract More than a century ago, Stokes (1819–1903)pointed out that the attenuation of surface waves could be exploited to measure viscosity. This paper provides the link between fluid viscosity and the attenuation of surface waves by invoking the conservation of energy. First we calculate the power loss per unit area due to viscous dissipation. Next we calculate the power loss per unit area as manifested in the decay of the wave amplitude. By equating these two quantities, we derive the relationship between the fluid viscosity and the decay coefficient of the surface waves in a transparent way. 1. Introduction Surface tension and gravity govern the propagation of surface waves on fluids while viscosity determines the wave attenuation. In this paper we focus on the relation between viscosity and attenuation of surface waves. More than a century ago, Stokes (1819–1903) pointed out that the attenuation of surface waves could be exploited to measure viscosity [1]. Since then, the determination of viscosity from the damping of surface waves has receivedmuch attention [2–10], particularly because the method presents the possibility of measuring viscosity noninvasively. In his attempt to obtain the functionalrelationship between viscosity and wave attenuation, Stokes observed that the harmonic solutions obtained by solving the Laplace equation for the velocity potential in the absence of viscosity also satisfy the linearized Navier–Stokes equation [11].
    [Show full text]
  • Lecture 18 Ocean General Circulation Modeling
    Lecture 18 Ocean General Circulation Modeling 9.1 The equations of motion: Navier-Stokes The governing equations for a real fluid are the Navier-Stokes equations (con­ servation of linear momentum and mass mass) along with conservation of salt, conservation of heat (the first law of thermodynamics) and an equation of state. However, these equations support fast acoustic modes and involve nonlinearities in many terms that makes solving them both difficult and ex­ pensive and particularly ill suited for long time scale calculations. Instead we make a series of approximations to simplify the Navier-Stokes equations to yield the “primitive equations” which are the basis of most general circu­ lations models. In a rotating frame of reference and in the absence of sources and sinks of mass or salt the Navier-Stokes equations are @ �~v + �~v~v + 2�~ �~v + g�kˆ + p = ~ρ (9.1) t r · ^ r r · @ � + �~v = 0 (9.2) t r · @ �S + �S~v = 0 (9.3) t r · 1 @t �ζ + �ζ~v = ω (9.4) r · cpS r · F � = �(ζ; S; p) (9.5) Where � is the fluid density, ~v is the velocity, p is the pressure, S is the salinity and ζ is the potential temperature which add up to seven dependent variables. 115 12.950 Atmospheric and Oceanic Modeling, Spring '04 116 The constants are �~ the rotation vector of the sphere, g the gravitational acceleration and cp the specific heat capacity at constant pressure. ~ρ is the stress tensor and ω are non-advective heat fluxes (such as heat exchange across the sea-surface).F 9.2 Acoustic modes Notice that there is no prognostic equation for pressure, p, but there are two equations for density, �; one prognostic and one diagnostic.
    [Show full text]
  • Modelling Viscous Free Surface Flow
    Chapter 8 Modelling Viscous Free Surface Flow Free surface flows, where a boundary of a fluid body is free to move constrained only by forces across the surface, are possibly the most commonly observed flow phenomenon, with the motion of the free surface readily allowing the observation of flow of the fluid. Flows that are commonly encountered by the layperson include the motion of the surface of a river, the waves on the surface of the ocean, and the more personal flow of that in a cup of tea. From an engineering perspective, flows of interest include the open channel flow of rivers and canals, the erosive forces of waves on the shoreline, and the wakes generated by ships when under way. When modelling these flows, the problems associated with the solution of the Navier–Stokes equations are compounded by the free motion of the surface of the fluid, with the boundaries of the flow domain being a function of the flow structure, and thus an unknown which must be calculated along with the flow field. The movement of the free surface therefore both aids the observation of the flow, and hinders it’s modelling. In this chapter an attempt is made to model the steady flow around a ships hull. The modelling of such a flow is of great interest to Naval Architects, with the ability to predict the resistance of ships allowing the design of more efficient hull forms, whilst the accurate modelling of the ships wake allows the design of hulls that create a smaller disturbance in confined waterways.
    [Show full text]
  • Waves and Weather
    Waves and Weather 1. Where do waves come from? 2. What storms produce good surfing waves? 3. Where do these storms frequently form? 4. Where are the good areas for receiving swells? Where do waves come from? ==> Wind! Any two fluids (with different density) moving at different speeds can produce waves. In our case, air is one fluid and the water is the other. • Start with perfectly glassy conditions (no waves) and no wind. • As wind starts, will first get very small capillary waves (ripples). • Once ripples form, now wind can push against the surface and waves can grow faster. Within Wave Source Region: - all wavelengths and heights mixed together - looks like washing machine ("Victory at Sea") But this is what we want our surfing waves to look like: How do we get from this To this ???? DISPERSION !! In deep water, wave speed (celerity) c= gT/2π Long period waves travel faster. Short period waves travel slower Waves begin to separate as they move away from generation area ===> This is Dispersion How Big Will the Waves Get? Height and Period of waves depends primarily on: - Wind speed - Duration (how long the wind blows over the waves) - Fetch (distance that wind blows over the waves) "SMB" Tables How Big Will the Waves Get? Assume Duration = 24 hours Fetch Length = 500 miles Significant Significant Wind Speed Wave Height Wave Period 10 mph 2 ft 3.5 sec 20 mph 6 ft 5.5 sec 30 mph 12 ft 7.5 sec 40 mph 19 ft 10.0 sec 50 mph 27 ft 11.5 sec 60 mph 35 ft 13.0 sec Wave height will decay as waves move away from source region!!! Map of Mean Wind
    [Show full text]
  • Basin Scale Tsunami Propagation Modeling Using Boussinesq Models: Parallel Implementation in Spherical Coordinates
    WCCE – ECCE – TCCE Joint Conference: EARTHQUAKE & TSUNAMI BASIN SCALE TSUNAMI PROPAGATION MODELING USING BOUSSINESQ MODELS: PARALLEL IMPLEMENTATION IN SPHERICAL COORDINATES J. T. Kirby1, N. Pophet2, F. Shi1, S. T. Grilli3 ABSTRACT We derive weakly nonlinear, weakly dispersive model equations for propagation of surface gravity waves in a shallow, homogeneous ocean of variable depth on the surface of a rotating sphere. A numerical scheme is developed based on the staggered-grid finite difference formulation of Shi et al (2001). The model is implemented using the domain decomposition technique in conjunction with the message passing interface (MPI). The efficiency tests show a nearly linear speedup on a Linux cluster. Relative importance of frequency dispersion and Coriolis force is evaluated in both the scaling analysis and the numerical simulation of an idealized case on a sphere. 1 INTRODUCTION The conventional models in the global-scale tsunami modeling are based on the shallow water equations and neglect frequency dispersion effects in wave propagation. Recent studies on tsunami modeling revealed that such tsunami models may not be satisfactory in predicting tsunamis caused by nonseismic sources (Løvholt et al., 2008). For seismic tsunamis, the frequency dispersion effects in the long distance propagation of tsunami fronts may become significant. The numerical simulations of the 2004 Indian Ocean tsunami by Glimsdal et al. (2006) and Grue et al. (2008) indicated the undular bores may evolve in shallow water, as the phenomenon evidenced in observations (Shuto, 1985). In the simulation for the same tsunami by Grilli et al. (2007), the dispersive effects were quantified by running the dispersive Boussinesq model FUNWAVE (Kirby et al., 1998) and the NSWE solver.
    [Show full text]
  • A Simple and Unified Linear Solver for Free-Surface and Pressurized
    water Technical Note A Simple and Unified Linear Solver for Free-Surface and Pressurized Mixed Flows in Hydraulic Systems Dechao Hu 1 , Songping Li 1, Shiming Yao 2 and Zhongwu Jin 2,* 1 School of Hydropower and Information Engineering, Huazhong University of Science and Technology, Wuhan 430074, China; [email protected] (D.H.); [email protected] (S.L.) 2 Yangtze River Scientific Research Institute, Wuhan 430010, China; [email protected] * Correspondence: [email protected]; Tel.: +86-027-8282-9873 Received: 19 August 2019; Accepted: 19 September 2019; Published: 23 September 2019 Abstract: A semi–implicit numerical model with a linear solver is proposed for the free-surface and pressurized mixed flows in hydraulic systems. It solves the two flow regimes within a unified formulation, and is much simpler than existing similar models for mixed flows. Using a local linearization and an Eulerian–Lagrangian method, the new model only needs to solve a tridiagonal linear system (arising from velocity-pressure coupling) and is free of iterations. The model is tested using various types of mixed flows, where the simulation results agree with analytical solutions, experiment data and the results reported by former researchers. Sensitivity studies of grid scales and time steps are both performed, where a common grid scale provides grid-independent results and a common time step provides time-step-independent results. Moreover, the model is revealed to achieve stable and accurate simulations at large time steps for which the CFL is greater than 1. In simulations of a challenging case (mixed flows characterized by frequent flow-regime conversions and a closed pipe with wide-top cross-sections), an artificial slot (A-slot) technique is proposed to cope with possible instabilities related to the discontinuous main-diagonal coefficients of the linear system.
    [Show full text]
  • Shallow Water Waves and Solitary Waves Article Outline Glossary
    Shallow Water Waves and Solitary Waves Willy Hereman Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, Colorado, USA Article Outline Glossary I. Definition of the Subject II. Introduction{Historical Perspective III. Completely Integrable Shallow Water Wave Equations IV. Shallow Water Wave Equations of Geophysical Fluid Dynamics V. Computation of Solitary Wave Solutions VI. Water Wave Experiments and Observations VII. Future Directions VIII. Bibliography Glossary Deep water A surface wave is said to be in deep water if its wavelength is much shorter than the local water depth. Internal wave A internal wave travels within the interior of a fluid. The maximum velocity and maximum amplitude occur within the fluid or at an internal boundary (interface). Internal waves depend on the density-stratification of the fluid. Shallow water A surface wave is said to be in shallow water if its wavelength is much larger than the local water depth. Shallow water waves Shallow water waves correspond to the flow at the free surface of a body of shallow water under the force of gravity, or to the flow below a horizontal pressure surface in a fluid. Shallow water wave equations Shallow water wave equations are a set of partial differential equations that describe shallow water waves. 1 Solitary wave A solitary wave is a localized gravity wave that maintains its coherence and, hence, its visi- bility through properties of nonlinear hydrodynamics. Solitary waves have finite amplitude and propagate with constant speed and constant shape. Soliton Solitons are solitary waves that have an elastic scattering property: they retain their shape and speed after colliding with each other.
    [Show full text]
  • Dispersion Coefficients for Coastal Regions
    (!1vL- L/0;z7 /- NUREG/CR-3149 PNL-4627 Dispersion Coefficients for Coastal Regions Prepared by B. L. MacRae, R. J. Kaleel, D. L. Shearer/ TRC TRC Environmental Consultants, Inc. Pacific Northwest Laboratory Operated by Battelle Memorial Institute Prepared for U.S. Nuclear Regulatory Commission NOTICE This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, or any of their employees, makes any warranty, expressed or implied, or assumes any legal liability of re ­ sponsibility for any third party's use, or the results of such use, of any information, apparatus, product or process disclosed in this report, or represents that its use by such third party would not infringe privately owned rights. Availability of Reference Materials Cited in NRC Publications Most documents cited in NRC publications will be available from one of the following sources: 1. The NRC Public Document Room, 1717 H Street, N.W. Washington, DC 20555 2 The NRC/GPO Sales Program, U.S. Nuclear Regulatory Commission, Washington, DC 20555 3. The National Technical Information Service, Springfield, VA 22161 Although the listing that follows represents the majority of documents cited in NRC publications, it is not intended to be exhaustive. Referenced documents available for inspection and copying for a fee from the NRC Public Docu­ ment Room include NRC correspondence and ir.ternal NRC memoranda; NRC Office of Inspection and Enforcement bulletins, circulars, information not1ces, inspection and investigation notices; licensee Event Reports; vendor reports and correspondence; Commission papers; and applicant and licensee documents and correspondence.
    [Show full text]
  • 2.016 Hydrodynamics Free Surface Water Waves
    2.016 Hydrodynamics Reading #7 2.016 Hydrodynamics Prof. A.H. Techet Fall 2005 Free Surface Water Waves I. Problem setup 1. Free surface water wave problem. In order to determine an exact equation for the problem of free surface gravity waves we will assume potential theory (ideal flow) and ignore the effects of viscosity. Waves in the ocean are not typically uni-directional, but often approach structures from many directions. This complicates the problem of free surface wave analysis, but can be overcome through a series of assumptions. To setup the exact solution to the free surface gravity wave problem we first specify our unknowns: • Velocity Field: V (x, y, z,t ) = ∇φ(x, y, z,t ) • Free surface elevation: η(x, yt, ) • Pressure field: p(xy,,z , t ) version 3.0 updated 8/30/2005 -1- ©2005 A. Techet 2.016 Hydrodynamics Reading #7 Next we need to set up the equations and conditions that govern the problem: • Continuity (Conservation of Mass): ∇=2φ 0 for z <η (Laplace’s Equation) (7.1) • Bernoulli’s Equation (given some φ ): 2 ∂φ 1 pp− a ∂t +∇2 φ + ρ + gz = 0 for z <η (7.2) • No disturbance far away: ∂φ ∂t , ∇→φ 0 and p = pa − ρ gz (7.3) Finally we need to dictate the boundary conditions at the free surface, seafloor and on any body in the water: (1) Pressure is constant across the free surface interface: p = patm on z =η . ⎧∂φ 1 2 ⎫ p =−ρ ⎨ − V − gz⎬ + c ()t= patm . (7.4) ⎩ ∂t 2 ⎭ Choosing a suitable integration constant, ct( ) = patm , the boundary condition on z =η becomes ∂φ 1 ρ{ + V 2 + gη} = 0.
    [Show full text]
  • Solitons Et Dispersion
    Solitons et dispersion Raphaël Côte Table des matières Table des matières i Introduction iii Preamble: Orbital stability of solitons vii Solitons . vii Orbital stability of solitons . ix 1 Multi-solitons 1 1 Construction of multi-solitons . 1 2 Families of multi-solitons . 8 2 Soliton resolution for wave type equations 11 1 Profile decomposition . 12 2 Linear energy outside the light cone . 14 3 The nonlinear wave equation in dimension 3 . 17 4 Equivariant wave maps to the sphere S2 ....................... 19 5 Back to the H˙ 1 × L2 critical wave equation: dimension 4 . 23 3 Blowup for the 1D semi linear wave equation 27 1 Blow up and blow up curve . 27 2 Description of the blow up . 28 3 Construction of characteristic points . 33 4 The Zakharov-Kuznetsov flow around solitons 37 1 Liouville theorems . 39 2 Asymptotic stability . 41 3 Multi-solitons . 42 5 Formation of Néel walls 45 1 A two-dimensional model for thin-film micromagnetism . 45 2 Compactness and optimality of Néel walls . 48 3 Formation of static Néel walls . 49 Index 53 Bibliography 55 References presented for the habilitation . 55 General references . 56 i Introduction objet de ce mémoire est de présenter quelques aspects de la dynamique des solutions L’ d’équations aux dérivées partielles dispersives. Le fil conducteur liant les différents travaux exposés ici est l’étude des solitons. Ce sont des ondes progressives ou stationnaires, des objets non linéaires qui préservent leur forme au cours du temps. Elles possèdent des propriétés de rigidité remarquables, et jouent un rôle prééminent dans la description de solutions génériques en temps long.
    [Show full text]