Modeling of unsteady aerodynamic loads at high angles of attack using a small amount of experimental or simulation data to construct predictive models for unknown states can greatly improve the efficiency of aircraft ...Modeling of unsteady aerodynamic loads at high angles of attack using a small amount of experimental or simulation data to construct predictive models for unknown states can greatly improve the efficiency of aircraft unsteady aerodynamic design and flight dynamics analysis.In this paper,aiming at the problems of poor generalization of traditional aerodynamic models and intelligent models,an intelligent aerodynamic modeling method based on gated neural units is proposed.The time memory characteristics of the gated neural unit is fully utilized,thus the nonlinear flow field characterization ability of the learning and training process is enhanced,and the generalization ability of the whole prediction model is improved.The prediction and verification of the model are carried out under the maneuvering flight condition of NACA0015 airfoil.The results show that the model has good adaptability.In the interpolation prediction,the maximum prediction error of the lift and drag coefficients and the moment coefficient does not exceed 10%,which can basically represent the variation characteristics of the entire flow field.In the construction of extrapolation models,the training model based on the strong nonlinear data has good accuracy for weak nonlinear prediction.Furthermore,the error is larger,even exceeding 20%,which indicates that the extrapolation and generalization capabilities need to be further optimized by integrating physical models.Compared with the conventional state space equation model,the proposed method can improve the extrapolation accuracy and efficiency by 78%and 60%,respectively,which demonstrates the applied potential of this method in aerodynamic modeling.展开更多
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 characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is p...The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.展开更多
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 pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and th...The pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and the three dimensional unsteady turbulent simulation concerning the rotor-stator interaction is executed by a Navier-Stoke solver embedded with k -ε turbulence model and with appropriate moving interface boundary conditions. Detecting points are distributed in the flow passage in different radial and circumferential positions to capture the static pressure fluctuation character for one cycle of the impeller. The time-domain spectrums show that the static pressure curves are periodic and have five peaks and five valleys. With the radius increasing, the pressure fluctuation peak-to-peak values in the impeller are increasing, and reach the maximum value on the interface. In the casing flow passage, those values are about 7% of local static pressure except some ones near the tongue. The values become decreasingly in the diffuser pipe. The frequency spectrums transformed by fast Fourier transform (FFT) show that the dominant frequency is approximate with the blade passing frequency, and the pressure fluctuations in impeller passage have high frequency content while those in casing ones have no such information.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54...In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54×106. The calculations were performed by means of solving the 2-D Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations with a k-ε turbulence model. The calculated results, produced flow structure drag and lift coefficients, as well as Strouhal numbers. The findings were in good agreement with previous published data, which also supplied us with a good understanding of the flow across cylinders of different high Reynolds numbers. Meanwhile, an effective measure was presented to control the lift force on a cylinder, which points the way to decrease the vortex induced vibration of marine structure in future.展开更多
An analysis of unsteady boundary layer flow and heat transfer over an exponentially shrinking porous sheet filled with a copper-water nanofluid is presented.Water is treated as a base fluid.In the investigation,non-un...An analysis of unsteady boundary layer flow and heat transfer over an exponentially shrinking porous sheet filled with a copper-water nanofluid is presented.Water is treated as a base fluid.In the investigation,non-uniform mass suction through the porous sheet is considered.Using Keller-box method the transformed equations are solved numerically.The results of skin friction coefficient,the local Nusselt number as well as the velocity and temperature profiles are presented for different flow parameters.The results showed that the dual non-similar solutions exist only when certain amount of mass suction is applied through the porous sheet for various unsteady parameters and nanoparticle volume fractions.The ranges of suction where dual non-similar solution exists,become larger when values of unsteady parameter as well as nanoparticle volume fraction increase.So,due to unsteadiness of flow dynamics and the presence of nanoparticles in flow field,the requirement of mass suction for existence of solution of boundary layer flow past an exponentially shrinking sheet is less.Furthermore,the velocity boundary layer thickness decreases and thermal boundary layer thickness increases with increasing of nanoparticle volume fraction in both non-similar solutions.Whereas,for stronger mass suction,the velocity boundary layer thickness becomes thinner for the first solution and the effect is opposite in the case of second solution.The temperature inside the boundary layer increases with nanoparticle volume fraction and decreases with mass suction.So,for the unsteadiness and for the presence of nanoparticles,the flow separation is delayed to some extent.展开更多
In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. I...In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.展开更多
Numerical simulations of unsteady cavitating flow around a NACA66-mod hydrofoil were performed using the partially-averaged Navier-Stokes method with different values of the resolution control parameters (f k =1.0-0....Numerical simulations of unsteady cavitating flow around a NACA66-mod hydrofoil were performed using the partially-averaged Navier-Stokes method with different values of the resolution control parameters (f k =1.0-0.2,fε =1).With decreasing fk,the predicted cavitating flow becomes unsteady as the time-averaged turbulent viscosity at the rear part of the attached cavity is gradually reduced.For fk =0.9 and 0.8,the cavity becomes unstable and its length dramatically expands and shrinks,but the calculation fails to predict the vapor cloud shedding behavior observed experimentally.With smaller fk less than 0.7,the cloud shedding behavior is simulated numerically and the predicted cavity shedding frequency increases.With fk =0.2,the whole cavitating flow evolution can be reasonably reproduced including the cavity growth/destabilization observed previously.The reentrant flow along the suction surface of the hydrofoil is the main trigger to cause the vapor cloud shedding.The wall pressure along the hydrofoil surface oscillates greatly due to the dynamic cavity shedding.Comparing the simulations and experiments,it is confirmed that for the PANS method,resolution control parameters of fk =0.2and fε =1 are recommended for numerical simulations of unsteady cavitating flows.Thus,the present study shows that the PANS method is an effective approach for predicting unsteady cavitating flow over hydrofoils.展开更多
Aerodynamic forces and power requirements in forward flight in a bumblebee (Bombus terrestris) were studied using the method of computational fluid dynamics. Actual wing kinematic data of free flight were used in th...Aerodynamic forces and power requirements in forward flight in a bumblebee (Bombus terrestris) were studied using the method of computational fluid dynamics. Actual wing kinematic data of free flight were used in the study (the speed ranges from 0 m/s to 4.5 m/s; advance ratio ranges from 0-0.66). The bumblebee employs the delayed stall mechanism and the fast pitching-up rotation mechanism to produce vertical force and thrust. The leading-edge vortex does not shed in the translatory phase of the half-strokes and is much more concentrated than that of the fruit fly in a previous study. At hovering and low-speed flight, the vertical force is produced by both the half-strokes and is contributed by wing lift; at medium and high speeds, the vertical force is mainly produced during the downstroke and is contributed by both wing lift and wing drag. At all speeds the thrust is mainly produced in the upstroke and is contributed by wing drag. The power requirement at low to medium speeds is not very different from that of hovering and is relatively large at the highest speed (advance ratio 0.66), i.e. the power curve is Jshaped. Except at the highest flight speed, storing energy elastically can save power up to 20%-30%. At the highest speed, because of the large increase of aerodynamic torque and the slight decrease of inertial torque (due to the smaller stroke amplitude and stroke frequency used), the power requirement is dominated by aerodynamic power and the effect of elastic storage of energy on power requirement is limited.展开更多
The heat transfer between an airway and the air flowing though it is an unsteady problem. The governing equation of unsteady heat transfer was solved using the method of separation of variables. The solution is an inf...The heat transfer between an airway and the air flowing though it is an unsteady problem. The governing equation of unsteady heat transfer was solved using the method of separation of variables. The solution is an infinite series including Bessel functions. The theoretical solution was analyzed by solving for the positive roots of the transcendental equation by iteration. The dimensionless surface temperature of the sur- rounding rock is only affected by the Bi number but not by the thermo-physical coefficients of the rock. The dimensionless coefficient of heat transfer, k, decreases with the Fo number similarly to the influence of the Bi number on k. A formula for determining the fully developed stage (FDS) suitable for unsteady heat transfer in the airway is proposed. The FDS from theoretical analysis occurs with Fo from 1.6 to 2. The ratio of excess temperature in the surrounding rock is independent of the initial conditions and only dependent on the Bi number and the relative position in the airway, at the FDS. The calculation error is large when using just the first term from the complete series when Fo is from 2 to 12. Five terms give a solution approximately equal to that found using the complete series. The first term could replace the complete series only when Fo is greater than 12. The FDS plays an important role in predication of the temperature field of the surrounding rock and in simplified calculations.展开更多
An analysis is presented for an unsteady boundary layer stagnation-point flow of a Newtonian fluid and the heat transfer towards a stretching sheet taking non-conventional partial slip conditions at the sheet.The self...An analysis is presented for an unsteady boundary layer stagnation-point flow of a Newtonian fluid and the heat transfer towards a stretching sheet taking non-conventional partial slip conditions at the sheet.The self-similar equations are obtained using similarity transformations and solved numerically by the shooting method.Effects of the parameters involved in the equations,especially velocity slip and thermal slip parameters on the velocity and temperature profiles,are analyzed extensively.It is revealed that due to the velocity and thermal slip parameters,the rate of heat transfer from the sheet and the wall skin friction change significantly.展开更多
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.展开更多
基金supported in part by the National Natural Science Foundation of China (No. 12202363)。
文摘Modeling of unsteady aerodynamic loads at high angles of attack using a small amount of experimental or simulation data to construct predictive models for unknown states can greatly improve the efficiency of aircraft unsteady aerodynamic design and flight dynamics analysis.In this paper,aiming at the problems of poor generalization of traditional aerodynamic models and intelligent models,an intelligent aerodynamic modeling method based on gated neural units is proposed.The time memory characteristics of the gated neural unit is fully utilized,thus the nonlinear flow field characterization ability of the learning and training process is enhanced,and the generalization ability of the whole prediction model is improved.The prediction and verification of the model are carried out under the maneuvering flight condition of NACA0015 airfoil.The results show that the model has good adaptability.In the interpolation prediction,the maximum prediction error of the lift and drag coefficients and the moment coefficient does not exceed 10%,which can basically represent the variation characteristics of the entire flow field.In the construction of extrapolation models,the training model based on the strong nonlinear data has good accuracy for weak nonlinear prediction.Furthermore,the error is larger,even exceeding 20%,which indicates that the extrapolation and generalization capabilities need to be further optimized by integrating physical models.Compared with the conventional state space equation model,the proposed method can improve the extrapolation accuracy and efficiency by 78%and 60%,respectively,which demonstrates the applied potential of this method in aerodynamic modeling.
基金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 (No. 90716011)
文摘The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.
基金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.
基金supported by National Outstanding Young Scientists Funds of China (Grand No.50825902)
文摘The pressure fluctuation in the flow passage of both impeller and casing is addressed on design condition. The initial conditions for the unsteady turbulent simulation are resulted from the steady calculations, and the three dimensional unsteady turbulent simulation concerning the rotor-stator interaction is executed by a Navier-Stoke solver embedded with k -ε turbulence model and with appropriate moving interface boundary conditions. Detecting points are distributed in the flow passage in different radial and circumferential positions to capture the static pressure fluctuation character for one cycle of the impeller. The time-domain spectrums show that the static pressure curves are periodic and have five peaks and five valleys. With the radius increasing, the pressure fluctuation peak-to-peak values in the impeller are increasing, and reach the maximum value on the interface. In the casing flow passage, those values are about 7% of local static pressure except some ones near the tongue. The values become decreasingly in the diffuser pipe. The frequency spectrums transformed by fast Fourier transform (FFT) show that the dominant frequency is approximate with the blade passing frequency, and the pressure fluctuations in impeller passage have high frequency content while those in casing ones have no such information.
基金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 (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.
基金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.
基金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 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.
基金Foundation item: Supported by Supported by the National Natural Science Foundation of China (Grant No. 51009070).
文摘In this paper, 2-D computational analyses were conducted for unsteady high Reynolds number flows around a smooth circular cylinder in the supercritical and upper-transition flow regimes, i.e. 8.21×104〈Re〈1.54×106. The calculations were performed by means of solving the 2-D Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations with a k-ε turbulence model. The calculated results, produced flow structure drag and lift coefficients, as well as Strouhal numbers. The findings were in good agreement with previous published data, which also supplied us with a good understanding of the flow across cylinders of different high Reynolds numbers. Meanwhile, an effective measure was presented to control the lift force on a cylinder, which points the way to decrease the vortex induced vibration of marine structure in future.
基金the National Board for Higher Mathematics (NBHM),Department of Atomic Energy,Government of India for the financial support in pursuing this workthe financial support from MOHE and the Research Management Center-UTM through FRGS and RUG vote number 4F109 and 02H80 for this research
文摘An analysis of unsteady boundary layer flow and heat transfer over an exponentially shrinking porous sheet filled with a copper-water nanofluid is presented.Water is treated as a base fluid.In the investigation,non-uniform mass suction through the porous sheet is considered.Using Keller-box method the transformed equations are solved numerically.The results of skin friction coefficient,the local Nusselt number as well as the velocity and temperature profiles are presented for different flow parameters.The results showed that the dual non-similar solutions exist only when certain amount of mass suction is applied through the porous sheet for various unsteady parameters and nanoparticle volume fractions.The ranges of suction where dual non-similar solution exists,become larger when values of unsteady parameter as well as nanoparticle volume fraction increase.So,due to unsteadiness of flow dynamics and the presence of nanoparticles in flow field,the requirement of mass suction for existence of solution of boundary layer flow past an exponentially shrinking sheet is less.Furthermore,the velocity boundary layer thickness decreases and thermal boundary layer thickness increases with increasing of nanoparticle volume fraction in both non-similar solutions.Whereas,for stronger mass suction,the velocity boundary layer thickness becomes thinner for the first solution and the effect is opposite in the case of second solution.The temperature inside the boundary layer increases with nanoparticle volume fraction and decreases with mass suction.So,for the unsteadiness and for the presence of nanoparticles,the flow separation is delayed to some extent.
文摘In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.
基金Supported by the National Natural Science Fund of China under Grant Nos 51179091 and 50976061, the Major NationalScientific Instrument and Equipment Development Project (No 2011YQ07004901), and the China Postdoctoral Science Foundation(No 2011M500314).
文摘Numerical simulations of unsteady cavitating flow around a NACA66-mod hydrofoil were performed using the partially-averaged Navier-Stokes method with different values of the resolution control parameters (f k =1.0-0.2,fε =1).With decreasing fk,the predicted cavitating flow becomes unsteady as the time-averaged turbulent viscosity at the rear part of the attached cavity is gradually reduced.For fk =0.9 and 0.8,the cavity becomes unstable and its length dramatically expands and shrinks,but the calculation fails to predict the vapor cloud shedding behavior observed experimentally.With smaller fk less than 0.7,the cloud shedding behavior is simulated numerically and the predicted cavity shedding frequency increases.With fk =0.2,the whole cavitating flow evolution can be reasonably reproduced including the cavity growth/destabilization observed previously.The reentrant flow along the suction surface of the hydrofoil is the main trigger to cause the vapor cloud shedding.The wall pressure along the hydrofoil surface oscillates greatly due to the dynamic cavity shedding.Comparing the simulations and experiments,it is confirmed that for the PANS method,resolution control parameters of fk =0.2and fε =1 are recommended for numerical simulations of unsteady cavitating flows.Thus,the present study shows that the PANS method is an effective approach for predicting unsteady cavitating flow over hydrofoils.
基金The project supported by the National Natural Science Foundation of China(10232010)the National Aeronautic Science fund of China(03A51049)
文摘Aerodynamic forces and power requirements in forward flight in a bumblebee (Bombus terrestris) were studied using the method of computational fluid dynamics. Actual wing kinematic data of free flight were used in the study (the speed ranges from 0 m/s to 4.5 m/s; advance ratio ranges from 0-0.66). The bumblebee employs the delayed stall mechanism and the fast pitching-up rotation mechanism to produce vertical force and thrust. The leading-edge vortex does not shed in the translatory phase of the half-strokes and is much more concentrated than that of the fruit fly in a previous study. At hovering and low-speed flight, the vertical force is produced by both the half-strokes and is contributed by wing lift; at medium and high speeds, the vertical force is mainly produced during the downstroke and is contributed by both wing lift and wing drag. At all speeds the thrust is mainly produced in the upstroke and is contributed by wing drag. The power requirement at low to medium speeds is not very different from that of hovering and is relatively large at the highest speed (advance ratio 0.66), i.e. the power curve is Jshaped. Except at the highest flight speed, storing energy elastically can save power up to 20%-30%. At the highest speed, because of the large increase of aerodynamic torque and the slight decrease of inertial torque (due to the smaller stroke amplitude and stroke frequency used), the power requirement is dominated by aerodynamic power and the effect of elastic storage of energy on power requirement is limited.
基金provided by the Postdoctoral Science Foundation of China (No. 2011M500974)Fundamental Research Funds for Central Universities (No. 2011QNA16)+1 种基金Postdoctoral Foundation of state key laboratory for GDUE (No. PD1101)Natural Science Foundation for Young Scholars of SMCE, CUMT(No. lj2010qnjj004)
文摘The heat transfer between an airway and the air flowing though it is an unsteady problem. The governing equation of unsteady heat transfer was solved using the method of separation of variables. The solution is an infinite series including Bessel functions. The theoretical solution was analyzed by solving for the positive roots of the transcendental equation by iteration. The dimensionless surface temperature of the sur- rounding rock is only affected by the Bi number but not by the thermo-physical coefficients of the rock. The dimensionless coefficient of heat transfer, k, decreases with the Fo number similarly to the influence of the Bi number on k. A formula for determining the fully developed stage (FDS) suitable for unsteady heat transfer in the airway is proposed. The FDS from theoretical analysis occurs with Fo from 1.6 to 2. The ratio of excess temperature in the surrounding rock is independent of the initial conditions and only dependent on the Bi number and the relative position in the airway, at the FDS. The calculation error is large when using just the first term from the complete series when Fo is from 2 to 12. Five terms give a solution approximately equal to that found using the complete series. The first term could replace the complete series only when Fo is greater than 12. The FDS plays an important role in predication of the temperature field of the surrounding rock and in simplified calculations.
文摘An analysis is presented for an unsteady boundary layer stagnation-point flow of a Newtonian fluid and the heat transfer towards a stretching sheet taking non-conventional partial slip conditions at the sheet.The self-similar equations are obtained using similarity transformations and solved numerically by the shooting method.Effects of the parameters involved in the equations,especially velocity slip and thermal slip parameters on the velocity and temperature profiles,are analyzed extensively.It is revealed that due to the velocity and thermal slip parameters,the rate of heat transfer from the sheet and the wall skin friction change significantly.
文摘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.