Complex Vortex Research Articles

Numerical simulation of start-up and shutdown transient processes of super large flow and high-power centrifugal pump

Abstract During start-up and shutdown processes, the transient flow of centrifugal pump is complex, which significantly affects the safe and stable operation of the entire pump station. In this paper, considering the effect of the transient process of the pump station, the start-up and shutdown transient processes of a large flow and high-power centrifugal pump was numerically simulated by use of STAR CCM+. The evolution of flow field in the volute and impeller was analyzed. The results show that the flow field of the impeller is filled with complex vortex structures of varying sizes at the beginning of the start-up process of a water pump. As the flow-rate increases, the fluid flow gradually tends to stabilize. During the shutdown process of the water pump, as the flow-rate decreases, the flow field vortices increase progressively, and the flow pattern is unstable. The research results provide a basis for the operation strategy of super large flow and high-power centrifugal pump.

Unsteady Interaction Between Purge Flow and Secondary Flow in High-Lift Low Pressure Turbine

Abstract The coupling effect between the complex vortex system and the rim purge flow in the endwall region of a high-lift low pressure turbine (LPT) will significantly influence the evolution of secondary flow and the corresponding losses. This paper focused on the unsteady interaction mechanism between the rim purge flow and the secondary flow inside the high-lift LPT under the periodic wake passing. The large eddy simulation (LES) method was used to reveal the influence mechanism of rim purge flow, rotor-stator cavity interaction and unsteady wakes on the secondary flow. Detailed experimental measurement was carried out for the flow field of high-lift LPT under the influences of static purge flow. The results showed that the rim purge flow significantly increased the overturning and underturning downstream of the endwall region, resulting in aggravating secondary loss. The rotating-rim increased the difference of the circumferential velocity between the mainstream and the purge flow, aggravating the Kelvin-Helmholtz (K-H) instability, inducing the K-H vortex structure with stronger energy amplitude, and further increasing the endwall flow loss. The incoming wakes weakened the energy amplitude of the K-H vortex at the rim purge outlet, thinned the thickness of the low-energy fluid at the blade leading edge. Additionally, the interaction between the wakes and the secondary vortices further suppressed the development of the secondary flow. Nevertheless, the incoming wakes caused additional mixing losses and made a negative impact on the overall aerodynamic performance of the high-lift LPT.

Analysis of the Flow Behavior and Pressure Fluctuation of a Pump Turbine with Splitter Blades in Part-Load Pump Mode

The internal flow of a pump turbine is unstable in part-load pump mode for small guide-vane openings, and the strong vibration caused by pressure pulsation is related to the safe and stable operation of the unit. A pump turbine with a six-splitter-blade runner was chosen for unsteady simulation analyses. A standard k-epsilon turbulence model was adopted to study the unsteady flow and pressure pulsation in part-load pump mode. The predicted results show that the flow in the draft tube and the runner with splitter blades was relatively stable and the flow of the blade-to-blade channel was symmetrical. When the inlet and outlet velocity distribution of the vanes was not uniform, a vortex began to form in the stay-vane domain. The reason for this vortex formation is explained, and it is pointed out that the existence of the vortex and backflow leads to uneven velocity distribution. The unsteady calculation results showed that the pressure-pulsation peak-to-peak amplitudes in the vaneless area and guide vanes were much higher than those of other monitor points because of rotor–stator interference between the rotating runner and the vanes. In addition, the pulsation characteristics of the monitor points at different circumferential positions in the vaneless region were quite different. In the vaneless area, the velocity gradient along the circumferential direction was very large, and there was a phenomenon of backflow. Also, the pressure pulsation was 0.2 times that of the runner rotational frequency, and the blade-passing frequency was a third-order frequency. At the outlet of the guide vane, the pressure pulsation was mainly of a low frequency with a complex vortex flow. Finally, the pressure pulsation began to decrease rapidly in the stay-vane region.

Mechanism analysis of the vortex ring and its effects on an axial descending rotor

A vortex ring is a flow phenomenon involving a complex toroidal vortex system. When a rotor is in descending flight, it may enter the vortex ring state (VRS), which will affect the rotor’s aerodynamic characteristics and even endanger it. In this paper, to clarify the aerodynamic mechanism by which the vortex ring exerts its effects on an axial descending rotor, an unsteady numerical simulation method for the rotor integrated with structured moving overset grids is proposed and validated using experimental data. This numerical simulation method is then applied to analyze the aerodynamic characteristics of the rotor in descending flight. The variations of the aerodynamic forces and the flow characteristics in the slipstream are analyzed to elucidate the physical mechanisms responsible for the relationships between the aerodynamic loads, flow field, and vortices when the rotor is in the VRS. The effects of the VRS cause a sharp drop in the average aerodynamic forces, which directly affects the safety and reliability of a rotor aircraft. In the slipstream of the descending rotor, a distinct vortex ring forms and moves upward as the velocity of descent increases. The most severe VRS occurs at a nondimensional velocity of descent of 1.13 when the center of the highly unsteady vortex ring is right at the blade tip, and this can be used as an indicator of the VRS. The physical mechanism by which the VRS exerts its effects can be attributed to the low pressure induced by strong vortices. The VRS decreases the pressure on the lower surface of the blade and increases the pressure on the upper surface, resulting in a reduction in the aerodynamic loads. In comparison with the hovering state, the VRS results in a much larger vortex strength and possesses a different vortex structure.

Large-Eddy Simulation and Aeroacoustic Prediction of Supercritical Airfoil Side-Edge Noise

Airfoil self-noise is investigated on a cantilever wing with a supercritical profile with a compressible wall-resolved large-eddy simulation. The Reynolds number based on the chord is 620,000, and the angle of attack is 5°. The aerodynamic results reveal the development of a complex vortex system at the side edge, including primary, secondary, and tertiary vortices that govern aerodynamic noise production. The wall-pressure coefficient and root-mean-square pressure coefficient contours highlight the side-edge shear layer and flow impingement of the primary vortex at the pressure-side edge as important noise-generation mechanisms. Dynamic mode decomposition shows that the dominant pressure modes on the airfoil exist along primary and secondary vortex impingement lines. Far-field acoustic predictions based on both solid and porous Ffowcs-Willliams and Hawkings’ analogies demonstrate good agreement with experimental results between 1.5 and 9 kHz. Low-frequency tonal humps observed in the spectra result from duct acoustic modes excited by airfoil self-noise, a consequence of the installation effect. For this configuration, airfoil side-edge noise did not dominate over airfoil surface noise levels, suggesting that at low to moderate angles of attack, airfoil side-edge noise is not a significant contributor to the overall acoustic emissions.

Study on the safety of a magnetic strainer at the deaerator outlet under two-phase and large-scale flow fields

Iron oxide corrosion exists in the deaerator, which causes serious damage to the water supply system. To absorb iron oxide impurities for the feed water quality, a magnetic strainer device has been designed at the deaerator outlet. In order to ensure the safe and stable operation of the deaerator, the computational fluid dynamics (CFD) method was used in this work to analyze the impact of complex outlet conditions on the gas-liquid two-phase flow field, pressure field, and structural safety within a deaerator. After the strainer is installed, turbulent kinetic energy (TKE) and turbulence intensity (TI) are not monotonically changing throughout strainer region and depend on the local flow state and geometry. TI relies more on turbulence disturbance and geometry. The overall flow near the strainer is accelerated, and the dimensionless velocity in the shear layer on both sides of the strainer center bracket increases by a maximum of 153.44%. Stable vortices and complex secondary vortices are formed near the strainer. The overall eddy is more fragmented, and eddy dissipation is prone to occur in the outlet region. However, this does not cause additional large turbulence losses. The impact of the strainer on the static pressure loss of the entire deaerator is very limited, and the reduction in static pressure difference does not exceed 0.04%. Finally, suggestions are made for the structural improvement and maintenance of the strainer based on the simulation results, which provide a reference for further safe and stable application of the deaerator magnetic strainer.

Combination of Advanced Actuator Line/Disk Model and High-Order Unstructured Finite Volume Solver for Helicopter Rotors

In the research field of rotorcraft aerodynamics, there are two fundamental challenges: resolving the complex vortex structures in rotor wakes and representing the moving rotor blades in the ambient airflow. In this paper, we address the first challenge by utilizing a third-order unstructured finite volume solver, which exhibits lower numerical dissipation than its second-order counterpart. This allows for sufficient resolution of small vortex structures on relatively coarse meshes. With this flow solver, the second challenge is addressed by modeling each rotor as an actuator disk (i.e., the actuator disk model (ADM)) or modeling each blade as an actuator line (i.e., the actuator line model (ALM)). Both of the two models are equipped with an improved tip loss correction, which is introduced in detail in the methodology section. In the section of numerical experiments, the numerical convergence properties of the two types of solvers have been compared in the case of two-dimensional infinite wing. In addition, the relationship between the ALM and the lifting line theory is discussed in the cases of fixed-wing calculations. Another goal of these cases is to validate the tip loss correction presented. The validation of the ALM/ADM and comparisons of computational efficiency are also demonstrated in simulations involving both hover and forward flight rotors. It was found that the combination of the third-order finite volume solver and the ALM/ADM with the improved tip loss correction presents an efficient way of performing the aerodynamic analysis of rotor-induced downwash flow.

Numerical analysis of aero-hydrodynamic wake flows of a floating offshore wind turbine subjected to atmospheric turbulence inflows

The influence of atmospheric turbulence inflow on coupled dynamic behaviors and wakes of floating offshore wind turbine (FOWT) is pronounced, especially as blade size increases. In this study, we explore the effects of atmospheric inflow on aero-hydrodynamics and wake characteristics of a spar-type FOWT. A well-validated in-house computational fluid dynamic (CFD) solver, FOWT-UALM-SJTU, is utilized to conduct the numerical simulations. To generate the realistic atmospheric inflow condition, the large eddy simulation (LES) with long-duration is employed. Two other cases involving uniform inflow and shear inflow are conducted to provide some comparable results. Compared to uniform inflow and shear inflow, the power and thrust of FOWT under atmospheric inflow display greater instability, with their power spectrums exhibiting higher intensity in high frequency region. The variation of yaw moment in atmospheric scenario is significantly drastic, leading to a remarkable response in platform yaw motion. Notably, the dominant frequency in this scenario is the blade passage frequency, as opposed to the incident wave frequency observed in the uniform and shear scenarios. The atmospheric inflow promotes the manifestation of wake breakdown and wake meandering, resulting in the faster wake recovery of FOWT. Besides, the meandering in the far wake is more significant compared to uniform and shear scenarios. The absence of vortex rings in atmospheric scenario is observed, along with the presence of more complex and smaller vortices in far wake. Our findings highlight the remarkable impacts of atmospheric inflow on dynamic responses and wake evolution of FOWT, and it is recommended that the atmospheric inflow be taken into account when assessing wake interactions between multiple FOWTs.

A numerical and experimental investigation into the mixing mechanism of hydrogen transverse jets into an air swirl flow

As a clean fuel with the advantages of abundant reserves, high calorific value, renewability, and zero carbon emissions, hydrogen has broad application prospects in the fields of energy and power. Moreover, the mixing characteristics of hydrogen and air play a crucial role in determining combustion performance. A novel mixing method of hydrogen transverse jets into an air swirl flow was investigated via numerical and experimental approaches. The Schlieren technique and high-speed photography were employed in the experiments. The effects of various swirl numbers and jet momentum flux ratios on the flow field structure, its transient characteristics, and mixing properties were studied. The research results indicate that the complex vortex structure in the mean flow field is jointly affected by the swirl number and the jet momentum flux ratio. An increase in the jet momentum flux ratio has distinct effects on the flow unsteadiness for different swirl numbers, and there exists a critical value of the jet momentum flux ratio that substantially affects the degree of mixing and a characteristic length suitable for normalization of the axial coordinates when describing the centerline concentration decay. This study provides a reference and basis for further research on combustion in air swirl flows of hydrogen transverse jets.

Jet-induced vortices of a row of distributed engines at vertical take-off condition

Abstract The aerodynamic performance of a wing model with a row of distributed engines are investigated at the vertical take-off condition. The engines are installed near the trailing edge of the wing. During vertical take-off, the jets exit from the engines and impinge perpendicularly to the ground, providing a thrust for the aircraft. Due to the ground effects, complex vortex structures are induced by the jets. The vortices are categorised into the spanwise vortices and the chordwise vortices. The underwing vortices can lead to low-pressure regions on the lower surface of the wing, resulting in an undesirable downward force. The underwing vortex structures are affected by the ratio of the engine distance to the engine diameter ( $S/D$ ). At a small $S/D$ = 1.10, the flow field is dominated by the spanwise vortices; at a large $S/D$ = 2.78, the flow field is dominated by the chordwise vortices. The range and strength of the spanwise vortices are affected by the vortices interaction. Competition mechanism exists between the range and strength effects, which results in the non-linear variation of the wing lift coefficient with engine spacing. The details of the flow physics underneath the wing and its mechanism on the lift of the wing during take-off are investigated.

Detached-eddy simulation of the vortex system on the high-lift common research model

A complex vortex system exists on the high-lift transport aircraft with components of the nacelle, pylon, high-lift devices, strakes, and brackets, and the investigation of the vortex system is helpful in profoundly understanding the aerodynamic characteristics of actual high-lift aircraft. To investigate the vortex system on the high-lift transport aircraft and evaluate the simulation capacity of detached-eddy simulation (DES) for the vortex system on the high-lift aircraft configuration, a high-resolution DES numerical investigation is carried out on a standard high-lift transport aircraft configuration, the National Aeronautics and Space Administration (NASA) high-lift common research model (CRM-HL). Reynolds-averaged Navier–Stokes simulation (RANS) is also performed to investigate the differences between DES and RANS in predicting time-averaged aerodynamic data on the CRM-HL. After experimental validation of time-averaged surface load and flow patterns, three vortex identification methods are adopted to extract the vortex structures of the flow field. It is found that the Rortex-based vortex identification method can more distinctly extract the fine vortex structures. The DES method can reliably predict the time-averaged surface load while capturing the complex vortex system on the CRM-HL under the concerned case of this work with high resolution.

Analysis of Wind Turbine Wake Dynamics by a Gaussian-Core Vortex Lattice Technique

The development and deployment of the next generation of wind energy systems calls for simulation tools that model the entire wind farm while balancing accuracy and computational cost. A full-system wind farm simulation must consider the atmospheric inflow, the wakes and consequent response of the multiple turbines, and the implementation of the appropriate farm-collective control strategies that optimize the entire wind farm’s output. In this article, we present a novel vortex lattice model that enables the effective representation of the complex vortex wake dynamics of the turbines in a farm subject to transient inflow conditions. This work extends the capabilities of our multi-physics suite, CODEF, to include the capability to simulate the wakes and the high-fidelity aeroelastic response of multiple turbines in a wind farm. Herein, we compare the results of our GVLM technique with the LiDAR measurements obtained at Sandia National Laboratories’ SWiFT facility. The comparison shows remarkable similarities between the simulation and field measurements of the wake velocity. These similarities demonstrate our model’s capabilities in capturing the entire wake of a wind turbine at a significantly reduced computational cost as compared to other techniques.

Numerical investigation of cavitation vortex dynamics in different cavitation patterns coupled implicit large eddy simulation and boundary data immersion method

The objective of this paper is to investigate the flow characteristics of different cavitation flow patterns around a NACA (National Advisory Committee for Aeronautics) 66 hydrofoil by applying the BDIM (boundary data immersion method) and ILES (implicit large eddy simulation) with an artificial code. Meanwhile, an artificial compressibility method is also employed to consider the effects of compressibility on cavitating flow. The results present that the numerical method can effectively capture different cavitation patterns, which agrees well with the previous experimental data. Subsequently, the detailed analysis of vortex structures and dynamics for the non-cavitation (σ = 3.0), sheet cavitation (σ = 2.0), and cloud cavitation (σ = 1.6) cases with the Liutex method and the vortex enstrophy transport equation have been investigated. When cavitation occurs, the degree of turbulence and the enstrophy in the flow field have been enhanced, due to the disturbance of the velocity field. For sheet cavitation, complex vortex structures appear in the attached cavity region with high-intensity enstrophy causing by the highly intense velocity and density gradient. As the cavitation pattern transits from the sheet cavitation to the cloud cavitation, more complex vortex structures can be observed in the cavitation region. Furthermore, the value and the fluctuation amplitude of enstrophy intensity increase significantly under the effect of reentrant jet. Analysis of the enstrophy transport equation indicates that the vortex stretching term and dilatation term for cloud cavitation increase relatively significantly with the movement of the reentrant flow and are highly dependent on the cavitation evolution. In addition, the region affected by the baroclinic torque also increases.

A numerical study on the unsteady flow in the guide vane of a water-jet propulsion mixed-flow pump

Mixing and dissipation of the complex vortices in the water-jet propulsion mixed-flow pump is an important source of flow loss. Furthermore, some vortices have strong unstable feature, which adversely affect vibration and noise. The unsteady flow in the guide vane of a mixed-flow pump is studied by numerical methods. Through the spectral results of hydrodynamic excitation and the instantaneous flow field of the mix-flow pump, the basic unsteady features of the unsteady flow are obtained. Peak frequency of approximately 1.68 rotation frequency is found as the dominant frequency of the excitation in the guide vane. Dynamic mode decomposition (DMD) is then used to extract the key flow modes in the guide vane. The results illustrate the dominant flow structure in the guide vane, including upstream sweeping structures, the separate vortex structures, and large-scale passage vortex structure. The upstream sweeping wake caused the fluctuation in the boundary layer and induced flow separation. Following the separation, passage vortices are induced which has the dominant frequency 1.68 rotation frequency.

Ion-Exchange Enabled Dual-Functional Swarms with Reconfigurability and Magnetic Controllability.

In nature, many organisms are capable of self-organizing into collective groups through local communications to perform complex tasks that individuals cannot complete. To date, the reported artificial microswarms either rely on toxic chemical reactions for communication or lack the hierarchical controllability and functionality, which is unfavorable for practical applications. To this end, this exploits the ion-exchange reaction enabled hierarchical swarm composed of cationic ion exchange resin and magnetic microspheres of internal information exchange. The swarm is reconfigurable under magnetic fields, generating ordered structures of controllable mobilities and even reversed hierarchy, able to navigate in confined and complex environments. Moreover, the swarm showsinteresting communications among each other, such as merging, splitting, and member exchange, forming multi-leader groups, living crystals, and complex vortices. Furthermore, the swarm functions as a dual-functional microreactor, which can load, transport, and release drugs in a pH-enhanced manner, as well as effectively degrade antibiotics via light-enhanced Fenton-like reaction in polluted water. The organized structure of the swarm greatly improves the drug loading/transport efficiency and the local concentration of catalysts for fast pollutant removal. This design lays the foundation for the design of dual-functional micro/nanorobots for intelligent drug delivery and advanced environmental remediation.

MODELLING THE EVOLUTION OF CASTING DEFECT CLOSURE BY RADIAL SHEAR ROLLING

This paper studies the behavior of a transverse through defect in an ingot during complex vortex flow of metal deformed by the radial shift rolling (RSR) method. The aim of the work was to investigate the applicability of the RSR method for the deformation treatment of oxide-dispersed strengthened steel (ODS-steel) ingots into a final semi-finished product. Vortex flow of metal promotes dispersion of alloying solid oxides and refinement of steel structure. However, the question of development or closure of possible casting defects of ODS-steel ingot remains open. For this purpose, an ODS-steel ingot with a through-hole modelling the defect was rolled according to the assumed process scheme. Also, to visualize the evolution of the defect in more detail during the rolling process, a simulation of thesame defect development in an aluminum bar has been performed in 2 mm steps. The results of both rolls showed good correlation. Based on many cross-sections of the experimental bars, detailed 3D models of the defect evolution were built. An average reduction of 50 % in defect volume was found. Welding of the outer areas occurs immediately, but the central areas remain unclosed, stretching proportionally with the drawing of the rod during rolling. For this reason, RSR cannot be used with large cavity defects in the axial zone, but this method is guaranteed to weld surface defects even at low reductions.

Effect of Archimedes number on the dynamics of free-falling perforated disks

The dynamics of perforated disks falling freely in a large expanse of viscous fluid at rest is investigated numerically. This complex fluid–structure interaction is solved via large eddy simulation. This numerical algorithm is verified and validated with available experimental results. The influence of Archimedes number expressing the ratio between the gravity-buoyancy and viscosity effects is discussed thoroughly, including kinematics and dynamics. Two critical Archimedes numbers are identified, Arcr1≈450 and Arcr2≈950, respectively. At these two critical Archimedes numbers, both kinematic and dynamic variables change trends. In this paper, we focus on the statistics of free-falling perforated disks. With the Archimedes number Ar increasing, the average angle of attack ⟨AoA⟩ and descent velocity ⟨Uz⟩ decrease gradually, and they arrive at a fixed value finally (here, ⟨·⟩ represents a time-average result); On the contrary, the other kinetic variables change violently when Ar is around 900, for example, terminal velocity ⟨Ut⟩. Additionally, phase differences of kinematic and dynamic variables are analyzed. A constant phase difference between the nutation angle θ and normal force FN is identified, about 66°, which is independent of Ar. Vortex structures are visualized using Q-criterion, and triangular vortex is omnipresent around holes. During the descent, a helical vortex always attaches to the perforated disk outer edge. With Ar increasing, complex vortex interaction appears, for example, merging and stretching. Some unusual behaviors in the numerical results are analyzed from the perspective of wake dynamics.

Modeling the Evolution of Casting Defect Closure in Ingots through Radial Shear Rolling Processing

This paper investigates the behavior of transverse defects under significant total strain in conditions of complex vortex metal flow implemented through the radial shear rolling (RSR) method. The aim of this study is to assess the applicability of RSR processing for the in-depth transformation of small ingots of special steel into bars, particularly for the manufacturing of structural elements in specialized construction projects such as nuclear power plants. Although a substantial total strain is anticipated to enhance the steel structure and contribute to defect closure, the question of the development or closure of potential casting defects remains unclear. To address this issue, model tests were conducted to simulate the implementation of RSR processing. Defect behavior data under very complex vortex metal flow and high strain conditions were obtained for the first time and have scientific merit. A small steel ingot with a 32 mm diameter, containing a simulated artificial defect in the form of a transverse through-hole with a 5 mm diameter, was employed. During rolling, the workpiece diameter was progressively reduced by 2 mm with each subsequent pass, reaching a final diameter of 20 mm. Additionally, to provide a more detailed visualization of the defect evolution process, the same defect was modeled in an aluminum bar over six passes, and changes in defect volume and shape were analyzed after each pass. A highly detailed 3D visualization of the actual defect evolution was achieved based on cross-sections from experimental workpieces. These data corresponded to the total strain levels obtained by finite element method (FEM) simulation. Notably, a consistent similarity was observed between the test results for both metals, revealing a reduction in defect volume of up to 67.7%. The deformational welding of defects in the outer sections, encompassing one-third of the rod’s radius, occurred in the initial passes. However, defects in the axial zone of the rods remained unclosed, lengthening and gradually decreasing proportionally to the elongation of the rod, akin to conventional rolling. Consequently, the radial shear rolling (RSR) method is unsuitable for ingots with substantial discontinuities in the axial zone post-casting. Nevertheless, the method ensures the total welding of defects located in the outer zones of the ingots, even with minor applied deformations and a slight decrease in the diameter of the deformed ingot. Such data were obtained for the first time and should contribute to future investigations in this field.

The vortex dynamics characteristics in a pump-turbine: A rigid vorticity analysis while varying guide vane openings in turbine mode

In off-design conditions, pump-turbines generate vortex structures of varying scales and intensities, jeopardizing energy conversion due to unstable flow phenomena. While off-design pump-turbine internal flow dynamics have been analyzed, understanding the spatiotemporal evolution of complex vortex structures within runner blade passages, which are prone to inducing hydraulic instabilities, has to be still addressed properly. This paper applies the rigid vorticity theory, also known as Liutex, to explore vortex characteristics in the pump-turbine under different guide vane openings during turbine mode operation. The investigation focuses on three main aspects: i) spatial distribution characteristics, ii) temporal evolution characteristics, and iii) vortex morphological features along with deformation mechanisms. Results showed that the suction surface side of the blade's trailing edge commonly experiences significant shear effects, while vortices tend to emerge on both pressure and suction surface sides. Under small Guide Vane Opening (GVO) condition, vortices in blade passages show periodic evolution, and Enstrophy of the Rigid Vorticity Dilatation Term (ERDT) predominantly contributes to the helical deformation of these vortices. Conversely, under large GVO condition, vortices maintain high helicity and remain stable. This study provides insights to reduce hydraulic instability in pump-turbines and optimize energy conversion efficiency. The proposed Reynolds-averaged transport equation for enstrophy of the rigid vorticity and the relative streamline coordinate system provide advanced tools for analyzing unstable flow structures.

Study on the Formation Characteristics and Disaster Mitigation Mechanisms of Rip Currents on Arc-Shaped Beach

Rip currents are fast offshore currents generated during the breaking process of waves propagating nearshore, posing a potential life safety threat to coastal bathers. This study utilizes a Boussinesq phase-resolving model to investigate the formation mechanism of rip currents at Dadonghai Beach, based on its actual topography, and explores the characteristics of rip current formation under various wave conditions, with an emphasis on analyzing vortices, the mean water level and the spatial distribution of average velocity. The results indicate that rip current formation is significantly influenced by wave height and period. The increase in wave height and period results in more intense rip currents and higher water level fluctuations on arc-shaped beaches and on both sides of the bay, leading to complex vortex distributions. An increase in the angle of wave incidence hinders rip current formation in arc-shaped beach areas but is favorable to the generation of deflection rips on both sides of the bay. Furthermore, an increase in bottom friction inhibits rip current formation. When the water depth decreases in the channels, rip currents transition into longshore currents. The findings of this research offer valuable scientific insights into the formation mechanisms of rip currents and contribute to their prediction and prevention.