Turbulent wind waves on a water current

Abstract. An analytical model of water waves generated by the wind over the water surface is presented. A simple modeling method of wind waves is described based on waves lengths diagram, azimuthal hodograph of waves velocities and others. Properties of the generated waves are described. The wave length and wave velocity are obtained as functions on azimuth of wave propagation and growth rate. Motionless waves dynamically trapped into the general picture of three dimensional waves are described. The gravitation force does not enter the three dimensional of turbulent wind waves. That is why these waves have turbulent and not gravitational nature. The Langmuir stripes are naturally modeled and existence of the rogue waves is theoretically proved.


Introduction
Great variety of wave types in liquids have been studied since Newton, Bernulli, Euler, Lagrange, Laplace, Cauchy, Airy, Lamb, Poincare, Lyapunov and many others.The role of the waves in the fluid dynamics theory is well understood.Remarkable advantages for the hundreds of years have been achieved in understanding of the wave properties.The main directions of research in the fluid wave motions cover (though not limited to) gravitational waves, internal waves, ship waves, tidal waves, acoustic, shock, solitary, convective waves, etc.
The investigation methods are not as diverse as the waves themselves.The basic tools of investigation are the statistics methods, approximation calculus and classical Euler equations (Lamb, 1932;Stocker, 1957).Not many specialists use the model of viscous fluid in their investigation (see, for example, Sajjadi, 1999).Sometimes Lagrange variables Correspondence to: P. B. Rutkevich (peter@d902.iki.rssi.ru)(Monin, 1999), and the conformal mapping theory are applied for the waves modeling (Nekrasov, 1961).
However, despite such an intense intellectual attack there still exist issues not well understood in the wave motions theory.One of them is the theory of wind waves on the water surface which is usually investigated based on the ideal fluid theory and neglecting the viscosity.In other cases numerical, experimental, semiempirical, probabilistic, statistic methods are used in investigations.All researches assume that the wind waves are gravitational.Therefore, nowadays there does not exist any analytical model of the wind waves with account of the fluid viscosity.
The idea of modeling of the wind waves neglecting the influence of viscosity between the water surface and the air wind seems hopeless since the wind is the generator of these waves.However, there is a strong reason to neglect the viscosity, and it is due to the internal structure of the Navier-Stokes equations.Indeed, in order to take into account analytically the wind influence on the water surface one has to assume that the location of the main masses of air driven by the wind is high enough as compared to the surface wavelengths.As a matter of fact this condition is always satisfied in nature.To model the wind waves in viscous fluid one should also assume that the waves do not affect the water current in the deep stream, existing in the sea.In other words undisturbed deep water flow in the first approximation should be located infinitely far from the water surface.And again, in the first approximation the fluid should be considered infinitely deep.And the reason why the wind waves cannot be described on the base of the Navier-Stokes equations is that the Navier-Stokes equations do not have a stationary solution for the fluid flow in the half-space under influence of stationary tangent stresses on its surface (Kochin et al., 1963;Kondrat'ev et al., 1990).
The idea that conventional Navier-Stokes equations are not complete has been mentioned before.For example, (Birkhoff, 1960) mentions that the NavierStokes equations Fig. 1.Air photograph of the sea surface at strong wind.The light Langmuir stripes and the dark stripes of wind wave troughs can be observed (Titov, 1969).
describe well enough incompressible flows of viscous fluids, to which belong usual gases and liquids at speeds considerably smaller speed of sound.However description of turbulent flows requires accurate use of statistical analysis at high mathematical level.Therefore, instead of use of the NavierStokes equations at formulation of obvious boundary problems (for example, the problems of stationary current) it is necessary to address directly to the physical validity for statement of the appropriate boundary problem.In our model we use the turbulent resistance to resolve this problem.
General structure of the surface waves generated by the wind can be understood even from a photo image of the sea surface Fig. 1 (Titov, 1969).One can notice distinct directions of bright Langmuir stripes and the direction at which the waves do not propagate.Motionless waves correspond to the zero wavelength in the corresponding direction.Therefore, there is a closed curved line the tangent to which is aligned in the direction of no wave propagation (since the waves are confined such a curve does not go to infinity).The chords lengths of such a line counted from the point of contact are equal to the length of wind waves.It is appropriate to call this curve as wavelength diagram along the azimuth of the chosen diagram chord.It is appropriate to call this curve azimuth wind waves velocity hodograph along the azimuth of the chosen hodograph chord.
The Langmuir stripes represent motionless waves since they are in fact motionless on the water surface.These waves "propagate" (i.e. are oriented) in the direction orthogonal to the Langmuir stripes.It is clear then that along the orientation of the motionless waves (orthogonal to the Langmuir stripes), the velocity of the wind eaves equals zero.
The velocity hodograph and the wavelength diagram of the wind waves really exist and will described later in this paper.The model of the wind waves is based on the Navier-Stokes equation with isotropic turbulent resistance.The presented analytical model of turbulent wind waves over the water current is remarkably simple, since the waves are described through elementary relations.Moreover, the model predicts existence of the new waves: motionless and barchan's waves.The model also describes the choppy sea the kind of waves that up to now have no analytical description.

Basic dispersion relations
As it has already been mentioned above, Navier-Stokes equations do not have a stable solution in a semi-space with tangent stress on the surface.However, if the deep water current and high air flow have different directions there should exist a stable solution, responsible for wind-generated waves on the surface.
We start from the Navier-Stokes equations with turbulent resistance (Gill, 1982;Lamb, 1932;Kondrat'ev et al., 1990;Struminskii, 1969): where V is the fluid velocity at infinity, and κ is turbulent resistance of the medium.Equation (1) has to be applied twice, first for the air in the upper half-space, and second for the water in the bottom half-space, and satisfy the boundary conditions between these two regions.
Introducing Cartesian coordinates Ox y z , with axis Oz directed against the gravity force g, and axis Ox directed along the direction of the constant air velocity vector V 1 at z →∞.Density, kinematic viscosity and turbulent resistance of the air denote as ρ 1 , ν 1 and κ 1 , respectively, finally velocity and hydrodynamic pressure of the air as v 1 (x , y , z , t ) and p 1 (x , y , z , t ).Similarly, water is described by the parameters ρ, ν and κ, (ρ ρ 1 ), velocity v (x , y , z , t ), pressure p (x , y , z , t ) and at z →−∞ moves with constant velocity V 0 =V 0 (i cos α + j sin α), where i and j are the unit vectors along Ox and Oy , respectively.The angle α is counted from vector V 1 towards V 0 (−π≤α≤π ).The water-air interface plane is z =ζ (x , y , t ).Therefore, to determine the parameters v 1 , p 1 , v , p and ζ we have Eq.( 1) with isotropic turbulent resistance for air: and water: supported by the fluxes continuity conditions at the air-water boundary z =ζ (i.e.equality of their velocities and stresses along the tangent direction), kinematic condition, that the interface points do not leave this interface and conditions at infinity.This system of equations with the mentioned boundary conditions has a stationary solution at the boundary, At large wind V 1 or water V 0 velocities the solution (2) loses its stability.Therefore we are looking for the solution of the problem different from (2): {x, y, z} = {x , y , z } κ/ν, t = t κ All the functions marked with tilde (∼) (except ζ ) depend on x, y, z and t.New variables (3 will convert initial equations into a dimensionless form, After linearization, the initial equations and limiting conditions (2) for new variables V , v, P , p and ζ satisfy the following equations: Omitting the terms of the order of O(ε), and having calculated ũ, ṽ and w in order to determine P , Ũ , Ṽ , W we obtain a nonhomogeneous boundary problem.We are not going to investigate this problem since the water waves are more interesting then those in the air over the water surface.We are looking for the solution of the matrix { ũ, ṽ, w, p, ζ } in the form where i = √ −1.After substitution into (4) one obtains the following differential equations for unknown variables u, v, w and p (for simplicity we used a vector representation for three variables and three equations): Here the prime stands for derivative with respect to z. Multiplying the first component of the vector row (6) by im, the second by in and add together the first two components of the obtained relation.Then taking into account the continuity equation we obtain Now calculate the derivative p and put it equal to the same derivative from the third row of the vector Eq. ( 6).This gives the equation for determination the component w (z<0): Multiplying the first component of the vector equation ( 6) by n, the second by m and subtracting the first from the second yields (z<0) The limiting conditions for the Eqs.( 8), ( 9) are obtained by substitution of the solution (5) into the limiting conditions (4) at ε=0: So (nu−mv) =0 (z=0) and nu−mv=0 (z=−∞).Excluding the parameter ζ from the two remained limiting conditions and using the continuity Eq. ( 6) and the Eq. ( 7), one obtains the limiting conditions for the Eq. ( 7), that is formulate the boundary problem for determining the amplitude of the vertical velocity w(z): The amplitudes of u, v, p and ζ are not necessary for obtaining the dispersion relation.These amplitudes can be derived later from the nonhomogeneous problem for the Eq. ( 9) after the component w is found.The dispersion relation connects real wavenumbers m and n with the complex frequency c, it is derived from the spectral problem ( 8),( 10), that can be solved approximately using the following idea.The equation coefficients (8) on account of Eq. ( 2) are linear functions of exp(z), (z<0).The parameter z=z √ κ/ν for finite z is large, since the turbulent resistance is strong during the wave formation.So the exponents in the Eq. ( 8) coefficients are the pronounced interface layers.For derivation of the dispersion relation these layers are taken into account by replacement the exponents with units.In the boundary conditions (10) which define the dispersion relation, the influence of the layers is taken into account more precisely.The so approximated spectral problem (8), (10) takes the form Solution of the Eq. ( 11) vanishing at z→−∞ has the form where ξ 1 , ξ 2 are the roots of characteristic equation: with negative sign of real parts.Substituting velocity (12) into the boundary conditions ( 11) at z=0, one gets a system of two linear homogeneous algebraic equations for variables C 1 and C 2 .The determinant of this system is must be equal to zero, thus giving following dispersion relation is obtained The most important is the situation when ξ 1 =ξ 2 since in the case of multiple roots of the Eq. ( 13) the disturbances (12) penetrate into water deeper than the disturbances of two different root.At ξ 1 =ξ 2 instead the relation ( 12) one has to use w=(C 1 +C 2 z) exp(ξ 1 z), (z < 0), which maximal value modulating the oscillations is attained bellow the surface, rather than on the surface as in the case of relations ( 12).The roots of Eq. ( 13) are equal at (q+k 2 ) 2 =4(qk 2 +iδB).Therefore the parameter Taking into account the parameter q definition (11) one obtains the complex oscillation frequency (5): (c=c r +ic i =b−i±(i−sgn(B)) √ ∈δ|B|) (sgnx=1 at x>0, sgn(x)=−1 for x<0).Hence the frequency c r and the growth rate c i of the waves (5), corresponding to the dispersion relation at Here one has to use simultaneously either the upper or the lower signs.One can see, that there is no gravity force g in the dispersion relations (15), thus the waves described by relation ( 15) can be called as non-gravitational.Analysis presented in the rest of this work is based on the relation (15).

Length and velocity of the wind waves
Introducing the azimuth θ of the wave vector k={m, n} with respect to the wind velocity V 1 , one obtains m=k cos θ, n=k sin θ, and the relations (11) for b and B can be rewriten as The waves defined by the expressions ( 17) can be called as κ-waves, since they have turbulent and not gravitational nature (i.e. the gravity force g does not appear in the formulas ( 17)).To describe the κ-waves let us abandon the tradition and investigate their properties not as function of the wave number k, but as a function of the growth rate, or, more precisely, as a function of parameter K=c i +1.The K dependence of k is presented in ( 17).The values K<1 correspond to the decaying waves, the value K=1 corresponds to the neutral waves, and at K>1 there exist non decaying waves.In order to obtain the dimensional velocity λK from Eq. ( 17) and the wavelength λ K of κ-waves the ratio c r /k sould be multiplied by the velocity scale √ κν, and the dimensional wavelength 2π/k by the length scale √ ν/κ: The formulas ( 18) and ( 19) describe the wave front velocity and the wave length at any direction θ with respect to the wave velocity.The formula ( 18) is convenient for description of the waves velocities at large values of the absolute value of dimensionless characteristics K of the growth rate c i .It follows, for example, from formulas ( 18) and ( 19) that at large values of |K| the decaying waves and the non decaying waves move at the same velocity symmetrically with respect to direction θ=(π/2)−γ , in which the velocity of the waves is maximal.Note that the wavelengths (19) of these waves are not symmetric with respect to this direction.If the parameter |K| does not increase then instead of (18) for the velocities description κ-waves it is more convenient to use another expression.Let us use in (18) the definition of the parameters A, C and the angles β, γ from (16).Then Here | K | is the maximum velocity of the κ-waves with respect to the azimuth θ at the fixed value K, corresponding to waves propagation in the direction θ=θ K .At K=0 the growth rate c i =−1.It follows from ( 15) and ( 16)that in this case ether k=0 (there is no wave motion), or θ=β or θ=β+π.The uncertainty zero divided by zero at K=0, θ=β, θ=β+π in the formula ( 19) is prescribed the value zero, since any κ-waves at K =0 at the directions θ=β, θ=β+π have zero lengths.The direction in which the developed water waves have zero lengths (that is not seen) are well known from the air photography observations (Titov, 1969).This direction corresponds to the orientation of the wind waves troughs.Formulas ( 16) define the β dependence of this orientation on the values and directions of the wind and water flow velocities.Note that the directions θ=β±π/2 of the waves propagation of the maximal azimuth wavelength and the directions θ=β, θ=β+π in which the wavelength equals zero do not depend on the growth rate.These directions do not depend even on air and water physical properties.At K=0 the velocity (21) and the length (19) of the κ-waves go to infinity.There takes place an internal resonance of the growth rate of the κ-waves with the turbulent resistance frequency (in this case c i =−κ).If one expresses the parameter K 2 with the use of L K and the formula (19), then the approach of the parameter K to zero means just increase of the maximal wind waves length.

Wave lengths diagram
The formula λ K =L K | sin(θ−β)| (see Eq. 19) describing the wind wave lengths for all values of the angle θ, means that the wave lengths diagram is a pair of equal circles.Let us mark an arbitrary pole O on the plane z=0 and fix a counting out direction V 1 of the azimuth θ.For simplicity the rays on the z-plane are marked by their azimuths.The circles of the wind wave lengths diagram touch the straight line θ=β, θ=β+π (Fig. 2).
The centers of the wave lengths diagram circles C and C are located on the straight line θ=β± 1 2 π at distance 1 2 L K from the diagram circles point of contact.
It follows from the formula ( 16) that the angle β does not depend of the wind waves growth rate and the tan β does not depend of the parameter d (see Eq. 2) that is, does not  depend neither of the viscosity, nor of their turbulent resistances.That is why the β angle is defined only by the air and water velocities V 1 , V 0 and the angle α between the vectors V 1 and V 0 .It means that the straight lines θ=β, θ=β+π and θ=β± 1 2 π, characterizing all (by the growth rate) wave length diagrams are obtained by measurements and simplest calculations.
The most important azimuth β(±π) property is that the wave length at this direction equals zero.This follows from the formula (19) and also from the wavelength diagram (Fig. 2).The zero wavelength direction θ=β is easily seen on the water waves air photograph at strong wind (Fig. 1).

The wave lengths similarity theorem
It was mentioned above that at K=1 the wind waves are neutral.In this case it follows from the formula (19) that the wind wavelengths at any value of parameter K are similar to the lengths λ 1 =L 1 | sin(θ−β)| of the neutral waves for any fixed direction θ.The stated theorem essentially takes into account the fact that the angle β does not depend of the K parameter.
The circle λ 1 =L 1 | sin(θ−β)| does not depend of the growth rate.That is why it also enters the skeleton of the wind wavelengths diagram.
Note also that the wind wave lengths have all the properties of circle chords.In particular the sum of wave lengths squared in two perpendicular directions equals to the squared maximal waves length:

Wind waves velocities azimuth hodograph
From formula λK = k sin(θ+θ K ) (see Eq. 21) one can see, that the wind waves velocities diagram is also a circle.We are going to call the wind waves velocities diagram azimuth hodograph of the wind waves.

Motionless wind waves
The azimuth hodograph circle of the wind waves touches the straight line θ=−θ K , θ=−θ K +π.Therefore at the direction θ=−θ K the wind waves velocity equals zero.
In the orthogonal direction θ=(π/2)sgn(K)−θ K (see Eq. 21) the zero velocity troughs are located.We will call the waves with zero velocity as motionless waves.In the motionless wave's troughs foam, driftweed and other minor material from the water surface is accumulated.Thereby the motionless waves are marked on the water surface as motionless light stripes called "the Langmuir stripes".The motionless waves are seen on the water surface air photographs at strong wind (Fig. 1).We do not consider in this paper the Langmuir circulations as it was made in the paper (Craik and Leibovich, 1976), though it should be noted, that our approach is more preferable, since we use both the water current and the wind currant and not just the surface stresses as in the paper (Craik and Leibovich, 1976).
If put aside all the damped neutral and non damped wind waves and leave only the motionless waves, the water surface will become crimped in the direction orthogonal to the Langmuir stripes.The amplitude of such a goffer diminishes at damped waves and increases up to total destruction at non damped waves.The goffer of the neutral wind waves is frozen-in the other waves until the wind blows.
As a result of the wind waves dynamic interaction the motionless Langmuir stripes are described here for the first time.
The center of the wind waves velocities hodograph is situated of the strait line θ=(π/2)sgn(K)−θ K at the distance 1 2 K from the point of contact of the hodograph with the line θ=−θ K , θ=−θ K +π (Fig. 3).
The axis θ=(π/2)sgn(K)−θ K of the wind waves velocities hodograph rotates at the parameter K change (see the formulas (21), defining the trigonometric functions of the angle θ K ).That is why there is no wave velocity similarity theorem.
Figure 3 describes one of the possible wind waves velocities azimuth hodographs with the center at the point C K .The vector K in Fig. 3 equals the maximal wind waves velocity for the given parameter K value.The vector OA in Fig. 3 equals the wind waves velocities at the direction θ (for the given parameter K value).The vector V 1 equals the wind velocity.All the azimuths are counted out from the vector V 1 on the z-plain.The velocity hodograph diameter directed along the maximal wind waves velocity K is called the main diameter of the hodograph.

σ -hodograph
Let us find the wind waves velocity along the direction θ=β.Let in the formula (18) θ=β and take into account formulas ( 16).Then Here velocity σ ≡λ K (β) is the velocity of the zero length waves that is the velocity of the wave ripples in the direction θ=β.σ -hodograph is the circle which main diameter equals to σ and is directed along the azimuth θ=β at 0≤α≤π and along the azimuth θ=β+π at (−π)<α<0.The σ -hodograph circle takes the form: Let us find the minimal value of the velocity K with respect to parameter K using the formulas ( 21) and ( 16).It turns out as a result that the minimum is obtained at and min This means that the minimal wind waves velocity with respect to K equals |σ |.
But the main property of the wind waves velocity of zero length is that the velocity σ does not depend of the growth rate1 .That is why since the value K σ corresponding to the main diameter of the σ -hodograph (see Eq. 24) belongs to a set of values of the parameter K, the main diameter σ of the σ -hodograph is a chord of any velocity hodograph of the wind waves.Now it is possible to formulate the main theorem of geometric dynamics of the non gravitational wind waves: The locus of all the centers of the main diameters of the wind waves velocity hodographs is a strait line passing perpendicular to the vector σ via its middle.The locus of all the main diameters ends of the wind waves velocity hodographs is a strait line passing perpendicular to the vector σ via their ends.

Consequence.
Since the vector σ is oriented along or against the ray θ=β, perpendicular to the vector σ centers line and the main diameters ends line are parallel to the axis θ=β±π/2 of the wavelengths diagram.
Therefore, the wavelengths diagram, the main σhodographs diameter and the two lines: the line of the centers of the main velocity hodographs diameters and the line of the ends of the same main diameters enter into the general skeleton of the wind waves dynamics.
The physical meaning of the theorem and its consequence is in the fact that the main σ -hodograph diameter can be measured by the value and direction2 .

Absolutely motionless waves, freak waves
Motionless wind waves in this approach appear in a natural way since the sine in the formula always reduces to zero at the direction θ=−θ K (see Eq. 21).
There appears a question whether wind waves with zero velocity exist at any azimuth value θ.To answer this question let us see at what values of the K parameter reduces to zero the coefficient K (K) in the formula (26).Applying once again to the formulas (21) we can see that there are two such values.Their expressions after some simplifications take the form: At these growth rate values the wind waves represent a dynamic barkhan-like water surface deformation.The height of the water barkhans diminishes if the values (27) are less than unity and increase up to overturning at K 0 >1 or K π >1.If K 0 =1 or K π =1, the wind born water barkhan are neutral.
The barkhan like waves are encountered in the vicinity of south-east Africa and more rearly along the Gulf Stream flow.In the case when the barkhan like wave's amplitude increases quickly (K>1), the absolutely motionless wave is called freak wave, rogue wave, rabid-dog wave and so on.The waves like this sometimes sink ships.In the south-east Africa region the rogue waves appear on the Agulhas Current when the wind is directly oriented along the water current (α=0).This Agulhas Current and the Agulhas Retroflection can give rise to immense rogue waves that can threaten supertankers.For that reason, mariners who successfully navigated the Cape of Good Hope frequently breathed a sigh of relief.A very well known is the New Year's Day wave, that struck the stationary Draupner oil platform in the North Sea on 1 January 1995.During this event, minor damage was inflicted on the platform, confirming that the reading was valid.From that time it was concidered that once thought to be only legendary, they are now known to be anatural phenomeon, not rare, but rarely encountered.On the Gulf Stream in the region of the Bermuda Islands, the rogue waves appear and gain height of 18 meters.In 1912 in the Bermuda Islands region a three masts sailing vessel was knocked out by a rogue waves of 30 m high.The next wave sank the vessel.

(−2δ)-hodograph
Let us assume in the formula (20) K=−2δ.Then obtain the velocity hodograph of the wind waves The main hodograph diameter ( 28) is oriented along the wind vector V .The radius of the velocity hodograph circle (28) equals half of the wind velocity.Since (−2δ)<1, the hodograph (28) (the wind velocity hodograph) belongs to the damped weaves set.
Despite the simplicity of the hodograph (28), it illustrates the main property of the considered waves: the wind waves propagate at many possible directions (in the considered case it is the direction (−π≤θ≤π)).The usual plane waves would propagate only in one direction (θ=0).

2-hodograph
Assume in the formula (20) K=2.Then obtain the velocity hodograph of wind waves The main hodograph diameter ( 29) is oriented along the velocity vector of the underwater current V 0 .The radius of the velocity hodograph circle (29) equals 1 2 V 0 .Since 2>1, the hodograph (29) (the underwater current velocity hodograph) belongs to the set of non damped waves.
Both hodographs ( 28) and ( 29) can be defined by means of measurements and simplest calculations of the angle β.Indeed, let us choose on the z plane an arbitrary pole O and draw two vectors from it: vector V of the air velocity and vector V 0 of the underwater current velocity.Consider the triangle OV V 0 .This triangle can always be plotted since its two sides (V and V 0 ) and the angle α between them can be measured.
Remind now that the locus of the main diameters of the wind waves velocities is a straight line.Since the vectors V and V 0 are the main diameters of some hodographs the whole line of the ends of the wind waves hodographs main diameters passes by the ends of the vectors V and V 0 that is, the line of the wind waves hodographs main diameters ends contains the triangle OV V 0 base opposite to the O apex.Now it is clear that the locus of the hodographs centers of the wind waves is the triangle's OV V 0 centerline parallel to the base V V 0 .
Therefore, one can conlude that the locuses of the centers of the wind waves velocities hodographs and the ends the same diameters are easily determined by means of wind and underwater current velocities measurements.
It is easily understood that the main diameter σ of the σ -hodograph is an oriented triangle OV V 0 altitude passed from the apex O.

The main hodographs diameters orientation sectors
The wind waves dynamics has one peculiarity.As we have just seen, the main diameter V of the hodograph (28) belongs to the region of non damped waves.The main diameter V 0 of the hodograph (29) belongs to the region of damped waves.What is the hodograph that is the boundary dividing the damped and the non damped waves?At first sight the answer is simple: the boundary is the hodograph λ1 = 1 sin(θ+θ 1 ) (see Eq. 21) of neutral waves (K=1).Indeed at K<1 the wind waves damp and at K>1 do not damp right up till destruction a cause of the wave amplitude growth.
So, if we move along the centers of the wind waves hodographs from the azimuth θ=β−π/2 to azimuth θ=β+π/2, the situation changes.Point is, that the hodograph K=±∞, that is, (by the Eq.18), the hodograph of the wind waves velocities of zero lengths (see Eq. 19) enters the wind waves dynamics.That is why the velocities hodograph (Eq.30) is located at the certain point C of the velocities hodograph centers line of the wind waves.Let us mark the center C of the hodograph (Eq.30) (K=∞) on the centers line of wind waves velocities hodographs.
If we move along the indicated line of the centers from the point C (K=∞) to the point K=1 we will pass the centers of the velocities hodographs of the non damped waves (1<K<∞).At K=1 we get to the center of the velocities hodograph of neutral waves.At the further decrease of the parameter K, at 0<K<1 we move along the hodograph centers of damped waves.At K=+0 we obtain an internal resonance of the growth rate of the waves with the turbulent resistance frequency c i =1 or in dimensional form c i =κ.This will take place on the damped waves.
Note, that at K=2>1 the hodograph (29) of the water current velocity enters the velocities hodographs set of non damped waves.The main hodograph (29) diameter is located on the ray θ=α.That is why at 0≤α≤π the velocities hodographs centers of the wind waves at K→+0 go upward along the ray3 θ=β+π/2.
Let us pass now from the point C (K=−∞) to the point K=−0.We will again pass along the hodographs centers of the wind waves velocities but this time along the velocities hodographs centers we go downward along the ray θ=β−π/2.
Thus the orientation sectors of the main diameters hodographs of the damped waves at +0<K<1 and −∞<K<−0 are divided by the orientation sector of the main hodographs diameters of the damped waves velocities (1<K<∞).In this sense the damped waves at +0<K<1 and at −∞<K<−0 differ from each other and the waves at +0<K<1 we will call weakly damped in contrast to strongly damped waves at −∞<K<−0 (see Fig. 5).Let us take an arbitrary pole O on the z-plane.Draw via the pole O a counting out vector V 1 of the azimuths θ.Than via pole O draw a skeleton θ=β, θ=β+π and θ=β±π/2 of the waves lengths diagram (Fig. 6).
Let us draw a vector σ the main diameter of the σ -hodograph from the pole O along the azimuth θ=β (0≤α≤π ).Along the azimuth θ=α from the pole O draw a vector V 0 of the velocity of the underwater current.Via the ends of the vectors V and V 0 draw a straight P Q -line of the ends of the main hodographs diameters of the wind waves velocities (Fig. 6).
Draw a line MN via the triangle OV V 0 centerline of the centers of the main hodographs diameters of the wind waves velocities.
Let us fix the parameter K , make for its value a waves lengths diagram and draw a tangent θ=−θ K , θ=−θ K +π to the velocities hodograph of the wind waves with the growth rate c i =K−1.
Draw a ray θ=θ K =(π/2)sgn(K)−θ K via the pole O up to the centers lines and the ends of he main wind waves velocities diameters intersection.
As a result make a main diameter of the velocities hodograph of the wind waves with the given growth rate.The hodograph center is in the intersection point of the line MN with the ray θ=θ K .The main hodograph diameter end of the velocities is in the intersection point of the line MN with the same ray θ=θ K (Fig. 6).
Knowing the ray θ=θ K let us define the maximum wind waves velocity K of the given growth rate and their length In the intersection point of the strait line θ=β±π/2 with the velocities find the velocity σ K of the waves of maximal length Having fixed some direction θ (Fig. 6), find the velocity OA and the wave length OB in this direction.
The non gravitational wind waves have one more interesting property: the waves of the equal velocity value can have different wave lengths.The waves of equal wave length can have different velocities.
For example, in the direction θ=β, the wave velocity equals σ , and the wave length equals zero.The waves symmetric with respect to the axis θ=θ K of the hodograph have the same velocity value OE, but not the same wave length OF .
The neutral waves (K=1) of the equal velocity and different wave length represent choppy sea.The problem of choppy sea has not been yet described analytically.Here this problem is solved with the use of a couple of circles.
The double lines on the Fig. 6 mark the Langmuir stripes that are evidently parallel to the main hodograph diameter K of the wind waves velocities at the given value of K.The lines L K B, B 0 O and others (Fig. 6) represent the troughs lines of the wind waves propagating with velocity λK along the azimuth θ.

Conclusions
A common practice in both theoretical and practical hydrodynamics is to consider sea waves as gravitational.Results presented in this work at least qualitatively coincide with both visual observations and with air-photography observations, therefore assure us that the main sea and river waves have not the gravitational but turbulent nature: the initial relations ( 14) do not contain the acceleration g.These relations have a profound hydrodynamic nature -cluster and microheterogeneous liquid structure and as the consequence their turbulent resistance -the main wind waves stimulator, since in its absence (κ=κ 1 =0) it is impossible to formulate the initial non disturbed adjacency current (2).It is quite natural, that though the approximately obtained dispersion relation ( 14), proves to be so informative.These relations not only explain all the main properties of the considered waves (variety, crestedness, propagation direction, velocity, wave length, choppy sea etc.), but also reveal new -motionless -waves that are marked on the water surface by driftweed, foam as light Langmuir stripes, and thus are observed on the water surface in windy weather as real motionless waves.The value of turbulent resiatance κ for air and water have to be measured experimentally.Recently, we obtained the value for κ in case of turbulent flow of water through a pipe (to be published).The value is κ=0.0231/s.Using this value for both water and air turbulent resistance, and reasonable estimation for wind velocity (10 m/s), water velocity (1 m/s), and α=30 • one obtains from (19) estimation of the wavelength of the Langmuir waves approximately 21 m.More precise values are subject of further experimental work.
Physically the multiple order roots of dispersion relation means more profound penetration of the influence of surface effects on the profound water layers.In this sense the turbulent part ξ 1 =ξ 2 of the dispersion relation is more important than the other parts of ( 14), including gravitational waves.The cause for this is the soundness of the description of the wind properties on the water surface resulting from these simple relations.The apparent three-dimensional wave chaos on the water surface in reality is fully regulated and in order to understand it one needs to know two angles: θ=β, θ=θ K and three circles: wave lengths diagrams and κ-hodograph of their velocities.Express information of the turbulent waves on the water surface is given only by three vectors: the wind velocity V 1 , the water current V 0 and the vector σ , direction θ=β (or θ=β+π) of which can be determined by visual observation and its magnitude is equal to the velocity of the wave front-line in this direction.That is why application of the obtained results into practice can be done on-the-fly by taking the wave characteristics (wave lengths, their velocities, propagation directions) on the monitors screens in ship deck-cabins, in dispatching offices of sea ports, drilling plants on shelves, lighthouses, etc.

Fig. 5 .
Fig. 5. Distribution of the sectors of the strongly-and weaklydamped waves