The aim of this study is to investigate CO2 two-phase nozzle flow in terms of both experimental and analytical aspects for the optimum design of two-phase flow nozzle of CO2 two-phase flow ejector. In the experiment, ...The aim of this study is to investigate CO2 two-phase nozzle flow in terms of both experimental and analytical aspects for the optimum design of two-phase flow nozzle of CO2 two-phase flow ejector. In the experiment, it is measured that the temperature profile in the stream-wise direction of a divergent-convergent nozzle through which CO2 in the supercritical pressure condition is blown down into the atmosphere. In the analysis, a one-dimensional model which assumes steady, adiabatic, frictionless, and equilibrium is proposed. In the convergent part of the nozzle the flow is treated as single-phase flow of liquid, whereas in the divergent part the flow is treated as separated two-phase flow with saturated condition. The analytical results indicate that the temperature and the pressure decrease rapidly in the divergent part, and the void fraction increases immediately near the throat. Although this analysis is quite simple, the analytical results can follow the experimental results well within this study.展开更多
Based on non-Darcian flow caused by non-Newtonian liquid, the theory of one-dimensional (1D) consolidation was modified to consider variation in the total vertical stress with depth and time. The finite difference met...Based on non-Darcian flow caused by non-Newtonian liquid, the theory of one-dimensional (1D) consolidation was modified to consider variation in the total vertical stress with depth and time. The finite difference method (FDM) was adopted to obtain numerical solutions for excess pore water pressure and average degree of consolidation. When non-Darcian flow is degenerated into Darcian flow, a comparison between numerical solutions and analytical solutions was made to verify reliability of finite difference solutions. Finally, taking into account the ramp time-dependent loading, consolidation behaviors with non-Darcian flow under various parameters were analyzed. Thus, a comprehensive analysis of 1D consolidation combined with non-Darcian flow caused by non-Newtonian liquid was conducted in this paper.展开更多
In order to use the second-order 5-point difference scheme mentioned to compute the solution of one dimension unsteady equations of the direct reflection of the strong plane detonation wave meeting a solid wall barrie...In order to use the second-order 5-point difference scheme mentioned to compute the solution of one dimension unsteady equations of the direct reflection of the strong plane detonation wave meeting a solid wall barrier, in this paper, we technically construct the difference schemes of the boundary and sub-boundary of the problem, and deduce the auto-analogue analytic solutions of the initial value problem, and at the same time, we present a method for the singular property of the initial value problem, from which we can get a satisfactory computation result of this difficult problem.The difference scheme used in this paper to deal with the discontinuity problems of the shock wave are valuable and worth generalization.展开更多
This paper describes a new method of calculation of one-dimensional steady compressible gas flows in channels with possible heat and mass exchange through perforated sidewalls. The channel is divided into small elemen...This paper describes a new method of calculation of one-dimensional steady compressible gas flows in channels with possible heat and mass exchange through perforated sidewalls. The channel is divided into small elements of a finite size for which mass, energy and momentum conservation laws are written in the integral form, assuming linear distribution of the parameters along the length. As a result, the calculation is reduced to finding the roots of a quadratic algebraic equation, thus providing an alternative to numerical methods based on differential equations. The advantage of this method is its high tolerance to coarse discretization of the calculation area as well as its good applicability for transonic flow calculations.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy ...This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy equation is solved with the closed form solutions for prescribed wall temperature and prescribed wall heat flux conditions.The wall temperature and heat flux have power dependence on both time and spatial distance.The solution domain,the velocity distribution,the flow field,and the temperature distribution in the fluids are studied for different controlling parameters.These parameters include the Prandtl number,the mass transfer parameter at the wall,the wall moving parameter,the time power index,and the spatial power index.It is found that two solution branches exist for certain combinations of the controlling parameters for the flow and heat transfer problems.The heat transfer solutions are given by the confluent hypergeometric function of the first kind,which can be simplified into the incomplete gamma functions for special conditions.The wall heat flux and temperature profiles show very complicated variation behaviors.The wall heat flux can have multiple poles under certain given controlling parameters,and the temperature can have significant oscillations with overshoot and negative values in the boundary layers.The relationship between the number of poles in the wall heat flux and the number of zero-crossing points is identified.The difference in the results of the prescribed wall temperature case and the prescribed wall heat flux case is analyzed.Results given in this paper provide a rare closed form analytical solution to the entire unsteady NS equations,which can be used as a benchmark problem for numerical code validation.展开更多
A new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation (SOR) iteration method and Jacobian matrix d...A new time-accurate marching scheme for unsteady flow calculations is proposed in the present work. This method is the combination of classical Successive Over-Relaxation (SOR) iteration method and Jacobian matrix diagonally dominant splitting method of LUSGS. One advantage of this algorithm is the second-order accuracy because of no factorization error. Another advantage is the low computational cost because the Jacobian matrices and fluxes are only calculated once in each physical time step. And, the SOR algorithm has better convergence property than Gauss-Seidel. To investigate its accuracy and convergency, several unsteady flow computational tests are carded out by using the proposed SOR algorithm. Roe's FDS scheme is used to discritize the inviscid flux terms. Unsteady computational results of SOR are compared with the experiment results and those of Gauss-Seidel, Results reveal that the numerical results agree well with the experimental data and the second-order accuracy can be obtained as the Gauss-Seidel for unsteady flow computations. The impact of SOR factor is investigated for unsteady computations by using different SOR factors in this algorithm to simulate each computational test. Different numbers of inner iterations are needed to converge to the same criterion for different SOR factors and optimal choice of SOR factor can improve the computational efficiency greatly.展开更多
The unsteady hydrodynamics of a biomimetic fin attached to a cylindrical body has been studied numerically using a computational fluid dynamic (CFD) simulator based on an in-house solver of the Navier-Stokes equations...The unsteady hydrodynamics of a biomimetic fin attached to a cylindrical body has been studied numerically using a computational fluid dynamic (CFD) simulator based on an in-house solver of the Navier-Stokes equations, combined with a recently developed multi-block, overset grid method. The fin-body CFD model is based on a mechanical pectoral fin device, which consists of a cylindrical body and an asymmetric fin and can mimic flapping, rowing and feathering motions of the pectoral fins in fishes. First the multi-block, overset grid method incorporated into the NS solver was verified through an extensive study of unsteady flows past a single fin undergoing rowing and feathering motion. Then unsteady flows past the biomimetic fin-body model undergoing the same motions were computed and compared with the measurements of forces of the mechanical pectoral fin, which shows good agreement in both time-varying and time-averaged hydrodynamic forces. The relationship between force generation and vortex dynamics points to the importance of the match in fin kinematics between power and recovery strokes and implies that an optimal selection of parameters of phase lags between and amplitudes of rowing and feathering motions can improve the performance of labriform propulsion in terms of either maximum force generation or minimum mechanical power.展开更多
Objective To develop the numerical method for the steady and unsteady wet steam condensing flow in turbine stage. Methods An Eulerian/Eulerian numerical model is used to describe the spontaneous condensation flow in t...Objective To develop the numerical method for the steady and unsteady wet steam condensing flow in turbine stage. Methods An Eulerian/Eulerian numerical model is used to describe the spontaneous condensation flow in the steam turbine. For the steady condensing flow computations, the mixing plane model was used. For the unsteady condensing flow computations, the sliding mesh method was used to simulate the rotor-stator interactions. Results The numerical results showed the obvious differences between non-condensing and condensing flows. The results also showed the unsteadiness effect due to rotor-stator interactions had a deep influence on the formation and growth process of water droplets. Conclusion The numerical methods presented in this paper are valid for the condensing flow in the turbine stage.展开更多
The unsteady laminar magnetohydrodynamics (MHD) boundary layer flow and heat transfer of nanofluids over an accelerating convectively heated stretching sheet are numerically studied in the presence of a transverse m...The unsteady laminar magnetohydrodynamics (MHD) boundary layer flow and heat transfer of nanofluids over an accelerating convectively heated stretching sheet are numerically studied in the presence of a transverse magnetic field with heat source/sink The unsteady governing equations are solved by a shooting method with the Runge-Kutta- Fehlberg scheme. Three different types of water based nanofluids, containing copper, aluminium oxide, and titanium dioxide, are taken into consideration. The effects of the pertinent parameters on the fluid velocity, the temperature, the entropy generation num- ber, the Bejan number, the shear stress, and the heat transfer rate at the sheet surface are graphically and quantitatively discussed in detail. A comparison of the entropy generation due to the heat transfer and the fluid friction is made with the help of the Bejan number. It is observed that the presence of the metallic nanoparticles creates more entropy in the nanofluid flow than in the regular fluid flow.展开更多
The present study refers to a cavitating Venturi type section geometry characterized by a convergent angle of 18°?and a divergent angle of about 8°?where the sheet cavity presents typical self-oscillation be...The present study refers to a cavitating Venturi type section geometry characterized by a convergent angle of 18°?and a divergent angle of about 8°?where the sheet cavity presents typical self-oscillation behavior with quasi-periodic vapor clouds shedding. This work is an extension of previous works concerning void ratio measurements and velocity fields using double optical probe and constitutes a complete analysis of the two-phase structure of unsteady cavitating flow. This paper provides a new method based on conditional and phase averaging technique with wall pressure signal to treat experimental data in order to evaluate more precisely time-averaged and rms values of the void ratio and instantaneous velocity fields. Conditional analysis shows a different behavior of the two-phase flow dynamics leading to highlight high void ratio events linked to the break-off cycle. Unsteady phase averaging of the optical probe signal gives the evolution of the void ratio at each studied location in the venturi and shows that the fluctuations close to the wall (where the re-entrant jet is predominant) are in phase with the upper part of the cavity instead of the thickness of the cavity which is unchanged.展开更多
Direct calculations of unsteady-state Weymouth equations for gas volumetric flow rate occur more frequently in the design and operation analysis of natural gas systems. Most of the existing gas pipelines design proced...Direct calculations of unsteady-state Weymouth equations for gas volumetric flow rate occur more frequently in the design and operation analysis of natural gas systems. Most of the existing gas pipelines design procedures are based on a particular friction factor and steady-state flow analysis. This paper examined the behavior of different friction factors and the need to develop model analysis capable of calculating unsteady-state gas flow rate in horizontal and inclined pipes. The results show different variation in flow rate with Panhandle A and Panhandle B attaining stability in accurate time with initial unsteadiness at the instance of flow. Chen and Jain friction factors have opposition to flow with high flow rate: The prediction also reveals that Colebrook-White degenerated to Nikuradse friction factor at high Reynolds number. The horizontal and inclined flow equations are considerably enhanced on the usage of different friction factors with the aid of Matlab to handle these calculations.展开更多
文摘The aim of this study is to investigate CO2 two-phase nozzle flow in terms of both experimental and analytical aspects for the optimum design of two-phase flow nozzle of CO2 two-phase flow ejector. In the experiment, it is measured that the temperature profile in the stream-wise direction of a divergent-convergent nozzle through which CO2 in the supercritical pressure condition is blown down into the atmosphere. In the analysis, a one-dimensional model which assumes steady, adiabatic, frictionless, and equilibrium is proposed. In the convergent part of the nozzle the flow is treated as single-phase flow of liquid, whereas in the divergent part the flow is treated as separated two-phase flow with saturated condition. The analytical results indicate that the temperature and the pressure decrease rapidly in the divergent part, and the void fraction increases immediately near the throat. Although this analysis is quite simple, the analytical results can follow the experimental results well within this study.
基金Supported by the National Natural Science Foundation of China (51109092,50878191)
文摘Based on non-Darcian flow caused by non-Newtonian liquid, the theory of one-dimensional (1D) consolidation was modified to consider variation in the total vertical stress with depth and time. The finite difference method (FDM) was adopted to obtain numerical solutions for excess pore water pressure and average degree of consolidation. When non-Darcian flow is degenerated into Darcian flow, a comparison between numerical solutions and analytical solutions was made to verify reliability of finite difference solutions. Finally, taking into account the ramp time-dependent loading, consolidation behaviors with non-Darcian flow under various parameters were analyzed. Thus, a comprehensive analysis of 1D consolidation combined with non-Darcian flow caused by non-Newtonian liquid was conducted in this paper.
文摘In order to use the second-order 5-point difference scheme mentioned to compute the solution of one dimension unsteady equations of the direct reflection of the strong plane detonation wave meeting a solid wall barrier, in this paper, we technically construct the difference schemes of the boundary and sub-boundary of the problem, and deduce the auto-analogue analytic solutions of the initial value problem, and at the same time, we present a method for the singular property of the initial value problem, from which we can get a satisfactory computation result of this difficult problem.The difference scheme used in this paper to deal with the discontinuity problems of the shock wave are valuable and worth generalization.
文摘This paper describes a new method of calculation of one-dimensional steady compressible gas flows in channels with possible heat and mass exchange through perforated sidewalls. The channel is divided into small elements of a finite size for which mass, energy and momentum conservation laws are written in the integral form, assuming linear distribution of the parameters along the length. As a result, the calculation is reduced to finding the roots of a quadratic algebraic equation, thus providing an alternative to numerical methods based on differential equations. The advantage of this method is its high tolerance to coarse discretization of the calculation area as well as its good applicability for transonic flow calculations.
基金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.
基金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.
基金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.
文摘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.
基金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.
文摘This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy equation is solved with the closed form solutions for prescribed wall temperature and prescribed wall heat flux conditions.The wall temperature and heat flux have power dependence on both time and spatial distance.The solution domain,the velocity distribution,the flow field,and the temperature distribution in the fluids are studied for different controlling parameters.These parameters include the Prandtl number,the mass transfer parameter at the wall,the wall moving parameter,the time power index,and the spatial power index.It is found that two solution branches exist for certain combinations of the controlling parameters for the flow and heat transfer problems.The heat transfer solutions are given by the confluent hypergeometric function of the first kind,which can be simplified into the incomplete gamma functions for special conditions.The wall heat flux and temperature profiles show very complicated variation behaviors.The wall heat flux can have multiple poles under certain given controlling parameters,and the temperature can have significant oscillations with overshoot and negative values in the boundary layers.The relationship between the number of poles in the wall heat flux and the number of zero-crossing points is identified.The difference in the results of the prescribed wall temperature case and the prescribed wall heat flux case is analyzed.Results given in this paper provide a rare closed form analytical solution to the entire unsteady NS equations,which can be used as a benchmark problem for numerical code validation.
基金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 unsteady hydrodynamics of a biomimetic fin attached to a cylindrical body has been studied numerically using a computational fluid dynamic (CFD) simulator based on an in-house solver of the Navier-Stokes equations, combined with a recently developed multi-block, overset grid method. The fin-body CFD model is based on a mechanical pectoral fin device, which consists of a cylindrical body and an asymmetric fin and can mimic flapping, rowing and feathering motions of the pectoral fins in fishes. First the multi-block, overset grid method incorporated into the NS solver was verified through an extensive study of unsteady flows past a single fin undergoing rowing and feathering motion. Then unsteady flows past the biomimetic fin-body model undergoing the same motions were computed and compared with the measurements of forces of the mechanical pectoral fin, which shows good agreement in both time-varying and time-averaged hydrodynamic forces. The relationship between force generation and vortex dynamics points to the importance of the match in fin kinematics between power and recovery strokes and implies that an optimal selection of parameters of phase lags between and amplitudes of rowing and feathering motions can improve the performance of labriform propulsion in terms of either maximum force generation or minimum mechanical power.
基金This work was supported by the National Natural Science Foundation of China (No50336050)
文摘Objective To develop the numerical method for the steady and unsteady wet steam condensing flow in turbine stage. Methods An Eulerian/Eulerian numerical model is used to describe the spontaneous condensation flow in the steam turbine. For the steady condensing flow computations, the mixing plane model was used. For the unsteady condensing flow computations, the sliding mesh method was used to simulate the rotor-stator interactions. Results The numerical results showed the obvious differences between non-condensing and condensing flows. The results also showed the unsteadiness effect due to rotor-stator interactions had a deep influence on the formation and growth process of water droplets. Conclusion The numerical methods presented in this paper are valid for the condensing flow in the turbine stage.
文摘The unsteady laminar magnetohydrodynamics (MHD) boundary layer flow and heat transfer of nanofluids over an accelerating convectively heated stretching sheet are numerically studied in the presence of a transverse magnetic field with heat source/sink The unsteady governing equations are solved by a shooting method with the Runge-Kutta- Fehlberg scheme. Three different types of water based nanofluids, containing copper, aluminium oxide, and titanium dioxide, are taken into consideration. The effects of the pertinent parameters on the fluid velocity, the temperature, the entropy generation num- ber, the Bejan number, the shear stress, and the heat transfer rate at the sheet surface are graphically and quantitatively discussed in detail. A comparison of the entropy generation due to the heat transfer and the fluid friction is made with the help of the Bejan number. It is observed that the presence of the metallic nanoparticles creates more entropy in the nanofluid flow than in the regular fluid flow.
文摘The present study refers to a cavitating Venturi type section geometry characterized by a convergent angle of 18°?and a divergent angle of about 8°?where the sheet cavity presents typical self-oscillation behavior with quasi-periodic vapor clouds shedding. This work is an extension of previous works concerning void ratio measurements and velocity fields using double optical probe and constitutes a complete analysis of the two-phase structure of unsteady cavitating flow. This paper provides a new method based on conditional and phase averaging technique with wall pressure signal to treat experimental data in order to evaluate more precisely time-averaged and rms values of the void ratio and instantaneous velocity fields. Conditional analysis shows a different behavior of the two-phase flow dynamics leading to highlight high void ratio events linked to the break-off cycle. Unsteady phase averaging of the optical probe signal gives the evolution of the void ratio at each studied location in the venturi and shows that the fluctuations close to the wall (where the re-entrant jet is predominant) are in phase with the upper part of the cavity instead of the thickness of the cavity which is unchanged.
文摘Direct calculations of unsteady-state Weymouth equations for gas volumetric flow rate occur more frequently in the design and operation analysis of natural gas systems. Most of the existing gas pipelines design procedures are based on a particular friction factor and steady-state flow analysis. This paper examined the behavior of different friction factors and the need to develop model analysis capable of calculating unsteady-state gas flow rate in horizontal and inclined pipes. The results show different variation in flow rate with Panhandle A and Panhandle B attaining stability in accurate time with initial unsteadiness at the instance of flow. Chen and Jain friction factors have opposition to flow with high flow rate: The prediction also reveals that Colebrook-White degenerated to Nikuradse friction factor at high Reynolds number. The horizontal and inclined flow equations are considerably enhanced on the usage of different friction factors with the aid of Matlab to handle these calculations.