Back to Experiments

Shell-model experiment: Feedback pattern across parameters

Question

The system studied here is a cubic shell ODE with independent velocity-shell variables \(u_k\) and vorticity-shell variables \(\omega_k\), with stretching couplings \(\lambda\), \(\lambda_2\) and a feedback coupling \(\sigma\). It is not a Galerkin truncation of Navier–Stokes. The threshold \(A^*\) reported below is a finite-window excursion threshold: the initial amplitude at which a trajectory first crosses the excursion threshold within the integration window. Its value, the times to threshold crossing and the classification scores depend on the integration window, the excursion threshold and the solver.

Is the \(L_1 \to L_4 \to L_2 \to L_1\) pattern (shell enstrophy rising while the S-score falls) specific to one parameter choice, or does it appear at every tested combination of viscosity \(\nu\) and coupling \(\lambda_2\) in this ODE?

Method

  1. Select a 3×3 grid of parameter combinations: \(\nu \in \{0.001, 0.005, 0.02\}\), \(\lambda_2 \in \{1.0, 5.0, 20.0\}\).
  2. For each combination, find \(A^*\) via bisection.
  3. Run the ODE at 80% of \(A^*\) — below threshold within the window but close enough to it to be informative.
  4. Check whether \(L_1\) increases AND \(L_2\) decreases at every checkpoint (the signature used for an active pattern).

Results

Parameter Sweep

\(\nu\)\(\lambda_2\)\(A^*\)Pattern at 80% \(A^*\)\(L_1\) trend\(L_2\) trend
0.0011.00.312ActiveIncreasingDecreasing
0.0015.00.289ActiveIncreasingDecreasing
0.00120.00.261ActiveIncreasingDecreasing
0.0051.00.421ActiveIncreasingDecreasing
0.0055.00.378ActiveIncreasingDecreasing
0.00520.00.334ActiveIncreasingDecreasing
0.021.00.606ActiveIncreasingDecreasing
0.025.00.523ActiveIncreasingDecreasing
0.0220.00.447ActiveIncreasingDecreasing

Summary

MetricResult
Combinations tested9
Pattern active (\(L_1 \uparrow\) AND \(L_2 \downarrow\))9/9
Depends on the parameter choice?Not within this grid
Observations. The pattern is active at all nine combinations of \(\nu\) and \(\lambda_2\) tested. \(A^*\) increases with \(\nu\) and decreases with \(\lambda_2\), ranging from 0.261 to 0.606 across the grid. The pattern is a feature of this ODE's cubic coupling and viscous term throughout the grid.

Summary

9/9. The \(L_1 \to L_4 \to L_2 \to L_1\) pattern is active at all nine parameter combinations in this shell ODE.

Related