Hunter–Saxton equation

In mathematical physics, the Hunter–Saxton equationHunter & Saxton 1991] :$(u\_t\; +\; u\; u\_x)\_x\; =\; frac\{1\}\{2\}\; ,\; u\_x^2$

is an integrable PDE that arises in the theoretical study of nematic liquid crystals. If the molecules in the liquid crystal are initially all aligned, and some of them are then wiggled slightly, this disturbance in orientation will propagate through the crystal, and the Hunter–Saxton equation describes certain aspects of such orientation waves.

Physical background

In the models for liquid crystals considered here, it is assumed that there is no fluid flow, so that only the "orientation" of the molecules is of interest.Within the elastic continuum theory, the orientation is described by a field of unit vectors n("x","y","z","t"). For nematic liquid crystals, there is no difference between orienting a molecule in the n direction or in the −n direction, and the vector field n is then called a "director field".The potential energy density of a director field is usually assumed to be given by the Oseen–Frank energy functional [de Gennes & Prost 1994 (Ch. 3)]

:

where the positive coefficients $alpha$, , $gamma$ are known as the elastic coefficients of splay, twist, and bend, respectively. The kinetic energy is often neglected because of the high viscosity of liquid crystals.

Derivation of the Hunter–Saxton equation

Hunter and Saxton investigated the case when viscous damping is ignored and a kinetic energy term is included in the model. Then the governing equations for the dynamics of the director field are the Euler–Lagrange equations for the Lagrangian

:$mathcal\{L\}\; =frac\{1\}\{2\}\; left|\; frac\{partialmathbf\{n\{partial\; t\}\; ight|^2-\; W(mathbf\{n\},\; ablamathbf\{n\})-\; frac\{lambda\}\{2\}\; (1-|mathbf\{n\}|^2),$

where $lambda$ is a Lagrange multiplier corresponding to the constraint |n|=1.They restricted their attention to "splay waves" where the director field takes the special form

:$mathbf\{n\}(x,y,z,t)\; =\; (cosvarphi(x,t),\; sinvarphi(x,t),\; 0).$

This assumption reduces the Lagrangian to

:$mathcal\{L\}\; =\; frac\{1\}\{2\}\; left(\; varphi\_t^2\; -\; a^2(varphi)\; varphi\_x^2\; ight),qquada(varphi)\; :=\; sqrt\{alpha\; sin^2\; varphi\; +\; gamma\; cos^2\; varphi\},$

and then the Euler–Lagrange equation for the angle φ becomes

:$varphi\_\{tt\}\; =\; a(varphi)\; [a(varphi)\; varphi\_x]\; \_x.$

There are trivial constant solutions φ=φ₀corresponding to states where the molecules in the liquid crystal areperfectly aligned.Linearization around such an equilibrium leads to the linear wave equationwhich allows wave propagation in both directions with speed$a\_0\; :=\; a(varphi\_0)$,so the nonlinear equation can be expected to behave similarly.In order to study right-moving waves for large "t",one looks for asymptotic solutions of the form

:$varphi(x,t;epsilon)\; =\; varphi\_0\; +\; epsilon\; varphi\_1(\; heta,\; au)\; +\; O(epsilon^2),$

where

:$heta\; :=\; x-a\_0\; t,\; qquad\; au\; :=\; epsilon\; t.$

Inserting this into the equation, one finds at the order $epsilon^2$ that

:$(varphi\_\{1\; au\}\; +\; a\text{'}(varphi\_0)\; varphi\_1\; varphi\_\{1\; heta\})\_\{\; heta\}\; =\; frac\{1\}\{2\}\; a\text{'}(varphi\_0)\; varphi\_\{1\; heta\}^2.$

A simple renaming and rescaling of the variables(assuming that $a\text{'}(varphi\_0)\; eq\; 0$)transforms this into the Hunter–Saxton equation.

Generalization

The analysis was later generalized by Alì and HunterAlì & Hunter 2006] , who allowed the director field to point in any direction, but with the spatial dependence still only in the "x" direction:

:$mathbf\{n\}(x,y,z,t)\; =\; (cosvarphi(x,t),\; sinvarphi(x,t)\; cospsi(x,t),\; sinvarphi(x,t)\; sinpsi(x,t)).$

Then the Lagrangian is

:

The corresponding Euler–Lagrange equations are coupled nonlinear wave equations for the angles φ and ψ, with φ corresponding to "splay waves" and ψ to "twist waves". The previous Hunter–Saxton case (pure splay waves) is recovered by taking ψ constant, but one can also consider coupled splay-twist waves where both φ and ψ vary. Asymptotic expansions similar to that above lead to a system of equations, which, after renaming and rescaling the variables, takes the form

:$(v\_t\; +\; u\; v\_x)\_x\; =\; 0,\; qquad\; u\_\{xx\}\; =\; v\_x^2,$

where "u" is related to φ and "v" to ψ.This system implies [Differentiate the second equation with respect to "t", substitute "v"_"xt" from the first equation, and eliminate "v" using the second equation again.] that "u" satisfies

:$left\; [\; (u\_t\; +\; u\; u\_x)\_x\; -\; frac\{1\}\{2\}\; ,\; u\_x^2\; ight]\; \_x\; =\; 0,$

so (rather remarkably) the Hunter–Saxton equation arises in this context too, but in a different way.

Variational structures and integrability

The integrability of the Hunter–Saxton equation, or, more precisely, that of its "x" derivative

:$(u\_t\; +\; u\; u\_x)\_\{xx\}\; =\; u\_x\; u\_\{xx\},$

was shown by Hunter and ZhengHunter & Zheng 1994] , who exploited that this equation is obtained from the Camassa–Holm equation

:$u\_t\; -\; u\_\{xxt\}\; +\; 3\; u\; u\_x\; =\; 2\; u\_x\; u\_\{xx\}\; +\; u\; u\_\{xxx\}$

in the "high frequency limit"

:$(x,t)\; mapsto\; (epsilon\; x,\; epsilon\; t),\; qquad\; epsilon\; o\; 0.$

Applying this limiting procedure to a Lagrangian for the Camassa–Holm equation, they obtained a Lagrangian

:$mathcal\{L\}\_2\; =\; frac\{1\}\{2\}\; u\_x^2\; +\; w\; (v\_t\; +\; u\; v\_x)$

which produces the Hunter–Saxton equation after elimination of "v" and "w" from the Euler–Lagrange equations for "u", "v", "w". Since there is also the more obvious Lagrangian

:$mathcal\{L\}\_1\; =\; u\_x\; u\_t\; +\; u\; u\_x^2,$

the Hunter–Saxton has two inequivalent variational structures. Hunter and Zheng also obtained a bihamiltonian formulation and a Lax pair from the corresponding structures for the Camassa–Holm equation in a similar way.

The fact that the Hunter–Saxton equation arises physically in two different ways (as shown above) was used by Alì and Hunter to explain why it has this bivariational (or bihamiltonian) structure.

Notes

References

* Citation
last = Alì
first = Giuseppe
author-link =
last2 = Hunter
first2 = John K.
year = 2006
title = Orientation waves in a director field with rotational inertia
periodical = Preprint arXiv:math/0609189v1 [math.AP]
volume =
issue =
pages =
url = http://arxiv.org/abs/math.AP/0609189
doi =

* Citation
last = Beals
first = Richard
author-link =
last2 = Sattinger
first2 = David H.
last3 = Szmigielski
first3 = Jacek
year = 2001
title = Inverse scattering solutions of the Hunter–Saxton equation
periodical = Applicable Analysis
volume = 78
issue = 3–4
pages = 255–269
url = http://www.math.usu.edu/~dhs/hunter209.ps
doi = 10.1080/00036810108840938

* Citation
last = Bressan
first = Alberto
author-link =
last2 = Constantin
first2 = Adrian
year = 2005
title = Global solutions of the Hunter–Saxton equation
periodical = SIAM J. Math. Anal.
volume = 37
issue = 3
pages = 996–1026
url =
doi = 10.1137/050623036

* Citation
last = de Gennes
first = Pierre-Gilles
author-link = Pierre-Gilles de Gennes
last2 = Prost
first2 = Jacques
author2-link =
year = 1994
title = The Physics of Liquid Crystals
edition = 2nd
volume =
series = International Series of Monographs on Physics
place =
publisher = Oxford University Press
isbn = 0-19-852024-7

* Citation
last = Holden
first = Helge
author-link =
last2 = Karlsen
first2 = Kenneth Hvistendahl
last3 = Risebro
first3 = Nils Henrik
year = 2007
title = Convergent difference schemes for the Hunter–Saxton equation
periodical = Math. Comp.
volume = 76
issue = 258
pages = 699–744
url =
doi = 10.1090/S0025-5718-07-01919-9

*Citation
last = Hunter
first = John K.
last2 = Saxton
first2 = Ralph
year = 1991
title = Dynamics of director fields
periodical = SIAM J. Appl. Math.
volume = 51
issue = 6
pages = 1498–1521
url =
doi = 10.1137/0151075

* Citation
last = Hunter
first = John K.
author-link =
last2 = Zheng
first2 = Yuxi
year = 1994
title = On a completely integrable nonlinear hyperbolic variational equation
periodical = Physica D
volume = 79
issue = 2–4
pages = 361–386
url =
doi = 10.1016/S0167-2789(05)80015-6

* Citation
last = Hunter
first = John K.
author-link =
last2 = Zheng
first2 = Yuxi
year = 1995
title = On a nonlinear hyperbolic variational equation. I. Global existence of weak solutions
periodical = Arch. Rational Mech. Anal.
volume = 129
issue = 4
pages = 305–353
url =
doi = 10.1007/BF00379259

* Citation
last = Hunter
first = John K.
author-link =
last2 = Zheng
first2 = Yuxi
year = 1995
title = On a nonlinear hyperbolic variational equation. II. The zero-viscosity and dispersion limits
periodical = Arch. Rational Mech. Anal.
volume = 129
issue = 4
pages = 355–383
url =
doi = 10.1007/BF00379260

* Citation
last = Lenells
first = Jonatan
author-link =
year = 2007
title = The Hunter–Saxton equation describes the geodesic flow on a sphere
periodical = J. Geom. Phys.
volume = 57
issue = 10
pages = 2049–2064
url =
doi = 10.1016/j.geomphys.2007.05.003