Large Eddy Simulation Method Research Articles

Re-entry rocket basic flow characteristics and thermal environment of different retro-propulsion modes

During the supersonic re-entry of multi-nozzle heavy rockets into the atmosphere, the basic flow state becomes increasingly complex due to the coupling effect between the retro-propulsion plumes and the freestream. A numerical method using the hybrid Reynolds-Averaged Navier-Stokes and Large Eddy Simulation (RES) method and discrete coordinate method is developed to accurately estimate the thermal environment. In addition, finite rate chemical kinetics is used to calculate the afterburning reactions. The numerical results agree well with wind tunnel data, which confirms the validity and accuracy of the numerical method. Computations are conducted for the heavy carrier rocket re-entry from 53.1 km to 39.5 km altitude with 180° angle of attack by using three different Supersonic Retro-Propulsion (SRP) modes. The numerical results reveal that these three SRP flow fields are all Short Penetration Models (SPM). As the re-entry altitudes decrease, both the plume-plume interaction and the plume-freestream interaction become weaker. The highest temperatures in the plume shear layers of the three SRP modes increase by 8.36%, 7.33% and 6.92% respectively after considering afterburning reactions, and all occur at a re-entry altitude of 39.5 km. As the rocket re-enters the atmosphere, the maximum heat flux on the rocket base plate of three SRP modes stabilizes at 290, 170 and 200 kW/m2 respectively, but the maximum heat flux on the side wall increases significantly. When the altitude declines to 39.5 km, the extreme heat flux of the three modes increase by 84.16%, 49.45% and 62.97% respectively compared to that at 53.1 km.

Research on wake characteristics and HCWF blade number matching of a new coupled propeller

For the purpose of wake energy recovery and hub vortex control in a new Coupled Propeller (CP), three hub cap with fins (HCWF) schemes with different blade numbers are designed for a specific multi-blade CP. A numerical investigation of the CP-HCWF system is conducted using the Large Eddy Simulation (LES) method. Firstly, the analysis focuses on the distribution of axial and tangential flow velocity characteristics of the CP. The wake spectral characterization analysis results reveal that the 10th-order component is equivalent to the 5th-order components within 0.7 R of the HCWF. However, beyond 0.7 R of the HCWF, the 10th-order component becomes more significant. Then, the influence of the HCWF blade number on the CP performance is further analyzed. Hydrodynamic efficiency, post-hub pressure recovery, wake energy recovery, and hub vortex elimination improve with increased HCWF blade numbers. The 10-blade HCWF demonstrates the best performance. Ultimately, the reliability of the numerical findings is verified through the model test. The efficiency variation with the fin number exhibits a consistent trend in the test and calculation results. It is suggested that the HCWF blade number is selected by considering the sum of the CP blades as a reference.

Study of hydrodynamic stress in cell culture bioreactors via lattice-Boltzmann CFD simulations supported by micro-probe shear stress method

Mammalian cell cultivation in pharmaceutical industry can last up to units of weeks and requires proper transport of nutrients and oxygen for cell growth and production. Given the long time period, cells experience flow fields from all bioreactor’s zones, where the energy dissipation rate (ε) varies substantially. Shear sensitive micro-probes with size comparable to cells and Kolmogorov eddies are used for the determination of the maximum hydrodynamic stress (τmax) in bioreactors. For the very first time, the micro-probe method is applied successfully not only to laboratory (3 L) and pilot scale (80 L and 200 L), but also to industrial production scale bioreactor (12,500 L) with Rushton turbine and pitched-blade (RT-PB) impeller configuration. Experimentally obtained data are used for the validation of comprehensive CFD scale-up study, using the Lattice-Boltzmann large eddy simulation (LB-LES) method. Besides τmax, this work also focuses on the study of mixing time and flow field attributes.

A wall-modeled immersed boundary/large eddy simulation method and its application to simulating heart valve flows

In this work, a wall-modeled immersed boundary (IB)/large eddy simulation (LES) method is extended to the simulation of moving-boundary flows. The used non-equilibrium algebraic wall model is based on an assumed velocity profile, the coefficients of which are determined from physical constraints provided by the full turbulent-boundary-layer equations. To implement the wall model in an IB method named the local domain-free discretization (DFD) method, a local coordinate system fixed on the moving body is introduced. Thus, wall modeling is transformed into a local two-dimensional problem and the complexity of implementation of the wall model is reduced. In the present LES-DFD method, the tangential velocity at an exterior dependent node is determined via wall shear stress prescribed by the wall model. To reduce computational cost for simulating an internal flow with moving boundaries, the stationary boundaries are handled by the body-fitted-grid method and the moving boundaries by the local DFD method. There is no need of an auxiliary grid for solving the non-equilibrium algebraic wall model. Therefore, the inbuilt advantage of an IB method can be retained when simulating moving-boundary problems, and the economy of equilibrium wall models can also be preserved. The present method is applied to simulating the pulsatile flows through a bileaflet mechanical heart valve implanted in a model aorta. The predicted results show an acceptable agreement with the referenced experimental measurements or numerical results at much higher resolution and the applicability of the non-equilibrium wall model to LES of complex moving-boundary flows is verified.

Autonomous large-eddy simulations of turbulence using eddy viscosity derived from the subgrid-scale similarity stress tensor

A previously developed method for large-eddy simulations (LES), based on spectral eddy-viscosity models, is generalised to the physical space representation. The method estimates the subgrid-scale (SGS) energy transfer using a similarity-type model expression for the SGS tensor obtained using Gaussian filtering of velocity fields advanced in the simulations. Following steps for the spectral space representation, the SGS transfer in the physical space is used to obtain a spatially varying eddy viscosity at each time step in LES. The computed eddy viscosity is employed to model the SGS stress tensor in the familiar Boussinesq form for use in LES. The method is tested in LES of isotropic turbulence at high Reynolds numbers where the inertial range dynamics is expected and for lower Reynolds number decaying turbulence under conditions of the classical Comte-Bellot and Corrsin experiments. In both cases the agreement with reference data is very good and the SGS transfer computed for the proposed eddy-viscosity model is highly correlated with the transfer computed for the similarity-type stress tensor.

A new probability-distribution scale synthetic eddy method for large eddy simulation of wind loads

The accurate generation of inflow turbulence in large eddy simulation (LES) plays a crucial role in assessing the wind loads on buildings. To faithfully recreate turbulent flow fields, the inflow turbulence must accurately replicate vortex structures. In this study, we propose a novel method for generating turbulent inflow called the probability-distribution scale synthetic eddy method (PDSSEM). PDSSEM utilizes a random generation of eddy positions, with the length of each eddy scale determined by its probability density function (PDF). This approach ensures that PDSSEM maintains consistency in the mean wind speed profile, turbulence integral length, turbulence intensity profile, and wind speed spectra of the atmospheric boundary layer (ABL). The PDSSEM method can be more easily applied compared with multi-scale synthetic eddy method (MSSEM), making it highly versatile and adaptable. The self-sustaining characteristics of the flow field generated by PDSSEM are validated, revealing abundant vortex structures and good agreement with targeted values for turbulence flow characteristics. Furthermore, comparing the wind pressures of a tall building with the experimental data, the accuracy and feasibility of PDSSEM are comprehensively demonstrated. Thus, it indicates the great potential of the proposed method for broader computational wind engineering applications, including wind load evaluation and flow pattern investigation.

A mass transfer-based LES modelling methodology for analyzing the movement of submarine sediment flows with extensive shear behavior

This study employs the large-eddy simulation (LES) method to investigate the interaction of submarine sediment flows with extensive shear behavior and ambient water. This method is validated with good accuracy by simulating the head velocity and evolution geometries of a submarine mud flow with non-Newtonian fluid characteristics and a submarine turbidity current with Newtonian fluid characteristics in inclined and regular lock-exchange flume experiments. This study finds that the violent mass transfer at the interface of the submarine sediment flow and seawater is caused by numerous eddies, which significantly alter the submarine sediment flow's geometry, resulting in varying degrees of undulation, including the deposition pattern of the submarine sediment flow tail. The acceleration zone, where the velocity of the submarine sediment flow increases significantly at this undulation, propels the sediment flow forward, supporting its long-distance transport. Furthermore, the turbulence and mass transport characteristics of submarine turbidity currents with low dynamic viscosity Newtonian fluid characteristics are stronger than those of submarine mud flows with high dynamic viscosity non-Newtonian fluid characteristics. Therefore, when submarine landslides develop into later stages, such as submarine turbidity currents with high velocity, large volume, and long run-out distance characteristics, more attention must be given to the mass transport process.

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

LES of autoignition process in a turbulent droplet-laden mixing layer – impact of an evaporation model and discretisation method

This paper focuses on the numerical modelling of the autoignition process in an ethanol-air mixture within a turbulent droplet-laden temporally evolving mixing layer, utilising the Large Eddy Simulation (LES) method. We compare results obtained by implementing two notably distinct discretisation methods (second-order Total Variation Diminishing – TVD, and fifth-order Weighted Essentially Non-Oscillatory – WENO) in the species and energy transport equations. Additionally, we examine three widely recognised evaporation models (the ‘ D 2 ’ model, the Abramzon–Sirignano model, and a non-equilibrium model based on the Langmuir-Knudsen law). For modelling the combustion process, we employ the Eulerian Stochastic Field method and a laminar chemistry model. The simulations encompass two droplet diameter distributions, one with low and the other with significant diameter non-uniformity. Furthermore, nine distinct initial conditions are considered. Our findings indicate that adjusting the discretisation scheme or the evaporation model can yield simulation outcomes differing to a substantial degree, similar to the impact of modifying the actual properties of the spray droplets or their initial distribution. This is particularly pronounced in the case of sprays featuring a narrowed droplet diameter distribution. We found that the Abramzon–Sirignano model results in the most rapid autoignition. The TVD scheme contributes to a more even spatial distribution of the evaporated fuel, consequently facilitating quicker evaporation and faster autoignition. Using the WENO scheme generates higher fuel vapor gradients near the droplets, leading to less intense evaporation and slower fuel transport in space.

Numerical study on the influence of riblet yaw-angle on drag reduction rate with different immersion heights

This paper uses the Large-eddy simulation method to investigate the influence of riblet yaw angle on drag reduction rate with different immersion heights, revealing that the vortex structure is gradually transformed from riblet-tip vortex pairs to boundary vortex forms with the increase of yaw angle, and analyzing the influence mechanism of immersion height on the critical yaw angle. The result shows that the drag-reduction rate of riblets decreases as the yaw angle increases, and the critical angle of the riblet varies with the different immersion heights. When the immersion height is 0.7, 0.5, and 0.3, the critical yaw angle is close to 30, 40, and 50°, respectively. As the immersion height increases, the cross-flow velocity at the riblet tip is more accessible to reach the critical velocity at which a boundary vortex can be formed, resulting in a smaller critical yaw angle for larger immersion heights. Furthermore, with the yaw angle increases, the normal diffusivity of normal vorticity and spanwise vorticity increase and decrease, respectively, which indicates that the lateral flow fluctuations decrease beyond a critical yaw angle, while the Reynolds stress component due to energy dissipation induced by the streamwise velocity decreases. Meanwhile, with the formation of a boundary vortex, the ability of the riblet to reduce viscous shear stresses increases due to the reduction of near-wall velocity gradient, and the emergence of a high-pressure region on the windward side of the riblet and a low-pressure region on the leeward side, consequently augmenting the differential pressure drag. This work provides valuable insights for designing riblets with heightened effectiveness at larger yaw angles, thereby extending their applicability under diverse flow conditions.

Effects of split injection strategy on combustion characteristics and NOx emissions performance in dual-fuel marine engine

In order to address environmental concerns and promote sustainable development, the maritime industry is facing a demand to reduce emissions from ship engines. To address this, a growing trend in the development of diesel natural gas dual-fuel marine engine involves implementing more flexible diesel injection strategies, with the aim of enhancing performance and minimizing harmful emissions. This study utilized a numerical investigation based on the Large Eddy Simulation (LES) method to examine pilot diesel injection strategies in a dual-fuel marine engine operating at a high degree of natural gas substitution rate. Specifically, the analysis focused on the start of injection (SOI) in the single injection strategy and the first start of injection (SOI1) in the split injection strategy for pilot diesel. The findings highlight distinct trade-off relationships between these two injection strategies regarding their impact on indicated specific fuel consumption (ISFC) and NOx emissions. Compared to the optimized single injection strategy (Case A: SOI = -15°CA ATDC), the optimized split injection strategy (Case B: SOI1 = -60°CA ATDC, SOI2 = -10°CA ATDC) increases thermal efficiency by 0.77 % while reducing NOx emissions by 36.3 %. The split injection strategy induces low-temperature reactions of n-heptane in the squish region, resulting in the generation of active free radicals such as CH2O and OH, thereby accelerating the propagation speed of the flame in the squish region. Meanwhile, the split injection strategy diminishes pilot diesel accumulation in the piston bowl region, narrows the distribution range of high-temperature areas surpassing 2400 K, consequently leading to a reduction in NOx emissions.

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.

A numerical study on hydrogen blending in natural gas pipeline by a T-Pipe

In order to study the flow blending and transporting process of hydrogen that injects into the natural gas pipelines, a three-dimensional T-pipe blending model is established and the flow characteristics are investigated systematically by the large eddy simulation (LES). Firstly, the mathematical formulation of hydrogen-methane blending process is provided and the LES method is introduced and validated by a benchmark gas blending model having experimental data. Subsequently, the T-pipe blending model is presented, and the effects of key parameters, such as the velocity of main pipe, hydrogen blending ratio, diameter of hydrogen injection pipeline, diameter of main pipe and operating pressure on the hydrogen-methane blending process, are studied systematically. The results show that, under certain conditions, the gas mixture will be stratified downstream of the blending point, with hydrogen at the top of the pipeline and methane at the bottom of the pipeline. In the no-stratified scenario, the mixing distance increases at lower hydrogen mixing ratio and larger diameter of the hydrogen injection pipe or the main pipe. Finally, based on the numerical results, the underlying physics of the stratification phenomenon during the blending process are explored and an indicator for stratification is proposed using the ratio between the Reynolds numbers of the natural gas and hydrogen.

Exploring Vortex–Flame Interactions and Combustion Dynamics in Bluff Body-Stabilized Diffusion Flames: Effects of Incoming Flow Velocity and Oxygen Content

An afterburner encounters two primary features: high incoming flow velocity and low oxygen concentration in the incoming airflow, which pose substantial challenges and contribute significantly to the deterioration of combustion performance. In order to research the influence of oxygen content on the dynamic combustion characteristics of the afterburner under various inlet velocities, the effect of oxygen content (14–23%) on the field structure of reacting bluff body flow, flame morphology, temperature pulsation, and pressure pulsation of the afterburner at different incoming flow velocities (0.1–0.2 Ma) was investigated in this study by using a large eddy simulation method. The results show that two different instability features, BVK instability and KH instability, are observed in the separated shear layer and wake, and are influenced by changes in the O2 mass fraction and Mach number. The oxygen content and velocity affected the oscillation amplitude of the downstream flow. As the O2 mass fraction decreases, the flame oscillation amplitude increases, the OH concentration in the combustion chamber decreases, and the flame temperature decreases. Additionally, the amplitude of the temperature pulsation in the bluff body flame was primarily influenced by the temperature intensity of the flame and BVK instability. Moreover, the pressure pulsation is predominantly affected by the dynamic characteristics of the flow field behind the bluff body. When the BVK instability dominated, the primary frequency of the pressure pulsation aligned with that of the temperature pulsation. Conversely, under the dominance of the KH instability, the temperature pulsation did not exhibit a distinct main frequency. At present, the influence of oxygen content and incoming flow rate on the combustion performance of the combustion chamber is not clear. The study of the effect of oxygen content on the combustion characteristics of the combustion chamber at different incoming flow rates provides a reference for improving the performance of the combustion chamber and enhancing the combustion stability.

Large eddy simulation of the cavity shedding characteristics of ventilated cavitation around the underwater vehicle

Ventilated cavitation is utilized for reducing drag and improving the headway of underwater vehicles. The ventilated cavitating flow around the axisymmetric body is simulated using the Large Eddy Simulation (LES) method in this study. Three cavity morphologies with different ventilation rates (Cq) are observed: Foamy Cavity (FC), Continuous Transparent and Foamy Cavity (CTFC) and Super Cavity (SC). The cavity shedding characteristic of the three different cavity morphologies is revealed. For partial cavity, the periodic cavity shedding is dominated by the local high pressure, re-entrant jet and vortex. The distributions of the turbulent kinetic energy and the pulsation enstrophy exhibit different migration mechanisms during the cavity shedding. The energy characteristic is reflected by three terms of the entropy generation. And the entropy production rate by direct dissipation occupies a dominant position, accounting for about 80% of the total entropy production, while the entropy production rate by viscous dissipation only accounts for 4%. For super cavity, gas leakage is dominated by a re-entrant jet induced by the intense adverse pressure gradient. And the results also demonstrated the velocity and vorticity pulsations around the test body wrapped in a super cavity are not intense.

Validation of numerical simulations and experiments on impulse characteristics induced by self-excited oscillation

The high-frequency pulse flow, equivalent to the natural frequency of rocks, is generated by a self-excited oscillating cavity to achieve resonance rock-breaking. The flow field and oscillating mechanism of the self-excited oscillating cavity were simulated using the large eddy simulation method of Computational Fluid Dynamics (CFD). A field-scale testing apparatus was developed to investigate the impulse characteristics and verify the simulation results. The results show that the fluid at the outlet at the tool is deflected due to the pulse oscillation of the fluid. The size and shape of low-pressure vortices constantly change, leading to periodic changes in fluid impedance within the oscillating cavity. The impulse frequency reaches its highest point when the length–diameter ratio is 0.67. As the length–diameter ratio increases, the tool pressure loss also increases. Regarding the cavity thickness, the impulse frequency of the oscillating cavity initially decreases, then increases, and finally decreases again. Moreover, both the impulse frequency and pressure loss increase with an increase in displacement. The numerical simulation findings align with the experimental results, thus confirming the validity of the theoretical model. This research provides theoretical guidance for the practical application of resonance rock-breaking technology.

Large eddy simulation of combustion characteristics during dual fuel switching process in gas turbine combustor

The present study aims to investigate the combustion characteristics, fuel switching strategy, and flame stability associated with fuel switching processes in a dual-fuel gas turbine combustor. Specifically, the Large Eddy Simulation method is employed in conjunction with a skeletal chemical reaction mechanism to simulate the fuel switching process. It is revealed that the presence of a local rich combustion zone, formed by an increased accumulation of fuel in the head region, leads to the displacement of the heat release position towards the downstream inner shear layer. Additionally, this condition causes a reduction in the recirculation zone area and a deterioration in the uniformity of fuel blending. Simultaneously, the combustor experiences low-frequency, high-amplitude pressure pulsations, as well as pulsations in the heat release rate. Furthermore, a comparative analysis is conducted to assess the discrepancies among three distinct fuel supply strategies employed during fuel switching. The fast-opening strategy is found to demonstrate superior performance indicators, as it accomplishes filling the combustor performance indicators within 80% of the allocated time due to its lower initial fuel changing rate.

Turbulence generation in the transitional wake flow behind a sphere

Turbulence generation in the transitional flow in the wake behind a sphere is studied with numerical simulations. The filtered Navier–Stokes equation and the large eddy simulation method are employed as the governing equation and the numerical method, respectively. The ΩR̃ vortex identification method is used to trace the evolution of vortices in the wake flow. The energy gradient theory is used to analyze the spike formation in the wake flow. The simulation results show that the vortex structure in the wake flow is the type of hairpin vortices, which is similar to that in a boundary layer flow. Ejection and sweep motions exist around the hairpin vortices. There are two most unstable regions in the wake where turbulence “burst” is first produced, one is near the center of the vortex head and the other is between the two vortex legs. There is a high-pressure zone above the vortex head due to the decrease in the streamwise velocity, and a soliton-like coherent structure exists in this area. The mechanism of turbulence generation in the wake is the discontinuity of the streamwise velocity, which makes the Navier–Stokes equation be singular. This singularity leads to the formation of the “negative spike” in the streamwise velocity. The amplitude of the “negative spike” reaches up to 60% of the incoming velocity, which is close to the situation in a boundary layer flow on a flat plate. It is concluded that the mechanism of turbulence generation in the wake flow is the same as that in the boundary layer flow.

Deep learning architecture for sparse and noisy turbulent flow data

The success of deep learning models in fluid dynamics applications will depend on their ability to handle sparse and noisy data accurately. This paper concerns the development of a deep learning model for reconstructing turbulent flow images from low-resolution counterparts encompassing noise. The flow is incompressible through a symmetric, sudden expansion featuring bifurcation, instabilities, and turbulence. The deep learning model is based on convolutional neural networks, in a high-performance, lightweight architecture. The training is performed by finding correlations between high- and low-resolution two-dimensional images. The study also investigates how to remove noise from flow images after training the model with high-resolution and noisy images. In such flow images, the turbulent velocity field is represented by significant color variations. The model's peak signal-to-noise ratio is 45, one of the largest achieved for such problems. Fine-grained resolution can be achieved using sparse data at a fraction of the time required by large-eddy and direct numerical simulation methods. Considering its accuracy and lightweight architecture, the proposed model provides an alternative when repetitive experiments are complex and only a small amount of noisy data is available.

Data-driven modal decomposition of R134a refrigerant cavitating flow in Venturi tube

This study utilized high-speed camera and large eddy simulation methods to explore the cavitating flow mechanisms and turbulence structures of R134a refrigerant inside a Venturi tube under varying cavitation numbers (CNs). Data-driven modal analysis approaches, proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD), were introduced to identify and extract the energy hierarchy and transient characteristics within the cavitating flow. The analysis of grayscale images indicated that the cavitating flow gradually transitioned from quasi-periodic to unsteady flow as the CN decreased, and the severity of cavitation correlates with lower peak frequencies. The POD analysis facilitated the extraction of coherent structures in the cavity's temporal evolution, and the results indicate that the quasi-ordering shedding and collapse of large-scale cavity clouds predominantly occur under low cavitation intensity conditions. As the CN increases, the influence of small-scale cavity shedding becomes more significant. The first 30 most energetic modes occupied over 75% of the entire energy, and they were used to reconstruct the cavitating flow, achieving good consistency with transient flow snapshots. Additionally, the DMD results of the cavitating flow yield three frequency spans, including several prominent characteristic frequencies. These spans are closely linked to the cavity cloud structures of varying scales, unveiling the structural characteristics of unsteady cavitating flow.