Vortex Ring Research Articles

Vortex ring beams in nonlinear P T-symmetric systems.

In this paper, we investigate a two-dimensional photonic array featuring a circular shape and an alternating gain and loss pattern. Our analysis revolves around determining the presence and resilience of optical ring modes with varying vorticity values. This investigation is conducted with respect to both the array's length and the strength of the non-Hermitian parameter. For larger values of the array's length, we observe a reduction in the stability domain as the non-Hermitian parameter increases. Interestingly, upon increasing the vorticity of the optical modes, full stability windows emerge for shorter lattice size regime.

Numerical analysis of turbulent nitrogen-diluted hydrogen flames stabilised by star-shaped bluff bodies

In this paper, we investigate the flow and flame dynamics within a combustion chamber supplied with nitrogen-diluted hydrogen fuel and air as oxidiser. The chamber is equipped with a bluff body that serves as a flame stabiliser. We focus primarily on the influence of the bluff body shape on the formation of large and small-scale vortical structures generated at the bluff body edge. We analyse how the geometry of the bluff body alters these structures and examine the impact of resulting modifications on flame characteristics, including key features such as flame length and radial extent, temperature and species field, and completeness of the combustion process. We consider four bluff body geometries, namely a typical cylindrical shape, a star-shaped wall, and a small triangular sharp and smooth wavy structure. The research employs the Large Eddy Simulation (LES) method, providing deep insight into the complex physics of unsteady flow. In the cylindrical bluff body configuration, the key mechanism driving the fuel-oxidiser mixing process is the periodic occurrence of large toroidal vortices shed at an exit of the oxidiser duct and on the bluff body edge. They exert the flow towards and away from a recirculation zone formed above the bluff body. When its geometry is irregular, the appearance of the large vortices is suppressed, however, the wall irregularities induce small-scale flow phenomena that enhance mixing in the shear layer formed between the recirculation zone and the oxidiser stream. The mixing efficiency is largely dependent on the shape of the bluff body wall, and it turns out that seemingly small changes (sharp vs. smooth small triangular parts), which in non-reactive regimes have only a minor effect on the flow dynamics, significantly change the flame length and fuel consumption rate. Compared to the cylindrical bluff body configuration, the small triangular parts cause an axial shift of the maximum temperature downstream and slow down fuel consumption. Contrary, large triangular elements in the star-shape bluff body clearly influence the flow both in non-reactive and reactive regimes. They cause the formation of local oxidiser streams directed towards the fuel-rich mixture. Consequently, the injected fuel burns at a rate similar to that observed in the case of the cylindrical bluff body, but the combustion process occurs at an average temperature approximately 200 K lower and its fluctuations are reduced. These factors may influence the thermal reduction of NOx. On the other hand, a negative aspect is the highest level of hydrogen in the closest vicinity of the bluff body, which in practical applications may accelerate the embrittlement of its surface.

Molding the acoustic cavity–analyzing the influence of toroidal vortex development on acoustic multi-bubble macrostructures under different ultrasonic horn tip diameters

The acoustic cavity structure typically experiences a sequence of transfigurations during its sinusoidal growth–collapse cycle. However, upon examining the cavity structure in aqueous bodies, it appears that the growth structure attained falls between two geometrical structures, namely, mushroom-like structure (MBS) and cone-like bubble structure (CBS), based on the actuated ultrasonic horn tip diameter. With the recurring observations of the emergence of proximal toroidal vortices, the present investigation conducts a numerical analysis exploring the vortex development under 3, 6, 13, 16, and 19 mm horn tips to establish a potential correlation between the vortex and the cavity structure. The study presents a computational fluid dynamic investigation to capture the nature of the vortex evolution, in terms of size and position, and its respective cavitation development. The first indicator of potential correlation was the equivalency of the vortex expansion–contraction frequency and the cavity's sub-harmonic frequency. It has been found that the cavity structure is molded into MBS by the presence of a symmetric locomotive vortex structure that extends up to 1.5 times the horn tip diameter. Meanwhile, CBS is observed to take shape in the presence of an eccentric locomotive vortex that attains a size within 0.2–0.6 times the horn tip diameter. The significance of the vortex size and position is also observed in the cavity's collapse, as the vortex appears to govern the ability of the cavity impinging jet to initialize the collapse phase.

Paraxial and ray approximations of acoustic vortex beams.

A compact analytical solution obtained in the paraxial approximation is used to investigate focused and unfocused vortex beams radiated by a source with a Gaussian amplitude distribution. Comparisons with solutions of the Helmholtz equation are conducted to determine bounds on the parameter space in which the paraxial approximation is accurate. A linear relation is obtained for the dependence of the vortex ring radius on the topological charge, characterized by its orbital number, in the far field of an unfocused beam and in the focal plane of a focused beam. For a focused beam, it is shown that as the orbital number increases, the vortex ring not only increases in radius but also moves out of the focal plane in the direction of the source. For certain parameters, it is demonstrated that with increasing orbital number, the maximum amplitude in a focused beam becomes localized along a spheroidal surface enclosing a shadow zone in the prefocal region. This field structure is described analytically by ray theory developed in the present work, showing that the spheroidal surface in the prefocal region coincides with a simple expression for the coordinates of the caustic surface formed in a focused vortex beam.

Visualization of a Toroidal Vortex Structure Forming in the Laboratory Ball Lightning of the Gatchina Discharge

Visualization of a Toroidal Vortex Structure Forming in the Laboratory Ball Lightning of the Gatchina Discharge

Fully resolved simulation of spherical and non-spherical particles in a turbulent channel flow

This paper investigates the behavior of finite-size particles in a turbulent channel flow using a custom direct numerical simulation solver within the FOAM-Extend framework. The solver integrates the cut-cell immersed boundary module with a Lagrangian particle-tracking subroutine capable of simulating the motion of both spherical and non-spherical particles. The study investigates the complex interactions between particles and turbulent structures, offering insights into how particle shape and orientation affect their behavior within the flow field. Additionally, it examines the collision dynamics of two spherical particles in a turbulent channel. The simulations reveal that particle shape significantly influences particle trajectories, rotation, and their interactions with turbulent structures. The Q-criterion visualization showed the creation of hairpin and vortex ring structures shed by the particles. Finally, the close proximity and collision of particles was shown to significantly modify the flow pattern and particle dynamics.

Numerical simulation investigating the impact of regulated underflow rate on the performance of a cyclone with split flow

A novel cyclone design called enhanced cyclone with split flow (ECSF) incorporates bypass flow and underflow to eliminate or suppress localized secondary flow observed in conventional cyclones. The present study aims to investigate the mechanism of ECSF in suppressing localized secondary flow and explore methods for regulating the bottom flow rate. The turbulent characteristics within the ECSF are obtained using Reynolds stress model, while the trajectory of particles is predicted using discrete phase model to derive the separation performance. The accuracy of the simulation model has been validated through experimental verification. The results demonstrated that the full separation efficiency of ECSF for PM10 and PM2.5 remained above 90 % and 70 %, respectively. Moreover, it exhibited an increasing trend with the rise in the full split ratio, eventually stabilizing at approximately 95 % and 77 %, correspondingly, when the full split ratio reached 20 %. The bypass flow and underflow effectively eliminated or suppressed the upper ash ring and processing vortex core (PVC) phenomenon, respectively. However, the exacerbation of the short-circuit flow phenomenon hindered further improvement in separation efficiency beyond a full diversion ratio of 20 %. When the full split ratio is relatively small, particles may accumulate in the discharge pipe, which can be resolved by introducing auxiliary airflow; however, it is crucial to ensure that the flow rate of auxiliary flow does not exceed 14.29 % of that of feed flow to avoid compromising the separation efficiency. Furthermore, the split flow induces an additional reversed vortex within the ECSF, which extends from the overflow pipe to conical section. Its structural characteristics are influenced by both back pressure of discharge outlet and diameter of underflow pipe. To ensure optimal performance of the ECSF, it is recommended that the underflow pipe diameter should not be smaller than that of the outlet in conical section.

Effect of the tip speed ratio on the wake characteristics of wind turbines using LBM‐LES

AbstractIn this paper, the wake characteristics of Zell 2000 wind turbine under different tip velocity ratios are studied by using the lattice Boltzmann method and large eddy simulation. The adaptive mesh refinement method is performed to capture the fine flow structure and wake characteristics development. In this paper, we mainly focus on the effect of the tip speed ratio on the flow structure and unsteady characteristics of wind turbine wake. The three‐dimensional flow vorticity structures, the section vorticity diagram, the pressure fluctuation of wake and the lift coefficient of wind turbine wake are utilized to explore the effect of the tip speed ratio on the unsteady physics mechanism of wind turbine wake. With the increase of the tip speed ratio, the distance between two adjacent vortex rings along the axial direction gradually decreases, as the position of the broken vortex circles gradually approaches the center of the blade, separated vortexes are rapidly generated, and the coherent structure appears closer to the wind turbine. A relationship is established between the tip speed ratios and the positions of the broken vortex circles. It is further found that the dominant frequency amplitude gradually increases with the increase of tip speed ratio and the pressure amplitude spectra of vortex increases with the decrease of the distance between the wake and the center of the blade axis. The above series of studies can provide significant physical insight into deep understanding the influence of the tip speed ratio on the wake characteristics of wind turbines.

Effect of stratification on the propagation of a cylindrical gravity current

Direct numerical simulations (DNSs) of three-dimensional cylindrical release gravity currents in a linearly stratified ambient are presented. The simulations cover a range of stratification strengths $0< S\leq 0.8$ (where $S=(\rho _b^*-\rho _0^*)/(\rho _c^*-\rho _0^*), \rho _b^*, \rho _0^*$ and $\rho _c^*$ are the dimensional density at the bottom of the domain, top of the domain and the dense fluid, respectively) at two different Reynolds numbers. A comparison between the stratified and unstratified cases illustrates the influence of stratification strength on the dynamics of cylindrical gravity currents. Specifically, the front velocity in the slumping phase decreases with increasing stratification strength whereas the duration of the slumping phase increases with increments of $S$ . The Froude number calculated in this phase shows a good agreement with models proposed by Ungarish & Huppert (J. Fluid Mech., vol. 458, 2002, pp. 283–301) and Ungarish (J. Fluid Mech., vol. 548, 2006, pp. 49–68), originally developed for planar gravity currents in a stratified ambient. In the inertial phase, the front velocity across cases with different stratification strengths adheres to a power-law scaling with an exponent of $-$ 1/2. Higher Reynolds numbers led to more frequent lobe splitting and merging, with lobe size diminishing as stratification strength increased. Strong interactions among inner vortex rings occurred during the slumping phase, leading to the early formation of hairpin vortices in weakly stratified cases, while strongly stratified cases exhibited delayed vortex formation and less turbulence.

Curvature conservation and conduction modulation for symmetric charged ferroelectric domain walls

Curvature conservation is a common feature for the physical properties with topological constraints, from the macroscopic cosmic string loops to the microscopic quantized vortex rings in the superfluid. Interestingly, ferroelectric domain-wall curvature is important to determine its conduction, which is the key parameter for designing domain-wall nanoelectronic devices. Here, we demonstrate the curvature conservation of the charged domain walls confined in ferroelectric topological structures with specific domain symmetry by combining piezoelectric force microscopy, conducting atomic force microscopy, and phase-field simulations. Significantly, the intrinsic conductivity of the charged domain walls exhibits an inverse functional relationship with the local domain-wall curvature in the ferroelectric BiFeO3 nano-islands, which can also be precisely and continuously modulated by applying an electric field. This work provides us with an insightful understanding of the conduction mechanism of charged domain walls confined in ferroelectric topological structures with a high degree of spatial symmetry, which paves the way for the application of charged domain walls as flexible conducting channels with continuously changing resistance states to be used for ferroelectric domain-wall memristors.

Interaction between an elliptic vortex ring and a bubble: effect of capture angle

Interaction between an elliptic vortex ring and a bubble: effect of capture angle

Rotor Performance Predictions for Urban Air Mobility: Single vs. Coaxial Rigid Rotors

This work details the development and validation of a methodology for high-resolution rotor models used in hybrid Blade Element Momentum Theory Unsteady Reynolds Averaged Navier–Stokes (BEMT-URANS) CFD. The methodology is shown to accurately predict single and coaxial rotor performance in a fraction of the time required by conventional CFD methods. The methodology has three key features: (1) a high-resolution BEMT rotor model enabling large reductions in grid size, (2) a discretized set of momentum sources to interface between the BEMT rotor model and the structured URANS flow solver, and (3) leveraging of the first two features to enable highly parallelized GPU-accelerated multirotor CFD simulations. The hybrid approach retains high-fidelity rotor inflow, wake propagation, and rotor–rotor interactional effects at a several orders of magnitude lower computational cost compared to conventional blade-resolved CFD while retaining high accuracy on steady rotor performance metrics. Rotor performance predictions of thrust and torque for both single and coaxial rotor configurations are compared to test the data that the authors obtained at the NASA Langley 14- by 22-ft. Subsonic Tunnel Facility. Simulations were run with both fully turbulent and free-transition airfoil performance tables to quantify the associated uncertainty. Single rotor thrust and torque were predicted on average within 4%. Coaxial thrust and power were predicted within an average of 5%. A vortex ring state (VRS) shielding phenomenon for coaxial rotor systems is also presented and discussed. The results support that this hybrid BEMT-URANS CFD methodology can be highly parallelized on GPU machines to obtain accurate rotor performance predictions across the full spectrum of possible UAM flight conditions in a fraction of the time required by conventional higher-fidelity methods. This strategy can be used to rapidly create look-up tables with hundreds to thousands of flight conditions using a three-dimensional multirotor CFD for UAM.

Current-Induced Creation of Topological Vortex Rings in a Magnetic Nanocylinder.

Vortex rings are ubiquitous topological structures in nature. In solid magnetic systems, their formation leads to intriguing physical phenomena and potential device applications. However, realizing these topological magnetic vortex rings and manipulating their topology on demand have still been challenging. Here, we theoretically show that topological vortex rings can be created by a current pulse in a chiral magnetic nanocylinder with a trench structure. The creation process involves the formation of a vortex ring street, i.e., a chain of magnetic vortex rings with an alternative linking manner. The created vortex rings can be bounded with monopole-antimonopole pairs and possess a rich and controllable linking topology (e.g., Hopf link and Solomon link), which is determined by the duration and amplitude of the current pulse. Our proposal paves the way for the realization and manipulation of diverse three-dimensional (3D) topological spin textures and could catalyze the development of 3D spintronic devices.

Pipe to annular flow conversion: Numerical study

Large Eddy Simulation (LES) of turbulent flow at the entrance to an annular pipe section has been performed. The flow over a single circular rod, coaxially immersed in a circular pipe is considered. The Reynolds number, based on mean velocity and hydraulic diameter of the annular section, is ReDh=40530. The study aims to gain insights into the turbulence characteristics of the primary and secondary separation bubbles, developing on a side surface of the rod. Key findings include tiny toroidal vortices in the primary bubble's front, forming a vortex train with intermittently shed vortex structures. Large-scale helical stripes have been reproduced at the entrance to the annular section. Tilting of the axially-oriented streaks was observed close to the rod surface. It resulted in a slight anisotropy of the normal Reynolds stresses. The LES data could refine classic Reynolds-averaged Navier-Stokes-based turbulence models.

Improving depth uncertainty in plenoptic camera-based velocimetry

This work describes the development of a particle tracking velocimetry (PTV) algorithm designed to improve three-dimensional (3D), three-component velocity field measurements using a single plenoptic camera. Particular focus is on mitigating the longstanding depth uncertainty issues that have traditionally plagued plenoptic particle image velocimetry (PIV) experiments by leveraging the camera’s ability to generate multiple perspective views of a scene in order to assist both particle triangulation and tracking. 3D positions are first estimated via light field ray bundling (LFRB) whereby particle rays are projected into the measurement volume using image-to-object space mapping. Tracking is subsequently performed independently within each perspective view, providing a statistical amalgamation of each particle’s predicted motion through time in order to help guide 3D trajectory estimation while simultaneously protecting the tracking algorithm from physically unreasonable fluctuations in particle depth positions. A synthetic performance assessment revealed a reduction in the average depth errors obtained by LFRB as compared to the conventional multiplicative algebraic reconstruction technique when estimating particle locations. Further analysis using a synthetic vortex ring at a magnification of − 0.6 demonstrated plenoptic-PIV capable of maintaining the equivalent of 0.1–0.15 voxel accuracy in the depth domain at a spacing to displacement ratio of 5.3–10.5, an improvement of 84–89% compared to plenoptic-PIV. Experiments were conducted at a spacing to displacement ratio of approximately 5.8 to capture the 3D flow field around a rotor within the rotating reference frame. The resulting plenoptic-PIV/PTV vector fields were evaluated with reference to a fixed frame stereoscopic-PIV (stereo-PIV) validation experiment. A systematic depth-wise (radial) component of velocity directed toward the wingtip, consistent with observations from prior literature and stereo-PIV experiments, was captured by plenoptic-PTV at magnitudes similar to the validation data. In contrast, the plenoptic-PIV did not discern any coherent indication of radial motion. Our algorithm constitutes a significant advancement in enhancing the functionality and versatility of single-plenoptic camera flow diagnostics by directly addressing the primary limitation associated with plenoptic imaging.Graphical abstract

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.

Monge–Ampère geometry and vortices

We introduce a new approach to Monge–Ampère geometry based on techniques from higher symplectic geometry. Our work is motivated by the application of Monge–Ampère geometry to the Poisson equation for the pressure that arises for incompressible Navier–Stokes flows. Whilst this equation constitutes an elliptic problem for the pressure, it can also be viewed as a non-linear partial differential equation connecting the pressure, the vorticity, and the rate-of-strain. As such, it is a key diagnostic relation in the quest to understand the formation of vortices in turbulent flows. We study this equation via an associated (higher) Lagrangian submanifold in the cotangent bundle to the configuration space of the fluid. Using our definition of a (higher) Monge–Ampère structure, we study an associated metric on the cotangent bundle together with its pull-back to the (higher) Lagrangian submanifold. The signatures of these metrics are dictated by the relationship between vorticity and rate-of-strain, and their scalar curvatures can be interpreted in a physical context in terms of the accumulation of vorticity, strain, and their gradients. We show explicity, in the case of two-dimensional flows, how topological information can be derived from the Monge–Ampère geometry of the Lagrangian submanifold. We also demonstrate how certain solutions to the three-dimensional incompressible Navier–Stokes equations, such as Hill’s spherical vortex and an integrable case of Arnol’d–Beltrami–Childress flow, have symmetries that facilitate a formulation of these solutions from the perspective of (higher) symplectic reduction.

An investigation of flow regimes and mixing in a novel arrow-shaped jet reactor

Fluid mixing plays a vital role in the industry of chemical engineering. In this work, a novel arrow-shaped jet reactor by varying the angle of fluid impingement in the inlet channel was investigated. Flow field visualization experiments using Planar Laser-Induced Fluorescence (PLIF) technique were carried out to study various flow regimes and mixing performance. Results show that differentiated flow patterns of engulfment flow regimes appear within the arrow-shaped jet reactor. Particularly, vortex ring shedding was found in the unsteady engulfing flow regime, which promoted the mixing significantly. Besides, the mechanism of the variation of fluid impinging angle on mixing in an arrow-shaped jet reactor was investigated by numerical simulation. The results indicate that the transformation from double-leg vortex to single-leg vortex and the increase of the volume of vortex clusters caused by vortex ring shedding are the main reasons for the improvement of the mixing of engulfment flow regime.

Velocity gradient analysis of a head-on vortex ring collision

We simulate the head-on collision between vortex rings with circulation Reynolds numbers of 4000 using an adaptive, multiresolution solver based on the lattice Green's function. The simulation fidelity is established with integral metrics representing symmetries and discretization errors. Using the velocity gradient tensor and structural features of local streamlines, we characterize the evolution of the flow with a particular focus on its transition and turbulent decay. Transition is excited by the development of the elliptic instability, which grows during the mutual interaction of the rings as they expand radially at the collision plane. The development of antiparallel secondary vortex filaments along the circumference mediates the proliferation of small-scale turbulence. During turbulent decay, the partitioning of the velocity gradients approaches an equilibrium that is dominated by shearing and agrees well with previous results for forced isotropic turbulence. We also introduce new phase spaces for the velocity gradients that reflect the interplay between shearing and rigid rotation and highlight geometric features of local streamlines. In conjunction with our other analyses, these phase spaces suggest that, while the elliptic instability is the predominant mechanism driving the initial transition, its interplay with other mechanisms, e.g. the Crow instability, becomes more important during turbulent decay. Our analysis also suggests that the geometry-based phase space may be promising for identifying the effects of the elliptic instability and other mechanisms using the structure of local streamlines. Moving forward, characterizing the organization of these mechanisms within vortices and universal features of velocity gradients may aid in modelling turbulent flows.

Fluid–acoustic interactions around an expanding pipe with orifice plates

Expanding pipes with orifice plates are often utilized as silencers for fluid machinery. However, intense tonal sounds can be generated from a flow through such expanding pipes. To clarify the mechanism of tonal sound from a flow through a circular expanding pipe with two orifice plates and the conditions for intense acoustic radiation, the flow and acoustic fields are directly solved based on the compressible Navier–Stokes equations. Phase-averaged flow fields indicate the occurrence of periodic vortex shedding in the free shear layers of the expanding pipe, resulting in acoustic radiation. The effects of the orifice radius and freestream Mach number on the acoustic radiation are focused on. The computational results demonstrate that vortex rings or spiral vortices are generated in the cavities formed by the orifice plates, where the primary vortical shape changes, depending on the freestream Mach number and orifice radius. The collision of the vortex ring and spiral vortex with the orifice plate or downstream edge of the expanding pipe leads to the occurrence of circumferentially in-phase and one-wavelength-mode pressure fluctuations, respectively. The orifice radius also affects the convective velocity of vortices and the position of the acoustic source, varying the frequency of the acoustic radiation. The findings of this research provide the first clarifications of fluid–acoustic interactions in an expanding pipe with orifice plates.