Tsunami Modeling

Total Page:16

File Type:pdf, Size:1020Kb

Tsunami Modeling Tsunami modeling Philip L-F. Liu Class of 1912 Professor School of Civil and Environmental Engineering Cornell University Ithaca, NY USA PASI 2013: Tsunamis and storm surges Valparaiso, Chile January 2-13, 2013 1/17/13 1 Mathematical models for tsunami waves 1/17/13 2 Model Types • Nonlinear Shallow Water Equations – Linear Shallow Water Equations • Boussinesq-type Equations • Computational Fluid Dynamics – RANS, SPH, LES, DNS 1/17/13 3 CFD models o Navier-Stokes equations where is the velocity vector and p the pressure. o Boundary conditions and initial conditions 1. On the free surface, zxyt = (, ,) the kinematic boundary condition Dz()∂∂∂ =++=0, or uvw Dt∂∂∂ t x y 1/17/13 4 2. On the free surface, the dynamic boundary condition requires: pp==atm 0, on z = ( xyt , , ) 3. On the seafloor, zhxyt = (, ,) , the kinematic boundary condition Dz() h∂∂∂ h h h =++=0, or uv w Dt∂∂∂ t x y If the time history of the seafloor displacement is prescribed, this boundary condition is linear. 1/17/13 5 o Initial conditions 1. If the time history of the seafloor displacement is prescribed, the initial conditions can be given as: ==0,u 0, at t =0. 2. If seafloor is stationary, the bottom boundary condition become ∂∂hh uv+== wzhxy, on ( , ) ∂∂xy The initial free surface displacement needs to be prescribed, (xyt , ,== 0)0 ( xy , ) with u( xyzt , , , == 0) 0, where 0 (, xy ) mimics the final seafloor displacement after an earthquake. 1/17/13 6 Approximated long-wave models o Long-wave assumptions o Depth-integrated equations (reduce 3D to 2D) • Boussinesq-type equations • Shallow-water equations 1/17/13 7 Brief Review on Tsunami Models • O(months-years) N-S model (3-D) HIGH Based on Navier-Stokes Equaons (or Euler Equaons for potenJal flow) Computational Cost e.g., Truchas (Los Alamos Lab, USA) Short • Boussinesq Model (H2D) O(days-months) - Based on Boussinesq-type Equaons for weakly Accuracy wave dispersive waves, e.g. FUNWAVE (Wei and Kirby, 1995; Kirby et al, 1998) CoulWave (Lyne] and Liu, 2004) • SWE Model (H2D) O(hours) based on Shallow Water Equaons for non-dispersive waves, e.g., COMCOT (Liu et al,1995), MOST (Titov et al, 1997) Anuga (Australia) Dispersion is neglected!! LOW TSUNAMI (ARSC, UAF) 1/17/13 8 Long wave approximaons Consider small amplitude single harmonic progressive wave: = Aeikx() t gkAcosh k ( z+ h ) uxzt(,,)= eikx() t , cosh kh igkAsinh k ( z+ h ) wxzt(,,)= eikx() t , cosh kh h/L>1/2 coshkz (+ h ) pg= cosh kh g Ckh== tanh kk h/L>1/2 1124 coshkz (+=+ h ) 1[ kz ( + h )] +[ kz ( + h )] + 2 24 1 3 sinhkz (+= h )[ kz ( + h )] +[ kz ( + h )] + 6 1/17/13 9 Shallow water wave equations • Assumptions: 1. Horizontal scales (i.e., wavelength) are much longer than the vertical scale (i.e., the water depth). Accurate only for very long waves, kh < 0.25 (wavelength > 25 water depths) 2. Ignore the viscous effects. • Consequences: 1. The vertical velocity, w, is much smaller than the horizontal velocity, (u,v). 2. The pressure is hydrostatic, pg = () z. 3. The leading horizontal velocity is uniform in water depth. 1/17/13 10 • Nonlinear shallow-water wave equations Conservation of mass: ∂∂ ∂ [++hhuhv] ( +) +( +) =0 ∂∂tx ∂ y Momentum equations: ∂u +uug +g =0 t ∂ • Linear shallow–water wave equations ∂ [ ++hh] g[u ] =0 t ∂ ∂u +g =0 ∂t 1/17/13 11 Special case: One-dimensional and constant depth ∂∂ ++=( hu) 0, ∂∂tx ∂∂∂uu ++ug =0. ∂∂txx ∂ which can be rewritten as ∂∂ ++(uC) ( u +20, C) = ∂∂tx ∂∂ +(uC) ( u 2 C) = 0. tx ∂∂ where . Cg=+( h) dx uC±=± are Riemann invariants along charateristics along u2 C . dt 1/17/13 12 Linear shallow water wave equation ∂∂ u +=h 0, ∂∂tx ∂∂u +=g 0. ∂∂tx which can be combined into a wave equaon ∂∂22 =gh 0 ∂∂tx22 The general soluJon is =±(),x Ct C = gh Wave form does not change! NO dispersion! 1/17/13 13 Boussinesq equations (Peregrine, 1967; Ngowu, 1993) o Pressure field is not hydrostatic: quadratic o Horizontal velocities are not vertically uniform 2 uxzt(,,)=++ Axt (,) zBxt (,) z Cxt (,) should be small compared to A(x,t) z Accurate for long and x u(x,z,t) intermediate depth waves, kh < 3 (wavelength > 2 water depths) Functions B, C lead to 3rd order spatial derivatives in model 1/17/13 14 High-Order Boussinesq Equations (Gobbi et al., 2000) 2 uxzt(,,)=+ Axt (,) z * Bxt (,) + z * Cxt (,) 34 ++zDxtzExt* ( , ) * ( , ) Should be small compared to B,C group z - Accurate for long, x intermediate, and moderately deep waves, kh < 6 (wavelength > 1 water depth) - Functions D, E lead to 5th order spatial derivatives in model 1/17/13 15 Derivation of depth-integrated wave equations . Normalize the conservation equations and the associated boundary conditions . Perturbation analysis h ;1 = 0 = L0 . Integrate the governing equations and employ o the irrotational condition o the boundary conditions 1/17/13 16 Boussinesq-type equations 2 OO()= 1,( )= 1 1/17/13 17 Boussinesq equaons OOkh()== (2 ); 11∂∂h +++(hhO)uu + 23 =( ) (,, 422 ) ∂∂tt 3 2 ∂∂u 2422h +++uu ( u) hhO( u) =(, ,) ∂∂tt2 Different Boussinesq equaon forms can be derived by using different Representave velocity or by subsJtuJng the higher order terms by an Order one relaon, e.g., ∂u = +O(,2 ) ∂t The resulJng equaons have slightly different dispersion characterisJcs. 1/17/13 18 Modeling the dispersion effects • To apply the depth-integrated equaons for shorter waves, higher order terms can be included. However, it more difficult to solve the high-order model – Momentum equaon: ∂∂uu ∂5 u +++=uC.... 0 ∂∂tx1 ∂ x5 – To solve consistently, numerical truncaon error (Taylor series error) for leading term must be less important than included terms. • For example: 2nd order in space finite difference: ∂ux(,) t ux (+ xt ,) ux ( xt ,) x2 ∂3 ux (,) t oo= o o ∂xx26∂ x3 • High-order model requires use of 6-point difference formulas (Dx6 accuracy). AddiJonally, Jme integraon would require a Dt6 accurate scheme. • It is difficult to specify boundary condiJons. 1/17/13 19 COULWAVE (Cornell University Long and Intermediate Wave Modeling) A Multiple Layers Approach • Divide water column into arbitrary layers – Each layer governed by an independent velocity profile, each in the same form as tradiJonal Boussinesq models: z x 2 uz1111()=+ A BzCz + z 1 uz11()= uz 21 () 2 uz2222()=+ A BzCz + z2 uz22()= uz 32 () 2 uz3333()=+ A BzCz + 1/17/13 20 COULWAVE (Cornell University Long and Intermediate Wave Modeling) • Regardless of # of layers, highest order of derivation is 3 • The more layers used, the more accurate the model • Any # of layers can be used – 1-Layer model = Boussinesq model – Numerical applications of 2-Layer model to be discussed • Location of layers will be optimized for good agreement with known, analytic properties of water waves 2 uz1111()=+ A BzCz + 2 uz2222()=+ A BzCz + 2 uz3333()=+ A BzCz + 1/17/13 21 COULWAVE . Governing equations: multi-layer, fully-nonlinear COULWAVE Linear theory high-order, single-layer result . Numerical scheme: predictor-corrector method o Predictor: explicit 3rd-order Adams-Bashforth o Corrector: implicit 4th-order Adams-Moulton 1/17/13 22 Linear Dispersion Properties of Multi-Layer Model • Compare phase and group velocity with linear theory 1.2 Boussinesq High-Order 1.15Equaons Boussinesq Two-Layer Model (One-Layer Equaons 1.1 Four-Layer Model) linear Three-Layer Model Model 1.05 C / C 1 0.95 5 10 15 20 25 30 35 kh Boussinesq High-Order Equaons Boussinesq 1.2 Equaons 1.15 Two-Layer Model 1.1 linear Three-Layer Model Four-Layer Model 1.05 Cg / Cg 1 0.95 5 10 15 20 25 30 35 1/17/13 23 kh General remarks on tsunami modeling by long-wave equations . Add-hoc dissipation mechanism o Bo]om fricJon: quadrac form o Wave breaking: parameterized model . Controlling wave mechanism o Near the earthquake source: linear, dispersive or non-dispersive waves o Deep ocean: linear, non-dispersive waves o Close to the shore: nonlinear, non-dispersive waves (Observaon from the wave gauge data, and, the analyJcal study of one-dimensional problem) 1/17/13 24 Tsunami modeling packages . Boussinesq-type equaons model o COULWAVE (Cornell University Long and Intermediate Wave Modeling Package) o An improved mulJ-layer approach . Shallow water equaons model o COMCOT (Cornell Mul-grid Coupled Tsunami Model) o Covers tsunami generaon, propagaon, and wave runup/rundown on coastal regions o If desired, physical frequency dispersion can be mimicked by the numerical dispersion o This is a more pracJcal choice for the tsunami simulaon 1/17/13 25 Numerical simulaon of tsunami waves . Governing equaons (COMCOT) o Linear model: deep ocean o Nonlinear model: shallower region 1/17/13 26 COMCOT: numerical scheme . Explicit leap-frog Finite Differencing Method center of sub-level grid . Nested grid center of parent grid volume flux of sub-level grid connecJng boundary volume flux of parent grid 1/17/13 27 Inputs (inial condions) for the tsunami wave modeling . Earthquake-generated seafloor displacement o Impulsive model • The seafloor deforms instantaneously and the enJre fault line ruptures simultaneously • The sea surface follows the seafloor deformaon instantaneously o Transient model • The seafloor deformaon and the rupture along the fault line are both described as transient processes • Require Jme histories of the horizontal and verJcal seafloor displacements 1/17/13 28 Fault plane models . Impulsive vs. transient models o Transient models are not always available o The speed of the rupture is orders of magnitude faster than the tsunami wave propagaon speed o Advanced impulsive models can provide detailed spaal variaons of the final seafloor displacement o No significant difference in the resulJng tsunami wave height and propagaon Jme from these two types of fault plane models • Analysis arguments: Kajiura (1963); Momoi (1964-5); Tuck and Hwang (1972) • Numerical experiments: Wang and Liu (2006) 1/17/13 29 Numerical simulaons of tsunami waves .
Recommended publications
  • Well-Balanced Schemes for the Shallow Water Equations with Coriolis Forces
    Well-balanced schemes for the shallow water equations with Coriolis forces Alina Chertock, Michael Dudzinski, Alexander Kurganov & Mária Lukáčová- Medvid’ová Numerische Mathematik ISSN 0029-599X Volume 138 Number 4 Numer. Math. (2018) 138:939-973 DOI 10.1007/s00211-017-0928-0 1 23 Your article is protected by copyright and all rights are held exclusively by Springer- Verlag GmbH Germany, part of Springer Nature. This e-offprint is for personal use only and shall not be self-archived in electronic repositories. If you wish to self-archive your article, please use the accepted manuscript version for posting on your own website. You may further deposit the accepted manuscript version in any repository, provided it is only made publicly available 12 months after official publication or later and provided acknowledgement is given to the original source of publication and a link is inserted to the published article on Springer's website. The link must be accompanied by the following text: "The final publication is available at link.springer.com”. 1 23 Author's personal copy Numer. Math. (2018) 138:939–973 Numerische https://doi.org/10.1007/s00211-017-0928-0 Mathematik Well-balanced schemes for the shallow water equations with Coriolis forces Alina Chertock1 · Michael Dudzinski2 · Alexander Kurganov3,4 · Mária Lukáˇcová-Medvid’ová5 Received: 28 April 2014 / Revised: 19 September 2017 / Published online: 2 December 2017 © Springer-Verlag GmbH Germany, part of Springer Nature 2017 Abstract In the present paper we study shallow water equations with bottom topog- raphy and Coriolis forces. The latter yield non-local potential operators that need to be taken into account in order to derive a well-balanced numerical scheme.
    [Show full text]
  • 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]
  • Part II-1 Water Wave Mechanics
    Chapter 1 EM 1110-2-1100 WATER WAVE MECHANICS (Part II) 1 August 2008 (Change 2) Table of Contents Page II-1-1. Introduction ............................................................II-1-1 II-1-2. Regular Waves .........................................................II-1-3 a. Introduction ...........................................................II-1-3 b. Definition of wave parameters .............................................II-1-4 c. Linear wave theory ......................................................II-1-5 (1) Introduction .......................................................II-1-5 (2) Wave celerity, length, and period.......................................II-1-6 (3) The sinusoidal wave profile...........................................II-1-9 (4) Some useful functions ...............................................II-1-9 (5) Local fluid velocities and accelerations .................................II-1-12 (6) Water particle displacements .........................................II-1-13 (7) Subsurface pressure ................................................II-1-21 (8) Group velocity ....................................................II-1-22 (9) Wave energy and power.............................................II-1-26 (10)Summary of linear wave theory.......................................II-1-29 d. Nonlinear wave theories .................................................II-1-30 (1) Introduction ......................................................II-1-30 (2) Stokes finite-amplitude wave theory ...................................II-1-32
    [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]
  • Fluid Simulation Using Shallow Water Equation
    Fluid Simulation using Shallow Water Equation Dhruv Kore Itika Gupta Undergraduate Student Graduate Student University of Illinois at Chicago University of Illinois at Chicago [email protected] [email protected] 847-345-9745 312-478-0764 ABSTRACT rivers and channel. As the name suggest, the main In this paper, we present a technique, which shows how characteristic of shallow Water flows is that the vertical waves once generated from a small drop continue to ripple. dimension is much smaller as compared to the horizontal Waves interact with each other and on collision change the dimension. Naïve Stroke Equation defines the motion of form and direction. Once the waves strike the boundary, fluids and from these equations we derive SWE. they return with the same speed and in sometime, depending on the delay, you can see continuous ripples in In Next section we talk about the related works under the the surface. We use shallow water equation to achieve the heading Literature Review. Then we will explain the desired output. framework and other concepts necessary to understand the SWE and wave equation. In the section followed by it, we Author Keywords show the results achieved using our implementation. Then Naïve Stroke Equation; Shallow Water Equation; Fluid finally we talk about conclusion and future work. Simulation; WebGL; Quadratic function. LITERATURE REVIEW INTRODUCTION In [2] author in detail explains the Shallow water equation In early years of fluid simulation, procedural surface with its derivation. Since then a lot of authors has used generation was used to represent waves as presented by SWE to present the formation of waves in water.
    [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]
  • 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]
  • Shallow-Water Equations and Related Topics
    CHAPTER 1 Shallow-Water Equations and Related Topics Didier Bresch UMR 5127 CNRS, LAMA, Universite´ de Savoie, 73376 Le Bourget-du-Lac, France Contents 1. Preface .................................................... 3 2. Introduction ................................................. 4 3. A friction shallow-water system ...................................... 5 3.1. Conservation of potential vorticity .................................. 5 3.2. The inviscid shallow-water equations ................................. 7 3.3. LERAY solutions ........................................... 11 3.3.1. A new mathematical entropy: The BD entropy ........................ 12 3.3.2. Weak solutions with drag terms ................................ 16 3.3.3. Forgetting drag terms – Stability ............................... 19 3.3.4. Bounded domains ....................................... 21 3.4. Strong solutions ............................................ 23 3.5. Other viscous terms in the literature ................................. 23 3.6. Low Froude number limits ...................................... 25 3.6.1. The quasi-geostrophic model ................................. 25 3.6.2. The lake equations ...................................... 34 3.7. An interesting open problem: Open sea boundary conditions .................... 41 3.8. Multi-level and multi-layers models ................................. 43 3.9. Friction shallow-water equations derivation ............................. 44 3.9.1. Formal derivation ....................................... 44 3.10.
    [Show full text]
  • Waves in Water from Ripples to Tsunamis and to Rogue Waves
    Waves in Water from ripples to tsunamis and to rogue waves Vera Mikyoung Hur (Mathematics) (Image from Flickr) The motion of a fluid can be very complicated as we know whenever we see waves break on a beach, (Image from the Internet) The motion of a fluid can be very complicated as we know whenever we fly in an airplane, (Image from the Internet) The motion of a fluid can be very complicated as we know whenever we look at a lake on a windy day. (Image from the Internet) Euler in the 1750s proposed a mathematical model of an incompressible fluid. @u + (u · r)u + rP = F, @t r · u = 0. Here, u(x; t) is the velocity of the fluid at the point x and time t, P(x; t) is the pressure, and F(x; t) is an outside force. The Navier-Stokes equations (adding ν∆u) allow the fluid to be viscous. It is concise and captures the essence of fluid behavior. The theory of fluids has provided source and inspiration to many branches of mathematics, e.g. Cauchy's complex function theory. Difficulties of understanding fluids are profound. e.g. the global-in-time solution of the Navier-Stokes equations in 3 dimensions is a Clay Millennium Problem! Waves, jets, drops come to mind when thinking of fluids. They involve one or more fluids separated by an unknown surface. In the mathematical community, they go by free boundary problems. (Image from the Internet) Free boundaries are mathematically challenging in their own right. They occur in many other situations, such as melting of ice stretching a membrane over an obstacle.
    [Show full text]
  • Modified Shallow Water Equations for Significantly Varying Seabeds Denys Dutykh, Didier Clamond
    Modified Shallow Water Equations for significantly varying seabeds Denys Dutykh, Didier Clamond To cite this version: Denys Dutykh, Didier Clamond. Modified Shallow Water Equations for significantly varying seabeds: Modified Shallow Water Equations. Applied Mathematical Modelling, Elsevier, 2016, 40 (23-24), pp.9767-9787. 10.1016/j.apm.2016.06.033. hal-00675209v6 HAL Id: hal-00675209 https://hal.archives-ouvertes.fr/hal-00675209v6 Submitted on 26 Apr 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Distributed under a Creative Commons Attribution - NonCommercial| 4.0 International License Denys Dutykh CNRS, Université Savoie Mont Blanc, France Didier Clamond Université de Nice – Sophia Antipolis, France Modified Shallow Water Equations for significantly varying seabeds arXiv.org / hal Last modified: April 26, 2016 Modified Shallow Water Equations for significantly varying seabeds Denys Dutykh∗ and Didier Clamond Abstract. In the present study, we propose a modified version of the Nonlinear Shal- low Water Equations (Saint-Venant or NSWE) for irrotational surface waves in the case when the bottom undergoes some significant variations in space and time. The model is derived from a variational principle by choosing an appropriate shallow water ansatz and imposing some constraints.
    [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]