This study investigates the turbulence-induced disturbances and stall precursor triggering mechanism in NACA65-18(10)cascade based on large eddy simulations.The results indicate that the disturbances exist under vario...This study investigates the turbulence-induced disturbances and stall precursor triggering mechanism in NACA65-18(10)cascade based on large eddy simulations.The results indicate that the disturbances exist under various operating conditions along the performance curve.The shear layer is the physical structure responsible for the generation,propagation,and dissipation of disturbances.When operating near stall,the separation on the suction surface intensifies,and strong unsteady backflow occurs at the trailing edge of the passage.Under the influence of inlet disturbances,unsteady behaviors between passages form specific phase differences,leading the entire system to oscillate in a first-order mode.As the flow develops from near-stall to stall,axial momentum decreases further,reducing the main flow’s ability to drive blockages downstream through convection.Consequently,the blockage accumulates during the circumferential propagation process until the stall onset.Based on the above mechanism,this study proposes factors describing the size of the backflow zone,shedding frequency,and convection velocity to characterize blockage dynamics,identifying critical values that represent the stall onset.展开更多
A numerical framework is developed for dust transport by coupling the regularized lattice Boltzmannmethod with large eddy simulation in an Eulerian-Eulerian formulation.The model incorporates wall-modeled LES through ...A numerical framework is developed for dust transport by coupling the regularized lattice Boltzmannmethod with large eddy simulation in an Eulerian-Eulerian formulation.The model incorporates wall-modeled LES through a wall shear stress treatment and adopts a recursive regularization strategy to enhance numerical stability for high Reynolds turbulence.Dust transport is modeled using an additional distribution function governed by an eddydiffusion-based paradigm.The model is validated using two benchmark cases:the inviscid Taylor-Green vortex with dust and a high-Reynolds-number turbulent channel flow.It is then applied to simulate dust transport in the atmospheric boundary layer,where large-scale coherent structures drive lateral dispersion and influence vertical transport and particle residence time.The effects of particle size are systematically examined,showing that larger particles exhibit stronger settling due to inertia,while smaller particles remain suspended through enhanced coupling with turbulence.A negative correlation between streamwise velocity and particle concentration is also identified,highlighting the interaction between turbulent flow structures and particulate transport.展开更多
At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures...At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures.A framework based on the interior penalty discontinuous Galerkin method and large eddy simulation approach was adopted in the present study.The performances of various subgrid models,including the Smagorinsky(SM)model,the dynamic Smagorinsky(DSM)model,the wall-adapting local-eddy-viscosity(WALE)model,and the VREMAN model,have been analyzed through flow simulations of the SD7003 airfoil at a Reynolds number of 60000.It turns out that the SM model fails to predict the emergence of LSB,even modified by the Van-Driest damping function.On the contrary,the best agreement is generally achieved by the WALE model in terms of flow separation,reattachment,and transition locations,together with the aerodynamic loads.Furthermore,the influence of numerical dissipation has also been discussed through the comparison of skin friction and resolved Reynolds stresses.As numerical dissipation decreases,the prediction accuracy of the WALE model degrades.Meanwhile,nonlinear variation could be observed from the performances of the DSM model,which could be attributed to the interaction between the numerical dissipation and the subgrid model.展开更多
The air-sea interactions with the submesoscale warm filament of the oceanic mixed layer are simulated by a coupled atmosphere-ocean model of the parallelized lager eddy simulation model.The results show that the warm ...The air-sea interactions with the submesoscale warm filament of the oceanic mixed layer are simulated by a coupled atmosphere-ocean model of the parallelized lager eddy simulation model.The results show that the warm core of the oceanic warm filament heats the bottom air of the atmospheric boundary layer,the rise of the bottom warm air results in the formation of the atmospheric warm filament.The variation in the width of the oceanic warm filament is generated by the change in the direction of the secondary circulations.The variation in the width of the atmospheric warm filament is created by that of the oceanic warm filament,because the direction of the secondary circulations of the atmospheric warm filament is invariable with time.The Coriolis effect results in the change in the direction of the secondary circulations for the oceanic warm filament.The secondary circulations of the atmospheric warm filament are produced by the rise of the bottom warm air caused by the oceanic warm filament,which leads to the unchanged direction of the secondary circulations.The thermal convection turbulence caused by the temperature difference of the ocean and atmosphere gradually weakens the structure of the oceanic and atmospheric warm filaments.展开更多
This paper investigates the influence of numerical methods and mesh resolution on the prediction accuracy of the aerodynamic behaviors of a 1/20 scaled generic high-speed train(HST)model.A thorough comparison is made ...This paper investigates the influence of numerical methods and mesh resolution on the prediction accuracy of the aerodynamic behaviors of a 1/20 scaled generic high-speed train(HST)model.A thorough comparison is made between partially averaged Navier-Stokes(PANS),large eddy simulation(LES),and wind tunnel experiments,covering aerodynamic forces,surface pressure,velocity distribution,and Reynolds stress and turbulent kinetic energy in the wake region.The Reynolds number for both simulations and experiments is set to 4.75×105.The results show that the PANS approach accurately predicts flow characteristics observed in experiments and fine LES calculations,even with a low resolution grid.PANS exhibits a distinct advantage over LES when grid resolutions are insufficient for resolving near wall flow structures around the HST,both in open-air conditions and crosswind environments.Additionally,grid refinement improves the predictive accuracy of the HST's aerodynamic performance,particularly in the presence of small yaw angle.展开更多
Coagulation and growth of nanoparticles subject to large coherent structures in a planar jet has been explored by using large eddy simulation. The particle field is obtained by employing a moment method to approximate...Coagulation and growth of nanoparticles subject to large coherent structures in a planar jet has been explored by using large eddy simulation. The particle field is obtained by employing a moment method to approximate the nanoparticle general dynamic equa- tion. An incompressible fluid containing particles of 1 nm in diameter is projected into a particle-free ambient. The results show that the coherent structures dominate the evolution of the nanoparticle number intensity, diameter and polydispersity distributions as the jet develops. In addition, the coherent structures act to increase the diffusion of particles, and the vortex rolling-up makes the particles distributing more irregularly while the vortex pairing causes particle distributions to become uniform. As the jet travels downstream, the time-averaged particle number concentration becomes lower in the jet core and higher in the outskirts, whereas the time- averaged particle mass over the entire flow field maintains unaltered, and the time-averaged particle diameter and geometric standard deviations grow and reach their maximum on the interface of the jet region and the ambient.展开更多
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 large eddy simulation (LES) of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the ...A large eddy simulation (LES) of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.展开更多
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.展开更多
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.展开更多
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.展开更多
Fourier neural operator(FNO)model is developed for large eddy simulation(LES)of three-dimensional(3D)turbulence.Velocity fields of isotropic turbulence generated by direct numerical simulation(DNS)are used for trainin...Fourier neural operator(FNO)model is developed for large eddy simulation(LES)of three-dimensional(3D)turbulence.Velocity fields of isotropic turbulence generated by direct numerical simulation(DNS)are used for training the FNO model to predict the filtered velocity field at a given time.The input of the FNO model is the filtered velocity fields at the previous several time-nodes with large time lag.In the a posteriori study of LES,the FNO model performs better than the dynamic Smagorinsky model(DSM)and the dynamic mixed model(DMM)in the prediction of the velocity spectrum,probability density functions(PDFs)of vorticity and velocity increments,and the instantaneous flow structures.Moreover,the proposed model can significantly reduce the computational cost,and can be well generalized to LES of turbulence at higher Taylor-Reynolds numbers.展开更多
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.展开更多
In this article, numerical investigation of the effects of different plasma actuation strengths on the film cooling flow characteristics has been conducted using large eddy simulation (LES). For this numerical resea...In this article, numerical investigation of the effects of different plasma actuation strengths on the film cooling flow characteristics has been conducted using large eddy simulation (LES). For this numerical research, the plasma actuator is placed downstream of the trailing edge of the film cooling hole and a phenomenological model is employed to provide the electric field generated by it, resulting in the body forces. Our results show that as the plasma actuation strength grows larger, under the downward effect of the plasma actuation, the jet trajectory near the cooling hole stays closer to the wall and the recirculation region observably reduces in size. Meanwhile, the momentum injection effect of the plasma actuation also actively alters the distributions of the velocity components downstream of the cooling hole. Consequently, the influence of the plasma actuation strength on the Reynolds stress downstream of the cooling hole is remarkable. Furthermore, the plasma actuation weakens the strength of the kidney shaped vortex and prevents the jet from lifting off the wall. Therefore, with the increase of the strength of the plasma actuation, the coolant core stays closer to the wall and tends to split into two distinct regions. So the centerline film cooling efficiency is enhanced, and it is increased by 55% at most when the plasma actuation strength is 10.展开更多
The turbulence enhancement by particle wake effect is studied by large eddy simulation (LES) of turbulent gas flows passing a single particle. The predicted time-averaged and root-mean-square fluctuation velocities ...The turbulence enhancement by particle wake effect is studied by large eddy simulation (LES) of turbulent gas flows passing a single particle. The predicted time-averaged and root-mean-square fluctuation velocities behind the particle are in agreement with the Reynolds-averaged Navier-Stokes modeling results and experimental results. A semi-empirical turbulence enhancement model is proposed by the present-authors based on the LES resuits. This model is incorporated into the second-order moment two-phase turbulence model for simulating vertical gas-particle pipe flows and horizontal gas-particle channel flows. The simulation results show that compared with the model not accounting for the particle wake effect, the present model gives simulation results for the gas turbulence modulation in much better agreement with the experimental results.展开更多
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.展开更多
Supersonic fuel film cooling is a promising way to simultaneously reduce the severe wall heat and friction load of the internal passage in a scramjet engine when operating at hypersonic conditions.Large eddy simulatio...Supersonic fuel film cooling is a promising way to simultaneously reduce the severe wall heat and friction load of the internal passage in a scramjet engine when operating at hypersonic conditions.Large eddy simulations were performed to investigate the cooling and wall friction characteristics of hydrogen and hydrocarbon films under inert and reactive circumstances.The results show that the essential difference of the turbulent state in the mixing layer contributes to the totally different behaviors of the cooling and wall friction reduction performances of the two fuel films.The turbulent transport processes between the hydrogen film and the mainstream are much weaker as compared to the case of hydrocarbon film,making inert hydrogen rather superior in cooling and friction reduction applications.Besides,the increase of wall temperature for hydrogen film under the inert case is mainly driven by the loss of hydrogen with high heat capacity instead of by direct heat addition.However,the film cooling performance severely deteriorates when the hydrogen film burns due to presence of severe heat release sources near the wall.On the other hand,combustion of hydrocarbon film in the boundary layer can remarkably improve its originally barely-satisfactory cooling and friction reduction performance to the level comparable to that of hydrogen film,due to the suppression of turbulent transport processes in the mixing layer and presence of heat absorption sources near the wall.Overall,the hydrogen film is more advantageous in friction reduction,while the hydrocarbon film is more suitable for cooling.展开更多
In this study, the large eddy simulation technique was applied on the flow in a baffled stirred tank driven by a Rushton turbine at Re=29000. The interaction between the rotating impeller and the static baffles was ac...In this study, the large eddy simulation technique was applied on the flow in a baffled stirred tank driven by a Rushton turbine at Re=29000. The interaction between the rotating impeller and the static baffles was accounted for by means of the improved inner-outer iterative algorithm. The sub-grid scale model was a conventional Smagorinsky model. The numerical solution of the governing equations was conducted in a cylindrical staggered grid. The momentum and the continuity equations were discretized using the finite difference method, with a third-order QUICK scheme used for convective terms. The phase-resolved predictions were compared with the experimental data of Wu and Patterson and good agreement was observed for both the mean and the turbulence quantities. They were much better than the Reynolds-averaged Navier-Stokes model including the Reynolds Stress Model for simulating the turbulence. The study also suggests the feasibility of LES in combination with the improved inner-outer iterative algorithm for the simulation of turbulent flow in stirred tanks.展开更多
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.展开更多
An investigation of large-eddy simulation(LES) for turbulent channel flow with buoyancy effects was performed by solving the resolved incompressible Navier-Stokes equations under the Boussinesq approximation. The Smag...An investigation of large-eddy simulation(LES) for turbulent channel flow with buoyancy effects was performed by solving the resolved incompressible Navier-Stokes equations under the Boussinesq approximation. The Smagorinsky eddy-viscosity model and Yoshizawa eddy-viscosity model were used to describe the unresolved subgrid scale (SGS) fluctuations respectively After some numerical testing, the latter was further simplified so that if can be used in the dynamic model closure. A LES code was developed for parallel computations by using the parallel technique, and aas run on the Dawn-1000 parallel computer. To demonstrate the viability and accuracy of the code, our results are compared with and found in good agreement with available LES results.展开更多
基金support of the National Natural Science Foundation of China(Nos.52322603 and U24A20141)the Science Center for Gas Turbine Project of China(No.P2023-B-Ⅱ-001-001)the Fundamental Research Funds for the Central Universities of China and the Beijing Nova Program of China(Nos.20220484074 and 20230484479)。
摘要This study investigates the turbulence-induced disturbances and stall precursor triggering mechanism in NACA65-18(10)cascade based on large eddy simulations.The results indicate that the disturbances exist under various operating conditions along the performance curve.The shear layer is the physical structure responsible for the generation,propagation,and dissipation of disturbances.When operating near stall,the separation on the suction surface intensifies,and strong unsteady backflow occurs at the trailing edge of the passage.Under the influence of inlet disturbances,unsteady behaviors between passages form specific phase differences,leading the entire system to oscillate in a first-order mode.As the flow develops from near-stall to stall,axial momentum decreases further,reducing the main flow’s ability to drive blockages downstream through convection.Consequently,the blockage accumulates during the circumferential propagation process until the stall onset.Based on the above mechanism,this study proposes factors describing the size of the backflow zone,shedding frequency,and convection velocity to characterize blockage dynamics,identifying critical values that represent the stall onset.
基金supported by the National Natural Science Foundation of China(Grant Nos.12572316 and 12388101)the Natural Science Foundation of Heilongjiang Province(Grant No.LH2021 A007).
摘要A numerical framework is developed for dust transport by coupling the regularized lattice Boltzmannmethod with large eddy simulation in an Eulerian-Eulerian formulation.The model incorporates wall-modeled LES through a wall shear stress treatment and adopts a recursive regularization strategy to enhance numerical stability for high Reynolds turbulence.Dust transport is modeled using an additional distribution function governed by an eddydiffusion-based paradigm.The model is validated using two benchmark cases:the inviscid Taylor-Green vortex with dust and a high-Reynolds-number turbulent channel flow.It is then applied to simulate dust transport in the atmospheric boundary layer,where large-scale coherent structures drive lateral dispersion and influence vertical transport and particle residence time.The effects of particle size are systematically examined,showing that larger particles exhibit stronger settling due to inertia,while smaller particles remain suspended through enhanced coupling with turbulence.A negative correlation between streamwise velocity and particle concentration is also identified,highlighting the interaction between turbulent flow structures and particulate transport.
基金This work was supported by the National Key R&D Program of China(Grant No.2022YFE0207000)the National Natural Science Foundation of China(Grant Nos.12372289,11972250,and 12102298)+1 种基金the China Postdoctoral Science Foundation(Grant No.2021M702443)Tianjin Natural Science Foundation(Grant No.22JCZDJC00910).
摘要At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures.A framework based on the interior penalty discontinuous Galerkin method and large eddy simulation approach was adopted in the present study.The performances of various subgrid models,including the Smagorinsky(SM)model,the dynamic Smagorinsky(DSM)model,the wall-adapting local-eddy-viscosity(WALE)model,and the VREMAN model,have been analyzed through flow simulations of the SD7003 airfoil at a Reynolds number of 60000.It turns out that the SM model fails to predict the emergence of LSB,even modified by the Van-Driest damping function.On the contrary,the best agreement is generally achieved by the WALE model in terms of flow separation,reattachment,and transition locations,together with the aerodynamic loads.Furthermore,the influence of numerical dissipation has also been discussed through the comparison of skin friction and resolved Reynolds stresses.As numerical dissipation decreases,the prediction accuracy of the WALE model degrades.Meanwhile,nonlinear variation could be observed from the performances of the DSM model,which could be attributed to the interaction between the numerical dissipation and the subgrid model.
基金The Scientific Research Fund for Doctor/Professor of Hezhou University under contract No.2024BSQD01the National Natural Science Foundation of China under contract No.42466001the project supported by Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)under contract No.SML2023SP240。
摘要The air-sea interactions with the submesoscale warm filament of the oceanic mixed layer are simulated by a coupled atmosphere-ocean model of the parallelized lager eddy simulation model.The results show that the warm core of the oceanic warm filament heats the bottom air of the atmospheric boundary layer,the rise of the bottom warm air results in the formation of the atmospheric warm filament.The variation in the width of the oceanic warm filament is generated by the change in the direction of the secondary circulations.The variation in the width of the atmospheric warm filament is created by that of the oceanic warm filament,because the direction of the secondary circulations of the atmospheric warm filament is invariable with time.The Coriolis effect results in the change in the direction of the secondary circulations for the oceanic warm filament.The secondary circulations of the atmospheric warm filament are produced by the rise of the bottom warm air caused by the oceanic warm filament,which leads to the unchanged direction of the secondary circulations.The thermal convection turbulence caused by the temperature difference of the ocean and atmosphere gradually weakens the structure of the oceanic and atmospheric warm filaments.
基金Project(2024YFB4303300)supported by the National Key R&D Program of ChinaProject(52202429)supported by the National Natural Science Foundation of ChinaProject(2023JJ40747)supported by the Natural Science Foundation of Hunan Province,China。
摘要This paper investigates the influence of numerical methods and mesh resolution on the prediction accuracy of the aerodynamic behaviors of a 1/20 scaled generic high-speed train(HST)model.A thorough comparison is made between partially averaged Navier-Stokes(PANS),large eddy simulation(LES),and wind tunnel experiments,covering aerodynamic forces,surface pressure,velocity distribution,and Reynolds stress and turbulent kinetic energy in the wake region.The Reynolds number for both simulations and experiments is set to 4.75×105.The results show that the PANS approach accurately predicts flow characteristics observed in experiments and fine LES calculations,even with a low resolution grid.PANS exhibits a distinct advantage over LES when grid resolutions are insufficient for resolving near wall flow structures around the HST,both in open-air conditions and crosswind environments.Additionally,grid refinement improves the predictive accuracy of the HST's aerodynamic performance,particularly in the presence of small yaw angle.
基金The project was supported by the National Natural Science Foundation of China (10372090)the Doctoral Program of Higher Education of China (20030335001)
摘要Coagulation and growth of nanoparticles subject to large coherent structures in a planar jet has been explored by using large eddy simulation. The particle field is obtained by employing a moment method to approximate the nanoparticle general dynamic equa- tion. An incompressible fluid containing particles of 1 nm in diameter is projected into a particle-free ambient. The results show that the coherent structures dominate the evolution of the nanoparticle number intensity, diameter and polydispersity distributions as the jet develops. In addition, the coherent structures act to increase the diffusion of particles, and the vortex rolling-up makes the particles distributing more irregularly while the vortex pairing causes particle distributions to become uniform. As the jet travels downstream, the time-averaged particle number concentration becomes lower in the jet core and higher in the outskirts, whereas the time- averaged particle mass over the entire flow field maintains unaltered, and the time-averaged particle diameter and geometric standard deviations grow and reach their maximum on the interface of the jet region and the ambient.
基金supported by the National Natural Science Foundation of China(Nos.51925808 and 52178516)the Natural Science Foundation of Hunan Province(Nos.2020JJ5745 and 2023JJ20073),China.
摘要Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the characteristics of inflow turbulence,can significantly impact the quality of predictions.In this study,we examined these boundary conditions within the context of the mountainous terrain around a long-span cable-stayed bridge using a wind tunnel experiment.Various sizes of computational domains and turbulent incoming wind velocities were used in large eddy simulations.The results show that when the height of the computational domain is five times greater than the height of the terrain model,there is minimal influence from the top wall on the wind field characteristics in this complex mountainous area.Expanding the length of the wake region of the computational domain has negligible effects on the wind fields.Turbulence in the inlet boundary reduces the length of the wake region on a leeward hill with a low slope,but has less impact on the mean wind velocity of steep hills.
基金supported by the National Natural Science Foundation of China (11302238, 11232011. and 11572331)support from the Strategic Priority Research Program (XDB22040104)+1 种基金the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences (QYZDJ-SSW-SYS002)the National Basic Research Program of China (973 Program 2013CB834100: Nonlinear science)
摘要A large eddy simulation (LES) of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.
基金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.
基金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 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.
基金supported by the National Natural Science Foundation of China(Nos.91952104,92052301,12172161,and 12161141017)National Numerical Windtunnel Project(No.NNW2019ZT1-A04)+4 种基金Shenzhen Science and Technology Program(No.KQTD20180411143441009)Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory(Guangzhou)(No.GML2019ZD0103)CAAI-Huawei Mind Spore open Fundand by Department of Science and Technology of Guangdong Province(No.2019B21203001)supported by Center for Computational Science and Engineering of Southern University of Science and Technology。
摘要Fourier neural operator(FNO)model is developed for large eddy simulation(LES)of three-dimensional(3D)turbulence.Velocity fields of isotropic turbulence generated by direct numerical simulation(DNS)are used for training the FNO model to predict the filtered velocity field at a given time.The input of the FNO model is the filtered velocity fields at the previous several time-nodes with large time lag.In the a posteriori study of LES,the FNO model performs better than the dynamic Smagorinsky model(DSM)and the dynamic mixed model(DMM)in the prediction of the velocity spectrum,probability density functions(PDFs)of vorticity and velocity increments,and the instantaneous flow structures.Moreover,the proposed model can significantly reduce the computational cost,and can be well generalized to LES of turbulence at higher Taylor-Reynolds numbers.
基金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.
摘要In this article, numerical investigation of the effects of different plasma actuation strengths on the film cooling flow characteristics has been conducted using large eddy simulation (LES). For this numerical research, the plasma actuator is placed downstream of the trailing edge of the film cooling hole and a phenomenological model is employed to provide the electric field generated by it, resulting in the body forces. Our results show that as the plasma actuation strength grows larger, under the downward effect of the plasma actuation, the jet trajectory near the cooling hole stays closer to the wall and the recirculation region observably reduces in size. Meanwhile, the momentum injection effect of the plasma actuation also actively alters the distributions of the velocity components downstream of the cooling hole. Consequently, the influence of the plasma actuation strength on the Reynolds stress downstream of the cooling hole is remarkable. Furthermore, the plasma actuation weakens the strength of the kidney shaped vortex and prevents the jet from lifting off the wall. Therefore, with the increase of the strength of the plasma actuation, the coolant core stays closer to the wall and tends to split into two distinct regions. So the centerline film cooling efficiency is enhanced, and it is increased by 55% at most when the plasma actuation strength is 10.
基金Supported by the Major Project of National Natural Science Foundation of China (No.10632070) the Postdoctoral ScienceFoundation (No.2004036239).
摘要The turbulence enhancement by particle wake effect is studied by large eddy simulation (LES) of turbulent gas flows passing a single particle. The predicted time-averaged and root-mean-square fluctuation velocities behind the particle are in agreement with the Reynolds-averaged Navier-Stokes modeling results and experimental results. A semi-empirical turbulence enhancement model is proposed by the present-authors based on the LES resuits. This model is incorporated into the second-order moment two-phase turbulence model for simulating vertical gas-particle pipe flows and horizontal gas-particle channel flows. The simulation results show that compared with the model not accounting for the particle wake effect, the present model gives simulation results for the gas turbulence modulation in much better agreement with the experimental results.
基金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.
基金supported by the National Natural Science Foundation of China(No.52176037).
摘要Supersonic fuel film cooling is a promising way to simultaneously reduce the severe wall heat and friction load of the internal passage in a scramjet engine when operating at hypersonic conditions.Large eddy simulations were performed to investigate the cooling and wall friction characteristics of hydrogen and hydrocarbon films under inert and reactive circumstances.The results show that the essential difference of the turbulent state in the mixing layer contributes to the totally different behaviors of the cooling and wall friction reduction performances of the two fuel films.The turbulent transport processes between the hydrogen film and the mainstream are much weaker as compared to the case of hydrocarbon film,making inert hydrogen rather superior in cooling and friction reduction applications.Besides,the increase of wall temperature for hydrogen film under the inert case is mainly driven by the loss of hydrogen with high heat capacity instead of by direct heat addition.However,the film cooling performance severely deteriorates when the hydrogen film burns due to presence of severe heat release sources near the wall.On the other hand,combustion of hydrocarbon film in the boundary layer can remarkably improve its originally barely-satisfactory cooling and friction reduction performance to the level comparable to that of hydrogen film,due to the suppression of turbulent transport processes in the mixing layer and presence of heat absorption sources near the wall.Overall,the hydrogen film is more advantageous in friction reduction,while the hydrocarbon film is more suitable for cooling.
基金Supported by the National Natural Science Foundation of China (No.20236050) and the State Key Development Program for Basic Research of China (No.2004CB217604).
摘要In this study, the large eddy simulation technique was applied on the flow in a baffled stirred tank driven by a Rushton turbine at Re=29000. The interaction between the rotating impeller and the static baffles was accounted for by means of the improved inner-outer iterative algorithm. The sub-grid scale model was a conventional Smagorinsky model. The numerical solution of the governing equations was conducted in a cylindrical staggered grid. The momentum and the continuity equations were discretized using the finite difference method, with a third-order QUICK scheme used for convective terms. The phase-resolved predictions were compared with the experimental data of Wu and Patterson and good agreement was observed for both the mean and the turbulence quantities. They were much better than the Reynolds-averaged Navier-Stokes model including the Reynolds Stress Model for simulating the turbulence. The study also suggests the feasibility of LES in combination with the improved inner-outer iterative algorithm for the simulation of turbulent flow in stirred tanks.
摘要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 project supported by the National Natural Science Foundation of Chinaby the Youngster Funding of Academia Sinica
摘要An investigation of large-eddy simulation(LES) for turbulent channel flow with buoyancy effects was performed by solving the resolved incompressible Navier-Stokes equations under the Boussinesq approximation. The Smagorinsky eddy-viscosity model and Yoshizawa eddy-viscosity model were used to describe the unresolved subgrid scale (SGS) fluctuations respectively After some numerical testing, the latter was further simplified so that if can be used in the dynamic model closure. A LES code was developed for parallel computations by using the parallel technique, and aas run on the Dawn-1000 parallel computer. To demonstrate the viability and accuracy of the code, our results are compared with and found in good agreement with available LES results.