Three-Dimensional Numerical Simulation of Near-Field Hydrodynamic Response and Suspended Sediment Transport Characteristics Around the Caofeidian Port Breakwaters

Breakwater construction at meso-tidal ports fundamentally alters near-field hydrodynamics and drives harbor sedimentation, yet the three-dimensional mechanisms linking entrance geometry to sediment flux remain poorly quantified. Here, we apply a validated Delft3D tidal–sediment coupled model to Caofeidian Port, Bohai Bay, comparing pre-construction baseline conditions against four entrance width scenarios (400, 300, 250, and 200 m). Breakwater enclosure reduces depth-averaged harbor velocities by 61.9–63.2% during spring tides, while generating tip-jet velocities of 1.41–1.53 m s−1 at the eastern breakwater head—exceeding pre-construction maxima by 14–18%. The eastern tip produces an ebb vortex (radius ~230 m; peak vorticity 0.034 s−1) approximately 34% larger and 62% more intense than its flood counterpart, driving vortex-assisted sediment recirculation toward the harbor interior despite ebb-dominant background velocities. Reynolds flux decomposition confirms that the eastern tip-vortex sector contributes ~39% of net sediment import (advective component: −0.7%), directly quantifying vortex-assisted recirculation as an independent transport mechanism. Bed shear stress falls below the critical erosion threshold (τce = 0.22 Pa) across 76.8% of the harbor area during spring tides (robust lower bound ~60% under wave-coupling correction), creating a structurally stable depositional interior, while the near-entrance zone sustains persistent tidal-cycle resuspension. Asymmetric tidal pumping—flood-phase open-sea SSC of 0.088 kg m−3 versus ebb-phase harbor SSC of 0.032–0.041 kg m−3—drives net spring-tide sediment import of 14.8 × 106 kg per cycle (wave-coupled upper bound: 17.8–19.2 × 106 kg per cycle). Entrance width reduction from 400 to 300 m achieves a favorable sedimentation-to-water exchange trade-off (marginal efficiency ratio 1.23), whereas further reduction to 200 m indicates onset of hydraulic choking. The marginal efficiency ratio declines sharply from 1.23 (400 → 300 m) to 1.03 (300 → 250 m) to 1.01 (250 → 200 m), indicating a hydraulic transition within the 250–300 m range that warrants targeted refinement in future studies.