To assess the aerodynamic performance and vibration characteristics of rotor blades during rotation,a study of unsteady blade surface forces is conducted in a low-speed axial flow compressor under a rotating coordinat...To assess the aerodynamic performance and vibration characteristics of rotor blades during rotation,a study of unsteady blade surface forces is conducted in a low-speed axial flow compressor under a rotating coordinate system.The capture,modulation,and acquisition of unsteady blade surface forces are achieved by using pressure sensors and strain gauges attached to the rotor blades,in conjunction with a wireless telemetry system.Based on the measurement reliability verification,this approach allows for the determination of the static pressure distribution on rotor blade surfaces,enabling the quantitative description of loadability at different spanwise positions along the blade chord.Effects caused by the factors such as Tip Leakage Flow(TLF)and flow separation can be perceived and reflected in the trends of static pressure on the blade surfaces.Simultaneously,the dynamic characteristics of unsteady pressure and stress on the blade surfaces are analyzed.The results indicate that only the pressure signals measured at the mid-chord of the blade tip can distinctly detect the unsteady frequency of TLF due to the oscillation of the low-pressure spot on the pressure surface.Subsequently,with the help of one-dimensional continuous wavelet analysis method,it can be inferred that as the compressor enters stall,the sensors are capable of capturing stall cell frequency under a rotating coordinate system.Furthermore,the stress at the blade root is higher than that at the blade tip,and the frequency band of the vibration can also be measured by the pressure sensors fixed on the casing wall in a stationary frame.While the compressor stalls,the stress at the blade root can be higher,which can provide valuable guidance for monitoring the lifecycle of compressor blades.展开更多
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.展开更多
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.展开更多
Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,bound...Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,boundary layer thickness and Stanton number demonstrate that the flow changes dramatically in the shock wave/turbulent boundary layer interaction area.It is found that the steady jet has no effect on suppressing flow separation unexpectedly,but increases its spatial scale instead.Instantaneous flow structures show that the turbulence amplification can be observed after the application of flow control,and abundant Görtler-like vorticities appear,but the strength of the main shock decreases.Analyzing the wall fluctuating pressure signals using weighted power spectral density,we found an interesting thing.That is,although the low-frequency oscillation phenomenon induced by separation shock is suppressed by the steady jet,wall fluctuating pressure beneath the jet shock is oscillating at a frequency lower than 0.1u∞/δref.Results of coherent and intermittency factor reveal that it is related to the backand-forth movement of the jet shock itself.展开更多
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.展开更多
Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous f...Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.展开更多
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.展开更多
Unsteady supersonic base flows around three afterbodies, cylindrical （Cy）, boattailed （BT） and three-step （MS）, are investigated in this paper. Reynolds-averaged Navier-Stokes （RANS） and two RANS/LES （large-e...Unsteady supersonic base flows around three afterbodies, cylindrical （Cy）, boattailed （BT） and three-step （MS）, are investigated in this paper. Reynolds-averaged Navier-Stokes （RANS） and two RANS/LES （large-eddy simulation） hybrid methods, detached eddy simulation （DES） and delayed-DES （DDES）, are used to predict the base flow characteristics around the baseline Cy afterbody. All the RANS and hybrid methods are based on the two-equation SST （shear-stress transport） model with compressible corrections （CC）. According to the comparison of measurements, both DES and DDES can produce more satisfactory results than RANS. RANS can only present the ＂stable＂ flow pat- terns, while the hybrid methods can demonstrate unsteady flow structures. DDES and DES results are little different from one another although the latter exhibits better agreement with the experiment. DES is taken to investigate the 5° BT and three-step afterbodies. The mean flow data and the instantaneous turbulent coherent structures are compared against available measurements.展开更多
Cavitation has a significant e ect on the flow fields and structural behaviors of a centrifugal pump. In this study, the unsteady flow and structural behaviors of a centrifugal pump are investigated numerically under ...Cavitation has a significant e ect on the flow fields and structural behaviors of a centrifugal pump. In this study, the unsteady flow and structural behaviors of a centrifugal pump are investigated numerically under di erent cavitation conditions. A strong two-way coupling fluid-structure interaction simulation is applied to obtain interior views of the e ects of cavitating bubbles on the flow and structural dynamics of a pump. The renormalization-group k-ε turbulence model and the Zwart–Gerbe–Belamri cavitation model are solved for the fluid side, while a transient structural dynamic analysis is employed for the structure side. The di erent cavitation states are mapped in the head-net positive suction head(H-NPSH) curves and flow field features inside the impeller are fully revealed. Results indicate that cavitating bubbles grow and expand rapidly with decreasing NPSH. In addition, the pressure fluctuations, both in the impeller and volute, are quantitatively analyzed and associated with the cavitation states. It is shown that influence of the cavitation on the flow field is critical, specifically in the super-cavitation state. The e ect of cavitation on the unsteady radial force and blade loads is also discussed. The results indicate that the averaged radial force increased from 8.5 N to 54.4 N in the transition progress from an onset cavitation state to a super-cavitation state. Furthermore, the structural behaviors, including blade deformation, stress, and natural frequencies, corresponding to the cavitation states are discussed. A large volume of cavitating bubbles weakens the fluid forces on the blade and decreases the natural frequencies of the rotor system. This study could enhance the understanding of the e ects of cavitation on pump flow and structural behaviors.展开更多
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.展开更多
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 fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases ...The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases are solved and the exact solutions are obtained by using the Weber transform and the Laplace transform for fractional calculus.展开更多
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.展开更多
An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS...An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS equations and physical analysis and relevant discussions are then presented.展开更多
Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panM...Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.展开更多
In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Backgroun...In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Background Oriented Schlieren technique (BOS) allows the measurement of the light deflection caused by density gradients in a compressible flow. For this purpose the distortion of the image of a background pattern observed through the flow is used. In order to increase the performance of the conventional Background Oriented Schlieren technique, the monochromatic background is replaced by a colored dot pattern. The different colors are treated separately using suitable correlation algorithms. Therefore, the precision and the spatial resolution can be highly increased. Furthermore a special arrangement of the different colored dot patterns in the background allows astigmatism in the region with high density gradients to be overcome. For the first time an algebraic reconstruction technique (ART) is then used to reconstruct the density field of unsteady flows around a spike-tipped model from CBOS measurements. The obtained images reveal the interaction between the free-stream flow and the high-pressure region in front of the model, which leads to large-scale instabilities in the flow.展开更多
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.展开更多
In previous studies,a theoretical model was developed after Acoustic Resonance(AR)was experimentally detected in a four-stage compressor,and AR inception was proposed to be triggered by an unknown sound source,which i...In previous studies,a theoretical model was developed after Acoustic Resonance(AR)was experimentally detected in a four-stage compressor,and AR inception was proposed to be triggered by an unknown sound source,which is a pressure perturbation of a specific frequency with a suitable circumferential propagation speed.The present paper,which is not dedicated to the simulation of acoustic field,aims to identify the specific sound source generated by the unsteady tip leakage flow using the unsteady Computational Fluid Mechanics(CFD)approach.After a comprehensive analysis of an Unsteady Reynolds Averaged Navier-Stokes(URANS)simulation,a pressure perturbation of non-integer multiple of rotor frequency is found at the blade tip.Since the essence of the tip leakage flow is a jet flow driven by the pressure difference between two sides of blade,a simplified tip leakage flow model is adopted using Large Eddy Simulation(LES)in order to simulate the jet flow through a tip clearance.It is found that the convection velocity of shedding vortices fits the expected propagation speed of the sound source,the frequency is also close to one of the dominating frequencies in the URANS simulation,and the resultant combination frequency coincides with the experimentally measured AR frequency.Since such a simplified model successfully captures the key physical mechanisms,it is concluded that this paper provides a piece of unambiguous evidence on the role of unsteady tip leakage vortex in triggering the AR inception of the multistage compressor.展开更多
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.展开更多
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.展开更多
基金funded by the National Natural Science Foundation of China(Nos.U24A20138 and No.52376039)the Beijing Natural Science Foundation,China(JQ24017)+1 种基金the National Science and Technology Major Project of China(Nos.J2019-II-0005-0025 and Y2022-II-0002-0005)the Special Fund for the Member of Youth Innovation Promotion Association of Chinese Academy of Sciences,China(No.2018173).
文摘To assess the aerodynamic performance and vibration characteristics of rotor blades during rotation,a study of unsteady blade surface forces is conducted in a low-speed axial flow compressor under a rotating coordinate system.The capture,modulation,and acquisition of unsteady blade surface forces are achieved by using pressure sensors and strain gauges attached to the rotor blades,in conjunction with a wireless telemetry system.Based on the measurement reliability verification,this approach allows for the determination of the static pressure distribution on rotor blade surfaces,enabling the quantitative description of loadability at different spanwise positions along the blade chord.Effects caused by the factors such as Tip Leakage Flow(TLF)and flow separation can be perceived and reflected in the trends of static pressure on the blade surfaces.Simultaneously,the dynamic characteristics of unsteady pressure and stress on the blade surfaces are analyzed.The results indicate that only the pressure signals measured at the mid-chord of the blade tip can distinctly detect the unsteady frequency of TLF due to the oscillation of the low-pressure spot on the pressure surface.Subsequently,with the help of one-dimensional continuous wavelet analysis method,it can be inferred that as the compressor enters stall,the sensors are capable of capturing stall cell frequency under a rotating coordinate system.Furthermore,the stress at the blade root is higher than that at the blade tip,and the frequency band of the vibration can also be measured by the pressure sensors fixed on the casing wall in a stationary frame.While the compressor stalls,the stress at the blade root can be higher,which can provide valuable guidance for monitoring the lifecycle of compressor blades.
基金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.
基金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(Grant Nos.12202488 and 12002377)the Natural Science Program of National University of Defense Technology(Grant No.ZK22-30)Independent Cultivation Project for Young Talents in College of Aerospace Science and Engineering.
文摘Direct numerical simulations of Mach 6 hypersonic flow over a 34°compression corner subject to steady jet are conducted.Distributions of skin friction coefficient,wall pressure,mean velocity and temperature,boundary layer thickness and Stanton number demonstrate that the flow changes dramatically in the shock wave/turbulent boundary layer interaction area.It is found that the steady jet has no effect on suppressing flow separation unexpectedly,but increases its spatial scale instead.Instantaneous flow structures show that the turbulence amplification can be observed after the application of flow control,and abundant Görtler-like vorticities appear,but the strength of the main shock decreases.Analyzing the wall fluctuating pressure signals using weighted power spectral density,we found an interesting thing.That is,although the low-frequency oscillation phenomenon induced by separation shock is suppressed by the steady jet,wall fluctuating pressure beneath the jet shock is oscillating at a frequency lower than 0.1u∞/δref.Results of coherent and intermittency factor reveal that it is related to the backand-forth movement of the jet shock itself.
基金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.
基金Supported by the National Natural Science Foundation of China(10372043,11172134)the Fundingof Jiangsu Innovation Program for Graduate Education(CXZZ11-0192)~~
文摘Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.
基金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.
基金supported by the National Natural Science Foundation of China (10502030 and 90505005)Innovation and Support Foundation of Chinese Astronautics
文摘Unsteady supersonic base flows around three afterbodies, cylindrical （Cy）, boattailed （BT） and three-step （MS）, are investigated in this paper. Reynolds-averaged Navier-Stokes （RANS） and two RANS/LES （large-eddy simulation） hybrid methods, detached eddy simulation （DES） and delayed-DES （DDES）, are used to predict the base flow characteristics around the baseline Cy afterbody. All the RANS and hybrid methods are based on the two-equation SST （shear-stress transport） model with compressible corrections （CC）. According to the comparison of measurements, both DES and DDES can produce more satisfactory results than RANS. RANS can only present the ＂stable＂ flow pat- terns, while the hybrid methods can demonstrate unsteady flow structures. DDES and DES results are little different from one another although the latter exhibits better agreement with the experiment. DES is taken to investigate the 5° BT and three-step afterbodies. The mean flow data and the instantaneous turbulent coherent structures are compared against available measurements.
基金Supported by National Natural Science Foundation of China(Grant Nos.51609212,51606167)China Postdoctoral Science Foundation(Grant No.2016M590546)Zhejiang Provincial Natural Science Foundation(Grant No.2016C31043)
文摘Cavitation has a significant e ect on the flow fields and structural behaviors of a centrifugal pump. In this study, the unsteady flow and structural behaviors of a centrifugal pump are investigated numerically under di erent cavitation conditions. A strong two-way coupling fluid-structure interaction simulation is applied to obtain interior views of the e ects of cavitating bubbles on the flow and structural dynamics of a pump. The renormalization-group k-ε turbulence model and the Zwart–Gerbe–Belamri cavitation model are solved for the fluid side, while a transient structural dynamic analysis is employed for the structure side. The di erent cavitation states are mapped in the head-net positive suction head(H-NPSH) curves and flow field features inside the impeller are fully revealed. Results indicate that cavitating bubbles grow and expand rapidly with decreasing NPSH. In addition, the pressure fluctuations, both in the impeller and volute, are quantitatively analyzed and associated with the cavitation states. It is shown that influence of the cavitation on the flow field is critical, specifically in the super-cavitation state. The e ect of cavitation on the unsteady radial force and blade loads is also discussed. The results indicate that the averaged radial force increased from 8.5 N to 54.4 N in the transition progress from an onset cavitation state to a super-cavitation state. Furthermore, the structural behaviors, including blade deformation, stress, and natural frequencies, corresponding to the cavitation states are discussed. A large volume of cavitating bubbles weakens the fluid forces on the blade and decreases the natural frequencies of the rotor system. This study could enhance the understanding of the e ects of cavitation on pump flow and structural behaviors.
基金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 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 project supported by the National Natural Science Foundation of China (10272067, 10426024)the Doctoral Program Foundation of the Education Ministry of China (20030422046)the Natural Science Foundation of Shandong University at Weihai.
文摘The fractional calculus is used in the constitutive relationship model of viscoelastic fluid. A generalized Maxwell model with fractional calculus is considered. Based on the flow conditions described, two flow cases are solved and the exact solutions are obtained by using the Weber transform and the Laplace transform for fractional calculus.
基金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.
文摘An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS equations and physical analysis and relevant discussions are then presented.
文摘Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.
文摘In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Background Oriented Schlieren technique (BOS) allows the measurement of the light deflection caused by density gradients in a compressible flow. For this purpose the distortion of the image of a background pattern observed through the flow is used. In order to increase the performance of the conventional Background Oriented Schlieren technique, the monochromatic background is replaced by a colored dot pattern. The different colors are treated separately using suitable correlation algorithms. Therefore, the precision and the spatial resolution can be highly increased. Furthermore a special arrangement of the different colored dot patterns in the background allows astigmatism in the region with high density gradients to be overcome. For the first time an algebraic reconstruction technique (ART) is then used to reconstruct the density field of unsteady flows around a spike-tipped model from CBOS measurements. The obtained images reveal the interaction between the free-stream flow and the high-pressure region in front of the model, which leads to large-scale instabilities in the flow.
基金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.
基金the National Natural Science Foundation of China(Nos.51976116 and 51976125)the Natural Science Fund of Shanghai,China(No.19ZR1425900)+1 种基金the Open Research Subject of Key Laboratory(Fluid Machinery and Engineering Research Base)of Sichuan Province,China(No.szjj2019-022)Alexander von Humboldt Foundation.
文摘In previous studies,a theoretical model was developed after Acoustic Resonance(AR)was experimentally detected in a four-stage compressor,and AR inception was proposed to be triggered by an unknown sound source,which is a pressure perturbation of a specific frequency with a suitable circumferential propagation speed.The present paper,which is not dedicated to the simulation of acoustic field,aims to identify the specific sound source generated by the unsteady tip leakage flow using the unsteady Computational Fluid Mechanics(CFD)approach.After a comprehensive analysis of an Unsteady Reynolds Averaged Navier-Stokes(URANS)simulation,a pressure perturbation of non-integer multiple of rotor frequency is found at the blade tip.Since the essence of the tip leakage flow is a jet flow driven by the pressure difference between two sides of blade,a simplified tip leakage flow model is adopted using Large Eddy Simulation(LES)in order to simulate the jet flow through a tip clearance.It is found that the convection velocity of shedding vortices fits the expected propagation speed of the sound source,the frequency is also close to one of the dominating frequencies in the URANS simulation,and the resultant combination frequency coincides with the experimentally measured AR frequency.Since such a simplified model successfully captures the key physical mechanisms,it is concluded that this paper provides a piece of unambiguous evidence on the role of unsteady tip leakage vortex in triggering the AR inception of the multistage compressor.
文摘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 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.
