Conservative and skew-symmetric forms of the incompressible Navier-Stokes equations in sigma-coordinates

This paper derives new conservative and skew-symmetric formulations of the incompressible Navier-Stokes equations in terrain-following sigma-coordinates that preserve the structural and energy-conserving properties typically found in Cartesian systems.

The Gist

Imagine you are trying to map out the movement of water in a complex, rocky riverbed or the way wind swirls over a jagged mountain range.

To do this, scientists use math equations called the Navier-Stokes equations. These equations are the "rules of the road" for how fluids move. However, there is a massive problem: most of these math rules were written for a flat, perfect world (like a swimming pool). When you try to apply them to a bumpy, irregular world (like a mountain range), the math "breaks." It gets messy, loses its balance, and can even cause computer simulations to "explode" with errors.

This paper is about a new way to rewrite those rules so they work perfectly on "bumpy" terrain without losing their fundamental logic.

Here is the breakdown of how they did it, using three simple analogies.

---

1. The "Sigma-Coordinate" Problem: The Stretching Map

Imagine you have a standard rectangular grid (like graph paper) laid over a mountain. The mountains don't fit into the squares; they poke through them. To fix this, scientists use something called $\sigma$-coordinates (Sigma-coordinates).

Think of this like taking a piece of stretchy spandex and pulling it over the mountains so that the "grid lines" follow the shape of the slopes. Now, the mountains are "inside" the grid.

The Catch: When you stretch the spandex, the squares become tilted diamonds. In math, this "stretching" creates "ghost forces" (metric-induced terms) that weren't there before. These ghost forces act like fake wind or fake gravity, confusing the computer and making the simulation unstable.

2. The "Conservative Form": The Strict Accountant

The researchers first created a "Conservative Form."

Think of a Strict Accountant. If you are tracking money in a business, a conservative accountant ensures that every single penny is accounted for. If \100 leaves one pocket, exactly \100 must appear in another. You can't just "lose" money to the air.

In fluid math, "money" is things like mass and momentum. The researchers found a way to rewrite the equations so that even though the "spandex" is stretched over mountains, the "Accountant" still tracks every drop of water perfectly. This is great for "shock-capturing"—like simulating a sudden wave hitting a rock.

3. The "Skew-Symmetric Form": The Perfect Athlete

Next, they created a "Skew-Symmetric Form." This is a much more sophisticated tool.

Think of a Professional Athlete on a trampoline. An athlete wants to maintain "energy stability." They want to jump, use energy, and land, but they don't want the trampoline to suddenly gain infinite energy and launch them into space (which is what happens when a computer simulation fails).

The "Skew-Symmetric" version is like a mathematical "energy stabilizer." It ensures that the energy in the system stays balanced. It doesn't create "fake energy" out of nowhere. This is crucial for studying turbulence—the chaotic, swirling eddies you see in a rushing stream. Without this stability, the computer would see a tiny swirl and accidentally turn it into a massive, fake hurricane.

---

Summary: Which tool should you use?

The paper concludes that there isn't just one "right" way; it depends on what you are trying to do:

• If you are simulating a sudden, violent event (like a flash flood hitting a canyon), use the Conservative Form (The Strict Accountant). It’s best at handling sudden changes and "shocks."
• If you are simulating complex, swirling motion (like the tiny eddies in a turbulent river), use the Skew-Symmetric Form (The Perfect Athlete). It’s best at keeping the energy stable so the simulation doesn't crash.

In short: They have given engineers a better "mathematical toolkit" to simulate the messy, bumpy real world with much higher accuracy and reliability.

Technical Summary

Technical Summary: Conservative and Skew-Symmetric Forms of the Incompressible Navier-Stokes Equations in $\sigma$-Coordinates

1. Problem Statement

Simulating incompressible fluid flows (such as atmospheric boundary layers or oceanic currents) over complex, irregular, or moving terrain presents a significant numerical challenge. Traditional Cartesian grids struggle with boundary fidelity, while body-conforming grids are computationally expensive and difficult to implement for evolving geometries.

The widely used $\sigma$-coordinate system (a terrain-following coordinate) rescales the vertical dimension to match topography. However, this transformation breaks the local orthogonality of the grid, introducing "metric-induced terms." In conventional numerical schemes, simply transforming Cartesian equations into $\sigma$-coordinates disrupts the fundamental mathematical structures of the governing equations—specifically, the conservative structure (essential for shock-capturing and mass/momentum conservation) and the skew-symmetric structure (essential for energy stability and boundedness).

2. Methodology

The authors derive new mathematical formulations for the incompressible Navier-Stokes equations within the $\sigma$-coordinate framework that restore these lost structural properties.

• Conservative Formulation:
  The authors start with the general conservation laws in Cartesian coordinates and apply the $\sigma$-transformation. By carefully managing the Jacobian ($J$) of the transformation and utilizing the properties of mixed partial derivatives, they derive a system where the mass and momentum equations are expressed in a strictly conservative form. They specifically apply this to the artificial compressibility method (a pseudo-compressibility approach) to analyze its eigenstructure.
• Skew-Symmetric Formulation:
  To ensure energy stability, the authors introduce a novel set of transformed variables:
  $$P \equiv \sqrt{Jp/\rho_o}, \quad U \equiv \sqrt{Ju}, \quad W \equiv \sqrt{Jw}$$
  By substituting these into the momentum equations, they construct a formulation where the convective terms are written in a skew-symmetric manner. This allows the use of the energy method (via summation-by-parts operators) to analyze the system's stability.
• Boundary Condition Analysis:
  The study employs characteristic-based boundary conditions. By diagonalizing the flux Jacobians, the authors determine exactly how many boundary conditions are required at the lateral, top, and bottom boundaries to ensure the system remains well-posed and energy-bounded.

3. Key Contributions

• Restoration of Mathematical Structure: The paper provides a rigorous derivation of both conservative and skew-symmetric forms that are specifically tailored for the $\sigma$-coordinate system, preventing the "structural inconsistencies" found in standard transformations.
• Eigenstructure Analysis: For the conservative form, the authors provide a complete analysis of the eigenvalues and eigenvectors, proving that the $\sigma$-transformation modifies the eigenstructure but maintains the necessary hyperbolic properties for pseudo-compressibility methods.
• Energy Stability Proof: The authors mathematically demonstrate that the new skew-symmetric formulation is:
  • Energy-conserving for the inviscid (Euler) equations.
  • Energy-bounded for the viscous (Navier-Stokes) equations.
• Physical Interpretation of Boundary Work: The study provides a detailed breakdown of how moving terrain ($z_t \neq 0$) performs work on the fluid through pressure, and how viscous stresses contribute to energy flux at the boundaries.

4. Results and Findings

• Conservative Form Utility: The derived conservative form is shown to be ideal for weakly compressible or pseudo-compressible modeling, as it preserves the ability to use Riemann solvers and capture sharp gradients.
• Skew-Symmetric Form Utility: The skew-symmetric form successfully handles the non-orthogonality of the $\sigma$-grid while ensuring that the kinetic energy evolution is governed solely by boundary terms and internal viscous dissipation.
• Boundary Requirements: The analysis identifies that for lateral boundaries, the number of required boundary conditions depends on the direction of the flow (inflow vs. outflow), and for the $\sigma$-boundaries, the number of conditions is dictated by the sign of the eigenvalues of the normal flux Jacobian.

5. Significance

This work is highly significant for the fields of computational fluid dynamics (CFD), meteorology, and oceanography.

By providing a mathematically sound way to handle terrain-following coordinates, the authors offer a roadmap for developing more stable and physically accurate numerical models. Researchers can now choose between a conservative formulation (for robust shock-capturing and high-speed/pseudo-compressible flows) and a skew-symmetric formulation (for high-fidelity, energy-stable simulations of turbulence), both of which are guaranteed to respect the underlying physics of the incompressible Navier-Stokes equations even over highly complex topography.