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.展开更多
For Francis turbines,frequent operations under extremely low load conditions result in significant noise and pressure fluctuation issues.These issues may cause vibration and fatigue damage to the unit,accompanied by d...For Francis turbines,frequent operations under extremely low load conditions result in significant noise and pressure fluctuation issues.These issues may cause vibration and fatigue damage to the unit,accompanied by difficulties in connecting to the grid and reductions in the power generation efficiency of renewable energy.However,there is limited research on the relationship between pressure fluctuations and the induced noise of Francis turbines during extreme operations.In the present study,an acoustic numerical simulation based on the Ffowcs Williams-Hawkings equation and large eddy simulation is used to analyze the acoustic performances of Francis turbines.In the current study,for evaluating the acoustic characteristics under such terrible conditions,the results of variable flow rate and guide vane opening conditions are compared.Results indicated that Francis turbine noise is mostly due to pressure fluctuations brought on by rotor-stator interference and corkscrew-shaped vortices.The blade passing frequency(BPF)of 130.00 Hz and the low frequency of 0.33 f_(n)(where f_(n)denotes the rotating frequency)are the key factors affecting pressure and noise fluctuations.The influence of low frequency is reduced as the flow rate rises,whereas the influence of BPF gradually increases.Besides,the hydrodynamic noise of Francis turbines is primarily low-frequency,with discrete and broad-band features.The rotating noise with distinct peak values and the turbulence noise produced by large-scale vortices(corkscrew-shaped vortices)make up the majority of low-frequency noise.Therefore,reducing pressure fluctuations is a key strategy for lowering flow-induced noise radiation.展开更多
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.展开更多
Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, signific...Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, significantly extended design freedom and knowledge must be available in the early design stages. This paper presents a high-fidelity framework for design and optimization of the liquid swirl injectors that are widely used in aerospace propulsion and power-generation systems. The framework assembles a set of techniques, including Design Of Experiment(DOE), high-fidelity Large Eddy Simulations(LES), machine learning, Proper Orthogonal Decomposition(POD)-based Kriging surrogate modeling(emulation), inverse problem optimization, and uncertainty quantification. LES-based simulations can reveal detailed spatiotemporal evolution of flow structures and flame dynamics in a high-fidelity manner, and identify important injector design parameters according to their effects on propellant mixing, flame stabilization, and thermal protection.For a given a space of design parameters, DOE determines the number of design points to perform LES-based simulations. POD-based emulations, trained by the LES database, can effectively explore the design space and deduce an optimal group of design parameters in a turn-around time that is reduced by three orders of magnitude. The accuracy of the emulated results is validated, and the uncertainty of prediction is quantified. The proposed design methodology is expected to profoundly extend the knowledge base and reduce the cost for initial design stages.展开更多
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.展开更多
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.展开更多
Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the c...Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the characteristics of inflow turbulence,can significantly impact the quality of predictions.In this study,we examined these boundary conditions within the context of the mountainous terrain around a long-span cable-stayed bridge using a wind tunnel experiment.Various sizes of computational domains and turbulent incoming wind velocities were used in large eddy simulations.The results show that when the height of the computational domain is five times greater than the height of the terrain model,there is minimal influence from the top wall on the wind field characteristics in this complex mountainous area.Expanding the length of the wake region of the computational domain has negligible effects on the wind fields.Turbulence in the inlet boundary reduces the length of the wake region on a leeward hill with a low slope,but has less impact on the mean wind velocity of steep hills.展开更多
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 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.展开更多
The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in...The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.展开更多
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 this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross fl...In this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross flow were studied through experiments and numerical simulations.The results indicate a strong dependency of the cylinder's vibration response on the fin parameters.The results indicate that the vibration response and wake structure of the cylinder are significantly influenced by the parameters of the fins.The introduction of a finned cylinder affects both its own vibration amplitude and frequency,as well as the downstream cylinder.The amplitudes of finned cylinders Ⅰ and Ⅲ are reduced by 57.8% and 59.9%,respectively,compared to the bare cylinder.This reduction helps to restrain vibration and diminishes the amplitudes of the downstream cylinder.Although finned cylinder Ⅱ slightly decreases its own vibration,it increases the amplitude of the downstream cylinder by 13.7%.The mean drag coefficient and the root mean square of the lift coefficient of the finned cylinder are higher than those of the bare cylinder when the finned cylinder is positioned upstream.Smaller pitch and larger equivalent diameter will lead to increased drag,resulting in enhanced vortex shedding in the wake,which amplifies the vibrations of the cylinder in that wake.The downstream of finned cylinder Ⅱ has the widest wake and higher vortex strength,and the dynamic load and vibration of the downstream cylinder are increased.The vortex intensity decays faster in the wake of finned cylinder Ⅲ,and the vibration of the downstream cylinder is weaker.展开更多
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.展开更多
Fractal theory provides a new strategy for equipment design.In this work,we propose a novel H-like fractal(HLF)impeller to improve the uniformity of the distribution of hydrodynamics in stirred tanks.The impellers are...Fractal theory provides a new strategy for equipment design.In this work,we propose a novel H-like fractal(HLF)impeller to improve the uniformity of the distribution of hydrodynamics in stirred tanks.The impellers are constructed by replacing two vertical blades or four legs with two or four H-like subblades by fractal iterations,respectively.Flow characteristics including velocity and turbulent kinetic energy(TKE)distributions,vortices,power number,are predicted by large eddy simulation.Compared with Rushton turbine(RT)impeller when H/T=1(or dual RTs when H/T=1.5,triple RTs when H/T=2),the HLF impeller can produce a flow field with more uniform distributions of larger velocities and TKE level.The impeller with more fractal iteration times can further improve the distribution uniformity of hydrodynamics in the case of high H/T.Power analysis shows that this is mainly due to the improved energy utilization efficiency by the fractal structure design.展开更多
The cutaneous ridges on dolphin skin have long been believed to effectively reduce friction drag, thereby contributing to overall drag reduction. However, since these skin ridges are oriented perpendicular to the swim...The cutaneous ridges on dolphin skin have long been believed to effectively reduce friction drag, thereby contributing to overall drag reduction. However, since these skin ridges are oriented perpendicular to the swimming direction, they also generate additional pressure drag, raising questions about the impact of the shape-induced pressure forces on swimming. Inspired by the microvibrations observed on dolphin skin, we hypothesize that the microstructure on dolphin skin is not static but dynamically oscillates in the form of Longitudinal Micro-Ultrasonic Waves (LMUWs). To explore this, we carried out a series of Computational Fluid Dynamics (CFD) simulations based on Large Eddy Simulation (LES) model to investigate the impact of pressure drag on the total drag acting on an oscillating skin surface under realistic turbulent flow conditions. The results indicate that the dynamic skin oscillations induce a new dynamic Stokes boundary layer, which has the potential to convert pressure drag into a negative force, thereby reducing total drag under the influence of traveling LMUW excitations. Furthermore, a relative velocity ξ, defined as the difference between the wave speed c and the external flow speed U, is introduced to evaluate the drag-reduction effect dominated by pressure. The findings reveal that pressure drag remains negative when ξ > 0. As ξ increases, the thrust effect induced by negative pressure becomes increasingly significant, ultimately counteracting friction drag and eliminating total drag. This pressure-dominated drag reduction mechanism thus demonstrates a novel strategy for the drag reduction technology and the potential of unveiling the mysteries behind dolphin swimming.展开更多
Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealize...Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealized numerical experiment using the large eddy simulation method was conducted to analyze the three-dimensional flow field structure and the mechanism for dense water subsidence.Specifically,a negative salt flux is set at the sea surface,in which salt flux enters the sea surface to simulate the icing and salting-out phenomena that occur at high latitudes.Results show that the mean-state 3D flow field of dense water subsidence exhibits a hollow conical distribution.The horizontal flow field is characterized by a cyclonic vortex that driven primarily by the pressure gradient and influenced by the Coriolis effect.Moreover,the inverse vertical pressure gradient generated by this vortex inhibits the sinking of the plume,leading to its off-axis deflection and the development of an anticyclonic precession.In addition,the impact of rotation on the structure of a sinking plume within a stratified environment is discussed.Both horizontal vortex intensity and cone angle of the hollow cone flow field are increased with increasing rotation rate,resulting in a decrease in the plume’s maximum sinking depth.Variances in rotation direction cause the horizontal vortex and sinking plumes of dense water in the northern and southern hemispheres to rotate in opposite directions.展开更多
基金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.
基金supported by the National Natural Science Foundation of China(Grant No.52179087).
文摘For Francis turbines,frequent operations under extremely low load conditions result in significant noise and pressure fluctuation issues.These issues may cause vibration and fatigue damage to the unit,accompanied by difficulties in connecting to the grid and reductions in the power generation efficiency of renewable energy.However,there is limited research on the relationship between pressure fluctuations and the induced noise of Francis turbines during extreme operations.In the present study,an acoustic numerical simulation based on the Ffowcs Williams-Hawkings equation and large eddy simulation is used to analyze the acoustic performances of Francis turbines.In the current study,for evaluating the acoustic characteristics under such terrible conditions,the results of variable flow rate and guide vane opening conditions are compared.Results indicated that Francis turbine noise is mostly due to pressure fluctuations brought on by rotor-stator interference and corkscrew-shaped vortices.The blade passing frequency(BPF)of 130.00 Hz and the low frequency of 0.33 f_(n)(where f_(n)denotes the rotating frequency)are the key factors affecting pressure and noise fluctuations.The influence of low frequency is reduced as the flow rate rises,whereas the influence of BPF gradually increases.Besides,the hydrodynamic noise of Francis turbines is primarily low-frequency,with discrete and broad-band features.The rotating noise with distinct peak values and the turbulence noise produced by large-scale vortices(corkscrew-shaped vortices)make up the majority of low-frequency noise.Therefore,reducing pressure fluctuations is a key strategy for lowering flow-induced noise radiation.
文摘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.
基金sponsored by the William RT Oakes Endowment of the Georgia Institute of Technology
文摘Engineering design is undergoing a paradigm shift from design for performance to design for affordability, operability, and durability, seeking multi-objective optimization. To facilitate this transformation, significantly extended design freedom and knowledge must be available in the early design stages. This paper presents a high-fidelity framework for design and optimization of the liquid swirl injectors that are widely used in aerospace propulsion and power-generation systems. The framework assembles a set of techniques, including Design Of Experiment(DOE), high-fidelity Large Eddy Simulations(LES), machine learning, Proper Orthogonal Decomposition(POD)-based Kriging surrogate modeling(emulation), inverse problem optimization, and uncertainty quantification. LES-based simulations can reveal detailed spatiotemporal evolution of flow structures and flame dynamics in a high-fidelity manner, and identify important injector design parameters according to their effects on propellant mixing, flame stabilization, and thermal protection.For a given a space of design parameters, DOE determines the number of design points to perform LES-based simulations. POD-based emulations, trained by the LES database, can effectively explore the design space and deduce an optimal group of design parameters in a turn-around time that is reduced by three orders of magnitude. The accuracy of the emulated results is validated, and the uncertainty of prediction is quantified. The proposed design methodology is expected to profoundly extend the knowledge base and reduce the cost for initial design stages.
基金supported by 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 National Natural Science Foundation of China(Grant No.50979059)
文摘Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.
基金supported by the National Natural Science Foundation of China(Nos.51925808 and 52178516)the Natural Science Foundation of Hunan Province(Nos.2020JJ5745 and 2023JJ20073),China.
文摘Large eddy simulations generally are used to predict 3D wind field characteristics in complex mountainous areas.Certain simulation boundary conditions,such as the height and length of the computational domain or the characteristics of inflow turbulence,can significantly impact the quality of predictions.In this study,we examined these boundary conditions within the context of the mountainous terrain around a long-span cable-stayed bridge using a wind tunnel experiment.Various sizes of computational domains and turbulent incoming wind velocities were used in large eddy simulations.The results show that when the height of the computational domain is five times greater than the height of the terrain model,there is minimal influence from the top wall on the wind field characteristics in this complex mountainous area.Expanding the length of the wake region of the computational domain has negligible effects on the wind fields.Turbulence in the inlet boundary reduces the length of the wake region on a leeward hill with a low slope,but has less impact on the mean wind velocity of steep hills.
文摘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.
基金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(10772082)~~
文摘The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.
基金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.
基金financially supported by the National Natural Science Foundation of China(22478286)。
文摘In this study,four types of spiral fins with varying parameters were mounted on an upstream cylinder,and the effects of spiral fins on the vibration response of heat exchange tubes and the vortex structure in cross flow were studied through experiments and numerical simulations.The results indicate a strong dependency of the cylinder's vibration response on the fin parameters.The results indicate that the vibration response and wake structure of the cylinder are significantly influenced by the parameters of the fins.The introduction of a finned cylinder affects both its own vibration amplitude and frequency,as well as the downstream cylinder.The amplitudes of finned cylinders Ⅰ and Ⅲ are reduced by 57.8% and 59.9%,respectively,compared to the bare cylinder.This reduction helps to restrain vibration and diminishes the amplitudes of the downstream cylinder.Although finned cylinder Ⅱ slightly decreases its own vibration,it increases the amplitude of the downstream cylinder by 13.7%.The mean drag coefficient and the root mean square of the lift coefficient of the finned cylinder are higher than those of the bare cylinder when the finned cylinder is positioned upstream.Smaller pitch and larger equivalent diameter will lead to increased drag,resulting in enhanced vortex shedding in the wake,which amplifies the vibrations of the cylinder in that wake.The downstream of finned cylinder Ⅱ has the widest wake and higher vortex strength,and the dynamic load and vibration of the downstream cylinder are increased.The vortex intensity decays faster in the wake of finned cylinder Ⅲ,and the vibration of the downstream cylinder is weaker.
基金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.
基金the financial support from the National Natural Science Foundation of China(22078058)。
文摘Fractal theory provides a new strategy for equipment design.In this work,we propose a novel H-like fractal(HLF)impeller to improve the uniformity of the distribution of hydrodynamics in stirred tanks.The impellers are constructed by replacing two vertical blades or four legs with two or four H-like subblades by fractal iterations,respectively.Flow characteristics including velocity and turbulent kinetic energy(TKE)distributions,vortices,power number,are predicted by large eddy simulation.Compared with Rushton turbine(RT)impeller when H/T=1(or dual RTs when H/T=1.5,triple RTs when H/T=2),the HLF impeller can produce a flow field with more uniform distributions of larger velocities and TKE level.The impeller with more fractal iteration times can further improve the distribution uniformity of hydrodynamics in the case of high H/T.Power analysis shows that this is mainly due to the improved energy utilization efficiency by the fractal structure design.
基金supported by the Japan Society for the Promotion of Science(JSPS KAKENHI No.23H01373).
文摘The cutaneous ridges on dolphin skin have long been believed to effectively reduce friction drag, thereby contributing to overall drag reduction. However, since these skin ridges are oriented perpendicular to the swimming direction, they also generate additional pressure drag, raising questions about the impact of the shape-induced pressure forces on swimming. Inspired by the microvibrations observed on dolphin skin, we hypothesize that the microstructure on dolphin skin is not static but dynamically oscillates in the form of Longitudinal Micro-Ultrasonic Waves (LMUWs). To explore this, we carried out a series of Computational Fluid Dynamics (CFD) simulations based on Large Eddy Simulation (LES) model to investigate the impact of pressure drag on the total drag acting on an oscillating skin surface under realistic turbulent flow conditions. The results indicate that the dynamic skin oscillations induce a new dynamic Stokes boundary layer, which has the potential to convert pressure drag into a negative force, thereby reducing total drag under the influence of traveling LMUW excitations. Furthermore, a relative velocity ξ, defined as the difference between the wave speed c and the external flow speed U, is introduced to evaluate the drag-reduction effect dominated by pressure. The findings reveal that pressure drag remains negative when ξ > 0. As ξ increases, the thrust effect induced by negative pressure becomes increasingly significant, ultimately counteracting friction drag and eliminating total drag. This pressure-dominated drag reduction mechanism thus demonstrates a novel strategy for the drag reduction technology and the potential of unveiling the mysteries behind dolphin swimming.
基金Supported by the National Natural Science Foundation of China(Nos.42250710152,42192562)the Southern Laboratory of Ocean Science and Engineering(Guangdong Zhuhai)(No.SML 2020 SP 007)。
文摘Convectively unstable processes caused by dense water subsidence are common occurrences in high-latitude oceanic regions,and significantly modulate mass and heat transport and mixing processes in the ocean.An idealized numerical experiment using the large eddy simulation method was conducted to analyze the three-dimensional flow field structure and the mechanism for dense water subsidence.Specifically,a negative salt flux is set at the sea surface,in which salt flux enters the sea surface to simulate the icing and salting-out phenomena that occur at high latitudes.Results show that the mean-state 3D flow field of dense water subsidence exhibits a hollow conical distribution.The horizontal flow field is characterized by a cyclonic vortex that driven primarily by the pressure gradient and influenced by the Coriolis effect.Moreover,the inverse vertical pressure gradient generated by this vortex inhibits the sinking of the plume,leading to its off-axis deflection and the development of an anticyclonic precession.In addition,the impact of rotation on the structure of a sinking plume within a stratified environment is discussed.Both horizontal vortex intensity and cone angle of the hollow cone flow field are increased with increasing rotation rate,resulting in a decrease in the plume’s maximum sinking depth.Variances in rotation direction cause the horizontal vortex and sinking plumes of dense water in the northern and southern hemispheres to rotate in opposite directions.
