This study presents a two-dimensional,transient model to simulate the flow and thermal behavior of CO_(2) within a fracturing wellbore.The model accounts for high-velocity flow within the tubing and radial heat exchan...This study presents a two-dimensional,transient model to simulate the flow and thermal behavior of CO_(2) within a fracturing wellbore.The model accounts for high-velocity flow within the tubing and radial heat exchange between the wellbore and surrounding formation.It captures the temporal evolution of temperature,pressure,flow velocity,and fluid density,enabling detailed analysis of phase transitions along different tubing sections.The influence of key operational and geological parameters,including wellhead pressure,injection velocity,inlet temperature,and formation temperature gradient,on the wellbore’s thermal and pressure fields is systematically investigated.Results indicate that due to intense convective transport by the high-speed CO_(2) flow,the temperature and velocity within the tubing are primarily governed by the inlet temperature and injection velocity,with relatively minor influence from radial heat transfer with the formation.The pressure,flow velocity,and density of CO_(2) within the tubing are strongly dependent on wellhead conditions.Frictional losses and well depth contribute to pressure variations,particularly in the horizontal section of the wellbore,where a noticeable pressurization effect increases the fluid density.During injection,liquid CO_(2) initially undergoes a rapid transition to a supercritical state,with the depth at which this phase change occurs stabilizing as injection progresses.展开更多
Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of un...Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.展开更多
The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic inter...The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.展开更多
The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized secon...The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
This paper presents a detailed investigation of unsteady supersonic flows around a typical two-body configuration, which consists of a capsule and a canopy. The cases with different trailing distances between the caps...This paper presents a detailed investigation of unsteady supersonic flows around a typical two-body configuration, which consists of a capsule and a canopy. The cases with different trailing distances between the capsule and canopy are simulated. The objective of this study is to examine the detailed effects of trailing distance on the flow fields and analyze the flow physics of the different flow modes around the parachute-like two-body model. The computational results show unsteady pulsating flow fields in the small trailing distance cases and are in reasonable agree- ment with the experimental data. As the trailing distance increases, this unsteady flow mode takes different forms along with the wake/shock and shock/shock interactions, and then gradually fades away and transits to oscillate mode, which is very different from the former. As the trailing distance keeps increasing, only the capsule wake/canopy shock interaction is present in the flow field around the two-body model, which reveals that the unsteady capsule shock/canopy shock interaction is a key mechanism for the pulsation mode.展开更多
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 new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation （SOR） iteration method and Jacobian matrix d...A new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation （SOR） iteration method and Jacobian matrix diagonally dominant splitting method of LUSGS. One advantage of this algorithm is the second-order accuracy because of no factorization error. Another advantage is the low computational cost because the Jacobian matrices and fluxes are only calculated once in each physical time step. And, the SOR algorithm has better convergence property than Gauss-Seidel. To investigate its accuracy and convergency, several unsteady flow computational tests are carded out by using the proposed SOR algorithm. Roe＇s FDS scheme is used to discritize the inviscid flux terms. Unsteady computational results of SOR are compared with the experiment results and those of Gauss-Seidel, Results reveal that the numerical results agree well with the experimental data and the second-order accuracy can be obtained as the Gauss-Seidel for unsteady flow computations. The impact of SOR factor is investigated for unsteady computations by using different SOR factors in this algorithm to simulate each computational test. Different numbers of inner iterations are needed to converge to the same criterion for different SOR factors and optimal choice of SOR factor can improve the computational efficiency greatly.展开更多
The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theor...The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 （2004））. By analyzing the distribution of the finite-time Lyapunov exponent （FTLE） along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.展开更多
Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studie...Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.展开更多
The unsteady hydrodynamics of a biomimetic fin attached to a cylindrical body has been studied numerically using a computational fluid dynamic (CFD) simulator based on an in-house solver of the Navier-Stokes equations...The unsteady hydrodynamics of a biomimetic fin attached to a cylindrical body has been studied numerically using a computational fluid dynamic (CFD) simulator based on an in-house solver of the Navier-Stokes equations, combined with a recently developed multi-block, overset grid method. The fin-body CFD model is based on a mechanical pectoral fin device, which consists of a cylindrical body and an asymmetric fin and can mimic flapping, rowing and feathering motions of the pectoral fins in fishes. First the multi-block, overset grid method incorporated into the NS solver was verified through an extensive study of unsteady flows past a single fin undergoing rowing and feathering motion. Then unsteady flows past the biomimetic fin-body model undergoing the same motions were computed and compared with the measurements of forces of the mechanical pectoral fin, which shows good agreement in both time-varying and time-averaged hydrodynamic forces. The relationship between force generation and vortex dynamics points to the importance of the match in fin kinematics between power and recovery strokes and implies that an optimal selection of parameters of phase lags between and amplitudes of rowing and feathering motions can improve the performance of labriform propulsion in terms of either maximum force generation or minimum mechanical power.展开更多
The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic diff...The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.展开更多
After the linear analytical method of unsteady flow theory is further improved,an innovative and faster algorithm is introduced.The water storage in a confined aquifer is derived from the water transmissivity coeffici...After the linear analytical method of unsteady flow theory is further improved,an innovative and faster algorithm is introduced.The water storage in a confined aquifer is derived from the water transmissivity coefficient and the water-pressure conductivity coefficient.The water transmissivity coefficient is approximated by a Taylorseries expansion of drawdown,and the water-pressure conductivity coefficient is obtained by the average drawdown.In this algorithm,the distance of the observation points from the pumping well must be short.When the distance is as short as the radius of the main pumping well,the data of the drawdown difference between the sidewall and the center of pumping well are difficult to measure,but the same results can be achieved based on the assumption that the drawdown difference approximates to the drawdown of the observation wells at a radial distance from the pumping well according to the algorithm.Without the help of charts,this algorithm is more concise and efficient,which has been verified by the test of water pumping project in Tianjin Binhai International Airport.展开更多
The study presents an experimental exploration into the mode transition of an overunder TBCC(Turbine-Based Combined Cycle)inlet,with a specific emphasis on the flow characteristics at off-design transition Mach number...The study presents an experimental exploration into the mode transition of an overunder TBCC(Turbine-Based Combined Cycle)inlet,with a specific emphasis on the flow characteristics at off-design transition Mach number.A systematic investigation was undertaken into the mode transition characteristics in both unthrottled and throttled conditions within a highspeed duct,employing high speed Schlieren and dynamic pressure acquisition systems.The results show that the high-speed duct faced flow oscillations primarily dictated by the separation bubble near the duct entrance during the downward rotation of splitter,leading to the duct’s unstart under the unthrottled condition.During the splitter’s reverse rotation,a notable hysteresis of unstart/restart of the high-speed duct was observed.Conversely,hysteresis vanishes when the initial flowfield nears the critical state owing to downstream throttling.Moreover,the oscillatory diversity,a distinctive characteristic of the high-speed duct,was firstly observed during the mode transition induced by throttling.The flow evolution was divided into four stages:an initial instability stage characterized by low-frequency oscillations below 255 Hz induced by shock train self-excitation oscillation and high-frequency oscillations around 1367 Hz caused by the movement of separation bubble.This stage is succeeded by the“big buzz”phase,comprised of pressure accumulation/release within the overflow-free duct and shock motion outside the duct to retain dynamic flow balance.The dominant frequency escalated with the increase of the internal contraction ratio in the range of 280 Hz to 400 Hz.This was followed by a high-frequency oscillation stage around 453 Hz dominated by a large internal contraction ratio with low pulsating energy,accompanied by a continuous supersonic overflow.Lastly,as the splitter gradually intersected the boundary layer of the first-stage compression surface,the capture area and the turbulence intensity of the incoming flow underwent a sudden shift,leading to a more diverse flow oscillation within the duct,manifested as various forms of mixed buzz.展开更多
To improve the performance of the positive displacement blower, it is imperative to understand the detailed internal flow characteristics or enable a visualization of flow status. However, the existing two-dimensional...To improve the performance of the positive displacement blower, it is imperative to understand the detailed internal flow characteristics or enable a visualization of flow status. However, the existing two-dimensional unsteady, three-dimensional steady or quasi-unsteady numerical simulation and theoretical analysis cannot provide the detailed flow information, which is unfavorable to improve the performance of positive displacement blower. Therefore, the unsteady flow characteristics in a three-lobe positive displacement blower are numerically investigated by solving the three-dimensional, unsteady, compressible Navier-Stokes equations coupled with RNG k-e turbulent model. In the numerical simulation, the dynamic mesh technique and overset mesh updating method are adopted. Due to the air being compressed in the process of the rotors rotating, the variation of the temperature field in the positive displacement blower is considered. By comparing the experimental measurements and the numerical results on the variation of flow rate with the outlet pressure, the maximum relative error of the flow rate is less than 2.15% even at the maximum outlet pressure condition, which means that the calculation model and numerical computational method used are effective. The numerical results show that in the intake region, the fluctuations of the inlet flow are greatly affected by the direction of the velocity vectors. In the exhaust region, the temperature changes significantly, which leads to the increase of the airflow pulsation. Through analysis on the velocity, pressure and temperature fields obtained from the numerical simulations, three-dimensional unsteady flow characteristics in the positive displacement blower are revealed. The studied results will provide useful reference for improving the performance and empirical correction in the design of the positive displacement blower.展开更多
An excited experiment system of 20# steel pipe was established with oil cylinder, 20# steel pipe, frequency converter, pump station and wave exciter generating unsteady flow artificially. The experimental results show...An excited experiment system of 20# steel pipe was established with oil cylinder, 20# steel pipe, frequency converter, pump station and wave exciter generating unsteady flow artificially. The experimental results showed that the 20# steel pipe could vibrate with the excitation of unsteady flows, and the vibration was periodic, instead of a harmonic one. Particles on the front and rear positions of pipe vibrated synchronously, and the vibration intensity of the pipe＇s two ends was greater than in the middle. System pressure and wave exciter＇s frequency had much influence upon pipe＇s amplitude. Pipe＇s vibration frequency was little affected by system pressure, and its value was close to the wave exciter＇s. Therefore, the active control of pipe＇s vibration can be realized by setting system pressure and adjusting frequency converter＇s frequency.展开更多
Left ventricular assist device( LVAD) in this study is a mechanical tool that is used to support blood flow in the patient with heart disease. It supports left ventricle by building up the pressure to the pump outlet ...Left ventricular assist device( LVAD) in this study is a mechanical tool that is used to support blood flow in the patient with heart disease. It supports left ventricle by building up the pressure to the pump outlet connected to the aorta. This pump was designed based on the magnetic driven centrifugal pump with a unique small washout hole constructed inside the impeller to generate the washout flow passage to prevent the stagnation at the region underneath and around the rotor. Computational fluid dynamics( CFD) was adopted in this study to assess the performance and optimize the design to avoid recirculation and high shear stress which is the main cause of stagnation and blood damage. Transient simulation was used for this study due to the asymmetric design of the washout hole and the complication of the bottom support of the impeller that has a risk of thrombosis,also,it was used to predict the variation of hydraulic performance caused by the rotation of the impeller and pulsed flow at the pump inlet. The simulation results show no excessive stress and no recirculation observed within the computational domain; in addition,the research result also provides information for further optimization and development to the pump.展开更多
Stall phenomena increase the complexity of the internal flow in centrifugal pump impellers.In order to tackle this problem,in the present work,a large eddy simulation(LES)approach is applied to determine the character...Stall phenomena increase the complexity of the internal flow in centrifugal pump impellers.In order to tackle this problem,in the present work,a large eddy simulation(LES)approach is applied to determine the characteristics of these unstable flows.Moreover,a vorticity identification method is used to characterize quantitatively the vortex position inside the impeller and its influencing area.By comparing the outcomes of the numerical simulations and experimental results provided by a Particle Image Velocimetry(PIV)technique,it is shown that an apparent“alternating stall”phenomenon exists inside the impeller when relatively small flow rate conditions are considered.The stall is generated near the suction side of the blade inlet,grows towards the high-pressure side of the blade in the circumferential direction,and gradually attenuates.As the flow rate decreases,the number of stalls remains unchanged,while the related influencing area and strength gradually increase and the circumferential velocity increases.展开更多
The expressions of fluid velocity under several conditions during gas atomization according to the unsteady flow mechanics are propsed, and the basic ways to maintain the stability of fluid flow is analyzed
A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusio...A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.展开更多
基金funded by National Natural Science Foundation of China(Mechanisms of proppant-carrying transport by magnetic cross-linked microparticle grids and their degradation patterns in CO_(2) fractured cracks).
文摘This study presents a two-dimensional,transient model to simulate the flow and thermal behavior of CO_(2) within a fracturing wellbore.The model accounts for high-velocity flow within the tubing and radial heat exchange between the wellbore and surrounding formation.It captures the temporal evolution of temperature,pressure,flow velocity,and fluid density,enabling detailed analysis of phase transitions along different tubing sections.The influence of key operational and geological parameters,including wellhead pressure,injection velocity,inlet temperature,and formation temperature gradient,on the wellbore’s thermal and pressure fields is systematically investigated.Results indicate that due to intense convective transport by the high-speed CO_(2) flow,the temperature and velocity within the tubing are primarily governed by the inlet temperature and injection velocity,with relatively minor influence from radial heat transfer with the formation.The pressure,flow velocity,and density of CO_(2) within the tubing are strongly dependent on wellhead conditions.Frictional losses and well depth contribute to pressure variations,particularly in the horizontal section of the wellbore,where a noticeable pressurization effect increases the fluid density.During injection,liquid CO_(2) initially undergoes a rapid transition to a supercritical state,with the depth at which this phase change occurs stabilizing as injection progresses.
基金supported by the National Natural Science Foundation of China(No.92252201)the Fundamental Research Funds for the Central Universitiesthe Academic Excellence Foundation of Beihang University(BUAA)for PhD Students。
文摘Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.
基金National Natural Science Foundation of China (506460210) Chinese Specialized Research Fund for the Doctoral Program of Higher Education (20060213007)Development Program for Outstanding Young Teachers in Harbin Institute of Technology (HITQNJS.2006.046)
文摘The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.
基金The project supported by the National Natural Science Foundation of China (10372007,10002003) and CNPC Innovation Fund
文摘The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
基金supported by the National Natural Science Foundation of China(No.11702332)
文摘This paper presents a detailed investigation of unsteady supersonic flows around a typical two-body configuration, which consists of a capsule and a canopy. The cases with different trailing distances between the capsule and canopy are simulated. The objective of this study is to examine the detailed effects of trailing distance on the flow fields and analyze the flow physics of the different flow modes around the parachute-like two-body model. The computational results show unsteady pulsating flow fields in the small trailing distance cases and are in reasonable agree- ment with the experimental data. As the trailing distance increases, this unsteady flow mode takes different forms along with the wake/shock and shock/shock interactions, and then gradually fades away and transits to oscillate mode, which is very different from the former. As the trailing distance keeps increasing, only the capsule wake/canopy shock interaction is present in the flow field around the two-body model, which reveals that the unsteady capsule shock/canopy shock interaction is a key mechanism for the pulsation mode.
文摘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.
基金National Natural Science Foundation of China (10032060)Aeronautical Basic Science Foundation of China (04A51040)
文摘A new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation （SOR） iteration method and Jacobian matrix diagonally dominant splitting method of LUSGS. One advantage of this algorithm is the second-order accuracy because of no factorization error. Another advantage is the low computational cost because the Jacobian matrices and fluxes are only calculated once in each physical time step. And, the SOR algorithm has better convergence property than Gauss-Seidel. To investigate its accuracy and convergency, several unsteady flow computational tests are carded out by using the proposed SOR algorithm. Roe＇s FDS scheme is used to discritize the inviscid flux terms. Unsteady computational results of SOR are compared with the experiment results and those of Gauss-Seidel, Results reveal that the numerical results agree well with the experimental data and the second-order accuracy can be obtained as the Gauss-Seidel for unsteady flow computations. The impact of SOR factor is investigated for unsteady computations by using different SOR factors in this algorithm to simulate each computational test. Different numbers of inner iterations are needed to converge to the same criterion for different SOR factors and optimal choice of SOR factor can improve the computational efficiency greatly.
基金supported by the National Natural Science Foundation of China(Nos.11372340 and 11732016)
文摘The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 （2004））. By analyzing the distribution of the finite-time Lyapunov exponent （FTLE） along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.
基金National Natural Science Foundation of China (51176013)Ph.D. Programs Foundation of Ministry of Education of China (20091101110014)National High-tech Research and Development Program of China (2007AA050502)
文摘Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.
文摘The unsteady hydrodynamics of a biomimetic fin attached to a cylindrical body has been studied numerically using a computational fluid dynamic (CFD) simulator based on an in-house solver of the Navier-Stokes equations, combined with a recently developed multi-block, overset grid method. The fin-body CFD model is based on a mechanical pectoral fin device, which consists of a cylindrical body and an asymmetric fin and can mimic flapping, rowing and feathering motions of the pectoral fins in fishes. First the multi-block, overset grid method incorporated into the NS solver was verified through an extensive study of unsteady flows past a single fin undergoing rowing and feathering motion. Then unsteady flows past the biomimetic fin-body model undergoing the same motions were computed and compared with the measurements of forces of the mechanical pectoral fin, which shows good agreement in both time-varying and time-averaged hydrodynamic forces. The relationship between force generation and vortex dynamics points to the importance of the match in fin kinematics between power and recovery strokes and implies that an optimal selection of parameters of phase lags between and amplitudes of rowing and feathering motions can improve the performance of labriform propulsion in terms of either maximum force generation or minimum mechanical power.
基金supported by the National Natural Science Foundation of China
文摘The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.
基金Supported by Major State Basic Research Development Program of China("973" Program,No.2010CB732106)
文摘After the linear analytical method of unsteady flow theory is further improved,an innovative and faster algorithm is introduced.The water storage in a confined aquifer is derived from the water transmissivity coefficient and the water-pressure conductivity coefficient.The water transmissivity coefficient is approximated by a Taylorseries expansion of drawdown,and the water-pressure conductivity coefficient is obtained by the average drawdown.In this algorithm,the distance of the observation points from the pumping well must be short.When the distance is as short as the radius of the main pumping well,the data of the drawdown difference between the sidewall and the center of pumping well are difficult to measure,but the same results can be achieved based on the assumption that the drawdown difference approximates to the drawdown of the observation wells at a radial distance from the pumping well according to the algorithm.Without the help of charts,this algorithm is more concise and efficient,which has been verified by the test of water pumping project in Tianjin Binhai International Airport.
基金funded by the National Natural Science Foundation of China(Nos.12025202,U20A2070 and 12172175)the National Science and Technology Major Project,China(No.J2019-Ⅱ-0014-0035)+2 种基金the Postdoctoral Fellowship Program of CPSF,China(No.GZB20230970)the Science Center for Gas Turbine Project,China(Nos.P2022-C-Ⅱ-002-001 and P2022-A-Ⅱ-002-001)the Young Scientific and Technological Talents Project of Jiangsu Association for Science and Technology,China(No.TJ-2021-052).
文摘The study presents an experimental exploration into the mode transition of an overunder TBCC(Turbine-Based Combined Cycle)inlet,with a specific emphasis on the flow characteristics at off-design transition Mach number.A systematic investigation was undertaken into the mode transition characteristics in both unthrottled and throttled conditions within a highspeed duct,employing high speed Schlieren and dynamic pressure acquisition systems.The results show that the high-speed duct faced flow oscillations primarily dictated by the separation bubble near the duct entrance during the downward rotation of splitter,leading to the duct’s unstart under the unthrottled condition.During the splitter’s reverse rotation,a notable hysteresis of unstart/restart of the high-speed duct was observed.Conversely,hysteresis vanishes when the initial flowfield nears the critical state owing to downstream throttling.Moreover,the oscillatory diversity,a distinctive characteristic of the high-speed duct,was firstly observed during the mode transition induced by throttling.The flow evolution was divided into four stages:an initial instability stage characterized by low-frequency oscillations below 255 Hz induced by shock train self-excitation oscillation and high-frequency oscillations around 1367 Hz caused by the movement of separation bubble.This stage is succeeded by the“big buzz”phase,comprised of pressure accumulation/release within the overflow-free duct and shock motion outside the duct to retain dynamic flow balance.The dominant frequency escalated with the increase of the internal contraction ratio in the range of 280 Hz to 400 Hz.This was followed by a high-frequency oscillation stage around 453 Hz dominated by a large internal contraction ratio with low pulsating energy,accompanied by a continuous supersonic overflow.Lastly,as the splitter gradually intersected the boundary layer of the first-stage compression surface,the capture area and the turbulence intensity of the incoming flow underwent a sudden shift,leading to a more diverse flow oscillation within the duct,manifested as various forms of mixed buzz.
基金Supported by Fundamental Research Funds for the Central Universities,China(Grant No.xjj20100073)Science and Technology Innovation Project of Shaanxi Province of China(Grant No.2011KTCL01-04)
文摘To improve the performance of the positive displacement blower, it is imperative to understand the detailed internal flow characteristics or enable a visualization of flow status. However, the existing two-dimensional unsteady, three-dimensional steady or quasi-unsteady numerical simulation and theoretical analysis cannot provide the detailed flow information, which is unfavorable to improve the performance of positive displacement blower. Therefore, the unsteady flow characteristics in a three-lobe positive displacement blower are numerically investigated by solving the three-dimensional, unsteady, compressible Navier-Stokes equations coupled with RNG k-e turbulent model. In the numerical simulation, the dynamic mesh technique and overset mesh updating method are adopted. Due to the air being compressed in the process of the rotors rotating, the variation of the temperature field in the positive displacement blower is considered. By comparing the experimental measurements and the numerical results on the variation of flow rate with the outlet pressure, the maximum relative error of the flow rate is less than 2.15% even at the maximum outlet pressure condition, which means that the calculation model and numerical computational method used are effective. The numerical results show that in the intake region, the fluctuations of the inlet flow are greatly affected by the direction of the velocity vectors. In the exhaust region, the temperature changes significantly, which leads to the increase of the airflow pulsation. Through analysis on the velocity, pressure and temperature fields obtained from the numerical simulations, three-dimensional unsteady flow characteristics in the positive displacement blower are revealed. The studied results will provide useful reference for improving the performance and empirical correction in the design of the positive displacement blower.
基金Funded by the National Natural Science Foundation of China(No. 50775154)the Natural Science Foundation for Young Scientists of Shanxi Province(No. 2010201025)
文摘An excited experiment system of 20# steel pipe was established with oil cylinder, 20# steel pipe, frequency converter, pump station and wave exciter generating unsteady flow artificially. The experimental results showed that the 20# steel pipe could vibrate with the excitation of unsteady flows, and the vibration was periodic, instead of a harmonic one. Particles on the front and rear positions of pipe vibrated synchronously, and the vibration intensity of the pipe＇s two ends was greater than in the middle. System pressure and wave exciter＇s frequency had much influence upon pipe＇s amplitude. Pipe＇s vibration frequency was little affected by system pressure, and its value was close to the wave exciter＇s. Therefore, the active control of pipe＇s vibration can be realized by setting system pressure and adjusting frequency converter＇s frequency.
文摘Left ventricular assist device( LVAD) in this study is a mechanical tool that is used to support blood flow in the patient with heart disease. It supports left ventricle by building up the pressure to the pump outlet connected to the aorta. This pump was designed based on the magnetic driven centrifugal pump with a unique small washout hole constructed inside the impeller to generate the washout flow passage to prevent the stagnation at the region underneath and around the rotor. Computational fluid dynamics( CFD) was adopted in this study to assess the performance and optimize the design to avoid recirculation and high shear stress which is the main cause of stagnation and blood damage. Transient simulation was used for this study due to the asymmetric design of the washout hole and the complication of the bottom support of the impeller that has a risk of thrombosis,also,it was used to predict the variation of hydraulic performance caused by the rotation of the impeller and pulsed flow at the pump inlet. The simulation results show no excessive stress and no recirculation observed within the computational domain; in addition,the research result also provides information for further optimization and development to the pump.
基金This research was funded by the Zhejiang Provincial Natural Science Foundation of China(Grant Nos.LGG21E090002,LY21E060004,LGG21E090003)National Natural Science Foundation of China(Grant No.51779226)the China Postdoctoral Science Foundation(Grant No.2021M691383).
文摘Stall phenomena increase the complexity of the internal flow in centrifugal pump impellers.In order to tackle this problem,in the present work,a large eddy simulation(LES)approach is applied to determine the characteristics of these unstable flows.Moreover,a vorticity identification method is used to characterize quantitatively the vortex position inside the impeller and its influencing area.By comparing the outcomes of the numerical simulations and experimental results provided by a Particle Image Velocimetry(PIV)technique,it is shown that an apparent“alternating stall”phenomenon exists inside the impeller when relatively small flow rate conditions are considered.The stall is generated near the suction side of the blade inlet,grows towards the high-pressure side of the blade in the circumferential direction,and gradually attenuates.As the flow rate decreases,the number of stalls remains unchanged,while the related influencing area and strength gradually increase and the circumferential velocity increases.
文摘The expressions of fluid velocity under several conditions during gas atomization according to the unsteady flow mechanics are propsed, and the basic ways to maintain the stability of fluid flow is analyzed
文摘A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.
