The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameteri...The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.展开更多
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.展开更多
The current research of large eddy simulation (LES) of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale (SGS) model to simulate turbulent flow, which aims at obtaining the flow fi...The current research of large eddy simulation (LES) of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale (SGS) model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry (PIV) and laser Doppler velocimetry (LDV). The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.展开更多
Eddy-damping quasinormal Markovian (EDQNM) theory is employed to calculate the resolved-scale spectrum and transfer spectrum, based on which we investigate the resolved-scale scaling law. Results show that the scali...Eddy-damping quasinormal Markovian (EDQNM) theory is employed to calculate the resolved-scale spectrum and transfer spectrum, based on which we investigate the resolved-scale scaling law. Results show that the scaling law of the resolved-scale turbulence, which is affected by several factors, is far from that of the full-scale turbulence and should be corrected. These results are then applied to an existing subgrid model to improve its performance. A series of simulations are performed to verify the necessity of a fixed scaling law in the subgrid modeling.展开更多
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.展开更多
The entrainment flux ratio Ae and the inversion layer (IL) thickness are two key parameters in a mixed layer model. Ae is defined as the ratio of the entrainment heat flux at the mixed layer top to the surface heat ...The entrainment flux ratio Ae and the inversion layer (IL) thickness are two key parameters in a mixed layer model. Ae is defined as the ratio of the entrainment heat flux at the mixed layer top to the surface heat flux. The IL is the layer between the mixed layer and the free atmosphere. In this study, a parameterization of Ae is derived from the TKE budget in the first- order model for a well-developed CBL under the condition of linearly sheared geostrophic velocity with a zero value at the surface. It is also appropriate for a CBL under the condition of geostrophic velocity remaining constant with height. LESs are conducted under the above two conditions to determine the coefficients in the parameterization scheme. Results suggest that about 43% of the shear-produced TKE in the IL is available for entrainment, while the shear-produced TKE in the mixed layer and surface layer have little effect on entrainment. Based on this scheme, a new scale of convective turbulence velocity is proposed and applied to parameterize the IL thickness, The LES outputs for the CBLs under the condition of linearly sheared geostrophic velocity with a non-zero surface value are used to verify the performance of the parameterization scheme. It is found that the parameterized Ae and IL thickness agree well with the LES outputs.展开更多
The flows past a number 3900 are simulated circular cylinder at Reynolds using large-eddy simulation (LES) and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volum...The flows past a number 3900 are simulated circular cylinder at Reynolds using large-eddy simulation (LES) and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volume scheme is used to discretize the incompressible Navier- Stokes equations. The dynamic global coefficient version of the Vreman's subgrid scale (SGS) model is used to com- pute the sub-grid stresses. Curie's integral of Lighthill's acoustic analogy is used to extract the sound radiated from the cylinder. The profiles of mean velocity and turbulent fluctua- tions obtained are consistent with the previous experimental and computational results. The sound radiation at far field exhibits the characteristic of a dipole and directivity. The sound spectra display the -5/3 power law. It is shown that Vreman's SGS model in company with dynamic procedure is suitable for LES of turbulence generated noise.展开更多
The so-called T-shaped reducing tees are typically used to divide,change and control(to a certain extent)the flow direction in pipe networks.In this study,the Ffowcs Williams–Hawkings(FW-H)equation and the Large Eddy...The so-called T-shaped reducing tees are typically used to divide,change and control(to a certain extent)the flow direction in pipe networks.In this study,the Ffowcs Williams–Hawkings(FW-H)equation and the Large Eddy Simulation(LES)methods are used to simulate the flow-induced noise related to T-shaped reducing tees under different inlet flow velocities and for different pipe diameter ratios.The results show that the maximum flow velocity,average flow velocity,and vorticity in the branch pipe increase gradually as the related diameter decreases.Strong vorticity and secondary flows are also observed in the branch pipe,and the associated violent pressure fluctuations are found to be the main sources of flow-induced noise.In particular,as the pipe diameter ratio decreases from 1 to 0.45,the Total Sound Pressure Level(TSPL)increases by 6.8,6.26,and 7.43 dB for values of the inlet flow velocity of 1,2,and 3 m/s,respectively.The distribution characteristics of the flow-induced noise in the frequency domain follow similar trends for different pipe diameter ratios.展开更多
The temperature fluctuation caused by thermal striping phenomena of hot and cold fluids mixing results in cyclical thermal stress fatigue failure of the pipe wall. Mean temperature difference between hot and cold flui...The temperature fluctuation caused by thermal striping phenomena of hot and cold fluids mixing results in cyclical thermal stress fatigue failure of the pipe wall. Mean temperature difference between hot and cold fluids was often used as thermal load in previous analysis of thermal fatigue failure, thereby the influences of the amplitude and frequency of temperature fluctuation on thermal fatigue failure were neglected. Based on the mechanism of flow and heat transfer which induces thermal fatigue, the turbulent mixing of hot and cold water in a tee junction is simulated with FLUENT platform by using the Large-eddy simulation(LES) turbulent flow model with the sub-grid scale(SGS) model of Smagorinsky-Lilly(SL) to capture the amplitude and frequency of temperature fluctuation. In a simulation case, hot water with temperature of 343.48 K and velocity of 0.15 m/s enters the horizontal main duct with the side length of 100 mm, while cold water with temperature of 296.78 K and velocity of 0.3 m/s enters the vertical branch duct with the side length of 50 mm. The numerical results show that the mean and fluctuating temperatures are in good agreement with the previous experimental data, which describes numerical simulation with high reliability and accuracy; the power spectrum density(PSD) on top wall is higher than that on bottom wall(as the frequency less than 1 Hz), while the PSD on bottom wall is relatively higher than that on top wall (as the frequency of 1-10Hz). The temperature fluctuations in full mixing region of the tee junction can be accurately captured by LES and can provide the theoretical basis for the thermal stress and thermal fatigue analyses.展开更多
Mixing processes of hot and cold fluids in a tee with and without sin- tered copper spheres are simulated by FLUENT using the large-eddy simulation (LES) turbulent flow model and the sub-grid scale (SGS) Smagorins...Mixing processes of hot and cold fluids in a tee with and without sin- tered copper spheres are simulated by FLUENT using the large-eddy simulation (LES) turbulent flow model and the sub-grid scale (SGS) Smagorinsky-Lilly (SL) model with buoyancy. Comparisons of numerical results of the two cases with and without sintered copper spheres show that the porous medium significantly reduces velocity and temper- ature fluctuations because the porous medium can effectively restrict the fluid flow and enhance heat transfer. The porous medium obviously increases the pressure drop in the main duct. The porous medium reduces the power spectrum density (PSD) of tempera- ture fluctuations in the frequency range from 1 Hz to 10 Hz.展开更多
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.展开更多
The aerodynamic performances and flow features of the capsule/rigid disk-gap-band(DGB)parachute system from the Mach number 1.8 to 2.2 are studied.We use the adaptive mesh refinement(AMR),the hybrid tuned center-diffe...The aerodynamic performances and flow features of the capsule/rigid disk-gap-band(DGB)parachute system from the Mach number 1.8 to 2.2 are studied.We use the adaptive mesh refinement(AMR),the hybrid tuned center-difference and weighted essentially non-oscillatory(TCD-WENO)scheme,and the large-eddy simulation(LES)with the stretched-vortex subgrid model.The simulations reproduce complex interaction of the flow structures,including turbulent wakes and bow shocks,as well as bow shocks and expansion waves.The results show that the calculated aerodynamic drag coefficient agrees well with the previou simulation.Both the aerodynamic drag coefficient and the aerodynamic drag oscillation of the parachute system decrease with the increase of the initial Mach number of the fluid.It is found that the position and angle of the bow shock ahead of the canopy change as the Mach number increases,which makes the flow inside the canopy and the turbulent wake behind the canopy more complex and unstable.展开更多
We analyze the error of large-eddy simulation(LES)in wall pressure fluctuation of a turbulent channel flow.To separate different sources of the error,we conduct both direct numerical simulations(DNS)and LES,and apply ...We analyze the error of large-eddy simulation(LES)in wall pressure fluctuation of a turbulent channel flow.To separate different sources of the error,we conduct both direct numerical simulations(DNS)and LES,and apply an explicit filter on DNS data to obtain filtered DNS(FDNS)data.The error of LES is consequently decomposed into two parts:The first part is the error of FDNS with respect to DNS,which quantifies the influence of the filter operation.The second part is the difference between LES and FDNS induced by the error of LES in velocity field.By comparing the root-mean-square value and the wavenumber-frequency spectrum of the wall pressure fluctuation,it is found that the inaccuracy of the velocity fluctuations is the dominant source that induces the error of LES in the wall pressure fluctuation.The present study provides a basis on future LES studies of the wall pressure fluctuation.展开更多
Turbulence is a phenomenon which presents peculiarities when it is experimented or simulated.This occurs due to its complexity and high sensibility to the inlet conditions of the turbulent flow fields,as well as the p...Turbulence is a phenomenon which presents peculiarities when it is experimented or simulated.This occurs due to its complexity and high sensibility to the inlet conditions of the turbulent flow fields,as well as the presence of a large range of time and length scales.A simplification for this situation is obtained with the use of approximations and turbulence models.In the present work,the Largeeddy Simulations methodology was applied,aiming the modeling of the previously mentioned complexity,which consists in using a filter to resolve the large scales while the remaining scales were determined by classical and dynamic Smagorinsky models.Three different approximations for the inlet conditions were applied:white noise,Random Flow Generation(RFG)and Synthetic Eddy Method(SEM).It was possible to realize that the use of the dynamic Smagorinky model and the RFG or SEM methodologies resulted in a better characterization of the studied flow.展开更多
A novel numerical framework is developed for large-eddy simulation(LES) of interactions among air, water, and solid bodies. The motions of air and water are solved on a fixed block-structured mesh, with the air–water...A novel numerical framework is developed for large-eddy simulation(LES) of interactions among air, water, and solid bodies. The motions of air and water are solved on a fixed block-structured mesh, with the air–water interface captured using the volume-of-fluid method. A new sub-grid scale stress model based on the vortex identifier is used to improve the robustness and efficiency of the simulation flows with air–water interface. The new framework is tested in the context of bow waves and Kelvin waves generated by a water-surface vehicle. Wave breaking at the bow of the vehicle is captured in LES. The LES results of wave geometry approaches the measurements progressively as the grid resolution is refined. The simulation results indicate that LES is a useful tool for studying wave dynamics of water-surface vehicles.展开更多
Large-eddy simulations (LES) based on the temporal approximate deconvolution model were performed for a forced homogeneous isotropic turbulence (FHIT) with polymer additives at moderate Taylor Reynolds number. Fin...Large-eddy simulations (LES) based on the temporal approximate deconvolution model were performed for a forced homogeneous isotropic turbulence (FHIT) with polymer additives at moderate Taylor Reynolds number. Finitely extensible nonlinear elastic in the Peterlin approximation model was adopted as the constitutive equation for the filtered conformation tensor of the polymer molecules. The LES results were verified through comparisons with the direct numerical simulation results. Using the LES database of the FHIT in the Newtonian fluid and the polymer solution flows, the polymer effects on some important parameters such as strain, vorticity, drag reduction, and so forth were studied. By extracting the vortex structures and exploring the flatness factor through a high-order correlation function of velocity derivative and wavelet analysis, it can be found that the small-scale vortex structures and small-scale intermittency in the FHIT are all inhibited due to the existence of the polymers. The extended self-similarity scaling law in the polymer solution flow shows no apparent difference from that in the Newtonian fluid flow at the currently simulated ranges of Reynolds and Weissenberg numbers.展开更多
The near-field dynamics of an aircraft wake is studied by means of temporal large-eddy simulations, with and without considering the effect of engine jets. In the absence of jets, the simulations showed the roll-up of...The near-field dynamics of an aircraft wake is studied by means of temporal large-eddy simulations, with and without considering the effect of engine jets. In the absence of jets, the simulations showed the roll-up of the initial vorticity sheet shed by the wing and the occurrence of short-wavelength instability in a pair of primary co-rotating vortices. The main consequence of the instability is the modification of the internal structure of the vortex, compared to the two-dimensional stable behavior. The presence of engine jets affects the roll-up of the vorticity sheet and causes an enlargement of the final merged vortex core compared to the case without jets.展开更多
The purpose of the present study was to establish a passive flow control method for a rectangular jet using two types of deflectors installed symmetrically inside a nozzle. This deflector in a rectangular nozzle gener...The purpose of the present study was to establish a passive flow control method for a rectangular jet using two types of deflectors installed symmetrically inside a nozzle. This deflector in a rectangular nozzle generates the rectangular coaxial jets. The effect of the slant angle of the deflectors on the flow characteristics and the spread of the rectangular jet was investigated experimentally and by large-eddy simulation. The experiment and the numerical simulation were carried out at a Reynolds number of 9000. The rectangular jet with no deflectors generates a vortex ring from the nozzle exit. The vortex ring collapses in the downstream region and the outline of the jet changes from rectangular to diamond-shaped as a result of the axis-switching phenomenon. The rectangular jet with divergent and convergent deflectors shows particularly noticeable changes in the flow characteristics and vortical structures, as compared to the case with no deflectors. In the case of the rectangular jet with divergent deflectors (slant angle of <i style="font-family:" font-size:13.3333px;white-space:normal;"="">α</i><span style="font-family:" font-size:13.3333px;white-space:normal;"=""> = 6<span style="white-space:nowrap;">°</span>), minor axis spread is promoted more than major axis spread, and axis switching occurs closer to the nozzle exit than that in the case of no deflectors. The outline of the jet also changes from lateral rectangular to vertical rectangular as a result of axis switching. On the other hand, in the case of a rectangular jet with convergent deflectors (</span><i style="font-family:" font-size:13.3333px;white-space:normal;"="">α</i><span style="font-family:" font-size:13.3333px;white-space:normal;"=""> = -6<span style="white-space:nowrap;">°</span>), minor axis spread is suppressed more than major axis spread, and axis switching occurs farther from the nozzle exit than that in the case with no deflectors. The outline of the jet does not change until the downstream region. Therefore, the spread and the axis-switching location for the rectangular jet can be controlled by the deflectors inside the rectangular nozzle.</span>展开更多
In this study, the development of a convective boundary layer (CBL) in the Badaujilin region was investigated by comparing the observation data of two cases. A deep neutral layer capped a CBL that occurred on 30 Aug...In this study, the development of a convective boundary layer (CBL) in the Badaujilin region was investigated by comparing the observation data of two cases. A deep neutral layer capped a CBL that occurred on 30 August 2009. This case was divided into five sublayers from the surface to higher atmospheric elevations: surface layer, mixed layer, inversion layer, neutral layer, and sub-inversion layer. The development process of the CBL was divided into three stages: S1, S2, and S3. This case was quite different from the development of the three-layer CBL observed on 31 August 2009 because the mixed layer of the five-layer CBL (CBL5) eroded the neutral layer during S2. The specific initial structure of the CBL5 was correlated to the synoptic background of atmosphere during nighttime. The three-stage development process of the CBL5 was confirmed by six simulations using National Center for Atmospheric Research (USA) large-eddy simulation (NCAR-LES), and some of its characteristics are presented in detail.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.42175099,42027804,42075073)the Innovative Project of Postgraduates in Jiangsu Province in 2023(Grant No.KYCX23_1319)+3 种基金supported by the National Natural Science Foundation of China(Grant No.42205080)the Natural Science Foundation of Sichuan(Grant No.2023YFS0442)the Research Fund of Civil Aviation Flight University of China(Grant No.J2022-037)supported by the National Key Scientific and Technological Infrastructure project“Earth System Science Numerical Simulator Facility”(Earth Lab)。
文摘The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.
基金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.
基金supported by National Natural Science Foundation of China (Grant Nos. 51139007, 51079151, 51079152)Research Fundfor the Doctoral Program of Higher Education of China (Grant No. 0100008110012)
文摘The current research of large eddy simulation (LES) of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale (SGS) model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry (PIV) and laser Doppler velocimetry (LDV). The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.
基金supported by the National Natural Science Foundation of China(11202013 and 51136003)the National Basic Research Program of China(2012CB720200)the Opening fundof State Key Laboratory of Nonlinear Mechanics
文摘Eddy-damping quasinormal Markovian (EDQNM) theory is employed to calculate the resolved-scale spectrum and transfer spectrum, based on which we investigate the resolved-scale scaling law. Results show that the scaling law of the resolved-scale turbulence, which is affected by several factors, is far from that of the full-scale turbulence and should be corrected. These results are then applied to an existing subgrid model to improve its performance. A series of simulations are performed to verify the necessity of a fixed scaling law in the subgrid modeling.
基金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.
基金sponsored by the National Natural Science Foundation of China(Grant No.40975004)the State Key Basic Program(973)Program(Grant No.2013CB430100)
文摘The entrainment flux ratio Ae and the inversion layer (IL) thickness are two key parameters in a mixed layer model. Ae is defined as the ratio of the entrainment heat flux at the mixed layer top to the surface heat flux. The IL is the layer between the mixed layer and the free atmosphere. In this study, a parameterization of Ae is derived from the TKE budget in the first- order model for a well-developed CBL under the condition of linearly sheared geostrophic velocity with a zero value at the surface. It is also appropriate for a CBL under the condition of geostrophic velocity remaining constant with height. LESs are conducted under the above two conditions to determine the coefficients in the parameterization scheme. Results suggest that about 43% of the shear-produced TKE in the IL is available for entrainment, while the shear-produced TKE in the mixed layer and surface layer have little effect on entrainment. Based on this scheme, a new scale of convective turbulence velocity is proposed and applied to parameterize the IL thickness, The LES outputs for the CBLs under the condition of linearly sheared geostrophic velocity with a non-zero surface value are used to verify the performance of the parameterization scheme. It is found that the parameterized Ae and IL thickness agree well with the LES outputs.
基金supported by the National Natural Science Foundation of China(Grant 11232011)
文摘The flows past a number 3900 are simulated circular cylinder at Reynolds using large-eddy simulation (LES) and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volume scheme is used to discretize the incompressible Navier- Stokes equations. The dynamic global coefficient version of the Vreman's subgrid scale (SGS) model is used to com- pute the sub-grid stresses. Curie's integral of Lighthill's acoustic analogy is used to extract the sound radiated from the cylinder. The profiles of mean velocity and turbulent fluctua- tions obtained are consistent with the previous experimental and computational results. The sound radiation at far field exhibits the characteristic of a dipole and directivity. The sound spectra display the -5/3 power law. It is shown that Vreman's SGS model in company with dynamic procedure is suitable for LES of turbulence generated noise.
基金supported by the Shandong Engineering Laboratory for High-Efficiency Energy Conservation and Energy Storage Technology&Equipment.
文摘The so-called T-shaped reducing tees are typically used to divide,change and control(to a certain extent)the flow direction in pipe networks.In this study,the Ffowcs Williams–Hawkings(FW-H)equation and the Large Eddy Simulation(LES)methods are used to simulate the flow-induced noise related to T-shaped reducing tees under different inlet flow velocities and for different pipe diameter ratios.The results show that the maximum flow velocity,average flow velocity,and vorticity in the branch pipe increase gradually as the related diameter decreases.Strong vorticity and secondary flows are also observed in the branch pipe,and the associated violent pressure fluctuations are found to be the main sources of flow-induced noise.In particular,as the pipe diameter ratio decreases from 1 to 0.45,the Total Sound Pressure Level(TSPL)increases by 6.8,6.26,and 7.43 dB for values of the inlet flow velocity of 1,2,and 3 m/s,respectively.The distribution characteristics of the flow-induced noise in the frequency domain follow similar trends for different pipe diameter ratios.
基金supported by National Natural Science Foundation of China(Grant No. 50906002)National Basic Research Program of China(973 Program, Grant No. 2011CB706900)Beijing Novel Program of China(Grant No. 2008B16)
文摘The temperature fluctuation caused by thermal striping phenomena of hot and cold fluids mixing results in cyclical thermal stress fatigue failure of the pipe wall. Mean temperature difference between hot and cold fluids was often used as thermal load in previous analysis of thermal fatigue failure, thereby the influences of the amplitude and frequency of temperature fluctuation on thermal fatigue failure were neglected. Based on the mechanism of flow and heat transfer which induces thermal fatigue, the turbulent mixing of hot and cold water in a tee junction is simulated with FLUENT platform by using the Large-eddy simulation(LES) turbulent flow model with the sub-grid scale(SGS) model of Smagorinsky-Lilly(SL) to capture the amplitude and frequency of temperature fluctuation. In a simulation case, hot water with temperature of 343.48 K and velocity of 0.15 m/s enters the horizontal main duct with the side length of 100 mm, while cold water with temperature of 296.78 K and velocity of 0.3 m/s enters the vertical branch duct with the side length of 50 mm. The numerical results show that the mean and fluctuating temperatures are in good agreement with the previous experimental data, which describes numerical simulation with high reliability and accuracy; the power spectrum density(PSD) on top wall is higher than that on bottom wall(as the frequency less than 1 Hz), while the PSD on bottom wall is relatively higher than that on top wall (as the frequency of 1-10Hz). The temperature fluctuations in full mixing region of the tee junction can be accurately captured by LES and can provide the theoretical basis for the thermal stress and thermal fatigue analyses.
基金supported by the National Natural Science Foundation of China (No. 50906002)the National Basic Research Program of China (No. 2011CB706900)+1 种基金the Research Fund for the Doctoral Program of Higher Education of China (No. 20090010110006)the Beijing Novel Program of China (No. 2008B16)
文摘Mixing processes of hot and cold fluids in a tee with and without sin- tered copper spheres are simulated by FLUENT using the large-eddy simulation (LES) turbulent flow model and the sub-grid scale (SGS) Smagorinsky-Lilly (SL) model with buoyancy. Comparisons of numerical results of the two cases with and without sintered copper spheres show that the porous medium significantly reduces velocity and temper- ature fluctuations because the porous medium can effectively restrict the fluid flow and enhance heat transfer. The porous medium obviously increases the pressure drop in the main duct. The porous medium reduces the power spectrum density (PSD) of tempera- ture fluctuations in the frequency range from 1 Hz to 10 Hz.
基金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.
基金Project supported by the National Natural Science Foundation of China(No.11372068)the National Basic Research Program of China(973 Program)(No.2014CB744104)。
文摘The aerodynamic performances and flow features of the capsule/rigid disk-gap-band(DGB)parachute system from the Mach number 1.8 to 2.2 are studied.We use the adaptive mesh refinement(AMR),the hybrid tuned center-difference and weighted essentially non-oscillatory(TCD-WENO)scheme,and the large-eddy simulation(LES)with the stretched-vortex subgrid model.The simulations reproduce complex interaction of the flow structures,including turbulent wakes and bow shocks,as well as bow shocks and expansion waves.The results show that the calculated aerodynamic drag coefficient agrees well with the previou simulation.Both the aerodynamic drag coefficient and the aerodynamic drag oscillation of the parachute system decrease with the increase of the initial Mach number of the fluid.It is found that the position and angle of the bow shock ahead of the canopy change as the Mach number increases,which makes the flow inside the canopy and the turbulent wake behind the canopy more complex and unstable.
基金This research is supported by the National Natural Science Foundation of China(NFSC)Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”(Grant 11988102)the National Key Project(Grant GJXM92579)Shizhao Wang acknowledges the support from the National Natural Science Foundation of China(Grant 11922214).
文摘We analyze the error of large-eddy simulation(LES)in wall pressure fluctuation of a turbulent channel flow.To separate different sources of the error,we conduct both direct numerical simulations(DNS)and LES,and apply an explicit filter on DNS data to obtain filtered DNS(FDNS)data.The error of LES is consequently decomposed into two parts:The first part is the error of FDNS with respect to DNS,which quantifies the influence of the filter operation.The second part is the difference between LES and FDNS induced by the error of LES in velocity field.By comparing the root-mean-square value and the wavenumber-frequency spectrum of the wall pressure fluctuation,it is found that the inaccuracy of the velocity fluctuations is the dominant source that induces the error of LES in the wall pressure fluctuation.The present study provides a basis on future LES studies of the wall pressure fluctuation.
文摘Turbulence is a phenomenon which presents peculiarities when it is experimented or simulated.This occurs due to its complexity and high sensibility to the inlet conditions of the turbulent flow fields,as well as the presence of a large range of time and length scales.A simplification for this situation is obtained with the use of approximations and turbulence models.In the present work,the Largeeddy Simulations methodology was applied,aiming the modeling of the previously mentioned complexity,which consists in using a filter to resolve the large scales while the remaining scales were determined by classical and dynamic Smagorinsky models.Three different approximations for the inlet conditions were applied:white noise,Random Flow Generation(RFG)and Synthetic Eddy Method(SEM).It was possible to realize that the use of the dynamic Smagorinky model and the RFG or SEM methodologies resulted in a better characterization of the studied flow.
基金supported by Lixing plan of Institute of Mechanics,Chinese Academy of Sciencesthe National Natural Science Foundation of China(91752119,11232011 and11572331)+1 种基金the support from the Strategic Priority Research Program(XDB22040104)the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences(QYZDJ-SSW-SYS002)
文摘A novel numerical framework is developed for large-eddy simulation(LES) of interactions among air, water, and solid bodies. The motions of air and water are solved on a fixed block-structured mesh, with the air–water interface captured using the volume-of-fluid method. A new sub-grid scale stress model based on the vortex identifier is used to improve the robustness and efficiency of the simulation flows with air–water interface. The new framework is tested in the context of bow waves and Kelvin waves generated by a water-surface vehicle. Wave breaking at the bow of the vehicle is captured in LES. The LES results of wave geometry approaches the measurements progressively as the grid resolution is refined. The simulation results indicate that LES is a useful tool for studying wave dynamics of water-surface vehicles.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51206033 and 51276046)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20112302110020)+2 种基金the China Postdoctoral Science Foundation(Grant No.2011M500652)the Heilongjiang Postdoctoral Science Foundation,China(Grant No.2011LBH-Z11139)the Natural Scientific Research Innovation Foundation in Harbin Institute of Technology,China(Grant No.HIT.NSRIF.2012070)
文摘Large-eddy simulations (LES) based on the temporal approximate deconvolution model were performed for a forced homogeneous isotropic turbulence (FHIT) with polymer additives at moderate Taylor Reynolds number. Finitely extensible nonlinear elastic in the Peterlin approximation model was adopted as the constitutive equation for the filtered conformation tensor of the polymer molecules. The LES results were verified through comparisons with the direct numerical simulation results. Using the LES database of the FHIT in the Newtonian fluid and the polymer solution flows, the polymer effects on some important parameters such as strain, vorticity, drag reduction, and so forth were studied. By extracting the vortex structures and exploring the flatness factor through a high-order correlation function of velocity derivative and wavelet analysis, it can be found that the small-scale vortex structures and small-scale intermittency in the FHIT are all inhibited due to the existence of the polymers. The extended self-similarity scaling law in the polymer solution flow shows no apparent difference from that in the Newtonian fluid flow at the currently simulated ranges of Reynolds and Weissenberg numbers.
文摘The near-field dynamics of an aircraft wake is studied by means of temporal large-eddy simulations, with and without considering the effect of engine jets. In the absence of jets, the simulations showed the roll-up of the initial vorticity sheet shed by the wing and the occurrence of short-wavelength instability in a pair of primary co-rotating vortices. The main consequence of the instability is the modification of the internal structure of the vortex, compared to the two-dimensional stable behavior. The presence of engine jets affects the roll-up of the vorticity sheet and causes an enlargement of the final merged vortex core compared to the case without jets.
文摘The purpose of the present study was to establish a passive flow control method for a rectangular jet using two types of deflectors installed symmetrically inside a nozzle. This deflector in a rectangular nozzle generates the rectangular coaxial jets. The effect of the slant angle of the deflectors on the flow characteristics and the spread of the rectangular jet was investigated experimentally and by large-eddy simulation. The experiment and the numerical simulation were carried out at a Reynolds number of 9000. The rectangular jet with no deflectors generates a vortex ring from the nozzle exit. The vortex ring collapses in the downstream region and the outline of the jet changes from rectangular to diamond-shaped as a result of the axis-switching phenomenon. The rectangular jet with divergent and convergent deflectors shows particularly noticeable changes in the flow characteristics and vortical structures, as compared to the case with no deflectors. In the case of the rectangular jet with divergent deflectors (slant angle of <i style="font-family:" font-size:13.3333px;white-space:normal;"="">α</i><span style="font-family:" font-size:13.3333px;white-space:normal;"=""> = 6<span style="white-space:nowrap;">°</span>), minor axis spread is promoted more than major axis spread, and axis switching occurs closer to the nozzle exit than that in the case of no deflectors. The outline of the jet also changes from lateral rectangular to vertical rectangular as a result of axis switching. On the other hand, in the case of a rectangular jet with convergent deflectors (</span><i style="font-family:" font-size:13.3333px;white-space:normal;"="">α</i><span style="font-family:" font-size:13.3333px;white-space:normal;"=""> = -6<span style="white-space:nowrap;">°</span>), minor axis spread is suppressed more than major axis spread, and axis switching occurs farther from the nozzle exit than that in the case with no deflectors. The outline of the jet does not change until the downstream region. Therefore, the spread and the axis-switching location for the rectangular jet can be controlled by the deflectors inside the rectangular nozzle.</span>
基金funded by National Basic Research Program of China (Grant Nos. 2009CB421402 and 2010CB950503)National Natural Science Foundation of China (GrantNo. 40975007)
文摘In this study, the development of a convective boundary layer (CBL) in the Badaujilin region was investigated by comparing the observation data of two cases. A deep neutral layer capped a CBL that occurred on 30 August 2009. This case was divided into five sublayers from the surface to higher atmospheric elevations: surface layer, mixed layer, inversion layer, neutral layer, and sub-inversion layer. The development process of the CBL was divided into three stages: S1, S2, and S3. This case was quite different from the development of the three-layer CBL observed on 31 August 2009 because the mixed layer of the five-layer CBL (CBL5) eroded the neutral layer during S2. The specific initial structure of the CBL5 was correlated to the synoptic background of atmosphere during nighttime. The three-stage development process of the CBL5 was confirmed by six simulations using National Center for Atmospheric Research (USA) large-eddy simulation (NCAR-LES), and some of its characteristics are presented in detail.