Kelvin-Helmholtz Billows Induced by Shear Instability along the North Passage of the Yangtze River Estuary , China

Kelvin-Helmholtz (K-H) instability plays a significant role in mixing. To investigate the existence of K-H instability along the North Passage of the Yangtze River Estuary, the non-hydrostatic model NHWAVE is utilized to simulate the fresh-salt water mixing process along the North Passage of the Yangtze River Estuary. Using high horizontal resolution, the structure of K-H billows have been successfully captured within the Lower Reach of the North Passage. The K-H instability occurs between the max flood and high-water slack. The duration and length scale of the K-H billows highly depends on the local interaction between fresh-water discharge and tide. The horizontal length scale of the instability is about 60 m, similar to the observations in other estuaries. In the vertical direction, the K-H billows exist within the pycnocline with length scale ranging from 6 to 7 m. The timescale of the billows is approximate 6 min. By analyzing the changes of potential energy during the mixing process, results show that the existence of K-H instability induces intense vertical mixing, which can greatly increase mixing efficiency in the North Passage of the Yangtze River Estuary.


Introduction
An estuary is a place where a freshwater river meets the salty open-sea.The hydrodynamics in estuaries are very complex due to the presence of multi-scale physical processes, such as river discharge, tides, winds, waves and offshore currents [1][2][3].The periodic mixing and stratification induced by the interaction of river discharge and tidal current plays an important role in mass transport and water quality in estuaries.Tidal-scale processes in fresh-salt water mixing are well understood through in situ observations, laboratory experiments and numerical simulations, however, relative little is know about the hydrodynamic processes on intra-tidal time scales.These small scale flow structures, such as turbulence and instabilities, are the dominant mechanism for the transition from organized flow to turbulence, which is critical to the understanding of the hydrodynamics and mass transport in estuarine areas [4][5][6][7].
The Kelvin-Helmholtz (K-H) instability is one of the intra-tidal scale process, which is a common phenomenon of sheared flows in ocean.The theoretical studies of K-H instability originate from Thomson [8] and von Helmholtz [9].These theories describe the development of the instability and the transition to turbulent mixing in fluids with velocity shear.K-H instability can occur when the destabilizing influence of the velocity shear overcomes the stabilizing effect of the buoyancy force.
Miles [10] and Howard [11] showed that a necessary condition for K-H instability in a parallel, stratified, inviscid flow, is that the gradient Richardson number (Ri = N 2 /S 2 , where N = −(g/ρ)(∂ρ/∂z) is the Brunt-Väisälä frequency, S = du/dz is the velocity shear, g is the gravity acceleration, ρ is density and u is a representative flow speed) is less than 0.25 somewhere in the flow.However, it has been demonstrated this criterion is not sufficient [12], because it is possible to have a stable shear layer when Ri < 0.25 at the pycnocline.Fringer and Street [13] found the interfacial wave can be stable with Ri < 0.25, but unstable perturbations occurred when Ri < 0.13.Barad and Fringer [14] used an adaptive numerical method to evaluate the critical Ri for instability, and the result showed a similar value of Ri < 0.1 is required for instability.Based on laboratory experiments, Fructus et al. [15] proposed a new criterion for instability, which is L x /λ > 0.86, where L x is the length of the region with Ri < 0.25 and λ is the wave width.Another alternative criterion for instability is based on the linear stability analysis with the Taylor-Goldstein equation.Troy and Koseff [16] used the equation to derive the criterion for instability, which requires σi T w > 5, where T w is the time the fluid spends in region with Ri < 0.25 and σi is the averaged growth rate of the instabilities in the region.
In situ observations of K-H instabilities are few due to the limitation of instruments in recording the spatial distribution of density with high resolution.The first observation of K-H instability was conducted by Woods [17].Through a dye-release experiment in the stratified thermocline of the Mediterranean Sea, the K-H billows associated with internal waves were successfully observed.However, most observations of K-H billows were obtained by echo sounder.Geyer and Smith [18], Geyer and Farmer [19] investigated the existence of K-H billows at the Fraser River Estuary, Canada.The images of K-H billows were most apparent during ebb tide.During the flood tide, there is some mixing but not sufficient to cause instability.Bourgault et al. [20] undertook a field experiment in the St. Lawrence Estuary, but instabilities were detected during both of ebb and flood tide.That means the occurrence of K-H instability depends on the local hydrodynamics and could vary due to the changes of upstream discharge and tides from open-sea.Although the occurrence of K-H instability varies in estuaries, it is consistent that the K-H instability occurs in the pycnocline rather than at the bottom and it is the main mechanism for mixing in the pycnocline [21].
The length scales of K-H billows vary significantly in estuaries.The observation conducted by Geyer and Farmer [19] showed K-H billows with a wavelength of approximate 10 m.In the St. Lawrence estuary, however, the observed length scale is 140-150 m and the vertical size is 10-25 m.Tedford et al. [22] analyzed observations made in the Fraser River Estuary.The length scale of the instabilities is impacted by both of discharge and tidal current.The horizontal length of the billows ranges from 20 to 65 m.More recently, Chang et al. [23] captured a train of K-H billows with a horizontal length scale of 200 m in the Kuroshio off southeastern Taiwan.
Although there is not a consistent length scale of the billows, the observations has indicated that the height to length ratio of K-H billows is about 0.1 [24].Thus, the typical resolution for well simulating the billows are at least 10-150 m in the horizontal direction and 1-15 m in the vertical direction.Furthermore, high resolution field measurements in the Connecticut River estuary [21,25] have proved the existence of secondary instabilities with wavelength of 0.5-1 m, which play an important role in the turbulent mixing within the pycnocline.That means accurate simulation of the K-H instabilities needs high spatial resolution with grid size of O(0.1 m)-O(1 m).Such small spatial grid is very time consuming in large scale ocean model.Most of the existing ocean models, such as ROMS [26], FVCOM [27] and Delft3D [28,29], are based on the hydrostatic assumption.The vertical accelerating terms are omitted in the momentum equation, hence the vertical mixing process cannot be well represented.
To well capture K-H instability during fresh-salt water mixing in estuaries, non-hydrostatic models are needed, in which the total pressure is divided into two parts: dynamic pressure and hydrostatic pressure.The dynamic pressure can be obtained by solving the pressure Poisson equation.Compared with hydrostatic models, non-hydrostatic models can simulate dispersive waves, strongly stratified fluids and flows over abrupt topography changes, in which the effects of vertical acceleration are not neglected.The free surface in non-hydrostatic models is assumed to be a single value function of the horizontal coordinate, thus, the surface elevation can be obtained by integrating the continuity equation over depth.Thus, the models are computationally more efficient with the lack of free surface tracking.However, the application of non-hydrostatic model in large-scale simulation is still very time consuming.In coastal and estuarine areas, high resolution non-hydrostatic simulations are limited to domains with lengths of O(100 km) and most of the simulations are based on two dimensional domain.Özgökmen et al. [30] studied turbulent mixing in the Red Sea Outflow Plume using a non-hydrostatic model.The spatial grid size is 73 m.This two dimensional model revealed the interaction between shear-induced mixing and internal waves, but the horizontal resolution was too coarse to capture the evolution of K-H billows.Wang et al. [31] applied the non-hydrostatic model SUNTANS [32] to resolve small scale intra-tidal dynamics for the Snohomish River Estuary.The results showed that the non-hydrostatic effect is not necessary for simulation of large spatial-temporal scales of estuarine flow, but is very important for small-scale flow structures.Vlasenko et al. [33] simulated the evolution of a river plume using a spatial resolution of 12.5 m.The model can predict the fine internal structures of the plume.It can be seen non-hydrostatic model can well simulate the detailed plume structure in estuary with spatial resolution of O(10m).However, finer grid is needed to investigate the structures of the instabilities.
The Yangtze River, with length of 6400 km, is the longest river in Asia.Its annual mean freshwater discharge is 28,630 m 3 /s according to measurements at Datong gauge station [34].The Yangtze River estuary (YRE), which is 120 km in length and 90 km in width at the mouth, is characterized by the three-order bifurcations and four outlets (Figure 1).The mass mixing and transport are very complicated due to the interactions of river discharge, tides, waves and wind in the estuary [35][36][37][38].Among the dominant mechanisms leading to mixing, such as bottom stirring, wind stirring and K-H instability, research on K-H instabilities are rare in this area.No observed images or numerical simulations of K-H billows have been published up to now.Pu et al. [39] calculated the gradient Richardson number(Ri) along the North Passage.The results showed the values of Ri is less than 0.25 at the interface within the middle and lower reaches, suggesting the occurrence of K-H instability there.Following the previous discussion the Richardson number is not condition sufficient for the KH instability to occur, thus they did not suggest its existence but rather the possibility that such an instability might occur.To explore whether K-H billows occur along the North Passage of the Yangtze River Estuary, a non-hydrostatic model is used to predict the mixing process in this area.The following parts are organized as: in Section 2, the governing equation, numerical schemes and model setup are briefly discussed.In Section 3, the results are compared with observations.The characteristics of K-H billows are analyzed in Section 4. Finally, concluding remarks are given in Section 5.

Study Area
In 2010, the Deep Navigation Channel (DNC) was constructed within the North Passage of YRE.The northern and southern dikes are about 60 km in length and the width of the DNC is around 400 m with a mean lowest water depth of 12.5 m.The width of the DNC is much smaller that the channel length.Due to the existence of the two dikes, the along channel velocity is much larger than the cross channel one, so flow within the DNC can be approximated as vertical two-dimensional.
In the present investigation, we focus on the vertical salt-fresh water mixing in intra-tidal scales.Observations from fieldwork conducted in April 2009, including sea level, velocities and salinity, provide boundary conditions and validation data for the model.The computational domain covers the North Passage beginning from the gauge station CS0 extending to CS4 close to the North Passage exit (Figure 1).

Numerical Model
The non-hydrostatic wave model (NHWAVE) has been implemented to solve the hydrodynamic and transport equation in this study.NHWAVE was developed by Ma et al. [40], which was originally utilized to simulate dispersive waves and tsunamis.Ma et al. [41] extended the application of NHWAVE by considering the baroclinic pressure forcing in the momentum equation.To decrease the time consumed in solving the Poisson equation, Shi et al. [42] where D is the total water depth.U = (Du, Dv, Dw) T , (u, v, w) are velocities in (x, y, z) directions.F, G, H represent fluxes terms.S h , S p , S ρ , S τ are source terms, which include the bottom slope term, dynamic pressure gradient, baroclinic forcing and turbulent mixing term derived from the Boussinesq approximation, respectively.ω is the vertical velocity in σ−coordinates, which can be given by To maintain the flux and source terms are automatically balanced in the steady state, the well balanced method [43,44], which can avoid artificial flow due to bottom slope, is adopted to treat the bed slope source term.Therefore, the source term can be written as where the (x * , y * , z * ) are the Cartesian coordinate system.η represents surface elevation.p is the dynamic pressure.r is the baroclinic pressure head, which can be calculated by where ρ 0 is the constant reference density and ρ servers as the dynamic density with expression as ρ = ρ m − ρ 0 .ρ m is the density of the mixture.

Numerical Approach
In the present version of NHWAVE, regular grids are used in horizontal directions, and the σ coordinate system is defined in vertical direction to represent bottom and surface geometry accurately.Following the Keller-box scheme proposed by Stelling and Zijlema [45], the velocities (u, v, w) are placed at the cell center, but the dynamic pressure (p) is defined at the velocity-facing cell faces.This grid configuration can exactly set the free surface pressure to be zero and reduce error in calculating of dynamic pressure.Using the Keller-box scheme, non-hydrostatic models can simulate highly dispersive waves well with only 3-5 vertical layers.
The second-order Runge-Kutta scheme is adopted for time stepping to obtain second-order accuracy in NHWAVE.The momentum equation (Equation ( 2)) is discretized by a second-order Godunov type finite volume method.The fluxes at the cell faces are estimated by HLL [46] or HLLC [47] approximate Riemann solvers.The pressure-correction projection method is adopted in NHWAVE to get the velocities at each time step, which devides the total pressure into two parts: dynamic pressure and hydrostatic pressure.In the first step, the intermediate values of velocities are calculated from the momentum equation without the dynamic term (S p ).In the second step, the dynamic pressure can be obtained by solving the pressure Poisson equation, then, the velocities are corrected by the dynamic term.Therefore, the velocities at kth stage in Runge-Kutta scheme can be expresses as where Assuming strong vertical mixing, the turbulence closure proposed by Pacanowski and Philander [48] was used.The ν z is calculated from where ν 0 = 10 −3 m 2 /s is the background dissipation parameters and ν 1 = 1.5 × 10 −2 m 2 /s, α = 5, n = 2 are adjustable parameters [49].Ri is the Richardson number.

Model Setup
Bathymetry was obtained from a field survey along North Passage in February 2009.The horizontal grid spacing is set to be 2.  [34] with consideration of baroclinic forcing.
At the surface, a free-slip boundary is applied.At the bottom, the effect of bottom friction is considered with the roughness height being set to be 0.001 m.The time step (∆t) is adaptive during the simulation using the Courant-Friedrichs-Lewy (CFL) criterion , min ∆y |v i,j,k | + gD i,j , min ∆σD i,j where C is the Courant number.

Model Results
The North Passage can be divided into three parts: the Upper Reach, the Middle Reach and the Lower Reach.The Upper Reach begins at CS0 with the length of 28.6 km.The Middle Reach is between 5km above and lower station CSW.The Lower Reach is from the lower end of the Middle Reach to station CS4.As shown in Figure 1, stations CS1, CSW and CS7 are gauge stations in the Upper, Middle and Lower Reaches, respectively.The numerical results are compared with the observed water level, velocity and salinity data at the three stations.
Figure 2 compares predicted sea level at the three gauge stations with observations.The model accurately capture the oscillating sea level.Both the phase and amplitude are simulated well in the whole channel.At station CS7, the model shows an excellent agreement with the observation.At the other two stations, a slight discrepancy occurs during the ebb tide.As the simulation is based on a vertical 2D domain, it cannot provide results representing the tidal distortion and asymmetry induced by the geometry of the North Passage.Velocity and salinity measurements were made at six different relative water depths h/H = 0.0, 0.2, 0.4, 0.6, 0.8 and 1.0 (where h is defined as the vertical distance from the bottom and H is the whole water depth.0.0 means the bottom and 1.0 refers to the free surface).Figures 3-5 shows the comparisons of velocity and salinity at CS1, CSW and CS7 respectively.Overall, the magnitudes of velocity are well predicted by the model.The model can predict the variations of vertical distributions of velocity with good accuracy.However, the comparison also shows some discrepancies near the bed.It might be induced by the constant bottom roughness used in the model.In the North Passage, the bottom friction can be impacted by the bottom material and topography.
Figures 3-5 also depicts the vertical distributions of salinity at CS1, CSW and CS7.The accurate prediction of salinity is a more difficult than that of surface water level and velocity.The model overestimate the salinity near bottom st CS7 from t = 9 to 11 h.This might due to the overestimation of the inflow from the right boundary during flood tide.Evidence can be found in Figure 5 that the model slightly overestimates the velocity during this period.These results show that the non-hydrostatic model can capture the qualitative features of the hydrodynamics in the North Passage over a full tidal cycle.Further improvements, such as more reasonable bottom roughness used in numerical modelling, are needed to be carried out for better prediction.However, these improvement methods are out the scope of this paper, and the discrepancies do not greatly impact the analysis of K-H instability.The discussion of the K-H billow characteristics will be focused on in the following sections.

Discussion
YRE is a mesotidal, partially-mixed estuary.The periodic stratification and mixing, induced by the interaction of tidal currents and freshwater discharge, can greatly impact the estuarine environment.During the process from stratification to mixing, K-H instability is a critical mechanism for turbulence development [21].However, both the spatial and temporal resolutions of the observations are too coarse to observe K-H instabilities in this area.Pu et al. [39] found that the K-H instability might exist within the middle and lower reaches of the North Passage.To validate the existence of K-H instability, the numerical results are used to analyze salinity and velocity dields within the North Passage.

The Existence of K-H Instability
The value of Ri < 0.25 is necessary for the occurrence of K-H instability, but it is not sufficient.Thus, this criterion can only be used to analyze the possibility of K-H instability, but cannot detect the occurrence of K-H instability.K-H billows are the direct manifestation of the existence of K-H instability.Therefore, the distributions of the salinity are plotted and used to evaluate the existence of K-H billows.Figure 6 shows the spatial distributions of salinity at t = 7.2, 8.2 and 9.2 h, respectively.Perturbations can be noticed, which occurs in the vicinity of the pynocline.The intensity of the perturbations increases as the saltwedge moves to upstream from t = 7.2 to 8.2 h, then, decreases from 8.2 to 9.2 h.The perturbations experience their strongest intensity at t = 8.2 h.The figure also shows the perturbations is not only vigorous enough to mix the fluid around pynocline, but also cause mixing in the whole vertical domain at x = 49 − 51 km.In the lower panel of Figure 6, the intensity of the perturbations have greatly decreased, indicating the duration of the perturbations of salinity is approximate 2 h.
Figure 7 is a zoom-in of the distribution of salinity at t = 7.2 h.The figure reveals a train of K-H billows around the pynocline at 50.8-51 km downstream of station CS0.The billows are approximately 50-60 m long and 6 m high.The generation of the billows is induced by the sheared flow within the density interface.The occurrence of K-H billows is associated with strong salinity gradients and shears across the pynocline.The corresponding values of N 2 , S 2 and Ri along the dashed line in Figure 7a are calculated and shown in Figure 7b-d.The value of S 2 confirms that the shear is strong at the positions where K-H billows exist.The maximum value of N 2 are also at the density interface.The generation of the billows is induced by the sheared flow within the density interface.

The Spatial and Temporal Scales of the K-H Billows
The horizontal size of the billows is estimated by calculating the spectral representation of the resolved flow field at t = 7.2 h.The Fast Fourier Transform (FFT) is employed to transform the vertical velocity distribution into the spectral field.The expression of FFT can be given by where k x = 2π/L is the wave number, L is wave length, w n is the Fourier coefficient, N is the half number of the grid in the x direction.To get the relative importance of the various spectral components, we normalize the spectrum by the following equation: The averaged wave number ( kx ) can be expressed as At t = 7.2 h, the K-H instability mainly occurs at x = 49-51 km.Therefore, we choose the streamwise distance ranging from 49 to 51 km to calculate the spectrum.The representative depth is chosen at the center of the K-H billows, which is at z = −5 m. Figure 8 depicts the normalized spectral density of the vertical velocity.It can be seen that the maximum value of the spectrum locates at k x = 0.1 m −1 , which corresponds to the billows with horizonal length of 62.8 m.The main energy is in wave numbers k x = 0 to 0.2 m −1 .The averaged wave number is 0.105 m −1 , thus, the mean horizontal size of the billows at x = −5 m is 59.8 m.To get the averaged billow size over depth, we calculated the averaged wave number from z = −2 to −9 m.The vertical distribution of averaged wave number is shown in Figure 9.The wave number, which shows a maximum at z = −7.8m, corresponds to the horizontal size of 56-61 m.The calculated length scale is consistent with the observations at the Fraser River Estuary conducted by Tedford et al. [22].In their observation, they also found that the length scale of the billows was impacted by incoming tidal waves and river discharges.To investigate the range of the length scale during our simulation, we calculate the averaged wave number at time between from 7 h and 10 h at 5 min intervals for x ∈ [49,51] km.The wave length of K-H billows ranges between 59 m and 83 m (See Figure 10).At the beginning of the occurrence of the K-H billows from t = 7-8.2h, the structure of the billows is relatively regular with wave length of approximate 60 m.Then, the wave length increases with decreasing intensity of the K-H instability.The maximum wave length is 83 m, which occurs at t = 9.8 h.
The generation and evolution of the billows can be seen in Figure 11.The figures are plotted in 150 s interval.In the first 150 s, forced by the sheared flow, K-H billows arise from the perturbations at the salinity interface.The structures of the billows are well described by the model.Then, the billows extend the shapes both in horizonal and vertical direction, but the braids between the billows become unobvious.Finally, in Figure 11d, the structures of the billows cannot be recognized, and the stratified flow is well mixed.Therefore, the K-H billows last for about 6 min, which is consistent with the observations of Chang et al. [23].

Mixing Efficiency
In this section, the effects of K-H instability on mixing efficiency is considered.The spatial distributions of the salinity, velocities and vorticity at t = 7.2 h are plotted in Figure 12.A trains of K-H billows are clearly shown in Figure 12a.There is an obvious shear layer in the vicinity of the pynocline.At the locations of the billows, the surface horizontal velocity seems to be smaller than the surrounding values.The distribution of vertical velocity shows a remarkable difference at the locations where the billows exist.The K-H instability has greatly accelerated the vertical exchange, and the maximum magnitude of w is around 0.08 m/s. Figure 12d depicts the distribution of spanwise vorticity (ω y ).It can be seen the large values of ω y are within the braids.According to Corcos and Sherman [50], the secondary instability of spanwise vorticity within the braid plays an essential role in mixing at high Reynolds number (Re).In this case, the Re can be calculated by where u is half the velocity difference across the shear layer, δ is the half-width of the shear zone, and ν is the kinematic viscosity of the fluid.We estimate the value of Re at t = 7.2 h, which is approximate 500,000.This value is similar to the observation which shows K-H billows in Geyer et al. [21].The result of intensive values of ω y within the braid is also consistent with the observations, which demonstrates the present numerical model can capture the evolution of K-H billows at high Re.
The amount of mixing rate is estimated by calculating the changes in gravitational potential energy.The relative mixing rate (M r ) is adopted, which is given by, where P PE = (g/A b ) V ρ(x, z, t)zdV is the potential energy per horizontal area [51,52].A b and V are the horizontal area and volume for the distribution of the density field ρ(x, z, t).The M r during the predicted tidal cycle is shown in Figure 13.Again, the values at x = 50, 48 and 46 km are shown to analyze the influence of K-H instability on mixing.It is shown that M r has a notable ebb-flood variation, especially at x = 50 and 48 km.During ebb tide, the changes of the potential energy at the three selected locations seem to be same.Due to the heavy salt water moves seaward, the potential energy decreases during ebb, which leads to a negative value of M r .After that the value of M r increases from the beginning of flood.The value of M r at x = 46 km is slightly larger than 0 during flood.However, the values show perturbations at the other two points during flood tide, which correspond to the existence of K-H billows.The value of M r during the existence of K-H billows is much larger than the corresponding values at other period.The maximum value of M r at x = 50 km occurs at t = 7.2-7.5 h with maximum value of M r = 5 × 10 −6 Jm −2 s −1 .The perturbations of M r confirm that the K-H instability can improve the mixing efficiency and the mixing mainly occurs in the lower reach due to the existence of K-H instability.

Conclusions
In this study, a high resolution non-hydrostatic model has been employed to explore the existence of K-H instability in the North Passage of YRE.The spacial grid is set to be 2.5 m to capture the evolution of K-H billows.
The model can capture the existence of K-H instability during flood tide.The instability occurs 48-51 km downstream of station CS0 (See Figure 1).The duration and intensity of K-H instability decreases from x = 51 to 48 km.The horizontal length scale of the K-H billows experiences a stable period from t = 7.2 to 8.2 h, with lengths of approximately 60 m, in agreement with the field observations by Tedford et al. [22].Therefore, the wavelength increases with the decrease of the K-H instability intensity.The K-H billows have vertical scales of 6-7 m and time scales of 6 min.The changes of potential energy is adopted to estimate the mixing efficiency.The mixing rate shows perturbations during the existence of K-H instability, which indicates intensive vertical mixing induced by K-H instability.The values of M r show a significant increase during the generation of the K-H instability because the K-H instability can greatly increase the mixing efficiency.

120. 8 Figure 1 .
Figure 1.Sketch of the Yangtze River Estuary.CX, HS and JDS represent the Changxing Island, Hengsha Island and Jiuduansha Shoal.Blue circles are the locations of the gauging stations in the North Passage.The red lines represent the Northern and Southern dikes.
5 m and 40 layers in the vertical direction were used.The simulation begins on 7:00 27 April 2009 UTC and ends on 19:30 27 April 2009 UTC.The model is forced by the observed velocities at the left boundary and by sea level at the right boundary.As there are only 5 gauge stations along the North Passage which cannot provide enough information to well represent the initial state, the initial conditions of the salinity and velocity fields are extracted by the model established by Lu et al.

Figure 2 .
Figure 2. Comparisons of sea level between predicted values (Solid lines) and observations (Cycles) at CS1, CSW and CS7 within the North Passage.

Figure 10 .
Figure 10.Time series of the mean horizontal length scale of the K-H billows.

Figure 11 .
Figure 11.Generation and evolution of the billows ((a) t = 7.117 h; (b-d) are taken 2.5, 5 and 7.5 min after (a); the dashed lines indicate the locations of the same billow).

Figure 13 .
Figure 13.Time series of the mixing rate (M r ) at x = 50 km (a), x = 48 km (b) and x = 46 km (c).
added a pressure decimation interpolation (PDI) method into NHWAVE, which can greatly improve the model's efficiency without significantly impacting the model's accuracy obviously.2.2.1.Governing EquationsNHWAVE is a parallelized model based on the incompressible Navier-Stokes equations in σ coordinate system.The continuity and momentum equations are written as u * represents the intermediate value at the first step.The pressure Poisson equation in (x, y, σ) coordinate system can be derived by substituting Equation (11) into the continuity equation (Equation (1)).