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.展开更多
A sub-grid scale(SGS) combustion model, which combines the artificial thickened flame(ATF) model with the flamelet generated manifold(FGM) tabulation method, is proposed. Based on the analysis of laminar flame structu...A sub-grid scale(SGS) combustion model, which combines the artificial thickened flame(ATF) model with the flamelet generated manifold(FGM) tabulation method, is proposed. Based on the analysis of laminar flame structures, two self-contained flame sensors are used to track the diffusion and reaction processes with different spatial scales in the flame front, respectively. The dynamic formulation for the proposed SGS combustion model is also performed. Large eddy simulations(LESs) of Bunsen flame F3 are used to evaluate the different SGS combustion models. The results show that the proposed SGS model has the ability in predicting the distributions of temperature and velocity reasonably, while the predictions for the distributions of some species need further improvement. The snapshots of instantaneous normalized progress variables reveal that the flame is more remarkably and severely wrinkled at the flame tip for flame F3.More satisfactory results obtained by the dynamic model indicate that it can preserve the premixed flame propagation characteristics better.展开更多
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.展开更多
Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling facto...Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds(FGM) method is employed to determine major scalars. This coupled sub-grid scale(SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence.The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.展开更多
A mixed subgrid-scale（SGS） model based on coherent structures and temporal approximate deconvolution（MCT） is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is...A mixed subgrid-scale（SGS） model based on coherent structures and temporal approximate deconvolution（MCT） is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation（LES） of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence（FHIT） with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation（DNS） results. Compared with the LES results using the temporal approximate deconvolution model（TADM） for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.展开更多
In this paper, large eddy simulation coupled with a dynamic subgrid scale (SGS) model is used to study turbulent channel flows with heat transfer. Some fundamental flow behaviors are discussed. Based on our calculate...In this paper, large eddy simulation coupled with a dynamic subgrid scale (SGS) model is used to study turbulent channel flows with heat transfer. Some fundamental flow behaviors are discussed. Based on our calculated results, the dynamic SGS model can reasonably predict some main characteristics of stratified turbulent channel flows. Our results are also in good agreement with theoretical analyses and previous calculated results.展开更多
A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a si...A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a sixth order compact differencing scheme utilizing the localized artificial diffusivity method for stabilizing shock waves and employing a wall model to enable the use of coarse mesh. In the validation study where a supersonic turbulent boundary layer flow over a flat plate is examined, both LES codes are well validated using velocity profile in the boundary layer given by the hot-wire anemometry and normal stress given by the laser Doppler anemometry. In particular, the sixth order compact differencing scheme gives closer agreements with these experimental data. Then, the validated LES codes are applied to solve the Mach number 2.5 supersonic turbulent flowfield involving the pseudo-shock waves. It is shown that typical features of unsteady flowfield of the pseudo-shock waves are well obtained by both schemes. Again, it is indicated that the sixth order compact differencing scheme gives closer agreements with the existing velocity data obtained by particle image velocimetry and pressure fluctuation data on the wall surface. Besides, the computational cost of the compact differencing scheme is found to be 1/7 of that for the upwind finite volume scheme, even though a wall model is solved at each grid point on the wall surface. Therefore, the obtained results in the present study allow recommending the sixth order compact differencing scheme with a wall model for simulating supersonic turbulent flowfield in an isolator involving the pseudo-shock waves.展开更多
Diesel spray is injected at high pressure. So, upper stream region of spray is high Weber number condition. However, even if the fuel is injected at high pressure, the downstream region of spray is corresponding to re...Diesel spray is injected at high pressure. So, upper stream region of spray is high Weber number condition. However, even if the fuel is injected at high pressure, the downstream region of spray is corresponding to relatively low Weber number condition. Thus, KH （Kelvin-Helmholtz） model modeled for high Weber number conditions and MTAB （modified Taylor analogy breakup） model are used for primary and secondary breakup processes respectively. This study is focused on the development of new hybrid breakup model The calculations are performed by LES （large eddy simulation） incorporated into KIVA code. LES of non-evaporating diesel spray are performed using KH ＆ RT （Rayleigh-Taylor） model, MTAB model and KH-MTAB model. Then, LES with these models were compared with experimental results. As the result, the availability of KH-MTAB model is showed. It is found that KH-MTAB is good agreement with experimental results of penetration and SMD （Sauter mean diameter） in relatively low density conditions.展开更多
Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained ...Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained using data with Mach number Ma=3.0 and Reynolds number Re=3000 was applied to situations with different Mach numbers and Reynolds numbers.The input variables of the neural network model were the filtered velocity gradients and temperature gradients at a single spatial grid point.The a priori test showed that the FCNN model had a correlation coefficient larger than 0.91 and a relative error smaller than 0.43,with much better reconstructions of SGS unclosed terms than the dynamic Smagorinsky model(DSM).In a posteriori test,the behavior of the FCNN model was marginally better than that of the DSM in predicting the mean velocity profiles,mean temperature profiles,turbulent intensities,total Reynolds stress,total Reynolds heat flux,and mean SGS flux of kinetic energy,and outperformed the Smagorinsky model.展开更多
An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating c...An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.展开更多
Large eddy simulations (LES) of mixing process in a stirred tank of 0.476m diameter with a 3-narrow blade hydrofoil CBY impeller were reported. The turbulent flow field and mixing time were calculated using LES with S...Large eddy simulations (LES) of mixing process in a stirred tank of 0.476m diameter with a 3-narrow blade hydrofoil CBY impeller were reported. The turbulent flow field and mixing time were calculated using LES with Sma-gorinsky-Lilly subgrid scale model. The impeller rotation was modeled using the sliding mesh technique. Better agree-ment of power demand and mixing time was obtained between the experimental and the LES prediction than that by the traditional Reynolds-averaged Navier-Stokes (RANS) approach. The curve of tracer response predicted by LES was in good agreement with the experimental. The results show that LES is a reliable tool to investigate the unsteady and quasi-periodic behavior of the turbulent flow in stirred tanks.展开更多
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.展开更多
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.展开更多
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.展开更多
The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordin...The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordinate system is adopted. Turbulence is modeled by a dynamic coherent eddy model. The σ-coordinate transformation is introduced to map the irregular physical domain with a wavy free surface and an uneven bottom onto a regular computational domain. The fractional step method is used to solve the filtered Navier–Stokes equations. Results presented include the distribution of velocity, the decay law of the mean velocity along the jet axis, self-similar characteristics and volume flux per unit width. In particular, the role of coherent structures on the momentum transfer along the jet centerline and the jet instantaneous characteristics in JONSWAP waves are a special focus of this research. The numerical results obtained are of great theoretical importance in understanding the behavior of turbulent jets in random wave environments.展开更多
The aim of this work is to investigate the flow instabilities in a baffled, stirred tank generated by a single Rushton turbine by means of large eddy simulation （LES）. The sliding mesh method was used for the coupli...The aim of this work is to investigate the flow instabilities in a baffled, stirred tank generated by a single Rushton turbine by means of large eddy simulation （LES）. The sliding mesh method was used for the coupling between the rotating and the stationary frame of references. The calculations were carried out on the ＂Shengcao-21C＂ supercomputer using a computational fluid dynamics （CFD） code CFX5. The flow fields predicted by the LES simulation and the simulation using standard κ-ε model were compared to the results from particle image velocimetry （PIV） measurements. It is shown that the CFD simulations using the LES approach and the standard κ-ε model agree well with the PIV measurements. Fluctuations of the radial and axial velocity are predicted at different frequencies by the LES simulation. Velocity fluctuations of high frequencies are seen in the impeller region, while low frequencies velocity fluctuations are observed in the bulk flow. A low frequency velocity fluctuation with a nondimensional frequency of 0.027Hz is predicted by the LES simulation, which agrees with experimental investigations in the literature. Flow circulation patterns predicted by the LES simulation are asymmetric, stochastic and complex, spanning a large portion of the tanks and varying with time, while circulation patterns calculated by the simulation using the standard κ-ε model are symmetric. The results of the present work give better understanding to the flow instabilities in the mechanically agitated tank. However, further analysis of the LES calculated velocity series by means of fast Fourier transform （FFT） and/or spectra analysis are recommended in future work in order to gain more knowledge of the complicated flow phenomena.展开更多
In this paper, the large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling lead...In this paper, the large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling leads to the increase of wave height, and then at a certain position, the wave will be breaking. The breaking wave is a powerful agent for generating turbulence, which plays an important role in most of the fluid dynamic processes throughout the surf zone, Such as transformation of wave energy, generation of near-shore current and diffusion of materials. So a proper numerical model for describing the turbulence effect is needed. In this paper, a revised Smagorinsky subgrid-scale model is used to describe the turbulence effect. The present study reveals that the coefficient of the Smagorinsky model for wave propagation or breaking simulation may be taken as a varying function of the water depth and distance away from the wave breaking point. The large eddy simulation model presented in this paper has been used to study the propagation of the solitary wave in constant water depth and the shoaling of the non-breaking solitary wave on a beach. The model is based on large eddy simulation, and to track free-surface movements, the Tokyo University Modified Marker and Cell (TUMMAC) method is employed. In order to ensure the accuracy of each component of this wave mathematical model, several steps have been taken to verify calculated solutions; with either analytical solutions or experimental data. For non-breaking waves, very accurate results are obtained for a solitary wave propagating over a constant depth and on a beach. Application of the model to cnoidal wave breaking in the surf zone shows that the model results are in good agreement with analytical solution and experimental data. From the present model results, it can be seen that the turbulent eddy viscosity increases from the bottom to the water surface in surf zone. In the eddy viscosity curve, there is a turn-point obviously, dividing water depth into two parts, in the upper part, the eddy viscosity becomes very large near the wave breaking position.展开更多
Viscous flow around a circular cylinder at a subcritical Reynolds number is investigated using a large eddy simulation (LES) coupled with the Smagorinsky subgrid-scale (SGS) model. A fractional-step method with a seco...Viscous flow around a circular cylinder at a subcritical Reynolds number is investigated using a large eddy simulation (LES) coupled with the Smagorinsky subgrid-scale (SGS) model. A fractional-step method with a second-order in time and a combined finite-difference/spectral approximations are used to solve the filtered three-dimensional incompressible Navier-Stokes equations. Calculations have been performed with and without the SGS model. Turbulence statistical behaviors and flow structures in the near wake of the cylinder are studied. Some calculated results, including the lift and drag coefficients, shedding frequency, peak Reynolds stresses, and time-average velocity profile, are in good agreement with the experimental and computational data, which shows that the Smagorinsky model can reasonably predict the global features of the flow and some turbulent statistical behaviors.展开更多
The effect of a cross-sectional exit plane on the downstream mixing characteristics of a circular turbulent jet is in- vestigated using large eddy simulation （LES）. The turbulent jet is issued from an orifice-type n...The effect of a cross-sectional exit plane on the downstream mixing characteristics of a circular turbulent jet is in- vestigated using large eddy simulation （LES）. The turbulent jet is issued from an orifice-type nozzle at an exit Reynolds number of 5 ×104. Both instantaneous and statistical velocity fields of the jet are provided. Results show that the rates of the mean velocity decay and jet spread are both higher in the case with the exit plate than without it. The existence of the plate is found to increase the downstream entrainment rate by about 10% on average over the axial range of 8-30de （exit diameter）. Also, the presence of the plate enables the formation of vortex rings to occur further downstream by 0.5-1 .Ode. A physical insight into the near-field jet is provided to explain the importance of the boundary conditions in the evolution of a turbulent jet. In addition, a method of using the decay of the centreline velocity and the half-width of the jet to calculate the entrainment rate is proposed.展开更多
Computational modeling methods have been increasingly employed to quantify aortic hemodynamic parameters that are challenging to in vivo measurements but important for the diagnosis/treatment of aortic disease.Althoug...Computational modeling methods have been increasingly employed to quantify aortic hemodynamic parameters that are challenging to in vivo measurements but important for the diagnosis/treatment of aortic disease.Although the presence of turbulence-like behaviors of blood flow in normal or diseased aorta has long been confirmed,the majority of existing computational model studies adopted the laminar flow assumption(LFA)in the treatment of sub-grid flow variables.So far,it remains unclear whether LFA would significantly compromise the reliability of hemodynamic simulation.In the present study,we addressed the issue in the context of a specific aortopathy,namely aortic dilation,which is usually accompanied by disturbed flow patterns.Three patient-specific aortas with treated/untreated dilation of the ascending segment were investigated,and their geometrical models were reconstructed from computed tomography angiographic images,with the boundary conditions being prescribed based on flow velocity information measured in vivo with the phase contrast magnetic resonance imaging technique.For the modeling of blood flow,apart from the traditional LFA-based method in which sub-grid flow dynamics is ignored,the large eddy simulation(LES)method capable of incorporating the dissipative energy loss induced by turbulent eddies at the sub-grid level,was adopted and taken as a reference for examining the performance of the LFA-based method.Obtained results showed that the simulated large-scale flow patterns with the two methods had high similarity,both agreeing well with in vivo measurements,although locally large between-method discrepancies in computed hemodynamic quantities existed in regions with high intensity of flow turbulence.Quantitatively,a switch from the LES to the LFAbased modeling method led to mild(<6%)changes in computed space-averaged wall shear stress metrics(i.e.,SA-TAWSS,SA-OSI)in the ascending aortic segment where intensive vortex evolution accompanied by high statistical Reynolds stress was observed.In addition,comparisons among the three aortas revealed that the treatment status of aortic dilation or the concomitant presence of aortic valve disease,despite its remarkable influence on flow patterns in the ascending aortic segment,did not significantly affect the degrees of discrepancies between the two modeling methods in predicting SA-TAWSS and SA-OSI.These findings suggest that aortic dilation per se does not induce strong flow turbulence that substantially negates the validity of LFA-based modeling,especially in simulating macro-scale hemodynamic features.展开更多
基金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.
基金Project supported by the National Natural Science Foundation of China(Nos.91441117 and 51576182)
文摘A sub-grid scale(SGS) combustion model, which combines the artificial thickened flame(ATF) model with the flamelet generated manifold(FGM) tabulation method, is proposed. Based on the analysis of laminar flame structures, two self-contained flame sensors are used to track the diffusion and reaction processes with different spatial scales in the flame front, respectively. The dynamic formulation for the proposed SGS combustion model is also performed. Large eddy simulations(LESs) of Bunsen flame F3 are used to evaluate the different SGS combustion models. The results show that the proposed SGS model has the ability in predicting the distributions of temperature and velocity reasonably, while the predictions for the distributions of some species need further improvement. The snapshots of instantaneous normalized progress variables reveal that the flame is more remarkably and severely wrinkled at the flame tip for flame F3.More satisfactory results obtained by the dynamic model indicate that it can preserve the premixed flame propagation characteristics better.
基金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 National Natural Science Foundation of China(Nos.91441117 and51576182)
文摘Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds(FGM) method is employed to determine major scalars. This coupled sub-grid scale(SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence.The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.
基金Project supported by the China Postdoctoral Science Foundation(Grant No.2011M500652)the National Natural Science Foundation of China(Grant Nos.51276046 and 51206033)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20112302110020)
文摘A mixed subgrid-scale（SGS） model based on coherent structures and temporal approximate deconvolution（MCT） is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation（LES） of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence（FHIT） with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation（DNS） results. Compared with the LES results using the temporal approximate deconvolution model（TADM） for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.
文摘In this paper, large eddy simulation coupled with a dynamic subgrid scale (SGS) model is used to study turbulent channel flows with heat transfer. Some fundamental flow behaviors are discussed. Based on our calculated results, the dynamic SGS model can reasonably predict some main characteristics of stratified turbulent channel flows. Our results are also in good agreement with theoretical analyses and previous calculated results.
文摘A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a sixth order compact differencing scheme utilizing the localized artificial diffusivity method for stabilizing shock waves and employing a wall model to enable the use of coarse mesh. In the validation study where a supersonic turbulent boundary layer flow over a flat plate is examined, both LES codes are well validated using velocity profile in the boundary layer given by the hot-wire anemometry and normal stress given by the laser Doppler anemometry. In particular, the sixth order compact differencing scheme gives closer agreements with these experimental data. Then, the validated LES codes are applied to solve the Mach number 2.5 supersonic turbulent flowfield involving the pseudo-shock waves. It is shown that typical features of unsteady flowfield of the pseudo-shock waves are well obtained by both schemes. Again, it is indicated that the sixth order compact differencing scheme gives closer agreements with the existing velocity data obtained by particle image velocimetry and pressure fluctuation data on the wall surface. Besides, the computational cost of the compact differencing scheme is found to be 1/7 of that for the upwind finite volume scheme, even though a wall model is solved at each grid point on the wall surface. Therefore, the obtained results in the present study allow recommending the sixth order compact differencing scheme with a wall model for simulating supersonic turbulent flowfield in an isolator involving the pseudo-shock waves.
文摘Diesel spray is injected at high pressure. So, upper stream region of spray is high Weber number condition. However, even if the fuel is injected at high pressure, the downstream region of spray is corresponding to relatively low Weber number condition. Thus, KH （Kelvin-Helmholtz） model modeled for high Weber number conditions and MTAB （modified Taylor analogy breakup） model are used for primary and secondary breakup processes respectively. This study is focused on the development of new hybrid breakup model The calculations are performed by LES （large eddy simulation） incorporated into KIVA code. LES of non-evaporating diesel spray are performed using KH ＆ RT （Rayleigh-Taylor） model, MTAB model and KH-MTAB model. Then, LES with these models were compared with experimental results. As the result, the availability of KH-MTAB model is showed. It is found that KH-MTAB is good agreement with experimental results of penetration and SMD （Sauter mean diameter） in relatively low density conditions.
基金Financial support provided by the National Natural Science Foundation of China(Grant Nos.11702042 and 91952104)。
文摘Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained using data with Mach number Ma=3.0 and Reynolds number Re=3000 was applied to situations with different Mach numbers and Reynolds numbers.The input variables of the neural network model were the filtered velocity gradients and temperature gradients at a single spatial grid point.The a priori test showed that the FCNN model had a correlation coefficient larger than 0.91 and a relative error smaller than 0.43,with much better reconstructions of SGS unclosed terms than the dynamic Smagorinsky model(DSM).In a posteriori test,the behavior of the FCNN model was marginally better than that of the DSM in predicting the mean velocity profiles,mean temperature profiles,turbulent intensities,total Reynolds stress,total Reynolds heat flux,and mean SGS flux of kinetic energy,and outperformed the Smagorinsky model.
基金the National Natural Science Foundation of China(50779069 and 90510007)the Start-up Scientific Research Foundation of China Agricultural University(2006021)the Beijing Natural Science Foundation(3071002).
文摘An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.
文摘Large eddy simulations (LES) of mixing process in a stirred tank of 0.476m diameter with a 3-narrow blade hydrofoil CBY impeller were reported. The turbulent flow field and mixing time were calculated using LES with Sma-gorinsky-Lilly subgrid scale model. The impeller rotation was modeled using the sliding mesh technique. Better agree-ment of power demand and mixing time was obtained between the experimental and the LES prediction than that by the traditional Reynolds-averaged Navier-Stokes (RANS) approach. The curve of tracer response predicted by LES was in good agreement with the experimental. The results show that LES is a reliable tool to investigate the unsteady and quasi-periodic behavior of the turbulent flow in stirred tanks.
基金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.
基金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.
文摘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.
基金supported by the National Natural Science Foundation of China (50679023, 50879019)Ph.D. Programs Foundation of Ministry of Education of China (20070294012)+2 种基金the National Science Fund for Distinguished Young Scholars (50925932)Outstanding Doctoral Dissertation Incubation Program of Hohai University (2010B18814)Qing Lan Project of Jiangsu Province, and 333 High-Level Talent Training Program of Jiangsu Province (2017-B08038)
文摘The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordinate system is adopted. Turbulence is modeled by a dynamic coherent eddy model. The σ-coordinate transformation is introduced to map the irregular physical domain with a wavy free surface and an uneven bottom onto a regular computational domain. The fractional step method is used to solve the filtered Navier–Stokes equations. Results presented include the distribution of velocity, the decay law of the mean velocity along the jet axis, self-similar characteristics and volume flux per unit width. In particular, the role of coherent structures on the momentum transfer along the jet centerline and the jet instantaneous characteristics in JONSWAP waves are a special focus of this research. The numerical results obtained are of great theoretical importance in understanding the behavior of turbulent jets in random wave environments.
文摘The aim of this work is to investigate the flow instabilities in a baffled, stirred tank generated by a single Rushton turbine by means of large eddy simulation （LES）. The sliding mesh method was used for the coupling between the rotating and the stationary frame of references. The calculations were carried out on the ＂Shengcao-21C＂ supercomputer using a computational fluid dynamics （CFD） code CFX5. The flow fields predicted by the LES simulation and the simulation using standard κ-ε model were compared to the results from particle image velocimetry （PIV） measurements. It is shown that the CFD simulations using the LES approach and the standard κ-ε model agree well with the PIV measurements. Fluctuations of the radial and axial velocity are predicted at different frequencies by the LES simulation. Velocity fluctuations of high frequencies are seen in the impeller region, while low frequencies velocity fluctuations are observed in the bulk flow. A low frequency velocity fluctuation with a nondimensional frequency of 0.027Hz is predicted by the LES simulation, which agrees with experimental investigations in the literature. Flow circulation patterns predicted by the LES simulation are asymmetric, stochastic and complex, spanning a large portion of the tanks and varying with time, while circulation patterns calculated by the simulation using the standard κ-ε model are symmetric. The results of the present work give better understanding to the flow instabilities in the mechanically agitated tank. However, further analysis of the LES calculated velocity series by means of fast Fourier transform （FFT） and/or spectra analysis are recommended in future work in order to gain more knowledge of the complicated flow phenomena.
基金This research project was supported by the National Natural Science Foundation of China and The Hong Kong Research Grants under contracts No. 59809006 and No. 59890200, also by the Science Foundation of Tianjin Municipality under contract No. 9837020
文摘In this paper, the large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling leads to the increase of wave height, and then at a certain position, the wave will be breaking. The breaking wave is a powerful agent for generating turbulence, which plays an important role in most of the fluid dynamic processes throughout the surf zone, Such as transformation of wave energy, generation of near-shore current and diffusion of materials. So a proper numerical model for describing the turbulence effect is needed. In this paper, a revised Smagorinsky subgrid-scale model is used to describe the turbulence effect. The present study reveals that the coefficient of the Smagorinsky model for wave propagation or breaking simulation may be taken as a varying function of the water depth and distance away from the wave breaking point. The large eddy simulation model presented in this paper has been used to study the propagation of the solitary wave in constant water depth and the shoaling of the non-breaking solitary wave on a beach. The model is based on large eddy simulation, and to track free-surface movements, the Tokyo University Modified Marker and Cell (TUMMAC) method is employed. In order to ensure the accuracy of each component of this wave mathematical model, several steps have been taken to verify calculated solutions; with either analytical solutions or experimental data. For non-breaking waves, very accurate results are obtained for a solitary wave propagating over a constant depth and on a beach. Application of the model to cnoidal wave breaking in the surf zone shows that the model results are in good agreement with analytical solution and experimental data. From the present model results, it can be seen that the turbulent eddy viscosity increases from the bottom to the water surface in surf zone. In the eddy viscosity curve, there is a turn-point obviously, dividing water depth into two parts, in the upper part, the eddy viscosity becomes very large near the wave breaking position.
基金The project supported by the National Science Fund for Distinguished Scholars (10125210)the Special Funds for Major State Basic Research Project (G1999032801)the National Natural Science Foundation of China (19772062)
文摘Viscous flow around a circular cylinder at a subcritical Reynolds number is investigated using a large eddy simulation (LES) coupled with the Smagorinsky subgrid-scale (SGS) model. A fractional-step method with a second-order in time and a combined finite-difference/spectral approximations are used to solve the filtered three-dimensional incompressible Navier-Stokes equations. Calculations have been performed with and without the SGS model. Turbulence statistical behaviors and flow structures in the near wake of the cylinder are studied. Some calculated results, including the lift and drag coefficients, shedding frequency, peak Reynolds stresses, and time-average velocity profile, are in good agreement with the experimental and computational data, which shows that the Smagorinsky model can reasonably predict the global features of the flow and some turbulent statistical behaviors.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11072005 and 10921202)the Fundamental Research Funds for the Central Universities,China(Grant No.3132013029)
文摘The effect of a cross-sectional exit plane on the downstream mixing characteristics of a circular turbulent jet is in- vestigated using large eddy simulation （LES）. The turbulent jet is issued from an orifice-type nozzle at an exit Reynolds number of 5 ×104. Both instantaneous and statistical velocity fields of the jet are provided. Results show that the rates of the mean velocity decay and jet spread are both higher in the case with the exit plate than without it. The existence of the plate is found to increase the downstream entrainment rate by about 10% on average over the axial range of 8-30de （exit diameter）. Also, the presence of the plate enables the formation of vortex rings to occur further downstream by 0.5-1 .Ode. A physical insight into the near-field jet is provided to explain the importance of the boundary conditions in the evolution of a turbulent jet. In addition, a method of using the decay of the centreline velocity and the half-width of the jet to calculate the entrainment rate is proposed.
基金The study was supported by the National Natural Science Foundation of China(Grant nos.11972231,11832003,81611530715)the China Postdoctoral Science Foundation(Grant no.2018M640385)the SJTU Medical-Engineering Cross-cutting Research Project(Grant no.YG2017MS45).
文摘Computational modeling methods have been increasingly employed to quantify aortic hemodynamic parameters that are challenging to in vivo measurements but important for the diagnosis/treatment of aortic disease.Although the presence of turbulence-like behaviors of blood flow in normal or diseased aorta has long been confirmed,the majority of existing computational model studies adopted the laminar flow assumption(LFA)in the treatment of sub-grid flow variables.So far,it remains unclear whether LFA would significantly compromise the reliability of hemodynamic simulation.In the present study,we addressed the issue in the context of a specific aortopathy,namely aortic dilation,which is usually accompanied by disturbed flow patterns.Three patient-specific aortas with treated/untreated dilation of the ascending segment were investigated,and their geometrical models were reconstructed from computed tomography angiographic images,with the boundary conditions being prescribed based on flow velocity information measured in vivo with the phase contrast magnetic resonance imaging technique.For the modeling of blood flow,apart from the traditional LFA-based method in which sub-grid flow dynamics is ignored,the large eddy simulation(LES)method capable of incorporating the dissipative energy loss induced by turbulent eddies at the sub-grid level,was adopted and taken as a reference for examining the performance of the LFA-based method.Obtained results showed that the simulated large-scale flow patterns with the two methods had high similarity,both agreeing well with in vivo measurements,although locally large between-method discrepancies in computed hemodynamic quantities existed in regions with high intensity of flow turbulence.Quantitatively,a switch from the LES to the LFAbased modeling method led to mild(<6%)changes in computed space-averaged wall shear stress metrics(i.e.,SA-TAWSS,SA-OSI)in the ascending aortic segment where intensive vortex evolution accompanied by high statistical Reynolds stress was observed.In addition,comparisons among the three aortas revealed that the treatment status of aortic dilation or the concomitant presence of aortic valve disease,despite its remarkable influence on flow patterns in the ascending aortic segment,did not significantly affect the degrees of discrepancies between the two modeling methods in predicting SA-TAWSS and SA-OSI.These findings suggest that aortic dilation per se does not induce strong flow turbulence that substantially negates the validity of LFA-based modeling,especially in simulating macro-scale hemodynamic features.
