The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and w...The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and wave fields are studied.The results show that the intense wind and wave fields further break the symmetries of submesoscale flow fields and suppress the levels of filament frontogenesis. The changes of secondary circulation directions—that is, the conversion between the convergence and divergence of the surface cross-filament currents with the downwelling and upwelling jets in the filament center—are associated with the inertial oscillation. The filament frontogenesis and frontolysis caused by the changes of secondary circulation directions may periodically sharpen and smooth the gradient of submesoscale flow fields.The lifecycle of the cold filament may include multiple stages of filament frontogenesis and frontolysis.展开更多
Shallow convection plays an important role in transporting heat and moisture from the near-surface to higher altitudes,yet its parameterization in numerical models remains a great challenge,partly due to the lack of h...Shallow convection plays an important role in transporting heat and moisture from the near-surface to higher altitudes,yet its parameterization in numerical models remains a great challenge,partly due to the lack of high-resolution observations.This study describes a large eddy simulation(LES)dataset for four shallow convection cases that differ primarily in inversion strength,which can be used as a surrogate for real data.To reduce the uncertainty in LES modeling,three different large eddy models were used,including SAM(System for Atmospheric Modeling),WRF(Weather Research and Forecasting model),and UCLA-LES.Results show that the different models generally exhibit similar behavior for each shallow convection case,despite some differences in the details of the convective structure.In addition to grid-averaged fields,conditionally sampled variables,such as in-cloud moisture and vertical velocity,are also provided,which are indispensable for calculation of the entrainment/detrainment rate.Considering the essentiality of the entraining/detraining process in the parameterization of cumulus convection,the dataset presented in this study is potentially useful for validation and improvement of the parameterization of shallow convection.展开更多
This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7...This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.展开更多
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.展开更多
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 vel展开更多
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.展开更多
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.展开更多
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.展开更多
A three dimensional numerical model in the σ coordinate system is developed to study the problem of waves. Turbulence effects are modeled by a subgrid scale (SGS) model with the concept of large eddy simulatio...A three dimensional numerical model in the σ coordinate system is developed to study the problem of waves. Turbulence effects are modeled by a subgrid scale (SGS) model with the concept of large eddy simulation (LES). The σ coordinate transformation is introduced to map the irregular physical domain of the wavy free surface and uneven bottom onto the regular computational domain of the shape of rectangular prism. The operator splitting method, which splits the solution procedure into the advection, diffusion, and propagation steps, is used to solve the modified Navier Stokes Equation. The model is used to simulate the propagation of solitary wave and wave passing over a submerged breakwater. Numerical results are compared with available analytical solutions and experimental data in terms of velocity profiles, free surface displacement, and energy conservation. Good agreement is obtained. The method is proved to be of high accuracy and efficiency in simulating surface wave propagation and wave structure interaction. It is suitable for the large and irregular physical domain, and requiring the non uniform grid system. The present work provides a foundation for further studies of random waves, wave structure interaction, wave discharge interaction, etc.展开更多
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 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.展开更多
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.展开更多
Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling facto...Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds(FGM) method is employed to determine major scalars. This coupled sub-grid scale(SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence.The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.展开更多
The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordin...The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordinate system is adopted. Turbulence is modeled by a dynamic coherent eddy model. The σ-coordinate transformation is introduced to map the irregular physical domain with a wavy free surface and an uneven bottom onto a regular computational domain. The fractional step method is used to solve the filtered Navier–Stokes equations. Results presented include the distribution of velocity, the decay law of the mean velocity along the jet axis, self-similar characteristics and volume flux per unit width. In particular, the role of coherent structures on the momentum transfer along the jet centerline and the jet instantaneous characteristics in JONSWAP waves are a special focus of this research. The numerical results obtained are of great theoretical importance in understanding the behavior of turbulent jets in random wave environments.展开更多
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.展开更多
This study uses a large eddy simulation (LES) model to investigate the turbulence processes in the ocean surface boundary layer at Zhangzi Island offshore. Field measurements at Zhangzi Island (39°N, 122°...This study uses a large eddy simulation (LES) model to investigate the turbulence processes in the ocean surface boundary layer at Zhangzi Island offshore. Field measurements at Zhangzi Island (39°N, 122°E) during July 2009 are used to drive the LES model. The LES results capture a clear diurnal cycle in the oceanic turbulence boundary layer. The process of the heat penetration and heat distribution characteristics are analyzed through the heat flux results from the LES and their differences between two diurnal cycles are discussed as well. Energy balance and other dynamics are investigated which show that the tide-induced shear production is the main source of the turbulence energy that balanced dissipation. Momentum flux near the surface shows better agreement with atmospheric data computed by the eddy correlation method than those computed by bulk formula.展开更多
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.展开更多
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.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 92158204, 41506001 and 42076019)a Project supported by the Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (Grant No. 311021005)。
文摘The variations of the frontogenetic trend of a cold filament induced by the cross-filament wind and wave fields are studied by a non-hydrostatic large eddy simulation. Five cases with different strengths of wind and wave fields are studied.The results show that the intense wind and wave fields further break the symmetries of submesoscale flow fields and suppress the levels of filament frontogenesis. The changes of secondary circulation directions—that is, the conversion between the convergence and divergence of the surface cross-filament currents with the downwelling and upwelling jets in the filament center—are associated with the inertial oscillation. The filament frontogenesis and frontolysis caused by the changes of secondary circulation directions may periodically sharpen and smooth the gradient of submesoscale flow fields.The lifecycle of the cold filament may include multiple stages of filament frontogenesis and frontolysis.
基金the National Key R&D Program of China(Grant No.2021YFC3000802)the National Natural Science Foundation of China(Grant No.42175165)the National Key Scientific and Technological Infrastructure project“Earth System Numerical Simulation Facility”(EarthLab).
文摘Shallow convection plays an important role in transporting heat and moisture from the near-surface to higher altitudes,yet its parameterization in numerical models remains a great challenge,partly due to the lack of high-resolution observations.This study describes a large eddy simulation(LES)dataset for four shallow convection cases that differ primarily in inversion strength,which can be used as a surrogate for real data.To reduce the uncertainty in LES modeling,three different large eddy models were used,including SAM(System for Atmospheric Modeling),WRF(Weather Research and Forecasting model),and UCLA-LES.Results show that the different models generally exhibit similar behavior for each shallow convection case,despite some differences in the details of the convective structure.In addition to grid-averaged fields,conditionally sampled variables,such as in-cloud moisture and vertical velocity,are also provided,which are indispensable for calculation of the entrainment/detrainment rate.Considering the essentiality of the entraining/detraining process in the parameterization of cumulus convection,the dataset presented in this study is potentially useful for validation and improvement of the parameterization of shallow convection.
文摘This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.
基金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.
基金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 vel
基金Project supported by the Boeing-COMAC Aviation Energy Conservation and Emissions Reduction Technology Center(AECER)
文摘A self-adaptive-grid method is applied to numerical simulation of the evolu- tion of aircraft wake vortex with the large eddy simulation (LES). The Idaho Falls (IDF) measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V^iis/il^i (BV) frequencies are com- puted with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, sec-ondary structures such as baroclinic vorticity and helical structures are also simulated. Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.
基金the 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 (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.
文摘A three dimensional numerical model in the σ coordinate system is developed to study the problem of waves. Turbulence effects are modeled by a subgrid scale (SGS) model with the concept of large eddy simulation (LES). The σ coordinate transformation is introduced to map the irregular physical domain of the wavy free surface and uneven bottom onto the regular computational domain of the shape of rectangular prism. The operator splitting method, which splits the solution procedure into the advection, diffusion, and propagation steps, is used to solve the modified Navier Stokes Equation. The model is used to simulate the propagation of solitary wave and wave passing over a submerged breakwater. Numerical results are compared with available analytical solutions and experimental data in terms of velocity profiles, free surface displacement, and energy conservation. Good agreement is obtained. The method is proved to be of high accuracy and efficiency in simulating surface wave propagation and wave structure interaction. It is suitable for the large and irregular physical domain, and requiring the non uniform grid system. The present work provides a foundation for further studies of random waves, wave structure interaction, wave discharge interaction, etc.
基金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.
基金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(Nos.91441117 and 51576182)
文摘A sub-grid scale(SGS) combustion model, which combines the artificial thickened flame(ATF) model with the flamelet generated manifold(FGM) tabulation method, is proposed. Based on the analysis of laminar flame structures, two self-contained flame sensors are used to track the diffusion and reaction processes with different spatial scales in the flame front, respectively. The dynamic formulation for the proposed SGS combustion model is also performed. Large eddy simulations(LESs) of Bunsen flame F3 are used to evaluate the different SGS combustion models. The results show that the proposed SGS model has the ability in predicting the distributions of temperature and velocity reasonably, while the predictions for the distributions of some species need further improvement. The snapshots of instantaneous normalized progress variables reveal that the flame is more remarkably and severely wrinkled at the flame tip for flame F3.More satisfactory results obtained by the dynamic model indicate that it can preserve the premixed flame propagation characteristics better.
基金Project supported by the National Natural Science Foundation of China(Nos.91441117 and51576182)
文摘Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds(FGM) method is employed to determine major scalars. This coupled sub-grid scale(SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence.The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.
基金supported by the National Natural Science Foundation of China (50679023, 50879019)Ph.D. Programs Foundation of Ministry of Education of China (20070294012)+2 种基金the National Science Fund for Distinguished Young Scholars (50925932)Outstanding Doctoral Dissertation Incubation Program of Hohai University (2010B18814)Qing Lan Project of Jiangsu Province, and 333 High-Level Talent Training Program of Jiangsu Province (2017-B08038)
文摘The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordinate system is adopted. Turbulence is modeled by a dynamic coherent eddy model. The σ-coordinate transformation is introduced to map the irregular physical domain with a wavy free surface and an uneven bottom onto a regular computational domain. The fractional step method is used to solve the filtered Navier–Stokes equations. Results presented include the distribution of velocity, the decay law of the mean velocity along the jet axis, self-similar characteristics and volume flux per unit width. In particular, the role of coherent structures on the momentum transfer along the jet centerline and the jet instantaneous characteristics in JONSWAP waves are a special focus of this research. The numerical results obtained are of great theoretical importance in understanding the behavior of turbulent jets in random wave environments.
基金The 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.
基金The National Basic Research Program of China under contract Nos 201 1CB403501 and 2012CB417402the Fund for Creative Research Groups by the National Natural Science Foundation of China under contract No.41121064+1 种基金the National Natural Science Foundation of China under contract Nos 41206015 and 41176016the Open Research Foundation for the State Key Laboratory of Satellite Ocean Environment Dynamics,Second Institute of Oceanography,State Oceanic Administration under contract No.SOED1210
文摘This study uses a large eddy simulation (LES) model to investigate the turbulence processes in the ocean surface boundary layer at Zhangzi Island offshore. Field measurements at Zhangzi Island (39°N, 122°E) during July 2009 are used to drive the LES model. The LES results capture a clear diurnal cycle in the oceanic turbulence boundary layer. The process of the heat penetration and heat distribution characteristics are analyzed through the heat flux results from the LES and their differences between two diurnal cycles are discussed as well. Energy balance and other dynamics are investigated which show that the tide-induced shear production is the main source of the turbulence energy that balanced dissipation. Momentum flux near the surface shows better agreement with atmospheric data computed by the eddy correlation method than those computed by bulk formula.
基金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.
基金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.