2018 Vol. 33, No. 8

Display Method:
Two-phase heat transfer characteristics of R141bin rectangular micro-channels
2018, 33(8): 1793-1800. doi: 10.13224/j.cnki.jasp.2018.08.001
Abstract:
An experimental study of two-phase flow heat transfer characteristics in rectangular micro-channels of 1mm and 0.5mm hydraulic diameter was conducted. R141b was used as the working fluid. In this experiment, the heat flux ranged from 1 to 16kW/m2, the mass flow rate ranges from 111.1kg/(m2·s) to 333.3kg/(m2·s), and the mass vapor quality ranged from 0 to 1. In order to know the main factor affecting heat transfer, the relationships of average heat transfer coefficient between heat flux, mass flow rate and mass vapor quality were analyzed. Result showed that, when the heat flux was low, the average heat transfer coefficient decreased as heat flux increase was affected by boiling heat transfer, and the opposite trend came when the mass flow rate became large. Under the condition of low heat flux, the average heat transfer coefficient changed obviously with the mass flow rate, and the mass flow rate had weak effect on the average heat transfer coefficient after the heat flux increased to a certain value. When the mass flow rate ranged from 111.1kg/(m2·s) to 333.3kg/(m2·s), the average heat transfer coefficient decreased as mass vapor quality increased.
Effect of thermal radiation on the heat transfer performance of a film cooling vane
2018, 33(8): 1801-1810. doi: 10.13224/j.cnki.jasp.2018.08.002
Abstract:
For a film-cooled first stage turbine vane, gas-solid-thermal coupled simulation was used to analyze the effect of radiation factors on the temperature distribution and cooling efficiency by comparing the temperature ratio and cooling efficiency with consideration/without consideration of the radiation heat flux over the vane. It showed that when the inlet blackbody temperature changed from 1200K to 1900K and the emissivity of vane surface changed from 0.3 to 0.7, the vane surface temperature increased obviously while considering the effect of radiation factors. When the inlet blackbody temperature was 1600K and the surface emissivity of the blade was 0.5, the temperature rise of the pressure surface was about 100K, and the maximum temperature(1350K) of the blade surface increased about 50K. The gas radiation caused about 5% temperature rise on the surface of suction side and trailing edge. The cooling efficiency decreased while considering the thermal radiation, the vane leading edge and pressure surface were not able to meet the cooling requirements even though the gas film holes were close together, the overall cooling efficiency reduced to 0.3 or less.
Structural strength simulation of film cooling vane after heat shock by thermal/flow/structure coupling
2018, 33(8): 1811-1820. doi: 10.13224/j.cnki.jasp.2018.08.003
Abstract:
Study on influence of coolant mass flux on the thermal stress around the hole provided a reference for the design of the reliability of the film cooling vane. Combining the finite element method and boundary element method, the pore size was changed, and the finite element model was established to obtain the maximum temperature, temperature imbalance degree and maximum thermal stress of the vane after thermal shock with transient thermal/flow/structure coupling technology. Research showed that increasing the coolant mass flux was beneficial to improve the cooling efficiency and reduce the vane temperature, but also increase the imbalance temperature of the vane, which may lead to more serious thermal stress concentration at the trailing edge. Increasing the pore diameter at the leading edge can increase the average cooling efficiency by 66%, helping to slow down the thermal stress in the pores. However, increasing the trailing air hole diameter had moderate effect on average cooling efficiency and thermal stress. Above all, numerical results were consistent with experiment and theory, showing the results were valuable reference for aero-engine turbine vane cooling designers.
Analysis of heat transfer characteristic of cylindrical foam porous block with high speed flow around
2018, 33(8): 1821-1829. doi: 10.13224/j.cnki.jasp.2018.08.004
Abstract:
Numerical investigation on heat transfer characteristic inside the cylindrical foam porous block attached to frontal surface of cylinder body with high speed flow around was performed by employing single-domain approach through integration of fluid and porous region. Radiation effect within porous region was considered based on Monte Carlo method. The influences of variable porous region length and porous drag characteristic on both the wave drag of entire model and the aero-heating effect of frontal porous region were analyzed. The results show that installation of certain length porous foam material with appropriate drag characteristic in front of the cylinder body can reduce both the wave drag of entire body and aero-heating effect of the frontal surface. Under simulation condition, the frontal foam porous with parameters (dimensionless length 1.0, viscous coefficient 0.2×107m-2 and inertial coefficient 200m-1) can reduce 13.5% of the wave drag and approximately 75% of the frontal surface average aero-heating heat flux. With the dimensionless length being kept unchanged, the wave drag was decreased with the decline of inertial coefficient of foam porous region, whereas the aero-heating flux of frontal surface was slightly increased.
Thermal cracking and heat transfer of hydrocarbon fuel with catalytic steam reforming
2018, 33(8): 1830-1837. doi: 10.13224/j.cnki.jasp.2018.08.005
Abstract:
For the thermal cracking and steam reforming reactions of hydrocarbon fuel in regenerative cooling channel, a numerical model coupling with flow, heat transfer and chemical reaction was established to investigate the heat absorbing and reaction properties of hydrocarbon fuel, taking thermo-physical properties under supercritical pressure into consideration. Results showed that, the calculated results were in good agreement with experimental data. The numerical model can predict the wall and fuel temperature distribution, fuel conversion and heat transfer deterioration phenomenon. Fuel with catalytic steam reforming had the advantages of heat sink promotion and lower outlet temperature. Thermal cracking reaction can be inhibited by the existence of steam reforming reactions. The analysis of fuel mass flow rate showed that residence time of the fuel decreased with the increasing mass flow rate, resulting in a lower conversion and chemical heat sink of the hydrocarbon fuel.
Effect of cavity on combustion characteristics of integrated strut flame stabilizer
2018, 33(8): 1838-1844. doi: 10.13224/j.cnki.jasp.2018.08.006
Abstract:
Effect of cavity on combustion efficiency and blow off characteristic of fuel injection/stabilization integrated strut flame stabilizer were experimentally investigated under the condition of incoming flow temperature 780-850℃, incoming flow Mach number 0.16, fuel air ratio 0.002-0.006. Numerical simulation was also used to assist analysis. The results show that integrated strut flame stabilizer with cavity could achieve stable and high efficiency flames under the condition of various fuel air ratios; combustion efficiency of integrated strut flame stabilizer without cavity is worse than that of the integrated strut flame stabilizer with cavity, the differences between the integrated strut flame stabilizer with cavity and without cavity decreases gradually with the increase of fuel air ratio; the flame stabilizer with cavity performs better on blow off than flame stabilizer without cavity. Cavity structure could promote the atomization and evaporation characteristic of fuel, and improve combustion characteristics of integrated strut flame stabilizer.
Unsteady cavitating flow of liquid hydrogenaround the ogive body
2018, 33(8): 1845-1854. doi: 10.13224/j.cnki.jasp.2018.08.007
Abstract:
In order to investigate the unsteady cavitating flow characteristics of the cryogenic fluids in the aircraft engine, Zwart cavitation model and large eddy simulation (LES) turbulence model were used to simulate the unsteady cavitating flow of liquid hydrogen around an ogive body. The results showed that compared with the experiments, the unsteady cavitating flow of liquid hydrogen around the ogive body could be accurately simulated by the numerical calculation model. The cavitation evolution process was divided into three stages: the growth of attached cavity, the development of large-scale cavity and the development of small-scale cavity. The analysis of flow field showed that the re-retrant jet was the main cause of the cavitation shedding and instability. Besides, the interaction between cavitation and vorticity was also analyzed. The vortex streching, vortex dilation and baroclinic torque terms were contributed to cavitation structures of leading and trailing parts of cavity, inside the cavity and at the interface between the cavity with mainstream and head of the re-retrant flow, respectively.
Numerical study on transport aircraft after-body flow separation control by spark jet
2018, 33(8): 1855-1863. doi: 10.13224/j.cnki.jasp.2018.08.008
Abstract:
Flow active control on transport aircraft after-body flow separation by the spark jet was numerically studied. The results show that the spark jet directly supplied the high-energy momentum to the boundary layer of the separation zone. Vortices induced by spark jet strengthened the mixing between the high-speed main flow and the low-speed boundary flow. It delayed the separation of the after-body flow and reduced the size of the separation zone. Under the same spark jet parameters, the effect of active control was better when the spark jet was led into the boundary layer at the position of the separation point or the downstream separation point. The average drag coefficient was reduced to 3.26%. Compared with the mass equivalent steady jet, pulsed jet and the pulsed jet had better performance on drag reduction than the steady jet. The spark jet was superior to the pulsed jet on drag reduction.
Numerical study on matching mechanism and efficiency ofpulse detonation combustor and impulse turbine
2018, 33(8): 1864-1871. doi: 10.13224/j.cnki.jasp.2018.08.009
Abstract:
As the turbine inlet was dominated by strong unsteady airflow when working with the pulse detonation combustor (PDC), the type of turbine was determined based on the loss mechanism and the aerodynamic design work of turbine components was carried out.The numerical method was used to calculate the efficiency of the designed turbine and explore the mechanism of the interaction between the PDC and the turbine.The results showed that partial intake, impulse turbine with less reaction degree was more suitable for PDC; a strong reflection shock wave was formed on the contraction section of nozzle,the leading edge and the pressure surface of the rotor blade, which cause energy loss; the efficiency of turbine was about 75% at the design point. The results can provide some reference for turbine design of pulse detonation turbine engine.
Effect of temperature ratio on film cooling characteristics at first-stage vane endwall
2018, 33(8): 1872-1879. doi: 10.13224/j.cnki.jasp.2018.08.010
Abstract:
Experimental investigations of the film cooling characteristics at a vane end-wall with ten rows of compound film-holes were carried out in a hot wind tunnel using infrared thermograph technique. Two kinds of mainstream-to-coolant temperature ratio(TR) of 1.64 and 2.68, and four kinds of coolant-to-mainstream blowing ratio(BR) of 0.6, 1.0, 1.5 and 2.0 were chosen. Comparing the results at high and low TR, it revealed that: (1) in the case of no coolant injection, with the increase of mainstream temperature, the point of maximum temperature at the end-wall leading edge was pushed to vane leading edge; (2) in the case of coolant injections, the overall cooling effectiveness near suction side was higher than pressure side. The area-averaged cooling effectiveness increased with BR, moreover, the largest increment was generated from BR of 1.0 to 1.5; (3) the cooling effectiveness increased with TR at the same BR, but the extent of variation was reduced with the increasing BR, i.e. comparing with TR of 1.64, the area-averaged overall effectiveness can be enlarged by 18.2% at BR of 0.5, but by 8.8% at BR of 2.0.
Physical properties of coke from the RP-3 kerosene
2018, 33(8): 1880-1885. doi: 10.13224/j.cnki.jasp.2018.08.011
Abstract:
Some morphology and physical properties of the coke from the RP-3 kerosene at ambient temperature 800K and 900K were experimentally studied. The empirical correlations of specific heat capacity and thermal conductivity were summed up. The results indicate that the coke is an accumulation of microcosmical spherical particles. There were pores and cracks on the surface of coke, which had an envelope density of approximately 1049kg/m3, a real density of approximately 1498kg/m3 and a porosity of approximately 29.9%. The coke had lower real density formed at ambient temperature 800K compared with that at ambient temperature 900K. Specific heat capacity and thermal conductivity of the coke were measured by the flash method. Specific heat capacity was approximately 1.1 to 2.2J/(g·K) and thermal conductivity was approximately 0.19 to 0.28W/(m·K), both of which increased with increasing temperature. The coke had higher specific heat capacity and thermal conductivity formed at ambient temperature 800K compared with that at 900K.
Stress intensity factor calculation and fracture evaluation based on weight function method
2018, 33(8): 1886-1895. doi: 10.13224/j.cnki.jasp.2018.08.012
Abstract:
An approximate and simple expression was presented to calculate the crack face displacements of collinear cracks and center cracks by use of only one reference stress intensity factor. The crack face displacement and its partial derivative determined by the present method were in good agreement with the exact solutions. Calculations of the weight function for the center crack were reduced to the simple quadrature of the correction function and the partial derivative of the crack face displacement. Based on the weight function method, the relationships between the stress intensity factor of the rotating blade with a center crack and the crack length, crack location, rotating speed and angular acceleration were studied. It showed that the stress intensity factor increased with the increasing crack length and the decreasing distance from the crack line to the rotating axis; the influence of the angular acceleration was negligible, but the increasing rotating speed made the stress intensity factors increase. The critical rotating speed of the blade was evaluated based on the fracture law, and decreased with the increasing crack length.
Fatigue life of turbine disk based on actively managed thermal loading method
2018, 33(8): 1896-1904. doi: 10.13224/j.cnki.jasp.2018.08.013
Abstract:
In order to reduce the failure risk of turbine disk and extend its life, actively managed thermal loading method was applied to the turbine disk with pre-set crack. Meanwhile, the crack stress intensity factors were calculated by general weight function method and the life of the turbine disk was analyzed. The influence of the energy distribution based on actively managed thermal loading method on the life of the turbine disk was studied. The relationship and change trend between the thermal boundary loading and the life of the turbine disk were explored. The reason and mechanism were analyzed by finite element numerical simulation. The consequence shows that actively managed thermal loading method can effectively optimize the temperature distribution of the turbine disk, reduce the stress near the crack, delay the expansion of the crack, and significantly improve its life and safety. When the thermal boundary loading coefficient was 0.05 and 0.10 respectively, the corresponding lives increase by 12.2% and 26.1%.
Aerothermoelastic flutter mechanisms of functionally graded curved panel in supersonic flows
2018, 33(8): 1905-1915. doi: 10.13224/j.cnki.jasp.2018.08.014
Abstract:
Aerothermoelastic flutter mechanisms and bifurcation characteristics of functionally graded curved panel (FGCP) in supersonic flows were investigated. Piston theory and Eckert reference enthalpy methods were used to simulate aerodynamic force and heating, respectively. 2-D heat conduction equation was solved and the impact of elevated temperature on physical properties was considered to build an aerothermal-aeroelastic two-way coupling model of FGCP. Using finite element method, bifurcation diagrams under different arch rise were presented, and the impact of arch rise on the bifurcation diagram was discussed in detail. Three kinds of behaviors such as: thermal buckling, chaos and regular vibration were discovered. Comparative studies of two different regular vibrations were conducted as the ratio of arch rise to panel thickness equaled to 1. It was found that, as Mach number increased to larger values, the regular vibration had higher frequency and different modes due to aerodynamic heating.
Test on mechanics of three dimensional four directional C/C composites at elevated temperature and oxidation environment
2018, 33(8): 1916-1922. doi: 10.13224/j.cnki.jasp.2018.08.015
Abstract:
In order to research the mechanics of 3D braided C/C composites at elevated temperature and oxidation environment, the tensile test and fatigue test of 3D four directional C/C composite plate specimens, which were covered with and without anti-oxidation coating, were carried out at 700℃.The tensile test results showed that: the tensile strength and elasticity modulus of 3D braided C/C composites without anti-oxidation coating separately descended to 70.33% and 58.57%, compared with the coated one with the effect of 1 hour thermal insulation at the temperature of 700℃. And if the time increased to 2 hour, the strength declined to 44.57% while the elasticity modulus was 38.99%. The total tensile strain of specimens without anti-oxidation coating increased significantly, and the stiffness did not suddenly drop in case of material access failure. The fatigue test results show that: the residual stiffness of 3D four directional C/C composites with anti-oxidation coating was increased first, then maintained, but suddenly dropped finally. The residual strength after 105 cycles at 83% stress level was 19.75% higher than the initial strength. The residual stiffness of 3D four directional C/C composites without anti-oxidation coating increased, and then decreased until the material was destroyed completely. The residual strength after 105 cycles was 20.40% lower than the initial strength at 75% stress level.
Topological optimization design of aero-engine blade considering bird strike
2018, 33(8): 1923-1932. doi: 10.13224/j.cnki.jasp.2018.08.016
Abstract:
In order to achieve the lightweight design in consideration of the phenomenon of bird strike, it is necessary to dig into this issue and put forward a concrete solution from the perspective of structural topological optimization. Deep theoretical analysis and study were conducted to deal with those key points such as the construction of topological optimization model, the structural topological optimization while considering a type of dynamic load, and multiple constraints and cases in the dynamic problem. The theoretical basis of how to apply topological optimization method into the lightweight design was then clarified, and a computing scheme called multi-step optimizing calculation was also proposed and applied based on a corresponding numerical simulation platform. The lightweight design was based on those above key points and methods. The final optimized structure of aero-engine blade achieved its lightweight while satisfying all the optimization constraints under the different bird strike working conditions. The result showed that it could acquire about 37.9% mass reduction of the aero-engine blade. Thus, it is of great practical significance. Furthermore, the proposed topological optimization method for dynamic cases offers broad application prospect in most lightweight designs of engineering structures.
Effects of dimples with different depths on loss performance of highly-loaded compressor cascade
2018, 33(8): 1933-1940. doi: 10.13224/j.cnki.jasp.2018.08.017
Abstract:
The CFD investigation was verified by an experiment to study the aerodynamic performance of a highly-loaded compressor cascade with dimples of different depth. Five dimple configurations with diameters of 0.2,0.3,0.4,0.5,0.6mm and a depth-diameter ratio of 0.25, were arranged in 42%-60% axial chord lengths of suction side respectively, and studied to validate their effectiveness in reducing total pressure loss. The dimples with different depth exhibited different performance in loss reduction. Result showed that, a 10.8% loss reduction can be required under the influence of dimples whose depth was 0.2mm. Three-dimensional spherical dimples can promote the transition of the boundary layer by enhancing the turbulent kinetic energy level. Given the eradication of the laminar separation bubble on the suction side, an improved aerodynamic performance of linear cascade can be required.
Coupling optimization design for large-turning tandem blade shape and relative position
2018, 33(8): 1941-1953. doi: 10.13224/j.cnki.jasp.2018.08.018
Abstract:
To improve the design quality of tandem blade, an automatic optimization system of tandem blade was developed based on an improved particle swarm optimization algorithm (IPSO), adaptive Kriging model and non-uniform rational B-splines(NURBS)method. The optimization system can be used to realize the coupling optimization for shape and relative position of tandem blade. Given that Particle swarm optimization (PSO) algorithm has the advantage of fast convergence speed and may also fall into local optimal solution, an improved particle swarm optimization algorithm was proposed. It can effectively balance the global and local searching ability of PSO by adaptively changing the inertia factor, learning factor, and the number of neighborhood particles. The artificial immune operator can effectively maintain the population diversity of PSO. In addition, NURBS method was used to parameterize tandem blade, and a perturbation method of NURBS control points was designed. It was proved that the improved expected improvement (EI) criterion can make Kriging more easily jump out of the local optimal solution. The optimization system was validated by optimizing a large-turning tandem blade. Results indicate that, as compared with original tandem blade, at design condition, the total pressure loss coefficient of the optimized blade decreased by 40.4%. Besides, the static pressure ratio of optimized blade was higher and the total pressure loss coefficient was smaller at all incidence conditions. The performance of optimized blade was largely improved at positive incidence. It also proved that the coupling optimization design method of tandem blade had a good application value.
Body force model for multistage axial compressor—theoretical method and simplified application
2018, 33(8): 1954-1963. doi: 10.13224/j.cnki.jasp.2018.08.019
Abstract:
Based on the time marching technique and the finite volume method, a three-dimensional numerical computational model was proposed through the combination of the body force model and elementary cascade method in order to predict the impact of the large-scale inlet distortion on multistage axial compressor aerodynamic performance and stability. Under the premise of axisymmetric inlet condition, the corresponding simplified computational model was applied to analyze the internal flow field of a four-stage axial compressor under uniform and tip radial total pressure distorted inlet in detail. The calculated result with clean inlet was consistent with the related experimental result. The calculated result with tip radial total pressure distorted inlet indicated that a significant radial mixing was caused by the radial distortion in the compressor passage,and the tip radial total pressure distortion could greatly worsen the compressor stability. The predicted stability boundary with tip radial inlet total pressure distortion gradually decreased with the increase of distortion intensity. The calculated results support the possibility of using the model to predict the radial inlet distortion impact.
Mechanism of stall development in a transonic axial compressor
2018, 33(8): 1964-1973. doi: 10.13224/j.cnki.jasp.2018.08.020
Abstract:
Three-dimensional multi-passage numerical simulations on the NASA Rotor 37 were conducted to investigate the evolution process of rotating stall in the transonic axial compressor. The development of the shock and the spillage flow and their impacts on the stall process were analyzed minutely. The results show that only one oblique shock was observed in the passage under peak efficiency condition, and the oblique shock evolved into a detached shock under stall condition. The tip leakage vortex broke down after interacting with the detached shock. With the axial development of the broken tip leakage vortex, because of the adverse pressure gradient, a significant vortex was formed in the middle of the passage. At the beginning of the stall process, the spillage flow appeared intermittently under the influence of the periodical evolution of the blockage region. Accompanied by the decrease of the flow rate, the blockage region extended gradually. At the same time, the spillage flow existed during the whole periodical evolution process of the blockage region, which can be treated as the onset of a rapid deterioration of the flow field.
An effective method for numerical simulation of individual components of pitching combined dynamic derivative
2018, 33(8): 1974-1980. doi: 10.13224/j.cnki.jasp.2018.08.021
Abstract:
Based on the field velocity approach, an effective method was presented for numerical simulating of floating motion, pure pitch motion and pitching oscillation, and then the lag of wash derivative and the pitch-damping derivative was directly calculated, without the mesh movement. Without the need to update the grid in real time, this approach reduced the time and memory required for computation, and avoided the negative volume. Numerical results were validated by comparison with the experiment results for NACA 0006 airfoil with gust response and NACA 0012 airfoil with forced oscillation. To validate the applicability for the present method farther, pitching combined dynamic derivative and its individual components were calculated from load history of the unsteady flow around a standard research configuration, known as the basic finner missile(BFM). Predicted results showed a good agreement with the available wind tunnel data, the maximum error was no more than 4%.
Roughness for wall turbulence extension model and flow numerical simulation
2018, 33(8): 1981-1989. doi: 10.13224/j.cnki.jasp.2018.08.022
Abstract:
Reynolds average Navier-Stokes equations were solved based on finite volume method to simulate the surface roughnesss influence on aerodynamic and thermodynamic characteristics. To account for the effect of wall roughness, the wall-value of ω of Knopps roughness extension for shear stress transport k-ω two equations turbulence model was modified to increase near wall eddy viscosity and total wall shear stress. Flows over flat plate and past NACA 652215 airfoil were simulated and the influences of roughness on skin friction coefficient, Stanton number and lift coefficient were analyzed. The predicted results showed the roughness increased the skin friction coefficient and Stanton number, shifted separation upstream and decreased the lift coefficient. The shift of velocity profile matched well with theoretical results when flow was in transitional regime compared with the roughness extension by Knopp, meanwhile, skin friction coefficient, Stanton number and life coefficent agreed better with experiments. The extension model allows for the simulation of flows over rough surfaces quickly and preciously at the same grid resolution requirements as for smooth walls.
Mechanism of aerodynamic interference of tankers jet flow on receiver
2018, 33(8): 1990-2005. doi: 10.13224/j.cnki.jasp.2018.08.023
Abstract:
The jet effect in aerial refueling was studied in hose-drogue after body (HDB), hose-drogue after pod (HDP) and flying boom (FB) refueling methods, based on structured multi-block grid and Reynolds-averaged Navier-Stokes(RANS) equations. This numerical method was verified by turbine powered simulation (TPS) and DLR-F6 cases. Compared with no jet situation, the mechanism of aerodynamic interference was explored in three different methods. Result showed that, the lift coefficient, drag coefficient and pitch down moment of receiver were larger than the case without jet flow. The receiver in FB refueling method was slightly affected, while the lateral characteristic in HDP refueling method was changed severely by jet. The increase of dynamic pressure, the change of local flow angle and the ejecting effect were main phenomena caused by jet flow, which resulted in remarkable pressure change on receiver. Moreover, the local flow angle was changed by acceleration effect and the wake vortex dissipated by jet.
Numerical investigation of influence of co-flow velocity on plane jet
2018, 33(8): 2006-2015. doi: 10.13224/j.cnki.jasp.2018.08.024
Abstract:
The characteristics of plane jet were investigated by large eddy simulation, at the outlet of plane jet Mach number of 0.9 and the co-flow velocity ratio of 0.1, 0.3 and 0.5, respectively. High resolution numerical method and Smagorinsky sub-grid scale model were used. The mean flow properties, fluctuations and the evolution of the vortical structures were analyzed. Result showed that with the growth of the co-flow velocity, the potential core length increased while the development of the shear layer was slowed down and the jet transition was delayed. The velocity distribution was self-similar, but possible similarity in the turbulence intensity was not found. Two-point space-time correlations of the velocity and pressure fluctuations in the shear layer were investigated. It was found that the increase of the co-flow velocity weakened the space-time correlations, but accelerated the spreads of fluctuations to the downstream. This study provides a basis for revealing the effect of co-flow velocity on the sound field.
Characteristic of helium cycle system parameters for pre-cooling air turbo rocket engine
2018, 33(8): 2016-2024. doi: 10.13224/j.cnki.jasp.2018.08.025
Abstract:
Pre-cooling air turbo rocket engine (PATR) is a typical scheme. A mathematical model of helium cycle system for PATR was established by considering the main factors, such as the structural features of components, geometry dimensions and working fluid properties. Numerical simulation showed that the thrust and specific impulse of PATR simulation were better, and the model could describe the thermodynamic process of helium cycle. In addition, the excess air coefficient of engine was a key parameter affecting the engine thrust and specific impulse significantly. Analysis results indicate that improving the highest pressure in helium system and decreasing the inlet helium temperature in air pre-cooler are effective methods to reduce the system pressure loss and helium compressor power. For every increase of the highest pressure with 1MPa, the average system pressure loss and the helium compressor power decreased 1.1% and 3.2%, respectively. Similarly,for every increase of air pre-cooler inlet helium temperature with 1K, the average system pressure loss and helium compressor power increased 0.086% and 2.3%, respectively.
Optimization for liquid film cooling based on orthogonal methodology and Kriging model in liquid rocket engine
2018, 33(8): 2025-2032. doi: 10.13224/j.cnki.jasp.2018.08.026
Abstract:
Transcritical methane film cooling steady flow field was calculated by numerical transfer study in the rocket engine chamber by FLUENT. Through the orthogonal methodology, cooling performance under the combined action of different diameters, axial angle, radial angle and shape of hole was simulated to choose the optimal geometric parameters. Based on the optimal geometric parameters, the optimal latin hypercube design was adopted to get sample points that can be used to establish Kriging model. Then the optimal transcritical liquid film mass flow, the partition ratio and location of cooling ring under multi-objective conditions can be calculated through genetic algorithm. The results showed that orthogonal methodology and Kriging model could resolve problems of high design cost and numerical noise in liquid rocket engine design. In the considered geometrical factors above, the significance order from big to small effects on cooling efficiency and non-uniformity degree was: shape, diameter, radial angle, axial angle; the optimal geometric parameter combination was 0.003mm diameter, 45° axial angle, 15° radial angle, and diffused shape. The Kriging model established precisely reflected the relation between objective function and the liquid film mass flow, partition ratio and location of cooling ring. The average cooling efficiency and specific impulse loss increased by 4.9% and 0.37%, respectively, and non-uniformity degree was reduced by 0.025, while the total objective function was increased by 184%. After the optimization, asymmetry of vortex pair made coolant spanwise well-distributed, while reversed vortex pair attenuated faster, and the property of liquid film attaching to the wall increased, thus improving the cooling efficiency.
Fault feature extraction method of rolling bearing based onspectral graph indices
2018, 33(8): 2033-2040. doi: 10.13224/j.cnki.jasp.2018.08.027
Abstract:
In order to more accurately extract nonlinear fault features of the rolling bearing vibration signals, the graph signal processing(GSP) was introduced into the field of mechanical fault diagnosis, and a rolling bearings fault feature extraction method based on spectral graph indices was proposed. In the proposed method, the rolling bearing vibration signal was converted into the path graph, from which a number of spectral graph indices were extracted. The spectral graph indices were sorted by the Fisher score (FS) algorithm, and some of the most sensitive spectral graph indices were chosen as rolling bearings fault feature parameters. The K-means clustering algorithm was used to identify different faults of rolling bearings. The application examples indicated that, when 1 to 5 optimal spectral graph indices, time-domain indices and frequency-domain indices were respectively selected to identify different bearing faults, the spectral graph indices had no identification errors, while both the time-domain indices and the frequency-domain indices had different numbers of identification errors. Therefore, the distinguishing ability to bearing faults of the spectral graph indices was superior to those of the time-domain indices and the frequency-domain indices.
Assessment of performance degradation for aero-engine based on denoising autoencoder
2018, 33(8): 2041-2048. doi: 10.13224/j.cnki.jasp.2018.08.028
Abstract:
Targeting the form and law of the aero-engines performance degradation, a degradation assessment method based on denoising autoencoder was proposed. On account of the collected six aero-engine condition monitoring parameters, denoising autoencoder and greedy layer-wise training algorithm were used to assess performance degradation in order to explore the deep influences of those parameters on engines performances and extract the data characteristics more conducive to the assessment. The comparison between the proposed algorithm and back propagation(BP) neural networks as well as support vector machine showed that the proposed method had high accuracy and robustness. The accuracy of the proposed method was 93.5% and the accuracy just reduced to 84.5% when signal-to-noise ratio was 10dB. The proposed method also can prevent the over-fitting of small samples in aero-engine condition monitoring.