Temporal mode direction numerical simulation has been done for the transition of a supersonic boundary layer on a flat plate with Mach number 4.5. Analysis of the result showed that during the breakdown process in lam...Temporal mode direction numerical simulation has been done for the transition of a supersonic boundary layer on a flat plate with Mach number 4.5. Analysis of the result showed that during the breakdown process in laminar-turbulent transition, the mechanism causing the mean flow profile to evolve swiftly from laminar to turbulent was that the modification of mean flow profile by the disturbances,when they became larger,led to a remarkable change of its stability characteristics. Though the most unstable T-S wave was of second mode for laminar flow, the first mode waves played the key role in the breakdown process in laminar-turbulent transition.展开更多
Based on Reynolds-averaged Navier-Stokes approach,a laminar-turbulence transition model is proposed in this study that takes into account the effects of different instability modes associated with the variations in Ma...Based on Reynolds-averaged Navier-Stokes approach,a laminar-turbulence transition model is proposed in this study that takes into account the effects of different instability modes associated with the variations in Mach numbers of compressible boundary layer flows.The model is based on k-ω-γ three-equation eddy-viscosity concept with k representing the fluctuating kinetic energy,ωthe specific dissipation rate and the intermittency factorγ.The particular features of the model are that:1)k includes the non-turbulent,as well as turbulent fluctuations;2)a transport equation for the intermittency factorγis proposed here with a source term set to trigger the transition onset;3)through the introduction of a new length scale normal to wall,the present model employs the local variables only avoiding the use of the integral parameters,like the boundary layer thicknessδ,which are often cost-ineffective with the modern CFD(Computational Fluid Dynamics)methods;4)in the fully turbulent region,the model retreats to the well-known k-ωSST(Shear Stress Transport)model.This model is validated with a number of available experiments on boundary layer transitions including the incompressible,supersonic and hypersonic flows past flat plates,straight/flared cones at zero incidences,etc.It is demonstrated that the present model can be successfully applied to the engineering calculations of a variety of aerodynamic flow transition.展开更多
The aero-heating of the rudder shaft region of a hypersonic vehicle is very harsh, as the peak heat flux in this region can be even higher than that at the stagnation point. Therefore, studying the aero-heating of the...The aero-heating of the rudder shaft region of a hypersonic vehicle is very harsh, as the peak heat flux in this region can be even higher than that at the stagnation point. Therefore, studying the aero-heating of the rudder shaft is of great significance for designing the thermal protection system of the hypersonic vehicle. In the wind tunnel test of the aero-heating effect, we find that with the increase of the angle of attack of the lifting body model, the increasement of the heat flux of the rudder shaft is larger under laminar flow conditions than that under turbulent flow conditions. To understand this, we design a wind tunnel experiment to study the effect of laminar/turbulent hypersonic boundary layers on the heat flux of the rudder shaft under the same wind tunnel freestream conditions. The experiment is carried out in the ?2 m shock tunnel(FD-14 A) affiliated to the China Aerodynamics Research and Development Center(CARDC). The laminar boundary layer on the model is triggered to a turbulent one by using vortex generators, which are 2 mm-high diamonds. The aero-heating of the rudder shaft(with the rudder) and the protuberance(without the rudder) are studied in both hypersonic laminar and turbulent boundary layers under the same freestream condition. The nominal Mach numbers are 10 and 12, and the unit Reynolds numbers are2.4 × 10~6 m^(-1) and 2.1 × 10~6 m-1. The angle of attack of the model is 20°, and the deflection angle of the rudder and the protuberance is 10°. The heat flux on the model surface is measured by thin film heat flux sensors, and the heat flux distribution along the center line of the lifting body model suggests that forced transition is achieved in the upstream of the rudder. The test results of the rudder shaft and the protuberance show that the heat flux of the rudder shaft is lower in the turbulent flow than that in the laminar flow, but the heat flux of the protuberance is the other way around,i.e., lower in the laminar flow than in the turbulent flow. The wind tunnel test resul展开更多
This paper describes a simplified transition model based on the recently developed correlation-based γ - Reot transition model. The transport equation of transition momentum thick- ness Reynolds number is eliminated ...This paper describes a simplified transition model based on the recently developed correlation-based γ - Reot transition model. The transport equation of transition momentum thick- ness Reynolds number is eliminated for simplicity, and new transition length function and critical Reynolds number correlation are proposed. The new model is implemented into an in-house com- putational fluid dynamics (CFD) code and validated for low and high-speed flow cases, including the zero pressure flat plate, airfoils, hypersonic flat plate and double wedge. Comparisons between the simulation results and experimental data show that the boundary-layer transition phenomena can be reasonably illustrated by the new model, which gives rise to significant improvements over the fully laminar and fully turbulent results. Moreover, the new model has comparable features of accuracy and applicability when compared with the original 3' - Reot model. In the meantime, the newly proposed model takes only one transport equation of intermittency factor and requires fewer correlations, which simplifies the original model greatly. Further studies, especially on separation- induced transition flows, are required for the improvement of the new model.展开更多
Transition prediction for boundary layers has always been one of the urgent problems waiting for a solution for the development of aero-space technology,yet there is no reliable and effective method due to the complex...Transition prediction for boundary layers has always been one of the urgent problems waiting for a solution for the development of aero-space technology,yet there is no reliable and effective method due to the complexity of the problem.The eN method has been regarded as an effective method for the transition prediction of boundary layers.However,it heavily relies on experiment or experience.And in cases with three-dimensional base flow,for instance,the boundary layer on a cone with angle of attack,the result of its application is not satisfactory.The authors have found its cause and proposed the method for its improvement,which did yield the fairly satisfactory result for a given test case,and also did not rely so much on experiment or experience.However,before people can really apply this method to practical problems,more test cases have to be studied.In this paper,more test cases for the application of the improved eN method to problems of transition prediction of supersonic and hypersonic boundary layers on cones with angle of attack will be studied.The results are compared with those obtained by experiments and/or direct numerical simulations,confirming that the improved eN method is effective and reliable.We also find that there may be more than one ZARF for each meridian plane,and which one should be chosen for the eN method has been clarified.展开更多
The effects induced by the presence of incoming wakes on the boundary layer developing over a high-lift low-pressure turbine profile have been investigated in a linear cascade at mid-span.The tested Reynolds number is...The effects induced by the presence of incoming wakes on the boundary layer developing over a high-lift low-pressure turbine profile have been investigated in a linear cascade at mid-span.The tested Reynolds number is 70000,typical of the cruise operating condition.The results of the investigations performed under steady and unsteady inflow conditions are analyzed.The unsteady investigations have been performed at the reduced frequency of f+=0.62,representative of the real engine operating condition.Profile aerodynamic loadings as well as boundary layer velocity profiles have been measured to survey the separation and transition processes.Spectral analysis has been also performed to better understand the phenomena associated with the transition process under steady inflow.For the unsteady case,a phase-locked ensemble averaging technique has been employed to reconstruct the time-resolved boundary layer velocity distributions from the hot-wire instantaneous signal output.The ensemble-averaging technique allowed a detailed analysis of the effects induced by incoming wakes-boundary layer interaction in separation suppression.Time-resolved results are presented in terms of mean velocity and unresolved unsteadiness time-space plots.展开更多
As boundary layer transition plays an important role in aerodynamic drag prediction, the proposal and study of transition prediction methods simulating the complex flow phenomena are prerequisite for aerodynamic desig...As boundary layer transition plays an important role in aerodynamic drag prediction, the proposal and study of transition prediction methods simulating the complex flow phenomena are prerequisite for aerodynamic design. In this paper, with the application of the linear stability theory based on amplification factor transport transition equations on the two-equation shear stress transport (SST) eddy-viscosity model, a new method, the SST-NTS-NCF model, is yielded. The new amplification factor transport equation for the crossflow instability induced transition is proposed to add to the NTS equation proposed by Coder, which simulates Tollmien-Schlichting wave transition. The turbulent kinetic energy equation is modified by introducing a new source term that simulates the transition process without the intermittency factor equation. Finally, coupled with these two amplification factor transport equations and SST turbulence model, a four-equation transition turbulence model is built. Comparisons between predictions using the new model and wind-tunnel experiments of NACA64(2) A015, NLF(2)-0415 and ONERA-D infinite swept wing and ONERAM6 swept wing validate the predictive quality of the new SST-NTS-NCF model. (C) 2016 Chinese Society of Aeronautics and Astronautics. Production and hosting by Elsevier Ltd.展开更多
The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The ...The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach numberMc = 0.8 andRe = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of A -vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.展开更多
In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach num...In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach number 5. For single roughness, the DNS results showed that both horseshoe vortices and hairpin vortices caused by shear layer instability can affect the boundary layer instability. The generation of the near-wall unstable structure is the key point of boundary layer transition behind the roughness element. At the downstream of the roughness element, the interaction between horseshoe vortices and hairpin vortices will spread in the spanwise direction.For double roughness elements, the effect of the spacing between roughness elements on the transition is studied. It is found that the case of higher spacing between roughness elements is more effective for inducing transition than the lower one. The interaction between two adjacent roughness elements can suppress the evolution of horseshoe vortices in the downstream and trigger the instability of shear layer. Thus, the transition will be suppressed accordingly.展开更多
Our previous study showed that the frictional drag decreases with increasing void fraction at Re〉1300, while it increases at Re 〈 1000. Decomposition of the Reynolds shear stress also implied that bubbles induce iso...Our previous study showed that the frictional drag decreases with increasing void fraction at Re〉1300, while it increases at Re 〈 1000. Decomposition of the Reynolds shear stress also implied that bubbles induce isotropy of turbulence. In order to confirm our previous analysis and to further investigate flow fields in the vicinity of bubbles, we analyze velocity fluctuations on the quadrant space in the streamwise and transverse directions (u′-v′ plane). Here, we focus on two specific Reynolds numbers (at Re≈900 and ≈1410, which are close to the laminar-to-turbulent transition regime) and discuss bubble effects on sweep (u′〉 0, v′〈 0 ) and ejection (u′〈 0, v′〉 0) events as a function of the Reynolds number. We also illustrate velocity fluctuations in the vicinity of an individual bubble and a swarm of bubbles on the u′- v′ coordinates. The results show that a bubble swarm suppresses the velocity fluctuations at Re≈1410.展开更多
Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20° cone angle (or 10° half-cone angle) is performed by using the mixed seventh- order up-wind biased finite difference sche...Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20° cone angle (or 10° half-cone angle) is performed by using the mixed seventh- order up-wind biased finite difference scheme and sixth-order central difference scheme. The free stream Mach number is 0.7 and free stream unit Reynolds number is 250000/inch. The characteristics of transition and turbulence of the sharp cone boundary layer are compared with those of the flat plate boundary layer. Statistics of fully developed turbulent flow agree well with the experimental and theoretical data for the turbulent flat-plate boundary layer flow. The near wall streak-like structure is shown and the average space between streaks (normalized by the local wall unit) keeps approximately invariable at different streamwise locations. The turbulent energy equation in the cylindrical coordinate is given and turbulent en-ergy budget is studied. The computed results show that the effect of circumferen-tial curvature on turbulence characteristics is not obvious.展开更多
The achievement of laminar flow in the boundary layer at high-speed cruise conditions may further, in addition to shock-wave control, reduce the drag and extend the range of military fighter aircraft. To this end, a f...The achievement of laminar flow in the boundary layer at high-speed cruise conditions may further, in addition to shock-wave control, reduce the drag and extend the range of military fighter aircraft. To this end, a further investigation on transitional boundary-layer flow of fighter wings is needed due to different configurations from the wings used on conventional transport aircraft. In this paper, wind tunnel experiments and numerical simulations were conducted on three-dimensional transition of thin diamond-shaped wings used on advanced fighter aircraft at tran/supersonic design points. A newly proposed correlation of crossflow transition which includes the effect of surface roughness was introduced into the c-Rehttransition model. Predicted results were in good agreement with flow visualizations. Results showed that the strength of the crossflow component grew rapidly around the leading edge because of the severe flow acceleration, just as the same as wings with a large aspect ratio. However, there seemed no regular pattern of instabilitydominance variation in span-wise for a diamond configuration. The dominance of different instability mechanisms strongly depended on the local pressure distribution. Hereby, the research recommended a ‘‘roof-like" shape of pressure distribution to suppress both crossflow and Tollmien-Schlichting(T-S) instabilities. Besides, a sharp suction peak with a serious pressure rise should be cut off to avoid stronger instabilities. Further discussions also revealed an independence of the unit Reynolds number when transition was triggered by T-S instabilities. Aerodynamic force comparisons indicated that further benefit on drag reduction could be expected by including the three-dimensional transition effect into a wing design process.展开更多
This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transp...This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transport equations, separately, to compute laminar kinetic energy, turbulent kinetic energy, and dissipation rate in a flow field. It needs neither correlations of intermittency factors nor knowledge of the transition onset. Two transition tests are carried out: flat plate boundary layer under zero ...展开更多
At the late stage of transitional boundary layers, the nonlinear evolution of the ring-like vortices and spike structures and their effects on the surrounding flow were studied by means of direct numerical simulation ...At the late stage of transitional boundary layers, the nonlinear evolution of the ring-like vortices and spike structures and their effects on the surrounding flow were studied by means of direct numerical simulation with high order accuracy. A spatial transition of the flat-plate boundary layers in the compressible flow was conducted. Detailed numerical results with high resolution clearly represented the typical vortex structures, such as ring-like vortices and so on, and induced ejection and sweep events. It was verified that the formation of spike structures in transitional boundary layers had close relationship with ring-like vortices. Especially, compared to the newly observed positive spike structure in the experiments, the same structure was found in the present numerical simulations, and the mechanism was also studied and analyzed.展开更多
Spatial mode direct numerical simulation has been applied to study the mechanism of breakdown in laminar-turbulent transition of a supersonic boundary layer on a fiat plate with Mach number 4.5. Analysis of the result...Spatial mode direct numerical simulation has been applied to study the mechanism of breakdown in laminar-turbulent transition of a supersonic boundary layer on a fiat plate with Mach number 4.5. Analysis of the result showed that, during the breakdown process in laminar-turbulent transition, the mechanism causing the mean flow profile to evolve swiftly from laminar to turbulent was that the modification of mean flow profile by the disturbance, when they became larger, leads to remarkable change of its stability characteristics. Though the most unstable T-S wave was of second mode for laminar flow, the first mode waves played the key role in the breakdown process in laminar-turbulent transition.展开更多
The thermal barrier coatings with NiCrAlY alloy bonding layer, NiCrAlY Y 2O 3 stabilized ZrO 2 transition layer and Y 2O 3 stabilized ZrO 2 ceramic layer are prepared on nickel alloy substrates using the plasma spray ...The thermal barrier coatings with NiCrAlY alloy bonding layer, NiCrAlY Y 2O 3 stabilized ZrO 2 transition layer and Y 2O 3 stabilized ZrO 2 ceramic layer are prepared on nickel alloy substrates using the plasma spray technique. The relationship among the composition, structure and property of the coatings are investiga ted by means of optical microscope, scanning electronic microscope and the experiments of thermal shock resistance cycling and high temperature oxidation resistance. The results show that the structure design of introdu cing a transition layer between Ni alloy substrate and ZrO 2 ceramic coating guarantees the high quality and properties of the coatings; ZrO 2 coatings doped with a little SiO 2 possesses better thermal shock resistance and more excellent hot corrosion resistance as compared with ZrO 2 coating materials without SiO 2 ;the improvement in performance of ZrO 2 coating doped with SiO 2 is due to forming more dense coating structure by self closing effects of the flaws and pores in the ZrO 2 coatings.展开更多
Experimental study was conducted for boundarylayers on a sharp 5° half-angle cone of 400mm length at angles of attack. The model was tested in the T-326 hypersonic wind tunnel (ITAM) at freestream Mach number M...Experimental study was conducted for boundarylayers on a sharp 5° half-angle cone of 400mm length at angles of attack. The model was tested in the T-326 hypersonic wind tunnel (ITAM) at freestream Mach number M = 5.95. Mean and fluctuation wall characteristics of the boundary layer are measured at 0°, 2°, 3° and 4° angles of attack for different stagnation pressures. Pulsation measurements are carried out by means of ALTP sensor. Pressure and temperature distributions along the model are obtained, and transition beginning and end locations have been found. Boundary layer stabilization with the increase of angle of attack and the decrease of stagnation pressure is observed. High frequency pulsations inherent to hypersonic boundary layer (second mode) have been detected.展开更多
Direct numerical simulation of spatially evolving compressible boundary layer over a blunt wedge is performed in this paper. The free-stream Mach number is 6 and the disturbance source produced by wall blowing and suc...Direct numerical simulation of spatially evolving compressible boundary layer over a blunt wedge is performed in this paper. The free-stream Mach number is 6 and the disturbance source produced by wall blowing and suction is located downstream of the sound-speed point. Statistics are studied and compared with the results in incompressible flat-plate boundary layer. The mean pressure gradient effects on the vortex structure are studied.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.90205021)the Foundation for the Author of National Excellent Doctoral Dissertation(FANEDD)of China(Grant No.200328)partially supported by Liu-Hui Center of Applied Mathematics,Nankai University and Tianjin University.
文摘Temporal mode direction numerical simulation has been done for the transition of a supersonic boundary layer on a flat plate with Mach number 4.5. Analysis of the result showed that during the breakdown process in laminar-turbulent transition, the mechanism causing the mean flow profile to evolve swiftly from laminar to turbulent was that the modification of mean flow profile by the disturbances,when they became larger,led to a remarkable change of its stability characteristics. Though the most unstable T-S wave was of second mode for laminar flow, the first mode waves played the key role in the breakdown process in laminar-turbulent transition.
基金Supported by the National Natural Science Foundation of China(Grant Nos.10232020 and 90505005)
文摘Based on Reynolds-averaged Navier-Stokes approach,a laminar-turbulence transition model is proposed in this study that takes into account the effects of different instability modes associated with the variations in Mach numbers of compressible boundary layer flows.The model is based on k-ω-γ three-equation eddy-viscosity concept with k representing the fluctuating kinetic energy,ωthe specific dissipation rate and the intermittency factorγ.The particular features of the model are that:1)k includes the non-turbulent,as well as turbulent fluctuations;2)a transport equation for the intermittency factorγis proposed here with a source term set to trigger the transition onset;3)through the introduction of a new length scale normal to wall,the present model employs the local variables only avoiding the use of the integral parameters,like the boundary layer thicknessδ,which are often cost-ineffective with the modern CFD(Computational Fluid Dynamics)methods;4)in the fully turbulent region,the model retreats to the well-known k-ωSST(Shear Stress Transport)model.This model is validated with a number of available experiments on boundary layer transitions including the incompressible,supersonic and hypersonic flows past flat plates,straight/flared cones at zero incidences,etc.It is demonstrated that the present model can be successfully applied to the engineering calculations of a variety of aerodynamic flow transition.
基金supported by the National Key Research and Development Program of China(No.2016YFA0401201)
文摘The aero-heating of the rudder shaft region of a hypersonic vehicle is very harsh, as the peak heat flux in this region can be even higher than that at the stagnation point. Therefore, studying the aero-heating of the rudder shaft is of great significance for designing the thermal protection system of the hypersonic vehicle. In the wind tunnel test of the aero-heating effect, we find that with the increase of the angle of attack of the lifting body model, the increasement of the heat flux of the rudder shaft is larger under laminar flow conditions than that under turbulent flow conditions. To understand this, we design a wind tunnel experiment to study the effect of laminar/turbulent hypersonic boundary layers on the heat flux of the rudder shaft under the same wind tunnel freestream conditions. The experiment is carried out in the ?2 m shock tunnel(FD-14 A) affiliated to the China Aerodynamics Research and Development Center(CARDC). The laminar boundary layer on the model is triggered to a turbulent one by using vortex generators, which are 2 mm-high diamonds. The aero-heating of the rudder shaft(with the rudder) and the protuberance(without the rudder) are studied in both hypersonic laminar and turbulent boundary layers under the same freestream condition. The nominal Mach numbers are 10 and 12, and the unit Reynolds numbers are2.4 × 10~6 m^(-1) and 2.1 × 10~6 m-1. The angle of attack of the model is 20°, and the deflection angle of the rudder and the protuberance is 10°. The heat flux on the model surface is measured by thin film heat flux sensors, and the heat flux distribution along the center line of the lifting body model suggests that forced transition is achieved in the upstream of the rudder. The test results of the rudder shaft and the protuberance show that the heat flux of the rudder shaft is lower in the turbulent flow than that in the laminar flow, but the heat flux of the protuberance is the other way around,i.e., lower in the laminar flow than in the turbulent flow. The wind tunnel test resul
基金supported by the State Key Development Program for Basic Research of China(No.2014CB340201)
文摘This paper describes a simplified transition model based on the recently developed correlation-based γ - Reot transition model. The transport equation of transition momentum thick- ness Reynolds number is eliminated for simplicity, and new transition length function and critical Reynolds number correlation are proposed. The new model is implemented into an in-house com- putational fluid dynamics (CFD) code and validated for low and high-speed flow cases, including the zero pressure flat plate, airfoils, hypersonic flat plate and double wedge. Comparisons between the simulation results and experimental data show that the boundary-layer transition phenomena can be reasonably illustrated by the new model, which gives rise to significant improvements over the fully laminar and fully turbulent results. Moreover, the new model has comparable features of accuracy and applicability when compared with the original 3' - Reot model. In the meantime, the newly proposed model takes only one transport equation of intermittency factor and requires fewer correlations, which simplifies the original model greatly. Further studies, especially on separation- induced transition flows, are required for the improvement of the new model.
基金Supported by the National Natural Science Foundation of China (Grant Nos. 10632050 and 90716007)the Special Foundation for the Authors of National Excel-lent Doctoral Dissertations (Grant No. 200328)the Foundation of Liu-Hui Center of Applied Mathematics of Nankai University and Tianjin University Contributed by ZHOU Heng
文摘Transition prediction for boundary layers has always been one of the urgent problems waiting for a solution for the development of aero-space technology,yet there is no reliable and effective method due to the complexity of the problem.The eN method has been regarded as an effective method for the transition prediction of boundary layers.However,it heavily relies on experiment or experience.And in cases with three-dimensional base flow,for instance,the boundary layer on a cone with angle of attack,the result of its application is not satisfactory.The authors have found its cause and proposed the method for its improvement,which did yield the fairly satisfactory result for a given test case,and also did not rely so much on experiment or experience.However,before people can really apply this method to practical problems,more test cases have to be studied.In this paper,more test cases for the application of the improved eN method to problems of transition prediction of supersonic and hypersonic boundary layers on cones with angle of attack will be studied.The results are compared with those obtained by experiments and/or direct numerical simulations,confirming that the improved eN method is effective and reliable.We also find that there may be more than one ZARF for each meridian plane,and which one should be chosen for the eN method has been clarified.
文摘The effects induced by the presence of incoming wakes on the boundary layer developing over a high-lift low-pressure turbine profile have been investigated in a linear cascade at mid-span.The tested Reynolds number is 70000,typical of the cruise operating condition.The results of the investigations performed under steady and unsteady inflow conditions are analyzed.The unsteady investigations have been performed at the reduced frequency of f+=0.62,representative of the real engine operating condition.Profile aerodynamic loadings as well as boundary layer velocity profiles have been measured to survey the separation and transition processes.Spectral analysis has been also performed to better understand the phenomena associated with the transition process under steady inflow.For the unsteady case,a phase-locked ensemble averaging technique has been employed to reconstruct the time-resolved boundary layer velocity distributions from the hot-wire instantaneous signal output.The ensemble-averaging technique allowed a detailed analysis of the effects induced by incoming wakes-boundary layer interaction in separation suppression.Time-resolved results are presented in terms of mean velocity and unresolved unsteadiness time-space plots.
基金supported by the National Basic Research Program of China (No. 2014CB744804)
文摘As boundary layer transition plays an important role in aerodynamic drag prediction, the proposal and study of transition prediction methods simulating the complex flow phenomena are prerequisite for aerodynamic design. In this paper, with the application of the linear stability theory based on amplification factor transport transition equations on the two-equation shear stress transport (SST) eddy-viscosity model, a new method, the SST-NTS-NCF model, is yielded. The new amplification factor transport equation for the crossflow instability induced transition is proposed to add to the NTS equation proposed by Coder, which simulates Tollmien-Schlichting wave transition. The turbulent kinetic energy equation is modified by introducing a new source term that simulates the transition process without the intermittency factor equation. Finally, coupled with these two amplification factor transport equations and SST turbulence model, a four-equation transition turbulence model is built. Comparisons between predictions using the new model and wind-tunnel experiments of NACA64(2) A015, NLF(2)-0415 and ONERA-D infinite swept wing and ONERAM6 swept wing validate the predictive quality of the new SST-NTS-NCF model. (C) 2016 Chinese Society of Aeronautics and Astronautics. Production and hosting by Elsevier Ltd.
文摘The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach numberMc = 0.8 andRe = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of A -vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.
基金the China Scholarship Council(CSC)the Aeronautics Science Foundation of China(No.20163252037)+2 种基金the China Postdoctoral Science Foundation(No.2017M610325)the Natural Science Foundation of Jiangsu Province of China(No.BK20170771)the Fundamental Research Funds for the Central Universities of China(No.NP2017202)for their support
文摘In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach number 5. For single roughness, the DNS results showed that both horseshoe vortices and hairpin vortices caused by shear layer instability can affect the boundary layer instability. The generation of the near-wall unstable structure is the key point of boundary layer transition behind the roughness element. At the downstream of the roughness element, the interaction between horseshoe vortices and hairpin vortices will spread in the spanwise direction.For double roughness elements, the effect of the spacing between roughness elements on the transition is studied. It is found that the case of higher spacing between roughness elements is more effective for inducing transition than the lower one. The interaction between two adjacent roughness elements can suppress the evolution of horseshoe vortices in the downstream and trigger the instability of shear layer. Thus, the transition will be suppressed accordingly.
文摘Our previous study showed that the frictional drag decreases with increasing void fraction at Re〉1300, while it increases at Re 〈 1000. Decomposition of the Reynolds shear stress also implied that bubbles induce isotropy of turbulence. In order to confirm our previous analysis and to further investigate flow fields in the vicinity of bubbles, we analyze velocity fluctuations on the quadrant space in the streamwise and transverse directions (u′-v′ plane). Here, we focus on two specific Reynolds numbers (at Re≈900 and ≈1410, which are close to the laminar-to-turbulent transition regime) and discuss bubble effects on sweep (u′〉 0, v′〈 0 ) and ejection (u′〈 0, v′〉 0) events as a function of the Reynolds number. We also illustrate velocity fluctuations in the vicinity of an individual bubble and a swarm of bubbles on the u′- v′ coordinates. The results show that a bubble swarm suppresses the velocity fluctuations at Re≈1410.
基金Supported by the National Natural Science Foundation of China (Grant Nos. 10632050 and 10502052) Informatization Construction of Knowledge Innovation Projects of the Chinese Academy of Sciences (Grant No. INF105-SCE)
文摘Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20° cone angle (or 10° half-cone angle) is performed by using the mixed seventh- order up-wind biased finite difference scheme and sixth-order central difference scheme. The free stream Mach number is 0.7 and free stream unit Reynolds number is 250000/inch. The characteristics of transition and turbulence of the sharp cone boundary layer are compared with those of the flat plate boundary layer. Statistics of fully developed turbulent flow agree well with the experimental and theoretical data for the turbulent flat-plate boundary layer flow. The near wall streak-like structure is shown and the average space between streaks (normalized by the local wall unit) keeps approximately invariable at different streamwise locations. The turbulent energy equation in the cylindrical coordinate is given and turbulent en-ergy budget is studied. The computed results show that the effect of circumferen-tial curvature on turbulence characteristics is not obvious.
基金supported by the National Natural Science Foundation of China (No.11372254)
文摘The achievement of laminar flow in the boundary layer at high-speed cruise conditions may further, in addition to shock-wave control, reduce the drag and extend the range of military fighter aircraft. To this end, a further investigation on transitional boundary-layer flow of fighter wings is needed due to different configurations from the wings used on conventional transport aircraft. In this paper, wind tunnel experiments and numerical simulations were conducted on three-dimensional transition of thin diamond-shaped wings used on advanced fighter aircraft at tran/supersonic design points. A newly proposed correlation of crossflow transition which includes the effect of surface roughness was introduced into the c-Rehttransition model. Predicted results were in good agreement with flow visualizations. Results showed that the strength of the crossflow component grew rapidly around the leading edge because of the severe flow acceleration, just as the same as wings with a large aspect ratio. However, there seemed no regular pattern of instabilitydominance variation in span-wise for a diamond configuration. The dominance of different instability mechanisms strongly depended on the local pressure distribution. Hereby, the research recommended a ‘‘roof-like" shape of pressure distribution to suppress both crossflow and Tollmien-Schlichting(T-S) instabilities. Besides, a sharp suction peak with a serious pressure rise should be cut off to avoid stronger instabilities. Further discussions also revealed an independence of the unit Reynolds number when transition was triggered by T-S instabilities. Aerodynamic force comparisons indicated that further benefit on drag reduction could be expected by including the three-dimensional transition effect into a wing design process.
文摘This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transport equations, separately, to compute laminar kinetic energy, turbulent kinetic energy, and dissipation rate in a flow field. It needs neither correlations of intermittency factors nor knowledge of the transition onset. Two transition tests are carried out: flat plate boundary layer under zero ...
基金supported by the National Natural Science Foundation of China (Grant No. 10772082)AFOSR (Grant No. FA9550-08-1-0201)
文摘At the late stage of transitional boundary layers, the nonlinear evolution of the ring-like vortices and spike structures and their effects on the surrounding flow were studied by means of direct numerical simulation with high order accuracy. A spatial transition of the flat-plate boundary layers in the compressible flow was conducted. Detailed numerical results with high resolution clearly represented the typical vortex structures, such as ring-like vortices and so on, and induced ejection and sweep events. It was verified that the formation of spike structures in transitional boundary layers had close relationship with ring-like vortices. Especially, compared to the newly observed positive spike structure in the experiments, the same structure was found in the present numerical simulations, and the mechanism was also studied and analyzed.
基金Project supported by the National Natural Science Foundation of China (No.90205021)the Special Foundation for Doctoral Dissertations Research (No.200328)the Science Foundation of Liuhui Center of Applied Mathematics, Nankai University and Tianjin University
文摘Spatial mode direct numerical simulation has been applied to study the mechanism of breakdown in laminar-turbulent transition of a supersonic boundary layer on a fiat plate with Mach number 4.5. Analysis of the result showed that, during the breakdown process in laminar-turbulent transition, the mechanism causing the mean flow profile to evolve swiftly from laminar to turbulent was that the modification of mean flow profile by the disturbance, when they became larger, leads to remarkable change of its stability characteristics. Though the most unstable T-S wave was of second mode for laminar flow, the first mode waves played the key role in the breakdown process in laminar-turbulent transition.
文摘The thermal barrier coatings with NiCrAlY alloy bonding layer, NiCrAlY Y 2O 3 stabilized ZrO 2 transition layer and Y 2O 3 stabilized ZrO 2 ceramic layer are prepared on nickel alloy substrates using the plasma spray technique. The relationship among the composition, structure and property of the coatings are investiga ted by means of optical microscope, scanning electronic microscope and the experiments of thermal shock resistance cycling and high temperature oxidation resistance. The results show that the structure design of introdu cing a transition layer between Ni alloy substrate and ZrO 2 ceramic coating guarantees the high quality and properties of the coatings; ZrO 2 coatings doped with a little SiO 2 possesses better thermal shock resistance and more excellent hot corrosion resistance as compared with ZrO 2 coating materials without SiO 2 ;the improvement in performance of ZrO 2 coating doped with SiO 2 is due to forming more dense coating structure by self closing effects of the flaws and pores in the ZrO 2 coatings.
基金the Russian Foundation for Basic Research (Grant 08-01-91956-NNIO)ADTP RNP 2.1.1/3963Program RAS (project 11/9)
文摘Experimental study was conducted for boundarylayers on a sharp 5° half-angle cone of 400mm length at angles of attack. The model was tested in the T-326 hypersonic wind tunnel (ITAM) at freestream Mach number M = 5.95. Mean and fluctuation wall characteristics of the boundary layer are measured at 0°, 2°, 3° and 4° angles of attack for different stagnation pressures. Pulsation measurements are carried out by means of ALTP sensor. Pressure and temperature distributions along the model are obtained, and transition beginning and end locations have been found. Boundary layer stabilization with the increase of angle of attack and the decrease of stagnation pressure is observed. High frequency pulsations inherent to hypersonic boundary layer (second mode) have been detected.
基金This work was supported by NKBRSF(Grant No.CG199032805)the National Natural Science Foundation of China(Grant Nos.90205025,19872069,170176033)the Informatization Construction of Knowledge Innovation Projects of the Chinese Academy of Sciences(Grant No.INF105-SCE).
文摘Direct numerical simulation of spatially evolving compressible boundary layer over a blunt wedge is performed in this paper. The free-stream Mach number is 6 and the disturbance source produced by wall blowing and suction is located downstream of the sound-speed point. Statistics are studied and compared with the results in incompressible flat-plate boundary layer. The mean pressure gradient effects on the vortex structure are studied.