Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, signific...Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, significantly extended design freedom and knowledge must be available in the early design stages. This paper presents a high-fidelity framework for design and optimization of the liquid swirl injectors that are widely used in aerospace propulsion and power-generation systems. The framework assembles a set of techniques, including Design Of Experiment(DOE), high-fidelity Large Eddy Simulations(LES), machine learning, Proper Orthogonal Decomposition(POD)-based Kriging surrogate modeling(emulation), inverse problem optimization, and uncertainty quantification. LES-based simulations can reveal detailed spatiotemporal evolution of flow structures and flame dynamics in a high-fidelity manner, and identify important injector design parameters according to their effects on propellant mixing, flame stabilization, and thermal protection.For a given a space of design parameters, DOE determines the number of design points to perform LES-based simulations. POD-based emulations, trained by the LES database, can effectively explore the design space and deduce an optimal group of design parameters in a turn-around time that is reduced by three orders of magnitude. The accuracy of the emulated results is validated, and the uncertainty of prediction is quantified. The proposed design methodology is expected to profoundly extend the knowledge base and reduce the cost for initial design stages.展开更多
Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering application...Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.展开更多
The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in...The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.展开更多
A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equ...A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.展开更多
A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation （LES）. The Idaho Falls （IDF） measurement of run 9 case is simulated numerica...A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation （LES）. The Idaho Falls （IDF） measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V^iis/il^i （BV） frequencies are com- puted with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, sec-ondary structures such as baroclinic vorticity and helical structures are also simulated. Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.展开更多
Large eddy simulations （LES） have been performed to investigate the flow and combustion fields in the scramjet of the German Aerospace Center （DLR）. Turbulent combustion is mod- eled by the tabulated thermo-chemis...Large eddy simulations （LES） have been performed to investigate the flow and combustion fields in the scramjet of the German Aerospace Center （DLR）. Turbulent combustion is mod- eled by the tabulated thermo-chemistry approach in combination with the presumed probability density function （PDF）. A/3-function is used to model the distribution of the mixture fraction, while two different PDFs, g-function （Model I） and //-function （Model II）, are applied to model the reaction progress. Temperature is obtained by solving filtered energy transport equation and the reaction rate of the progress variable is rescaled by pressure to consider the effects of compressibil- ity. The adaptive mesh refinement （AMR） technique is used to properly capture shock waves, boundary layers, shear layers and flame structures. Statistical results of temperature and velocity predicted by Model II show better accuracy than that predicted by Model I. The results of scatter points and mixture fraction-conditional variables indicate the significant differences between Model I and Model II. It is concluded that second moment information in the presumed PDF of the reaction progress is very important in the simulation of supersonic combustion. It is also found that an unstable flame with extinction and ignition develops in the shear layers of bluff body and a fuel- rich partially premixed flame stabilizes in the central recirculation bubble.展开更多
The noise induced by the fluctuant saturated steam flow under 250 °C in a stop-valve was numerically studied.The simulation was carried out using computational fluid dynamics（CFD） and ACTRAN.The acoustic field ...The noise induced by the fluctuant saturated steam flow under 250 °C in a stop-valve was numerically studied.The simulation was carried out using computational fluid dynamics（CFD） and ACTRAN.The acoustic field was investigated with Lighthill＇s acoustic analogy based on the properties of the flow field obtained using a large-eddy simulation that employs the LES-WALE dynamic model as the sub-grid-scale model.Firstly,the validation of mesh was well conducted,illustrating that two million elements were sufficient in this situation.Secondly,the treatment of the steam was deliberated,and conclusions indicate that when predicting the flow-induced noise of the stop-valve,the steam can be treated as incompressible gas at a low inlet velocity.Thirdly,the flow-induced noises under different inlet velocities were compared.The findings reveal it has remarkable influence on the flow-induced noises.Lastly,whether or not the heat preservation of the wall has influence on the noise was taken into account.The results show that heat preservation of the wall had little influence.展开更多
Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated ac...Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated accuracy has been evaluated using the LES verification and validation method.Various cavitation numbers are simulated to study important flow characteristics in the sphere wake,e.g.periodic cavity growth/contraction,interactions between the cloud and sheet cavitations and the vortex structure evolution.The spectral characteristics of the wake for typical cloud cavitation conditions were classified as the periodic cavitation mode,high Strouhal number mode and low Strouhal number mode.Main frequency distributions in the wake were analyzed and different dominant flow structures were identified for each of the three modes.Further,the cavitation and vortex relationship was also studied,which is an important issue associated with complex cavitating sphere wakes.Three types of cavitating vortex structures alternate,which indicates that three different cavity shedding regimes may exist in the wake.Analysis of vorticity transport equation shows a significant vorticity increase at the cavitation closure region and in the vortex cavitation region.This study provides a physical perspective to further understand the flow mechanisms in cavitating sphere wakes.展开更多
Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The...Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The modeling and analysis show that the LES model can simulate the planetary boundary layer (PBL) with a uniform underlying surface under various stratifications very well. Then, similar to the description of a forest canopy, the drag term on momentum and the production term of TKE by subgrid city buildings are introduced into the LES equations to account for the area-averaged effect of the subgrid urban canopy elements and to simulate the meteorological fields of the urban boundary layer (UBL). Numerical experiments and comparison analysis show that: (1) the result from the LES of the UBL with a proposed formula for the drag coefficient is consistent and comparable with that from wind tunnel experiments and an urban subdomain scale model; (2) due to the effect of urban buildings, the wind velocity near the canopy is decreased, turbulence is intensified, TKE, variance, and momentum flux are increased, the momentum and heat flux at the top of the PBL are increased, and the development of the PBL is quickened; (3) the height of the roughness sublayer (RS) of the actual city buildings is the maximum building height (1.5-3 times the mean building height), and a constant flux layer (CFL) exists in the lower part of the UBL.展开更多
In this paper,the Eulerian Stochastic Field(ESF)model in the Transported Probability Density Function(TPDF)class model is combined with the Flamelet Generated Manifolds(FGM)model.This method solves the joint probabili...In this paper,the Eulerian Stochastic Field(ESF)model in the Transported Probability Density Function(TPDF)class model is combined with the Flamelet Generated Manifolds(FGM)model.This method solves the joint probability density function transport equation by ESF method that considers the interaction mechanism between flame and turbulence with high precision.At the same time,by making use of the advantage of the FGM model,this model is able to incorporate the detailed chemical reaction mechanism(GRI 3.0)with acceptable computational cost.The new model has been implemented in the open source CFD suite-Open FOAM.Validation of the model has been carried out by simulating the Sandia flame series(three turbulent piloted methane jet flames)issued by the National Laboratory of the United States.The accuracy and advancement of the ESF/FGM turbulent combustion model are verified by comparing the LES results of the new model with the rich experimental data as well as the RANS results.The results demonstrate that the model has a strong ability in capturing combustion phenomena such as extinction and re-ignition in turbulent flame,which is essential in the accurate prediction of the combustion process in real combustion devices,for example,aircraft engines.展开更多
In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake d...In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.展开更多
A hybrid method is presented to numerically investigate the wind turbine aerodynamic characteristics.The wind turbine blade is replaced by an actuator line model.Turbulence is treated using a dynamic one-equation subg...A hybrid method is presented to numerically investigate the wind turbine aerodynamic characteristics.The wind turbine blade is replaced by an actuator line model.Turbulence is treated using a dynamic one-equation subgrid-scale model in large eddy simulation.Detailed information on the basic characteristics of the wind turbine wake is obtained and discussed.The rotor aerodynamic performance agrees well with the measurements.The actuator line method large-eddy simulation(ALM-LES)technique demonstrates its high potential in providing accurate load prediction and high resolution of turbulent fluctuations in the wind turbine wakes and the interactions within a feasible cost.展开更多
It is highly attractive to develop an efficient and flexible large eddy simulation(LES)technique for high-Reynolds-number atmospheric boundary layer(ABL)simulation using the low-order numerical scheme on a relatively ...It is highly attractive to develop an efficient and flexible large eddy simulation(LES)technique for high-Reynolds-number atmospheric boundary layer(ABL)simulation using the low-order numerical scheme on a relatively coarse grid,that could reproduce the logarithmic profile of the mean velocity and some key features of large-scale coherent structures in the outer layer.In this study,an improved near-wall correction scheme for the vertical gradient of the resolved streamwise velocity in the strain-rate tensor is proposed to calculate the eddy viscosity coefficient in the subgrid-scale(SGS)model.The LES code is realized with a second-order finite-difference scheme,the scale-dependent dynamic SGS stress model,the equilibrium wall stress model,and the proposed correction scheme.Very-high-Reynolds-number ABL flow simulation under the neutral stratification condition is conducted to assess the performance of the method in predicting the mean and fluctuating characteristics of the rough-wall turbulence.It is found that the logarithmic profile of the mean streamwise velocity and some key features of large-scale coherent structures can be reasonably predicted by adopting the proposed correction method and the low-order numerical scheme.展开更多
Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics h...Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.展开更多
Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the c...Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the characteristics of inflow turbulence,can significantly impact the quality of predictions.In this study,we examined these boundary conditions within the context of the mountainous terrain around a long-span cable-stayed bridge using a wind tunnel experiment.Various sizes of computational domains and turbulent incoming wind velocities were used in large eddy simulations.The results show that when the height of the computational domain is five times greater than the height of the terrain model,there is minimal influence from the top wall on the wind field characteristics in this complex mountainous area.Expanding the length of the wake region of the computational domain has negligible effects on the wind fields.Turbulence in the inlet boundary reduces the length of the wake region on a leeward hill with a low slope,but has less impact on the mean wind velocity of steep hills.展开更多
A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large edd...A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.展开更多
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.展开更多
Introducing the surface properties [initial vortex, ground temperature and surface momentum impact height （SMIH）] for the boundary conditions, dust-devil-scale large eddy simulations （LES） were carried out. Given ...Introducing the surface properties [initial vortex, ground temperature and surface momentum impact height （SMIH）] for the boundary conditions, dust-devil-scale large eddy simulations （LES） were carried out. Given three parameters of initial vortex, ground temperature and the SMIH based on Sinclair＇s observation, the dust devil physical characteristics, such as maximum tangential velocity, updraft velocity, pressure drop in the inner core region, and even reverse flow at the top of the core region, are predicted, and are found to be close to the observations, thus demonstrating the ability of the simulation. The physical characteristics of different modeled dust devils are reproduced and compared to the values predicted by Renno et al.＇ theory. Even for smaller temperature differences or weaker buoyancy, severe dust devils may be formed by strong incipient vortices. It is also indicated that SMIH substantially affects the near-surface shape of terrestrial dust devils.展开更多
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.展开更多
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.展开更多
基金sponsored by the William RT Oakes Endowment of the Georgia Institute of Technology
文摘Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, significantly extended design freedom and knowledge must be available in the early design stages. This paper presents a high-fidelity framework for design and optimization of the liquid swirl injectors that are widely used in aerospace propulsion and power-generation systems. The framework assembles a set of techniques, including Design Of Experiment(DOE), high-fidelity Large Eddy Simulations(LES), machine learning, Proper Orthogonal Decomposition(POD)-based Kriging surrogate modeling(emulation), inverse problem optimization, and uncertainty quantification. LES-based simulations can reveal detailed spatiotemporal evolution of flow structures and flame dynamics in a high-fidelity manner, and identify important injector design parameters according to their effects on propellant mixing, flame stabilization, and thermal protection.For a given a space of design parameters, DOE determines the number of design points to perform LES-based simulations. POD-based emulations, trained by the LES database, can effectively explore the design space and deduce an optimal group of design parameters in a turn-around time that is reduced by three orders of magnitude. The accuracy of the emulated results is validated, and the uncertainty of prediction is quantified. The proposed design methodology is expected to profoundly extend the knowledge base and reduce the cost for initial design stages.
基金supported by the National Natural Science Foundation of China(Grant No.50979059)
文摘Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.
基金Supported by the National Natural Science Foundation of China(10772082)~~
文摘The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.
基金Project supported by the National Natural Science Foundation of China(Nos.91752118,11672305,11232011,and 11572331)the Strategic Priority Research Program(No.XDB22040104)the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences(No.QYZDJ-SSWSYS002)
文摘A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.
基金Project supported by the Boeing-COMAC Aviation Energy Conservation and Emissions Reduction Technology Center(AECER)
文摘A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation （LES）. The Idaho Falls （IDF） measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V^iis/il^i （BV） frequencies are com- puted with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, sec-ondary structures such as baroclinic vorticity and helical structures are also simulated. Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.
基金financial support by the National Natural Science Foundation of China (Nos. 51176178, 50936005)
文摘Large eddy simulations （LES） have been performed to investigate the flow and combustion fields in the scramjet of the German Aerospace Center （DLR）. Turbulent combustion is mod- eled by the tabulated thermo-chemistry approach in combination with the presumed probability density function （PDF）. A/3-function is used to model the distribution of the mixture fraction, while two different PDFs, g-function （Model I） and //-function （Model II）, are applied to model the reaction progress. Temperature is obtained by solving filtered energy transport equation and the reaction rate of the progress variable is rescaled by pressure to consider the effects of compressibil- ity. The adaptive mesh refinement （AMR） technique is used to properly capture shock waves, boundary layers, shear layers and flame structures. Statistical results of temperature and velocity predicted by Model II show better accuracy than that predicted by Model I. The results of scatter points and mixture fraction-conditional variables indicate the significant differences between Model I and Model II. It is concluded that second moment information in the presumed PDF of the reaction progress is very important in the simulation of supersonic combustion. It is also found that an unstable flame with extinction and ignition develops in the shear layers of bluff body and a fuel- rich partially premixed flame stabilizes in the central recirculation bubble.
文摘The noise induced by the fluctuant saturated steam flow under 250 °C in a stop-valve was numerically studied.The simulation was carried out using computational fluid dynamics（CFD） and ACTRAN.The acoustic field was investigated with Lighthill＇s acoustic analogy based on the properties of the flow field obtained using a large-eddy simulation that employs the LES-WALE dynamic model as the sub-grid-scale model.Firstly,the validation of mesh was well conducted,illustrating that two million elements were sufficient in this situation.Secondly,the treatment of the steam was deliberated,and conclusions indicate that when predicting the flow-induced noise of the stop-valve,the steam can be treated as incompressible gas at a low inlet velocity.Thirdly,the flow-induced noises under different inlet velocities were compared.The findings reveal it has remarkable influence on the flow-induced noises.Lastly,whether or not the heat preservation of the wall has influence on the noise was taken into account.The results show that heat preservation of the wall had little influence.
基金This work was financially supported by the National Natural Science Foundation of China(Grants 51822903 and 11772239)the Natural Science Foundation of Hubei Province(Grant 2018CFA010).The numerical calculations were done on the supercomputing system in the Supercomputing Center of Wuhan University.
文摘Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated accuracy has been evaluated using the LES verification and validation method.Various cavitation numbers are simulated to study important flow characteristics in the sphere wake,e.g.periodic cavity growth/contraction,interactions between the cloud and sheet cavitations and the vortex structure evolution.The spectral characteristics of the wake for typical cloud cavitation conditions were classified as the periodic cavitation mode,high Strouhal number mode and low Strouhal number mode.Main frequency distributions in the wake were analyzed and different dominant flow structures were identified for each of the three modes.Further,the cavitation and vortex relationship was also studied,which is an important issue associated with complex cavitating sphere wakes.Three types of cavitating vortex structures alternate,which indicates that three different cavity shedding regimes may exist in the wake.Analysis of vorticity transport equation shows a significant vorticity increase at the cavitation closure region and in the vortex cavitation region.This study provides a physical perspective to further understand the flow mechanisms in cavitating sphere wakes.
基金The research was supported by the National Natural Science Foundation of China under Grant Nos.40333027 and 40075004.
文摘Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The modeling and analysis show that the LES model can simulate the planetary boundary layer (PBL) with a uniform underlying surface under various stratifications very well. Then, similar to the description of a forest canopy, the drag term on momentum and the production term of TKE by subgrid city buildings are introduced into the LES equations to account for the area-averaged effect of the subgrid urban canopy elements and to simulate the meteorological fields of the urban boundary layer (UBL). Numerical experiments and comparison analysis show that: (1) the result from the LES of the UBL with a proposed formula for the drag coefficient is consistent and comparable with that from wind tunnel experiments and an urban subdomain scale model; (2) due to the effect of urban buildings, the wind velocity near the canopy is decreased, turbulence is intensified, TKE, variance, and momentum flux are increased, the momentum and heat flux at the top of the PBL are increased, and the development of the PBL is quickened; (3) the height of the roughness sublayer (RS) of the actual city buildings is the maximum building height (1.5-3 times the mean building height), and a constant flux layer (CFL) exists in the lower part of the UBL.
基金supported by the National Natural Science Foundation of China(No.51706241).
文摘In this paper,the Eulerian Stochastic Field(ESF)model in the Transported Probability Density Function(TPDF)class model is combined with the Flamelet Generated Manifolds(FGM)model.This method solves the joint probability density function transport equation by ESF method that considers the interaction mechanism between flame and turbulence with high precision.At the same time,by making use of the advantage of the FGM model,this model is able to incorporate the detailed chemical reaction mechanism(GRI 3.0)with acceptable computational cost.The new model has been implemented in the open source CFD suite-Open FOAM.Validation of the model has been carried out by simulating the Sandia flame series(three turbulent piloted methane jet flames)issued by the National Laboratory of the United States.The accuracy and advancement of the ESF/FGM turbulent combustion model are verified by comparing the LES results of the new model with the rich experimental data as well as the RANS results.The results demonstrate that the model has a strong ability in capturing combustion phenomena such as extinction and re-ignition in turbulent flame,which is essential in the accurate prediction of the combustion process in real combustion devices,for example,aircraft engines.
基金Project supported by the National Key Research and Development Program of China(Nos.2019YFE0192600,2017YFE0132000,and 2019YFB1503700)the National Natural Science Foundation of China(Nos.51761135012 and 11872248)。
文摘In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.
基金funded jointly by the National Basic Research Program of China(″973″Program)(No.2014CB046200)the Jiangsu Provincial Natural Science Foundation(No.BK20140059)+2 种基金the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe National Natural Science Foundation of China(No.11172135)the EU Seventh Framework Program(No.FP7-PEOPLE-2010-IRSES-269202)
文摘A hybrid method is presented to numerically investigate the wind turbine aerodynamic characteristics.The wind turbine blade is replaced by an actuator line model.Turbulence is treated using a dynamic one-equation subgrid-scale model in large eddy simulation.Detailed information on the basic characteristics of the wind turbine wake is obtained and discussed.The rotor aerodynamic performance agrees well with the measurements.The actuator line method large-eddy simulation(ALM-LES)technique demonstrates its high potential in providing accurate load prediction and high resolution of turbulent fluctuations in the wind turbine wakes and the interactions within a feasible cost.
基金Project supported by the National Natural Science Foundation of China(No.11490551)the Fundamental Research Funds for the Central Universities(Nos.lzujbky-2016-k13 and lzujbky-2018-k07)
文摘It is highly attractive to develop an efficient and flexible large eddy simulation(LES)technique for high-Reynolds-number atmospheric boundary layer(ABL)simulation using the low-order numerical scheme on a relatively coarse grid,that could reproduce the logarithmic profile of the mean velocity and some key features of large-scale coherent structures in the outer layer.In this study,an improved near-wall correction scheme for the vertical gradient of the resolved streamwise velocity in the strain-rate tensor is proposed to calculate the eddy viscosity coefficient in the subgrid-scale(SGS)model.The LES code is realized with a second-order finite-difference scheme,the scale-dependent dynamic SGS stress model,the equilibrium wall stress model,and the proposed correction scheme.Very-high-Reynolds-number ABL flow simulation under the neutral stratification condition is conducted to assess the performance of the method in predicting the mean and fluctuating characteristics of the rough-wall turbulence.It is found that the logarithmic profile of the mean streamwise velocity and some key features of large-scale coherent structures can be reasonably predicted by adopting the proposed correction method and the low-order numerical scheme.
文摘Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.
基金supported by the National Natural Science Foundation of China(Nos.51925808 and 52178516)the Natural Science Foundation of Hunan Province(Nos.2020JJ5745 and 2023JJ20073),China.
文摘Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the characteristics of inflow turbulence,can significantly impact the quality of predictions.In this study,we examined these boundary conditions within the context of the mountainous terrain around a long-span cable-stayed bridge using a wind tunnel experiment.Various sizes of computational domains and turbulent incoming wind velocities were used in large eddy simulations.The results show that when the height of the computational domain is five times greater than the height of the terrain model,there is minimal influence from the top wall on the wind field characteristics in this complex mountainous area.Expanding the length of the wake region of the computational domain has negligible effects on the wind fields.Turbulence in the inlet boundary reduces the length of the wake region on a leeward hill with a low slope,but has less impact on the mean wind velocity of steep hills.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51479116 and 11272213)
文摘A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.
基金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.
文摘Introducing the surface properties [initial vortex, ground temperature and surface momentum impact height （SMIH）] for the boundary conditions, dust-devil-scale large eddy simulations （LES） were carried out. Given three parameters of initial vortex, ground temperature and the SMIH based on Sinclair＇s observation, the dust devil physical characteristics, such as maximum tangential velocity, updraft velocity, pressure drop in the inner core region, and even reverse flow at the top of the core region, are predicted, and are found to be close to the observations, thus demonstrating the ability of the simulation. The physical characteristics of different modeled dust devils are reproduced and compared to the values predicted by Renno et al.＇ theory. Even for smaller temperature differences or weaker buoyancy, severe dust devils may be formed by strong incipient vortices. It is also indicated that SMIH substantially affects the near-surface shape of terrestrial dust devils.
基金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(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.
