Spatial symmetry invariance of solution of Kolmogorov flow

This paper presents a mathematical proof demonstrating that the spatial symmetry of two-dimensional Kolmogorov flow solutions is preserved for all time, thereby validating clean numerical simulation (CNS) results while exposing the rapid symmetry-breaking artifacts caused by numerical noise in direct numerical simulation (DNS).

The Gist

Here is an explanation of the paper, translated into simple language with creative analogies.

The Big Picture: A Mathematical "Mirror" Test

Imagine you have a perfectly symmetrical pattern drawn on a piece of paper. If you fold the paper in half, the left side matches the right side perfectly. Now, imagine you start a simulation (a computer model) of how that pattern moves and changes over time.

The Core Question: If the pattern starts perfectly symmetrical, does the true physics of the universe guarantee it stays symmetrical forever? Or will it eventually become messy and asymmetrical just because of the way it moves?

The Author's Answer: Shijun Liao, a mathematician, proves that if the starting pattern is symmetrical, the true mathematical solution must stay symmetrical forever. It is a law of the universe for this specific type of fluid flow.

However, there is a catch: Most computer simulations are lying to us. They break the symmetry quickly, not because the physics changed, but because the computer is making tiny mistakes.

---

The Characters in the Story

1. The Kolmogorov Flow: Think of this as a giant, swirling bathtub of water in a square box. It's being pushed by a specific force (like a hand pushing the water in a rhythmic pattern). It's a classic test case for turbulence (chaotic, swirling water).
2. The "Butterfly Effect": You know the idea that a butterfly flapping its wings can cause a tornado weeks later? In this paper, it means that even the tiniest, invisible error in a computer calculation can grow into a massive, visible mistake very quickly.
3. DNS (Direct Numerical Simulation): This is the "Standard Computer Simulation." It's like trying to draw a perfect circle using a pencil and a ruler, but your hand shakes just a tiny bit. Over time, the circle gets wobbly and ugly.
4. CNS (Clean Numerical Simulation): This is the "Super-Precision Simulation." It's like using a laser-guided robotic arm to draw that circle. It eliminates the shaking so effectively that the circle stays perfect for a very long time.

---

The Problem: The "Broken Mirror"

In the past, scientists used DNS (the standard method) to simulate this swirling water. They started with a perfectly symmetrical pattern.

• What happened? Within a very short time, the symmetry broke. The left side looked different from the right side.
• The Conclusion: Scientists thought, "Ah, the turbulence naturally destroys symmetry! The flow is chaotic and unpredictable."

Liao's Discovery:
Liao ran the same simulation using CNS (the super-precise method).

• What happened? The symmetry never broke. The left side stayed perfectly matched to the right side for a very long time.
• The Realization: The standard simulations (DNS) weren't showing us the truth of the fluid; they were showing us the noise of the computer. The "chaos" they saw was actually just computer errors growing so big they looked like real physics.

The Mathematical Proof: The "Infinite Chain"

How did Liao prove this? He didn't just run a simulation; he wrote a mathematical proof.

Imagine the movement of the water is a story being told one word at a time.

1. The Start: The story begins with a symmetrical sentence (the initial condition).
2. The Next Word: To know what happens next, you look at the current sentence. If the current sentence is symmetrical, the rules of physics (the Navier-Stokes equations) say the next word must also be symmetrical.
3. The Chain: Liao proved that if the first word is symmetrical, the second must be, the third must be, and so on, forever.

He showed that as long as the math works (which it does for this specific flow), the symmetry is locked in. It is impossible for the true solution to lose its symmetry.

The Verdict: Who is Right?

• The Standard Simulations (DNS): They are like a broken mirror. They show a distorted image because the "glass" (the computer's math) is dirty with tiny numerical errors. These errors act like a tiny push that grows into a giant shove, destroying the symmetry.
• The Clean Simulations (CNS): They are like a perfect mirror. They show that the symmetry is actually preserved.
• The Math Proof: This is the judge. It confirms that the CNS is right and the DNS is wrong (in this specific context).

Why Does This Matter?

This paper is a wake-up call for the world of fluid dynamics.

1. We might have been wrong about Chaos: We thought turbulence was so chaotic that it destroyed all patterns instantly. Liao suggests that some patterns (like symmetry) are actually very robust and should survive, but our computers are too "noisy" to see it.
2. The "Ultra-Chaos" Concept: The paper suggests that turbulence is "ultra-chaotic." This means that two starting points that are almost identical (differing by a tiny, invisible amount) will end up looking completely different.
   • Analogy: Imagine two identical twins starting a race. If one trips on a microscopic grain of sand (a tiny numerical error), they might end up running in opposite directions.
3. A New Tool: The author argues that we need to use "Clean Numerical Simulation" (CNS) to find the "mathematical truths" of turbulence, rather than relying on standard methods that are polluted by computer noise.

In a Nutshell

The paper proves that if you start with a symmetrical swirl of water, it should stay symmetrical forever.

If your computer simulation shows it becoming messy, it's not the water that's messy; it's your computer. The standard way of doing these simulations is so full of tiny errors that it creates fake chaos. By using a cleaner, more precise method, we can see the true, beautiful symmetry that the math promises.

Technical Summary

Here is a detailed technical summary of the paper "Mathematical proof about spatial symmetry of solutions of the two-dimensional Kolmogorov flow" by Shijun Liao.

1. Problem Statement

The paper addresses a fundamental discrepancy in the numerical simulation of turbulent flows governed by the Navier-Stokes (NS) equations. Specifically, it investigates the two-dimensional (2D) Kolmogorov flow under periodic boundary conditions.

• The Conflict: Traditional Direct Numerical Simulations (DNS) of chaotic NS turbulence often lose the spatial symmetry present in the initial conditions very quickly. In contrast, Clean Numerical Simulations (CNS) preserve this symmetry for much longer durations.
• The Core Question: Is the loss of symmetry in DNS a physical phenomenon of the turbulence itself, or is it an artifact of numerical errors (noise)?
• Hypothesis: The author posits that the NS equations possess an inherent mathematical property where spatial symmetry (rotation and/or translation) is preserved for all time $t > 0$ if the initial condition is smooth and symmetric. Therefore, any simulation that loses this symmetry is mathematically incorrect due to numerical pollution.

2. Methodology

The paper employs a rigorous mathematical proof based on the properties of the Navier-Stokes equations and Taylor series expansions, rather than relying solely on numerical experimentation.

• Governing Equations: The study focuses on the dimensionless NS equation for 2D Kolmogorov flow expressed via a stream-function $\psi$:
  $$\frac{\partial}{\partial t}(\nabla^2\psi) + \psi_x\nabla^2\psi_y - \psi_y\nabla^2\psi_x - \frac{1}{Re}\nabla^4\psi + n_K \cos(n_K y) = 0$$
  where $n_K$ is an even integer, and periodic boundary conditions apply.
• Symmetry Definitions: The initial condition $\psi(x, y, 0)$ is assumed to possess:
  • Rotational Symmetry: $\psi(x, y, 0) = \psi(2\pi - x, 2\pi - y, 0)$
  • Translational Symmetry: $\psi(x, y, 0) = \psi(x + \pi, y + \pi, 0)$
• Analytical Approach:
  1. Taylor Series Expansion: The solution $\psi(x, y, t)$ is expanded as a Taylor series around an arbitrary time $t_0$.
  2. Inductive Proof: The author proves that if the solution and its time derivatives up to order $m$ satisfy the symmetry at $t_0$, then the $(m+1)$-th derivative also satisfies the symmetry.
  3. Coordinate Transformation: The proof utilizes coordinate transformations ($x' = 2\pi - x$, etc.) to demonstrate that the differential operators ($\nabla^2, \nabla^4$) and the forcing term ($\cos(n_K y)$) remain invariant under the symmetry operations, provided $n_K$ is even.
  4. Extension to Global Time: By chaining local time intervals where the Taylor series converges (assuming a non-zero radius of convergence for all $t \geq 0$), the symmetry is extended to all $t > 0$.

3. Key Contributions

The paper makes three primary theoretical contributions:

1. Theorem 1 (Inductive Symmetry): Proves that if a solution to the 2D Kolmogorov flow possesses spatial symmetry at a specific time $t_0$, all higher-order time derivatives of the solution at that time also possess the same symmetry.
2. Theorem 2 (Local Preservation): Establishes that if the solution is symmetric at $t_0$ and the Taylor series has a non-zero radius of convergence, the symmetry is preserved throughout the interval $[t_0, t_0 + \rho)$.
3. Theorem 3 (Global Preservation): The central result. It states that for smooth initial conditions with spatial symmetry, the solution to the 2D Kolmogorov flow retains this symmetry for all $t > 0$, provided the Taylor series convergence radius remains non-zero.

Comparison of Simulation Methods:

• CNS (Clean Numerical Simulation): Uses high-order methods and arbitrary-precision arithmetic to reduce truncation and round-off errors to negligible levels. CNS results agree with Theorem 3, preserving symmetry.
• DNS (Direct Numerical Simulation): Uses standard double-precision arithmetic. DNS results violate Theorem 3 by losing symmetry quickly, confirming that DNS results are "polluted" by numerical noise that acts as an artificial stochastic disturbance.

4. Results

• Validation of CNS: The mathematical proof rigorously confirms previous CNS findings (from references [1, 2]) that 2D Kolmogorov flow maintains spatial symmetry over long time intervals.
• Invalidation of Standard DNS: The proof demonstrates that standard DNS results, which lose symmetry, are mathematically inconsistent with the exact NS equations for this specific flow. The loss of symmetry is attributed to the exponential growth of numerical noise (the "butterfly effect" applied to numerical errors).
• Ultra-Chaos and Non-Uniqueness: The paper discusses the concept of "ultra-chaos," where tiny disturbances (like $10^{-20}$) can lead to macroscopic differences in both trajectory and statistical properties.
  • It supports Conjecture A: The NS equations may admit distinct smooth global solutions from almost identical initial conditions.
  • It supports Conjecture B: Specifically for 2D turbulent Kolmogorov flow, distinct global solutions exist for initial conditions differing by an arbitrarily small $\delta$.

5. Significance

• Rigorous Benchmarking: The paper provides a mathematical "ground truth" (Theorem 3) to evaluate the reliability of numerical simulations. It argues that spatial symmetry is a critical metric for validating turbulence simulations.
• Critique of DNS: It challenges the reliability of traditional DNS for chaotic systems, suggesting that standard numerical noise is not just a minor error but a fundamental contaminant that alters the physical nature of the solution (turning a symmetric flow into a non-symmetric one).
• Role of CNS: It elevates CNS from a numerical technique to a tool for discovering mathematical truths about the NS equations, such as the potential non-uniqueness of smooth solutions and the nature of turbulence as "ultra-chaos."
• Implications for Turbulence Theory: The findings suggest that turbulence is not just chaotic in trajectory but in statistical properties as well, and that the "randomness" observed in turbulence may be partially an artifact of numerical noise in traditional simulations.

In conclusion, Liao proves that the 2D Kolmogorov flow is mathematically bound to preserve its initial spatial symmetry. Consequently, any simulation that fails to do so is numerically flawed, providing strong theoretical backing for the use of Clean Numerical Simulations in studying turbulence.