The rapidly increasing penetration of electric vehicles(EVs) in modern metropolises has been witnessed during the past decade, inspired by financial subsidies as well as public awareness of climate change and environm...The rapidly increasing penetration of electric vehicles(EVs) in modern metropolises has been witnessed during the past decade, inspired by financial subsidies as well as public awareness of climate change and environment pro-tection. Integrating charging facilities, especially highpower chargers in fast charging stations, into power distribution systems remarkably alters the traditional load flow pattern, and thus imposes great challenges on the operation of distribution network in which controllable resources are rare. On the other hand, provided with appropriate incentives, the energy storage capability of electric vehicle offers a unique opportunity to facilitate the integration of distributed wind and solar power generation into power distribution system. The above trends call for thorough investigation and research on the interdependence between transportation system and power distribution system. This paper conducts a comprehensive survey on this line of research. The basic models of transportation system and power distribution system are introduced,especially the user equilibrium model, which describes the vehicular flow on each road segment and is not familiar to the readers in power system community. The modelling of interdependence across the two systems is highlighted.Taking into account such interdependence, applications ranging from long-term planning to short-term operation are reviewed with emphasis on comparing the description of traffic-power interdependence. Finally, an outlook of prospective directions and key technologies in future research is summarized.展开更多
In order to investigate the micro-process and inner mechanism of rock failure under impact loading, the laboratory tests were carried out on an improved split Hopkinson pressure bar (SHPB) system with synchronized m...In order to investigate the micro-process and inner mechanism of rock failure under impact loading, the laboratory tests were carried out on an improved split Hopkinson pressure bar (SHPB) system with synchronized measurement devices including a high-speed camera and a dynamic strain meter. The experimental results show that the specimens were in the state of good stress equilibrium during the post failure stage even when visible cracks were forming in the specimens. Rock specimens broke into strips but still could bear the external stress and keep force balance. Meanwhile, numerical tests with particle flow code (PFC) revealed that the failure process of rocks can be described by the evolution of micro-fractures. Shear cracks emerged firstly and stopped developing when the external stress was not high enough. Tensile cracks, however, emerged when the rock specimen reached its peak strength and played an important role in controlling the ultimate failure during the post failure stage.展开更多
As one of the most important components of river mechanics,sediment transport capacity of sediment-laden flows has attracted much attention from many researchers working on river mechanics and hydraulic engineering. B...As one of the most important components of river mechanics,sediment transport capacity of sediment-laden flows has attracted much attention from many researchers working on river mechanics and hydraulic engineering. Based on the time-averaged equation for a turbulent energy equilibrium in solid and liquid two-phase flow,an expression for the efficiency coefficient of suspended load movement was derived for the two-dimensional,steady,uniform,fully-developed turbulent flow. A new structural expression of sediment transport capacity was achieved. Using 115 runs of flume experimental data,which were obtained through two kinds of sediment transport experiments in the state of equilibrium,in combination with the basic rheological and sediment transporting characteristics of hyperconcentrated flow,the main parameters in the structural expression of sediment transport capacity were calibrated,and a new formula of sediment transport capacity for hyperconcentrated flow was developed. A large amount of field data from the Yellow River,Wuding River,and Yangtze River,etc. were adopted to verify the new formula and good agreement was obtained. These results above contribute to an improved theoretical system of river mechanics and a reliable tool for management of rivers carrying high concentration of sediments.展开更多
With discretized particle velocity space,a multi-scale unified gas-kinetic scheme for entire Knudsen number flows has been constructed based on the kinetic model in one-dimensional case[J.Comput.Phys.,vol.229(2010),pp...With discretized particle velocity space,a multi-scale unified gas-kinetic scheme for entire Knudsen number flows has been constructed based on the kinetic model in one-dimensional case[J.Comput.Phys.,vol.229(2010),pp.7747-7764].For the kinetic equation,to extend a one-dimensional scheme to multidimensional flow is not so straightforward.The major factor is that addition of one dimension in physical space causes the distribution function to become two-dimensional,rather than axially symmetric,in velocity space.In this paper,a unified gas-kinetic scheme based on the Shakhov model in two-dimensional space will be presented.Instead of particle-based modeling for the rarefied flow,such as the direct simulation Monte Carlo(DSMC)method,the philosophical principal underlying the current study is a partial-differential-equation(PDE)-based modeling.Since the valid scale of the kinetic equation and the scale of mesh size and time step may be significantly different,the gas evolution in a discretized space is modeled with the help of kinetic equation,instead of directly solving the partial differential equation.Due to the use of both hydrodynamic and kinetic scales flow physics in a gas evolution model at the cell interface,the unified scheme can basically present accurate solution in all flow regimes from the free molecule to the Navier-Stokes solutions.In comparison with the DSMC and Navier-Stokes flow solvers,the current method is much more efficient than DSMC in low speed transition and continuum flow regimes,and it has better capability than NS solver in capturing of non-equilibrium flow physics in the transition and rarefied flow regimes.As a result,the current method can be useful in the flow simulation where both continuum and rarefied flow physics needs to be resolved in a single computation.This paper will extensively evaluate the performance of the unified scheme fromfreemolecule to continuum NS solutions,and fromlow speedmicro-flow to high speed non-equilibrium aerodynamics.The test cases clearly demonstrate tha展开更多
This paper explores the application of noncooperative game theory together with the concept of Nash equilibrium to the investigation of some basic problems on multi-scale structure, especially the meso-scale structure...This paper explores the application of noncooperative game theory together with the concept of Nash equilibrium to the investigation of some basic problems on multi-scale structure, especially the meso-scale structure in the multi-phase complex systems in chemical engineering. The basis of this work is the energy-minimization-multi-scale (EMMS) model proposed by Li and Kwauk (1994) and Li, et al. (2013) which identifies the multi-scale structure as a result of 'compromise-in-competition between dominant mechanisms' and tries to solve a multi-objective optimization problem. However, the existing methods often integrate it into a problem of single objective optimization, which does not clearly reflect the 'compromise-in-competition' mechanism and causes heavy computation burden as well as uncertainty in choosing suitable weighting factors. This paper will formulate the compromise in competition mechanism in EMMS model as a noncooperative game with constraints, and will describe the desired stable system state as a generalized Nash equilibrium. Then the authors will investigate the game theoretical approach for two typical systems in chemical engineering, the gas-solid fluidiza- tion (GSF) system and turbulent flow in pipe. Two different cases for generalized Nash equilibrinm in such systems will be well defined and distinguished. The generalize Nash equilibrium will be solved accurately for the GSF system and a feasible method will be given for turbulent flow in pipe. These results coincide with the existing computational results and show the feasibility of this approach, which overcomes the disadvantages of the existing methods and provides deep insight into the mechanisms of multi-scale structure in the multi-phase complex systems in chemical engineering.展开更多
A method is presented to calculate the resistance of a high-speed displacement ship taking the effect of sinkage and trim and viscosity of fluid into account.A free surface flow field is evaluated by solving Reynolds ...A method is presented to calculate the resistance of a high-speed displacement ship taking the effect of sinkage and trim and viscosity of fluid into account.A free surface flow field is evaluated by solving Reynolds averaged Navier-Stokes(RANS) equations with volume of fluid(VoF) method.The sinkage and trim are computed by equating the vertical force and pitching moment to the hydrostatic restoring force and moment.The software Fluent,Maxsurf and MATLAB are used to implement this method.With dynamic mesh being used,the position of a ship is updated by the motion of "ship plus boundary layer" grid zone.The hull factors are introduced for fast calculating the running attitude of a ship.The method has been applied to the ship model INSEAN2340 for different Froude numbers and is found to be efficient for evaluating the flow field,resistance,sinkage and trim.展开更多
This paper studies the local receptivity of the Mack-mode instability to localized unsteady blowing and suction(UBS)in a chem-ical non-equilibrium(CNE)hypersonic boundary layer.The five-species CNE model is employed,a...This paper studies the local receptivity of the Mack-mode instability to localized unsteady blowing and suction(UBS)in a chem-ical non-equilibrium(CNE)hypersonic boundary layer.The five-species CNE model is employed,and the receptivity efficiency is formulated by use of the residual theorem.Compared with the results for the calorically perfect gas(CPG)model,we find that the real-gas effect enhances the receptivity efficiency remarkably in the majority of the second-mode frequency band,and the enhancement is mainly attributed to the modification of the base flow due to the CNE effect,which is akin to the cold-wall effect in hypersonic boundary layers.Combined with the destabilizing role of the CNE effect on the Mack second mode,it is concluded that the CNE effect would lead to a greater linearly accumulated perturbation amplitude,implying premature of transition to turbulence in a high-enthalpy hypersonic boundary layer subject to localized perturbations.展开更多
The Rotation and Curvature(RC)correction is an important turbulence model modifi-cation approach,and the Spalart-Allmaras model with the RC correction(SA-RC)has been exten-sively studied and used.As a multiplier of th...The Rotation and Curvature(RC)correction is an important turbulence model modifi-cation approach,and the Spalart-Allmaras model with the RC correction(SA-RC)has been exten-sively studied and used.As a multiplier of the modelling equation’s production term,the rotation function f_(r1)should have a cautiously designed value range,but its limit varies in different models and flow solvers.Therefore,the need of restriction is discussed theoretically,and the common range of f_(r1)is explored in Burgers vortexes.Afterwards,the SA-RC model with different limits is tested numerically.Negative f_(r1)always appears in the SA-RC model,and the difference between simula-tion results brought by the limits is not negligible.A lower limit of 0 enhances turbulence produc-tion,and therefore the vortex structures are dissipated faster and shrink in size,while an upper limit plays an opposite role.Considering that the lower limit of 0 usually promotes the simulation accu-racy and fixes the numerical defect,whereas the upper limit worsens the predictive performance in most cases,it is recommended to limit f_(r1)non-negative while utilizing the SA-RC model.In addi-tion,the RC-corrected model has a better prediction of the attached flow near curved walls,while the SA-Helicity model largely improves the simulation accuracy of three-dimensional large-scale vortices.The model combining both corrections has the potential to become more adaptive and more accurate.展开更多
The circumferentially averaged equation of the inlet flow radial equilibrium in axial compressor was deduced. It indicates that the blade inlet radial pressure gradient is closely related to the radial component of th...The circumferentially averaged equation of the inlet flow radial equilibrium in axial compressor was deduced. It indicates that the blade inlet radial pressure gradient is closely related to the radial component of the circumferential fluctuation(CF) source item. Several simplified cascades with/without aerodynamic loading were numerically studied to investigate the effects of blade bowing on the inlet flow radial equilibrium. A data reduction program was conducted to obtain the CF source from three-dimensional(3D) simulation results. Flow parameters at the passage inlet were focused on and each term in the radial equilibrium equation was discussed quantitatively. Results indicate that the inviscid blade force is the inducement of the inlet CF due to geometrical asymmetry. Blade bowing induces variation of the inlet CF, thus changes the radial pressure gradient and leads to flow migration before leading edge(LE) in the cascades. Positive bowing drives the inlet flow to migrate from end walls to mid-span and negative bowing turns it to the reverse direction to build a new equilibrium. In addition, comparative studies indicate that the inlet Mach number and blade loading can efficiently impact the effectiveness of blade bowing on radial equilibrium in compressor design.展开更多
Surface heterogeneous catalysis in a high-enthalpy dissociated environment leads to a remarkable enhancement of aerodynamic heating into the thermal protection surface of hypersonic aircraft.To more accurately predict...Surface heterogeneous catalysis in a high-enthalpy dissociated environment leads to a remarkable enhancement of aerodynamic heating into the thermal protection surface of hypersonic aircraft.To more accurately predict this catalytic heating,a kinetic catalytic model was constructed.This model involved four elementary reactions,the rates of which were determined on mean-field approximation and surface steady-state reaction assumption.By coupling this model into the viscous wall boundary condition of computational fluid dynamics(CFD)solver,the influences of metal material catalytic properties on heat and mass transfer into thermal protection materials were numerically investigated.Numerical results showed that atomic oxygen recombination catalyzed by surface material accounts for a major contribution to aerodynamic heating and thus variation in recombination rates from different materials leads to the significant difference in surface heat fluxes.From a comparative analysis of various materials,the catalytic activity increases from the inert platinum(Pt)to nickel(Ni)and finally to the active copper(Cu).As a result,the catalytic heating on Cu surface was more than twice of that on Pt surface.Further parametrical research revealed that the proper layout of inert material at the nose of aircraft could prevent stagnation catalytic heating from thermal damage by carrying near-wall dissociated atoms from the stagnation zone downstream.The material-relied heterogeneous catalysis mechanism in this study provides some technical support for the thermal protection system design of hypersonic aircraft.展开更多
基金support by the Young Elite Scientists Program of CSEE (No. JLB-2018-95)the National Natural Science Foundation of China (No. 51621065, No. U1766203)+1 种基金the support by FEDER funds through COMPETE 2020by Portuguese funds through FCT, under SAICT-PAC/0004/2015 (No. POCI-01-0145-FEDER-016434), 02/SAICT/2017 (No. POCI-01-0145-FEDER-029803) and UID/EEA/50014/2019 (No. POCI-01-0145-FEDER-006961)
文摘The rapidly increasing penetration of electric vehicles(EVs) in modern metropolises has been witnessed during the past decade, inspired by financial subsidies as well as public awareness of climate change and environment pro-tection. Integrating charging facilities, especially highpower chargers in fast charging stations, into power distribution systems remarkably alters the traditional load flow pattern, and thus imposes great challenges on the operation of distribution network in which controllable resources are rare. On the other hand, provided with appropriate incentives, the energy storage capability of electric vehicle offers a unique opportunity to facilitate the integration of distributed wind and solar power generation into power distribution system. The above trends call for thorough investigation and research on the interdependence between transportation system and power distribution system. This paper conducts a comprehensive survey on this line of research. The basic models of transportation system and power distribution system are introduced,especially the user equilibrium model, which describes the vehicular flow on each road segment and is not familiar to the readers in power system community. The modelling of interdependence across the two systems is highlighted.Taking into account such interdependence, applications ranging from long-term planning to short-term operation are reviewed with emphasis on comparing the description of traffic-power interdependence. Finally, an outlook of prospective directions and key technologies in future research is summarized.
基金Project(2015CB060200)supported by the National Basic Research and Development Program of ChinaProjects(51322403,51274254)supported by the National Natural Science Foundation of China
文摘In order to investigate the micro-process and inner mechanism of rock failure under impact loading, the laboratory tests were carried out on an improved split Hopkinson pressure bar (SHPB) system with synchronized measurement devices including a high-speed camera and a dynamic strain meter. The experimental results show that the specimens were in the state of good stress equilibrium during the post failure stage even when visible cracks were forming in the specimens. Rock specimens broke into strips but still could bear the external stress and keep force balance. Meanwhile, numerical tests with particle flow code (PFC) revealed that the failure process of rocks can be described by the evolution of micro-fractures. Shear cracks emerged firstly and stopped developing when the external stress was not high enough. Tensile cracks, however, emerged when the rock specimen reached its peak strength and played an important role in controlling the ultimate failure during the post failure stage.
基金the National Natural Science Foundation of China (Grand No. 10672024)the Scientific Research Foundation for the Returned Overseas Chinese Scholars, Ministry of Education of China
文摘As one of the most important components of river mechanics,sediment transport capacity of sediment-laden flows has attracted much attention from many researchers working on river mechanics and hydraulic engineering. Based on the time-averaged equation for a turbulent energy equilibrium in solid and liquid two-phase flow,an expression for the efficiency coefficient of suspended load movement was derived for the two-dimensional,steady,uniform,fully-developed turbulent flow. A new structural expression of sediment transport capacity was achieved. Using 115 runs of flume experimental data,which were obtained through two kinds of sediment transport experiments in the state of equilibrium,in combination with the basic rheological and sediment transporting characteristics of hyperconcentrated flow,the main parameters in the structural expression of sediment transport capacity were calibrated,and a new formula of sediment transport capacity for hyperconcentrated flow was developed. A large amount of field data from the Yellow River,Wuding River,and Yangtze River,etc. were adopted to verify the new formula and good agreement was obtained. These results above contribute to an improved theoretical system of river mechanics and a reliable tool for management of rivers carrying high concentration of sediments.
基金supported by Hong Kong Research Grant Council 621709 and 621011,National Natural Science Foundation of China(Project No.10928205)National Key Basic Research Program(2009CB724101)J.C.Huang was supported by National Science Council of Taiwan through grant no.NSC 100-2221-E-019-048-MY3。
文摘With discretized particle velocity space,a multi-scale unified gas-kinetic scheme for entire Knudsen number flows has been constructed based on the kinetic model in one-dimensional case[J.Comput.Phys.,vol.229(2010),pp.7747-7764].For the kinetic equation,to extend a one-dimensional scheme to multidimensional flow is not so straightforward.The major factor is that addition of one dimension in physical space causes the distribution function to become two-dimensional,rather than axially symmetric,in velocity space.In this paper,a unified gas-kinetic scheme based on the Shakhov model in two-dimensional space will be presented.Instead of particle-based modeling for the rarefied flow,such as the direct simulation Monte Carlo(DSMC)method,the philosophical principal underlying the current study is a partial-differential-equation(PDE)-based modeling.Since the valid scale of the kinetic equation and the scale of mesh size and time step may be significantly different,the gas evolution in a discretized space is modeled with the help of kinetic equation,instead of directly solving the partial differential equation.Due to the use of both hydrodynamic and kinetic scales flow physics in a gas evolution model at the cell interface,the unified scheme can basically present accurate solution in all flow regimes from the free molecule to the Navier-Stokes solutions.In comparison with the DSMC and Navier-Stokes flow solvers,the current method is much more efficient than DSMC in low speed transition and continuum flow regimes,and it has better capability than NS solver in capturing of non-equilibrium flow physics in the transition and rarefied flow regimes.As a result,the current method can be useful in the flow simulation where both continuum and rarefied flow physics needs to be resolved in a single computation.This paper will extensively evaluate the performance of the unified scheme fromfreemolecule to continuum NS solutions,and fromlow speedmicro-flow to high speed non-equilibrium aerodynamics.The test cases clearly demonstrate tha
基金supported by the National Natural Science Foundation of China under Grant Nos.11688101,91634203,61304159by the National Center for Mathematics and Interdisciplinary Sciences
文摘This paper explores the application of noncooperative game theory together with the concept of Nash equilibrium to the investigation of some basic problems on multi-scale structure, especially the meso-scale structure in the multi-phase complex systems in chemical engineering. The basis of this work is the energy-minimization-multi-scale (EMMS) model proposed by Li and Kwauk (1994) and Li, et al. (2013) which identifies the multi-scale structure as a result of 'compromise-in-competition between dominant mechanisms' and tries to solve a multi-objective optimization problem. However, the existing methods often integrate it into a problem of single objective optimization, which does not clearly reflect the 'compromise-in-competition' mechanism and causes heavy computation burden as well as uncertainty in choosing suitable weighting factors. This paper will formulate the compromise in competition mechanism in EMMS model as a noncooperative game with constraints, and will describe the desired stable system state as a generalized Nash equilibrium. Then the authors will investigate the game theoretical approach for two typical systems in chemical engineering, the gas-solid fluidiza- tion (GSF) system and turbulent flow in pipe. Two different cases for generalized Nash equilibrinm in such systems will be well defined and distinguished. The generalize Nash equilibrium will be solved accurately for the GSF system and a feasible method will be given for turbulent flow in pipe. These results coincide with the existing computational results and show the feasibility of this approach, which overcomes the disadvantages of the existing methods and provides deep insight into the mechanisms of multi-scale structure in the multi-phase complex systems in chemical engineering.
基金the National Natural Science Foundation of China (No.50879090)the Advanced Research Program of GAD of the P.L.A (No.7131005)
文摘A method is presented to calculate the resistance of a high-speed displacement ship taking the effect of sinkage and trim and viscosity of fluid into account.A free surface flow field is evaluated by solving Reynolds averaged Navier-Stokes(RANS) equations with volume of fluid(VoF) method.The sinkage and trim are computed by equating the vertical force and pitching moment to the hydrostatic restoring force and moment.The software Fluent,Maxsurf and MATLAB are used to implement this method.With dynamic mesh being used,the position of a ship is updated by the motion of "ship plus boundary layer" grid zone.The hull factors are introduced for fast calculating the running attitude of a ship.The method has been applied to the ship model INSEAN2340 for different Froude numbers and is found to be efficient for evaluating the flow field,resistance,sinkage and trim.
基金supported by the National Science Foundation of China(Grant Nos.U20B2003,11988102,12302292).
文摘This paper studies the local receptivity of the Mack-mode instability to localized unsteady blowing and suction(UBS)in a chem-ical non-equilibrium(CNE)hypersonic boundary layer.The five-species CNE model is employed,and the receptivity efficiency is formulated by use of the residual theorem.Compared with the results for the calorically perfect gas(CPG)model,we find that the real-gas effect enhances the receptivity efficiency remarkably in the majority of the second-mode frequency band,and the enhancement is mainly attributed to the modification of the base flow due to the CNE effect,which is akin to the cold-wall effect in hypersonic boundary layers.Combined with the destabilizing role of the CNE effect on the Mack second mode,it is concluded that the CNE effect would lead to a greater linearly accumulated perturbation amplitude,implying premature of transition to turbulence in a high-enthalpy hypersonic boundary layer subject to localized perturbations.
基金supported by the National Natural Science Foundation of China(Nos.51976006,51790513)the Aeronautical Science Foundation of China(No.2018ZB51013)+1 种基金the National Science and Technology Major Project,China(2017-II-003-0015)the Open Fund from State Key Laboratory of Aerodynamics,China(No.SKLA2019A0101).
文摘The Rotation and Curvature(RC)correction is an important turbulence model modifi-cation approach,and the Spalart-Allmaras model with the RC correction(SA-RC)has been exten-sively studied and used.As a multiplier of the modelling equation’s production term,the rotation function f_(r1)should have a cautiously designed value range,but its limit varies in different models and flow solvers.Therefore,the need of restriction is discussed theoretically,and the common range of f_(r1)is explored in Burgers vortexes.Afterwards,the SA-RC model with different limits is tested numerically.Negative f_(r1)always appears in the SA-RC model,and the difference between simula-tion results brought by the limits is not negligible.A lower limit of 0 enhances turbulence produc-tion,and therefore the vortex structures are dissipated faster and shrink in size,while an upper limit plays an opposite role.Considering that the lower limit of 0 usually promotes the simulation accu-racy and fixes the numerical defect,whereas the upper limit worsens the predictive performance in most cases,it is recommended to limit f_(r1)non-negative while utilizing the SA-RC model.In addi-tion,the RC-corrected model has a better prediction of the attached flow near curved walls,while the SA-Helicity model largely improves the simulation accuracy of three-dimensional large-scale vortices.The model combining both corrections has the potential to become more adaptive and more accurate.
基金supported by the National Natural Science Foundation of China (Nos.51236001,51006005)the National Basic Research Program of China (No. 2012CB720201)Beijing Natural Science Foundation (No. 3151002)
文摘The circumferentially averaged equation of the inlet flow radial equilibrium in axial compressor was deduced. It indicates that the blade inlet radial pressure gradient is closely related to the radial component of the circumferential fluctuation(CF) source item. Several simplified cascades with/without aerodynamic loading were numerically studied to investigate the effects of blade bowing on the inlet flow radial equilibrium. A data reduction program was conducted to obtain the CF source from three-dimensional(3D) simulation results. Flow parameters at the passage inlet were focused on and each term in the radial equilibrium equation was discussed quantitatively. Results indicate that the inviscid blade force is the inducement of the inlet CF due to geometrical asymmetry. Blade bowing induces variation of the inlet CF, thus changes the radial pressure gradient and leads to flow migration before leading edge(LE) in the cascades. Positive bowing drives the inlet flow to migrate from end walls to mid-span and negative bowing turns it to the reverse direction to build a new equilibrium. In addition, comparative studies indicate that the inlet Mach number and blade loading can efficiently impact the effectiveness of blade bowing on radial equilibrium in compressor design.
基金financial support of the National Key Research and Development Plan of China through the project(No.2019YFA0405202)National Natural Science Foundation of China through the project(No.12072361)。
文摘Surface heterogeneous catalysis in a high-enthalpy dissociated environment leads to a remarkable enhancement of aerodynamic heating into the thermal protection surface of hypersonic aircraft.To more accurately predict this catalytic heating,a kinetic catalytic model was constructed.This model involved four elementary reactions,the rates of which were determined on mean-field approximation and surface steady-state reaction assumption.By coupling this model into the viscous wall boundary condition of computational fluid dynamics(CFD)solver,the influences of metal material catalytic properties on heat and mass transfer into thermal protection materials were numerically investigated.Numerical results showed that atomic oxygen recombination catalyzed by surface material accounts for a major contribution to aerodynamic heating and thus variation in recombination rates from different materials leads to the significant difference in surface heat fluxes.From a comparative analysis of various materials,the catalytic activity increases from the inert platinum(Pt)to nickel(Ni)and finally to the active copper(Cu).As a result,the catalytic heating on Cu surface was more than twice of that on Pt surface.Further parametrical research revealed that the proper layout of inert material at the nose of aircraft could prevent stagnation catalytic heating from thermal damage by carrying near-wall dissociated atoms from the stagnation zone downstream.The material-relied heterogeneous catalysis mechanism in this study provides some technical support for the thermal protection system design of hypersonic aircraft.