2015 Vol. 30, No. 6

Display Method:
Characteristics of entropy generation of Blasius flow over permeable wall vertical flat plate
ZHAO Guo-chang, SHAN Long, SONG Li-ping, KONG Jing-ru, ZHAO Heng, LIU Yu, MAO Xiao-dong
2015, 30(6): 1281-1292. doi: 10.13224/j.cnki.jasp.2015.06.001
Abstract:
The entropy generation in the Blasius flow boundary layer over permeable wall vertical flat plate and the ratio of entropy generation due to heat transfer and total entropy generation (Bejan number) were investigated in detail. The coupled momentum and energy partial differential equations (PDEs) were transformed into nonlinear ordinary differential equations (ODEs) using similarity transformation, which were then solved numerically using the Runge-Kutta method and the similarity solutions of the dimensionless flow function and dimensionless temperature were obtained. The effects of dimensionless parameters on the entropy generation and Bejan number were investigated. The results are as follows: for suction, the maximum entropy generation appears on the wall and entropy generation decreases monotonically with the increase of similarity variables, the greater suction speed means the greater entropy generation. For injection, the greater injection speed at the wall means the smaller entropy generation. Near the wall, entropy generation increases with the variable viscosity parameter increases. However this effect away from the wall is reversed. Increases in either the injection/suction parameter or the variable viscosity parameter will increase the Bejan number when it's away from the wall and decrease the Bejan number when near the wall. Increases in either the radiation parameter or Reynolds number will decrease entropy generation and Bejan number, while increases in the Biot number will increase entropy generation and Bejan number. The entropy generation due to heat transfer changes with Biot number dramatically when Biot number is within the range of 0 to 2. The total entropy generation is dominated by the entropy generation due to flow when Biot number is small, and dominated by the entropy generation due to heat transfer when Biot number is large. The entropy generation will increase as the slip parameter increases, while the Bejan number decreases. The effect of Grashof number on entropy generation is relatively large, and entropy generation will have more dramatic changes due to changes of similarity variables when Grashof number is larger. Furthermore, entropy generation due to flow is smaller and Bejan number is bigger with a larger slip parameter.
Experiment on effect of inlet airflow on performance of a reverse-flow combustor
ZHAO Shuo, HUANG Zhang-fang, CHENG Sheng-jun
2015, 30(6): 1293-1297. doi: 10.13224/j.cnki.jasp.2015.06.002
Abstract:
In order to understand the effect of different inlet airflow angles on the performance of a reverse-flow combustor, performance experiments of reverse-flow combustor were carried out. The results show that, different inlet airflow angles have different effects on the outlet temperature distribution of reverse-flow combustor. Compared with the 30 degree vaned axial diffuser and the vaneless axial diffuser in the reverse-flow combustor, the oulet temperature distribution factor in the former is only half of the latter, and the radial temperature distribution factor is also much smaller in the former. Also, different inlet airflow angles have bigger effect on ignition performance.
Thermo-elastic coupling analysis of round-hole flat-plate film-cooling
WANG Zhan, ZHANG Chao, LIU Jian-jun
2015, 30(6): 1298-1306. doi: 10.13224/j.cnki.jasp.2015.06.003
Abstract:
The thermo-elastic coupling property of the basic cooling configuration of flat-plate with round-hole was analyzed. The factors that may influence the steady state thermal stress were studied. The thermal stress distribution of the round-hole flat-plate film-cooling was predicted and analyzed using the multi-field coupling method. The influences of the injection angle, compound angle and blowing ratio on thermal stress were systematically investigated. The results show that the overall thermal stress in the flat-plate is much lower than that of the place around the film holes, and the thermal stress only concentrates near the leading edge line and trailing edge line of the film holes, while the thermal stress remains low on the lateral sides of the film holes. The inner force imposed on the film holes is perpendicular to the axis of film hole. The higher blowing ratio means higher temperature gradient near the film holes, and more intensive thermal stress concentration near the leading edge line and tailing edge line. The smaller injection angle means higher thermal stress near the film holes, and more serious thermal stress concentration on the leading edge point of the film hole exit and the trailing edge point of the film hole inlet. The compound angle is beneficial to weakening the thermal stress concentration.
Heat transfer and flow characteristics of low aspect ratio pin-fins in cooling channel
WANG Long-fei, WANG Song-tao, LU Shao-peng, LUO Lei, WEN Feng-bo
2015, 30(6): 1307-1318. doi: 10.13224/j.cnki.jasp.2015.06.004
Abstract:
Numerical simulation was conducted to investigate the heat transfer and flow characteristics of low aspect ratio cylinder pin-fins arrays in turbine blade trailing edge cooling channel, and analyze the action process of inlet Reynolds numbers and pin-fins spacing on cooling channel heat transfer and flow characteristics. The result shows that, with the increase of inlet Reynolds numbers, the channel wall heat transfer performance is improved effectively while the augmentation ability is gradually weakened, and the array pressure loss coefficient is reduced. For the two kinds of spacing, streamwise spacing is the main factor affecting the wall heat transfer performance, and with the decrease of streamwise spacing, the cooling channel heat transfer performance becomes better while reducing pressure loss coefficient; spanwise spacing is the dominant factor affecting cooling channel flow loss, which has an inverse relation between each other. In calculation of the cooling channel, average heat transfer performance of pin-fins is about 1.8 times of the endwall, and endwall heat transfer weight is approximately 0.824 times of the heat transfer area ratio, while this weight is nothing to do with inlet Reynolds numbers.
Calculation of surface heat transfer coefficient based on Reynolds-averaged Navier-Stokes equations
HOU Shuo, CAO Yi-hua
2015, 30(6): 1319-1327. doi: 10.13224/j.cnki.jasp.2015.06.005
Abstract:
Reynolds-averaged Navier-Stokes (RANS) equations governing two-dimensional external flow were numerically solved using finite volume method. Local surface heat transfer coefficients on surface of smooth and rough NACA0012 airfoil and cylinder were calculated. The effects of near wall mesh spacing, turbulent and surface roughness model on the numerical calculation results were discussed. The features of convective heat transfer in laminar and turbulent boundary layers can be distinguished by shear stress transport(SST)turbulent model, moreover, the transition can be predicted using SST turbulent model. The reasonable heat transfer coefficients on rough walls are calculated using Spalart-Allmaras (S-A) extended model, but the convection heat transfer coefficients in laminar boundary layer can not be calculated effectively by S-A turbulent model neglecting transition function. The near wall mesh spacing should approach order of 10-5 in viscous sub-layer for obtaining accurate heat transfer coefficients on both smooth and rough walls. The curves of surface heat transfer coefficients are close to experimental results if appropriate near wall mesh spacing, turbulent and surface roughness models are used. The effect of compressibility and viscous dissipation is neglected by incompressible RANS equations, but can be incorporated into overall thermal analysis in the form of adiabatic heating term.
Performance of continuously rotating detonation combustor
LIU Qian, ZHENG Hong-tao, LI Zhi-ming, CHEN Xi
2015, 30(6): 1328-1336. doi: 10.13224/j.cnki.jasp.2015.06.006
Abstract:
Using two-dimensional compressible Euler equations, the flow field in continuously rotating detonation combustor (CRDC) with equivalent hydrogen/air was numerically analyzed. The effects of fuel inflow total pressure at different azimuthal sizes of CRDC on pressure ratio, combustion efficiency, volume heat release rate and NOx emissions were examined. The results show that with the increase of fuel inflow total pressure, when azimuthal size of CRDC is 150mm, pressure ratio remains unchanged, and combustion efficiency and volume heat release rate grow linearly; when azimuthal size of CRDC is 250mm, detonation wave number redoubles, and volume heat release rate still grow linearly, but pressure ratio and combustion efficiency turn down due to the increase of fuel inflow total pressure. Compared with conventional isobaric combustor, CRDC has a higher performance. In the calculation region of 250mm×40mm, the pressure ratio goes up to 2.706, greatly improving the potential power of gas; volume heat release rate is 2.18 ×1010W/m3; NOx emissions is 0.692mg/m3, thereby ensuring low emissions requirements.
Influence of fuel stage ratio on pressure oscillation frequency of LESS combustor
QIN Hao, TANG Guan-qiong, LIN Yu-zhen, LI Ji-bao
2015, 30(6): 1337-1343. doi: 10.13224/j.cnki.jasp.2015.06.007
Abstract:
Influence of fuel stage ratio on pressure oscillation frequency of low emissions with stirred swirls (LESS) combustor was tested. The pressure oscillation frequency dropped with the decrease of fuel stage ratio (from 32.7% to 5.21%) under the combustor inlet pressure of 1.14MPa to 2.77MPa, combustor inlet temperature of 645K to 808K, total fuel air ratio of 0.0242 to 0.0303 and combustor pressure drop of 3.14% to 3.64%. A homogeneous straight pipe model and CFD numerical method were adopted to analyze the influence of fuel stage ratio on pressure oscillation frequency. The results show that changes of original position, axial length and mean temperature of flame zone caused by fuel stage ratio are the main reasons contributing to the changes of pressure oscillation frequency.
Micro bubble behavior and droplet oscillations in pendant droplet combustion
LU Qi-zhi, GONG Jing-song, PIAO Ying, YIN Ya-jun, LI Hui, FAN Hong-yu
2015, 30(6): 1344-1349. doi: 10.13224/j.cnki.jasp.2015.06.008
Abstract:
High-speed microscopic camera technology was applied to observe and research RP-3 micro-scale pendant droplet combustion process, the micro bubble nucleation, growth and coalescence and overflow behavior process and droplet oscillations phenomenon. The results show that:①The location of the micro bubbles generated inside droplet.②The volume expansion of the droplet caused by the micro bubble in incipient evaporation.③The droplet oscillates caused by the bubble behavior in combustion.④The droplets stability combust (not burst) phenomenon. The role of surface tension is analyzed in pendant droplets combustion. It is established that the critical condition of stability combust combustion for pendant droplet, the critical characteristic valume within 2.5-5.0μL.
Investigation on structural parameters of airfoil shaped strut flame stabilizer
LIU Guang-hai, LIU Yu-ying
2015, 30(6): 1350-1356. doi: 10.13224/j.cnki.jasp.2015.06.009
Abstract:
Structural parameters of airfoil shaped strut flame stabilizer were designed based on the algorithm of NACA airfoil parameters. The impacts of maximum thickness position and thickness of trailing edge of airfoil shaped stabilizers on the total pressure recovery coefficient and length of recirculation zone were analyzed by numerical simulation. The results show that total pressure recovery coefficient of strut flame stabilizer increases firstly, then decreases sharply with the maximum thickness position of airfoil shaped stabilizer moving towards to trailing edge under the conditions of flow Mach number from 0.15 to 0.30. However, the length of recirculation zone shows a contrary trend. The smaller trailing edge thickness of trailing edge of airfoil shaped stabilizer means the greater maximum total pressure recovery coefficient. The length of recirculation zone is proportional to 1.5 power of thickness of trailing edge.
Reverse design method of hypersonic aircraft inlet compression surface
YUE Xiao-kui, QI Kai-hua
2015, 30(6): 1357-1365. doi: 10.13224/j.cnki.jasp.2015.06.010
Abstract:
For the problem of hypersonic aircraft inlet compression surface, a reverse design method of constant intensity and isentropic compression surface was proposed firstly based on the investigation of the reverse algorithm of the shock wave theory, then inlet compression surfaces of hypersonic flight vehicle were designed using this method with the flow parameters of inlet exit, next performance analysis and numerical simulation were conducted. The results show that: the reverse design method is effective and time-saving as compared with the traditional method. Theoretical calculation shows that both the geometry parameters and the incoming flow parameters of two compression surfaces are most consistent, the total pressure recovery coefficient of the isentropic compression increases by 6.3% to the constant intensity compression at designed point. The numerical simulation has verified the effectiveness of the proposed method and studied the performance parameters of two compression surfaces under off-design condition, while comparing the internal wave system distribution and airflow uniformity of the inlets under different Mach numbers.
Unsteady numerical simulation of propeller/wing interaction
QIAO Yu-hang, MA Dong-li, LI Zhi
2015, 30(6): 1366-1373. doi: 10.13224/j.cnki.jasp.2015.06.011
Abstract:
Unsteady numerical simulation was conducted to study distributed propeller slipstream problem of solar unmanned aerial vehicle based on structured/unstructured hybrid grid CFD methods. Wing was segmented to eliminate the influence of finite span three-dimensional effects and simplify the research object, and the sliding mesh and transition model were applied in the numerical simulation; then the conditions of propeller in front of or behind the wing were discussed and compared with pure wing and pure propeller. The result shows that propeller in front of or behind the wing makes wing lift, drag and nose-down pitching moment increase and fluctuate periodically, while propeller behind the wing has less impact on the wing. The wing causes propeller thrust, absorbed power and efficiency increase and fluctuate periodically, while non-uniformity of flow field results in vibration of peopeller.
Analyze on longitudinal controllability and stability of OSR helicopter
XU Ming, LI Jian-bo, HAN Dong
2015, 30(6): 1374-1381. doi: 10.13224/j.cnki.jasp.2015.06.012
Abstract:
In order to research the effect of rotor speed, blade spring stiffness and horizontal tail area on the controllability and stability of optimum speed rotor (OSR) helicopter, the flight dynamics model was established, and the induced velocity of rotor and aerodynamic interference of helicopter was calculated by free wake model. Finally the differential method was used to calculate the aerodynamic derivatives of the helicopter and control the derivatives matrix. The result shows that the controllability of the helicopter decreases with the decreasing of the rotor speed, but the change of stability about angle of attack was opposite. However the increase of blade spring stiffness in root can not only improve helicopter controllability, but also reduce the stability. The longitudinal stability of the OSR helicopter was decreasing without horizontal tail, and the handling qualities of sample helicopter are put at the third level in ADS-33E, but the period of instability oscillation cycle was longer, and the controllability of the helicopter was enhanced after cancelling horizontal tail. Thus, the scheme without horizontal tail for OSR helicopter is still controlled completely through the flight control system.
Fluid-structure interaction study of slide tail plate of single expansion ramp nozzle
GUO Shuai, XU Jing-lei, GU Rui
2015, 30(6): 1382-1390. doi: 10.13224/j.cnki.jasp.2015.06.013
Abstract:
The present research focused on fluid-structure interaction effects on an adjustable single expansion ramp nozzle (SERN). Firstly, the accuracy of the computational fluid dynamics/computational structure dynamics (CFD/CSD) serially loosely coupled algorithm was validated. Then, the numerical procedure was adopted to investigate the interaction between the slide tail plate mounted at the end of SERN upper ramp and the surrounding flow-field. The fluid-structure interaction responses with different slide tail plate lengths (150mm, 180mm) and thicknesses (5mm, 7mm) were researched. The numerical results indicate that, the vibration of the slide tail plate, which is caused by the fluid-structure interaction effects, leads to the aerodynamic parameters fluctuation in nozzle.For the slide tail plate length of 180mm, the nozzle lift is reduced by 1.38%, and the pitch moment difference is increased by 16.9%. Under the same slide tail plate thickness (5mm), shorter slide tail plate possesses lower aerodynamic performance (with a decrease of 9.9% in lift) and better dynamic characteristic (with a decrease of 44.7% in response time). Under the same slide tail plate length (180mm), better aerodynamic and dynamic performances can be achieved when the slide tail plate thickness is increased, namely, the lift is increased by 0.98% when steady state is reached, the pitch moment difference is decreased by 10.5%, and the response time is reduced by 40.4%.
Empirical mixing model of transverse gaseous jet in supersonic flow
YANG Yi-xin, WANG Zhen-guo, SUN Ming-bo, WANG Hong-bo
2015, 30(6): 1391-1399. doi: 10.13224/j.cnki.jasp.2015.06.014
Abstract:
In order to obtain empirical correlations of jet penetration for engineering application, a empirical mixing model for the analysis of transverse gaseous jet in supersonic flow was investigated by modeling the jet centerline, penetration boundary and injection's spanwise distributional characteristics. Several experimental cases were used to verify the performance of the present model. Results indicate that, within the momentum flux ratio of 0.5-5.3, the present model agreed well with the experiment image both in the jet centerline and penetration boundary, and in the area near jet centerline, the model predicted a injection mole fraction with an error less than 5%, which shows a higher precision in contrast with the existing models. With the increasing of the momentum flux ratio, the counter-rotating vortex pairs are strengthened in the near region downstream of the injector, which leads to some inaccuracy of spanwise mole fraction distribution calculated by the present model. This indicates that the present model validity is affected by the momentum flux ratio and downstream distance.
Effect of entry shape on starting characteristics of inward turning inlets
WANG De-peng, TIAN Fang-chao, ZHANG Qi-fan, TAN Hui-jun, BU Huan-xian, LI Guang-sheng
2015, 30(6): 1400-1406. doi: 10.13224/j.cnki.jasp.2015.06.015
Abstract:
In order to investigate the low Mach number starting characteristics of hypersonic inward turning inlets, four kinds of inlets with rectangular, square, elliptical and round entry shapes were designed. Simulation method was used to investigate their unstarting/restarting flow-field structures. Results show that, the inlets with rectangular and elliptical entry shape have similar restarting Mach number while the inlet with square and round entry shape have the similar and higher restarting Mach number. Further analysis shows that the difference of restarting capability may be caused by the different aspect ratios of the entry shape, therefore, a series of inward turning inlets with different aspect ratios were designed, and their unstarting/restarting flow-field structures were investigated. Results show that the aspect ratio of the entry plays an important role in inlet starting characteristics of inward turning inlet. The larger aspect ratio means the easier inlet restarts.
Analysis and design of stratospheric airship propulsion system
LIU Xin-qiang, HE Wei-liang
2015, 30(6): 1407-1413. doi: 10.13224/j.cnki.jasp.2015.06.016
Abstract:
For the problem of propeller and motor matching, brushless direct current (DC) motor model and propeller model were constructed. In order to research the control response problems of propulsion system at 0 to 10m/s cruising speed of airship and the simulation model of control system of propeller and motor was designed in Simulink by combining matching design and motor speed control. The interface of propeller and motor matching was compiled in graphical user interface (GUI). Propeller and motor matching software was developed. The matching design curves for operation at 0 to 10m/s cruising speed of airship, the working range of propeller from 535 to 1071r/min and working efficiency of the propulsion system from 0.558 to 0.593, were obtained, which can be used as the reference basis for matching design between high altitude propeller and motor.
Geometric modeling and 3-D numerical simulation for gooseneck of compressor
TU Qiu-ye, CHEN Jie, JIANG Ping, YAN Hong-ming, CAI Yuan-hu
2015, 30(6): 1414-1422. doi: 10.13224/j.cnki.jasp.2015.06.017
Abstract:
The employed S-shaped polynomial curves as the basis of the gooseneck of compressor geometric modeling. The method used the relative location of S-shaped inflection point of inner wall curve, the peak point and its value of the area distribution rate to control the gooseneck geometrics modeling. A series models were constructed from the baseline gooseneck modeling to analyze with 3-D numerical simulation. The results indicate that the peak value of area distribution rate has the most important effects on the total pressure loss of gooseneck. An relative Mach number can be obtained by changing the peak value of area distribution in order to reduce the flow separation near the outer wall and after the struts.Meanwhile inner walls with large curvature near the entrance of gooseneck modeling, and area distributions rate with an appropriate position of peak point can improve the streamline of outer wall. For example, the total pressure loss of gooseneck reached its minimum at relative location of inflection point of inner wall of 0.18, and relative location of peak point of area distribution rate of 0.20 at relative Mach number of 0.65.
Discussion on several problems of axial flow fan/compressor aerodynamic performance experiment
XIANG Hong-hui, HOU Min-jie, GE Ning, LIU Zhi-gang, WEN Zhen-rong
2015, 30(6): 1423-1431. doi: 10.13224/j.cnki.jasp.2015.06.018
Abstract:
In order to improve the rationality of experimental process control and the authenticity of experimental results assessment of aero-engine fan/compressor, several important technology problems were summarized by analyzing the aerodynamic performance experimental research situation of the domestic aero-engine compression components. Through theoretical analysis and partial experimental validation, these technology problems were discussed in depth. The analysis results show that the requirements of experimental facility exhaust system ability are different for fan aerodynamic performance experiment and compressor aerodynamic performance experiment; compared with the method of rising rotating speed collection, decreasing rotating speed collection can optimize the fan/compressor experimental process; the inherent effects of measuring probes on fan/compressor experimental results can be reduced through different methods; the stagger angle positioning and full stroke adjustment accuracy of the adjustable guide vane/stator vane are the main influential factors to the repeatability of aerodynamic performance experimental results of variable geometry multistage compressors.
Aerodynamic calculation of wind turbine wheel based on hybrid panel viscous-vortex particle method
HU Hao, SONG Xiao-yong, GU Bo, SONG Xian-cheng
2015, 30(6): 1432-1439. doi: 10.13224/j.cnki.jasp.2015.06.019
Abstract:
To accurately and rapidly predict the aerodynamic performance of large wind turbine wheel, a hybrid panel viscous-vortex particle (HPVP) method was established to calculate the wind turbine aerodynamic performance, and corresponding calculational program was developed, the model experiments in controlled conditions (MEXICO) wheel was chosen to validate the accuracy and reliability of HPVP method by comparison of calculational result, test data and CFD method. The results show that the pressure distribution of the main working zone of blade is predicted accurately in HPVP method. Compared with CFD method, when the flow separation is small, the aerodynamic calculation in HPVP method is faster than that of CFD method, and the calculational time by HPVP method is one thousandth of CFD method. Additionally, some other flow details of flow field of the wind turbine wheel can be also showed by HPVP method including the pressure distribution of blade.
Experiment on the instability process in a one-stage transonic axial compressor
PAN Tian-yu, HE Lei, LI Zhi-ping, LI Qiu-shi
2015, 30(6): 1440-1447. doi: 10.13224/j.cnki.jasp.2015.06.020
Abstract:
To fully understand the instability process and reveal the occurrence, development and essence of the unstable flow phenomena of one-stage transonic axial compressor with high loading, test investigation was carried out with both steady and dynamic measurements. The analysis was performed by the original signal, the low pass filtering and the FFT (fast Fourier transform) method. The results show that the large amplitude, low frequency axisymmetric axial disturbance appears firstly at the stator hub during the instability process. At that time, there is a decrease in the performance below 50% span of the compressor. But the compressor can still be operated. Meanwhile, considering the features of this disturbance such as low frequency, axisymmetricity, large amplitude, this kind of unstable flow phenomena, named as partial surge, is identified. Then, this disturbance spreads axially and radially with the further decrease of the mass flow. Next, this disturbance extend to the full span. Finally, it leads to the rotating stall cells at the rotor tip region, and the indicating the compressor reached the instability boundary.
Numerical simulation of leakage characteristic of assembled endwall slot of turbine guide vane
LIU Ya-fei, MAO Jun-kui, XU Rui, ZHAO Zun-sheng, HE Yi-hong, LI Wei
2015, 30(6): 1448-1459. doi: 10.13224/j.cnki.jasp.2015.06.021
Abstract:
Numerical simulations were applied to study the leakage flow of assembled endwall slot of turbine guide vane of an aero-engine. The effects of leakage flow on the structure of flow field of cascade passage and performance parameters of the cascade were also discussed. The leakage flow pressures, the slot width and position of the slot were changed to investigate the rules of leakage flow rate and energy loss coefficient varying with these parameters and discover the effect rules. It is found that the leakage can appear due to different pressures between cooling air and gas, spiral vortex system is formed in the middle position of the slot at the same time, and the leakage flow affects the secondary flow around the endwall significantly. The effect decreases along the direction of blade height and the largest region is up to 44.44% of the blade height. With the increase of leakage flow pressure, the slot width and the angle between the slot and major axis of engine, the energy loss coefficient of the cascade and the leakage flow rate increase monotonically. Numerical results show that the energy loss coefficient of the cascade increases by 14%-62% because of the leakage flow through the slot of the endwall in the parameter range of present study.
Numerical investigation of parameters transmission on blade rows interface of contra-rotating compressor
DENG Xiao-ming, GAO Li-min, GAO Lei, CAI Yu-tong
2015, 30(6): 1460-1469. doi: 10.13224/j.cnki.jasp.2015.06.022
Abstract:
Five steady blade rows interface methods were adopted to investigate an extraordinary turbomachinery—contra-rotating compressor. The whole performance, near design condition and near stall condition flow field of contra-rotating compressor were analyzed in combination with experimental data, so as to make a quantitative comparison of the effect on blade row interface parameters transmission among different interface methods. Results validate: (1) the whole compressor efficiency shows obvious difference among different interface methods, while the pressure ratio agrees with experimental data very well. (2) Under the near design condition, the parameters transmission on different interfaces are uniform for conservative coupling by pitchwise rows method. Local conservative coupling method and matching mixing plane methods behave better along downstream interface, while non reflecting 1-D method and full non matching frozen rotor method perform just the opposite. (3) Under the near stall condition, the high entropy loss area surrounding blade leading edge of non reflecting 1-D method and full non matching frozen rotor method expands much bigger. Mainstream and overflow interface within the blade tip twists into an arc closer to the upstream blade row interface.
CFD simulation of cryogenic propellant tank pressure variation
CHEN Liang, LIANG Guo-zhu, WEI Yi, HU Wei
2015, 30(6): 1470-1477. doi: 10.13224/j.cnki.jasp.2015.06.023
Abstract:
In order to predict the pressure variation and investigate the interrelation among the physical processes in cryogenic propellant tanks during ground parking, a 2-D axial symmetry volume of fluid (VOF) computational fluid dynamic (CFD) model was established, including a liquid propellant phase and a mixture gas phase. The propellant phase change model was cited based on the assumption of thermodynamic equilibrium. The difference between the simulated pressure rise rate and the test data of a test liquid hydrogen tank was 9.1%. The simulation results of the pressurized stage of liquid hydrogen tank and liquid oxygen tank during ground parking indicate that the heating on the ullage by the heat leakage from the inner surface is the key factor of the pressure rise in the liquid hydrogen tank, and the liquid hydrogen vaporization has comparatively less influence. However, the pressure in the liquid oxygen tank varies as the liquid oxygen phase transition changes at the beginning of the propellant pressurized parking stage. The pressure variations in these two tanks are both under the influence of the propellant phase change affected by the different change trends of the propellant convection induced by the heat transfer between the gas and the liquid.
Design of slosh baffle in large cryogenic tank under harmonic excitation
JIN Yong-hua, LI Yan-zhong, WANG Lei, LIU Zhan, LI Cui
2015, 30(6): 1478-1485. doi: 10.13224/j.cnki.jasp.2015.06.024
Abstract:
A fixed construction slosh baffle arranged in the whole tank was designed for a large cryogenic tank based on the horizontal ring slosh baffle and vertical slosh baffle. The effects of this slosh baffle at different fill levels were simulated numerically, using CFD software Flow-3D under a harmonic excitation. Comparisons with a horizontal ring slosh baffle at the liquid surface show that the effects of these two baffles are similar, both of them can significantly reduce the wave range of the center of liquid mass at different fill levels. The wave range can even be decreased by 64% at 50% fill level. During the slosh, the pressure imposed on the single board of the fixed construction slosh baffle is smaller than that of horizontal ring slosh baffle at the liquid surface, while the tank wall of the former undergoes higher pressure. The effect of horizontal ring slosh baffle at the liquid surface is more obvious than the vertical slosh baffle when they have the same board areas. The wave range of the center of liquid mass in the tank decreases as the fill level increases.
Nonlinear fatigue damage of nitrate ester plasticized polyether propellant for strain-control
GAO Yan-bin, XU Jin-sheng, CHEN Xiong, HU Shao-qing
2015, 30(6): 1486-1491. doi: 10.13224/j.cnki.jasp.2015.06.025
Abstract:
In order to explore the nonlinear fatigue damage characteristics of nitrate ester plasticized polyether (NEPE) propellant, a nonlinear fatigue damage model was proposed based on damage mechanics and visco-elasticity theory, which took into account the stiffness degradation in the fatigue process. With dynamic momechanical analysis (DMA) apparatus, NEPE propellant fatigue tests were conducted under five kinds of strain levels, the model parameters were obtained and the presented model was verified. A morphology analyze was conducted with charge coupled device (CCD), the effect of strain loading history on the mechanical properties of the NEPE propellant was analyzed. The result shows that the presented model is capable of describing micro-mechanically the nonlinear fatigue damage evolution of NEPE propellant properties; the repeated tensile leads to dewetting between the particles and matrix propellant. At the initial stage of the cycle phases, the velocity of dewetting is fast. With the increased number of cycles, the extent of dewetting tend to become stable.
Loading method on fatigue test of propeller hub
YANG Lei, SU Jun, WANG Hong-bin
2015, 30(6): 1492-1497. doi: 10.13224/j.cnki.jasp.2015.06.026
Abstract:
To resolve the problem of interference between aerodynamic load and centrifugal load in fatigue test of propeller hub, a loading method was designed based on the analysis of load on the propeller hub and its testability. Using the test principle of full bridge circuit flexibly, the loading method was formed on the basis of strain output control with the processes of static and dynamic calibration. Besides, the correlation between the calibration load and non-calibration load during the processes of load calibration and the loading error resulted from the correlation were analyzed. The interference between these two loads was reduced by means of strain linear superposition. Thus aerodynamic load and centrifugal load are applied accurately with the error of the loading system less than 5%.
Magnetic-elastic coupled model of giant magnetostrictive actuator with flexure-hinges amplifier
YANG Zhao-shu, HE Zhong-bo, BAI Hong-bai, LI Dong-wei, CUI Xu
2015, 30(6): 1498-1506. doi: 10.13224/j.cnki.jasp.2015.06.027
Abstract:
To extend limited stroke of giant magnetostrictive actuator (GMA), a mechanical amplifier with flexure-hinges was embodied, the actual amplification stroke ratio and dynamic behavior of the mechanical amplifier with flexure-hingers were studied using theoretical material mechanics, Lagrange's dynamical equation and finite element method. The magnetic-elastic coupled model of GMA with flexure-hinges amplifier was established based on the constituted model of Jiles-Atherton(J-A) model and quadratic model; the coupled model was simulated in the co-simulation environment of software AMESim&Simulink. A prototype of GMA with flexure-hinges amplifer was manufactured, followed by a series of dynamic performance tests. Both simulation and test result shows that the stroke of GMA with flexure-hinges amplifier reaches 0.4mm; the frequency bandwidth is above 100Hz and the effectiveness of the coupled model is validated.
A method of structural dynamic characteristics analysis for multistage bladed disks
CHEN Yu-di, WANG Yan-rong, XU Zhu-shan
2015, 30(6): 1507-1511. doi: 10.13224/j.cnki.jasp.2015.06.028
Abstract:
Concerning different blade numbers in each stage, the commercial finite element solvers developed on the traditional cyclic symmetry method can not be applied directly to the structural dynamic analysis of multistage bladed disks. A method was developed to deal with this issue by applying corresponding cyclic symmetry boundary constraints on each stage's sector model, and coupling constraints were also applied on the inter-stage's nodes to ensure continuous displacement field.To achieve a better accuracy, appropriate number of sectors for each stage was selected to consider the coupling effects to the greatest degree. The method could be easily implemented in the commercial finite element solvers since the cyclic symmetry boundary conditions for each stage and the inter-stage coupling constraints could be handled by the displacement degrees of freedom respectively. The modal analysis got by this method shows a good agreement with the full circle method and the frequency difference is below 3%. The influence of coupling effect on the multistage bladed disks structure was presented and with the increase of nodal number the coupling effects decrease.
Stiffness design, simulation and test of laminated spherical elastomeric bearing
CHEN Gao-sheng, YANG Yan
2015, 30(6): 1512-1519. doi: 10.13224/j.cnki.jasp.2015.06.029
Abstract:
The computational formulations of compression stiffness and torsion stiffness of every rubber layer of laminated spherical elastomeric bearing were deduced in theory, and the theoretical results and simulated results of compression stiffness and torsion stiffness of every rubber layer of laminated spherical elastomeric bearing were compared through examples. The results show that the theoretical result of compression stiffness is basically in accord with the simulated result, the maximum error is 3.98%. There are some differences from simulated result and theoretical result for torsion stiffness in every rubber layer; the simulated result of torsion stiffness is bigger than theoretical result for the first half of all the rubber layers from the little end, the maximum error is 33.3%; the simulated result of torsion stiffness is smaller than theoretical result for the second half of all the rubber layers, the maximum error is 32.8%. The nonlinear characteristic of compression stiffness is not evident for the alterable pressure, and that of torsion stiffness is more evident with the increase of the torsion angle. But the equivalent compression stiffness and equivalent torsion stiffness got from theoretical method and finite element simulation method are in good agreement with the test results.
Modeling and analysis of critical speed for tail drive shaft of helicopter considering space maneuvering flight
NI De, ZHU Ru-peng, LU Feng-xia, BAO He-yun, JIN Guang-hu, LI Miao-miao
2015, 30(6): 1520-1528. doi: 10.13224/j.cnki.jasp.2015.06.030
Abstract:
The universal dynamical model of tail drive shaft of helicopter considering maneuvering flight was established. The reasons for maneuvering flight affecting the critical speed of tail drive shaft of helicopter were analyzed. It is pointed out that only a portion of the rotational motion components can affect the critical speed of tail drive shaft, including: pitching angular velocity, yawing angular velocity, rolling angular velocity and angular acceleration. The influence law of the pitching angular velocity, yawing angular velocity and rolling angular velocity on the first three order critical speeds of drive shaft was discussed through a numerical example. The results reveal that the pitching angular velocity reduces the critical speeds of tail drive shaft, but only has a greater impact on the first order critical speed of tail drive shaft in a backward whirl, and the yawing angular velocity has the same influence. The positive rolling angular velocity increases the critical speed of tail drive shaft in a forward whirl, and reduces the critical speed of tail drive shaft in a backward whirl; but in the case of negative roll maneuver, there is a reverse case compared with the positive roll maneuver that the backward rolling angular velocity reduces the critical speeds of tail drive shaft in a forward whirl and increases the critical speeds of shaft in a backward whirl, before the boundary of the rolling angular velocity equals to a special value. The rolling angular velocities corresponding to the boundary of variation curve of the first, second and third critical speeds are respectively 1.86, 27.41 and 124.5rad/s.
Fault diagnosis for actuator and its sensor of aero-engine fuel system
JIANG Jie, LI Qiu-hong, ZHANG Gao-qian, JU Xin-xing
2015, 30(6): 1529-1536. doi: 10.13224/j.cnki.jasp.2015.06.031
Abstract:
A single fault diagnosis and location method of aero-engine actuator and its sensor was proposed. The method was based on actuator model and aero-engine inverse model. A three-order transfer function model of actuator was established according to its closed-loop construction. The aero-engine inverse model was built based on two parallel BP (back propagation) neural network including a steady inverse module and a dynamic compensation module, so the fuel mass flow estimation can be achieved based on the outputs of aero-engine. Fault diagnosis according to the deviation between the actuator model output and the sensor output was adopted. Fault location according to the deviation between the aero-engine inverse model output and the sensor output was used. The simulation results on an aero-engine and its fuel system actuator models show that the fault diagnosis system can diagnose amplitude in more than 1.6% of actuator and its sensor fault and locate it accurately at the steady and dynamic states of aero-engine, which verifies the effectiveness of the proposed method.