The characteristics of turbulent boundary layer over streamwise aligned drag reducing riblet surface under zero-pressure gradient are investigated using particle image velocimetry. The formation and distribution of la...The characteristics of turbulent boundary layer over streamwise aligned drag reducing riblet surface under zero-pressure gradient are investigated using particle image velocimetry. The formation and distribution of large-scale coherent structures and their effect on momentum partition are analyzed using two-point correlation and probability density function. Compared with smooth surface, the streamwise riblets reduce the friction velocity and Reynolds stress in the turbulent boundary layer, indicating the drag reduction effect. Strong correlation has been found between the occurrence of hairpin vortices and the momentum distribution. The number and streamwise length scale of hairpin vortices decrease over streamwise riblet surface. The correlation between number of uniform momentum zones and Reynolds number remains the same as smooth surface.展开更多
The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parab...The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations (PSE). Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.展开更多
Bent inlet pipes are often used in centrifugal compressors due to limited installation space,and an understanding of the effect on compressor stability is essential for safety and durability.This paper firstly investi...Bent inlet pipes are often used in centrifugal compressors due to limited installation space,and an understanding of the effect on compressor stability is essential for safety and durability.This paper firstly investigates flow instability behaviors in two compressors,one with a straight inlet pipe and the other with an S-shaped bent pipe.In detail,it analyzes the resulting flow fields,instability evolution paths and surge boundaries.The results show that the S-shaped pipe obviously affects the flow field at high mass flow rates,while reverse flow mainly influences the flow field at low mass flow rates.Reverse flow first occurs at certain flow passages with a high pressure difference that is predominantly decided by the volute rather than the S-shaped bent pipe.In addition,centrifugal compressors can tolerate reverse flow to some extent so that surge would not occur immediately if reverse flow occurs unless the reverse flow region extends circumferentially and radially to a sufficiently large size.Since the S-shaped pipe is not dominant in the creation and extension of reverse flow,it does not exacerbate the stability of the central compressor to a great extent.Last but not least,the S-shaped pipe is noted to delay the occurrence of surge at 90%rotating speed,which suggests the possibility of improving compressor stability with bent inlet pipes.This result differs from the conventional understanding that inlet distortion usually deteriorates compressor stability and emphasizes the particularity of centrifugal compressors.展开更多
Two techniques that improve the aerodynamic performance of wind turbine airfoils are described. The airfoil $809, designed specially for wind turbine blades, and the airfoil FX60-100, having a higher lift-drag ratio, ...Two techniques that improve the aerodynamic performance of wind turbine airfoils are described. The airfoil $809, designed specially for wind turbine blades, and the airfoil FX60-100, having a higher lift-drag ratio, are selected to verify the flow control techniques. The flow deflector, fixed at the leading edge, is employed to control the boundary layer separation on the airfoil at a high angle of attack. The multi-island genetic algorithm is used to optimize the parameters of the flow deflector. The results indicate that the flow deflector can suppress the flow separation, delay the stall, and enhance the lift. The characteristics of the blade tip vortex, the wake vortex, and the surface pressure distributions of the blades are analyzed. The vortex diffuser, set up at the blade tip, is employed to control the blade tip vortex. The results show that the vortex diffuser can increase the total pressure coefficient of the core of the vortex, decrease the strength of the blade tip vortex, lower the noise, and improve the efficiency of the blade.展开更多
The nonlinear evolution problem in nonparallel boundary layer stability was studied. The relative parabolized stability equations of nonlinear nonparallel boundary layer were derived. The developed numerical method, w...The nonlinear evolution problem in nonparallel boundary layer stability was studied. The relative parabolized stability equations of nonlinear nonparallel boundary layer were derived. The developed numerical method, which is very effective, was used to study the nonlinear evolution. of T-S disturbance wave at finite amplitudes. Solving nonlinear equations of different modes by using predictor-corrector and iterative approach, which is uncoupled between modes, improving computational accuracy by using high order compact differential scheme, satisfying normalization condition I determining tables of nonlinear terms at different modes, and implementing stably the spatial marching, were included in this method. With different initial amplitudes, the nonlinear evolution of T-S wave was studied. The nonlinear nonparallel results of examples compare with data of direct numerical simulations (DNS) using full Navier-Stokes equations.展开更多
Formation and evolution of secondary streamwise vortices in the compressible transitional boundary layers over a flat plate are studied using a direct numerical simu- lation method with high-order accuracy and highly ...Formation and evolution of secondary streamwise vortices in the compressible transitional boundary layers over a flat plate are studied using a direct numerical simu- lation method with high-order accuracy and highly effective non-reflecting characteristic boundary conditions. Generation and development processes of the secondary streamwise vortices in the complicated transitional boundary flow are clearly analyzed based on the of numerical results, and the effects on the formation of the ring-like vortex that is vital to the boundary layer transition are explored. A new mechanism forming the ring-like vortex through the mutual effect of the primary and secondary streamwise vortices is expressed.展开更多
The parallel acceleration of well-developed serial codes for numerical simulations of fluid dynamic problems is implemented in this paper. To solve the flow field,the lattice Boltzmann method( LBM) is used. The OpenAC...The parallel acceleration of well-developed serial codes for numerical simulations of fluid dynamic problems is implemented in this paper. To solve the flow field,the lattice Boltzmann method( LBM) is used. The OpenACC( OpenACCelerator) application programming interface,which is a relatively new standard for parallel computation,is adopted to perform the acceleration. By parallelizing the computation intensive loops and decreasing some meaningless data movement,the serial codes based on LBM are obviously accelerated. Several benchmark problems are considered and then simulated through the OpenACC programming model to assess performance and computational efficiency. After optimization,a reasonable speedup is obtained in comparison with the original version.展开更多
In order to assess the influence of high vibrational temperatures typically found in shock-tunnel-generated flows on model surface heat flux, measurements are compared to simulations using the DSMC (direct simulation...In order to assess the influence of high vibrational temperatures typically found in shock-tunnel-generated flows on model surface heat flux, measurements are compared to simulations using the DSMC (direct simulation Monte Carlo) method. The two limiting cases of either full accommodation of molecular vibrational energy to the surface temperature or none at all are simulated for three shock tunnel conditions using Nitrogen (N2) as a test gas. The conditions mainly differ by an order of magnitude in density, however, all three comparisons to the corresponding measurements suggest that vibrational surface relaxation happens only to a small extent. While model surface pressure is not affected, heat fluxes differ by up to 20%, depending on the modelling. Furthermore, lower flow densities generally lead to higher differences.展开更多
This paper presents a novel evasion guidance law for hypersonic morphing vehicles,focusing on determining the optimized wing's unfolded angle to promote maneuverability based on an intelligent algorithm.First,the ...This paper presents a novel evasion guidance law for hypersonic morphing vehicles,focusing on determining the optimized wing's unfolded angle to promote maneuverability based on an intelligent algorithm.First,the pursuit-evasion problem is modeled as a Markov decision process.And the agent's action consists of maneuver overload and the unfolded angle of wings,which is different from the conventional evasion guidance designed for fixed-shape vehicles.The reward function is formulated to ensure that the miss distances satisfy the prescribed bounds while minimizing energy consumption.Then,to maximize the expected cumulative reward,a residual learning method is proposed based on proximal policy optimization,which integrates the optimal evasion for linear cases as the baseline and trains to optimize the performance for nonlinear engagement with multiple pursuers.Therefore,offline training guarantees improvement of the constructed evasion guidance law over conventional ones.Ultimately,the guidance law for online implementation includes only analytical calculations.It maps from the confrontation state to the expected angle of attack and the unfolded angle while retaining high computational efficiency.Simulations show that the proposed evasion guidance law can utilize the change of unfolded angle to extend the maximum overload capability.And it surpasses conventional maneuver strategies by ensuring better evasion efficacy and higher energy efficiency.展开更多
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.展开更多
Geometric and working condition uncertainties are inevitable in a compressor,deviating the compressor performance from the design value.It’s necessary to explore the influence of geometric uncertainty on performance ...Geometric and working condition uncertainties are inevitable in a compressor,deviating the compressor performance from the design value.It’s necessary to explore the influence of geometric uncertainty on performance deviation under different working conditions.In this paper,the geometric uncertainty influences at near stall,peak efficiency,and near choke conditions under design speed and low speed are investigated.Firstly,manufacturing geometric uncertainties are analyzed.Next,correlation models between geometry and performance under different working conditions are constructed based on a neural network.Then the Shapley additive explanations(SHAP)method is introduced to explain the output of the neural network.Results show that under real manufacturing uncertainty,the efficiency deviation range is small under the near stall and peak efficiency conditions.However,under the near choke conditions,efficiency is highly sensitive to flow capacity changes caused by geometric uncertainty,leading to a significant increase in the efficiency deviation amplitude,up to a magnitude of-3.6%.Moreover,the tip leading-edge radius and tip thickness are two main factors affecting efficiency deviation.Therefore,to reduce efficiency uncertainty,a compressor should be avoided working near the choke condition,and the tolerances of the tip leading-edge radius and tip thickness should be strictly controlled.展开更多
A fast hybrid algorithm based on gridless method coupled with finite volume method (FVM) is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm i...A fast hybrid algorithm based on gridless method coupled with finite volume method (FVM) is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm is improved to the level of finite volume method for most parts of the flow filed arc covered with grid cells. Moreover, the hybrid method is flexible to deal with the configurations as clouds of points are used to cover the region adjacent to the bodies. Mirror satellites and mirror grid cells arc introduced to the interface to accomplish data communication between the different parts of the flow field. The Euler Equations arc spatially discretized with finite volume method and gridless method in mesh and clouds of points respectively, and an explicit four-stage Runge-Kutta scheme is utilized to reach the steady-state solution. Internal flows in channels and external flows over airfoils arc investigated with hybrid method, and the solutions arc comparad to those using pure finite volume method and pure gridless method. Numerical examples show that the hybrid algorithm captures the shock waves accurately, and it is as efficient as fmite volume method.展开更多
High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is investigated in combination with discont...High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is investigated in combination with discontinuous Galerkin(DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry representation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomerated high-order mesh works well with DG method in dealing with flow problems of curved geometries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant efficiency improvements compared to its explicit counterparts.展开更多
A coupled immersed boundary-lattice Boltzmann method (IB-LBM) is introduced to solve biomimetic problems. Compared to the conventional IB-LBM, the strict satisfaction of no-slip boundary condition is implemented in ...A coupled immersed boundary-lattice Boltzmann method (IB-LBM) is introduced to solve biomimetic problems. Compared to the conventional IB-LBM, the strict satisfaction of no-slip boundary condition is implemented in the current method. Consequently, the phenomenon of flow penetration that is frequently observed in the conventional IB-LBM is fully prevented, and subsequently the force on the boundary can be calculated more accurately. This feature is of importance for the simulation of biomimetic problems. Moreover, by applying the relationship between the velocity correction and forcing term, the boundary force can be calculated easily. Several biomimetic problems are then simulated. Based on the good agreement between the current results and those in the literature, it may be concluded that the present IB-LBM has the capabilitv to handle various biomimetic oroblems.展开更多
The article represents the discussion of several separate directions of investigations,which are performed by TsAGI flight vehicles aerodynamics specialists at the time.There are some major trends of classical layout ...The article represents the discussion of several separate directions of investigations,which are performed by TsAGI flight vehicles aerodynamics specialists at the time.There are some major trends of classical layout of route aircraft and also peculiarities of some prospective flight vehicles.Also there are some hypersonic vehicles aerodynamics questions examined along with problems of creation of civil supersonic transport aircraft.There is a description given for well-known and some newer methods of flow control for drag reduction.展开更多
An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. I...An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.展开更多
基金supported by the National Natural Science Foundation of China (Nos. 11721202 and 11672020)
文摘The characteristics of turbulent boundary layer over streamwise aligned drag reducing riblet surface under zero-pressure gradient are investigated using particle image velocimetry. The formation and distribution of large-scale coherent structures and their effect on momentum partition are analyzed using two-point correlation and probability density function. Compared with smooth surface, the streamwise riblets reduce the friction velocity and Reynolds stress in the turbulent boundary layer, indicating the drag reduction effect. Strong correlation has been found between the occurrence of hairpin vortices and the momentum distribution. The number and streamwise length scale of hairpin vortices decrease over streamwise riblet surface. The correlation between number of uniform momentum zones and Reynolds number remains the same as smooth surface.
文摘The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations (PSE). Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.
基金co-supported by the Tsinghua University"Shuimu Tsinghua Scholar"Programthe National Science and Technology Major Project+2 种基金China(No.2017-II-0004-0016)the National Natural Science Foundation of China(No.51876097)the IHI Corporation,Yokohama,Japan。
文摘Bent inlet pipes are often used in centrifugal compressors due to limited installation space,and an understanding of the effect on compressor stability is essential for safety and durability.This paper firstly investigates flow instability behaviors in two compressors,one with a straight inlet pipe and the other with an S-shaped bent pipe.In detail,it analyzes the resulting flow fields,instability evolution paths and surge boundaries.The results show that the S-shaped pipe obviously affects the flow field at high mass flow rates,while reverse flow mainly influences the flow field at low mass flow rates.Reverse flow first occurs at certain flow passages with a high pressure difference that is predominantly decided by the volute rather than the S-shaped bent pipe.In addition,centrifugal compressors can tolerate reverse flow to some extent so that surge would not occur immediately if reverse flow occurs unless the reverse flow region extends circumferentially and radially to a sufficiently large size.Since the S-shaped pipe is not dominant in the creation and extension of reverse flow,it does not exacerbate the stability of the central compressor to a great extent.Last but not least,the S-shaped pipe is noted to delay the occurrence of surge at 90%rotating speed,which suggests the possibility of improving compressor stability with bent inlet pipes.This result differs from the conventional understanding that inlet distortion usually deteriorates compressor stability and emphasizes the particularity of centrifugal compressors.
基金Project supported by the National Basic Research Program of China (973 Program) (No. 2007CB714600)
文摘Two techniques that improve the aerodynamic performance of wind turbine airfoils are described. The airfoil $809, designed specially for wind turbine blades, and the airfoil FX60-100, having a higher lift-drag ratio, are selected to verify the flow control techniques. The flow deflector, fixed at the leading edge, is employed to control the boundary layer separation on the airfoil at a high angle of attack. The multi-island genetic algorithm is used to optimize the parameters of the flow deflector. The results indicate that the flow deflector can suppress the flow separation, delay the stall, and enhance the lift. The characteristics of the blade tip vortex, the wake vortex, and the surface pressure distributions of the blades are analyzed. The vortex diffuser, set up at the blade tip, is employed to control the blade tip vortex. The results show that the vortex diffuser can increase the total pressure coefficient of the core of the vortex, decrease the strength of the blade tip vortex, lower the noise, and improve the efficiency of the blade.
文摘The nonlinear evolution problem in nonparallel boundary layer stability was studied. The relative parabolized stability equations of nonlinear nonparallel boundary layer were derived. The developed numerical method, which is very effective, was used to study the nonlinear evolution. of T-S disturbance wave at finite amplitudes. Solving nonlinear equations of different modes by using predictor-corrector and iterative approach, which is uncoupled between modes, improving computational accuracy by using high order compact differential scheme, satisfying normalization condition I determining tables of nonlinear terms at different modes, and implementing stably the spatial marching, were included in this method. With different initial amplitudes, the nonlinear evolution of T-S wave was studied. The nonlinear nonparallel results of examples compare with data of direct numerical simulations (DNS) using full Navier-Stokes equations.
基金supported by the National Natural Science Foundation of China (No. 10772082)AFOSR(No. FA9550-08-1-0201)
文摘Formation and evolution of secondary streamwise vortices in the compressible transitional boundary layers over a flat plate are studied using a direct numerical simu- lation method with high-order accuracy and highly effective non-reflecting characteristic boundary conditions. Generation and development processes of the secondary streamwise vortices in the complicated transitional boundary flow are clearly analyzed based on the of numerical results, and the effects on the formation of the ring-like vortex that is vital to the boundary layer transition are explored. A new mechanism forming the ring-like vortex through the mutual effect of the primary and secondary streamwise vortices is expressed.
基金Sponsored by the Research Fund of State Key Laboratory of Mechanics and Control of Mechanical Structures(Nanjing University of Aeronautics and Astronautics)(Grant No.MCMS-0117G01)
文摘The parallel acceleration of well-developed serial codes for numerical simulations of fluid dynamic problems is implemented in this paper. To solve the flow field,the lattice Boltzmann method( LBM) is used. The OpenACC( OpenACCelerator) application programming interface,which is a relatively new standard for parallel computation,is adopted to perform the acceleration. By parallelizing the computation intensive loops and decreasing some meaningless data movement,the serial codes based on LBM are obviously accelerated. Several benchmark problems are considered and then simulated through the OpenACC programming model to assess performance and computational efficiency. After optimization,a reasonable speedup is obtained in comparison with the original version.
文摘In order to assess the influence of high vibrational temperatures typically found in shock-tunnel-generated flows on model surface heat flux, measurements are compared to simulations using the DSMC (direct simulation Monte Carlo) method. The two limiting cases of either full accommodation of molecular vibrational energy to the surface temperature or none at all are simulated for three shock tunnel conditions using Nitrogen (N2) as a test gas. The conditions mainly differ by an order of magnitude in density, however, all three comparisons to the corresponding measurements suggest that vibrational surface relaxation happens only to a small extent. While model surface pressure is not affected, heat fluxes differ by up to 20%, depending on the modelling. Furthermore, lower flow densities generally lead to higher differences.
基金This work was supported by the National Natural Science Foundation of China(No.52202438).
文摘This paper presents a novel evasion guidance law for hypersonic morphing vehicles,focusing on determining the optimized wing's unfolded angle to promote maneuverability based on an intelligent algorithm.First,the pursuit-evasion problem is modeled as a Markov decision process.And the agent's action consists of maneuver overload and the unfolded angle of wings,which is different from the conventional evasion guidance designed for fixed-shape vehicles.The reward function is formulated to ensure that the miss distances satisfy the prescribed bounds while minimizing energy consumption.Then,to maximize the expected cumulative reward,a residual learning method is proposed based on proximal policy optimization,which integrates the optimal evasion for linear cases as the baseline and trains to optimize the performance for nonlinear engagement with multiple pursuers.Therefore,offline training guarantees improvement of the constructed evasion guidance law over conventional ones.Ultimately,the guidance law for online implementation includes only analytical calculations.It maps from the confrontation state to the expected angle of attack and the unfolded angle while retaining high computational efficiency.Simulations show that the proposed evasion guidance law can utilize the change of unfolded angle to extend the maximum overload capability.And it surpasses conventional maneuver strategies by ensuring better evasion efficacy and higher energy efficiency.
基金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.
基金supported by the National Science and Technology Major Project,China(No.2017-II-0004-0016)。
文摘Geometric and working condition uncertainties are inevitable in a compressor,deviating the compressor performance from the design value.It’s necessary to explore the influence of geometric uncertainty on performance deviation under different working conditions.In this paper,the geometric uncertainty influences at near stall,peak efficiency,and near choke conditions under design speed and low speed are investigated.Firstly,manufacturing geometric uncertainties are analyzed.Next,correlation models between geometry and performance under different working conditions are constructed based on a neural network.Then the Shapley additive explanations(SHAP)method is introduced to explain the output of the neural network.Results show that under real manufacturing uncertainty,the efficiency deviation range is small under the near stall and peak efficiency conditions.However,under the near choke conditions,efficiency is highly sensitive to flow capacity changes caused by geometric uncertainty,leading to a significant increase in the efficiency deviation amplitude,up to a magnitude of-3.6%.Moreover,the tip leading-edge radius and tip thickness are two main factors affecting efficiency deviation.Therefore,to reduce efficiency uncertainty,a compressor should be avoided working near the choke condition,and the tolerances of the tip leading-edge radius and tip thickness should be strictly controlled.
基金Aeronautical Science Foundation of China (02A52002), National Natural Science Foundation of China(10372043)
文摘A fast hybrid algorithm based on gridless method coupled with finite volume method (FVM) is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm is improved to the level of finite volume method for most parts of the flow filed arc covered with grid cells. Moreover, the hybrid method is flexible to deal with the configurations as clouds of points are used to cover the region adjacent to the bodies. Mirror satellites and mirror grid cells arc introduced to the interface to accomplish data communication between the different parts of the flow field. The Euler Equations arc spatially discretized with finite volume method and gridless method in mesh and clouds of points respectively, and an explicit four-stage Runge-Kutta scheme is utilized to reach the steady-state solution. Internal flows in channels and external flows over airfoils arc investigated with hybrid method, and the solutions arc comparad to those using pure finite volume method and pure gridless method. Numerical examples show that the hybrid algorithm captures the shock waves accurately, and it is as efficient as fmite volume method.
基金co-supported by the Aeronautical Science Foundation of China (No. 20152752033)the National Natural Science Foundation of China (No. 11272152)the Open Project of Key Laboratory of Aerodynamic Noise Control
文摘High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is investigated in combination with discontinuous Galerkin(DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry representation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomerated high-order mesh works well with DG method in dealing with flow problems of curved geometries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant efficiency improvements compared to its explicit counterparts.
基金supported by the National Natural Science Foundation of China(11272153)
文摘A coupled immersed boundary-lattice Boltzmann method (IB-LBM) is introduced to solve biomimetic problems. Compared to the conventional IB-LBM, the strict satisfaction of no-slip boundary condition is implemented in the current method. Consequently, the phenomenon of flow penetration that is frequently observed in the conventional IB-LBM is fully prevented, and subsequently the force on the boundary can be calculated more accurately. This feature is of importance for the simulation of biomimetic problems. Moreover, by applying the relationship between the velocity correction and forcing term, the boundary force can be calculated easily. Several biomimetic problems are then simulated. Based on the good agreement between the current results and those in the literature, it may be concluded that the present IB-LBM has the capabilitv to handle various biomimetic oroblems.
文摘The article represents the discussion of several separate directions of investigations,which are performed by TsAGI flight vehicles aerodynamics specialists at the time.There are some major trends of classical layout of route aircraft and also peculiarities of some prospective flight vehicles.Also there are some hypersonic vehicles aerodynamics questions examined along with problems of creation of civil supersonic transport aircraft.There is a description given for well-known and some newer methods of flow control for drag reduction.
基金co-supported by the National High Technology Research and Development Program of China(No.2015AA015303)the National Natural Science Foundation of China(No.11272152)+1 种基金the Aeronautical Science Foundation of China(No.20152752033)the Open Project of Key Laboratory of Aerodynamic Noise Control
文摘An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.