**1.0 The Periodic Navier-Stokes Matrix Differential Equation**

Begin with the Navier-Stokes equation [1] describing the time evolution of an incompressible flow’s Fourier modes

(1)

where and are discrete wavenumbers with for some even integer and where . Included in (1) are , a known viscosity , the incompressibility condition , the projection tensor

, (2)

and

. (3)

Simultaneously consider all wavenumbers by writing (1) in consecutive trios of rows for . This forms the periodic Navier-Stokes matrix differential equation

(4)

where

(5)

and

, (6)

with

, (7)

, (8)

(9)

and

. (10)

**2.0 Power ****Series Solution**

Expand as a time series

(11)

with unknown coefficients and more simply . Substituting (11) into (1) and performing the time differentiation gives

(12)

with

. (13)

Matching coefficients in (12) when gives a recursion formula for the unknown flow coefficients

. (14)

It begins with a known initial flow with an assumed magnitude of unity. As equations (11) and (14) describe the continuous time series solution to the periodic Navier-Stokes equation.

**3.0 The Recursion Formula**

Reducing the dimensionality of the quantities in (14) to a scalar, beginning with and ignoring the wavenumber factors quickly reveals that this recursion formula creates coefficients geometric in nature , and therefore time series coefficients admitting only a finite time interval of convergence. Extensive numerical simulations evaluating for increasing values of , using (14) and mostly white-noise water flows with a paltry and without the no-slip boundaries, corroborate this universally-predicted, power-law blow-up behavior. This is far from a proof though.

**4.0 The Immortal Matrix Exponential Solution**

There is a way out of this apparent hopelessness. It involves restricting wavenumber portions of the solution to self-similar vector coefficients over all time orders and parallel to the initial flow. Numerical simulations show that solution coefficient vectors for successive time orders are highly correlated.

As a means of its mathematical introduction, notice that back substituting in (14) for and keeping just the first term involving etc., yields

, (15)

which when the time variable is added, represents the series expansion of . This is stable for non-positive real eigenvalues which is certainly the case for viscous and turns out to be the case for all convective . The convective eigenvalues are purely imaginary (or zero) and come in conjugate pairs. Its eigenvectors form a unitary basis. Neither are proven here.

The viscous part of (15) is well understood, but the convective part is not, and, along with (4) suggests the universally convergent matrix exponential solution [2]

. (16)

It arises when commute. If

(17)

is the eigendecomposition of , commutativity requires that eigenvectors are time-order independent

. (18)

In view of (3) and (10) one way of achieving (18) is to require self-similar

, (19)

for scalar .

Constraint (19) requires that the projection matrix project flow vectors parallel to initial flow rather than just onto the plane perpendicular to wavenumber. This can be achieved by

(20)

rather than (14). Most importantly, (20) maintains incompressibility because the initial flow by assumption is incompressible.

**Remarks**

The above sections describe a purely mathematical invention and construction of a regular solution to the periodic Navier-Stokes equation (1), and, as such, should stand on its own merit. Is it universal? No, because it assumes self-similar wavenumber solution coefficient vectors that are parallel to the initial flow. It simply explains the required assumption for the matrix exponential solution. It is motivated by the discovery through numerical simulations that solution coefficient vectors are highly correlated.

**References**

[1] Pope SB 2003 *Turbulent Flows* eq 6.146 Cambridge University Press NY NY [2] arxiv.org/abs/math/0610086