The formation of a recirculating subsurface core in an internal solitary wave (ISW) of depression, shoaling over realistic bathymetry, is explored through fully nonlinear and nonhydrostatic two-dimensional simulations. The computational approach is based on a high-resolution/accuracy deformed spectral multidomain penalty-method flow solver, which employs the recorded bathymetry, background current, and stratification profile in the South China Sea. The flow solver is initialized using a solution of the fully nonlinear Dubreil-Jacotin-Long equation. During shoaling, convective breaking precedes core formation as the rear steepens and the trough decelerates, allowing heavier fluid to plunge forward, forming a trapped core. This coreformation mechanism is attributed to a stretching of a near-surface background vorticity layer. Since the sign of the vorticity is opposite to that generated by the propagating wave, only subsurface recirculating cores can form. The onset of convective breaking is visualized, and the sensitivity of the core properties to changes in the initial wave, near-surface background shear, and bottom slope is quantified. The magnitude of the near-surface vorticity determines the size of the convective-breaking region, and the rapid increase of local bathymetric slope accelerates core formation. If the amplitude of the initial wave is increased, the subsequent convective-breaking region increases in size. The simulations are guided by field data and capture the development of the recirculating subsurface core. The analyzed parameter space constitutes a baseline for future three-dimensional simulations focused on characterizing the turbulent flow engulfed within the convectively unstable ISW.

The pressure field inside a porous bed induced by the passage of an Internal Solitary Wave (ISW) of depression is examined using high‐accuracy numerical simulations. The velocity and density fields are obtained by solving the Dubreil‐Jacotin‐Long Equation, for a two‐layer, continuously stratified water column. The total wave‐induced pressure across the surface of the bed is computed by vertically integrating for the hydrostatic and nonhydrostatic contributions. The bed is assumed to be a continuum composed of either sand or silt, with a small amount of trapped gas. Results show variations in pore‐water pressure penetrating deeper into more conductive materials and remaining for a prolonged period after the wave has passed. In order to quantify the potential for failure, the vertical pressure gradient is compared against the buoyant weight of the bed. The pressure gradient exceeds this weight for weakly conductive materials. Failure is further enhanced by a decrease in bed saturation, consistent with studies in surface‐wave induced failure. In deeper water, the ISW‐induced pressure is stronger, causing failure only for weakly conductive materials. The pressure associated with the free‐surface displacement that accompanies ISWs is significant, when the water depth is less than 100 m, but has little influence when it is greater than 100 m, where the hydrostatic pressure due to the pycnocline displacement is much larger. Since the pore‐pressure gradient reduces the specific weight of the bed, results show that particles are easier for the flow to suspend, suggesting that pressure contributes to the powerful resuspension events observed in the field.

No abstract

No abstract