transitjon in subcritjcal shear fmows invariant solutjons
play

Transitjon in subcritjcal shear fmows Invariant solutjons and the - PowerPoint PPT Presentation

Transitjon in subcritjcal shear fmows Invariant solutjons and the edge of chaos. Ashley P. Willis 1 , Rich Kerswell 2 , Predrag Cvitanovi 3 , Yohann Duguet 4 1 School of Mathematjcs and Statjstjcs, University of Sheffjeld. 2 DAMTP,


  1. Transitjon in subcritjcal shear fmows – Invariant solutjons and the edge of chaos. Ashley P. Willis 1 , Rich Kerswell 2 , Predrag Cvitanović 3 , Yohann Duguet 4 1 School of Mathematjcs and Statjstjcs, University of Sheffjeld. 2 DAMTP, University of Cambridge. 3 School of Physics, Georgia Tech. (chaosbook.org) 4 LIMSI-CNRS, Orsay, France.

  2. Localised turbulence in a pipe Peixinho & Mullin, PRL Illuminated fmakes Simulatjon Axial vortjcity

  3. APPROACH: TURBULENCE AS A CHAOTIC DYNAMICAL SYTSTEM

  4. Trajectory in phase space, structured by stable/unstable manifolds of the equilibrium points.

  5. Knoll and Keyes (2004) Jacobian-free Newton-Krylov Viswanath (2007) → GMRES - Only involves evaluatjons of F ( x ). - No preconditjoner necessary!

  6. Knoll and Keyes (2004) Jacobian-free Newton-Krylov Viswanath (2007)

  7. Knoll and Keyes (2004) Jacobian-free Newton-Krylov - Hookstep Viswanath (2007) GMRES: Stabilise Newton: Don’t take too large step δx ...

  8. Code at openpipefmow.org (non-problem specifjc: black box) arxiv:1908.06730

  9. PROBLEMS 1. Where to get startjng guess x 0 !? 2. Look for recurrences: ‘small’ how small? what norm!?

  10. Osborne Reynolds’ Experiments, 1883 Observed importance of combination Re = LU / ν L, diameter U, mean axial flow ν, kinematic viscosity “The only idea I had formed before commencing the experiments, was that at some critical velocity the motion must become unstable, so that any disturbance from perfectly steady motion would result in eddies.” i.e. surprised to not find critical flow rate for linear instability “…the steady motion breaks down suddenly… for disturbances of the magnitude that cause it to break down… while it is stable for a smaller disturbance…” i.e. finite amplitude disturbance required to trigger turbulence

  11. Turbulent frictjon / Laminar frictjon (drag). Turbulent frictjon Laminar frictjon Turbulent 10 2 10 slope approx 0.75 Intermittent Laminar 1 10 3 10 4 10 5 10 6 Re ( fmow rate ) Adapted from Nikuradse (1950) / Blasius (1913)

  12. Subcritjcal instability, nonlinearity important || u’ || || u’ || ‘turbulence’ ? Saddle-node bifurcatjon Re Re laminar-turbulent linearly stable Disconnected from boundary laminar state. or ‘edge of chaos’ Q: how to fjnd states?

  13. Shear Flows U 2L U channel fmow Couetue fmow Reynolds number Re = LU /  Kinematjc viscosity  L 2U pipe fmow

  14. Stability of shear-fmows turbulence linear Re observed instability Pipe fmow 1720 inf.? 950 5772 Channel fmow ASBL 367 54370 312 inf. Couetue fmow

  15. Pipe fmow Re=81 lower bound, energy stability theory Re≈2000 T urbulence observed laminar Flow rate Re>770 (Faisst & Kerswell 2003) (Pringle & Kerswell 2007) lowest fjnite-amplitude solutjon lowest fjnite-amplitude solutjon m=2, m=3 m=1 Travelling Waves (TWs)

  16. Travelling Waves (TWs) / Vortex-Wave Interactjon (VWI) state / ‘Exact’ Coherent Structures (ECS) / Invariant Solutjons [Boundary-layer: ] Hall & Smith (1991), via asymptotjc theory. [Plane-couetue:] Walefge (1998), via contjnuatjon from Taylor-Couetue

  17. Travelling Wave solutjons (TWs) S2 S3 S4 S5 Streaks near Slower core walls Faisst & Eckhardt (2003) Wedin & Kerswell (2004) [Pipe:] via (painful) contjnuatjon from system with body force

  18. Discovery of TWs. (Self sustaining cycle completed ‘by hand’!) END: unforced 3D state, self-sustained Faisst & Eckhardt (2003) Wedin & Kerswell (2004) Reduce force 3D instability Force 2D streaks 0 Level of forcing - A lot can go wrong! START: - State really linked to dynamics? least-stable 2D eigen mode  Force

  19. Trajectory in phase space, structured by stable/unstable manifolds of the equilibrium points. Dimension of the space N → inf., In simulatjons, N = O(10 5 -10 6 ), dimension of unstable manifolds n = O(10)

  20. IF n=1 : Timestepping + Bisectjon between ICs . Turbulence Unstable TW becomes aturactor within ‘edge’ Laminar stable point Laminar turbulent boundary calculated by bisectjon : Skufca, Yorke & Eckardt (2006) for a reduced model of shear fmow Schneider, Eckhardt & Yorke (2007) for a short periodic pipe Itano & Toh (2000) for channel fmow.

  21. IF n>1 : Timestepping + Bisectjon between ICs . Turbulence Chaotjc aturactor within ‘edge’ Chaos within edge much Laminar milder than turbulence stable point → good candidates for Newton search Duguet, W. & Kerswell 2008,10 JFM long pipe, localised coherent structures within laminar-turbulent boundary

  22. Discovery of many (spatially periodic) TWs solutions for pipe flow

  23. Pufg-like invariant solutjons Avila, Mellibovsky, Rolland & Hof 2013 Exact localised periodic orbits found in m=2 + mirror space periodic solutjons Chantry, Willis & Kerswell 2014 Exact localised periodic orbits connected to periodic TWs via spatjal subharmonic L=2 π /α bifurcatjon. localised solutjon

  24. ‘Edge tracking’: Avila, Mellibovsky, Rolland & Hof 2013 Exact localised solution found in m=2 + mirror space 100s of simulations!

  25. TWs disconnected from laminar state. How to find them? A A Re Re Re = const. Re (t) = Re 0 + ĸ.(A 0 – A(t)) (Willis, Duguet, Omel’Chenko & Wolfrum, 2017, JFM)

  26. Pipe simulation: L = 2π /1.25 R, m = 2, no S&R etc. IC Re (t) = Re 0 + ĸ.(A 0 – A(t)) increasing ĸ

  27. Re (t) = Re 0 + ĸ.(A 0 – A(t)) increasing ĸ

  28. TW in ‘controlled’ and ‘uncontrolled’ system const. ĸ

  29. stabilized RPO unstable UPO

  30. ‘Method’: 1. Find a ‘suitable’ amplitude measure A 2. Link control parameter to A(t), e.g. Re (t) = Re 0 + к (A 0 – A(t)) 3. Increase slowly к = к(t) → reductjon in A(t) 4. Fix к if hit a stable point / orbit! (Willis, Duguet, Omel’Chenko & Wolfrum, 2017, JFM)

  31. Key points: • TWs (and POs) are weakly unstable solutjons of the N-S equatjons • Some are found in the laminar-turbulent boundary → transitjon • With hindsight, we could have found them yonks ago! (Willis, Duguet, Omel’Chenko & Wolfrum, 2017, JFM)

  32. Recurrent cycles (periodic orbits, POs) in turbulence

  33. Need to go into moving frame... yellow ,  2 = -0.3 blue , u z = -0.1

  34. Which phase speed !?

  35. ‘Slicing’

  36. ‘Slicing’

  37. Fourier ‘Slicing’ (Dynamics of TW just goes around circle)

  38. Slice vs Poincaré sectjon

  39. Recurrence plot (sliced dynamics)

  40. Sliced pipe

  41. Recurrence plot ‘Compensatory’ norm. ‘Crap’ norm ?

  42. Recurrence plot DMD/Koopman analysis ? Page & Kerswell (arXiv:1906.01310) Machine learning ? Page & Kerswell CSC methods ? Marensi & Willis

  43. Summary ● Jacobian-free Newton-Krylov (- Hookstep): workhorse for dynamical systems approach. arxiv:1908.06730 ● Relatjve equilibria (TWs) and relatjve periodic orbits (RPOs) embedded in laminar-turbulent boundary. ● Gettjng intjtal guesses for JFNK main issue. ● Bisectjon / ‘Surfjng’ edge. ● ‘Slicing’ (symmetry reductjon) ● RPOs embedded in turbulence → proxy for turbulence ● Norm problem.

Recommend


More recommend