Long-wave Instability in Anisotropic Double-Diffusion Jean-Luc Thiffeault Institute for Fusion Studies and Department of Physics University of Texas at Austin and Neil J. Balmforth Department of Theoretical Mechanics University of Nottingham November 1997 APS DFD-97
Overview Want to capture asymptotic dynamics near Takens Bogdanov bifurcation. Problem: typical scaling leads to Hamiltonian (and thus conservative) equation, which obviously does not capture a lot of the dynamics. Try using different scaling, but then get unremovable resonant terms. Solution: extend parameter space to allow removal of resonant terms. Raises codimension, but asymptotic. APS DFD-97 1
Takens Bogdanov TB bifurcation occurs when two modes become unstable at the same parameter values. Equations for the reduced dynamics near this bifurcation point capture more of the diverse behaviour of the system than simple steady or Hopf bifurcation. For double-diffusive convection in long-wave theory such a bifurcation is present. Problem: the reduced equations contain terms of differing order in the standard asymptotic expansion parameter. The asymptotic theory fails to collect a dissipative nonlinear term; the amplitude equations are Hamiltonian to leading order (Childress and Spiegel 1981). APS DFD-97 2
Possible solutions Normal form theory: not available for extended (continuum of excited modes) systems. Reconstitution: Not asymptotic, so hard to judge validity. May be flawed in some cases (Clune, Depassier, and Knobloch, 1994). Nonlocal averaging: 1988). Difficult to solve (Pismen, Alternative route: if more parameters were available, could remove resonant terms at the cost of augmenting the codimension of the bifurcation. To introduce needed extra parameters, we choose anisotropic double-diffusion as our system. (possible transport model for ocean, astrophysics, tokamak plasmas) APS DFD-97 3
Illustration of Procedure Normal form for three real marginal modes: F = G Ġ = H Ḣ = η H ν G λ F + a F 2 + b G 2 + c F G + d F H Assuming strongly damped mode ( η >> ν, λ ) we should recover the two-mode normal form. One way to do this (Spiegel et al) is to use the scaling t = t/δ, λ = δ 2 λ, ν = δ2 ν, F = δ 2 F, G = δ 3 Ḡ. This leads to a Hamiltonian equation, not two-mode normal form as one would expect. If instead of rescaling the amplitudes one rescales the nonlinear terms t = t/δ, λ = δ 2 λ, ν = δ ν, a = δ2 ā, c = δ c, we recover two-mode normal form, at the cost of raising the codimension. APS DFD-97 4
Model Equations The equations for anisotropic double-diffusion are σ 1 d dt 2 ψ = R x Θ S x Σ + (D 2 + 2 x) 2 ψ, d dt Θ = xψ + (D 2 + Λ 2 x)θ, Le d dt Σ = Le xψ + (D 2 + Ξ 2 x)σ; with no-slip, fixed-flux boundary conditions ψ = Dψ = 0, D Θ = D Σ = 0, z = 0 and 1 Fixed flux favors convection cells that are as large as the system will permit. Use this to define small parameter ɛ. Scaling: x = ɛ X, t = ɛ 4 T, ψ = ɛ φ X APS DFD-97 5
Order ɛ 0 and ɛ 2 The fixed flux boundary conditions give at order ɛ 0. Θ 0 = Θ 0 (X, T ), Σ 0 = Σ 0 (X, T ) At order ɛ 2, we get the solvability condition (linear at this order): ( 1 720 R 0 Λ 0 1 720 S 0 1 720 Le 0 R 0 1 720 Le 0 S 0 Ξ 0 ) ( Θ0XX Σ 0XX ) = 0. The requirement that the matrix have zero eigenvalues means that its trace and determinant must vanish. This is obtained by letting R 0 = 720 Λ 2 0, S 0 = 720 Λ 0 Ξ 0 /Le 0 Le 2 0 Ξ 2 0 Λ 0 Ξ 0 /Le 0, The eigenvector for the matrix is parametrized by Σ 0 = (Le 0 Λ 0 /Ξ 0 )Θ 0 (it only has one). APS DFD-97 6
Order ɛ 4 Get two solvability conditions again, this time involving T : Θ 0T = Σ 0T = (Le 0 Λ 0 /Ξ 0 )Θ 0T = Must be compatible since Θ 0T and Σ 0T are related. This is not satisfied automatically; this is why we now make use of the extra parameters. By letting Le 0 = 1 5(Λ 0 + Ξ 0 ) = 11(1 + 0 ) R 2 Λ 0 Ξ 0 S 2 = 720Λ 0(Λ 2 Ξ 2 + Le 2 Ξ 0 ) Λ 0 Ξ 0 the two become compatible. codimension by three. This increases the APS DFD-97 7
Marginal Stability Curves 7400 7200 7000 R 6800 6600 6400 6200 0 1 2 3 4 5 k 0 = 0.1 (long-dashed), 1.2727 (solid), 5 (dashed) APS DFD-97 8
Order ɛ 6 We get a solvability condition involving only the ɛ 2 integration constants g(x, T ) := Σ 2,0 (X, T ) Λ 0 Ξ 0 Θ 2,0 (X, T ) at this order. After rescaling to eliminate some parameters we have the coupled system f T = g XX + α f XX + f XXXX + ( fx 3 ) X g T = λ f XX + κ f XXXX γ f XXXXXX + β g XX ρ g XXXX + ξ ( fx 3 )X + ( fx 2 ) g X X + η ( f X fxx 2 ) ζ ( f 3 ) X X XXX We fixed Le 0, 0, and Λ 2. However, we are left with enough parameters to vary independently all the coefficients except η and ζ. APS DFD-97 9
Captures Turnaround 2 1.75 1.5 1.25 1 0.75 0.5 0.25 Supercritical 2 4 6 8 10 Steady state Blue A f 3 Red A f 3 + B f 5 Critical 2 1.75 1.5 1.25 1 0.75 0.5 0.25-10 -5 5 10 Subcritical 3 2.5 2 1.5 1 0.5-20 -10 10 20 APS DFD-97 10
Numerical Solution 0.3 0.2 f 0.1 phase 0.0 0 100 200 300 t 4 2 0-2 -4 0 100 200 300 t Traveling waves stable. APS DFD-97 11
Conclusions For anisotropic double-diffusion in long-wave theory, we have shown that an extended system equation can be asymptotically derived. The equation contains several known equations as limits (Chapman&Proctor 1980, Childress&Spiegel 1981, Knobloch 1989). Compare reconstituted result. Explore numerical solutions. Make connection with physics. APS DFD-97 12