Mathematical modeling of mixed convective MHD Falkner-Skan squeezed Sutterby multiphase flow with non-Fourier heat flux theory and porosity

In a wide variety of mechanical and industrial applications,e.g.,space cooling,nuclear reactor cooling,medicinal utilizations(magnetic drug targeting),energy generation,and heat conduction in tissues,the heat transfer phenomenon is involved.Fourier’s law of heat conduction has been used as the foundation for predicting the heat transfer behavior in a variety of real-world contexts.This model’s production of a parabolic energy expression,which means that an initial disturbance would immediately affect the system under investigation,is one of its main drawbacks.Therefore,numerous researchers worked on such problem to resolve this issue.At last,this problem was resolved by Cattaneo by adding relaxation time for heat flux in Fourier’s law,which was defined as the time required to establish steady heat conduction once a temperature gradient is imposed.Christov offered a material invariant version of Cattaneo’s model by taking into account the upper-connected derivative of the Oldroyd model.Nowadays,both models are combinedly known as the Cattaneo-Christov(CC)model.In this attempt,the mixed convective MHD Falkner-Skan Sutterby nanofluid flow is addressed towards a wedge surface in the presence of the variable external magnetic field.The CC model is incorporated instead of Fourier’s law for the examination of heat transfer features in the energy expression.A two-phase nanofluid model is utilized for the implementation of nano-concept.The nonlinear system of equations is tackled through the bvp4c technique in the MATLAB software 2016.The influence of pertinent flow parameters is discussed and displayed through different sketches.Major and important results are summarized in the conclusion section.Furthermore,in both cases of wall-through flow(i.e.,suction and injection effects),the porosity parameters increase the flow speed,and decrease the heat transport and the influence of drag forces.

Characterization of Passive Films Formed on As-received and Sensitized AISI 304 Stainless Steel

The current research of corrosion resistance of stainless steels mainly focuses on characterization of the passive films by point defect mode and mixed-conduction model. The corrosion resistance of the passive films formed on asreceived and sensitized AISI304 stainless steel in borate bu er solution were evaluated in this paper. The degree of sensitization and corrosion resistance of AISI304 stainless steels was evaluated by double loop electrochemical potentiodynamic reactivation and electrochemical impedance spectroscopy. The passive films formed on the stainless steels were studied by XPS technique. It was found that as-received specimen had higher pitting corrosion potential and corrosion resistance than sensitized one. The Mott-Schottky results showed that sensitized stainless steel had more defects in the passive film than as-received one. The compositions of the passive films were mainly Cr and Fe oxides according to XPS results.

Investigation of high-speed abrasion behavior of an abradable seal rubber in aero-engine fan application

Abradable seal rubber has been widely used in aero-engine fans to improve their efficiency by reducing the clearance between rotating and stationary components. To investigate the high-speed scraping behavior between a vulcanized silicone rubber and a Ti-6Al-4V fan blade and evaluate the abradable performance of seal rubber, abrasion tests were conducted at a blade tip velocity of 50–300 m/s with an incursion rate of 100 lm/s. The influences of the blade tip velocity on the wear mechanism and interaction forces were specially analyzed. It is shown that abrasive wear and pattern wear are the predominant wear mechanisms, and pattern wear can be seen as the emergence and propagation of cracks. With an increase of the blade tip velocity, both of the final incursion depth and wear mass loss of seal rubber exhibit growth trends. The gradual changes of rubbing forces with an increase of rubbing time are the characteristic of abrasive wear, and force curves with unstable mutations are a reflection of pattern wear. At a constant incursion rate of 100 lm/s, the maximum values of interaction forces decrease first and then grow with an increase of the blade tip velocity, and the blade tip velocity of 150 m/s becomes the cut-off point between abrasive wear and pattern wear.

Recent development of casing treatments for aeroengine compressors

Casing treatment is a mature stabilization technique which has been widely applied on aero-engines for modern aircrafts and turbo-chargers for automobiles. After the investigations of half century since the 1960 s, this technique has been well developed for various configurations with different effectiveness. From the perspective of stabilization mechanism, this paper roughly categorizes the configurations of casing treatment into two types: traditional ones which work by affecting the flow structure of blade tip region; a novel one named as Stall Precursor-Suppressed(SPS) casing treatment. The effectiveness of both types will be demonstrated for their applications on axial compressors and centrifugal compressors with uniform or distorted inlet. And the stabilization mechanism of casing treatments for regular types and SPS one will also be explained respectively. In addition, this review will summarize the methodologies of casing treatments with the numerical simulations for regular grooved configurations and the eigenvalue approach for SPS casing treatment.Looking forward to the future of compressor stabilization, casing treatment technique will still exist as a general and inexpensive option, and the exploration for its effectiveness and mechanism will be deeper with the development of computational fluid dynamics and advanced measurement techniques.

Effects of tooth bending damage on the leakage performance and rotordynamic coefficients of labyrinth seals

Tooth bending damage resulting from an intense impact by the rotor sometimes occurs in the transient operation.To investigate the influence of after-damage clearance and tooth bending length on the leakage performance and rotordynamic coefficients of labyrinth seals,three tooth bending damages were taken into consideration,including the unbent tooth damage(abbreviated as Unbent),the partial tooth bending damage(abbreviated as Pbent)and the complete tooth bending damage(abbreviated as Cbent).The transient CFD solution was utilized to calculate the leakage flow rates and rotordynamic coefficients of labyrinth seals with clearances of 0.3,0.4,0.5,0.6 mm for three tooth bending damages.The obtained result shows that the Unbent tooth damage leaks least while the Pbent tooth bending damage leaks most,and an increase of 6.1%for Cbent tooth bending damage and an increase of 19.4%for Pbent tooth bending damage are discovered at the tooth clearance of 0.6 mm in comparison with the Unbent tooth damage.Compared to the Unbent tooth damage,the effective damping for Pbent tooth bending damage and Cbent tooth bending damage is lower and drops by 9.7%–33.6%and 8.5%–22.6%respectively at the tooth clearance of 0.6 mm,suggesting that Pbent tooth bending damage or Cbent tooth bending damage tends to weaken the seal stability when compared to the Unbent tooth damage.

Experimental study on NO_x emission correlation of fuel staged combustion in a LPP combustor at high pressure based on NO-chemiluminescence

Experimental investigations on NOx emissions of a single-cup,Lean Premixed Prevaporized(LPP),module combustor were carried out at elevated inlet temperature and pressure up to810 K and 2.0 MPa,close to the real operating conditions of aero-engine combustors.This LPP combustor adopts centrally staged fuel injections which could produce separated stratified swirling spray flame.In the NOx emissions measurements,the ranges of dome equivalence ratio and fuel stage ratio were from 0.55 to 0.58 and 8%to 24%,respectively.The optical diagnosis on separated stratified swirling spray flame were carried out with fuel stage ratio changing from 15%to 30%.Therefore,NO*and OH*chemiluminescence images were obtained.The results show that NOx emissions increase with the increase of the fuel stage ratio.And from the chemiluminescence images,the main flame and pilot flame are found weakly coupled.The pilot flame plays a significant role in NOx emission production because of its higher adiabatic flame temperature.Based on the results of chemiluminescence optical tests,a new NOx emission prediction model is proposed based on the Lefebvre’s single flame model The estimate of local equivalence ratio of the pilot stage’s nonpremixed flame is modified considering the characteristics of spray combustion,and a"PLUS"emission prediction model suitable for separated stratified swirling spray flame is obtained.Compared to the experimental data,the"PLUS"model exhibits a good prediction in a range of±13%of deviation.

A dual timescale model for micro-mixing and its application in LES-TPDF simulations of turbulent nonpremixed flames

The numerical simulation of modern aero-engine combustion chamber needs accurate description of the interaction between turbulence and chemical reaction mechanism. The Large Eddy Simulation(LES) method with the Transported Probability Density Function(TPDF) turbulence combustion model is promising in engineering applications. In flame region, the impact of chemical reaction should be considered in TPDF molecular mixing model. Based on pioneer research, three new TPDF turbulence-chemistry dual time scale molecular mixing models were proposed tentatively by adding the chemistry time scale in molecular mixing model for nonpremixed flame. The Aero-Engine Combustor Simulation Code(AECSC) which is based on LES-TPDF method was combined with the three new models. Then the Sandia laboratory's methane-air jet flames: Flame D and Flame E were simulated. Transient simulation results show that all the three new models can predict the instantaneous combustion flow pattern of the jet flames. Furthermore,the average scalar statistical results were compared with the experimental data. The simulation result of the new TPDF arithmetic mean modification model is the closest to the experimental data:the average error in Flame D is 7.6% and 6.6% in Flame E. The extinction and re-ignition phenomena of the jet flames especially Flame E were captured. The turbulence time scale and the chemistry time scale are in different order in the whole flow field. The dual time scale TPDF combustion model has ability to deal with both the turbulence effect and the chemistry reaction effect, as well as their interaction more accurately for nonpremixed flames.

Experimental study on film cooling performance of imperfect holes

An experimental study is made to investigate the film cooling performance of imperfect holes due to in-hole blockage over a flat plate. A specifically pyramid-shaped element is used to simulate the in-hole blockage. Six in-hole blockage orientations（such as leading-inlet, leading-middle,leading-exit, trailing-inlet, trailing-middle and trailing-exit） and four blocking ratios（ranging from 0.1 to 0.4） are taken into considerations. Based on the experimental results, the influences of in-hole blockage on the film cooling effectiveness and discharge coefficient under typical blowing ratios are analyzed. It is confirmed that the in-hole blockage results in a reduction of discharge coefficient related to the perfect film cooling holes, especially for the leading-exit and trailing-inlet orientations with a big blocking ratio. However, in the view of film cooling effectiveness, the in-hole blockage shows complicated affecting roles. In general, except for the leading-exit orientation, the in-hole blockages produce detrimental influence on the film cooling effectiveness.

An experimental method for squealer tip flow field considering relative casing motion

Squealer tip is widely used in turbines to reduce tip leakage loss.In typical turbine environment,the squealer tip leakage flow is affected by multiple factors such as the relative casing motion and the wide range of variable incidence angles.The development of experimental methods which can accurately model the real turbine environment and influencing factors is of great significance to study the squealer tip leakage flow mechanism.In the present paper,a low-speed turbine cascade test facility which can model the relative casing motion and wide range of variable incidence angles(-25°to 55°)is built.Based on the similarity criteria,a high-low speed similarity transformation method of the turbine cascade is established by considering the thickness of the turbine blade.A combined testing method of Particle Image Velocimetry(PIV)and local pressure measurement is proposed to obtain the complex flow structures within the tip cavity.The results show that the experimental method can successfully model the relative casing motion and the wide range of variable incidence angles.The low-speed cascade obtained by the similarity transformation can model the high-speed flow accurately.The measurement technique developed can obtain the complex flow field and successfully capture the scraping vortex within the squealer tip.

Containment of soft wall casing wrapped with Kevlar fabric

Aramid fabrics have been commonly used in the civil turbofan engine fan blade containment system for its excellent performance. To investigate the behavior and capability of soft wall containment casing, a series of fan blade released tests were conducted on the high-speed spin tester.The soft wall casing was fabricated by wrapping multiple layers of Kevlar49 plain woven fabric around a thin steel ring. Casings with different inner metal ring and outer fabric layers number were compared. The method of using the explicit dynamic software LS-DYNA to establish the finite element analysis model for the quantitative analysis of the containment process was developed and conducted. The simulation results are in good agreement with the test results. It is shown that the containment process of the soft wall casing can be divided into three impact stages. The casing with low-stiffness inner metal ring will get severe overall deformation and lose the structural integrity when it suffers the blade impact. Kevlar fabric layers will appear large bulge on outside surface and absorb the most impact dynamic energy of the high speed released fan blade. By summing up the results of the test and simulation, an empirical critical equation was derived to describe the relationship between the released blade dynamic energy and the Kevlar fabric thickness.

Two freedom linear parameter varying μ synthesis control for flight environment testbed

To solve the problem of robust servo performance of Flight Environment Testbed(FET)of Altitude Ground Test Facilities(AGTF) over the whole operational envelope, a two-degree-offreedom μ synthesis method based on Linear Parameter Varying(LPV) schematic is proposed, and meanwhile a new structure frame of μ synthesis control on two degrees of freedom with double integral and weighting functions is presented, which constitutes a core support part of the paper. Aimed at the problem of reference command's rapid change, one freedom feed forward is adopted, while another freedom output feedback is used to meet good servo tracking as well as disturbance and noise rejection; furthermore, to overcome the overshoot problem and acquire dynamic tuning,the integral is introduced in inner loop, and another integral controller is used in outer loop in order to guarantee steady errors; additionally, two performance weighting functions are designed to achieve robust specialty and control energy limit considering the uncertainties in system. As the schedule parameters change over large flight envelope, the stability of closed-loop LPV system is proved using Lyapunov inequalities. The simulation results show that the relative tracking errors of temperature and pressure are less than 0.5% with LPV μ synthesis controller. Meanwhile, compared with non-LPV μ synthesis controller in large uncertainty range, the proposed approach in this research can ensure robust servo performance of FET over the whole operational envelope.

Transfer element method with application to acoustic design of aeroengine nacelle

In the present survey, various methods for the acoustic design of aeroengine nacelle are first briefly introduced along with the comments on their advantages and disadvantages for practi- cal application, and then detailed analysis and discussion focus on a kind of new method which is called ＂transfer element method＂ （TEM） with emphasis on its application in the following three problems： turbomachinery noise generations, sound transmission in ducts and radiation from the inlet and outlet of ducts, as well as the interaction between them. In the theoretical frame of the TEM, the solution of acoustic field in an infinite duct with stator sound source or liner is extended to that in a finite domain with all knows and unknowns on the interface plane, and the relevant acoustic field is solved by setting up matching equation. In addition, based on combining the TEM with the boundary element method （BEM） by establishing the pressure and its derivative con- tinuum conditions on the inlet and outlet surface, the sound radiation from the inlet and outlet of ducts can also be investigated. Finally, the effects of various interactions between the sound source and acoustic treatment have been discussed in this survey. The numerical examples indicate that it is quite important to consider the effect of such interactions on sound attenuation during the acoustic design of aeroengine nacelle.

Research on vibration suppression of a mistuned blisk by a piezoelectric network

The work aims to provide a further investigation of the dynamic characteristics of an integral bladed disk（also called ‘blisk＇） with a Parallel Piezoelectric Network（PPN）. The PPN is constructed by parallelly interconnecting the piezoelectric patches distributed in the blisk. Two kinds of PPN are considered, namely mono-periodic PPN and bi-periodic PPN. The former has a piezoelectric patch in each sector, and the later has one patch every few sectors. The vibration suppression performance of both kinds of PPN has been studied through modal analysis, forced response analysis, and statistical analysis. The research results turn out that the PPN will only affect mechanical frequencies near the electrical frequency clusters slightly, and the bi-periodic PPN will make the nodal diameter spectrum of the modes more complex, but the amplitude corresponding to the new nodal diameter component is much smaller than that of the nodal diameter component corresponding to the mono-periodic system. The mechanical coupling between the blades and the disk plays an important role in the damping effect of the PPN, and it should be paid attention to in applications. The mono-periodic PPN can effectively suppress the amplitude magnification of the forced response induced by the mistuning of the blisk; meanwhile, it can mitigate the vibration localization of the mistuned electromechanical system. If piezoelectric patches are set only in part of the sectors, the bi-periodic PPN still has a vibration suppression ability, but the effect is related to the number and spatial distribution of the piezoelectric patches.

Containment ability and groove depth design of U type protection ring

High-energy rotor uncontained failure can cause catastrophic damage effects to aircraft systems if not addressed in design. In this paper, numerical simulations of three high-energy rotor disk fragments impacting on U type protection rings are carried out using LS-DYNA. Protection rings with the same mass and different groove depths are designed to study the influence of the groove depth. Simulation results including kinetic energy and impact force variation of single frag- ment are presented. It shows that the groove depth infects both the axial containment ability of the protection ring and the transfer process of energy. The depth of groove ought to be controlled to an appropriate value to meet both the requirement of axial containment and higher safety factor. Ver- ification test on high-speed spin tester has been conducted and shows that protection ring with appropriate U structure can resist the impact of the disk burst fragments. The ring is inflated from a circular to an oval-triangle shape. The corresponding simulation shows good agreement with the test.

Effect of high temperature on compression property and deformation recovery of ceramic fiber reinforced silica aerogel composites

Ceramic fiber reinforced silica aerogel composites are novel insulation materials in the thermal protection field for hypersonic vehicles. Before the aerogel composites are applied in load-bearing structures, it is necessary to investigate their mechanical properties including load-bearing and deformation recovery capabilities. High temperature from service conditions will have important effects on the mechanical properties of thermal protection materials. In this paper, compression tests including loading and unloading stages were conducted for ceramic fiber reinforced silica aerogel composites at room temperature and elevated temperatures(300℃, 600℃ and 900℃). Influences of thermal exposure to high temperature and high temperature service environment on the compression property and deformation recovery were both investigated. Scanning electron microscopy(SEM), Fourier transform infrared spectroscopy(FT-IR) and X-ray diffraction(XRD) were applied to help understand the mechanisms of mechanical property variations. The experimental results show that the compression modulus and strength both increase with the increasing thermal exposure temperature and testing temperature,but the deformation recovery capability decreases. The micro structure changes caused by thermal sintering are considered as the main reason for the property variations.Viscous flow and matter transport due to high temperature resulted in the fusion of aerogel particles. This made the particle skeleton thicker and stronger, which led to higher stiffness and strength of the composites. However, matrix cracks induced by the formation and fracture of larger pores made unrecoverable deformation more serious. In the tests at elevated temperatures,the aggregation of aerogel particles in a fused state got more severe because of the addition of mechanical load. As a result, the degradation of deformation recovery capability became more significant.

Effect of bristle pack position on the rotordynamic characteristics of brush-labyrinth seals at various operating conditions

Over the last few decades,the research on the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals is not sufficient.To this end,two kinds of brush-labyrinth seals for the bristle pack element installed upstream of the labyrinth teeth named BSU and installed downstream of the labyrinth teeth called BSD were used to investigate the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals.Using the numerical model combining the porous medium model and the whirling rotor method,the rotordynamic characteristics of the BSU and BSD at various operating conditions including four kinds of pressure ratios,five kinds of inlet preswirl speeds and four kinds of rotor spinning speeds were conducted.The obtained results show that the effects of operating conditions on rotordynamic coefficients for the different seal configurations are different.The direct stiffness,cross-coupled stiffness and direct damping of the BSU are lower than those of the BSD.The rotordynamic coefficients of the BSU are more insensitive to the operating conditions variation.From the perspective of the seal stability,the BSU is a better brush-labyrinth seal configuration at high pressure ratio,high positive preswirl or high rotor spinning speed conditions.While in the case of low pressure ratio,low positive preswirl or low rotor spinning speed conditions,the BSD is a better choice.

Dynamic characteristics of a cyclic-periodic structure with a piezoelectric network

This paper deals with a cyclic-periodic structure with a piezoelectric network. In such a system, there is not only mechanical connection but also electrical connection between adjacent periodic sectors. The objective is to learn whether the presence of a piezoelectric network would change the dynamic characteristics of the system. The background of the research is about vibration reduction of a bladed disk in an aero-engine, and the system is simulated by a lumped parameter model. The dynamic equations of the system are derived, and then the analytical solution corresponding to the eigenvalue problem is given. The vibration responses to single traveling wave excitations （EO excitations） and multiple traveling wave excitations （NEO excitations） are studied. The results show that the presence of a piezoelectric network would change the natural frequencies of the system compared with those of the system with the piezoelectric shunt circuit. The forced response is sensitive to the connection type and the elements of the network. An energy analysis of the electro-mechanical coupling system has been performed to understand its dynamic behavior, and the following conclusion is obtained： a vibration reduction to excitations whose primary har- monic component is not zero can be achieved by a parallel piezoelectric network, while a reduction to other excitations should be based on a series piezoelectric network.

Turbulent characteristics and rotation correction of wall function in rotating channel with high local rotation parameter

The turbulent fluctuation and the rotation correction of wall function law are investigated in the entrance section of a rotating channel. The one-dimensional hot wire probe and the X-type probe are utilized to measure the boundary layer at four streamwise stations. Through the analysis on the boundary layer near the leading side and trailing side, it is found that the turbulent fluctuation is promoted in the trailing side whereas suppressed in the leading side. This difference is attributed to the Coriolis instability near the trailing side. In addition, considering the local rotation parameter Rc, whose maximum absolute value is 0.014, is larger than that in previous research, whose maximum value is 0.007, the whole process of the relaminarization is captured. To understand this phenomenon better, the effects of the generation term and the Coriolis term in the transport equation of the Reynolds stress are discussed. In addition, the rotation correction of the viscous-Coriolis region and the Coriolis region are discussed, a new revising method for the wall function is proposed.

Effect of rotor-mounted protrusion on sealing performance and flow structure in rotor-stator cavity

This paper presents an experimental investigation on the effect of protrusion radial position and height on the sealing performance and flow structure in the rotor-stator cavity. The rotormounted protrusions are assembled at three radial positions and are set to three heights. The cavity is equipped with three rim seals: a radial seal, an axial seal and a seal with double fins on the stator.The annulus Reynolds number is set at 4:39 ×10;and the rotational Reynolds number ranges from 7:51×10;to 1:20×10;. Heat and mass transfer analogy is applied. Pressure and CO;concentration are measured. The experimental results show that in cavities with different rim seals, radial distributions of the sealing efficiency, pressure and swirl ratio are basically the same. The sealing performance is improved by protrusions compared with the cavity without protrusion and improves with the increase of protrusion radial position and height. The effect of protrusion increases with the increase of the rotational Reynolds number. The windage loss and the flow resistance introduced by protrusions are investigated. It is found that induced windage loss and flow resistance decrease with the increase of protrusion radial position but increase with the protrusion height.

Adaptive sliding mode control for limit protection of aircraft engines

In practice, some sensors of aircraft engines naturally fail to obtain an acceptable measurement for control propose, which will severely degrade the system performance and even deactivate the limit protection function. This paper proposes an adaptive strategy for the limit protection task under unreliable measurement. With the help of a nominal system, an online estimator with gradient adaption law and low-pass filter is devised to evaluate output uncertainty.Based on the estimation result, a sliding mode controller is designed by defining a sliding surface and deriving a control law. Using Lyapunov theorem, the stability of the online estimator and the closed-loop system is detailedly proven. Simulations based on a reliable turbofan model are presented, which verify the stability and effectiveness of the proposed method. Simulation results show that the online estimator can operate against the measurement noise, and the sliding controller can keep relevant outputs within their limits despite slow-response sensors.