For Francis turbines,frequent operations under extremely low load conditions result in significant noise and pressure fluctuation issues.These issues may cause vibration and fatigue damage to the unit,accompanied by d...For Francis turbines,frequent operations under extremely low load conditions result in significant noise and pressure fluctuation issues.These issues may cause vibration and fatigue damage to the unit,accompanied by difficulties in connecting to the grid and reductions in the power generation efficiency of renewable energy.However,there is limited research on the relationship between pressure fluctuations and the induced noise of Francis turbines during extreme operations.In the present study,an acoustic numerical simulation based on the Ffowcs Williams-Hawkings equation and large eddy simulation is used to analyze the acoustic performances of Francis turbines.In the current study,for evaluating the acoustic characteristics under such terrible conditions,the results of variable flow rate and guide vane opening conditions are compared.Results indicated that Francis turbine noise is mostly due to pressure fluctuations brought on by rotor-stator interference and corkscrew-shaped vortices.The blade passing frequency(BPF)of 130.00 Hz and the low frequency of 0.33 f_(n)(where f_(n)denotes the rotating frequency)are the key factors affecting pressure and noise fluctuations.The influence of low frequency is reduced as the flow rate rises,whereas the influence of BPF gradually increases.Besides,the hydrodynamic noise of Francis turbines is primarily low-frequency,with discrete and broad-band features.The rotating noise with distinct peak values and the turbulence noise produced by large-scale vortices(corkscrew-shaped vortices)make up the majority of low-frequency noise.Therefore,reducing pressure fluctuations is a key strategy for lowering flow-induced noise radiation.展开更多
Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, signific...Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, significantly extended design freedom and knowledge must be available in the early design stages. This paper presents a high-fidelity framework for design and optimization of the liquid swirl injectors that are widely used in aerospace propulsion and power-generation systems. The framework assembles a set of techniques, including Design Of Experiment(DOE), high-fidelity Large Eddy Simulations(LES), machine learning, Proper Orthogonal Decomposition(POD)-based Kriging surrogate modeling(emulation), inverse problem optimization, and uncertainty quantification. LES-based simulations can reveal detailed spatiotemporal evolution of flow structures and flame dynamics in a high-fidelity manner, and identify important injector design parameters according to their effects on propellant mixing, flame stabilization, and thermal protection.For a given a space of design parameters, DOE determines the number of design points to perform LES-based simulations. POD-based emulations, trained by the LES database, can effectively explore the design space and deduce an optimal group of design parameters in a turn-around time that is reduced by three orders of magnitude. The accuracy of the emulated results is validated, and the uncertainty of prediction is quantified. The proposed design methodology is expected to profoundly extend the knowledge base and reduce the cost for initial design stages.展开更多
Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering application...Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.展开更多
Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the c...Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the characteristics of inflow turbulence,can significantly impact the quality of predictions.In this study,we examined these boundary conditions within the context of the mountainous terrain around a long-span cable-stayed bridge using a wind tunnel experiment.Various sizes of computational domains and turbulent incoming wind velocities were used in large eddy simulations.The results show that when the height of the computational domain is five times greater than the height of the terrain model,there is minimal influence from the top wall on the wind field characteristics in this complex mountainous area.Expanding the length of the wake region of the computational domain has negligible effects on the wind fields.Turbulence in the inlet boundary reduces the length of the wake region on a leeward hill with a low slope,but has less impact on the mean wind velocity of steep hills.展开更多
The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in...The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.展开更多
Fractal theory provides a new strategy for equipment design.In this work,we propose a novel H-like fractal(HLF)impeller to improve the uniformity of the distribution of hydrodynamics in stirred tanks.The impellers are...Fractal theory provides a new strategy for equipment design.In this work,we propose a novel H-like fractal(HLF)impeller to improve the uniformity of the distribution of hydrodynamics in stirred tanks.The impellers are constructed by replacing two vertical blades or four legs with two or four H-like subblades by fractal iterations,respectively.Flow characteristics including velocity and turbulent kinetic energy(TKE)distributions,vortices,power number,are predicted by large eddy simulation.Compared with Rushton turbine(RT)impeller when H/T=1(or dual RTs when H/T=1.5,triple RTs when H/T=2),the HLF impeller can produce a flow field with more uniform distributions of larger velocities and TKE level.The impeller with more fractal iteration times can further improve the distribution uniformity of hydrodynamics in the case of high H/T.Power analysis shows that this is mainly due to the improved energy utilization efficiency by the fractal structure design.展开更多
In this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross fl...In this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross flow were studied through experiments and numerical simulations.The results indicate a strong dependency of the cylinder's vibration response on the fin parameters.The results indicate that the vibration response and wake structure of the cylinder are significantly influenced by the parameters of the fins.The introduction of a finned cylinder affects both its own vibration amplitude and frequency,as well as the downstream cylinder.The amplitudes of finned cylinders Ⅰ and Ⅲ are reduced by 57.8% and 59.9%,respectively,compared to the bare cylinder.This reduction helps to restrain vibration and diminishes the amplitudes of the downstream cylinder.Although finned cylinder Ⅱ slightly decreases its own vibration,it increases the amplitude of the downstream cylinder by 13.7%.The mean drag coefficient and the root mean square of the lift coefficient of the finned cylinder are higher than those of the bare cylinder when the finned cylinder is positioned upstream.Smaller pitch and larger equivalent diameter will lead to increased drag,resulting in enhanced vortex shedding in the wake,which amplifies the vibrations of the cylinder in that wake.The downstream of finned cylinder Ⅱ has the widest wake and higher vortex strength,and the dynamic load and vibration of the downstream cylinder are increased.The vortex intensity decays faster in the wake of finned cylinder Ⅲ,and the vibration of the downstream cylinder is weaker.展开更多
The cutaneous ridges on dolphin skin have long been believed to effectively reduce friction drag, thereby contributing to overall drag reduction. However, since these skin ridges are oriented perpendicular to the swim...The cutaneous ridges on dolphin skin have long been believed to effectively reduce friction drag, thereby contributing to overall drag reduction. However, since these skin ridges are oriented perpendicular to the swimming direction, they also generate additional pressure drag, raising questions about the impact of the shape-induced pressure forces on swimming. Inspired by the microvibrations observed on dolphin skin, we hypothesize that the microstructure on dolphin skin is not static but dynamically oscillates in the form of Longitudinal Micro-Ultrasonic Waves (LMUWs). To explore this, we carried out a series of Computational Fluid Dynamics (CFD) simulations based on Large Eddy Simulation (LES) model to investigate the impact of pressure drag on the total drag acting on an oscillating skin surface under realistic turbulent flow conditions. The results indicate that the dynamic skin oscillations induce a new dynamic Stokes boundary layer, which has the potential to convert pressure drag into a negative force, thereby reducing total drag under the influence of traveling LMUW excitations. Furthermore, a relative velocity ξ, defined as the difference between the wave speed c and the external flow speed U, is introduced to evaluate the drag-reduction effect dominated by pressure. The findings reveal that pressure drag remains negative when ξ > 0. As ξ increases, the thrust effect induced by negative pressure becomes increasingly significant, ultimately counteracting friction drag and eliminating total drag. This pressure-dominated drag reduction mechanism thus demonstrates a novel strategy for the drag reduction technology and the potential of unveiling the mysteries behind dolphin swimming.展开更多
Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealize...Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealized numerical experiment using the large eddy simulation method was conducted to analyze the three-dimensional flow field structure and the mechanism for dense water subsidence.Specifically,a negative salt flux is set at the sea surface,in which salt flux enters the sea surface to simulate the icing and salting-out phenomena that occur at high latitudes.Results show that the mean-state 3D flow field of dense water subsidence exhibits a hollow conical distribution.The horizontal flow field is characterized by a cyclonic vortex that driven primarily by the pressure gradient and influenced by the Coriolis effect.Moreover,the inverse vertical pressure gradient generated by this vortex inhibits the sinking of the plume,leading to its off-axis deflection and the development of an anticyclonic precession.In addition,the impact of rotation on the structure of a sinking plume within a stratified environment is discussed.Both horizontal vortex intensity and cone angle of the hollow cone flow field are increased with increasing rotation rate,resulting in a decrease in the plume’s maximum sinking depth.Variances in rotation direction cause the horizontal vortex and sinking plumes of dense water in the northern and southern hemispheres to rotate in opposite directions.展开更多
A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equ...A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.展开更多
A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation （LES）. The Idaho Falls （IDF） measurement of run 9 case is simulated numerica...A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation （LES）. The Idaho Falls （IDF） measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V^iis/il^i （BV） frequencies are com- puted with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, sec-ondary structures such as baroclinic vorticity and helical structures are also simulated. Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.展开更多
Large eddy simulations （LES） have been performed to investigate the flow and combustion fields in the scramjet of the German Aerospace Center （DLR）. Turbulent combustion is mod- eled by the tabulated thermo-chemis...Large eddy simulations （LES） have been performed to investigate the flow and combustion fields in the scramjet of the German Aerospace Center （DLR）. Turbulent combustion is mod- eled by the tabulated thermo-chemistry approach in combination with the presumed probability density function （PDF）. A/3-function is used to model the distribution of the mixture fraction, while two different PDFs, g-function （Model I） and //-function （Model II）, are applied to model the reaction progress. Temperature is obtained by solving filtered energy transport equation and the reaction rate of the progress variable is rescaled by pressure to consider the effects of compressibil- ity. The adaptive mesh refinement （AMR） technique is used to properly capture shock waves, boundary layers, shear layers and flame structures. Statistical results of temperature and velocity predicted by Model II show better accuracy than that predicted by Model I. The results of scatter points and mixture fraction-conditional variables indicate the significant differences between Model I and Model II. It is concluded that second moment information in the presumed PDF of the reaction progress is very important in the simulation of supersonic combustion. It is also found that an unstable flame with extinction and ignition develops in the shear layers of bluff body and a fuel- rich partially premixed flame stabilizes in the central recirculation bubble.展开更多
The noise induced by the fluctuant saturated steam flow under 250 °C in a stop-valve was numerically studied.The simulation was carried out using computational fluid dynamics（CFD） and ACTRAN.The acoustic field ...The noise induced by the fluctuant saturated steam flow under 250 °C in a stop-valve was numerically studied.The simulation was carried out using computational fluid dynamics（CFD） and ACTRAN.The acoustic field was investigated with Lighthill＇s acoustic analogy based on the properties of the flow field obtained using a large-eddy simulation that employs the LES-WALE dynamic model as the sub-grid-scale model.Firstly,the validation of mesh was well conducted,illustrating that two million elements were sufficient in this situation.Secondly,the treatment of the steam was deliberated,and conclusions indicate that when predicting the flow-induced noise of the stop-valve,the steam can be treated as incompressible gas at a low inlet velocity.Thirdly,the flow-induced noises under different inlet velocities were compared.The findings reveal it has remarkable influence on the flow-induced noises.Lastly,whether or not the heat preservation of the wall has influence on the noise was taken into account.The results show that heat preservation of the wall had little influence.展开更多
Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated ac...Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated accuracy has been evaluated using the LES verification and validation method.Various cavitation numbers are simulated to study important flow characteristics in the sphere wake,e.g.periodic cavity growth/contraction,interactions between the cloud and sheet cavitations and the vortex structure evolution.The spectral characteristics of the wake for typical cloud cavitation conditions were classified as the periodic cavitation mode,high Strouhal number mode and low Strouhal number mode.Main frequency distributions in the wake were analyzed and different dominant flow structures were identified for each of the three modes.Further,the cavitation and vortex relationship was also studied,which is an important issue associated with complex cavitating sphere wakes.Three types of cavitating vortex structures alternate,which indicates that three different cavity shedding regimes may exist in the wake.Analysis of vorticity transport equation shows a significant vorticity increase at the cavitation closure region and in the vortex cavitation region.This study provides a physical perspective to further understand the flow mechanisms in cavitating sphere wakes.展开更多
Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The...Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The modeling and analysis show that the LES model can simulate the planetary boundary layer (PBL) with a uniform underlying surface under various stratifications very well. Then, similar to the description of a forest canopy, the drag term on momentum and the production term of TKE by subgrid city buildings are introduced into the LES equations to account for the area-averaged effect of the subgrid urban canopy elements and to simulate the meteorological fields of the urban boundary layer (UBL). Numerical experiments and comparison analysis show that: (1) the result from the LES of the UBL with a proposed formula for the drag coefficient is consistent and comparable with that from wind tunnel experiments and an urban subdomain scale model; (2) due to the effect of urban buildings, the wind velocity near the canopy is decreased, turbulence is intensified, TKE, variance, and momentum flux are increased, the momentum and heat flux at the top of the PBL are increased, and the development of the PBL is quickened; (3) the height of the roughness sublayer (RS) of the actual city buildings is the maximum building height (1.5-3 times the mean building height), and a constant flux layer (CFL) exists in the lower part of the UBL.展开更多
In this paper,the Eulerian Stochastic Field(ESF)model in the Transported Probability Density Function(TPDF)class model is combined with the Flamelet Generated Manifolds(FGM)model.This method solves the joint probabili...In this paper,the Eulerian Stochastic Field(ESF)model in the Transported Probability Density Function(TPDF)class model is combined with the Flamelet Generated Manifolds(FGM)model.This method solves the joint probability density function transport equation by ESF method that considers the interaction mechanism between flame and turbulence with high precision.At the same time,by making use of the advantage of the FGM model,this model is able to incorporate the detailed chemical reaction mechanism(GRI 3.0)with acceptable computational cost.The new model has been implemented in the open source CFD suite-Open FOAM.Validation of the model has been carried out by simulating the Sandia flame series(three turbulent piloted methane jet flames)issued by the National Laboratory of the United States.The accuracy and advancement of the ESF/FGM turbulent combustion model are verified by comparing the LES results of the new model with the rich experimental data as well as the RANS results.The results demonstrate that the model has a strong ability in capturing combustion phenomena such as extinction and re-ignition in turbulent flame,which is essential in the accurate prediction of the combustion process in real combustion devices,for example,aircraft engines.展开更多
In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake d...In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.展开更多
A hybrid method is presented to numerically investigate the wind turbine aerodynamic characteristics.The wind turbine blade is replaced by an actuator line model.Turbulence is treated using a dynamic one-equation subg...A hybrid method is presented to numerically investigate the wind turbine aerodynamic characteristics.The wind turbine blade is replaced by an actuator line model.Turbulence is treated using a dynamic one-equation subgrid-scale model in large eddy simulation.Detailed information on the basic characteristics of the wind turbine wake is obtained and discussed.The rotor aerodynamic performance agrees well with the measurements.The actuator line method large-eddy simulation(ALM-LES)technique demonstrates its high potential in providing accurate load prediction and high resolution of turbulent fluctuations in the wind turbine wakes and the interactions within a feasible cost.展开更多
It is highly attractive to develop an efficient and flexible large eddy simulation(LES)technique for high-Reynolds-number atmospheric boundary layer(ABL)simulation using the low-order numerical scheme on a relatively ...It is highly attractive to develop an efficient and flexible large eddy simulation(LES)technique for high-Reynolds-number atmospheric boundary layer(ABL)simulation using the low-order numerical scheme on a relatively coarse grid,that could reproduce the logarithmic profile of the mean velocity and some key features of large-scale coherent structures in the outer layer.In this study,an improved near-wall correction scheme for the vertical gradient of the resolved streamwise velocity in the strain-rate tensor is proposed to calculate the eddy viscosity coefficient in the subgrid-scale(SGS)model.The LES code is realized with a second-order finite-difference scheme,the scale-dependent dynamic SGS stress model,the equilibrium wall stress model,and the proposed correction scheme.Very-high-Reynolds-number ABL flow simulation under the neutral stratification condition is conducted to assess the performance of the method in predicting the mean and fluctuating characteristics of the rough-wall turbulence.It is found that the logarithmic profile of the mean streamwise velocity and some key features of large-scale coherent structures can be reasonably predicted by adopting the proposed correction method and the low-order numerical scheme.展开更多
Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics h...Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.52179087).
文摘For Francis turbines,frequent operations under extremely low load conditions result in significant noise and pressure fluctuation issues.These issues may cause vibration and fatigue damage to the unit,accompanied by difficulties in connecting to the grid and reductions in the power generation efficiency of renewable energy.However,there is limited research on the relationship between pressure fluctuations and the induced noise of Francis turbines during extreme operations.In the present study,an acoustic numerical simulation based on the Ffowcs Williams-Hawkings equation and large eddy simulation is used to analyze the acoustic performances of Francis turbines.In the current study,for evaluating the acoustic characteristics under such terrible conditions,the results of variable flow rate and guide vane opening conditions are compared.Results indicated that Francis turbine noise is mostly due to pressure fluctuations brought on by rotor-stator interference and corkscrew-shaped vortices.The blade passing frequency(BPF)of 130.00 Hz and the low frequency of 0.33 f_(n)(where f_(n)denotes the rotating frequency)are the key factors affecting pressure and noise fluctuations.The influence of low frequency is reduced as the flow rate rises,whereas the influence of BPF gradually increases.Besides,the hydrodynamic noise of Francis turbines is primarily low-frequency,with discrete and broad-band features.The rotating noise with distinct peak values and the turbulence noise produced by large-scale vortices(corkscrew-shaped vortices)make up the majority of low-frequency noise.Therefore,reducing pressure fluctuations is a key strategy for lowering flow-induced noise radiation.
基金sponsored by the William RT Oakes Endowment of the Georgia Institute of Technology
文摘Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, significantly extended design freedom and knowledge must be available in the early design stages. This paper presents a high-fidelity framework for design and optimization of the liquid swirl injectors that are widely used in aerospace propulsion and power-generation systems. The framework assembles a set of techniques, including Design Of Experiment(DOE), high-fidelity Large Eddy Simulations(LES), machine learning, Proper Orthogonal Decomposition(POD)-based Kriging surrogate modeling(emulation), inverse problem optimization, and uncertainty quantification. LES-based simulations can reveal detailed spatiotemporal evolution of flow structures and flame dynamics in a high-fidelity manner, and identify important injector design parameters according to their effects on propellant mixing, flame stabilization, and thermal protection.For a given a space of design parameters, DOE determines the number of design points to perform LES-based simulations. POD-based emulations, trained by the LES database, can effectively explore the design space and deduce an optimal group of design parameters in a turn-around time that is reduced by three orders of magnitude. The accuracy of the emulated results is validated, and the uncertainty of prediction is quantified. The proposed design methodology is expected to profoundly extend the knowledge base and reduce the cost for initial design stages.
基金supported by the National Natural Science Foundation of China(Grant No.50979059)
文摘Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.
基金supported by the National Natural Science Foundation of China(Nos.51925808 and 52178516)the Natural Science Foundation of Hunan Province(Nos.2020JJ5745 and 2023JJ20073),China.
文摘Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the characteristics of inflow turbulence,can significantly impact the quality of predictions.In this study,we examined these boundary conditions within the context of the mountainous terrain around a long-span cable-stayed bridge using a wind tunnel experiment.Various sizes of computational domains and turbulent incoming wind velocities were used in large eddy simulations.The results show that when the height of the computational domain is five times greater than the height of the terrain model,there is minimal influence from the top wall on the wind field characteristics in this complex mountainous area.Expanding the length of the wake region of the computational domain has negligible effects on the wind fields.Turbulence in the inlet boundary reduces the length of the wake region on a leeward hill with a low slope,but has less impact on the mean wind velocity of steep hills.
基金Supported by the National Natural Science Foundation of China(10772082)~~
文摘The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.
基金the financial support from the National Natural Science Foundation of China(22078058)。
文摘Fractal theory provides a new strategy for equipment design.In this work,we propose a novel H-like fractal(HLF)impeller to improve the uniformity of the distribution of hydrodynamics in stirred tanks.The impellers are constructed by replacing two vertical blades or four legs with two or four H-like subblades by fractal iterations,respectively.Flow characteristics including velocity and turbulent kinetic energy(TKE)distributions,vortices,power number,are predicted by large eddy simulation.Compared with Rushton turbine(RT)impeller when H/T=1(or dual RTs when H/T=1.5,triple RTs when H/T=2),the HLF impeller can produce a flow field with more uniform distributions of larger velocities and TKE level.The impeller with more fractal iteration times can further improve the distribution uniformity of hydrodynamics in the case of high H/T.Power analysis shows that this is mainly due to the improved energy utilization efficiency by the fractal structure design.
基金financially supported by the National Natural Science Foundation of China(22478286)。
文摘In this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross flow were studied through experiments and numerical simulations.The results indicate a strong dependency of the cylinder's vibration response on the fin parameters.The results indicate that the vibration response and wake structure of the cylinder are significantly influenced by the parameters of the fins.The introduction of a finned cylinder affects both its own vibration amplitude and frequency,as well as the downstream cylinder.The amplitudes of finned cylinders Ⅰ and Ⅲ are reduced by 57.8% and 59.9%,respectively,compared to the bare cylinder.This reduction helps to restrain vibration and diminishes the amplitudes of the downstream cylinder.Although finned cylinder Ⅱ slightly decreases its own vibration,it increases the amplitude of the downstream cylinder by 13.7%.The mean drag coefficient and the root mean square of the lift coefficient of the finned cylinder are higher than those of the bare cylinder when the finned cylinder is positioned upstream.Smaller pitch and larger equivalent diameter will lead to increased drag,resulting in enhanced vortex shedding in the wake,which amplifies the vibrations of the cylinder in that wake.The downstream of finned cylinder Ⅱ has the widest wake and higher vortex strength,and the dynamic load and vibration of the downstream cylinder are increased.The vortex intensity decays faster in the wake of finned cylinder Ⅲ,and the vibration of the downstream cylinder is weaker.
基金supported by the Japan Society for the Promotion of Science(JSPS KAKENHI No.23H01373).
文摘The cutaneous ridges on dolphin skin have long been believed to effectively reduce friction drag, thereby contributing to overall drag reduction. However, since these skin ridges are oriented perpendicular to the swimming direction, they also generate additional pressure drag, raising questions about the impact of the shape-induced pressure forces on swimming. Inspired by the microvibrations observed on dolphin skin, we hypothesize that the microstructure on dolphin skin is not static but dynamically oscillates in the form of Longitudinal Micro-Ultrasonic Waves (LMUWs). To explore this, we carried out a series of Computational Fluid Dynamics (CFD) simulations based on Large Eddy Simulation (LES) model to investigate the impact of pressure drag on the total drag acting on an oscillating skin surface under realistic turbulent flow conditions. The results indicate that the dynamic skin oscillations induce a new dynamic Stokes boundary layer, which has the potential to convert pressure drag into a negative force, thereby reducing total drag under the influence of traveling LMUW excitations. Furthermore, a relative velocity ξ, defined as the difference between the wave speed c and the external flow speed U, is introduced to evaluate the drag-reduction effect dominated by pressure. The findings reveal that pressure drag remains negative when ξ > 0. As ξ increases, the thrust effect induced by negative pressure becomes increasingly significant, ultimately counteracting friction drag and eliminating total drag. This pressure-dominated drag reduction mechanism thus demonstrates a novel strategy for the drag reduction technology and the potential of unveiling the mysteries behind dolphin swimming.
基金Supported by the National Natural Science Foundation of China(Nos.42250710152,42192562)the Southern Laboratory of Ocean Science and Engineering(Guangdong Zhuhai)(No.SML 2020 SP 007)。
文摘Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealized numerical experiment using the large eddy simulation method was conducted to analyze the three-dimensional flow field structure and the mechanism for dense water subsidence.Specifically,a negative salt flux is set at the sea surface,in which salt flux enters the sea surface to simulate the icing and salting-out phenomena that occur at high latitudes.Results show that the mean-state 3D flow field of dense water subsidence exhibits a hollow conical distribution.The horizontal flow field is characterized by a cyclonic vortex that driven primarily by the pressure gradient and influenced by the Coriolis effect.Moreover,the inverse vertical pressure gradient generated by this vortex inhibits the sinking of the plume,leading to its off-axis deflection and the development of an anticyclonic precession.In addition,the impact of rotation on the structure of a sinking plume within a stratified environment is discussed.Both horizontal vortex intensity and cone angle of the hollow cone flow field are increased with increasing rotation rate,resulting in a decrease in the plume’s maximum sinking depth.Variances in rotation direction cause the horizontal vortex and sinking plumes of dense water in the northern and southern hemispheres to rotate in opposite directions.
基金Project supported by the National Natural Science Foundation of China(Nos.91752118,11672305,11232011,and 11572331)the Strategic Priority Research Program(No.XDB22040104)the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences(No.QYZDJ-SSWSYS002)
文摘A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.
基金Project supported by the Boeing-COMAC Aviation Energy Conservation and Emissions Reduction Technology Center(AECER)
文摘A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation （LES）. The Idaho Falls （IDF） measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V^iis/il^i （BV） frequencies are com- puted with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, sec-ondary structures such as baroclinic vorticity and helical structures are also simulated. Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.
基金financial support by the National Natural Science Foundation of China (Nos. 51176178, 50936005)
文摘Large eddy simulations （LES） have been performed to investigate the flow and combustion fields in the scramjet of the German Aerospace Center （DLR）. Turbulent combustion is mod- eled by the tabulated thermo-chemistry approach in combination with the presumed probability density function （PDF）. A/3-function is used to model the distribution of the mixture fraction, while two different PDFs, g-function （Model I） and //-function （Model II）, are applied to model the reaction progress. Temperature is obtained by solving filtered energy transport equation and the reaction rate of the progress variable is rescaled by pressure to consider the effects of compressibil- ity. The adaptive mesh refinement （AMR） technique is used to properly capture shock waves, boundary layers, shear layers and flame structures. Statistical results of temperature and velocity predicted by Model II show better accuracy than that predicted by Model I. The results of scatter points and mixture fraction-conditional variables indicate the significant differences between Model I and Model II. It is concluded that second moment information in the presumed PDF of the reaction progress is very important in the simulation of supersonic combustion. It is also found that an unstable flame with extinction and ignition develops in the shear layers of bluff body and a fuel- rich partially premixed flame stabilizes in the central recirculation bubble.
文摘The noise induced by the fluctuant saturated steam flow under 250 °C in a stop-valve was numerically studied.The simulation was carried out using computational fluid dynamics（CFD） and ACTRAN.The acoustic field was investigated with Lighthill＇s acoustic analogy based on the properties of the flow field obtained using a large-eddy simulation that employs the LES-WALE dynamic model as the sub-grid-scale model.Firstly,the validation of mesh was well conducted,illustrating that two million elements were sufficient in this situation.Secondly,the treatment of the steam was deliberated,and conclusions indicate that when predicting the flow-induced noise of the stop-valve,the steam can be treated as incompressible gas at a low inlet velocity.Thirdly,the flow-induced noises under different inlet velocities were compared.The findings reveal it has remarkable influence on the flow-induced noises.Lastly,whether or not the heat preservation of the wall has influence on the noise was taken into account.The results show that heat preservation of the wall had little influence.
基金This work was financially supported by the National Natural Science Foundation of China(Grants 51822903 and 11772239)the Natural Science Foundation of Hubei Province(Grant 2018CFA010).The numerical calculations were done on the supercomputing system in the Supercomputing Center of Wuhan University.
文摘Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated accuracy has been evaluated using the LES verification and validation method.Various cavitation numbers are simulated to study important flow characteristics in the sphere wake,e.g.periodic cavity growth/contraction,interactions between the cloud and sheet cavitations and the vortex structure evolution.The spectral characteristics of the wake for typical cloud cavitation conditions were classified as the periodic cavitation mode,high Strouhal number mode and low Strouhal number mode.Main frequency distributions in the wake were analyzed and different dominant flow structures were identified for each of the three modes.Further,the cavitation and vortex relationship was also studied,which is an important issue associated with complex cavitating sphere wakes.Three types of cavitating vortex structures alternate,which indicates that three different cavity shedding regimes may exist in the wake.Analysis of vorticity transport equation shows a significant vorticity increase at the cavitation closure region and in the vortex cavitation region.This study provides a physical perspective to further understand the flow mechanisms in cavitating sphere wakes.
基金The research was supported by the National Natural Science Foundation of China under Grant Nos.40333027 and 40075004.
文摘Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The modeling and analysis show that the LES model can simulate the planetary boundary layer (PBL) with a uniform underlying surface under various stratifications very well. Then, similar to the description of a forest canopy, the drag term on momentum and the production term of TKE by subgrid city buildings are introduced into the LES equations to account for the area-averaged effect of the subgrid urban canopy elements and to simulate the meteorological fields of the urban boundary layer (UBL). Numerical experiments and comparison analysis show that: (1) the result from the LES of the UBL with a proposed formula for the drag coefficient is consistent and comparable with that from wind tunnel experiments and an urban subdomain scale model; (2) due to the effect of urban buildings, the wind velocity near the canopy is decreased, turbulence is intensified, TKE, variance, and momentum flux are increased, the momentum and heat flux at the top of the PBL are increased, and the development of the PBL is quickened; (3) the height of the roughness sublayer (RS) of the actual city buildings is the maximum building height (1.5-3 times the mean building height), and a constant flux layer (CFL) exists in the lower part of the UBL.
基金supported by the National Natural Science Foundation of China(No.51706241).
文摘In this paper,the Eulerian Stochastic Field(ESF)model in the Transported Probability Density Function(TPDF)class model is combined with the Flamelet Generated Manifolds(FGM)model.This method solves the joint probability density function transport equation by ESF method that considers the interaction mechanism between flame and turbulence with high precision.At the same time,by making use of the advantage of the FGM model,this model is able to incorporate the detailed chemical reaction mechanism(GRI 3.0)with acceptable computational cost.The new model has been implemented in the open source CFD suite-Open FOAM.Validation of the model has been carried out by simulating the Sandia flame series(three turbulent piloted methane jet flames)issued by the National Laboratory of the United States.The accuracy and advancement of the ESF/FGM turbulent combustion model are verified by comparing the LES results of the new model with the rich experimental data as well as the RANS results.The results demonstrate that the model has a strong ability in capturing combustion phenomena such as extinction and re-ignition in turbulent flame,which is essential in the accurate prediction of the combustion process in real combustion devices,for example,aircraft engines.
基金Project supported by the National Key Research and Development Program of China(Nos.2019YFE0192600,2017YFE0132000,and 2019YFB1503700)the National Natural Science Foundation of China(Nos.51761135012 and 11872248)。
文摘In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.
基金funded jointly by the National Basic Research Program of China(″973″Program)(No.2014CB046200)the Jiangsu Provincial Natural Science Foundation(No.BK20140059)+2 种基金the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe National Natural Science Foundation of China(No.11172135)the EU Seventh Framework Program(No.FP7-PEOPLE-2010-IRSES-269202)
文摘A hybrid method is presented to numerically investigate the wind turbine aerodynamic characteristics.The wind turbine blade is replaced by an actuator line model.Turbulence is treated using a dynamic one-equation subgrid-scale model in large eddy simulation.Detailed information on the basic characteristics of the wind turbine wake is obtained and discussed.The rotor aerodynamic performance agrees well with the measurements.The actuator line method large-eddy simulation(ALM-LES)technique demonstrates its high potential in providing accurate load prediction and high resolution of turbulent fluctuations in the wind turbine wakes and the interactions within a feasible cost.
基金Project supported by the National Natural Science Foundation of China(No.11490551)the Fundamental Research Funds for the Central Universities(Nos.lzujbky-2016-k13 and lzujbky-2018-k07)
文摘It is highly attractive to develop an efficient and flexible large eddy simulation(LES)technique for high-Reynolds-number atmospheric boundary layer(ABL)simulation using the low-order numerical scheme on a relatively coarse grid,that could reproduce the logarithmic profile of the mean velocity and some key features of large-scale coherent structures in the outer layer.In this study,an improved near-wall correction scheme for the vertical gradient of the resolved streamwise velocity in the strain-rate tensor is proposed to calculate the eddy viscosity coefficient in the subgrid-scale(SGS)model.The LES code is realized with a second-order finite-difference scheme,the scale-dependent dynamic SGS stress model,the equilibrium wall stress model,and the proposed correction scheme.Very-high-Reynolds-number ABL flow simulation under the neutral stratification condition is conducted to assess the performance of the method in predicting the mean and fluctuating characteristics of the rough-wall turbulence.It is found that the logarithmic profile of the mean streamwise velocity and some key features of large-scale coherent structures can be reasonably predicted by adopting the proposed correction method and the low-order numerical scheme.
文摘Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.
