No Arabic abstract
Cloud cavitation is related to an intrinsic instability where clouds are shed periodically. The shedding process is initiated either by the motion of a liquid re-entrant jet or a condensation shock. Cloud cavitation in nozzles interacts with the flow field in the nozzle, the mass flow and the spray break-up, and causes erosion damage. For nozzle geometries cloud shedding and the associated processes have not yet been studied in detail. In this paper, we investigate the process of cloud cavitation shedding, the re-entrant jet and condensation shocks in a scaled-up generic step nozzle with injection into gas using implicit Large-Eddy Simulations (LES). For modeling of the cavitating liquid we employ a barotropic equilibrium cavitation model, embedded in a homogeneous multi-component mixture model. Full compressibility of all components is taken into account to resolve the effects of collapsing vapor structures. We carry out simulations of two operating points exhibiting different cavitation regimes. The time-resolved, three-dimensional simulation results cover several shedding cycles and provide deeper insight into the flow field. Our results show that at lower cavitation numbers, shedding is initiated by condensation shocks, which has not yet been reported for nozzle flows with a constant cross-section. We analyze the cavitation dynamics and the shedding cycles of both operating points. Based on our observations we propose modifications to established schematics of the cloud shedding process. Additionally, we analyze the near-wall upstream flow in and underneath the vapor sheet and possible driving mechanism behind the formation of the re-entrant jet.
We present results of large eddy simulations of a cavitating nozzle flow and injection into gas, investigating the interactions of cavitation in the nozzle, primary jet breakup, mass-flow rates, and gas entrainment. During strong cavitation, detached vapor structures can reach the nozzle outlet, leading to partial entrainment of gas from the outflow region into the nozzle. The gas entrainment can affect cavitation dynamics, mass-flow rates, and jet breakup. Moreover, the implosion of detached vapor structures induces pressure peaks that on the one hand amplify turbulent fluctuations and subsequently can enhance jet breakup and on the other hand can damage walls in the proximity and thus lead to cavitation erosion. Our numerical setup is based on a reference experiment, in which liquid water is discharged into ambient air through a step nozzle. The cavitating liquid and the non-condensable gas phase are modeled with a barotropic homogeneous mixture model while for the numerical model a high-order implicit large eddy approach is employed. Full compressibility of all components is taken into account, enabling us to capture the effects of collapsing vapor structures. Two operating points covering different cavitation regimes and jet characteristics are investigated. Special emphasis is placed on studying the effects of cavitation on the mass flow and the jet as well as the impact of partial gas entrainment. Therefore, frequency analyses of the recorded time-resolved signals are performed. Furthermore, the dynamics and intensities of imploding vapor structures are assessed.
A nonlocal subgrid-scale stress (SGS) model is developed based on the convolution neural network (CNN), a powerful supervised data-driven approach. The CNN is an ideal approach to naturally consider nonlocal spatial information in prediction due to its wide receptive field. The CNN-based models used here only take primitive flow variables as input, then the flow features are automatically extracted without any $priori$ guidance. The nonlocal models trained by direct numerical simulation (DNS) data of a turbulent channel flow at $Re_{tau}=178$ are accessed in both the $priori$ and $posteriori$ test, providing physically reasonable flow statistics (like mean velocity and velocity fluctuations) closing to the DNS results even when extrapolating to a higher Reynolds number $Re_{tau}=600$. In our model, the backscatter is also predicted well and the numerical simulation is stable. The nonlocal models outperform local data-driven models like artificial neural network and some SGS models, e.g. the Smagorinsky model in actual large eddy simulation (LES). The model is also robust since stable solutions can be obtained when examining the grid resolution from one-half to double of the spatial resolution used in training. We also investigate the influence of receptive fields and suggest using the two-point correlation analysis as a quantitative method to guide the design of nonlocal physical models. To facilitate the combination of machine learning (ML) algorithms to computational fluid dynamics (CFD), a novel heterogeneous ML-CFD framework is proposed. The present study provides the effective data-driven nonlocal methods for SGS modelling in the LES of complex anisotropic turbulent flows.
We have performed Coherent-vorticity Preserving Large-Eddy simulations of a trefoil knot-shaped vortex, inspired by the experiments of Kleckner and Irvine. The flow parameter space is extended in the present study, including variations of the circulation Reynolds numbers in the range Re = 2000 - 200000, where Re = 20000 is the value used in the experiments. The vortex line corresponding to the trefoil knot is defined using a parametric equation and the Biot-Savart law is employed to initialize the velocity field. The CvP LES computation displays a good qualitative match with the experiment. In particular, the vortex entanglement process is accurately represented as well as the subsequent separation of the main vortex in two distinct structures - a small and a large vortex - with different self-advection speeds that have been quantified. The small vortex propagates faster than the large oscillatory vortex which carries an important amount of vorticity. The advection velocity of the vortex before bursting is found to be independent of the Reynolds number. The low Reynolds number computation leads to a decrease of the separated vortices velocity after bursting, compared to the higher Reynolds computations. The computation of energy spectra emphasizes intense energy transfers from large to small scales during the bursting process. The evolution of volume-averaged enstrophy shows that the bursting leads to the creation of small scales that are sustained a long time in the flow, when a sufficiently large Reynolds number is considered (Re>20000). The low Reynolds number case Re = 2000 hinders the generation of small scales during the bursting process and yields essentially laminar dynamics. The onset of background turbulence due to the entanglement process can be observed at Re = 200000
This study aims to quantify how turbulence in a channel flow mixes momentum in the mean sense. We applied the macroscopic forcing method to direct numerical simulation (DNS) of a turbulent channel flow at Re$_tau$=180 using two different forcing strategies that are designed to separately assess the anisotropy and nonlocality of momentum mixing. In the first strategy, the leading term of the Kramers-Moyal expansion of the eddy viscosity operator is quantified where the macroscopic forcing is employed to reveal all 81 tensorial coefficients that essentially represent the local-limit eddy viscosity. The results indicate: (1) eddy viscosity has significant anisotropy, (2) Reynolds stresses are generated by both mean strain rate and mean rotation rate tensors, and (3) the local-limit eddy viscosity generates asymmetric Reynolds stress tensors. In the second strategy, the eddy viscosity operator is considered as an integration kernel representing the nonlocal influence of mean gradients on the Reynolds stresses. Considering the average of this kernel in the homogeneous directions, the macroscopic forcing is designed to reveal the nonlocal effects in the wall-normal direction for all 9 components of the Reynolds stresses. Our results indicate that while the shear component of the Reynolds stress is reasonably controlled by the local mean gradients, other components of the Reynolds stress are highly nonlocal. These two analyses provide accurate verification data for quantitative testing of anisotropy and nonlocality effects in turbulence closure models.
The flow structure obtained when Localized Arc Filament Plasma Actuators (LAFPA) are employed to control the flow issuing from a perfectly expanded Mach 1.3 nozzle is elucidated by visualizing coherent structures obtained from Implicit Large-Eddy Simulations. The computations reproduce recent experimental observations at the Ohio State University to influence the acoustic and mixing properties of the jet. Eight actuators were placed on a collar around the periphery of the nozzle exit and selectively excited to generate various modes, including first and second mixed (m = +/- 1 and m = +/- 2) and axisymmetric (m = 0). In this fluid dynamics video http://ecommons.library.cornell.edu/bitstream/1813/13723/2/Alljoinedtotalwithmodetextlong2-Datta%20MPEG-1.m1v, http://ecommons.library.cornell.edu/bitstream/1813/13723/3/Alljoinedtotalwithmodetextlong2-Datta%20MPEG-2.m2v}, unsteady and phase-averaged quantities are displayed to aid understanding of the vortex dynamics associated with the m = +/- 1 and m = 0 modes excited at the preferred column-mode frequency (Strouhal number 0.3). The unsteady flow in both contains a broad spectrum of coherent features. For m = +/- 1, the phase-averaged flow reveals the generation of successive distorted elliptic vortex rings with axes in the flapping plane, but alternating on either side of the jet axis. This generates a chain of structures where each interacts with its predecessor on one side and its successor on the other. Through self and mutual interaction, the leading segment of each loop is pinched and passes through the previous ring before rapidly breaking up, and the mean jet flow takes on an elliptic shape. The m = 0 mode exhibits relatively stable roll-up events, with vortex ribs in the braid regions connecting successive large coherent structures.