Conservative interaction of N three-dimensional internal waves controlled by the Navier-Stokes equations has been studied in , where the existence theorem is proved for a partial solution of the correspondent boundary-value problem via the Stationary Kinematic Euler-Fourier (SKEF) functions. However, an extreme sophistication of the partial solution derived with the help of experimental and theoretical programming in Maple doesn’t permit development of a general solution for propagation and interaction of wave groups of internal waves.
To overcome this challenge, a structural approach to this problem has been developed in , where the differential algebra of SKEF structures is studied and various applications of the SKEF structures for solving the Laplace equation and the Helmholtz equations are considered. Here, the method of Decomposition in Invariant Structures (DIS) is systematically developed to solve various kinematic (linear) problems of computing the Helmholtz decomposition.
The objective of the current paper is a continuation of the structural approach towards development of a cascade of dynamic (nonlinear) invariant structures, which permit the computation of the exact general three-dimensional solution for propagation and interaction of J wave groups with M internal waves per group.
The contents of this paper are following. Experimental Deterministic Scalar Kinematic (DSK) structures, which are the revised SKEF structures of , are systematically considered and complemented by theoretical DSK structures in Section 2. The experimental and theoretical DSK structures, which are used to solve the kinematic problem for scalar fluid-dynamic variables, generate experimental and theoretical Deterministic Vector Kinematic (DVK) structures, which are studied in Section 3 and are utilized to find vector variables of the kinematic problem. To compute scalar variables of the dynamic problem, experimental and theoretical Deterministic Scalar Dynamic (DSD) structures are defined with the help of the relevant DSK structures in Section 4 and the differential algebra of the DSD structures is treated, as well. In Section 5, experimental and theoretical DSK and DVK structures are combined to produce and develop experimental and theoretical Deterministic Vector Dynamic (DVD) structures for vector variables of the dynamic problem. The section concludes with the Helmholtz decomposition of directional derivatives of the DVK structures. The Helmholtz decomposition of the Navier-Stokes problem into the Archimedean, Stokes, and Navier problems is tackled in Section 6, where the general solution of the Stokes problem is expanded in the theoretical DSK and DVK structures. The general solution for deterministic chaos of J wave groups with M internal waves is completed in Section 7, where the dynamic variables are decomposed in the theoretical DSD and DVD structures. Section 8 deals with the orthogonal and non-orthogonal decompositions of the kinetic energy between wave groups and internal waves. Here, exponential oscillons and pulsons, which are novel nonlinear conservative solutions for propagation and interaction of various constituents of the kinetic energy and the total pressure, are introduced and visualized. The most interesting properties of the DSK, DVK, DSD, and DVD structures are discussed in Section 9 together with important properties of the general solution. Directions of further exploration of the deterministic chaos of internal waves are also outlined there.
2. Scalar Kinematic Structures
Scalar kinematic structures have been constructed in  and later developed in  as the SKEF structures. In a modified notation, the non-orthogonal normalized structural basis of four Deterministic Scalar Kinematic (DSK) structures becomes
where is an index of a structural wave group generated by the DSK structure , is an index of a harmonic wave, M is a total number of harmonic waves in the structural wave group, functional coefficients , which do not depend on the independent Cartesian coordinates (x, y, z) and time t, are functional amplitudes of a harmonic variable v(x, y, z, t). The functional amplitudes are normalized by
In the general solution, and in partial solutions.
In Equation (1), the DSK functions are defined as follows:
where are wave parameters related by the Pythagorean identity
In Equations (3), the dependent Cartesian coordinates
include x- and y-components of the velocity of a Cartesian frame moving with the mth wave and coordinates of the origin (x, y, z) = (0, 0, 0) of the motionless frame in the mth moving frame at the initial moment t = 0. In an upper domain,
a decay parameter
and, in a lower domain,
All wave parameters , , and of the DSK functions do not depend on (x, y, z) and t, as well.
If , the DSK structures (1) are reduced to the DSK functions (3) as
It is a straightforward matter to show completeness of the DSK structures (1) with respect to differentiation in (x, y, z) of any order   since
The first derivatives of the DSK structures in (x, y) are covariant since they are proportional to costructures in the x- and y-directions, respectively. The first derivatives of the DSK structures with respect to z are invariant because they are proportional to themselves.
The first spatial derivatives of the DSK structure are visualized in Figure 1. Differentiation in x and y moves the DSK structure from one corner of a differentiation rectangle to another one, while differentiation in z does not change a location of . For each invariant structure ,
Figure 1. A differentiation diagram of the first spatial derivatives of the DSK and DVK structures.
there is a costructure in the x-direction , which is located on the same horizontal leg at a distance given by the coefficient of a relevant derivative. There is also a costructure in the y-direction , which is located on the same vertical leg at a distance , which is provided by the coefficient of a correspondent derivative. Differentiation in z is invariant with the invariable coefficient since differentiation in z doesn't change the location of . Thus, the differentiation diagram demonstrates transformation of the DSK structures produced by differentiation, which is shown by blue arrows for derivatives in x, by green arrows for derivatives in y, and red arrows for derivatives in z. The length of arrows of derivatives in (x, y, z) are proportional to differentiation scales , which are shown with colors corresponding to those of arrows.
By a tedious calculation, it may be shown that the differentiation table of the first spatial derivatives of the theoretical DSK structures takes the following form:
Connections between the theoretical DSK structures and the experimental DSK structures (1) and the corresponding values of the sign parameters and are
The table of definitions (13) is complemented by a table of oscillation properties of the DSK costructures and the sign parameters of higher orders
The experimental DSK structures present a fixed notation of the corner points of the differentiation rectangle in Figure 1 and the theoretical DSK structures represent a floating notation of the vertices of the differentiation diagram. The differentiation table (12) shows equivalence of the theoretical DSK structures with respect to the spatial differentiation of the first order since each of the theoretical DSK structures may be used to obtain the experimental derivatives (11). The described comparison of the theoretical and experimental results shows a quadrality of the theory, i.e. there are four equivalent theoretical ways of explaining the experimental results. It may be shown that the quadrality of the theoretical DSK structures holds with respect to the second spatial derivatives, the Laplacians, and the first temporal derivatives, as well. For the aim of brevity, further theoretical results will be demonstrated mainly for the theoretical DSK structure that is sufficient to describe experimental results.
Using the theoretical table (12) of the first spatial derivatives, we find the second derivatives of the theoretical DSK structure in (x, y, z)
The differentiation diagram immediately explains invariance of the repeated second spatial derivatives. The second-order differentiation moves the DSK structure from any corner to another corner of the differentiation rectangle transforming to a costructure and then returns in back both in the x- and y-directions restoring the original value of . Similar to the physical oscillation, this effect of differentiation may be called a scalar structural oscillation.
Since the first derivative of the DSK structure with respect to z is invariant, the coefficient of the second derivative in z is positive because
The coefficients of the second derivative of in (x, y) are negative since the first derivatives in (x, y) are covariant.
In agreement with the differentiation diagram (Figure 1), the second derivative of in (x, y) becomes a costructure in (x, y), which is located at an opposite vertex of the differentiation rectangle to that of . The second derivatives of in (x, z) and (y, z) become the costructures and in the x- and y-directions, respectively, since differentiation in z is invariant.
Summation of the repeated second derivatives of the differentiation table (15) of the second spatial derivatives shows that the theoretical DSK structure is harmonic, i.e.
because of the Pythagorean identity (4), where the Laplacian
With the help of the chain rule, definitions (3) and (5), and the spatial derivatives (12), we find
Expanding (19) for and using (13) yields that the experimental DSK structures are also closed with respect to temporal differentiation of any order. In view of Figure 1, we observe that the first derivative of the DSK structure in t is a superposition of the costructures and in the x- and y-directions, which are located in the adjacent corner points of the differentiation rectangle. Whereas structural amplitudes of the costructures depend on products of the propagation velocities and the wave numbers in the (x, y) plane.
3. Vector Kinematic Structures
Four experimental Deterministic Vector Kinematic (DVK) structures are defined as gradients of the DSK structures (1) by
where the gradients are calculated in accordance with (11).
If , the DVK structures (20) are reduced to DVK functions, which are set as follows:
Application of the gradient to the theoretical DSK structures in the view of the differentiation table (12) results in the following definitions of the theoretical DVK structures:
Comparison of the definitions of the theoretical DVK structures (22) with the definition of the experimental DVK structures (20) again reveals the quadrality of theoretical formulas. The quadrality of the theoretical DVK structures is also confirmed by the tables of the divergences, the curls, the first spatial derivatives, the second spatial derivatives, the Laplacians and the first temporal derivatives. For the purpose of conciseness, further theoretical results will be shown mostly for the theoretical DVK structure that is sufficient for explanation of experimental results.
Calculation of the divergence of the theoretical DVK structure
shows that the theoretical and experimental DVK structures are divergence-free due to the Pythagorean identity (4).
Similarly, we find the curl of as follows:
Thus, the DVK structures are curl-free, as well.
Computation of a differentiation table of the first spatial derivatives of the theoretical DVK structures yields
In agreement with (25) and (12), the differentiation diagram of the DVK structures is identical to that of the DSK structures shown in Figure 1 because the differentiation tables of the DVK and DSK structures coincide up to the substitution for all i and m. This property of the invariant DVK and DSK structures may be called a scalar-vector structural invariance. As we will find out later, the scalar-vector structural invariance also holds for the second spatial derivatives, the Laplacians, and the first temporal derivatives.
Namely, we obtain the second spatial derivatives of in the following form:
The repeated second spatial derivatives of the DVK structures are invariant and the mixed second spatial derivatives are covariant, similar to the DSK structures.
Harmonicity of the DVK structures may be readily obtained by summation of the repeated second spatial derivatives (26) and the parametric identity (4) as
The DVK structures are closed regarding the spatial differentiation of any order due to (25). Completeness of the DVK structures with respect to temporal differentiation follows from
Again, the first temporal derivative of the DVK structures is a superposition of the first spatial derivatives in the x- and y-directions, while structural amplitudes are the same as for the DSK structures in (19).
4. Scalar Dynamic Structures
Define 16 experimental Deterministic Scalar Dynamic (DSD) structures as the following products of DSK structures and :
where and are indices of internal waves and M is a number of internal waves per wave group.
Therefore, 16 theoretical DSD structures are defined as products
where and are indices of wave groups.
Taking into account (12) and Equation (12) with substitution , we find first derivatives of in (x, y, z) as follows:
Expansion of (31) for all and usage of (13) yield completeness of the experimental DSD structures with respect to spatial differentiation of any order. The first derivative of the DSD structure with respect to z is invariant and with respect to (x, y) is covariant.
We then take second derivatives of (31) to obtain the repeated second derivatives of the theoretical DSD structures
and the following mixed second derivatives of the theoretical DSD structures:
Thus, Equations (32) demonstrate the invariance of the repeated second derivatives of the DSD structures in z and a partial invariance of the repeated second derivatives in (x, y). Equations (33) show the covariance of the mixed second derivatives in (x, y) and a partial covariance of the mixed second derivatives in (x, z) and (y, z).
Summation of the repeated second derivatives (32) establishes anharmonicity of the DSD structures since
By definition, the Laplacian of the product of scalar functions and
If the scalar functions are harmonic, then
where and are the gradients of f and g, respectively,
Substitution of the theoretical DSK structures and in (36) and simplification by (12) and Equation (12) with substitution gives the following Laplacian of the theoretical DSD structures:
where the scalar product of the DVK structures
Therefore, the non-orthogonality of the DVK structures is a cause of anharmonicity of the DSD structures.
Equations (31)-(34), (38)-(39) for the theoretical DSD structures have been verified by the differentiation tables of the experimental DSD structures, while each theoretical formula corresponds to a table of 16 experimental formulas.
5. Vector Dynamic Structures
Define 32 experimental Deterministic Vector Dynamic (DVD) structures as 16 products of four DSK structures with index m (m-DSK structures) and four DVK structures with index n (n-DVK structures) and 16 products of four m-DVK structures and four n-DSK structures :
So, 32 correspondent theoretical DVD structures are defined as products
First, we use the vector definitions of the derivative of the n-DVK structure in the direction of the m-DVK structure and the derivative of the m-DVK structure in the direction of the n-DVK structure
to compute with the help of (12), (25) with substitution , (25), and (12) with the subsequent directional derivatives through the theoretical DVD structures (41) in the vector form:
Alternatively, the component definitions of the directional derivatives (42)
after substituting (12), (15) with (15), and (12) with yield the directional derivatives of the theoretical DVK structures via the theoretical DSD structures (30) in the following component form:
The vector definitions, like (42), describe non-orthogonal vector decompositions of the fields of directional derivatives, anticommutators, and commutators, e.g. (43). Alternatively, a list notation of the component definition of any vector field
may be always written in the form of the orthogonal vector decomposition
Typically, the non-orthogonal vector decompositions of vector fields, for example (43), are simpler than correspondent orthogonal vector decompositions, for instance (45). However, the fundamental definitions of the differential operators of curl, divergence, and gradient are given through the orthogonal vector decompositions. We will consider both non-orthogonal and orthogonal vector decompositions of the vector fields generated by the directional derivatives, anticommutators, and commutators since the non-orthogonal vector decompositions do not form a closed set of equations. The list notation of the orthogonal vector decomposition will be used further for the aim of robustness of symbolic computation.
Summation and subtraction of the theoretical directional derivatives in the vector form (43) produce an anticommutator and a commutator of the directional derivatives of the theoretical DVK structures in terms of the theoretical DVD structures, respectively:
The directional derivatives (43), the anticommutators and the commutators (48) of the DVK structures in the vector form are decomposed in the DVD structures (41). The directional derivatives, the anticommutator, and the commutator of the DVK structures in the component form (45) are expanded via the DSD structures (30).
Second, the component definition of the gradient of the dot product of the m-DVK structure and the n-DVK structure becomes
We substitute (12), (15) with, (15), and (12) with for the theoretical DSK structures in (49) and simplify it to find the following component form of the gradient:
Subtraction of the directional derivatives of the theoretical DVK structures (45) from (50) gives a vector form of the gradient of the dot product
The gradient of the dot product of the DVK structures in the component form (50) is decomposed through the DSD structures (30). The gradient of the dot product of the DVK structures in the vector form (51) is expressed via the anticommutator of the directional derivatives of the DVK structures (48).
Third, the component definition of the cross product of two gradients and (37) yields
We then find the following definition of the curl of the cross product of and
Usage of the harmonicity conditions for f and g
results in the component definition of the curl of the cross product of two harmonic gradients and:
Substitution of the theoretical DSK structures and in (55) and simplification by (12), (15) with, (15), and (12) with returns
Eventually, we add to (56) and subtract from (56) the relevant directional derivatives of (45) to get
The curl of the cross product of the DVK structures in the component form (56) is expanded in terms of the DSD structures (30). The curl of the cross product of the DVK structures in the vector form (57) is expressed via the commutator of the directional derivatives of the DVK structures (48).
Fourth, solving (51) with respect to the anticommutator yields
So, the scalar potential of the Helmholtz decomposition (58) of the anticommutator field is the dot product of the m-DVK and n-DVK structures.
Finding the commutator from (57) gives
Thus, the vector potential of the Helmholtz decomposition (59) of the commutator field is equal to the negative cross product of the m-DVK and n-DVK structures.
Multiplication of the Helmholtz decompositions (58)-(59) by 1/2, addition, and subtraction produce
is a scalar potential and
is a vector potential of the Helmholtz decompositions (60)-(61) of the fields of the directional derivatives of the m-DVK and n-DVK structures.
It was also verified that the theoretical expansions in terms of the DVD and DSD structures (43), (45), (48), (50)-(51), (56)-(61) completely describe the correspondent expansions via the experimental DVD and DSD structures.
6. The Stokes Field
Internal waves of a Newtonian fluid with a constant density and a constant kinematic viscosity in a field of gravity are governed by the momentum conservation law
and the mass conservation law 
The quasi-scalar Dirichlet problem for the Navier-Stokes Equations (64)-(65) may be set for the z-component of the velocity field on the upper and lower boundaries of the upper domain U (6)
and the upper and lower boundaries of the lower domain L (8), e.g.  ,
An admissible decomposition of the boundary function in DSK boundary structures will be considered in Section 6 later. The conditions at infinity are satisfied due to (7) and (9). Configuration of the upper and lower domains of internal waves, which are generated through the boundary conditions (66)-(67) by surface waves propagating in a generation domain, is shown in Figure 2.
From viewpoint of the Fundamental Theorem of Vector Analysis  the Dirichlet problem (66)-(67) for the Navier-Stokes Equations (64)-(65) may be treated as a construction of the Helmholtz decomposition for the Archimedean field, the Stokes field, and the Navier field:
The Archimedean, Stokes, and Navier fields are decomposed using the correspondent scalar potentials and as follows:
where stands in fluid dynamics for the hydrostatic pressure, for the kinematic pressure, and for the dynamic pressure.
Summation of (68)-(71) yields the Helmholtz representation of the momentum conservation law (64)
where, which is termed in fluid dynamics a total pressure, is a scalar potential of the superposition of the Archimedean, Stokes, and Navier fields.
The problem of finding the scalar potential of the Archimedean field has a general solution 
Figure 2. Configuration of the upper domain (6) and the lower domain (8) of the Dirichlet problem (66)-(67) for the Navier-Stokes Equations (64)-(65).
where is a reference pressure of the Cauchy integral .
The problem of calculating the kinematic velocity field and the scalar potential of the Stokes field is equivalent to the following quasi-scalar Dirichlet problem for the Stokes equations:
It will be called afterwards the Stokes problem.
The problem of computing the scalar potential of the Navier field for the given field
will be later referred to as the Navier problem.
Decomposition of the Navier-Stokes problem (64)-(67) into the superposition of the Archimedean problem (68), (71), (73), the Stokes problem (69), (71), (74), and the Navier problem (70), (71), (75) is justified by the structure of the total pressure (72). From the mathematical standpoint, an exact solution of the Stokes problem (74) requires a balance of terms, which are linear in, but an exact solution of the Navier problem (75) necessitates a balance of terms, which are quadratic in. As to the fluid-dynamic viewpoint, the Navier problem models the inertial regime of intermittency of flows at high Reynolds numbers and the Stokes problem describes the viscous regime of generation/decay of flows at low Reynolds numbers.
To construct general solutions for the kinematic velocity and the kinematic pressure of the Stokes problem (74), we start from the following definitions of anormalized DSK structures:
where the DSK functions are given by (3) and functional amplitudes are wave parameters, which do not depend on (x, y, z, t).
Using a quadratic norm, which is given by a weight
the anormalized amplitudes are reduced to the normalized amplitudes of (1) by
Substitution of (78) in definitions (76) yields
Relationships between anormalized DVK structures and the normalized ones (20) are similar due to the scalar-vector structural invariance since
The general solution for the kinematic velocity is decomposed in four wave groups with M internal waves per group using the anormalized DVK structures in the following form:
In agreement with (81) and (80),
and the general solution in the normalized DVK structures after changing the order of summation reads
The velocity field (84) delivers the general solution of a direct problem of construction of the velocity field with given weights of quartets of the experimental DVK structures.
For an inverse problem of reconstruction of the velocity field from data of a laboratory experiment, a linear norm, which is set via a weight
may be used to transform the quartet weights into quartet probabilities
So, the general solution of the inverse problem may be written as follows:
where equiprobability of the experimental DVK structures, which are shown in Figure 1, is explained by quadrality of the theoretical DVK structures.
In agreement with the Cauchy integral, we construct the field of kinematic pressure via general terms in the following form:
Substitution of the temporal derivative (19) gives the general solution for of the direct problem
and the general solution for of the inverse problem
Substituting the velocity field and the pressure field in the Partial Differential Equations (PDEs) of the Stokes problem (74) and using the temporal derivative of the DVK structures (28), the Laplacian of the DVK structures (27), the gradient of the DSK structures (22), and the divergence of the DVK structures (23) yield that PDES of (74) are satisfied identically both for the direct solutions (84), (90) and the inverse solutions (88), (91).
To construct the DSK boundary structures, primarily, we set the DSK boundary functions for and as
So, the DSK boundary functions (92) are compatible with the DSK functions (3), where are wave numbers of boundary waves,
are the boundary Cartesian coordinates, are propagation velocities, and are initial coordinates of the boundary waves.
In accordance with (1), the DSK boundary structures are defined by
where functional amplitudes are normalized, i.e.
All boundary parameters do not depend on (x, y, z, t).
Eventually, we define the boundary function of (66)-(67) via the DSK boundary structures in the following form:
In agreement with (22), the z-component of the general solution (84) for the velocity field becomes
If and only if
then the dependent Cartesian coordinates, the scale multiplier in the z-direction, the theoretical DSK function, and the theoretical DSK structure are related with the correspondent boundary variables as follows:
In this case, the z-component of the velocity field (97) is reduced to the boundary function (96) and the Dirichlet boundary condition of the Stokes problem (74) is fulfilled exactly.
7. The Navier Field
A rectangular summation of all vector elements of a square matrix
with, which is represented by a double sum, may be always reduced to a triangular summation of diagonal terms and pairs of terms above and beneath the diagonal by the following formula:
Indeed, the triangular summation is also valid for a square matrix with scalar elements.
Primarily, we expand the velocity field (84) into J wave groups
where the velocity fields and of the ith and jth wave groups, respectively, are
Substitution of expansions (102) in the Navier field (70), reduction by the rectangular summation, and transformation by the triangular summation (101) in wave groups with substitutions yields
The directional derivative that describes propagation of the ith wave group in view of (103) may be written using the rectangular and triangular summations through the DVK structures as follows:
Likewise, the directional derivatives that express interaction of the ith and jth wave groups via the DVK structures become
Secondly, we start potentialization of the directional derivative (105) substituting the following Helmholtz decompositions:
which are derived from (60) and (61) by substitutions and. The result is
To proceed with potentialization of the directional derivatives (106), the necessary Helmholtz decompositions of the relevant directional derivatives
are obtained by substituting in (60) and (61) and. Substitution of (109) in (106) gives
We then use (108) and (110) to get from (104) the Navier field in the potentialized form
Thirdly, the Navier field is reduced using conversion of the triangular summation in wave groups into the rectangular summation as
By decomposition of rectangular sums into products of sums and conversion of the triangular summation with respect to internal waves into the rectangular summation, is transformed into the following completely rectangular form:
Comparison of (113) with the definition of yields the dynamic pressure of the direct problem
and, together with (86), the dynamic pressure of the inverse problem
Using connection (38) between the Laplacian of the DSD structure and the scalar product of the DVK structures, we get
Equation (116) demonstrates a mathematical meaning of the dynamic pressure that is a source of the stationary diffusion of the weighted superposition of all DSD structures with a diffusion coefficient.
Potentialization of the Navier field is possible since internal vortex forces, which are described by the vector potentials of the Helmholtz decompositions (107) and (109), counterbalance each other, in agreement with Newton’s third law. On the contrary, external potential forces, which are expressed via the scalar potentials of (107) and (109), superpose together to form the gradient of the dynamic pressure, in accord with Newton’s second law.
Transformation of the Navier field (113) into the following product of two double sums
shows that the dynamic pressure and the kinetic energy have the same magnitude as since in agreement with (71), (113), (114), (117), (22), and (12), an orthogonal decomposition of in the DSD structures becomes
A non-orthogonal distribution of the kinetic energy between the wave groups, which is obtained using the triangular summation with respect to wave groups, has the following form:
where is the kinetic energy of propagation of the ith wave group and is the kinetic energy of interaction between the ith and jth wave groups.
We then find a non-orthogonal distribution of the propagation energy between internal waves, making the triangular summation (101) in regard to internal waves, to get
where is the kinetic energy of propagation of the mth internal wave of the ith wave group and is the kinetic energy of interaction between the mth and nth internal waves of the ith wave group.
Similarly, the triangular summation of gives the following non-orthogonal distribution of the interaction energy between internal waves:
where is the kinetic energy of interaction between the mth waves of the ith and jth groups, is the kinetic energy of interaction between the mth wave of the ith group and the nth wave of the jth group, and is the kinetic energy of interaction between the nth wave of the ith group and the mth wave of the jth group.
Substituting the computed general terms (120)-(121) in (119), we conclude with the non-orthogonal decomposition of the kinetic energy
into the following seven constituents. First, the cumulative kinetic energy of propagation of J wave groups
Second, the cumulative kinetic energy of interaction between J wave groups
Third, the cumulative kinetic energy of propagation of the mth internal waves of J wave groups
Fourth, the cumulative kinetic energy of interaction between the mth and nth internal waves of J wave groups
Fifth, the cumulative kinetic energy of interaction between the mth internal waves for J distinct wave groups
Sixth, the cumulative kinetic energy of interaction between the mth wave of the ith group and the nth wave of the jth group for J distinct wave groups
Seventh, the cumulative kinetic energy of interaction between the nth wave of the ith group and the mth wave of the jth group for J distinct wave groups
The kinetic energy and the cumulative kinetic energy of propagation demonstrate the pulsatory topology since the exact, nonlinear, three-dimensional (3-D), periodic solutions of the Navier-Stokes Equations (64)-(65) are strictly positive for all (x, y, z, t), and may be called exponential pulsons of elevation. Topology of the exponential pulsons is the same as of the solitons on shallow water, the solitary waves on shallow water with uniform and linear vorticity  , the solitary waves generated by crossed electric and magnetic fields , and the pulsatory waves of the Korteweg-de Vries equation ,
Figure 3. The kinetic energy and its constituents: (a)—(118); (b)—(125); (c)—(126) with wave parameters given by (130).
which are exact nonlinear solutions for propagation and conservative interaction of aperiodic one-dimensional (1-D) waves. In view of (118), the corresponding solutions for the dynamic pressure and its propagation constituent are given by exponential pulsons of depression since they are strictly negative for all (x, y, z, t).
Figure 4. The cumulative kinetic energy of interaction of distinct wave groups: (a)—(127); (b)—(128); (c)—(129) with wave parameters given by (130).
The cumulative kinetic energies of interaction exhibit the oscillatory topology since the exact, nonlinear, 3-D, periodic solutions of the Navier-Stokes Equations (64)-(65) alternate signs for various (x, y, z, t) and might be called exponential oscillons. Topology of the exponential oscillons is the same as of the solitons on deep water , which are exact nonlinear solutions for propagation and conservative interaction of aperiodic 1-D waves. The related solutions for the interactive constituents of the dynamic pressure and the total pressure are represented by the exponential oscillons, as well. Animations of visualizing propagation and conservative interaction of the exponential pulsons and oscillons, which are various superpositions of the time-dependent, 3-D, periodic DSD structures (29)-(30), confirm conservation of both the pulsatory and oscillatory topologies for all (x, y, z, t).
The most interesting properties of the scalar and vector kinematic structures that are studied in Section 2 and Section 3, respectively, are the scalar and vector structural oscillations, the scalar-vector duality, the quadrality of the theoretical DSK and DVK structures, and the equiprobability of the experimental DSK and DVK structures. The differentiation diagram in Figure 1 is actually an image of a mathematical two-dimensional dice, where all vertices have exactly the same probability despite the various scales of differentiation in the x- and y-directions.
The success of the scalar and vector dynamic structures treated in Section 4 and Section 5, correspondingly, is explained by their exact similarity to the Navier-Stokes Equations (64)-(65) as the momentum conservation law (64) represents a superposition of kinematic terms and dynamic terms with an algebraic nonlinearity. Further attractive features of the DSD and DVD structures are the absence of a formal restriction on the functional amplitudes, the structural amplitudes, and the quartet weights and the global convergence the DSD and DVD structures for all (x, y, z, t). For instance, the general solution of fluid dynamics in the Boussinesq-Rayleigh series converges only locally and the radius of convergence diminishes with the Reynolds number .
The large number of the experimental DSD and DVD structures, 16 and 32, respectively, is required for completeness of the expansions of the directional derivative of the global velocity and the gradient of the total pressure since otherwise the vortical forces produced by the vector potential (63) of the Helmholtz decomposition may not compensate each other. Therefore, the experimental and theoretical programming in Maple that facilitates computation and verification of the numerous large arrays of scalar and vector terms is essential for development of the general solution for J wave groups with M internal waves. Since the exact general solution is not affected by viscous dissipation it may serve as the 3-D model of conservative propagation and interaction of kinetic energy of internal waves in ocean and atmosphere via the exponential oscillons and pulsons visualized in Figure 3 and Figure 4.
Research in the functional amplitudes of the DSK structures started in  has revealed existence of three types of wave lattices: 1) slanted up and down, 2) stepped up and down with rectangular steps, and 3) stepped up and down with smooth steps. It is interesting to explore the effect of the wave parameters on the Eulerian and Lagrangian properties of the DSD structures, as well. By the well-known conjecture , chaos of physical systems is explained by superposing of a large number of deterministic solutions. Therefore, it also looks appealing to study stochastic properties of the general solution.
The support of CAAM is gratefully acknowledged. The author would like to thank reviewers and editor for valuable comments which have significantly improved the paper.
 Miroshnikov, V.A. (2017) Harmonic Wave Systems: Partial Differential Equations of the Helmholtz Decomposition. Scientific Research Publishing, USA.
 Miroshnikov, V.A. (2002) The Boussinesq-Rayleigh Approximation for Rotational Solitary Waves on Shallow Water with Uniform Vorticity. Journal of Fluid Mechanics, 456, 1-32.
 Miroshnikov, V.A. (1995) Solitary Wave on the Surface of a Shear Stream in Crossed Electric and Magnetic Fields: The Formation of a Single Vortex. Magnetohydrodynamics, 31, 149-165.
 Miroshnikov, V.A. (2014) Interaction of Two Pulsatory Waves of the Korteweg-de Vries Equation in a Zigzag Hyperbolic Structure. American Journal of Computational Mathematics, 4, 254-270.