An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional ...An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional anisotropic Gaussian kernel,this method effectively addresses the severe aerodynamic load fluctuations commonly associated with traditional Virtual Blade Method(VBM)due to turbulent flow around blade elements.To manage the issues of regional shape and grid cell quantity variations caused by virtual blade sweeping,a universal hybrid grid generation strategy is established without body-fitted and disk interpolation grids,which enhances the computational stability at both blade elements and blade edges.Aerodynamic numerical simulations of helicopter are performed using this method,focusing on rotor/fuselage interaction dominated by rotor wake motion and fuselage blockage effects,Blade-Vortex Interaction(BVI)induced by tip vortices,and maneuvering flights involving collective pitch ramp increases.The results indicate that the advanced ASM demonstrates reliability and robustness in the simulation of complex unsteady flow fields around helicopter.Under similar computational accuracy,the advanced ASM improves computational efficiency by nearly 40 times compared to the oversetgrid-based full Blade-Resolved(B-R)method,and by 6 times compared to the VBM.It shows significant advantages when applied to complex full-aircraft interaction and maneuvering flight conditions that require substantial computational resources.展开更多
The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under th...The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.展开更多
Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the ...Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.展开更多
Computational Fluid Dynamics(CFD)simulation is one of the important methods to study the performance and influencing factors of turbine flow sensors.According to the working characteristics of the turbine flow sensor,...Computational Fluid Dynamics(CFD)simulation is one of the important methods to study the performance and influencing factors of turbine flow sensors.According to the working characteristics of the turbine flow sensor,the passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed in this paper.The reliability of the simulation method is verified by steady-state experiments and unsteady-state experiments.The results show that the trend of meter factor with flow rate acquired from the simulation is close to the experimental results,and the deviation between the simulation result and the experiment result is low with a maximum deviation of 2.88%.In the unsteady simulation study,the impeller speed changes with the inlet velocity of the turbine flow sensor,which has a good follow-up.The passive simulation method can be used to predict the dynamic performance of the turbine flow sensor.展开更多
Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate f...Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.展开更多
The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the...The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.展开更多
The distribution of remaining oil is often described qualitatively. The remaining oil distributed in the whole reservoir is calculated according to the characteristics of the space distribution of the saturation of re...The distribution of remaining oil is often described qualitatively. The remaining oil distributed in the whole reservoir is calculated according to the characteristics of the space distribution of the saturation of remaining oil. Logging data are required to accomplish this. However, many such projects cannot be completed. Since the old study of remaining oil distribution could not be quantified efficiently, the ＂dynamic two-step method＂ is presented. Firstly, the water cut of every flow unit in one well at one time is calculated according to the comprehensive water cut of a single well at one time. Secondly, the remaining oil saturation of the flow unit of the well at one time is calculated based on the water cut of the flow unit at a given time. The results show that ＂dynamic two-step method＂ has characteristics of simplicity and convenience, and is especially suitable for the study of remaining oil distribution at high water-cut stage. The distribution of remaining oil presented banding and potato form, remaining oil was relatively concentrated in faultage neighborhood and imperfect well netting position, and the net thickness of the place was great. This proposal can provide an effective way to forecast remaining oil distribution and enhance oil recovery, especially applied at the high water-cut stage.展开更多
The present paper reviews the recent developments of a high⁃order⁃spectral method(HOS)and the combination with computational fluid dynamics(CFD)method for wave⁃structure interactions.As the numerical simulations of wa...The present paper reviews the recent developments of a high⁃order⁃spectral method(HOS)and the combination with computational fluid dynamics(CFD)method for wave⁃structure interactions.As the numerical simulations of wave⁃structure interaction require efficiency and accuracy,as well as the ability in calculating in open sea states,the HOS method has its strength in both generating extreme waves in open seas and fast convergence in simulations,while computational fluid dynamics(CFD)method has its advantages in simulating violent wave⁃structure interactions.This paper provides the new thoughts for fast and accurate simulations,as well as the future work on innovations in fine fluid field of numerical simulations.展开更多
This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology co...This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology considers the variables of debris flow depth and velocity and introduces the intensity index IDV(IDV = DV) to evaluate three different levels of debris flow impact intensity. The primary focus of this study is to investigate the dynamic response and failure mechanism of RC-frame buildings under debris flow impact, including structural failure patterns, impact force and column displacement. The results show that under a highintensity impact, a gradual collapse process of the RCframe building can be observed, and the damage mode of the frame column reflects shear failure or plastic hinge failure mechanism. First, the longitudinal infill walls are damaged owing to their low out-of-plane flexural capacity;the critical failure intensity index IDV value is approximately 7.5 m2/s. The structure cannot withstand debris flows with an intensity index IDV greater than 16 m2/s, and it is recommended that the peak impact force should not exceed 2100 k N. The impact damage ability of debris flow on buildings mostly originates from the impact force of the frontal debris flow, with the impact force of the debris flow body being approximately 42% lower than that of the debris flow head. Finally, a five-level classification system for evaluating the damage status of buildings is proposed based on the numerical simulation and investigation results of the disaster site.展开更多
The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, ...The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.展开更多
The main purpose of this paper is to develop a gridless method for unsteady flow simulation. A quadrantal point infilling strategy is developed to generate point and combine clouds of points automatically. A point-mov...The main purpose of this paper is to develop a gridless method for unsteady flow simulation. A quadrantal point infilling strategy is developed to generate point and combine clouds of points automatically. A point-moving algorithm is introduced to ensure the clouds of points following the movements of bodyboundaries. A dual time method for solving the two-dimenslonal Euler equations in Arbitrary Lagrangian-Eulerian （ALE） formulation is presented. Dual time method allows the real-time step to be chosen on the basis of accuracy rather than stability. It also permits the acceleration techniques, which are commonly used to speed up steady flow calculations, to be used when marching the equations in pseudo time. The spatial derivatives, which are used to estimating the inviscid flux, are directly approximated by using local least-squares curve method. An explicit multistage Runge-Kutta algorithm is used to advance the flow equations in pseudo time. In order to accelerate the solution to convergence, local time stepping technique and residual averaging are employed. The results of NACA0012 airfoil in transonic steady flow are presented to verify the accuracy of the present spatial discretization method. Finally, two AGARD standard test cases in which NACA0012 airfoil and NACA64A010 airfoil oscillate in transonic flow are simulated. The computational results are compared with the experimental data to demonstrate the validity and practicality of the presented method.展开更多
A number of studies have found that abnormal changes of dynamic derivatives occurred at very low reduced frequencies, but its inducement mechanism is not very clear. This paper has researched the abnormal changes and ...A number of studies have found that abnormal changes of dynamic derivatives occurred at very low reduced frequencies, but its inducement mechanism is not very clear. This paper has researched the abnormal changes and analyzed the influence on some parameters by solving the unsteady flow around forced oscillation airfoils based on Navier-Stokes equations. Results indicate that when the reduced frequency approaches to zero, the dynamic derivatives obtained by the numerical method will diverge. We have also proven it in theory that this phenomenon is not physical but completely caused by numerical singularity. Furthermore, the abnormal phenomenon can be effectively mitigated by using the time spectral method to solve the aerodynamic forces and the integral method to obtain the dynamic derivatives. When the reduced frequency is in the range of 0.001–0.01, the dynamic derivative maintains nearly unchanged for the whole speed region. This study can provide a reference for the reasonable choice of the reduced frequency in calculations and experiments of dynamic derivatives.展开更多
A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid ...A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid dynamics(CFD) approach, in which the unsteady, noncompressible Reynolds Averaged Navier-Stokes(RANS) method is used. The main focus of the study is to analyze the tower shadow effect on the aerodynamic performance of the wind turbine under different inlet flow conditions. Subsequently, the finite element model is established by considering fluid/structure interactions to study the structural stress, displacement, strain distributions and flow field information of the structure under the uniform wind speed. Finally, the fluid-structure interaction model is established by considering turbulent wind and the tower shadow effect. The variation rules of the dynamic response of the one-way and two-way fluid-structure interaction(FSI) models under different wind speeds are analyzed, and the numerical calculation results are compared with those of the centralized mass model. The results show that the tower shadow effect and structural deformation are the main factors affecting the aerodynamic load fluctuation of the wind turbine, which in turn affects the aerodynamic performance and structural stability of the blades. The structural dynamic response of the coupled model shows significant similarity, while the structural displacement response of the former exhibits less fluctuation compared with the conventional centralized mass model. The one-way fluid-structure interaction(FSI)model shows a higher frequency of stress-strain and displacement oscillations on the blade compared with the two-way FSI model.展开更多
A passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed according to the working principle to study the dynamic characteristics of the turbine flow sensors.This simulat...A passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed according to the working principle to study the dynamic characteristics of the turbine flow sensors.This simulation method controls the six degrees of freedom of the impeller using the user-defined functions(UDF)program so that it can only rotate under the impact of fluid.The impeller speed can be calculated in real-time,and the inlet speed can be set with time to obtain the dynamic performance of the turbine flow sensors.Based on this simulation method,three turbine flow sensors with different diameters were simulated,and the reliability of the simulation method was verified by both steady-state and unsteady-state experiments.The results show that the trend of meter factor with flow rate acquired from the simulation is close to the experimental results.The deviation between the simulation and experiment results is low,with a maximum deviation of 2.88%.In the unsteady simulation study,the impeller speed changed with the inlet velocity of the turbine flow sensor,showing good tracking performance.The passive simulation method can be used to predict the dynamic performance of the turbine flow sensor.展开更多
The knowledge of bubble profiles in gas-liquid two-phase flows is crucial for analyzing the kinetic processes such as heat and mass transfer, and this knowledge is contained in field data obtained by surface-resolved ...The knowledge of bubble profiles in gas-liquid two-phase flows is crucial for analyzing the kinetic processes such as heat and mass transfer, and this knowledge is contained in field data obtained by surface-resolved computational fluid dynamics (CFD) simulations. To obtain this information, an efficient bubble profile reconstruction method based on an improved agglomerative hierarchical clustering (AHC) algorithm is proposed in this paper. The reconstruction method is featured by the implementations of a binary space division preprocessing, which aims to reduce the computational complexity, an adaptive linkage criterion, which guarantees the applicability of the AHC algorithm when dealing with datasets involving either non-uniform or distorted grids, and a stepwise execution strategy, which enables the separation of attached bubbles. To illustrate and verify this method, it was applied to dealing with 3 datasets, 2 of them with pre-specified spherical bubbles and the other obtained by a surface-resolved CFD simulation. Application results indicate that the proposed method is effective even when the data include some non-uniform and distortion.展开更多
Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with...Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with the theory of heat-transfer and hydrodynamics by means of 3-dimensional finite element method. From simulation results, it is found that the variation of temperature and velocityis mainly influenced by the casting temperature of aluminum, rotating speed of Castex wheel and flow of cooling water. Among theseinfluencing factors, the casting temperature distributes most to the length of liquid phase metal. Moreover, the faster the metal solidifies,the higher the metal there moves with the overall trend of descending from the bottom of the wheel to the shoe wall as well as from sidewalls to the center of wheel groove. In comparison with the practical value, the simulation is reliable.展开更多
Liquid-liquid dispersion is often performed in stirred tanks,which are valued for their ease of operation,high droplet generation rate and effective droplet dispersion.Many relevant simulations use the Eulerian-Euleri...Liquid-liquid dispersion is often performed in stirred tanks,which are valued for their ease of operation,high droplet generation rate and effective droplet dispersion.Many relevant simulations use the Eulerian-Eulerian method,combining population balance equations with statistical models to forecast droplet breakage.Conversely,the Eulerian-Lagrangian(E-L)method provides precise tracking of individual droplets,which is crucial for simulating dispersion processes.However,E-L simulation faces challenges in integrating droplet breakage effectively.To address this issue,our research introduces a probabilistic approach for droplet breakages.It assumes that a longer time increases the likelihood of breakup;a droplet breaks if the calculated probability exceeds a random value from 0 to 1.Consequently,the simulated breakage frequency becomes independent of the Lagrangian time step.The Sauter mean diameter and droplet size distribution can be accurately predicted by this probabilistic approach.By closely monitoring droplet motion,we reveal the complexity of droplet trajectories and the detailed patterns of circulation in stirred tanks.These insights contribute to a deeper understanding of liquidliquid dispersion dynamics.展开更多
1.Introduction Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)is a powerful tool for simulating dense gas-solid reacting flows,which is essential in combustion,metallurgy,and waste management.Traditional...1.Introduction Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)is a powerful tool for simulating dense gas-solid reacting flows,which is essential in combustion,metallurgy,and waste management.Traditional methods face challenges in CFD-DEM modeling of dense gas-solid flows due to multi-scale characteristics,limiting resolution and creating simulation bottlenecks.By integrating fluid dynamics and particle behavior,it optimizes industrial processes.This review highlights advancements,applications,and challenges,emphasizing its role in sustainable engineering.展开更多
基金co-supported by the Foundation of the State Key Laboratory of Aerodynamics(No.RAL202203)the National Key Laboratory Foundation(No.6142202202)the China Postdoctoral Science Foundation(No.2024M754133)。
文摘An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional anisotropic Gaussian kernel,this method effectively addresses the severe aerodynamic load fluctuations commonly associated with traditional Virtual Blade Method(VBM)due to turbulent flow around blade elements.To manage the issues of regional shape and grid cell quantity variations caused by virtual blade sweeping,a universal hybrid grid generation strategy is established without body-fitted and disk interpolation grids,which enhances the computational stability at both blade elements and blade edges.Aerodynamic numerical simulations of helicopter are performed using this method,focusing on rotor/fuselage interaction dominated by rotor wake motion and fuselage blockage effects,Blade-Vortex Interaction(BVI)induced by tip vortices,and maneuvering flights involving collective pitch ramp increases.The results indicate that the advanced ASM demonstrates reliability and robustness in the simulation of complex unsteady flow fields around helicopter.Under similar computational accuracy,the advanced ASM improves computational efficiency by nearly 40 times compared to the oversetgrid-based full Blade-Resolved(B-R)method,and by 6 times compared to the VBM.It shows significant advantages when applied to complex full-aircraft interaction and maneuvering flight conditions that require substantial computational resources.
基金financially supported by the National Natural Science Foundation of China (No. 42102149)the Fundamental Research Funds for the Central Universities (No. 2462021YXZZ005)。
文摘The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.
基金Project supported by the National Natural Science Foundation of China(No.52109068)the Water Conservancy Technology Project of Jiangsu Province of China(No.2022060)。
文摘Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.
基金supported by National Natural Science Foundation of China,China(62173122)Hebei Key Project of Natural Science Foundation(F2021201031)Beijing-Tianjin-Hebei Collaborative Innovation Community Construction Project(20540301D)
文摘Computational Fluid Dynamics(CFD)simulation is one of the important methods to study the performance and influencing factors of turbine flow sensors.According to the working characteristics of the turbine flow sensor,the passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed in this paper.The reliability of the simulation method is verified by steady-state experiments and unsteady-state experiments.The results show that the trend of meter factor with flow rate acquired from the simulation is close to the experimental results,and the deviation between the simulation result and the experiment result is low with a maximum deviation of 2.88%.In the unsteady simulation study,the impeller speed changes with the inlet velocity of the turbine flow sensor,which has a good follow-up.The passive simulation method can be used to predict the dynamic performance of the turbine flow sensor.
基金Aeronautical Science Foundation of China (99A52007)
文摘Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.
基金supported by Scientific Research Foundation for Returned Scholars,Ministry of Education of China
文摘The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.
文摘The distribution of remaining oil is often described qualitatively. The remaining oil distributed in the whole reservoir is calculated according to the characteristics of the space distribution of the saturation of remaining oil. Logging data are required to accomplish this. However, many such projects cannot be completed. Since the old study of remaining oil distribution could not be quantified efficiently, the ＂dynamic two-step method＂ is presented. Firstly, the water cut of every flow unit in one well at one time is calculated according to the comprehensive water cut of a single well at one time. Secondly, the remaining oil saturation of the flow unit of the well at one time is calculated based on the water cut of the flow unit at a given time. The results show that ＂dynamic two-step method＂ has characteristics of simplicity and convenience, and is especially suitable for the study of remaining oil distribution at high water-cut stage. The distribution of remaining oil presented banding and potato form, remaining oil was relatively concentrated in faultage neighborhood and imperfect well netting position, and the net thickness of the place was great. This proposal can provide an effective way to forecast remaining oil distribution and enhance oil recovery, especially applied at the high water-cut stage.
基金National Natural Science Foundation of China(Grant No.51879159)the National Key Research and Development Program of China(Grant Nos.2019YFB1704200 and 2019YFC0312400)+2 种基金the Chang Jiang Scholars Program(Grant No.T2014099)the Shanghai Excellent Academic Leaders Program(Grant No.17XD1402300)the Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(Grant No.2016-23/09).
文摘The present paper reviews the recent developments of a high⁃order⁃spectral method(HOS)and the combination with computational fluid dynamics(CFD)method for wave⁃structure interactions.As the numerical simulations of wave⁃structure interaction require efficiency and accuracy,as well as the ability in calculating in open sea states,the HOS method has its strength in both generating extreme waves in open seas and fast convergence in simulations,while computational fluid dynamics(CFD)method has its advantages in simulating violent wave⁃structure interactions.This paper provides the new thoughts for fast and accurate simulations,as well as the future work on innovations in fine fluid field of numerical simulations.
基金supported by the National Natural Science Foundation of China (Grant No. 41877524, No. 42172320, No. 41971214)。
文摘This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology considers the variables of debris flow depth and velocity and introduces the intensity index IDV(IDV = DV) to evaluate three different levels of debris flow impact intensity. The primary focus of this study is to investigate the dynamic response and failure mechanism of RC-frame buildings under debris flow impact, including structural failure patterns, impact force and column displacement. The results show that under a highintensity impact, a gradual collapse process of the RCframe building can be observed, and the damage mode of the frame column reflects shear failure or plastic hinge failure mechanism. First, the longitudinal infill walls are damaged owing to their low out-of-plane flexural capacity;the critical failure intensity index IDV value is approximately 7.5 m2/s. The structure cannot withstand debris flows with an intensity index IDV greater than 16 m2/s, and it is recommended that the peak impact force should not exceed 2100 k N. The impact damage ability of debris flow on buildings mostly originates from the impact force of the frontal debris flow, with the impact force of the debris flow body being approximately 42% lower than that of the debris flow head. Finally, a five-level classification system for evaluating the damage status of buildings is proposed based on the numerical simulation and investigation results of the disaster site.
基金Supported by the Aeronautical Science Foundation of China(20061453020)Foundation for Basic Research of Northwestern Polytechnical University(03)~~
文摘The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.
文摘The main purpose of this paper is to develop a gridless method for unsteady flow simulation. A quadrantal point infilling strategy is developed to generate point and combine clouds of points automatically. A point-moving algorithm is introduced to ensure the clouds of points following the movements of bodyboundaries. A dual time method for solving the two-dimenslonal Euler equations in Arbitrary Lagrangian-Eulerian （ALE） formulation is presented. Dual time method allows the real-time step to be chosen on the basis of accuracy rather than stability. It also permits the acceleration techniques, which are commonly used to speed up steady flow calculations, to be used when marching the equations in pseudo time. The spatial derivatives, which are used to estimating the inviscid flux, are directly approximated by using local least-squares curve method. An explicit multistage Runge-Kutta algorithm is used to advance the flow equations in pseudo time. In order to accelerate the solution to convergence, local time stepping technique and residual averaging are employed. The results of NACA0012 airfoil in transonic steady flow are presented to verify the accuracy of the present spatial discretization method. Finally, two AGARD standard test cases in which NACA0012 airfoil and NACA64A010 airfoil oscillate in transonic flow are simulated. The computational results are compared with the experimental data to demonstrate the validity and practicality of the presented method.
基金supported by the National Science Foundation for Distinguished Young Scholars of China (No.11622220)the Programme of Introducing Talents of Discipline to Universities (No.B17037)
文摘A number of studies have found that abnormal changes of dynamic derivatives occurred at very low reduced frequencies, but its inducement mechanism is not very clear. This paper has researched the abnormal changes and analyzed the influence on some parameters by solving the unsteady flow around forced oscillation airfoils based on Navier-Stokes equations. Results indicate that when the reduced frequency approaches to zero, the dynamic derivatives obtained by the numerical method will diverge. We have also proven it in theory that this phenomenon is not physical but completely caused by numerical singularity. Furthermore, the abnormal phenomenon can be effectively mitigated by using the time spectral method to solve the aerodynamic forces and the integral method to obtain the dynamic derivatives. When the reduced frequency is in the range of 0.001–0.01, the dynamic derivative maintains nearly unchanged for the whole speed region. This study can provide a reference for the reasonable choice of the reduced frequency in calculations and experiments of dynamic derivatives.
基金supported by the National Natural Science Foundation of China(Grant No.52078010)Beijing Natural Science Foundation(Grant No.JQ19029).
文摘A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid dynamics(CFD) approach, in which the unsteady, noncompressible Reynolds Averaged Navier-Stokes(RANS) method is used. The main focus of the study is to analyze the tower shadow effect on the aerodynamic performance of the wind turbine under different inlet flow conditions. Subsequently, the finite element model is established by considering fluid/structure interactions to study the structural stress, displacement, strain distributions and flow field information of the structure under the uniform wind speed. Finally, the fluid-structure interaction model is established by considering turbulent wind and the tower shadow effect. The variation rules of the dynamic response of the one-way and two-way fluid-structure interaction(FSI) models under different wind speeds are analyzed, and the numerical calculation results are compared with those of the centralized mass model. The results show that the tower shadow effect and structural deformation are the main factors affecting the aerodynamic load fluctuation of the wind turbine, which in turn affects the aerodynamic performance and structural stability of the blades. The structural dynamic response of the coupled model shows significant similarity, while the structural displacement response of the former exhibits less fluctuation compared with the conventional centralized mass model. The one-way fluid-structure interaction(FSI)model shows a higher frequency of stress-strain and displacement oscillations on the blade compared with the two-way FSI model.
基金The National Natural Science Foundation of China(No.62173122)the Hebei Key Project of Natural Science Foundation(No.F2021201031)。
文摘A passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed according to the working principle to study the dynamic characteristics of the turbine flow sensors.This simulation method controls the six degrees of freedom of the impeller using the user-defined functions(UDF)program so that it can only rotate under the impact of fluid.The impeller speed can be calculated in real-time,and the inlet speed can be set with time to obtain the dynamic performance of the turbine flow sensors.Based on this simulation method,three turbine flow sensors with different diameters were simulated,and the reliability of the simulation method was verified by both steady-state and unsteady-state experiments.The results show that the trend of meter factor with flow rate acquired from the simulation is close to the experimental results.The deviation between the simulation and experiment results is low,with a maximum deviation of 2.88%.In the unsteady simulation study,the impeller speed changed with the inlet velocity of the turbine flow sensor,showing good tracking performance.The passive simulation method can be used to predict the dynamic performance of the turbine flow sensor.
基金Projects(51634010,51676211) supported by the National Natural Science Foundation of ChinaProject(2017SK2253) supported by the Key Research and Development Program of Hunan Province,China
文摘The knowledge of bubble profiles in gas-liquid two-phase flows is crucial for analyzing the kinetic processes such as heat and mass transfer, and this knowledge is contained in field data obtained by surface-resolved computational fluid dynamics (CFD) simulations. To obtain this information, an efficient bubble profile reconstruction method based on an improved agglomerative hierarchical clustering (AHC) algorithm is proposed in this paper. The reconstruction method is featured by the implementations of a binary space division preprocessing, which aims to reduce the computational complexity, an adaptive linkage criterion, which guarantees the applicability of the AHC algorithm when dealing with datasets involving either non-uniform or distorted grids, and a stepwise execution strategy, which enables the separation of attached bubbles. To illustrate and verify this method, it was applied to dealing with 3 datasets, 2 of them with pre-specified spherical bubbles and the other obtained by a surface-resolved CFD simulation. Application results indicate that the proposed method is effective even when the data include some non-uniform and distortion.
文摘Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with the theory of heat-transfer and hydrodynamics by means of 3-dimensional finite element method. From simulation results, it is found that the variation of temperature and velocityis mainly influenced by the casting temperature of aluminum, rotating speed of Castex wheel and flow of cooling water. Among theseinfluencing factors, the casting temperature distributes most to the length of liquid phase metal. Moreover, the faster the metal solidifies,the higher the metal there moves with the overall trend of descending from the bottom of the wheel to the shoe wall as well as from sidewalls to the center of wheel groove. In comparison with the practical value, the simulation is reliable.
基金support from the National Key Research and Development Program of China,China(2023YFE0106600)the National Natural Science Foundation of China,China(22421003,22178354,21925805)funding from FFG(Austria)under project“ABATE”(903872).
文摘Liquid-liquid dispersion is often performed in stirred tanks,which are valued for their ease of operation,high droplet generation rate and effective droplet dispersion.Many relevant simulations use the Eulerian-Eulerian method,combining population balance equations with statistical models to forecast droplet breakage.Conversely,the Eulerian-Lagrangian(E-L)method provides precise tracking of individual droplets,which is crucial for simulating dispersion processes.However,E-L simulation faces challenges in integrating droplet breakage effectively.To address this issue,our research introduces a probabilistic approach for droplet breakages.It assumes that a longer time increases the likelihood of breakup;a droplet breaks if the calculated probability exceeds a random value from 0 to 1.Consequently,the simulated breakage frequency becomes independent of the Lagrangian time step.The Sauter mean diameter and droplet size distribution can be accurately predicted by this probabilistic approach.By closely monitoring droplet motion,we reveal the complexity of droplet trajectories and the detailed patterns of circulation in stirred tanks.These insights contribute to a deeper understanding of liquidliquid dispersion dynamics.
文摘1.Introduction Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)is a powerful tool for simulating dense gas-solid reacting flows,which is essential in combustion,metallurgy,and waste management.Traditional methods face challenges in CFD-DEM modeling of dense gas-solid flows due to multi-scale characteristics,limiting resolution and creating simulation bottlenecks.By integrating fluid dynamics and particle behavior,it optimizes industrial processes.This review highlights advancements,applications,and challenges,emphasizing its role in sustainable engineering.
