Spatio-Temporal Chaos in Pattern-Forming Systems: Defects and Bursts

Transcription

1 Spatio-Temporal Chaos in Pattern-Forming Systems: Defects and Bursts with Santiago Madruga, MPIPKS Dresden Werner Pesch, U. Bayreuth Yuan-Nan Young, New Jersey Inst. Techn. DPG Frühjahrstagung supported by NASA and DOE

2 Pattern Formation Symmetry-breaking instabilities patterns Rayleigh-Bénard convection of a fluid layer heated from below T T+dT (Plapp & Bodenschatz, 1996) (Bodenschatz, de Bruyn, Ahlers, Cannell, 1991)

3 Spatio-Temporal Chaos Spiral-Defect Chaos Küppers-Lortz Chaos (Morris, Bodenschatz, Cannell, Ahlers, 1993) (Hu, Ecke, Ahlers, 1997) Small Prandtl number: large-scale flows Rotation: Küppers-Lortz Instability

4 Goals: Quantitative characterization of spatio-temporally chaotic states and of transitions between them (HR & Madruga, 2006) Origin of temporal and spatial chaos and the mechanisms maintaining it Microscopic equations Reduction Order parameters Symmetries Macroscopic equations

5 Non-Boussinesq Convection T T+dT (Bodenschatz, de Bruyn, Ahlers, Cannell, 1991)

6 Hexagons in Rotating Systems Weakly nonlinear description of hexagons without rotation with rotation γ=0 Rolls γ>0 Rolls Steady Hexagons Steady Hexagons Whirling Hexagons µ µ (J.W. Swift, 1984; Soward, 1985) Hopf bifurcation to whirling hexagons weak non-boussinesq effects: coupled Ginzburg-Landau equations CGL defect chaos (Echebarria & HR, 2000)

7 Rotating Non-Boussinesq Convection - Navier-Stokes equation for momentum V i ( ( )) 1 Vi t V i + V j j Pr ρ = i p + δ i3 ( 1 + γ 1 ( 2z + Θ R ) ) Θ + + j [νρ - Heat equation and continuity equation - Weakly temperature-dependent fluid parameters ( i ( V j ρ ) + j( V )] i ρ ) + 2Ωǫ ij3 V j T T 0 ρ(t) = 1 γ 0 R (1 + γ T T 0 1 R ) ν(t) = 1 + γ T T 0 2 R...

8 Rotating Non-Boussinesq Convection - Navier-Stokes equation for momentum V i ( ( )) 1 Vi t V i + V j j Pr ρ = i p + δ i3 ( 1 + γ 1 ( 2z + Θ R ) ) Θ + + j [νρ - Heat equation and continuity equation - Weakly temperature-dependent fluid parameters ( i ( V j ρ ) + j( V )] i ρ ) + 2Ωǫ ij3 V j T T 0 ρ(t) = 1 γ 0 R (1 + γ T T 0 1 R ) ν(t) = 1 + γ T T 0 2 R... Fully nonlinear hexagon solution linear stability analysis Direct numerical simulation of chaotic states Interpretation: Reduction to complex Ginzburg-Landau equation

9 Reentrant Hexagons in Non-Rotating Convection Stability of hexagons with respect to amplitude perturbations Red. Rayleigh Number ε Stable Hexagons Stable Rolls h=1.8 mm Unstable Rolls Reentrant Hexagons Unstable Hexagons Mean Temperature T 0 Non-Boussinesq effects increase with decreasing mean temperature T 0

10 Reentrant Hexagons in Non-Rotating Convection Stability of hexagons with respect to amplitude perturbations Red. Rayleigh Number ε Stable Hexagons Stable Rolls h=1.8 mm Unstable Rolls Reentrant Hexagons Unstable Hexagons Mean Temperature T 0 Strong non-boussinesq effects (low mean temperature): - no instability of hexagons to rolls - coexistence of stable hexagons and rolls

11 Reentrant Hexagons in Non-Rotating Convection Stability of hexagons with respect to amplitude perturbations Red. Rayleigh Number ε Stable Hexagons Stable Rolls h=1.8 mm Unstable Rolls Reentrant Hexagons Unstable Hexagons Mean Temperature T 0 Weak non-boussinesq effects

12 Reentrant Hexagons in Non-Rotating Convection Stability of hexagons with respect to amplitude perturbations Red. Rayleigh Number ε Stable Hexagons Stable Rolls h=1.8 mm Unstable Rolls Reentrant Hexagons Unstable Hexagons Mean Temperature T 0 Weak non-boussinesq effects Intermediate non-boussinesq effects

13 Reentrant Hexagons in Non-Rotating Convection Stability of hexagons with respect to amplitude perturbations Red. Rayleigh Number ε Stable Hexagons Stable Rolls h=1.8 mm Unstable Rolls Reentrant Hexagons Unstable Hexagons Mean Temperature T 0 Weak non-boussinesq effects Intermediate non-boussinesq effects Strong non-boussinesq effects

14 Weak Non-Boussinesq Effects water: thickness h = 4.92mm mean temperature T 0 = 12 o C critical temperature difference T c = 6.4 o C rotation rate Ω = 65 ( 1 Hz) Hopf bifurcation at ǫ = 0.07 Amplitude (a.u.) steady oscillating Oscillation Amplitude Red. Rayleigh Number ε ǫ = 0.2

15 Description within CGL Framework Extract oscillation amplitude H(X, T) wavevectors near q n frequencies near ω H v x (x, t, z = 0) = 3 n=1 ( R + [ ]) e 2πni/3 H(X, T)e iωht + c.c. t H = µh + (1 + ib 1 ) 2 H (b 3 i) H H 2 ~ q 2 ~ q 3 ~ q 1 exp (i q n x) +... Complex Ginzburg-Landau equation: the universal description of weakly nonlinear oscillations

16 Description within CGL Framework Extract oscillation amplitude H(X, T) wavevectors near q n frequencies near ω H v x (x, t, z = 0) = 3 n=1 ( R + [ ]) e 2πni/3 H(X, T)e iωht + c.c. t H = µh + (1 + ib 1 ) 2 H (b 3 i) H H 2 ~ q 2 ~ q 3 ~ q 1 exp (i q n x) +... CGL Defect Chaos (H r ) b S 2 Defect Chaos Frozen Vortices Phase Chaos b 3 T L BF Stable Plane Waves (Chaté & Manneville, 1996)

17 Description within CGL Framework Extract oscillation amplitude H(X, T) wavevectors near q n frequencies near ω H v x (x, t, z = 0) = 3 n=1 ( R + [ ]) e 2πni/3 H(X, T)e iωht + c.c. t H = µh + (1 + ib 1 ) 2 H (b 3 i) H H 2 ~ q 2 ~ q 3 ~ q 1 exp (i q n x) +... Extract coefficients b 1 and b 3 from direct Navier-Stokes simulations Bistability: Plane waves defect chaos b L Defect Chaos Stable Plane Waves BF b 3 T (Chaté & Manneville, 1996)

18 Description within CGL Framework Extract oscillation amplitude H(X, T) wavevectors near q n frequencies near ω H v x (x, t, z = 0) = 3 n=1 ( R + [ ]) e 2πni/3 H(X, T)e iωht + c.c. t H = µh + (1 + ib 1 ) 2 H (b 3 i) H H 2 ~ q 2 ~ q 3 ~ q 1 exp (i q n x) +... Navier-Stokes Simulation ( H ) 1 L Defect Chaos BF T b Stable Plane Waves b 3 (Chaté & Manneville, 1996)

19 Intermediate Non-Boussinesq Effects h = 4.92mm Ω = 65 Pr = 8.7 T 0 = 14 o C T c = 8.3 o C Ω = 65 Hopf bifurcation backward Hysteresis and bistability ǫ = 0.5 of steady and oscillating hexagons Amplitude (a.u.) steady hexagons oscillating hex ε H Red. Rayleigh ε restabilization of steady hexagons at larger ǫ fluctuating localized domains of whirling hexagons

20 Quintic Complex Ginzburg-Landau Equations Quintic Ginzburg-Landau equation t H = µh + (d r + id i ) 2 H (c r + ic i )H H 2 (g r + ig i )H H 4 Extract coefficients d = i, c = i, g = i Demodulation Navier-Stokes quintic Ginzburg-Landau

21 Localization Mechanism CGL coupled to phase modes t H = (µ Qδ ϕ)h + d 2 H ch H 2 gh H 4 Two contributions: 1. Wavenumber selection by front H = Re iψ k = ψ large c i, g i : gradients in the oscillation magnitude R induce wavevector k diffusion d r : oscillation amplitude damped H k (Bretherton and Spiegel, 1983;...,Coullet and Kramer, 2004)

22 Localization Mechanism CGL coupled to phase modes t H = (µ Qδ ϕ)h + d 2 H ch H 2 gh H 4 Two contributions: 1. Wavenumber selection by front 2. Compression ϕ of underlying hexagon pattern Pattern Compression ϕ Growthrate σ A ε=0.2 B ε= Wavenumber q

23 Strong Non-Boussinesq Effects h = 4.6mm Ω = 65 Pr = 8.7 T 0 = 12 C T c = 10 C no Hopf bifurcation to whirling hexagons only side-band instabilities Reduced Rayleigh (R-R c )/R c steady hexagons linearly stable Wavenumber q ǫ = 1.0 whirling destroys order of hexagonal lattice

24 Strong Non-Boussinesq Effects h = 4.6mm Ω = 65 Pr = 8.7 T 0 = 12 C T c = 10 C no Hopf bifurcation to whirling hexagons only side-band instabilities Reduced Rayleigh (R-R c )/R c steady hexagons linearly stable Wavenumber q ǫ = 0.87 whirling destroys order of hexagonal lattice

25 Defect Statistics ǫ = 1 Triangulation Heptagons Pentagons Defect Statistics Broader than squared Poisson: Correlations Relative Frequency N

26 Conclusions Defects and bursts in rotating non-boussinesq convection Whirling hexagons - reduction to CGL - 2d CGL defect chaos Localized bursts of oscillations quintic CGL: retracting fronts, collapse compression of lattice Whirling chaos whirling penta-hepta defects Phys. Rev. Lett. 96 (2006) , J. Fluid Mech. 548 (2006) 341, New J. Phys. 5 (2003)