Edinburgh Meeting Notes 1

These are notes from a meeting entitled Advanced Numerical Studies in Nonlinear PDEs in Edinburgh, Scotland.

Walter Craig (McMaster): Water Wave Interactions

I’m an analyst but I’m going to talk about numerics and experiments as well as analysis. We will discuss the problem of water waves and then I’ll talk about two specific settings in which the theory has led to good and quite elegant numerics and the numerics have started to answer some questions.

(joint work with P. Guyenne and C. Sulem)


  • Free surface water waves
  • Hamiltonian PDEs
  • Periodic Traveling wave patterns
  • Solitary wave Interactions
  • The KdV scaling limit

Free surface water waves

Euler’s equations of hydrodynamics, incompressible and irrotational flow. This is therefore given as a potential flow. The irrotational assumption is really an oceanographers assumption. Of course, there is vorticity but we follow the models of oceanographers.

The fluid domain is $-h < y < \eta (x,t)$. So, the domain is changing. Free surface boundary conditions hold on $y = \eta (x,t)$.

Zakharov’s Hamiltonian

  • The energy functional
$$ H = K+P $$ $$ K = \int_x \int_{-h}^{\eta(x)} \frac{1}{2} |\nabla \phi|^2 dy dx. $$ $$ P = \int_x \frac{g}{2} \eta^2 dx. $$

This could also include surface tension effects.

  • Zakharov’s choice of variables, $$ z = ( \eta(x), \xi(x) = \phi(x, \eta(x))), $$ for which we consider $\phi = \phi[\eta, \xi] (x,y)$.
  • Express the energy in terms of $\xi$ and $\eta$. This involves the Dirichlet-Neumann operator $G(\eta)$.

Dirichlet-Neumann operator

  • Laplace’s equation on the fluid domain: $\Delta \phi = 0$ subject to bottom Neumann boundary condition. Free surface boundary data $\phi (x, \eta(x)) = \xi(x)$, for which the D-N operator is given by $$ \xi(x) \longmapsto \phi(x,y) \longmapsto N \cdot \nabla \phi (1+ |\nabla_x \eta|^2)^{1/2} := G(\eta) \xi(x). $$
  • In these coordinates, we can rewrite the boundary conditions in a new (and nicer) form. This reexpresses the water wave problem as a Hamiltonian system in Darboux coordinates.

Hamiltonian PDEs

  • KdV is a Hamiltonian PDE with a different symplectic structure.
  • Other Hamiltonian PDEs
    • shallow water equations
    • Boussinesq
    • KP
    • NLS
    • Dysthe equation
Many of these problems arise in scaling limits of the water wave problem.

Lemma (Properties of D-N operator):

  1. $G(\eta) \geq 0$ and $G(\eta) 1 = 0$.
  2. $G(\eta)^* = G(\eta)$ Hermitian Symmetric
  3. $G(\eta): H^1_\xi \rightarrow L^2_\xi$ is analytic in $\eta$ for $\eta \in C^1$. There is an operator valued power series expansion of $G(\eta)$ (using a theorem of Christ-Journé 1987).
  4. Some explicit calculations of the Taylor expansion (I couldn’t keep up….)
  5. Conservation Laws
    • Mass: $M = \int \eta dx$ (He shows the calculation using properties of $G$.)
    • Momentum: Similar calculation
    • Energy: Easy since the commutator of $H$ with itself vanishes.
  6. Taylor expansion of the Hamiltonian
  7. Linearized equations; comparison with the harmonic oscillator.

Periodic Traveling wave patterns

  • Can I find traveling wave solutions? $$ \eta(x,t) = \eta ( x-tc); \xi(x,t) = \xi(x - tc) $$
  • Spatially periodic, $\Gamma \subset {\mathbb{R}^{d-1}}$. $$ \eta(x + \gamma, \cdot) = \eta(x, \cdot), \xi(x+\gamma, \cdot) = \xi(x, \cdot), ~ \forall \gamma \in \Gamma. $$
On such domains, we can use the Fourier tranform.

Rk: Notice this is a mathematician imposing a period rather than the physics making that selection. More can be said in this direction, but let’s proceed this way.

Rk: These (time independent) traveling wave patterns can be imagined to emerge in transient interactions in seas. The nonlinear actions create large amplitudes and this might be related to the phenomena of freak waves.

Equations for traveling waves.

Periodic traveling wave patterns are critical points of the Hamltonian on the variety $I = const$, with Lagrange multiplier $c \in {\mathbb{R}^{d-1}}.$

This leads to a bifurcation problem.

brief history (dimension $d=2$):

  • Levi-Civita 1925; existence of traveling waves
  • Struik 1926; traveling waves case
  • Zeidler 1971
  • Beale 1979
  • Jones-Toland 1985
brief history (dimension $d=3$):
  • Reeder Shinbrot 1981
  • Sun 1986
  • Craig-Nicholls 2000
  • Iooss-Plotnikov-Toland 2000 (small divisor problem)
He shows a picture from the wave tank at Penn State. He then shows some numerics which are trying to model those observations and they look beautiful.

Kuksin Question: Stability of these patterns?

Craig Answer: This is a very good question. I don’t know results like that. This is related to Benjamin-Feir. McLean showed instability for $d = 3$. Some further discussion….We need the Bloch theory of stability for these wave patterns. This appears to be difficult analytically so might need some numerical studies at first. There are instability zones….

Solitary wave Interactions

Solitary waves in 2-dimensions (Friedreichs-Hyers 1954, Amick-Fraenkel-toland 1980s)
  • Head-on collisions of solitons.
The numerics reveal some inelasticity in the collision. We’d like to understand those. If we make the amplitude of the solitons bigger, the dispersive ripples are more visible.

The KdV scaling limit

Titi’s Question: Can we reduce to the surface equations including rotation?
Craig’s Answer: Yes and No. You can make a rotation depending purely on y and impose that. Then it is reducible. But this is rather artificial. There is stuff that happens in the middle which is not a surface effect. Therefore, this problem requires a more complete analysis of the Euler equation and will not collapse to a system on the surface.

Sergey Nazarenko (Warwick): Assumptions, Techniques, Cahllenges in Wave Turbulence

This is not so much about new result. Instead, this is an attempt by a physicist trying to explain wave turbulence ideas being explored by physicists to mathematicians. My view is that there is a lot of interesting work to be done. Lots of open problems….

What is wave turublence?

He shows a picture of a relatively calm seashore from Nice. He emphasizes there is a wide range fo scales in these problems. WT is a statistical system of nonlinear waves.


  • Water waves
  • Waves in rotating and stratified fluids (internal and inertial waves, Rossby waves)
  • Plasma waves
  • Waves in Bose-Einstein condensates
  • Kelvin waves on quantized vortex filnments
  • MHD turbulence in interstellar turbulence and solar wind
  • Nonlinear optics
  • Solids: phonons, spin waves. Kinetics of phonons in weakly anharmonic crystals is a first example of study in tis direction (1920s). I didn’t catch the name….
He shows a picutre of a wave take of Lukaschuk.

Waves in fusion plasmas. Shows a picture of a Tokamak. Drift wave turbulence causes anomalous heat and particle loss - major problem for fusion. The devices have grown larger and larger basically to carry out the confinement for a longer period of time.

MHD turbulence in astrophysics. He shows some data from the Ulysses/Swoops (los alamos) solar wind studies.

Bose Einstein Condensates Nazarenko-Onorato 2006:

  • Inverse cascade - condensation
  • Condensate strongly affects WT
Quantum Turbulence (see Lvov et. al. 2007) (Superfluid turbulence)
  • Kelvin waves on quantized vortex filaments
  • Interaction with hydro eddies (vortex bundles) is important
  • Kelvin Wave Turbulence
Optical Turbulence
  • Bortolozzo et. al 2008
  • This project studies nonlinear corrections (coming from the optical physics) which are included beyond the 1d NLS model.
Kuksin Question: Which corrections? Can you write them down?
Nazareknko: Something like a DNLS correction…not so clear.

Ingredients in the approach

He writes $NLS_3^\pm (T^d)$ and comments that this is a physically reasonable model but we are really interested in the study in infinite space with finite energy density.

He reexpresses the NLS equation in Fourier language.

Set of wave modes: amplitudes and phases.

N-mode joint probability density function. Some notation….probability…sectors in the wave modes setting.

Random Phase (RP) and Random Phase Amplitude (RPA) systems


All phases are independent random variables such that uniformly distributed on $S^1$.


  1. All amplitudes and all phases are independent random variables.
  2. All phases are uniformly distributed on $S^1$.
Note: RPA does not mean Gaussian. Nevertheless, we have obtained successful closures without assuming the Gaussian statistics.

Frog Jumps!

  • expanding in small nonlinearity
  • Assuming RP at $t=0$.
  • Taking limit of a large box followed by the limit of small nonlinearity.
(The order of these steps is important.)

Evolution of joint PDF? We can derive the evolution equation under these assumptions. The derivation is rather systematic, in fact it is perhaps rigorous.

Mathematical Challenges:

  • WT is formally derived for $t=0$.
  • Does it work at the long time of nonlinear evolution?
  • Does RPA survive over this time?
  • Adding forcing and dissipation: will WT describe the steady state?
Hmmmm….This RPA condition at $t=0$ reminds me a bit of the assumption of product wave function in the QMB theory. The dynamics in the Hartree derivation might drive the multiparticle wave function away from the product case. Here we have a dynamic that might drive us away from the RPA condition.
Evolution of 1-mode PDF.

Kinetic equation (Hasselmann 1962).

Kolmogorov-Zakharov state.

  • Explained a steady state spectrum corresponding to energy cascade.
  • Exact solution of the asymptotic closure.
Numerics and Analysis of KE.
  • What is the role of KZ solutions with respect to the thermodynamic Rayleigh-Jeans state?
  • Similar issues for the classical Boltzmann equation.
Zakharov was awarded the 2003 Dirac Medal for “putting the theory of wave turbulence on a firm mathematical ground”! What is it that we want to do?

Gregor Tanner (Nottingham): A wave chaos approach towards describing the vibro-acoustic response of engineering structures

(Joint work with D. Chappel, Stefano Gianai, Hanya Ben Hamdin, Dmitrii Maksimov)

This talk is more directed toward engineering applications. inuTech is an industrial collaborator.


  • Introduction - the need for numerical short wavelenght methods in vibroacoustics
  • From wave equations to the Liouville equation
  • Solving the Liouville equation - a boundary integral approach (Dynamical Energy Analysis - DEA)
  • Tackling the Midfrequency problem - hybrid methods
  • Numerical results
Aim: predicting wave intensity distributions for the vibro-acoustical response of mechanical structures. Think of a car. Companies like Bombardier and Airbus use these methods. It is a difficult problem. You want these structures to be quite and with no noise in the interior.

Where is the problem?


  • Low frequencies - wavelength around the size of the object
  • Finite Element method
  • Boundary elemtn method
  • plane wave methods
High Frequencies:
  • Ray tracing
  • Statistical energy analysis
Midfrequency problems:
  • Structures with large variations in the local wavelength. (Large variations in the stiffness of components, ie body frame and side panels.)
  • Hybrid methods. Try to connect exact numerical methods with the statistical methods.

Short wavelength approximations - from wave chaos to statistical methods

  • Wave chaos -short wavelength asymptotics
    • Keller
    • Gutzwiller
    • Berry
    • Bogolmolny
    • Smilansky
  • Nonlinear dynamics - thermodynamic formalism
    • Ruelle
    • Arnold
    • Sinai
    • Eckmann
    • Cvitanovic - chaosbook.org
  • Wave transport - statistical methods in vibro-acoustics
    • Lyon - SEA (1967 paper)
    • Langley - WIA
    • Heron
    • Weaver - diffusion equation
    • Le Bot - radiative transformation
Linear wave equation. WKB ansatz. Hamiltoninan equations for the amplitude and phase. Characteristics of JH; nonlinear ODE; Liouville equation (linear).

Linear wave –> WKB –> HJ equation –> Liouville Equation

Think of polygonal billiards, not necessarily convex.

We want to understand the influence of a source (transfmitting at frequency $\omega$) at one location on the wave amplitude at another point. He writes this as a green’s function $G(r, r_0, \omega)$.

Small wavelength limit, so low frequency waves.

Write things as sums over all paths.

Perron-Frobenius operator…

Tandem Satellite images of coastline of Madagascar

I started wondering about connections between these ideas and quantum ergodicity…

Typical Wave Function in a stadium Billiard

Bouncing Ball Modes


I had a nice conversation with Gregor after the break. I learned from him about microlasers. The idea is to build a circular region out of a lasing material. We energize the material somehow with hopes to excite the whispering gallery mode. The laser light propagates near the boundary but can be arranged to exit the medium by raising the curvature at a specific location. These appear to be rather hard to control to create a unidirectional beam. Since the losses take place all along the boundary, there is very little power in the output beam. Some web searching revealed an advance made by the Capasso group at Harvard. Elliptic Notched Microlaser Cavity Drawing Elliptic notched microlaser cavity SEM photograph. Microlaser Cavity (Artistic Rendition) Schematic Image Artistic Rendering

David Dritschel (St. Andrews): CLAM, The Combined Lagrangian Advection Method

(Many many collaborators)

Courbet's "The Wave"

I’ll be speaking a bit about a numerical method. I’ll focus mostly on the results we’ve obtained to understand the large scale atmospheres, like Jupiter and perhaps also the ocean.

The numerical method (CLAM) emerges from a Lagrangian method from the 50s for studying fluid dynamics. Zabusky then built from these developments to develop new methods in plasmas. We’ve been extending these ideas to treat certain geophysical fluid flows.

The atmosphere and the oceans are extremely complex, turbulent flows. Accurate computer simulation is immensely difficult to achieve. However, much of this difficulty is inherent in the computational methods employed:

  • None take direct advantage of the natural inherent Lagrangian advection of dynamical, chemical and biological tracers. (Exploit Lagrangian Descriptions.)
  • None seek to separate slow vortical (eddying) and fast wave-like motions and use appropriate, optimal numerical methods for each. (Slow Rossby waves interacting with fast waves inertial-gravity waves.)
We can build the mathematical theory of the separation into the numerical methods and this will lead to better predictions.

Contour Advection (CASL) Dritschel & Ambaum 1997

geostropic and hydrostatic balances are basic features for describing atmospheric wave dynamics.

This talk reminds me somehow of Bourgain’s high/low method for proving low regularity GWP.

The idea is to use the advection of the vorticity to resolve some (especially relevant) sub-grid scales.

Dugald Duncan: IDE equation

  • full IDE equation and how it looks like, where it arises
  • Linear part of the IDEbehaviour and approximation
  • the full problem - behanviour and approximation
  • examples
$$ u_t = \sigma \int_\Omega J(x-y) [u(y,t) - u(x,t)]dy + f(u) dx ~ \forall x \in \Omega, t>0. $$ Typically, $f(u) = u - u^3$. This should be contrasted with the Allen-Cahn equation

$$ u_t = \sigma \Delta u + f(u) dx ~ \forall x \in \Omega, t>0. $$

There are no spatial derivatives. Therefore, there are now boundary conditions. Instead, this is some kind of integral dynamical equation. It is similar to the Cahn-Allen equation.

This equation is also related to sandpiles, neurons, phase transitions.

Other variations recently: Rossi, Perez-Llanos, Andreu, Mazon, Toledo et. al. They study a nonlocal version of the $p$-laplacian.

Linear IDE:

  • Ignore the nonlinear reaction term for now and take $\sigma \geq 0$ and $\Omega \subset {\mathbb{R}}$: $$ u_t = Lu.$$
  • L is a linear operator - partly a convolution: $$ Lu = \int_\Omega J(x-y) [u(y,t) - u(x,t)]dy = J * u $$ …ack slide changed….
Discontinuities don’t move. The solution collapses to the average value. Ther eis acomparison principle.

Snapshots of linear behavior.

He does a Fourier analysis of the behavior of plane waves. Instead of having an $\omega^2$, we have $$\hat{J} (\omega) - \hat{J} (0) \thicksim \frac{\omega^2}{2} \frac{d^2}{d \omega^2} {\widehat{J}} (\omega).$$

Peter Bates and Paul Fife did some of the earliest analysis on this equation.