Secondary Instability Research Articles

An Extreme Predecessor Rain Event in Central China Amplified by Upper-Level Jet Streak

Abstract A comprehensive study is conducted on a predecessor rain event (PRE) closely associated with the upper-level jet (ULJ) in central China during July 2021, utilizing the Weather Research and Forecasting Model (WRF). The ULJ and precipitation patterns are reasonably reproduced in the control experiment. In the sensitivity experiment (NOULJ), the ULJ is largely suppressed by the piecewise potential vorticity inversion method. The comparison of the two simulations indicates that the ULJ promotes the development of secondary circulation, lower-level jet (LLJ), moisture transport, and atmospheric instability. Synergistic effects among these processes amplify this PRE. The strong updrafts in the precipitation region are closely associated with the inverse secondary circulation that develops on the south side of ULJ entrance. The upper-level mass divergence associated with the ULJ enhances the lower-level convergence, promoting the development of lower-level vorticity and LLJ. The LLJ facilitates moisture transport into the precipitation region and increasing atmospheric instability. In contrast, the precipitation in NOULJ is decreased by 34% due to the weakening of favorable conditions.

Influence of spanwise wall vibration on non-modal perturbations subject to freestream vortical disturbances in hypersonic boundary layers

In this paper, we study the effect of lateral wall vibrations on the excitation and evolution of non-modal perturbations in hypersonic boundary layers subject to low-frequency freestream vortical disturbances (FSVDs). A novel, high-efficiency numerical approach, combining the harmonic weakly nonlinear Navier–Stokes and nonlinear parabolised stability equation approaches, is developed, which is sufficient to accommodate both the rapid distortion of the perturbation in the leading-edge vicinity and the nonlinear development of finite-amplitude high-order harmonics in the downstream region. The boundary-layer response to low-frequency FSVDs shows a longitudinal streaky structure, for which the temperature perturbation shows much greater magnitude than the streamwise velocity perturbation. The lateral vibration induces a Stokes layer solution for the spanwise velocity perturbation, which interacts with the FSVD-induced perturbations and leads to a suppression of the non-modal perturbation and an enhancement of the downstream modal perturbation. The new perturbations excited by the FSVD–vibration interaction strengthen as the vibration intensifies, and they could become comparable with the FSVD-induced perturbations in downstream locations at a high vibration intensity, indicating a remarkable modification of the streaky structure and its instability property. Secondary instability (SI) analyses based on the streaky base flow indicate that the vibration could enhance or suppress the SI modes, depending on their initial phases over the vibration period. Overall, the average effect is that the low-frequency and high-frequency SI modes are stabilised and destabilised by the vibration, respectively. Since the high-frequency SI modes undergo higher amplifications, the subsequent bypass transition is likely to be promoted by relatively strong vibrations.

Convective Heat Transfer in Uniformly Accelerated and Decelerated Turbulent Pipe Flows

This study presents a detailed investigation of the temporal evolution of the Nusselt number (Nu) in uniformly accelerated and decelerated turbulent pipe flows under a constant heat flux using direct numerical simulations. The influence of different acceleration and deceleration rates on heat transfer is systematically studied, addressing a gap in the previous research. The simulations confirm several key experimental findings, including the presence of three distinct phases in the Nusselt number temporal response—delay, recovery, and quasi-steady phases—as well as the characteristics of thermal structures in unsteady pipe flow. In accelerated flows, the delay in the turbulence response to changes in velocity results in reduced heat transfer, with average Nu values up to 48% lower than those for steady-flow conditions at the same mean Reynolds number. Conversely, decelerated flows exhibit enhanced heat transfer, with average Nu exceeding steady values by up to 42% due to the onset of secondary instabilities that amplify turbulence. To characterize the Nu response across the full range of acceleration and deceleration rates, a new model based on a hyperbolic tangent function is proposed, which provides a more accurate description of the heat transfer response than previous models. The results suggest the potential to design unsteady periodic cycles, combining slow acceleration and rapid deceleration, to enhance heat transfer compared to steady flows.

Crossflow-Induced Breakdown and Transition Correlation for a Hypersonic Swept Plate Flow

The crossflow-induced transition for a Mach 6 flow over a swept plate with a 35 mm nose radius and a sweep angle of 45 deg is investigated using direct numerical simulations (DNSs). In order to shed light on the complete transition process, the evolution of a stationary crossflow wave is first simulated. Then, unsteady wall blowing/suction perturbations are introduced to trigger transition to turbulence. The results show that both the type I and III secondary modes are excited, and they subsequently undergo a linear stage of amplification before breakdown. The type II mode is undetected despite its amplification predicted by the two-dimensional eigenvalue stability approach. Overall, the type I mode achieves a dominant amplitude and plays a key role in the transition. Furthermore, a transition correlation method is proposed based on the dominant secondary instability and threshold amplitude concept, and in order to verify and calibrate it, wall perturbations with different amplitudes are introduced to mimic the varying intensities of the background noise. The transition locations predicted by the most amplified secondary mode agree well with those by DNSs, validating the secondary-instability-based criterion for hypersonic crossflow transition.

Influence of large-scale free-stream turbulence on bypass transition in air and organic vapour flows

The free-stream turbulence (FST) induced transition in perfect and non-ideal gas zero-pressure-gradient flat-plate boundary layers is investigated by means of large-eddy simulations. The study focuses on the influence of large incoming disturbances over the laminar-to-turbulent transition, by comparing two different integral length scales $L_f$ , which differ by a factor of seven, at different FST intensities $T_u$ . High-subsonic dense-gas boundary layers of the organic vapour Novec649, representative of organic Rankine cycle applications, are compared with air flows at Mach numbers 0.1 and 0.9. Compressibility and non-ideal gas effects are shown to be of minor importance in comparison to the influence of the FST integral length scale $L_f$ . An increase of the inlet turbulent intensity always promotes transition, whereas an increase of $L_f$ has a double effect on the transition onset. At $T_u=2.5\,\%$ , increasing $L_f$ promotes the transition, while it tends to delay transition for an FST intensity of 4 %. Larger FST integral scales tend to increase the spanwise distance between laminar streaks generated in the boundary layer. Two competing transition scenarios are observed. When the incoming turbulence intensity and length scale are moderate, the classical bypass route consists in the linear non-modal growth of streaks, which then experience secondary instabilities (sinuous or varicose) and lead to the generation of turbulent spots. The second scenario is characterized by the appearance of $\Lambda$ -shaped structures near the inlet, which are further stretched to hairpin vortices before breaking down to turbulence. Spot inceptions can therefore occur at earlier locations than the streak growth. We are then faced with a competition between the classical bypass transition and nonlinear response mechanisms that ‘bypass’ this route. The present case at high $L_f$ and low $T_u$ is an example of a competing scenario, but even for the higher $T_u$ and $L_f$ conditions, only approximately one-third of turbulent spots are due to the $\Lambda$ -shaped events. The nonlinear alternative route has strong similarities with scenarios described previously in the literature in the presence of leading edge effects or due to passing wakes. Such a path is governed by the turbulence intensity, but also by the integral length scale, with both parameters playing a critical role in the generation of the $\Lambda$ -shaped structures near the inlet. This alternative mechanism is found to be robust under varying flow and thermodynamic conditions.

Baroclinic nonlinear saturation and secondary instability of current-undercurrent meanders

Baroclinic nonlinear saturation and secondary instability of current-undercurrent meanders

Observation of an Extraordinary Type V Solar Radio Burst: Nonlinear Evolution of the Electron Two-Stream Instability

Solar type V radio bursts are associated with type III bursts. Several processes have been proposed to interpret the association, electron distribution, and emission. We present the observation of a unique type V event observed by e-CALLISTO on 7 May 2021. The type V radio emission follows a group of U bursts. Unlike the unpolarized U bursts, the type V burst is circularly polarized, leaving room for a different emission process. Its starting edge drifts to higher frequency four times slower than the descending branch of the associated U burst. The type V processes seem to be ruled by electrons of lower energy. The observations conform to a coherent scenario where a dense electron beam drives the two-stream instability (causing type III emission) and, in the nonlinear stage, becomes unstable to another instability, previously known as the electron firehose instability (EFI). The secondary instability scatters some beam electrons into velocities perpendicular to the magnetic field and produces, after particle loss, a trapped distribution prone to electron cyclotron masering (ECM). A reduction in beaming and the formation of an isotropic halo are predicted for electron beams continuing to interplanetary space, possibly observable by Parker Solar Probe and Solar Orbiter.

The Observational and Numerical Analysis of the Rayleigh–Taylor Instability beneath a Hedgerow Prominence

Using Swedish 1 m Solar Telescope Crisp Imaging Spectro-Polarimeter 6563 Å (Hα) observations and Mancha3D simulations, we analyze the formation and evolution of falling knots beneath a hedgerow prominence. By comparing the observed knot widths and kinematics to those of a parametric survey of simulations, we estimate the range of magnetic field values and characteristic wavelengths to test if the magnetic Rayleigh–Taylor instability (MRTI) can provide a physically meaningful explanation. We recover observational parameters using a novel semiautomated method and find knot velocities with a mean of −9.68 km s−1 and a mean width of 614 km. Our simulations survey a range of critical wavelengths, λ c , of 100 to 500 km, and magnetic field strengths, B 0, of 1 to 20 G, finding the closest match to observations around λ c = 300 km, and B 0 = 2 to 6 G. As both the observational and simulated values match expected values, we conclude that the MRTI can provide a physically meaningful explanation of this observation. Additionally, we also predict that the Daniel K. Inouye Solar Telescope will be able to observationally recover secondary instabilities on the leading edge of the falling mass through applying a point-spread function to an example from the simulated results.

Flow and coherent structures generated by a circular array of rigid, emerged cylinders in a shallow channel

The study investigates the flow structure, dynamics of the large-scale coherent structures, drag forces and sediment entrainment mechanisms generated by a circular array of diameter D containing rigid emerged circular cylinders of diameter d placed in a smooth-bed channel of depth h under strong shallow flow conditions (D/h = 20, d/D = 0.0125 and 0.025). Eddy resolving simulations are conducted with different values of the solid volume fraction 0.025 ≤ SVF ≤ 0.1 and non-dimensional frontal area per unit volume (2.56 ≤ aD ≤ 10.2, where for the present configuration, a = SVF(1/d)4/ ${\rm \pi}$ ) and a fixed channel Reynolds number (Reh = 10 000). The flow conditions are such that a vortex street (VS)-type of shallow wake is expected to form for a solid cylinder of diameter D. Findings are compared with previous results obtained for cases with moderately shallow conditions 2.5 ≤ D/h ≤ 3.5 and with the limiting case of flow past a solid cylinder with D/h = 20. For moderately shallow conditions, the core of the main horseshoe vortex (HV) forming around the array occupies a small fraction of the water column and its coherence is the largest in front of the array. By contrast, for very shallow conditions, multiple HVs form around the array and their cores occupy a large fraction of the water column. Moreover, the coherence of most of these HVs peaks close to the sides of the array. Only for relatively large aD values, the main HV extends over the whole upstream face of the array, similar to the limiting case of a solid cylinder. For sufficiently high aD, a secondary instability is present inside the near wake that leads to the formation of parallel horizontal vortices in the vicinity of the wake roller vortices. The changes in the wake structure with decreasing SVF and aD are qualitatively similar to those observed for cases with moderately shallow flow conditions, with the antisymmetric wake shedding mode being suppressed for aD ≤ 2.5. However, the Strouhal number associated with the shedding of wake rollers, which is still close to 0.2 for a solid cylinder with D/h = 20, can be as high as 0.4 for aD = 2.5. The paper also discusses how the steady wake and total wake lengths and the strength of the bleeding flow vary with the SVF. Simulation results show that the capacity of the flow to entrain sediment inside and around the array peaks for SVF = 0.05 and aD = 5.2, with sediment entrainment monotonically increasing with the SVF outside the array and monotonically decreasing inside the array. Despite the differences in the flow structure next to and inside the array, the variation of the mean, time-averaged streamwise drag coefficient of the solid cylinders ${\bar{C}_d}$ with aD is close to that observed for arrays with moderately shallow flow conditions. The combined drag coefficient for the array decreases with increasing flow shallowness, with the decay being stronger for relatively large values of aD.

Observations and Model of Subauroral Sporadic E Layer Irregularities Driven by Turning Shears and Dynamic Instability

AbstractObservations of coherent scatter from patchy sporadic E layers in the subauroral zone made with a 30‐MHz coherent scatter radar imager are presented. The quasiperiodic (QP) echoes are similar to what has been observed at middle latitudes but with some differences. The echoes arise from bands of scatterers aligned mainly northwest to southeast and propagating to the southwest. A notable difference from observations at middle latitudes is the appearance of secondary irregularities or braids oriented obliquely to the primary bands and propagating mainly northward along them. We present a spectral simulation of the patchy layers that describes neutral atmospheric dynamics with the incompressible Navier Stokes equations and plasma dynamics with an extended MHD model. The simulation is initialized with turning shears in the form of an Ekman spiral. Ekman‐type instability deforms the sporadic E layer through compressible and incompressible motion. The layer ultimately exhibits both the QP bands and the braids, consequences mainly of primary and secondary neutral dynamic instability. Vorticity due to dynamic instability is an important source of structuring in the sporadic E layer.

Direct Numerical Simulation of Tandem-Wing Aerodynamic Interactions at Low Reynolds Number

A three-dimensional direct numerical simulation was conducted to investigate the vortex–wing interactions through two NACA 0012 stationary wings placed in tandem at a low Reynolds number Re=10,000. The aerodynamic characteristics and three-dimensional flow structures were analyzed for the tandem wings. The back wing disturbed by the upstream vortices gained an evident increase in aerodynamic performance, where the advantage is related to the suppression of the large-scale vortex formation near the trailing edge. The Liutex method was applied to visualize the vortical structures for investigating the three-dimensional evolution and instability when interacting with the back wing. The upstream wake triggered dual-secondary vortices and intensified the secondary instability on the back wing. The induced vortices contributed to the lift enhancement because they provided an extra low-pressure region when propagating downstream on the suction side of the back wing. Because of the three-dimensional destabilization, the vortex interaction in the evolution process accelerated the transition and injection of the high-momentum flow into the boundary layer attached to the back wing, energizing the turbulent boundary layer and eliminating the large-scale separation near the trailing edge. This study provided a new perspective on the enhanced aerodynamic performance of tandem layout.

Quantum vortex stability in draining fluid flows

Quantum vortices with more than a single circulation quantum are usually unstable and decay into clusters of smaller vortices. One way to prevent the decay is to place the vortex at the center of a convergent (draining) fluid flow, which tends to force vortices together. It is found that while the primary splitting instability is suppressed in this way (and completely quenched for strong enough flows) a secondary instability can emerge in circular trapping geometries. This behavior is related to an instability of rotating black holes when superradiantly amplified waves are confined inside a reflective cavity. The end state of the secondary instability is dramatic, manifesting as a shock wave that propagates round the circular wall and nucleates many more vortices. Published by the American Physical Society 2024

The Effect of Domain Length and Initialization Noise on Direct Numerical Simulation of Shear Stratified Turbulence

Direct numerical simulation (DNS) has been employed with success in a variety of oceanographic applications, particularly for investigating the internal dynamics of Kelvin–Helmholtz (KH) billows. However, it is difficult to relate these results directly with observations of ocean turbulence due to the significant scale differences involved (ocean shear layers are typically on the order of tens to hundreds of meters in thickness, compared to DNS studies, with layers on the order of one to tens of centimeters). As efforts continue to inform our understanding of geophysical-scale turbulence by extrapolating DNS results, it is important to understand the impact of model setup and initial conditions on the resulting turbulent quantities. Given that geophysical-scale measurements, whether through microstructures or other techniques, can only provide estimates of averaged TKE quantities (e.g., TKE dissipation or buoyancy flux), it may be necessary to compare mean turbulent quantities derived from DNS (i.e., across one or more complete billow evolutions) with ocean measurements. In this study, we analyze the effect of domain length and initial velocity noise on resulting turbulent quantities. Domain length is important, as dimensions that are not integer multiples of the natural KH billow wavelength may compress or stretch the billows and impact their energetics. The addition of random noise in the initial velocity field is often used to trigger turbulence and suppress secondary instabilities; however, the impact of noise on the resulting turbulent energetics is largely unknown. In this study, we conclude that domain lengths on the order of 1.5 times the natural wavelength or less can affect the resulting turbulent energetics by a factor of two or more. We also conclude that increasing the amplitude of random initial velocity noise decreases the resulting turbulent energetics, but that different realizations of the random noise field may have an even greater impact than amplitude. These results should be considered when designing a DNS experiment.

Streamwise-Elongated Sinusoidal Roughness Elements with Enhanced Laminarizing Effect on Three-Dimensional Boundary Layer

As a laminar flow control device for delaying the crossflow-induced transition of a three-dimensional boundary layer, sinusoidal roughness elements (SREs) are placed in a Falkner–Skan–Cooke boundary layer, and the resultant laminarizing effect is numerically investigated in comparison with discrete roughness elements (DREs). Because SREs are elongated in the streamwise direction and designed to avoid flow tripping, the critical height of SREs is much higher than that of DREs. Moreover, the wake flow behind SREs efficiently generates and sustains crossflow vortices that are not dangerously unstable against secondary instabilities but able to strongly distort the mean crossflow profile into a less unstable one. By measuring this mean flow distortion by SREs and DREs, the laminarizing effect is compared among them. It is shown that the effect of SREs is higher than that of DREs and can be enhanced by choosing the appropriate height, angle, and wavelength depending on the local boundary-layer profile.

Development of a full aperture backscatter system for the Orion laser.

A full aperture backscatter system (FABS) is currently in development on the Orion laser at AWE to measure scattered light from the stimulated Brillouin scattering (SBS) and stimulated Raman scattering (SRS) processes. The light is to be collected through the full aperture of the final optic assembly and traverse back down the beam path, with fractions of this light being directed to an optics table. By measuring the energy of this backscattered light, it is possible to gain insight into some of the laser-plasma instabilities that are present on the laser system and should give an indication of some of the scattered light losses due to the SBS and SRS processes. The uncontrolled scattered light can lead to hotter electrons, which then preheat the target causing a degradation in compression and can inhibit ignition in inertial confinement fusion experiments, as well as secondary instabilities whereby the scattered light may in itself cause further LPIs, such as two-ion decay or the Langmuir decay instability. The FABS diagnostic on Orion is planned to enhance the plasma diagnostics suite available and provide quantitative guidance on increasing the energy coupling. Current progress includes the characterization of filters and, hence, a broadband xenon lamp to be used in measuring the transmission efficiency of the optics chain, desktop alignment of the backscatter optics, and characterization of the streak cameras.

Breakdown mechanisms induced by stationary crossflow vortices in hypersonic three-dimensional boundary layers

This study investigates the crossflow breakdown of a Mach 6 flow over a swept flat plate by direct numerical simulation (DNS) considering three cases with different spanwise wavenumbers of stationary vortices. Transition in these cases is initiated by the linear and nonlinear evolution of these vortices, followed by secondary instabilities and breakdown due to type-I, type-II modes, and wall blowing/suction perturbations, respectively. The results showed that amplified secondary instabilities significantly distort the mean flow, causing a steep rise in the wall friction coefficient. Fourier analysis shows that, in this fast-varying flow region, the low-frequency disturbances undergo significantly greater amplifications than high-frequency disturbances. Moreover, the stability characteristics of the time- and spanwise-averaged mean flow were examined to elucidate the breakdown mechanisms. It was found that the unstable region initially contracts to a lower frequency band and then expands significantly in the spanwise wavenumber range at low frequencies. This suggests the significant amplifications of low-frequency disturbances, consistent with the observations from DNS. These amplified low-frequency disturbances, in turn, modify the mean flow, leading to the final breakdown. The presented mechanisms, highlighting the critical role of low-frequency disturbances in the breakdown process, are likely to be universally relevant across various parameter regimes.

Concurrent Joint Contact in Anterior Cruciate Ligament Injury induces cartilage micro-injury and subchondral bone sclerosis, resulting in knee osteoarthritis.

Anterior Cruciate Ligament (ACL) injury initiates post-traumatic osteoarthritis (PTOA) via two distinct processes: initial direct contact injury of the cartilage surface during ACL injury, and secondary joint instability due to the ACL deficiency. Using the well-established Compression-induced ACL rupture method (ACL-R) and a novel Non-Compression ACL-R model, we aimed to reveal the individual effects of cartilage compression and joint instability on PTOA progression after ACL injury in mice. Twelve-week-old C57BL/6J male were randomly divided to three experimental groups: Compression ACL-R, Non-Compression ACL-R, and Intact. Following ACL injury, we performed joint laxity testing and microscopic analysis of the articular cartilage surface at 0 days, in vivo optical imaging of matrix-metalloproteinase (MMP) activity at 3 and 7 days, and histological and microCT analysis at 0, 7, 14, and 28 days. The Compression ACL-R group exhibited a significant increase of cartilage roughness immediately after injury compared with the Non-Compression group. At 7 days, the Compression group exhibited increased MMP-induced fluorescence intensity and MMP-13 positive cell ratio of chondrocytes. Moreover, histological cartilage degeneration was observable in the Compression group at the same time point. Sclerosis of tibial subchondral bone in the Compression group was more significantly developed than in the Non-Compression group at 28 days. Both Compression and Non-Compression ACL injury initiated PTOA progression due to joint instability. However, joint contact during ACL rupture also caused initial micro-damage on the cartilage surface and initiated early MMP activity, which could accelerate PTOA progression compared to ACL injury without concurrent joint contact.

Weakly nonlinear analysis of rotating magnetoconvection with anisotropic thermal diffusivity effect

The influence of anisotropic thermal diffusive coefficient on the stability of the horizontal fluid planar layer rotating about its vertical axis and permeated by the horizontal homogeneous magnetic field is studied. The linear stability analysis is carried out using the normal mode method. The stationary cross/oblique and parallel modes are calculated for different ranges of control parameters arising in the system. The SA (Stratification Anisotropy) parameter, α (the ratio of horizontal and vertical thermal diffusivities), plays a key role in deciding the boundaries between these modes and their instability regions when there is a combination of high and low rotation with weak and strong magnetic fields. The obtained isotropic results coincide with those obtained by pioneers in the literature. The weakly nonlinear behavior of the stationary convective motion in the vicinity of primary instability threshold is studied using the two-dimensional Landau–Ginzburg (LG) equation with cubic nonlinearity. This equation derived using the multiple scale analysis is similar to the one obtained in the literature having different relaxation time, nonlinear coefficient, and coherence lengths. These coefficients are used to study the heat transfer rate. In the case of high rotation, Nusselt number gets decreased from atmospheric (α < 1) to oceanic (α > 1) SA types. The domain for secondary instability of Eckhaus is obtained using the spatiotemporal LG equation and it is observed that the Eckhaus instability region decreases with increasing α.

Large eddy simulation of shock wave/boundary layer interactions in a transonic compressor cascade

In this paper, large eddy simulation (LES) was performed to investigate the shock wave/boundary layer interaction (SBLI) phenomenon in transonic compressor cascades with a chord Reynolds number of 2.12 × 106. A comprehensive analysis was conducted on both the SBLI structures inherent to the transonic compressor cascade and the coherent vortex structures within the boundary layer. The underlying mechanisms of the shock-induced boundary layer transition and the shock low-frequency unsteadiness in the transonic compressor cascade were elucidated through spectral and dynamic mode decomposition (DMD) analysis. The results revealed that boundary layer separation induced by the SBLI cannot reattach, leading to the formation of large-scale coherent vortex structures. Spectral analysis revealed that the shock-induced boundary layer transition in the transonic compressor cascade was dominated by inviscid Kelvin–Helmholtz (K–H) and secondary instability mechanisms, characterized by a dimensionless Strouhal number of 0.06. Additionally, pressure signals showed the variations in sub-frequency from the separated shear layer to the main flow. The oscillation amplitude of the shock foot was significantly greater than that of the shock main body, and the oscillation frequency of the shock foot was consistent with the sub-frequency. The oscillation frequency of the shock main body coincided with that of the compression ramp and flat plate configurations. Finally, DMD modal analysis indicated that high-frequency modes were correlated with turbulent fluctuations in the boundary layer, while medium- and low-frequency modes corresponded to shedding motion in the separated shear layer and low-frequency motion of the shock. This work promotes the understanding of the complex flow mechanisms of SBLI in the transonic compressor cascades.

Flow dynamics and mixing of jet in crossflow with cylindrical cavity

The flow dynamics and mixing characteristics of an air jet issued from a cylindrical cavity in an air crossflow are numerically studied using a large-eddy-simulation technique. The cavity, aligned concentrically with the jet, is located beneath the crossflow wall. The jet-to-crossflow velocity ratio is 4, and the Reynolds number for the jet flow is 1.39×104 based on its diameter and centerline velocity. The presence of the cavity significantly influences the early evolution of the jet and its interaction with the crossflow. Complex vortical structures are observed. Notably, windward vortices on the jet surface increase in size, accompanied by a reduction in the Strouhal number. For a deep cavity, these vortices break down and result in small vortical tubes in the jet streamwise direction due to secondary instability. Also examined are leeward shear-layer vortices, hanging vortices, wake vortices, and the recirculating flow within the cavity. Their roles in the mixing between the jet fluid and the crossflow are identified. The cavity enhances mixing. The effect is significant in the near field but diminishes in the far field. By adjusting the cavity radius and depth, it is determined that the cavity depth exercises a more profound impact on jet evolution and mixing than the cavity radius. The most substantial influence occurs when a narrow and deep cavity is implemented. These findings may serve as guidelines for optimizing cavity design for effective modulation of jet behaviors.