J. Ocean Eng. Technol. Search

CLOSE


J. Ocean Eng. Technol. > Volume 38(5); 2024 > Article
Shin: General Solutions to the Navier-Stokes Equations for Incompressible Flow

Abstract

Waves are mainly generated by wind via the transfer of wind energy to the water through friction. When the wind subsides, the waves transition into swells and eventually dissipate. Friction plays a crucial role in the generation and dissipation of waves. Numerous wave theories have been developed based on the assumption of inviscid flow, but these theories are inadequate in explaining the transformation of waves into swells. This study addressed these limitations by analytically deriving general solutions to the Navier–Stokes equations. By expressing the velocity field as the product of a solution to the Helmholtz equation and a time-dependent univariate function, the Navier–Stokes equations are decomposed into an ordinary differential equation and the Euler equations, which are solved using tensor calculus. This paper provides solutions for viscous flow with shear currents when applied to the water wave problem. These solutions were validated through their application to the vorticity equation. The decay modulus of water waves was compared with experimental data, showing a significant degree of concordance. In contrast to other wave theories, this study clarified the process through which waves evolve into swells.

1. Introduction

Water waves are considered one of the most fascinating natural phenomena in the marine environment. Waves in the open ocean are generally generated by winds. The transmission of energy from the wind to the water is facilitated by friction between these two fluid media. This energy is transmitted as a wave across the water surface. Consequently, friction plays a pivotal role in the circulation of fluids and is a crucial factor in wave generation. The energy transferred from the wind to waves diminishes as the wind intensity decreases, leading to the transformation of waves into swells that eventually dissipate. Therefore, friction also contributes to the dissipation of waves. Nevertheless, the mathematical description of water waves has been conducted under the assumption of irrotational, inviscid flow (Chaplin, 1980; Chappeleakimr, 1961; De, 1955; Dean, 1965; Fenton, 1988; Kishida and Sobey, 1988; Korteweg and de Vries, 1895; Rienecker and Fenton, 1981; Shin, 2016, 2019; Stokes, 1847; Stoke, 1880; Skjelbreia and Hendrickson, 1960; Vanden-Broeck and Schwartz., 1979) or under the assumption of rotational, inviscid flow (Rankine, 1863; Chen and Zou, 2019; Constantin, 2005; Dalrymple, 1974; Henry; 2008; Nwogu, 2009; Shin, 2022, 2023). Hence, these studies cannot account for the natural dissipation of waves. In a few exceptional cases, viscosity is considered in the linearized Navier-Stokes equations, e.g., water waves over a viscous mud bottom and seiche phenomena in coastal regions. The equation (Dean and Dalrymple, 1984) for long waves with bottom friction is analogous to that of a damped harmonic oscillator. The amplitude decreases exponentially with time, and the angular frequency also changes, providing an important clue for solving the Navier–Stokes (NS) equations. NS equations are a set of nonlinear partial differential equations that govern the flow of Newtonian fluids. These equations are valuable because of their ability to elucidate numerous physical phenomena. The structural design process commences by evaluating the external appearance of the design object and the loads on the surface generated by fluid flow. The NS equations are used for these determinations. Despite the wide range of practical applications, the existence of smooth solutions in three dimensions has not been proven definitively. Theoretical solutions remain limited for most cases, with only a few exceptions. Analytical approaches are deemed almost impossible for solving these problems. Thus, only numerical and experimental methods are allowed. This study applied a theoretical approach with tensor calculus to address these limitations. General solutions were introduced without considering the boundary conditions and initial conditions.
NS equations are the Euler equations with the addition of the gradient of the viscous stress tensor. The viscous stress tensor is proportional to the strain rate, making the gradient of viscous stress a linear function of the velocity field. The convective acceleration is the only nonlinear term included in the Euler equation. Shin (2022, 2023) reported the complete solutions to the Euler equations. They transformed the Euler equations into Bernoulli’s principle when the velocity field was described by the solutions of the Helmholtz equation. Using the methodology used in the studies conducted by Shin (2022, 2023) with the solution for a damped harmonic oscillator, the findings were extended to cover general three-dimensional viscous flow in the current study.
When the velocity field is described by solutions of the Helmholtz equation, the gradient of the viscous stress tensor for incompressible flow is proportional to the velocity field. The product of a solution to the Helmholtz equation and any time function is also a solution to the Helmholtz equation. By substituting the product into the NS equations, these equations can be decomposed into the Euler equations and a linear ordinary differential equation for the time function. Based on Shin’s (2022) research, this study solved the Euler equations for general flow. Some exact solutions to NS equations exist. The Taylor–Green vortex (Taylor and Green, 1937) and Antonio’s solution (Antuono, 2020) are two interesting examples. The Taylor–Green vortex was an unsteady flow of a decaying vortex in Cartesian coordinates. Antonio’s solution is considered a three-dimensional generalization of the classic two-dimensional Taylor–Green vortex. For verification, the Taylor–Green vortex and Antonio’s solutions are special cases of the solutions presented in this study. Substituting the results into the vorticity equations confirmed that the vorticity equations were satisfied.
This study adjusted Shin’s results (2022, 2023) to comply with NS equations. The process through which waves evolve into swells was elucidated. The modified wave theory was verified by comparing the decay modulus with the experimental data of Shiau and Rumer (1974) and the theoretical solution presented by Dean and Dalrymple (1984). The comparison results showed a high level of agreement.

2. Navier-Stokes Equations for Incompressible Flow

Tensor notation was used for simplicity, and the following rule was adopted when using tensor notation. Whenever the subscript appeared twice as a letter in a grouping, the expressions of this grouping were summed, with the repeated index taking on successively the values 1, 2, and 3. The non-repeated letter subscripts or numerical subscripts remain fixed during the process (Dym and Shames, 1973). The partial derivatives of a tensor are denoted using commas and indices as ()xi=(),i where xi is the position vector; the conservation of mass is presented as follows:
(1)
vi,i=0
where vi is the velocity field, which is a vector function of position, xi, and time, t. The conservation of momentum is presented as follows:
(2)
v˙i+vjvi,j=p,iρ+μvi,jj+fi
A time derivative will be denoted with the ‘over-dot’ notation. Therefore, i is the local acceleration. ρ is the density of the fluid, which is a constant; p = p(xi, t) is pressure field; μ is kinematic viscosity, which is a constant; fi is body force per unit mass, which is a conservative force. The NS equations for an incompressible fluid are composed of Eqs. (1) and (2). The vorticity is defined as follows:
(3)
Ωidef=ipqvq,p
where ipq is the Levi–Civita symbol. Taking the inner product of the vorticity and the Levi-Civita symbol results in
(4)
vq,p=vp,q+Ωiipq
Taking the inner product of the velocity vp and the deformation tensor vq,p, the convective acceleration in Eq. (2) is expressed as
(5)
vpvq,p=vpvp,q+vpΩiipq
Replacing p = j, q = i, and i = k, Eq. (5) can be presented as follows because vpΩiipq = vpΩiqip =−vpΩiqpi :
(6)
vjvi,j=vjvj,iijkvjΩk
Substituting Eq. (6) into Eq. (2), the alternative form of the momentum equation is represented as follows:
(7)
v˙i+vjvj,i+p,iρfi=ijkvjΩk+μvi,jj
Because fi is a conservative force, it has potential energy. The first term on the right-hand side represents the outer product of the velocity and vorticity. Therefore, the term is zero when the vorticity and the velocity are parallel to each other.

3. Decomposition of the Navier-Stokes Equations

A definitive clue to solving the NS equations can be found in the solution of a single-degree-of-freedom damped harmonic oscillator or long waves with bottom friction (Dean and Dalrymple, 1984). In the oscillator, the damping force is proportional to the velocity of the object because it is due to fluid viscosity: −cż, where c is the damping coefficient. Newton’s second law for a damped harmonic oscillator is presented as follows:
(8)
z¨+2ξωz˙+ω2z=0
where ω is the undamped angular frequency of the oscillator, and ζ is the damping ratio. The displacement is expressed as
(9)
z(t)=Aeξωtsin(1ξ2ωt+φ)
where A is the amplitude, ϕ is the phase, and e ξωt is the exponential decay of the amplitude. Eq. (9) indicates that the damping force reduces the natural frequency, decreases the amplitude, and eventually eliminates the oscillation. Considering that the gradient of viscous stress in Eq. (7) arises from a Newtonian fluid, similar to the oscillator, the solutions to the NS equations are presented in a similar form to Eq. (9). If the velocity field is a solution to the Helmholtz equation, i.e., vi,jj = bvi, the gradient of viscous stress is presented as μvi,jj = μbvi, where b is a constant with the unit of inverse area. Therefore, like the damped harmonic oscillator, the gradient of viscous stress is proportional to the velocity field. Eq. (9) is based on the assumption that the velocity field for viscous flow is represented as follows:
(10)
vi(t,x1,x2,x3)def=F(t)ui(t,x1,x2,x3)
where F(t) is a function of time, which is a dimensionless quantity, and ui is a velocity field for inviscid flow to satisfy the Helmholtz equation. Examples of ui include the water particle velocities presented by Shin (2022, 2023).
Theorem 1. When the velocity fields are represented as a product of a time function and solutions of the Helmholtz equation, the NS equations can be decomposed into the Euler equations and an ordinary equation for the time function.
Proof. Substituting Eq. (10) into Eq. (1), the conservation of mass is presented as follows:
(11)
ui,i=0
The gradient of viscous stress is presented as μvi,jj = μbvi = μbF (t)ui, and the local acceleration is presented as follows:
(12)
v˙i=F˙ui+Fu˙i
Substituting Eq. (12) and μvi,jj = μbF (t)ui into Eq. (7), Eq. (7) is converted to the Euler momentum equation, provided that the following condition is satisfied:
(13)
F˙ui=μbFui
Because ui ≠ 0, the solution to Eq. (13) can be presented as follows:
(14)
F(t)=eμbt
The constant, b, must be negative because the flow must eventually dissipate. Eq. (14) shows similarities to the exponential decay of the amplitude, as expressed in Eq. (9). Eq. (7) is presented as follows:
(15)
Fu˙i+F2ujuj,i+p,iρfi=F2ijkujΩ¯k
where the vorticity, Ωkdef=F(t)Ω¯k. Therefore, Ω¯k=kpquq,p.
Theorem 2. When the velocity field is a solution of the Helmholtz equation, then the vorticity is also a solution of the Helmholtz equation.
Proof. Ωi,jj =(ipqvq,p ),jj =(ipqvq,jj),p =(ipqbvq ),p = b(ipqvq,p )= bΩi
As a result, the velocity field and the vorticity field were determined using the Helmholtz equation. The next step was to determine the pressure field using Eq. (15). “ 0 < F (t) < 1”, and it has no dimensions. Hence, it is considered a non-dimensional perturbation parameter. Using the perturbation method, Eq. (15) can be solved in the following chapters for two cases: ui is steady flow and ui is unsteady flow.

4. Pressure Field for Steady Flow

When the velocity field, ui is independent of time, Eq. (10) is presented as follows:
(16)
vi(t,x1,x2,x3)=F(t)ui(x1,x2,x3)
Therefore, the local acceleration in Eq. (15) is zero. The pressure field is presented as the sum of the dynamic pressure and static pressure. The pressure field must be presented by a quadratic function with respect to F(t) because the dynamic pressure is proportional to the square of the velocity, and static pressure is proportional to the body force, as follows.
(17)
pdef=ρS(x1,x2,x3)+F2(t)ρD(x1,x2,x3)
The first and second terms mean the static and dynamic pressures, respectively. Substituting Eqs. (16)(17) into Eq. (15), the result is presented as
(18)
{S,ifi}_{D,i+ujuj,iijkujΩ¯k}F2=0
Eq. (18) is decomposed with two equations because Eq. (18) is valid for all t, and F (t) ≠ 0, using the perturbation method, as follows. The zero-order perturbation equation is presented as follows:
(19)
S,i=fi
The second-order perturbation equation is presented as
(20)
D,i+ujuj,iijkujΩ¯k=0
Integrating Eq. (19), the potential energy of the body force is determined as follows:
(21)
S=fixi
The static pressure field is determined using Eq. (21). ui is a solution of the Helmholtz equation. Hence, uj,pp = buj. Therefore, the last term of Eq. (20) is presented as follows:
(22)
ijkujΩ¯k=1bijkuj,ppΩ¯k
Using Eq. (4),
(23)
uj,pp=(up,j+Ω¯rrpj),p=up.pj+Ω¯r,prpj
Using Eq. (11), up,pj = 0, and uj,pp=Ω¯r,prpj. Substituting it into Eq. (22),
(24)
ijkujΩ¯k=1bijkrpjΩ¯r,pΩ¯k=1b(δkrδipδkpδri)Ω¯r,pΩ¯k=1b(Ω¯kΩ¯k,iΩ¯kΩ¯i,k)
Where δij is the Kronecker delta. Since divergence of curl is zero, we have Ω¯k,k=0. Therefore, Ω¯kΩ¯i,k=(Ω¯kΩ¯i),k. Let’s introduce a quantity, E,idef=(Ω¯kΩ¯i),k. Therefore, E=(Ω¯kΩ¯i),k dxi is a type of kinetic energy. Because Ω¯k is also a solution of the Helmholtz equation and is determined by the velocity field, the energy, E can be easily determined. Substituting Eq. (24) into Eq. (20), the result is presented as follows.
(25)
D,i+ujuj,i1b(Ω¯kΩ¯k,iE,i)=0
Bernoulli’s principle can be expressed as Eq. (26) by integrating Eq. (25) with respect to xi :
(26)
D+12ujujr2bΩ¯kΩ¯k=Q(t)
The dynamic pressure field is determined using Eq. (26). The second term represents the kinetic energy associated with irrotational motion. In contrast, the third term represents the kinetic energy related to the vortex. The total kinetic energy is the linear sum of the two terms because the constant, b is negative, which is consistent with the law of conservation of energy. The ratio, r, is defined as
(27)
rdef=12EΩ¯kΩ¯k
Ω¯kΩ¯k2E0 because the kinetic energy is a positive-definite function and the constant, b, is negative. Therefore 0≤r≤1.
The energy E is referred to as an anti-vortex energy because it reduces the kinetic energy due to the vortex like frictional energy. Two cases of anti-vortex energy are considered as follows:

4.1 Anti-Vortex Energy Zero Case

Ω¯1=Ω¯2 and Ω¯3,3=0 for 2-dimensional flow. Hence,
(28)
Ω¯kΩ¯i,k=0
Therefore, r = 1 for two-dimensional flow. This condition is not generally valid in three-dimensional flow. Examples include the Beltrami flows considered in the next section.

4.2 Beltrami Flow

Beltrami flows are flows in which the vorticity vector and the velocity vector are parallel to each other. The outer product of the velocity and the vorticity is zero, i.e., ijkujΩ¯k=0. Applying it to Eq. (24), Ω¯kΩ¯k,i=Ω¯kΩ¯i,k. Therefore, r = 0 for a Beltrami flow. In this case, Eq. (26) represents Bernoulli’s principle for irrotational flow, and Beltrami flow is like irrotational flow.

5. Pressure Field for Unsteady Flow

When the velocity field, ui is unsteady flow, the velocity field can be expressed as Eq. (10). A typical example is progressive water waves (Shin, 2022, 2023). Using the moving coordinate system proposed by Dean (1965), an unsteady flow can be converted into a steady flow, as presented in the previous Chapter. The moving coordinate system is defined as
(29)
Xidef=xiGi(t)
Gi (t) has a length unit. In a damped harmonic oscillator, angular frequency is modified with a damping ratio, “ ζ,” as shown in Eq. (9). Therefore, the function Gi (t) shows that viscosity also affects the angular frequency of water waves. The first term in Eq. (15) must be a quadratic term of F (t), it is assumed that the function Gi (t) is defined as
(30)
G˙j(t)def=CjF(t)
where Cj is the phase speed for water waves or a reference speed for other flow, which is a constant. The relative velocity with respect to the moving coordinate system is presented as
(31)
Ui=uiCi
Ci,j = 0 because Ci is a constant. Hence,
(32)
Xixj=δij
Using Eqs. (29) and (30),
(33)
X˙i=CiF(t)
Using the coordinate system, the velocity field in Eq. (10) can be expressed as
(34)
vi(t,x1,x2,x3)=F(t)ui(X1,X2,X3)
The pressure field can be presented by a quadratic function with respect to F (t) as follows:
(35)
p(t,x1,x2,x3)=ρS(x1,x2,x3)+F2ρD(X1,X2,X3)
where the velocity field ui is a solution of the Helmholtz equation. Using Eqs. (29), (31), and (33), the local acceleration in Eq. (15) is presented as
(36)
uit=uiXjXjt=FCjuiXj=FCjUiXj
The convective acceleration in Eq. (15) is presented as
(37)
ujujxi=(Uj+Cj)UjXi
Differentiating Eq. (35), the gradient of the pressure is presented as
(38)
pxi=F2(t)ρDXi+ρSxi
Substituting Eqs. (36)(38) into Eq. (15), the result is expressed as follows:
(39)
{Sxifi}+{DXi+UjUjXi1b(Ω¯kΩ¯kxiΩ¯kΩ¯ixk)}F2=0
The zero-order perturbation equation is presented as the same in Eq. (19), and the second-order perturbation equation is presented as
(40)
DXi+UjUjXi1b(Ω¯kΩ¯kxiΩ¯kΩ¯ixk)=0
Because dXi = dxi, integrating the above equation with respect to dxi gives Bernoulli’s principle as follows:
(41)
D+12UjUjr2bΩ¯kΩ¯k=Q(t)
The dynamic pressure field was determined using Eq. (41). Substituting Eq. (14) into Eq. (30) and integrating the result results in
(42)
Gj(t)=Cjμbeμbt+B
where B is an integral constant. For inviscid flow, Gj(t)=Cjt as proposed by Dean (1965); hence, there is a condition to determine the constant.
(43)
limμ0Gj(t)=Cjt
The integral constant is determined as
(44)
B=Cjμb

6. Results and Verification

Some exact solutions to the NS equations exist for specific conditions (Antuono, 2020; Ethier and Steinman, 1994; Taylor and Green, 1937; Wang, 1991). For verification, the Taylor–Green vortex and Antonio’s solution were proven to be specific cases of the general solutions presented in this study.

6.1 General Form of the Velocity Field

A function, f(x)= eax, where “a” is a complex number in this section, satisfies that f”(x)= a2f(x). Therefore, the general solution to the Helmholtz equation is presented with a linear combination, i.e., ui=q=1NAqeaqx1G1(t)ecqx2G2(t)edqx3G3(t), where aq, cq, and dq are complex numbers in this section, and b=aq2+cq2+dq2 to satisfy ui,jj = bui . The coefficients Aq are determined to satisfy Eq. (11), boundary conditions, and initial conditions. The vorticity is also presented with the same form. Three samples are presented in the following sections.

6.2 Vorticity Transport Equation

Differencing Eq. (2) with respect to xq and taking the inner product of the result and the Levi-Civita symbol, rqi with respect to the indices i and q, the vorticity transport equation can be derived as follows.
(45)
Ω˙r+vjΩr,j=Ωjvr,j+μΩr,jj
The left-hand side is the material derivative of the vorticity vector. The first term on the right-hand side describes the stretching or turning of vorticity caused by the flow velocity gradients. Eq. (45) can be used as a criterion to judge the validity of the solution because the velocity field has already been determined in the previous section. Substituting Eq. (10) into Eq. (45) and using Theorem 2, the result is as follows.
(46)
F(t)Ω¯˙r+F2ujΩ¯r,j=F2Ω¯jur,j
The unsteady case is only considered because unsteady flow contains steady flow. When uj is unsteady, using the moving coordinate system, Eq. (46) can be expressed as
(47)
UjΩ¯r,j=Ω¯jUr,j
The first term on the right-hand side in Eq. (47) represents vortex stretching and turning, which is zero for two-dimensional flows. Therefore, Eq. (47) is presented for two-dimensional flows as follows since Ω¯1=Ω¯2=0.
(48)
UjΩ¯3,j=0
A stream function can be defined as U1 = ψ,2 and U2 =−ψ,1. Therefore, Ω3 = ψ,ii . Ω3 = ψ,ii = because the stream function is also a solution of Helmholtz equation. Substituting it into Eq. (48), gives
(49)
U1ψ,1+U2ψ,2=U1U2+U2U1=0
Therefore, Eq. (47) is satisfied for two-dimensional flows. Ω¯j,j=0 and Uj,j = 0. Hence, Eq. (47) is expressed as
(50)
(UjΩ¯r),j=(Ω¯rUr),j
For Beltrami flow, i.e., Ω¯i=φUi, where ϕ is a scalar field, Eq. (50) is satisfied. Therefore, Eq. (47) is satisfied for Beltrami flow. As a result, the solution presented in this study satisfies the vorticity equation.

6.3 Taylor–Green Vortex

The Taylor–Green vortex solution is used to test and validate the temporal accuracy of NS algorithms (Chorin, 1968; Kim and Moin, 1985). Taylor–Green vortex solution (Taylor and Green, 1973) is presented as follows. In the domain 0≤x1, x2 ≤2π, the velocity field is expressed as
(51)
v1=sinx1cosx2F(t)
(52)
v2=cosx1sinx2F(t)
where F(t)= e− 2μt . The pressure field is expressed as
(53)
p=ρ4(cos2x1+cos2x2)F2(t)
and the vorticity is expressed as
(54)
Ω3=2sinx1sinx2F(t)
Eqs. (51)(54) were cited from Taylor and Green (1937). Therefore, v1,jj =− 2v1 and v2,jj =− 2v2. The velocity is a solution of the homogeneous Helmholtz equation. Therefore, this study for b =−2 gives the Tayor–Green vortex. Ω3,jj =−2Ω3, which is a solution of the Helmholtz equation. Because the flow is 2-D, Ω¯1=0, Ω¯2=0, and Ω¯3=Ω¯3(x1,x2). Therefore, Ω¯kΩ¯i,k=0, and r = 1. Substituting the velocity and vorticity into Eq. (26), the dynamic pressure for b =− 2 is calculated as follows:
D=14(cos2x1+cos2x2)+Q-14
For Q = 1/4, the above gives Eq. (53). Therefore, the Taylor–Green vortex is a special case for steady flow in Sec. 4.1.

6.4 Viscous Three-Dimensional Periodic Solution

Antuono’s (2020) solutions for three-dimensional problems are defined on a three-dimensional torus 0≤x1,x2,x3 L and are characterized by positive and negative helicity, respectively. The solution with positive helicity is expressed as
vi=Uo[sin(kxi-π3)cos(kxj+π3)sin(kxl+π2)-cos(kxl-π3)sin(kxi+π3)sin(kxj+π2)]e-3μk2t
If i = 1, then j = 2, and l = 3; if i = 2, then j = 3, and l = 1; if i = 3, then j = 1, and l = 2. The solution is a sample for steady flow in Sec. 4.2. The velocity field is presented with a linear combination of three trigonometric functions with respect to three variables, x1, x2 and x3 . The velocity field is a solution of the Helmholtz equation for b =− 3k2 because the linear combination is a solution of the Helmholtz equation. The vorticity field is parallel to the velocity because both solutions belong to the class of the Beltrami flow presented in Sec. 4.2. For the case with positive helicity, the solution is Ωi = ϕui where φ=-3k, and the dynamic pressure is given by Eq. (26) for r = 0. Antuono’s (2020) solutions are also special cases of the solutions presented in Sec. 4.2.

6.5 Two-Dimensional Progressive Water Waves

Using the conventional coordinate system considered in Shin (2022, 2023), the phase is defined as β = kxωt where k is the wave number. Therefore, kG corresponds to ωt = k Ct for 2-D inviscid flow. Using Eqs. (42) and (44),
(55)
kG(t)=ωμb(eμbt1)
Using power series expansion, kG = ωt + O((μbt)2), where O stands for Big-O notation. Therefore, for inviscid flow, Eq. (55) gives ωt as defined in Shin (2022, 2023). “ μb” has a unit of angular frequency. Thus, a dimensionless quantity called the “damping ratio” is defined as ξdef=μb/ω. Using the damping ratio, the phase for viscous flow is represented as
(56)
β=kxωt{1eξωtξωt}
Multiplying e ξωt by the velocity field in Shin (2022, 2023) gave a solution to the NS equation. Multiplying e− 2ξωt to the dynamic pressure in Shin (2022, 2023) gives the pressure field for viscous flow. In the conventional coordinate system, the kinematic free surface boundary condition (KFSBC) is presented as
(57)
v2=ηt+v1ηx1
where v1 and v2 are horizontal and vertical velocities, respectively, for viscous flow, and η is the wave profile. The KFSBC stipulates that there must be no flow crossing any fluid interface. The free surface, which is a fluid interface, remains unaffected by the viscosity. Consequently, the kinematic boundary condition should remain unaltered by viscosity. Using Eq. (10), the KFSBC is presented as
(58)
eξωtu2=ηt+eξωtu2ηx1
where u1 and u2 are the horizontal and vertical velocities, respectively, for inviscid flow. The dimensionless elevation of the wave profile is defined as γ = k(η +h) in Shin (2022, 2023), where h is the water depth. Using the phase defined in Eq. (56),
(59)
ηt=ηββt=βtddβ(γkh)=ωkeξωtdγdβ
and
(60)
ηx1=ηββx1=βx1ddβ(γkh)=dγdβ
Substituting Eqs. (59) and (60) into Eq. (58), results in Shin (2022, 2023) because e ξωt ≠0. The KFSBC in Shin (2022, 2023) is still valid for viscous flow. Applying the dynamic boundary condition on the free surface to Eq. (35), the wave height for viscous flow is expressed as
(61)
H=Hoe2ξωt
where Ho is the wave height at t = 0. Waves emanating from a remote source are referred to as swells, whereas waves produced in the nearby vicinity are termed seas. Specifically, waves formed due to the wind-induced forced oscillation of water are generally identified as seas. In contrast, those originating from the natural oscillation of water post the cessation of wind are labeled as swells. Eqs. (56) and (61) indicate that the height of swell waves decreases when wind energy no longer influences wave formation, leading to a longer wave period and an increase in wavelength over time. The quantity “b” represents the strength of vorticity, σ defined in Shin (2022). b =− ()2 because the length is normalized by the wavelength in Shin (2022, 2023). Using ξω =− μb, the damping ratio is determined as follows:
(62)
ξ=μ(kσ)2ω
A logarithmic decrease is used to determine the damping ratio. Using Eq. (61), the damping ratio can be expressed as
(63)
ξ=14πln(H1H2)
where H1, and H2 represent the wave height of any two successive peaks. Although the strength of vorticity, σ, is determined using the average water particle velocities on the free surface and the seabed in Shin (2022, 2023), it can also be determined using Eqs. (62) and (63), which indicate that the currents considered in Shin (2022, 2023) are generated by the same source as the waves, but they differ from the Stokes drift. This study has proven that viscosity is the source of circulation and current in fluids.
Dean and Dalrymple (1984) provide the decay modulus for long waves with frictional damping as follows:
(64)
α=πP1/4
where P is the Proudman number, Pdef=μ2/gk2h5. Shiau and Rumer (1974) carried out a series of experiments to examine the decay of seiches in a basin. The experiments were conducted in very shallow water (0.15 < h < 8.5 cm). The theoretical values of α in Eq. (64) were comparable to the experimental data in Fig. 1.
The velocity fields for viscid flows, vi, were checked. The velocity fields for inviscid flows, ui, were compared with experimental data by Le Méhauté et al. (1968) in Shin (2022, 2023). Therefore, according to Eq. (10), the function, F(t) was checked instead of the velocity field, vi, as follows. Eq. (61) is valid for long waves because Shin (2022, 2023) is valid even for ultra-shallow water waves. α = 2ξωT = 4πξ from Eq. (61) because the decay modulus is defined by the relative reduction in amplitude over one wave period. Using Eq. (62), the decay modulus is expressed as
(65)
α=4π(kh)2σ2P
Note that ω2 = gk2h for shallow water waves and μ=Pgk2h5. Using Eq. (11) in Shin (2023), Eq. (65) is expressed as
(66)
α=4π{cos1(1u¯s1u¯b)}2P
where us is the average speed in the horizontal direction over a period on still water and ub is the average speed in the horizontal direction over a period on the seabed. s = us/C and b = ub/C, where C is the celerity. Upon comparing Eq. (64) and Eq. (66), two distinct differences emerge between these equations. Eq. (64) is based on long-wave theory. Owing to the relatively small water depth in relation to the wavelength, variations in velocity in the vertical direction were disregarded, and only the average velocity was considered in Eq. (64). Consequently, the decay modulus in Eq. (64) is not dependent on the velocity distribution. On the other hand, the decay modulus in Eq. (66) is closely associated with the velocity distribution. Dean and Dalrymple (1965) used a quadratic friction law where the shear stress is directly proportional to the square of the velocity. The decay modulus in Eq. (64) was calculated assuming that bottom shear stress is presented by τ12=ρωμu1, where u1 is the horizontal velocity. Therefore, the shear stress is directly proportional to μ. The present study used a Newtonian fluid model in which shear stress is linearly related to the strain rate, expressed as τ12 = ρμ(u1,2 + u2,1), where x1 is the horizontal direction and x2 is the vertical direction. Therefore, the shear stress is directly proportional to μ . Consequently, this leads to a change in the exponent of P from 1/4 to 1/2. When the same constitutive equation and the average concept considered in the long wave theory are used, the average shear stress gradient is as follows.
(67)
τ12,2=ρωμ{cos1(1u¯s1u¯b)}u1h
Substituting Eq. (67) into the right side of Eq. (13), the decay modulus is presented as
(68)
α=4π{cos1(1u¯s1u¯b)}P1/4
For the velocity ratio, 1-U¯s1-U¯b=cos(14)=0.97, Eqs. (64) and (68) are equivalent to one another, as shown in Fig.1. Eqs. (64) and (68) are compared in Fig. 1. The comparison results indicated a significant level of concordance between Eq. (64) and Eq. (68). In contrast, this indicates a distinction between Eq. (64) and Eq. (66). The distinction elucidates the rationale behind Dean and Dalrymple’s (1965) choice to use a quadratic friction law instead of a Newtonian fluid model. Eqs. (67) and (68) express the applicability of this study even for non-Newtonian fluids, provided that the linearization of shear stresses is achieved using the technique described by Dean and Dalrymple (1984).

7. Conclusions

NS equations play a crucial role in elucidating various scientific and engineering phenomena. They are used in modeling ocean currents, water waves, and other water flows. They help with the design of the structure. Despite the wide range of practical applications, the equations are mainly solved using numerical or experimental methods, while there are some exact solutions to specific problems.
General solutions to NS equations for incompressible flow are presented in this study. The solutions were calculated by combining the solutions of the Helmholtz equation and the solution of a single-degree-of-freedom damped harmonic oscillator. They are analytical exact solutions to NS equations.
The velocity fields and the vorticity fields were presented with solutions of the Helmholtz equations. When the velocity fields are expressed as the product of the solution of the Helmholtz equation and a univariate function of time, the NS momentum equation is decomposed into the Euler momentum equations and an ordinary differential equation for the univariate function of time. The Euler momentum equations were converted to Bernoulli’s principle to determine the pressure field. The time function was determined by solving the ordinary differential equation. Consequently, the general solutions were derived through the linear combinations of the solutions to the Helmholtz equation. The solutions presented in this study are functions of class C because solutions of the Helmholtz equations are smooth (Lin, 1988).
The Taylor–Green vortex solution and the solutions proposed by Antuono (2020) are particular cases of these solutions, and these solutions satisfy the vorticity equations. The solutions presented by Shin (2022, 2023) were enhanced to integrate viscid water waves with shear currents.
As time elapses or as a swell moves away from its generating area, its height decreases while the period and wavelength increase. Eventually, the swell dissipates entirely because of the removal of the energy source. This phenomenon was substantiated mathematically in the present study. The decay modulus was compared with the experimental data from Shiau and Rumer (1974) and the solution proposed by Dean and Dalrymple (1984). The agreement was excellent. Unlike Dean and Dalrymple (1984), the study showed that the decay modulus is influenced by the velocity gradient in the depth direction.
The technique presented in this study can also be applied to non-Newtonian fluid, provided that the linearization of shear stresses is achieved through the technique described by Dean and Dalrymple (1984).

Conflict of Interest

The author declares that they have no conflicts of interest.

Fig. 1.
Decay modulus versus Proudman number [Experiment data from Shiau and Rumer (1974)].
ksoe-2024-051f1.jpg

References

BAntuono, M. (2020). Tri-periodic fully three-dimensional analytic solutions for the Navier–Stokes equations. Journal of Fluid Mechanics, 890, A23. https://doi.org/10.1017/jfm.2020.126
crossref
Chaplin, J. R. (1980). Developments of stream-function wave theory. Coastal Engineering, 3, 179-205. https://doi.org/10.1016/0378-3839(79)90020-6
crossref
Chappelear, J. E. (1961). Direct numerical calculation of wave properties. Journal of Geophysical Research, 66(2), 501-508. https://doi.org/10.1029/JZ066i002p00501
crossref
Chen, H., & Zou, Q. (2019). Effects of following and opposing vertical current shear on nonlinear wave interactions. Applied Ocean Research, 89, 23-35. https://doi.org/10.1016/j.apor.2019.04.001
crossref
Chorin, A. J. (1968). Numerical solution of the Navier–Stokes equations. Mathematics of Computation, 22, 745-762.
crossref
Dym, C. L., & Shames, I. H. (1973). Solid Mechanics: A variational approach. McGRAW-Hill.

Constantin, A. (2005). A Hamiltonian formulation for free surface water waves with non-vanishing vorticity. Journal of Nonlinear Mathematical Physics, 12(1), 202-211. https://doi.org/10.2991/jnmp.2005.12.s1.17
crossref
Dalrymple, R. A. (1974). A finite amplitude wave on a linear shear current. Journal of Geophysical Research, 79(30), 4498-4504. https://doi.org/10.1029/JC079i030p04498
crossref
De, S. C. (1955). Contributions to the theory of Stokes waves. Proceedings of the Cambridge Philosophical Society, 51(4), 713-736.
crossref
Dean, R. G. (1965). Stream function representation of nonlinear ocean waves. Journal of Geophysical Research, 70(18), 4561-4572. https://doi.org/10.1029/JZ070i018p04561
crossref
Dean, R. G., & Dalrymple, R. A. (1984). Water wave mechanics for engineers and scientists. Prentice-Hall, Inc.

Ethier, C. R., & Steinman, D. A. (1994). Exact fully 3D Navier–Stokes solutions for benchmarking. International Journal for Numerical Methods in Fluids. 19(5), 369-375. https://doi.org/10.1002/fld.1650190502
crossref
Henry, D. (2008). On Gerstner’s water wave. Journal of Nonlinear Mathematical Physics, 15(2), 87-95. https://doi.org/10.2991/jnmp.2008.15.s2.7
crossref
Fenton, J. D. (1988). The numerical solution of steady water wave problems. Computers and Geosciences, 14(3), 357-368. https://doi.org/10.1016/0098-3004(88)90066-0
crossref
Kim, J., & Moin, P. (1985). Application of a fractional-step method to incompressible Navier–Stokes equations. Journal of computational physics, 59(2), 308-323. https://doi.org/10.1016/0021-9991(85)90148-2
crossref
Kishida, N., & Sobey, R. J. (1988). Stokes theory for waves on a linear shear current. Journal of engineering mechanics, 114(8), 1317-1334. https://doi.org/10.1061/(ASCE)0733-9399(1988)114:8(1317
crossref
Korteweg, D. J., & de Vries, G. (1895). XLI. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 39(240), 422-443. https://doi.org/10.1080/1478644950862073
crossref
Le Méhauté, B., Divoky, D., & Lin, A. (1968). Shallow Water Waves: A Comparison of Theories and Experiments. Proceedings of 11th Conference of Coastal Engineeing, London, UK, 86-107. https://doi.org/10.1061/9780872620131.00
crossref
Lin, T.-C. (1988). Smoothness results of single and double layer solutions of the Helmholtz equations. Journal of Integral Equations and Applications, 1(1), 83-121. https://www.jstor.org/stable/26162915
crossref
Nwogu, O. (2009). Interaction of finite-amplitude waves with vertically sheared current fields. Journal of fluid mechanics, 627, 179-213. https://doi.org/10.1017/S0022112009005850
crossref
Rankine, WJ. M. (1863). VI. On the exact form of waves near the surface of deep water. Philosophical Transactions of Royal of London, 153, 127-138. https://doi.org/10.1098/rstl.1863.0006
crossref
Rienecker, M. M., & Fenton, J. D. (1981). A Fourier approximation method for steady water waves. Journal of Fluid Mechanics, 104(1), 119-137.
crossref
Shiau, J. C., & Rummer, R. R Jr. (1974). Decay of mass oscillations in rectangular basins. Journal of the Hydraulics Division, 100(1), 119-136. https://doi.org/10.1061/JYCEAJ.000384
crossref
Shin, J. (2016). Analytical approximation in deep water waves. Journal of Advanced Research in Ocean Engineering, 2(1), 1-11. https://doi.org/10.5574/JAROE.2016.2.1.001
crossref
Shin, J. (2019). A regression analysis result for water waves on irrotational flow over a horizontal bed. International Journal of Offshore and Polar Engineering, 29(4), 461-466. https://doi.org/10.17736/ijope.2019.hc17
crossref
Shin, J. (2022). Solution for water waves with a shear current and vorticity over a flat bed. International Journal of Offshore and Polar Engineering, 32(4), 418-423. https://doi.org/10.17736/ijope.2022.hc29
crossref
Shin, J. (2023). Numerical method for calculating fourier coefficients and properties of water waves with shear current and vorticity in finite depth. Journal of Ocean Eng Technol, 37(6), 256-265. https://doi.org/10.26748/KSOE.2023.034
crossref
Stokes, G. G. (1847). On the theory of oscillatory waves. Transactions of the Cambridge Philosophical Society, 8, 441-473.
crossref
Stokes, G. G. (1880). In: Stokes G. G, ed. (2009). Mathematical and Physical Papers vol.1: Supplement to a paper on the theory of oscillatory waves. p 314-326 Cambridge University Press: https://doi.org/10.1017/CBO9780511702242.016

Skjelbreia, L., & Hendrickson, J. (1960). Fifth order gravity wave theory. Coastal Engineering Proceedings, 1(7), 10. https://doi.org/10.9753/icce.v7.10
crossref
Taylor, G. I., & Green, A. E. (1937). Mechanism of the production of small eddies from large ones. Proceedings of the Royal Society A -Mathematical and Physical Sciences, 158(895), 499-521. https://doi.org/10.1098/rspa.1937.0036
crossref
Vanden-Broeck, J. M., & Schwarts, L. W. (1979). Numerical computation of steep gravity waves in shallow water. Physics of Fluids, 22(10), 1868-1871. https://doi.org/10.1063/1.862492
crossref
Wang, C. Y. (1991). Exact solutions of the steady-state Navier–Stokes equations. Annual Review of Fluid Mechanics, 23, 159-177. https://doi.org/10.1146/annurev.fl.23.010191.001111
crossref
TOOLS
Share :
Facebook Twitter Linked In Google+ Line it
METRICS Graph View
  • 0 Crossref
  •   Scopus
  • 174 View
  • 25 Download


ABOUT
BROWSE ARTICLES
ARTICLE CATEGORY

Browse all articles >

PUBLICATION ETHICS
FOR CONTRIBUTORS
Editorial Office
President Office BD Rm. 1302, 13 Jungang-daero 180beon-gil, Dong-gu, Busan 48821, Republic of Korea
Tel: +82-51-759-0656    Fax: +82-51-759-0656    E-mail: ksoehj@ksoe.or.kr                

Copyright © 2024 by The Korean Society of Ocean Engineers.

Developed in M2PI

Close layer
prev next