Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics h...Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.展开更多
The current research of large eddy simulation （LES） of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale （SGS） model to simulate turbulent flow, which aims at obtaining the flow fi...The current research of large eddy simulation （LES） of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale （SGS） model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry （PIV） and laser Doppler velocimetry （LDV）. The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.展开更多
Realistic hydrological response is sensitive to the spatial variability of landscape properties. For a grid-based distributed rainfall-runoff model with a hypothesis of a uniform grid, the high-frequency information w...Realistic hydrological response is sensitive to the spatial variability of landscape properties. For a grid-based distributed rainfall-runoff model with a hypothesis of a uniform grid, the high-frequency information within a grid cell will be gradually lost as the resolution of the digital elevation model(DEM) grows coarser. Therefore, the performance of a hydrological model is usually scale-dependent. This study used the GridXinanjiang(GXAJ) model as an example to investigate the effects of subgrid variability on hydrological response at different scales. With the aim of producing a more reasonable hydrological response and spatial description of the landscape properties, a new distributed rainfall-runoff model integrating the subgrid variability(the GXAJSV model) was developed. In this model, the topographic index is used as an auxiliary variable correlated with the soil storage capacity. The incomplete beta distribution is suggested for simulating the probability distribution of the soil storage capacity within the raster grid. The Yaogu Basin in China was selected for model calibration and validation at different spatial scales.Results demonstrated that the proposed model can effectively eliminate the scale dependence of the GXAJ model and produce a more reasonable hydrological response.展开更多
The subgrid-scale effects on particle motion were investigated in forced isotropic turbulence by DNS and priorLES methods.In the DNS field,the importance of Kolmogorov scaling to preferential accumulation was validate...The subgrid-scale effects on particle motion were investigated in forced isotropic turbulence by DNS and priorLES methods.In the DNS field,the importance of Kolmogorov scaling to preferential accumulation was validated by comparing the radial distribution functions under various particle Stokes numbers.The prior-LES fields were generated by filtering the DNS data.The subgrid-scale Stokes number(StSGS)is a useful tool for determining the effects of subgrid-scale eddies on particle motion.The subgrid-scale eddies tend to accumulate particles with StSGSb 1 and disperse particles with 1 b StSGSb 10.For particles with StSGS?1,the effects of subgrid-scale eddies on particle motion can be neglected.In order to restore the subgrid-scale effects,the Langevin-type stochastic model with optimized parameters was adopted in this study.This model is effective for the particles with StSGS N 1 while has an adverse impact on the particles with StSGSb 1.The results show that the Langevin-type stochastic model tends to smooth the particle distribution in the isotropic turbulence.展开更多
In this paper, large eddy simulation coupled with a dynamic subgrid scale (SGS) model is used to study turbulent channel flows with heat transfer. Some fundamental flow behaviors are discussed. Based on our calculate...In this paper, large eddy simulation coupled with a dynamic subgrid scale (SGS) model is used to study turbulent channel flows with heat transfer. Some fundamental flow behaviors are discussed. Based on our calculated results, the dynamic SGS model can reasonably predict some main characteristics of stratified turbulent channel flows. Our results are also in good agreement with theoretical analyses and previous calculated results.展开更多
A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution di...A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.展开更多
A mixed subgrid-scale（SGS） model based on coherent structures and temporal approximate deconvolution（MCT） is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is...A mixed subgrid-scale（SGS） model based on coherent structures and temporal approximate deconvolution（MCT） is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation（LES） of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence（FHIT） with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation（DNS） results. Compared with the LES results using the temporal approximate deconvolution model（TADM） for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.展开更多
We review the previous attempts of rational subgrid-scale （SGS） modelling by employing theKolmogorov equation of filtered quantities. Aiming at explaining and solving the underlyingproblems in these models, we ...We review the previous attempts of rational subgrid-scale （SGS） modelling by employing theKolmogorov equation of filtered quantities. Aiming at explaining and solving the underlyingproblems in these models, we also introduce the recent methodological investigations for therational SGS modelling technique by defining the terms of assumption and restriction. Thesemethodological works are expected to provide instructive criterions for not only the rational SGSmodelling, but also other types of SGS modelling practices.展开更多
Vortices play a fundamental role in fluid dynamics,but mathematically defining them remains elusive.While many vortex identification methods are scalar-valued,vortices are inherently rotational,vector-based phenomena....Vortices play a fundamental role in fluid dynamics,but mathematically defining them remains elusive.While many vortex identification methods are scalar-valued,vortices are inherently rotational,vector-based phenomena.Liutex,as a vector quantity,addresses these limitations by accurately capturing the local rotational characteristics of fluid elements while remaining independent of shear influences.This unique property makes Liutex particularly well-suited for vortex identification and the quantitative analysis of turbulent flows.This paper explores the statistical analysis of Liutex in various turbulence regimes and proposes an objective Liutex-based vortex identification method.The objective method is rooted in the statistical properties of Liutex.Furthermore,the paper investigates the performance of Liutex-based subgrid models in large eddy simulation(LES).The effectiveness of these models is evaluated by comparing their performance in different flow conditions,such as decaying homogeneous isotropic turbulence and turbulent channel flows,against conventional models.Results demonstrate that the inclusion of Liutex significantly enhances the ability of subgrid models to accurately capture flow structures.Importantly,the new model maintains the same form regardless of whether strong or weak shear is present,ensuring robustness and consistency in both vortex identification and turbulence modeling.These findings highlight the significant potential of Liutex to improve turbulence modeling in both theoretical and practical contexts,with ongoing research aimed at further refining its theoretical foundations and expanding its application in more complex flow scenarios.展开更多
In numerical weather prediction(NWP),the parameterization of orographic drag plays an important role in representing subgrid orographic effects.The subgrid orographic parameters are the key input to the parameterizati...In numerical weather prediction(NWP),the parameterization of orographic drag plays an important role in representing subgrid orographic effects.The subgrid orographic parameters are the key input to the parameterization of orographic drag.Currently,the subgrid orographic parameters in most NWP models were produced based on elevation datasets generated many years ago,with a coarse resolution and low quality.In this paper,using the latest high-quality elevation data and considering the applicable scale range of the subgrid orographic parameters,we construct the orographic parameters,including the subgrid orographic standard deviation,anisotropy,orientation,and slope,that are required as input to the orographic gravity wave drag(OGWD)parameterization.Finally,we introduce the newly constructed orographic parameters into the Yin-He Global Spectral Model(YHGSM),optimize the description of the orographic effect in the model,and improve the simulation of two typical heavy rainfall events in Beijing and Henan.展开更多
Mesoscale eddies play a pivotal role in deciphering the intricacies of ocean dynamics and the transport of heat,salt,and nutrients.Accurate representation of these eddies in ocean models is essential for improving mod...Mesoscale eddies play a pivotal role in deciphering the intricacies of ocean dynamics and the transport of heat,salt,and nutrients.Accurate representation of these eddies in ocean models is essential for improving model predictions.In this study,we propose a convolutional neural network(CNN)that combines data-driven techniques with physical principles to develop a robust and interpretable parameterization scheme for mesoscale eddies in ocean modeling.We use a highresolution reanalysis dataset to extract subgrid eddy momentum and then applying machine learning algorithms to identify patterns and correlations.To ensure physical consistency,we have introduced conservation of momentum constraints in our CNN parameterization scheme through soft and hard constraints.The interpretability analysis illustrate that the pre-trained CNN parameterization shows promising results in accurately solving the resolved mean velocity and effectively capturing the representation of unresolved subgrid turbulence processes.Furthermore,to validate the CNN parameterization scheme offline,we conduct simulations using the Massachusetts Institute of Technology general circulation model(MITgcm)ocean model.A series of experiments is conducted to compare the performance of the model with the CNN parameterization scheme and high-resolution simulations.The offline validation demonstrates the effectiveness of the CNN parameterization scheme in improving the representation of mesoscale eddies in the MITgcm ocean model.Incorporating the CNN parameterization scheme leads to better agreement with high-resolution simulations and a more accurate representation of the kinetic energy spectra.展开更多
We propose a hybrid quantum-classical method,the quantum-enriched large eddy simulation(QELES),for simulating turbulence.The QELES combines the large-scale motion of the large eddy simulation(LES)and the subgrid motio...We propose a hybrid quantum-classical method,the quantum-enriched large eddy simulation(QELES),for simulating turbulence.The QELES combines the large-scale motion of the large eddy simulation(LES)and the subgrid motion of the incompressible Schrodinger flow(ISF).The ISF is a possible way to be simulated on a quantum computer,and it generates subgrid scale turbu-lent structures to enrich the LES field.The enriched LES field can be further used in turbulent combustion and multi-phase flows in which the subgrid scale motion plays an important role.As a conceptual study,we perform the simulations of ISF and LES separately on a classical computer to simulate decaying homogeneous isotropic turbulence.Then,the QEI ES velocity is obtained by the time matching and the spectral blending methods.The QEL ES achieves significant improvement in predicting the energy spectrum,probaility density functions of velocity and vorticity components,and velocity structure functions,and reconstructs coherent small-scales vortices in the direct numerical simulation(DNS).On the other hand,the vortices in the QELES are less elongated and tangled than those in the DNS,and the magnitude of the third-order structure function in the QELES is less than that in the DNS,due to the diferent constitutive relations in the viscous flow and ISE.展开更多
We consider a new subgrid eddy viscosity method based on pressure projection and extrapolated trapezoidal rule for the transient Navier-Stokes equations by using lowest equal-order pair of finite elements. The scheme ...We consider a new subgrid eddy viscosity method based on pressure projection and extrapolated trapezoidal rule for the transient Navier-Stokes equations by using lowest equal-order pair of finite elements. The scheme stabilizes convection dominated problems and ameliorates the restrictive inf-sup compatibility stability. It has some attractive fea- tures including parameter free for the pressure stabilized term and calculations required for higher order derivatives. Moreover, it requires only the solutions of the linear system arising from an Oseen problem per time step and has second order temporal accuracy. The method achieves optimal accuracy with respect to solution regularity.展开更多
The Ensemble Kalman Filter (EnKF) is well known and widely used in land data assimilation for its high precision and simple operation. The land surface models used as the forecast operator in a land data assimilation ...The Ensemble Kalman Filter (EnKF) is well known and widely used in land data assimilation for its high precision and simple operation. The land surface models used as the forecast operator in a land data assimilation system are usually designed to consider the model subgrid-heterogeneity and soil water thawing and freezing. To neglect their effects could lead to some errors in soil moisture assimilation. The dual EnKF method is employed in soil moisture data assimilation to build a soil moisture data as- similation framework based on the NCAR Community Land Model version 2.0 (CLM 2.0) in considera- tion of the effects of the model subgrid-heterogeneity and soil water thawing and freezing: Liquid volumetric soil moisture content in a given fraction is assimilated through the state filter process, while solid volumetric soil moisture content in the same fraction and solid/liquid volumetric soil moisture in the other fractions are optimized by the parameter filter. Preliminary experiments show that this dual EnKF-based assimilation framework can assimilate soil moisture more effectively and precisely than the usual EnKF-based assimilation framework without considering the model subgrid-scale heteroge- neity and soil water thawing and freezing. With the improvement of soil moisture simulation, the soil temperature-simulated precision can be also improved to some extent.展开更多
The concept of vortex is crucial in both understanding and modeling of turbulence.For large eddy simulation(LES),the effect of small-scale eddies onto the large scales or the resolved flow field is modeled by subgrid ...The concept of vortex is crucial in both understanding and modeling of turbulence.For large eddy simulation(LES),the effect of small-scale eddies onto the large scales or the resolved flow field is modeled by subgrid stress models.Even though the rotating motions of fluids,i.e.,vortices or eddies are central in developing turbulent models,vortex identification methods are seldom used in these models.In this study,we develop a new subgrid model based on the Liutex vector,a new quantity introduced to decompose fluid motions into rigid rotation,pure shear and stretching,and thus identify vortices.The new model is then applied in a decaying homogeneous isotropic turbulence(DHIT)and a turbulent channel flow at Reynolds number Reτ=180.It is shown that the new model can predict accurate energy spectra compared with experiments in DHIT and give a well-matched velocity profile in turbulent channel flow without changing the form of the model.Future directions include improvement of the Liutex based model,for example developing anisotropic subgrid models,and its applications in various turbulent flows.展开更多
In recent years, the Lattice Boltzmann Method (LBM) has developed into an alternative and promising numerical scheme for simulating fluid flows and modeling physics in fluids. In order to propose LBM for high Reynolds...In recent years, the Lattice Boltzmann Method (LBM) has developed into an alternative and promising numerical scheme for simulating fluid flows and modeling physics in fluids. In order to propose LBM for high Reynolds number fluid flow applications, a subgrid turbulence model for LBM was introduced based on standard Smagorinsky subgrid model and Lattice Bhatnagar-Gross-Krook (LBGK) model. The subgrid LBGK model was subsequently used to simulate the two-dimensional driven cavity flow at high Reynolds numbers. The simulation results including distribution of stream lines, dimensionless velocities distribution, values of stream function, as well as location of vertex center, were compared with benchmark solutions, with satisfactory agreements.展开更多
文摘Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.
基金supported by National Natural Science Foundation of China (Grant Nos. 51139007, 51079151, 51079152)Research Fundfor the Doctoral Program of Higher Education of China (Grant No. 0100008110012)
文摘The current research of large eddy simulation （LES） of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale （SGS） model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry （PIV） and laser Doppler velocimetry （LDV）. The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.
基金supported by the Graduate Students Scientific Research Innovation Plan of Jiangsu Province(Grant No.CXZZ12_0243)the China Scholarship Council(Grant No.201306710013)+1 种基金the State Major Project of Water Pollution Control and Management(Grant No.2014ZX07101-011)the Special Fund for Public Welfare of Ministry of Water Resources of China(Grant No.201501007)
文摘Realistic hydrological response is sensitive to the spatial variability of landscape properties. For a grid-based distributed rainfall-runoff model with a hypothesis of a uniform grid, the high-frequency information within a grid cell will be gradually lost as the resolution of the digital elevation model(DEM) grows coarser. Therefore, the performance of a hydrological model is usually scale-dependent. This study used the GridXinanjiang(GXAJ) model as an example to investigate the effects of subgrid variability on hydrological response at different scales. With the aim of producing a more reasonable hydrological response and spatial description of the landscape properties, a new distributed rainfall-runoff model integrating the subgrid variability(the GXAJSV model) was developed. In this model, the topographic index is used as an auxiliary variable correlated with the soil storage capacity. The incomplete beta distribution is suggested for simulating the probability distribution of the soil storage capacity within the raster grid. The Yaogu Basin in China was selected for model calibration and validation at different spatial scales.Results demonstrated that the proposed model can effectively eliminate the scale dependence of the GXAJ model and produce a more reasonable hydrological response.
基金Supported by the National Natural Science Foundation of China(51761125011).
文摘The subgrid-scale effects on particle motion were investigated in forced isotropic turbulence by DNS and priorLES methods.In the DNS field,the importance of Kolmogorov scaling to preferential accumulation was validated by comparing the radial distribution functions under various particle Stokes numbers.The prior-LES fields were generated by filtering the DNS data.The subgrid-scale Stokes number(StSGS)is a useful tool for determining the effects of subgrid-scale eddies on particle motion.The subgrid-scale eddies tend to accumulate particles with StSGSb 1 and disperse particles with 1 b StSGSb 10.For particles with StSGS?1,the effects of subgrid-scale eddies on particle motion can be neglected.In order to restore the subgrid-scale effects,the Langevin-type stochastic model with optimized parameters was adopted in this study.This model is effective for the particles with StSGS N 1 while has an adverse impact on the particles with StSGSb 1.The results show that the Langevin-type stochastic model tends to smooth the particle distribution in the isotropic turbulence.
文摘In this paper, large eddy simulation coupled with a dynamic subgrid scale (SGS) model is used to study turbulent channel flows with heat transfer. Some fundamental flow behaviors are discussed. Based on our calculated results, the dynamic SGS model can reasonably predict some main characteristics of stratified turbulent channel flows. Our results are also in good agreement with theoretical analyses and previous calculated results.
文摘A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.
基金Project supported by the China Postdoctoral Science Foundation(Grant No.2011M500652)the National Natural Science Foundation of China(Grant Nos.51276046 and 51206033)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20112302110020)
文摘A mixed subgrid-scale（SGS） model based on coherent structures and temporal approximate deconvolution（MCT） is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation（LES） of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence（FHIT） with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation（DNS） results. Compared with the LES results using the temporal approximate deconvolution model（TADM） for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.
基金supported by the National Natural Science Foundation of China (11772032, 11572025, and 51420105008)
文摘We review the previous attempts of rational subgrid-scale （SGS） modelling by employing theKolmogorov equation of filtered quantities. Aiming at explaining and solving the underlyingproblems in these models, we also introduce the recent methodological investigations for therational SGS modelling technique by defining the terms of assumption and restriction. Thesemethodological works are expected to provide instructive criterions for not only the rational SGSmodelling, but also other types of SGS modelling practices.
基金supported by the Jiangsu Shuangchuang Project(Grant No.JSSCTD202209)the National Science Foundation of the Jiangsu Higher Education Institutions of China(Grant No.22KJB130011).
文摘Vortices play a fundamental role in fluid dynamics,but mathematically defining them remains elusive.While many vortex identification methods are scalar-valued,vortices are inherently rotational,vector-based phenomena.Liutex,as a vector quantity,addresses these limitations by accurately capturing the local rotational characteristics of fluid elements while remaining independent of shear influences.This unique property makes Liutex particularly well-suited for vortex identification and the quantitative analysis of turbulent flows.This paper explores the statistical analysis of Liutex in various turbulence regimes and proposes an objective Liutex-based vortex identification method.The objective method is rooted in the statistical properties of Liutex.Furthermore,the paper investigates the performance of Liutex-based subgrid models in large eddy simulation(LES).The effectiveness of these models is evaluated by comparing their performance in different flow conditions,such as decaying homogeneous isotropic turbulence and turbulent channel flows,against conventional models.Results demonstrate that the inclusion of Liutex significantly enhances the ability of subgrid models to accurately capture flow structures.Importantly,the new model maintains the same form regardless of whether strong or weak shear is present,ensuring robustness and consistency in both vortex identification and turbulence modeling.These findings highlight the significant potential of Liutex to improve turbulence modeling in both theoretical and practical contexts,with ongoing research aimed at further refining its theoretical foundations and expanding its application in more complex flow scenarios.
基金Supported by the National Natural Science Foundation of China(42375158 and 41875121).
文摘In numerical weather prediction(NWP),the parameterization of orographic drag plays an important role in representing subgrid orographic effects.The subgrid orographic parameters are the key input to the parameterization of orographic drag.Currently,the subgrid orographic parameters in most NWP models were produced based on elevation datasets generated many years ago,with a coarse resolution and low quality.In this paper,using the latest high-quality elevation data and considering the applicable scale range of the subgrid orographic parameters,we construct the orographic parameters,including the subgrid orographic standard deviation,anisotropy,orientation,and slope,that are required as input to the orographic gravity wave drag(OGWD)parameterization.Finally,we introduce the newly constructed orographic parameters into the Yin-He Global Spectral Model(YHGSM),optimize the description of the orographic effect in the model,and improve the simulation of two typical heavy rainfall events in Beijing and Henan.
基金The National Key Research and Development Program of China under contract No.2021YFC3101602the National Natural Science Foundation of China under contract Nos 42176017 and 41976019.
文摘Mesoscale eddies play a pivotal role in deciphering the intricacies of ocean dynamics and the transport of heat,salt,and nutrients.Accurate representation of these eddies in ocean models is essential for improving model predictions.In this study,we propose a convolutional neural network(CNN)that combines data-driven techniques with physical principles to develop a robust and interpretable parameterization scheme for mesoscale eddies in ocean modeling.We use a highresolution reanalysis dataset to extract subgrid eddy momentum and then applying machine learning algorithms to identify patterns and correlations.To ensure physical consistency,we have introduced conservation of momentum constraints in our CNN parameterization scheme through soft and hard constraints.The interpretability analysis illustrate that the pre-trained CNN parameterization shows promising results in accurately solving the resolved mean velocity and effectively capturing the representation of unresolved subgrid turbulence processes.Furthermore,to validate the CNN parameterization scheme offline,we conduct simulations using the Massachusetts Institute of Technology general circulation model(MITgcm)ocean model.A series of experiments is conducted to compare the performance of the model with the CNN parameterization scheme and high-resolution simulations.The offline validation demonstrates the effectiveness of the CNN parameterization scheme in improving the representation of mesoscale eddies in the MITgcm ocean model.Incorporating the CNN parameterization scheme leads to better agreement with high-resolution simulations and a more accurate representation of the kinetic energy spectra.
基金supported by the National Natural Science Foundation of China (Grant Nos.11925201,and 11988102)the National Key R&D Program of China (Grant No.2020YFE0204200)the Xplore Prize.Numerical simulations were carried out on the TH-2A supercomputer in Guangzhou,China.
文摘We propose a hybrid quantum-classical method,the quantum-enriched large eddy simulation(QELES),for simulating turbulence.The QELES combines the large-scale motion of the large eddy simulation(LES)and the subgrid motion of the incompressible Schrodinger flow(ISF).The ISF is a possible way to be simulated on a quantum computer,and it generates subgrid scale turbu-lent structures to enrich the LES field.The enriched LES field can be further used in turbulent combustion and multi-phase flows in which the subgrid scale motion plays an important role.As a conceptual study,we perform the simulations of ISF and LES separately on a classical computer to simulate decaying homogeneous isotropic turbulence.Then,the QEI ES velocity is obtained by the time matching and the spectral blending methods.The QEL ES achieves significant improvement in predicting the energy spectrum,probaility density functions of velocity and vorticity components,and velocity structure functions,and reconstructs coherent small-scales vortices in the direct numerical simulation(DNS).On the other hand,the vortices in the QELES are less elongated and tangled than those in the DNS,and the magnitude of the third-order structure function in the QELES is less than that in the DNS,due to the diferent constitutive relations in the viscous flow and ISE.
基金Acknowledgments. The work is supported by the Natural Science Foundation of China （No. 10671154 and No. 11071184） and the National Basic Research Program （No. 2005CB321703）. It is also supported by Sichuan Science and Technology Project （No. 05GG006-006-2） and Science Research Foundation of UESTC.
文摘We consider a new subgrid eddy viscosity method based on pressure projection and extrapolated trapezoidal rule for the transient Navier-Stokes equations by using lowest equal-order pair of finite elements. The scheme stabilizes convection dominated problems and ameliorates the restrictive inf-sup compatibility stability. It has some attractive fea- tures including parameter free for the pressure stabilized term and calculations required for higher order derivatives. Moreover, it requires only the solutions of the linear system arising from an Oseen problem per time step and has second order temporal accuracy. The method achieves optimal accuracy with respect to solution regularity.
基金the National Natural Science Foundation of China (Grant No. 40705035)the Knowledge Innovation Project of Chinese Academy of Sciences (Grant Nos. KZCX2-YW-217 and KZCX2-YW-126-2)the National Basic Research Program of China (Grant No.2005CB321704)
文摘The Ensemble Kalman Filter (EnKF) is well known and widely used in land data assimilation for its high precision and simple operation. The land surface models used as the forecast operator in a land data assimilation system are usually designed to consider the model subgrid-heterogeneity and soil water thawing and freezing. To neglect their effects could lead to some errors in soil moisture assimilation. The dual EnKF method is employed in soil moisture data assimilation to build a soil moisture data as- similation framework based on the NCAR Community Land Model version 2.0 (CLM 2.0) in considera- tion of the effects of the model subgrid-heterogeneity and soil water thawing and freezing: Liquid volumetric soil moisture content in a given fraction is assimilated through the state filter process, while solid volumetric soil moisture content in the same fraction and solid/liquid volumetric soil moisture in the other fractions are optimized by the parameter filter. Preliminary experiments show that this dual EnKF-based assimilation framework can assimilate soil moisture more effectively and precisely than the usual EnKF-based assimilation framework without considering the model subgrid-scale heteroge- neity and soil water thawing and freezing. With the improvement of soil moisture simulation, the soil temperature-simulated precision can be also improved to some extent.
基金Project supported by the National Science Foundation of the Jiangsu Higher Education Institutions of China(Grant No.22KJB130011)the Supercomputing Center in Yancheng(Grant No.FW(W)20221001).
文摘The concept of vortex is crucial in both understanding and modeling of turbulence.For large eddy simulation(LES),the effect of small-scale eddies onto the large scales or the resolved flow field is modeled by subgrid stress models.Even though the rotating motions of fluids,i.e.,vortices or eddies are central in developing turbulent models,vortex identification methods are seldom used in these models.In this study,we develop a new subgrid model based on the Liutex vector,a new quantity introduced to decompose fluid motions into rigid rotation,pure shear and stretching,and thus identify vortices.The new model is then applied in a decaying homogeneous isotropic turbulence(DHIT)and a turbulent channel flow at Reynolds number Reτ=180.It is shown that the new model can predict accurate energy spectra compared with experiments in DHIT and give a well-matched velocity profile in turbulent channel flow without changing the form of the model.Future directions include improvement of the Liutex based model,for example developing anisotropic subgrid models,and its applications in various turbulent flows.
文摘In recent years, the Lattice Boltzmann Method (LBM) has developed into an alternative and promising numerical scheme for simulating fluid flows and modeling physics in fluids. In order to propose LBM for high Reynolds number fluid flow applications, a subgrid turbulence model for LBM was introduced based on standard Smagorinsky subgrid model and Lattice Bhatnagar-Gross-Krook (LBGK) model. The subgrid LBGK model was subsequently used to simulate the two-dimensional driven cavity flow at high Reynolds numbers. The simulation results including distribution of stream lines, dimensionless velocities distribution, values of stream function, as well as location of vertex center, were compared with benchmark solutions, with satisfactory agreements.
