A global optimization approach to turbine blade design based on hierarchical fair competition genetic algorithms with dynamic niche （HFCDN-GAs） coupled with Reynolds-averaged Navier-Stokes （RANS） equation is prese...A global optimization approach to turbine blade design based on hierarchical fair competition genetic algorithms with dynamic niche （HFCDN-GAs） coupled with Reynolds-averaged Navier-Stokes （RANS） equation is presented. In order to meet the search theory of GAs and the aerodynamic performances of turbine, Bezier curve is adopted to parameterize the turbine blade profile, and a fitness function pertaining to optimization is designed. The design variables are the control points＇ ordinates of characteristic polygon of Bezier curve representing the turbine blade profile. The object function is the maximum lift-drag ratio of the turbine blade. The constraint conditions take into account the leading and trailing edge metal angle, and the strength and aerodynamic performances of turbine blade. And the treatment method of the constraint conditions is the flexible penalty function. The convergence history of test function indicates that HFCDN-GAs can locate the global optimum within a few search steps and have high robustness. The lift-drag ratio of the optimized blade is 8.3% higher than that of the original one. The results show that the proposed global optimization approach is effective for turbine blade.展开更多
A hybrid approach coupled with a surface panel method for the propeller and a Reynolds averaged Navier-Stokes（RANS） model for the hull with the propeller body forces are presented for predicting the self-propulsion ...A hybrid approach coupled with a surface panel method for the propeller and a Reynolds averaged Navier-Stokes（RANS） model for the hull with the propeller body forces are presented for predicting the self-propulsion performance and the effective wake field of underwater vehicles. To achieve a high accuracy and simplicity, a radial basis function（RBF） based approach is proposed for mapping the force field from the blade surface panels to the RANS model. The effective wake field is evaluated in two ways, i.e., by extrapolation from the flat planes upstream of the propeller disk, and by direct computation in a curved surface upstream of and parallel to the blade leading edges. The hull-propeller system of a real propeller geometry is further simulated with the sliding mesh model to numerically verify the hybrid approach. Numerical simulations are conducted for the fully appended SUBOFF submarine model. The high accuracy of the RBF-based interpolation scheme is confirmed, and the effective wake fraction predicted by the hybrid approach is found consistent with that obtained by the sliding mesh model. The effective wake fractions predicted by the two methods are, respectively, 4.6% and 3% larger than the nominal one.展开更多
A numerical approach based on the solution of the Reynolds-averaged Navier-Stokes(RANS) equations using the shear-stress transport(SST) turbulence model has been employed to investigate the hydrodynamic performance an...A numerical approach based on the solution of the Reynolds-averaged Navier-Stokes(RANS) equations using the shear-stress transport(SST) turbulence model has been employed to investigate the hydrodynamic performance and flow of tunnel thrusters.The flow passages between adjacent blades are discretized with prismatic cells so that the boundary layer flow is resolved down to the viscous sub-layer.The hydrodynamic performances predicted by the quasi-steady approach agree well with the experimental data for three impellers covering a range of blade area and pitch.Through analysis of the flow field,the reason why the hub of impeller also contributes to thrust which can amount to 40%—60% of the impeller thrust,and the mechanism of the impeller inducing an axial force on the hull are elucidated.展开更多
The roll motions are influenced by significant viscous effects such as the flow separation.The 3D simulations of free decay roll motions for the ship model DTMB 5512 are carried out by Reynold averaged NavierStokes(RA...The roll motions are influenced by significant viscous effects such as the flow separation.The 3D simulations of free decay roll motions for the ship model DTMB 5512 are carried out by Reynold averaged NavierStokes(RANS) method based on the dynamic mesh technique.A new moving mesh technique is adopted and discussed in details for the present simulations.The purpose of the research is to obtain accurate numerical prediction for roll motions with their respective numerical/modeling errors and uncertainties.Errors and uncertainties are estimated by performing the modern verification and validation(V&V) procedures.Simulation results for the free-floating surface combatant are used to calculate the linear,nonlinear damping coefficients and resonant frequencies including a wide range of forward speed.The present work can provide a useful reference to calculate roll damping by computational fluid dynamics(CFD) method and simulate a general ship motions in waves.展开更多
In this study,the performance of a twin-screw propeller under the influence of the wake field of a fully appended ship was investigated using a coupled Reynolds-averaged Navier–Stokes(RANS)/boundary element method(BE...In this study,the performance of a twin-screw propeller under the influence of the wake field of a fully appended ship was investigated using a coupled Reynolds-averaged Navier–Stokes(RANS)/boundary element method(BEM)code.The unsteady BEM is an efficient approach to predicting propeller performance.By applying the time-stepping method in the BEM solver,the trailing vortex sheet pattern of the propeller can be accurately captured at each time step.This is the main innovation of the coupled strategy.Furthermore,to ascertain the effect of the wake field of the ship with acceptable accuracy,a RANS solver was developed.A finite volume method was used to discretize the Navier–Stokes equations on fully unstructured grids.To simulate ship motions,the volume of the fluid method was applied to the RANS solver.The validation of each solver(BEM/RANS)was separately performed,and the results were compared with experimental data.Ultimately,the BEM and RANS solvers were coupled to estimate the performance of a twin-screw propeller,which was affected by the wake field of the fully appended hull.The proposed model was applied to a twin-screw oceanography research vessel.The results demonstrated that the presented model can estimate the thrust coefficient of a propeller with good accuracy as compared to an experimental self-propulsion test.The wake sheet pattern of the propeller in open water(uniform flow)was also compared with the propeller in a real wake field.展开更多
Reynolds-Averaged Navier-Stokes(RANS) Computational Fluid Dynamics(CFD) has been widely used in compressor design and analysis. However, reasonable prediction of compressor flow and its impact on compressor performanc...Reynolds-Averaged Navier-Stokes(RANS) Computational Fluid Dynamics(CFD) has been widely used in compressor design and analysis. However, reasonable prediction of compressor flow and its impact on compressor performance remains challenging. In this study, Menter’s Shear Stress Transport(SST) model and its variants, as well as the ω-based Reynolds stress model(Stress-BSL) are assessed. For a single rotor(Rotor 67), under the peak efficiency operating condition, all studied turbulence models predict its performance with reasonable accuracy;under the off-design conditions, SST with Helicity correction(SST-Helicity) shows superiority in predicting the effect of flow on the spanwise distribution of aerodynamic parameters. For Darmstadt’s 1.5-stage transonic axial compressor, SST-Helicity outperforms SST, SST with the Quadratic Constitutive Relation(SST-QCR) and Stress-BSL in predicting the performance as well as the spanwise distribution of aerodynamic parameters. At the design rotating speed, the stall margin given by SST-Helicity(20.90%) is the closest to the experimental measurement(24.81%), which is more than twice that by SST(8.71%) and 1.5 times that by SST-QCR(14.14%). This paper demonstrates that SSTHelicity model, together with a good quality and sufficiently refined grid, can capture the compressor flow features with reasonable accuracy, which results in a credible prediction of compressor performance and stage matching.展开更多
Rotor-stator cavities are frequently encountered in engineering applications such as gas turbine engines.They are usually subject to an external hot mainstream crossflow which in general is highly swirled under the ef...Rotor-stator cavities are frequently encountered in engineering applications such as gas turbine engines.They are usually subject to an external hot mainstream crossflow which in general is highly swirled under the effect of the nozzle guide vanes.To avoid hot mainstream gas ingress,the cavity is usually purged by a stream of sealing flow.The interactions between the external crossflow,cavity flow,and sealing flow are complicated and involve all scales of turbulent unsteadiness and flow instability which are beyond the resolution of the Reynolds-average approach.To cope with such a complex issue,a wall-modeled large-eddy simulation(WMLES)approach is adopted in this study.In the simulation,a 20°sector model is used and subjected to a uniform pre-swirled external crossflow and a stream of radial sealing flow.It is triggered by a convergent Reynoldsaveraged Navier-Stokes(RANS)result in which the shear stress transport(SST)turbulent model is used.In the WMLES simulation,the Smagoringsky sub-grid scale(SGS)model is applied.A scalar transportation equation is solved to simulate the blending and transportation process in the cavity.The overall flow field characteristics and deviation between RANS and WMLES results are discussed first.Both RANS and WMLES results show a Batchelor flow mode,while distinct deviation is also observed.Deviations in the small-radius region are caused by the insufficiency of the RANS approach in capturing the small-scale vortex structures in the boundary layer while deviations in the large-radius region are caused by the insufficiency of the RANS approach in predicting the external crossflow ingestion.The boundary layer vortex and external ingestion are then discussed in detail,highlighting the related flow instabilities.Finally,the large-flow structures induced by external flow ingress are analyzed using unsteady pressure oscillation signals.展开更多
Different turbulence closures were used to predict the flow interaction between the wakes created by compressor outlet guide vanes(OGVs) and a downstream annular pre-diffuser.Two statistical turbulence models were tes...Different turbulence closures were used to predict the flow interaction between the wakes created by compressor outlet guide vanes(OGVs) and a downstream annular pre-diffuser.Two statistical turbulence models were tested based on the classical Reynolds-averaged Navier-Stokes(RANS) approach.Both high-Re and low-Re(Launder-Sharma) versions of the k-ε model were applied to a selected test problem for OGV wake/diffuser flows.The test problem was specifically chosen because experimentally determined inlet conditions and both profile and performance data were available to validate predictions.A preliminary study was also reported of the more advanced large eddy simulation(LES) approach.The LES sub-grid-scale(SGS) model was the basic Smagorinsky eddy viscosity assumption,with a Van-Driest damping function for improved capture of near-wall viscous behaviour.Comparison between the two RANS models showed little difference in terms of velocity contours at OGV trailing edge and diffuser exit.In terms of overall diffuser performance(static pressure recovery and total pressure loss coefficients),the high-Re model was shown to agree well with experimental data.The preliminary LES study indicates the highly unsteady character of the OGV wake flow,but requires improved treatment of inlet conditions.展开更多
On the basis of the transient water wave(TWW) theory,focused wave is generated in the circulating water channel.Numerical simulation of the focused wave is carried out by solving the Reynolds averaged NavierStokes(RAN...On the basis of the transient water wave(TWW) theory,focused wave is generated in the circulating water channel.Numerical simulation of the focused wave is carried out by solving the Reynolds averaged NavierStokes(RANS) equations.The dynamic grid technique is adopted to simulate the motion of the wave maker,and the volume of fluid(VOF) method is used to capture the free surface of the wave.The simulation results are compared with the measured data,and good agreement is obtained.For quantitative estimation of the numerical simulation error and uncertainty,the uncertainty analysis method recommended by the International Towing Tank Conference(ITTC) procedure is performed for the simulation results of the surface elevations at different positions.Both grid-convergence and time-step-size convergence studies are conducted using three types of grids and time step sizes.The simulation results are all monotonously convergent in the verification procedure,and the validations of the simulated surface elevations with the positions at 3.5,4.0 and 4.5 m are all achieved by comparing with the validation uncertainty.It is found that the numerical simulation errors caused by the grid and time-step-size in the convergence studies have the same order of magnitude.In addition,the numerical errors and uncertainties for the surface elevations at different positions are compared and discussed in detail.This paper presents the first attempt to carry out the uncertainty analysis of the simulation of focused wave,and the effectiveness of the proposed verification and validation procedures in the uncertainty analysis is demonstrated.展开更多
Shock wave-boundary layer interactions(SWBLI)are observed in several practical high-speed internal flows,such as compressor blades,turbine cascades,nozzles and so on.Shock induced oscillations(SIO),aerodynamic instabi...Shock wave-boundary layer interactions(SWBLI)are observed in several practical high-speed internal flows,such as compressor blades,turbine cascades,nozzles and so on.Shock induced oscillations(SIO),aerodynamic instabilities so-called buffet flows,flutter,aeroacoustic noise and vibration are the detrimental consequences of this unsteady shockboundary layer interactions.In the present study,a numerical computation has been performed to investigate the compressible flow characteristics around a 12%thick biconvex circular arc airfoil in a two dimensional channel.Reynolds averaged Navier-Stokes equations with two equation k-ωshear stress transport(SST)turbulence model have been applied for the computational analysis.The flow field characteristics has been studied from pressure ratio(ratio of back pressure,pb to inlet total pressure,p01)of 0.75 to 0.65.The present computational results have been compared and validated with the available experimental data.The results showed that the internal flow field characteristics such as shock wave structure,its behavior(steady or unsteady)and the corresponding boundary layer interaction are varied with pressure ratio.Self-excited shock oscillation was observed at certain flow conditions.Moreover,the mode of unsteady shock oscillation and its frequency are varied significantly with change of pressure ratio.展开更多
Numerical works have been conducted to investigate the effect of nozzle geometries on the discharge coefficient.Several contoured converging nozzles with finite radius of curvatures,conically converging nozzles and co...Numerical works have been conducted to investigate the effect of nozzle geometries on the discharge coefficient.Several contoured converging nozzles with finite radius of curvatures,conically converging nozzles and conical divergent orifices have been employed in this investigation.Each nozzle and orifice has a nominal exit diameter of 12.7x10^(-3)m.A 3rd order MUSCL finite volume method of ANSYS Fluent 13.0 was used to solve the Reynolds-averaged Navier-Stokes equations in simulating turbulent flows through various nozzle inlet geometries.The numerical model was validated through comparison between the numerical results and experimental data.The results obtained show that the nozzle geometry has pronounced effect on the sonic lines and discharge coefficients.The coefficient of discharge was found differ from unity due to the non-uniformity of flow parameters at the nozzle exit and the presence of boundary layer as well.展开更多
Turbofan engine intakes are designed to provide separation-free flow at the fan faceover a wide range of operating conditions. But at some off-design conditions, like at high flightspeeds and high angles of attack (Ao...Turbofan engine intakes are designed to provide separation-free flow at the fan faceover a wide range of operating conditions. But at some off-design conditions, like at high flightspeeds and high angles of attack (AoA), the aero engine intake may encounter flow separation.This boundary layer separation inside the nacelle inlet of an aircraft engine can lead to a largenumber of undesirable outcomes like reduction in fan efficiency, engine stall and high levels ofstress on the fan blades. Active flow control is a promising solution to reduce inlet boundarylayer separation and the associated fan-face flow distortion at such off-design conditions. Byblowing pressurized air into the intake near the separation point, the boundary layer is ener-gized and separation can be controlled. This study investigates the applicability of lip blowing,an active flow control technique, to control intake separation and flow distortion at the fan-face.First, intake separation was triggered in a 3D CFD model based on the NASA CommonResearch Model (CRM) using high AoA cases at cruise condition (Mach number 0.85, Massflow capture ratio w0.7) and the features of separated flow were analyzed. Thereafter, activeflow control was introduce to the intake in the form of two types of lip blowing, direct andpitched blowing. The efficacy of lip blowing at achieving separation control in an ultra highbypass ratio turbofan engine intake has been established through this study. The present paperalso examines the significance of blowing parameters like the type of blowing, blowing pres-sure ratio, and blowing slot dimension, at different angles of attack to identify the critical con-trol parameters. Our research successfully establishes proof of concept by demonstrating the feasibility of using lip blowing for separation control in aero-intakes, via numerical modelling.Furthermore, this study also provides crucial insights regarding the important variables to beconsidered for future experimental studies, and also for detailed studies covering a wider展开更多
Conical plug nozzle and truncated conical plug nozzle are advanced rocket nozzles suitable for use as altitude compensating nozzles.In this study flow through the conical plug and truncated conical plug nozzles are nu...Conical plug nozzle and truncated conical plug nozzle are advanced rocket nozzles suitable for use as altitude compensating nozzles.In this study flow through the conical plug and truncated conical plug nozzles are numerically simulated to first validate with experimental data and then to compare the performance when a base bleed is introduced.The numerical analysis has considered two-dimensional axisymmetric models.Reynolds-averaged NavierStokes equations are solved with two equation shear stress transport k-ω turbulence model.For the validation of the plug nozzle,flow features and wall pressure along the length of the nozzle is taken for different nozzle pressure ratios.For the validation of truncated plug nozzle,flow features and base pressures at various nozzle pressure ratios are compared.The base bleed is taken as 2%of the inlet mass flow rate.The comparison of results shows that the introduction of base bleed helps to compensate for the loss of thrust due to conical plug nozzle truncation.展开更多
Effects of spontaneous condensation of moist air on the shock wave dynamics around butterfly valves in transonic flows are investigated by experimental and numerical simulations.Two symmetric valve disk shapes namely-...Effects of spontaneous condensation of moist air on the shock wave dynamics around butterfly valves in transonic flows are investigated by experimental and numerical simulations.Two symmetric valve disk shapes namely-a flat rectangular plate and a mid-plane cross-section of a prototype butterfly valve have been studied in the present research.Results showed that in case with spontaneous condensation,the root mean square of pressure oscillation(induced by shock dynamics)is reduced significantly with those without condensation for both shapes of the valves.Moreover,local aerodynamic moments were reduced in case with condensation which is considered to be beneficial in torque requirement in case of on/off applications of valves as flow control devices.However,total pressure loss was increased with spontaneous condensation in both the valves.Furthermore,the disk shape of a prototype butterfly valve showed better aerodynamic performances compared to flat rectangular plate profile in respect of total pressure loss and vortex shedding frequency in the wake region.展开更多
This article presents an approach which employs a commercial Reynolds-Averaged Navier-Stokes（RANS）solver to predict the steady wake field and loading distributions for a rim driven thruster.Four different cases of p...This article presents an approach which employs a commercial Reynolds-Averaged Navier-Stokes（RANS）solver to predict the steady wake field and loading distributions for a rim driven thruster.Four different cases of propeller blades are chosen to be calculated with the presented method.The propeller blade radial circulation and chordwise circulation density distributions are analyzed.The maximum radial circulation is found at the blade tip,which is different from conventional shaft driven propeller.The numerical results indicate that there is no tip leakage vortex in rim driven propulors.But there exist the tip joint vortex and the root region vortex.Bollard characteristics are calculated by taking rim surface effect into account.From the predicted results the second case in this paper is selected as the final one to perform hydrodynamic experiment.The calculation results with empirical rim surface corrections are compared with the measurement.It shows that the developed numerical method can well predict hydrodynamic performances of the rim driven thruster.展开更多
The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the...The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.展开更多
The separated turbulent flow around a circular cylinder is investigated using Large-Eddy Simulation （LES）, Detached-Eddy Simulation （DES, or hybrid RANS/LES methods）, and Unsteady Reynolds-Averaged Navier-Stokes ...The separated turbulent flow around a circular cylinder is investigated using Large-Eddy Simulation （LES）, Detached-Eddy Simulation （DES, or hybrid RANS/LES methods）, and Unsteady Reynolds-Averaged Navier-Stokes （URANS）. The purpose of this study is to examine some typical simulation approaches for the prediction of complex separated turbulent flow and to clarify the capability of applying these approaches to a typical case of the separated turbulent flow around a circular cylinder. Several turbulence models, i.e. dynamic Sub-grid Scale （SGS） model in LES, the DES-based Spalart-Allmaras （S-A） and κ-ω Shear-Stress- Transport （SST） models in DES, and the S-A and SST models in URANS, are used in the calculations. Some typical results, e.g., the mean pressure and drag coefficients, velocity profiles, Strouhal number, and Reynolds stresses, are obtained and compared with previous computational and experimental data. Based on our extensive calculations, we assess the capability and performance of these simulation approaches coupled with the relevant turbulence models to predict the separated turbulent flow.展开更多
基金This project is supported by National Natural Science Foundation of China (No,50776056)National Hi-tech Research and Development Program of China (863 Program,No.2006AA05Z250).
文摘A global optimization approach to turbine blade design based on hierarchical fair competition genetic algorithms with dynamic niche （HFCDN-GAs） coupled with Reynolds-averaged Navier-Stokes （RANS） equation is presented. In order to meet the search theory of GAs and the aerodynamic performances of turbine, Bezier curve is adopted to parameterize the turbine blade profile, and a fitness function pertaining to optimization is designed. The design variables are the control points＇ ordinates of characteristic polygon of Bezier curve representing the turbine blade profile. The object function is the maximum lift-drag ratio of the turbine blade. The constraint conditions take into account the leading and trailing edge metal angle, and the strength and aerodynamic performances of turbine blade. And the treatment method of the constraint conditions is the flexible penalty function. The convergence history of test function indicates that HFCDN-GAs can locate the global optimum within a few search steps and have high robustness. The lift-drag ratio of the optimized blade is 8.3% higher than that of the original one. The results show that the proposed global optimization approach is effective for turbine blade.
基金Project supported by the National Basic Research Development Program of China(973 Program,Grant No.613134)
文摘A hybrid approach coupled with a surface panel method for the propeller and a Reynolds averaged Navier-Stokes（RANS） model for the hull with the propeller body forces are presented for predicting the self-propulsion performance and the effective wake field of underwater vehicles. To achieve a high accuracy and simplicity, a radial basis function（RBF） based approach is proposed for mapping the force field from the blade surface panels to the RANS model. The effective wake field is evaluated in two ways, i.e., by extrapolation from the flat planes upstream of the propeller disk, and by direct computation in a curved surface upstream of and parallel to the blade leading edges. The hull-propeller system of a real propeller geometry is further simulated with the sliding mesh model to numerically verify the hybrid approach. Numerical simulations are conducted for the fully appended SUBOFF submarine model. The high accuracy of the RBF-based interpolation scheme is confirmed, and the effective wake fraction predicted by the hybrid approach is found consistent with that obtained by the sliding mesh model. The effective wake fractions predicted by the two methods are, respectively, 4.6% and 3% larger than the nominal one.
文摘A numerical approach based on the solution of the Reynolds-averaged Navier-Stokes(RANS) equations using the shear-stress transport(SST) turbulence model has been employed to investigate the hydrodynamic performance and flow of tunnel thrusters.The flow passages between adjacent blades are discretized with prismatic cells so that the boundary layer flow is resolved down to the viscous sub-layer.The hydrodynamic performances predicted by the quasi-steady approach agree well with the experimental data for three impellers covering a range of blade area and pitch.Through analysis of the flow field,the reason why the hub of impeller also contributes to thrust which can amount to 40%—60% of the impeller thrust,and the mechanism of the impeller inducing an axial force on the hull are elucidated.
基金the National Natural Science Foundation of China(No.51579147)
文摘The roll motions are influenced by significant viscous effects such as the flow separation.The 3D simulations of free decay roll motions for the ship model DTMB 5512 are carried out by Reynold averaged NavierStokes(RANS) method based on the dynamic mesh technique.A new moving mesh technique is adopted and discussed in details for the present simulations.The purpose of the research is to obtain accurate numerical prediction for roll motions with their respective numerical/modeling errors and uncertainties.Errors and uncertainties are estimated by performing the modern verification and validation(V&V) procedures.Simulation results for the free-floating surface combatant are used to calculate the linear,nonlinear damping coefficients and resonant frequencies including a wide range of forward speed.The present work can provide a useful reference to calculate roll damping by computational fluid dynamics(CFD) method and simulate a general ship motions in waves.
文摘In this study,the performance of a twin-screw propeller under the influence of the wake field of a fully appended ship was investigated using a coupled Reynolds-averaged Navier–Stokes(RANS)/boundary element method(BEM)code.The unsteady BEM is an efficient approach to predicting propeller performance.By applying the time-stepping method in the BEM solver,the trailing vortex sheet pattern of the propeller can be accurately captured at each time step.This is the main innovation of the coupled strategy.Furthermore,to ascertain the effect of the wake field of the ship with acceptable accuracy,a RANS solver was developed.A finite volume method was used to discretize the Navier–Stokes equations on fully unstructured grids.To simulate ship motions,the volume of the fluid method was applied to the RANS solver.The validation of each solver(BEM/RANS)was separately performed,and the results were compared with experimental data.Ultimately,the BEM and RANS solvers were coupled to estimate the performance of a twin-screw propeller,which was affected by the wake field of the fully appended hull.The proposed model was applied to a twin-screw oceanography research vessel.The results demonstrated that the presented model can estimate the thrust coefficient of a propeller with good accuracy as compared to an experimental self-propulsion test.The wake sheet pattern of the propeller in open water(uniform flow)was also compared with the propeller in a real wake field.
文摘Reynolds-Averaged Navier-Stokes(RANS) Computational Fluid Dynamics(CFD) has been widely used in compressor design and analysis. However, reasonable prediction of compressor flow and its impact on compressor performance remains challenging. In this study, Menter’s Shear Stress Transport(SST) model and its variants, as well as the ω-based Reynolds stress model(Stress-BSL) are assessed. For a single rotor(Rotor 67), under the peak efficiency operating condition, all studied turbulence models predict its performance with reasonable accuracy;under the off-design conditions, SST with Helicity correction(SST-Helicity) shows superiority in predicting the effect of flow on the spanwise distribution of aerodynamic parameters. For Darmstadt’s 1.5-stage transonic axial compressor, SST-Helicity outperforms SST, SST with the Quadratic Constitutive Relation(SST-QCR) and Stress-BSL in predicting the performance as well as the spanwise distribution of aerodynamic parameters. At the design rotating speed, the stall margin given by SST-Helicity(20.90%) is the closest to the experimental measurement(24.81%), which is more than twice that by SST(8.71%) and 1.5 times that by SST-QCR(14.14%). This paper demonstrates that SSTHelicity model, together with a good quality and sufficiently refined grid, can capture the compressor flow features with reasonable accuracy, which results in a credible prediction of compressor performance and stage matching.
基金This work is supported by the National Natural Science Foundation of China(No.5212201273)the National Science and Technology Major Project of China(No.J2019-III-0003)The CFX software and computation resource supplied by Beijing Super Cloud Computing Center,China are acknowledged.
文摘Rotor-stator cavities are frequently encountered in engineering applications such as gas turbine engines.They are usually subject to an external hot mainstream crossflow which in general is highly swirled under the effect of the nozzle guide vanes.To avoid hot mainstream gas ingress,the cavity is usually purged by a stream of sealing flow.The interactions between the external crossflow,cavity flow,and sealing flow are complicated and involve all scales of turbulent unsteadiness and flow instability which are beyond the resolution of the Reynolds-average approach.To cope with such a complex issue,a wall-modeled large-eddy simulation(WMLES)approach is adopted in this study.In the simulation,a 20°sector model is used and subjected to a uniform pre-swirled external crossflow and a stream of radial sealing flow.It is triggered by a convergent Reynoldsaveraged Navier-Stokes(RANS)result in which the shear stress transport(SST)turbulent model is used.In the WMLES simulation,the Smagoringsky sub-grid scale(SGS)model is applied.A scalar transportation equation is solved to simulate the blending and transportation process in the cavity.The overall flow field characteristics and deviation between RANS and WMLES results are discussed first.Both RANS and WMLES results show a Batchelor flow mode,while distinct deviation is also observed.Deviations in the small-radius region are caused by the insufficiency of the RANS approach in capturing the small-scale vortex structures in the boundary layer while deviations in the large-radius region are caused by the insufficiency of the RANS approach in predicting the external crossflow ingestion.The boundary layer vortex and external ingestion are then discussed in detail,highlighting the related flow instabilities.Finally,the large-flow structures induced by external flow ingress are analyzed using unsteady pressure oscillation signals.
文摘Different turbulence closures were used to predict the flow interaction between the wakes created by compressor outlet guide vanes(OGVs) and a downstream annular pre-diffuser.Two statistical turbulence models were tested based on the classical Reynolds-averaged Navier-Stokes(RANS) approach.Both high-Re and low-Re(Launder-Sharma) versions of the k-ε model were applied to a selected test problem for OGV wake/diffuser flows.The test problem was specifically chosen because experimentally determined inlet conditions and both profile and performance data were available to validate predictions.A preliminary study was also reported of the more advanced large eddy simulation(LES) approach.The LES sub-grid-scale(SGS) model was the basic Smagorinsky eddy viscosity assumption,with a Van-Driest damping function for improved capture of near-wall viscous behaviour.Comparison between the two RANS models showed little difference in terms of velocity contours at OGV trailing edge and diffuser exit.In terms of overall diffuser performance(static pressure recovery and total pressure loss coefficients),the high-Re model was shown to agree well with experimental data.The preliminary LES study indicates the highly unsteady character of the OGV wake flow,but requires improved treatment of inlet conditions.
基金the National Natural Science Foundation of China(No.51379124)the Knowledge-Based Ship-Design Hyper-Integrated Platform(KSHIP) 2nd Term(No.GKZY010004)
文摘On the basis of the transient water wave(TWW) theory,focused wave is generated in the circulating water channel.Numerical simulation of the focused wave is carried out by solving the Reynolds averaged NavierStokes(RANS) equations.The dynamic grid technique is adopted to simulate the motion of the wave maker,and the volume of fluid(VOF) method is used to capture the free surface of the wave.The simulation results are compared with the measured data,and good agreement is obtained.For quantitative estimation of the numerical simulation error and uncertainty,the uncertainty analysis method recommended by the International Towing Tank Conference(ITTC) procedure is performed for the simulation results of the surface elevations at different positions.Both grid-convergence and time-step-size convergence studies are conducted using three types of grids and time step sizes.The simulation results are all monotonously convergent in the verification procedure,and the validations of the simulated surface elevations with the positions at 3.5,4.0 and 4.5 m are all achieved by comparing with the validation uncertainty.It is found that the numerical simulation errors caused by the grid and time-step-size in the convergence studies have the same order of magnitude.In addition,the numerical errors and uncertainties for the surface elevations at different positions are compared and discussed in detail.This paper presents the first attempt to carry out the uncertainty analysis of the simulation of focused wave,and the effectiveness of the proposed verification and validation procedures in the uncertainty analysis is demonstrated.
基金The present work has been carried out with computa-tional resource support from Higher Education Quality Enhancement Project(HEQEP)AIF(2nd Round)-Sub-Project CP 2099UGC,MoE,Government of Bangladesh(Contract no.28/2012).
文摘Shock wave-boundary layer interactions(SWBLI)are observed in several practical high-speed internal flows,such as compressor blades,turbine cascades,nozzles and so on.Shock induced oscillations(SIO),aerodynamic instabilities so-called buffet flows,flutter,aeroacoustic noise and vibration are the detrimental consequences of this unsteady shockboundary layer interactions.In the present study,a numerical computation has been performed to investigate the compressible flow characteristics around a 12%thick biconvex circular arc airfoil in a two dimensional channel.Reynolds averaged Navier-Stokes equations with two equation k-ωshear stress transport(SST)turbulence model have been applied for the computational analysis.The flow field characteristics has been studied from pressure ratio(ratio of back pressure,pb to inlet total pressure,p01)of 0.75 to 0.65.The present computational results have been compared and validated with the available experimental data.The results showed that the internal flow field characteristics such as shock wave structure,its behavior(steady or unsteady)and the corresponding boundary layer interaction are varied with pressure ratio.Self-excited shock oscillation was observed at certain flow conditions.Moreover,the mode of unsteady shock oscillation and its frequency are varied significantly with change of pressure ratio.
文摘Numerical works have been conducted to investigate the effect of nozzle geometries on the discharge coefficient.Several contoured converging nozzles with finite radius of curvatures,conically converging nozzles and conical divergent orifices have been employed in this investigation.Each nozzle and orifice has a nominal exit diameter of 12.7x10^(-3)m.A 3rd order MUSCL finite volume method of ANSYS Fluent 13.0 was used to solve the Reynolds-averaged Navier-Stokes equations in simulating turbulent flows through various nozzle inlet geometries.The numerical model was validated through comparison between the numerical results and experimental data.The results obtained show that the nozzle geometry has pronounced effect on the sonic lines and discharge coefficients.The coefficient of discharge was found differ from unity due to the non-uniformity of flow parameters at the nozzle exit and the presence of boundary layer as well.
文摘Turbofan engine intakes are designed to provide separation-free flow at the fan faceover a wide range of operating conditions. But at some off-design conditions, like at high flightspeeds and high angles of attack (AoA), the aero engine intake may encounter flow separation.This boundary layer separation inside the nacelle inlet of an aircraft engine can lead to a largenumber of undesirable outcomes like reduction in fan efficiency, engine stall and high levels ofstress on the fan blades. Active flow control is a promising solution to reduce inlet boundarylayer separation and the associated fan-face flow distortion at such off-design conditions. Byblowing pressurized air into the intake near the separation point, the boundary layer is ener-gized and separation can be controlled. This study investigates the applicability of lip blowing,an active flow control technique, to control intake separation and flow distortion at the fan-face.First, intake separation was triggered in a 3D CFD model based on the NASA CommonResearch Model (CRM) using high AoA cases at cruise condition (Mach number 0.85, Massflow capture ratio w0.7) and the features of separated flow were analyzed. Thereafter, activeflow control was introduce to the intake in the form of two types of lip blowing, direct andpitched blowing. The efficacy of lip blowing at achieving separation control in an ultra highbypass ratio turbofan engine intake has been established through this study. The present paperalso examines the significance of blowing parameters like the type of blowing, blowing pres-sure ratio, and blowing slot dimension, at different angles of attack to identify the critical con-trol parameters. Our research successfully establishes proof of concept by demonstrating the feasibility of using lip blowing for separation control in aero-intakes, via numerical modelling.Furthermore, this study also provides crucial insights regarding the important variables to beconsidered for future experimental studies, and also for detailed studies covering a wider
文摘Conical plug nozzle and truncated conical plug nozzle are advanced rocket nozzles suitable for use as altitude compensating nozzles.In this study flow through the conical plug and truncated conical plug nozzles are numerically simulated to first validate with experimental data and then to compare the performance when a base bleed is introduced.The numerical analysis has considered two-dimensional axisymmetric models.Reynolds-averaged NavierStokes equations are solved with two equation shear stress transport k-ω turbulence model.For the validation of the plug nozzle,flow features and wall pressure along the length of the nozzle is taken for different nozzle pressure ratios.For the validation of truncated plug nozzle,flow features and base pressures at various nozzle pressure ratios are compared.The base bleed is taken as 2%of the inlet mass flow rate.The comparison of results shows that the introduction of base bleed helps to compensate for the loss of thrust due to conical plug nozzle truncation.
文摘Effects of spontaneous condensation of moist air on the shock wave dynamics around butterfly valves in transonic flows are investigated by experimental and numerical simulations.Two symmetric valve disk shapes namely-a flat rectangular plate and a mid-plane cross-section of a prototype butterfly valve have been studied in the present research.Results showed that in case with spontaneous condensation,the root mean square of pressure oscillation(induced by shock dynamics)is reduced significantly with those without condensation for both shapes of the valves.Moreover,local aerodynamic moments were reduced in case with condensation which is considered to be beneficial in torque requirement in case of on/off applications of valves as flow control devices.However,total pressure loss was increased with spontaneous condensation in both the valves.Furthermore,the disk shape of a prototype butterfly valve showed better aerodynamic performances compared to flat rectangular plate profile in respect of total pressure loss and vortex shedding frequency in the wake region.
文摘This article presents an approach which employs a commercial Reynolds-Averaged Navier-Stokes（RANS）solver to predict the steady wake field and loading distributions for a rim driven thruster.Four different cases of propeller blades are chosen to be calculated with the presented method.The propeller blade radial circulation and chordwise circulation density distributions are analyzed.The maximum radial circulation is found at the blade tip,which is different from conventional shaft driven propeller.The numerical results indicate that there is no tip leakage vortex in rim driven propulors.But there exist the tip joint vortex and the root region vortex.Bollard characteristics are calculated by taking rim surface effect into account.From the predicted results the second case in this paper is selected as the final one to perform hydrodynamic experiment.The calculation results with empirical rim surface corrections are compared with the measurement.It shows that the developed numerical method can well predict hydrodynamic performances of the rim driven thruster.
基金supported by Scientific Research Foundation for Returned Scholars,Ministry of Education of China
文摘The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.
基金supported by the National Natural Science Foundation of China (Grant No. 90405007)the Hundred Talents Program of the Chinese Academy of SciencesProgram for Changjiang Scholars and Innovative Research Team in University.
文摘The separated turbulent flow around a circular cylinder is investigated using Large-Eddy Simulation （LES）, Detached-Eddy Simulation （DES, or hybrid RANS/LES methods）, and Unsteady Reynolds-Averaged Navier-Stokes （URANS）. The purpose of this study is to examine some typical simulation approaches for the prediction of complex separated turbulent flow and to clarify the capability of applying these approaches to a typical case of the separated turbulent flow around a circular cylinder. Several turbulence models, i.e. dynamic Sub-grid Scale （SGS） model in LES, the DES-based Spalart-Allmaras （S-A） and κ-ω Shear-Stress- Transport （SST） models in DES, and the S-A and SST models in URANS, are used in the calculations. Some typical results, e.g., the mean pressure and drag coefficients, velocity profiles, Strouhal number, and Reynolds stresses, are obtained and compared with previous computational and experimental data. Based on our extensive calculations, we assess the capability and performance of these simulation approaches coupled with the relevant turbulence models to predict the separated turbulent flow.
基金Project supported by the Key Scientific and Technological Project of Jilin Province(No.20170204066GX)the Green Design Platform Construction Project for High-end Earthwork Machinery,Ministry of Industry and Information Technology,China(No.[2017]327)+2 种基金the Open Foundation of Key Laboratory of Road Construction Technology and Equipment(Chang’an University)Ministry of Education of China(No.310825171104)the Advanced Manufacturing Projects of Government and University Co-construction Program of Jilin Province,China(No.SXGJSF2017-2)
