摘要
昆虫和鸟类的飞行性能极为高超,其飞行过程给予我们很多的启示。根据生物仿生学的原理,同时借助于空气动力学技术的发展,各种新型的变形飞行器应运而生。如今,高强度固定翼飞行器均对特定的飞行条件进行了优化设计,但在非设计状态往往难以达到较高的飞行效率。结构设计的困难、气动弹性的不稳定性、机翼蒙皮材料的性能都限制了飞行器的飞行包线进一步扩大,为此科研工作者提出了机翼形状控制[1](wing shape control)的概念。例如,洛克希德·马丁公司提出的折叠翼无人战斗飞行器(UCAV),具有长时间的巡航能力,能够通过折叠机翼快速进入俯冲和战斗状态,战斗结束后又可转换为长时间的巡航状态。这种构型有两处可以折叠,能够提供内翼折叠130度、最多200%的翼面积变形。变形使得巡航/侦察的功效和俯冲/突防能力集于一体,相比于常规作战平台显著提高了作战能力。本文的目的就是希望建立一套较为完善的柔性变形体数值模拟方法,数值模拟变形飞机机翼折叠/展开过程的非定常流动,研究在机翼变形过程中变形飞机的动态气动特性。
变形飞机机翼的折叠/展开运动是一个强非定常、非线性过程,其计算网格必须随机翼的运动而发生大变形。因此对变形飞机流场的数值模拟,需要建立合适的动态网格生成方法,并建立高效的可压缩流非定常计算方法。由此,本文设计了如下的组织结构:
在第一章的引言中,介绍了关于飞行器的多种变形方案以及国内外研究进展,对混合网格生成方法、动态网格生成方法和非定常数值计算方法做了简要的回顾,最后介绍了本文的工作。
第二章重点介绍了静态/动态混合网格生成方法,包含了动态混合网格生成的整体方案,静态混合网格生成方法,Delaunay背景网格生成方法,以及基于Delaunay背景网格插值法的局部网格重构等问题,并给出了本文数值计算所采用的动态混合网格的生成实例。
第三章介绍了基于上述动态混合网格的非定常计算方法,包括流动控制方程,空间和时间离散方法,动网格计算的几何守恒率,边界条件的处理等。在章节末尾给出了定常数值验证结果和典型非定常数值验证结果,计算结果与文献数据的吻合证明了本文采用的计算方法有较高的计算精度和计算效率。
第四章研究了局部翼面主动变形对跨声速翼型气动特性的影响,试图根据“流动滚轴效应”的思想,通过物面的动态变形,改变跨声速条件下主导翼型升力系数的上翼面激波出现的位置,从而提高翼型升力系数。为此,我们将局部行波壁模型和局部主动振动模型加载至NACA0012翼型和超临界翼型RAE2822上,考察了不同参数的计算状态。计算结果表明,局部主动变形模型带来的翼面局部变形在大部分情况下导致上翼面激波前移,从而降低了翼型的升力系数。由此表明,局部主动振动模型应用于跨声速条件下后移激波的方法是不成功的。
第五章研究局部翼面主动变形对低速翼型气动特性的影响,探讨了在低速多段翼型30P30N的主翼段上表面加载局部主动变形模型后翼型的气动特性,主要关注升力特性的影响,包括局部左行波壁模型、右行波壁模型、左行合成波壁模型、右行合成波壁模型和局部主动振动模型对翼型升力特性的影响。对于翼型在失速状态下加入前述5个变形模型,翼型都有不同程度的增升作用。本章讨论的各种局部主动变形模型应用于翼型的失速控制对三维机翼的失速控制有一定的参考价值。
第六章进行了折叠翼变形飞机动态气动特性及非定常流动机理的初步研究,本章首先参考文献资料设计了简化的折叠翼变形飞机模型,然后利用前述方法生成了非结构动态网格,最后利用第三章的非定常数值计算方法进行了无粘流数值模拟。数值计算结果表明,变形飞机机翼在折叠/展开过程中,飞机的升力特性变化较大,上翼面的激波位置在机翼折叠/展开过程中均会发生前移,从而降低了升力系数。
最后,在结束语中,对现有工作进行了总结,指出了存在的不足和今后的研究方向。
The flying capability of insects and birds is so excellent that we can benefit abundantly from their flying behavior. According to the principles of bionics, lots of new morphing aircrafts have been designed with the development of aerodynamic technology. As one follows the evolution of the aircraft, however, fixed geometry, high stiffness wings emerge as the dominant trait. As a result, air vehicles are optimized for specific flight conditions, such as the cruising condition. Flight dynamiscists have realized that changes in geometry of flight vehicles could provide significant improvement in flight efficiency. Therefore, researchers brought up the concept of wing shape control[1]. As an example, there comes the morphing aircraft with folding wings proposed by the Lockheed Martin Company. It is an unmanned combat air vehicle (UCAV), capable of long range cruising with unfolding wings, and fast diving with folded wings, and transition back for the long range cruising. The configuration has two foldable wing-sections that allow radical morphing of span and wetted area. The inner wing can fold approximately 130 degrees from the unfolded to the folded configuration, while the outer wing keeps the original horizontal status to achieve the 200% wing area changing. Comparing with conventional aircrafts, configuration morphing improves significantly the mission performance. The goal of this thesis is to simulate the unsteady flows and study the dynamic aerodynamic characteristics of the morphing aircraft during the movement of its wing’s folding/unfolding.
The folding/unfolding movement of morphing aircraft is of strongly unsteady and nonlinear, so its computational grids have to be deformed with the wing morphing. Therefore, to simulate the unsteady flows over morphing aircraft, dynamic grid technique and corresponding efficient unsteady flow solver should be set up firstly. So this thesis is organized as follows:
In the preface as chapter one of this thesis, the research background is discussed firstly in brief. Then the progress in hybrid mesh generation techniques, dynamic grid generation methods and the numerical methods for unsteady flows are reviewed. And the work of this paper is presented finally.
In the second chapter, we introduce the static/dynamic hybrid mesh generation technique, including the fundamental strategy of dynamic hybrid mesh generation, static hybrid mesh generation, Delaunay graph generation and the coupling dynamic hybrid grid generation method of Delaunay graph mapping and local remeshing. At the end of this chapter, some examples of grid generation are shown, which are used in the following numerical simulations of this thesis.
In the third chapter, an unsteady flow solver based on the hybrid dynamic grid is proposed, including the governing equations, spatial and time discretization scheme, geometric conservation law (GCL) for moving grid simulations, and boundary conditions. Some typical steady/unsteady test cases are simulated to validate the unsteady flow solver and the dynamic hybrid grid approach. The numerical results demonstrate the efficiency and accuracy of present method.
In the fourth chapter, the unsteady flows over transonic morphing airfoils are simulated by adding some morphing models to the original airfoils (NACA0012 and RAE2822). From the idea of flow-roller effect in the boundary layer, we wish that these models could push backward the position of shock-wave on the leeward of the airfoils to increase the lift coefficient. Unfortunately, most of the results show that the morphing models pull the shock wave on the airfoil leeward forward, resulting in the decreasing of the lift coefficient. So these morphing models are not suitable for transonic airfoils with shock wave dominating the flow fields.
In the fifth chapter, the morphing models are adopted on a typical subsonic three-element airfoil (30P30N) to find the helpful effects on the airfoil aerodynamic characteristics. By comparing the numerical results of left-travel wave model, right-travel wave model, left-travel compound wave model, right-travel compound wave model and local oscillation model, we find that the lift coefficients of 30P30N airfoil after stall can increase more or less.
In the sixth chapter, numerical simulations over a morphing aircraft with foldable wings are carried out, and the dynamic aerodynamic characteristics are analyzed. Firstly, we design a simplified model of morphing aircraft according to the references, and then dynamic unstructured grids are generated by the dynamic hybrid grid generation technique presented in chapter two. After that, numerical simulations are undertaken with the unsteady flow solver, here the inviscid cases are considered only. The computational results demonstrate that the aerodynamic characteristics of the morphing aircraft change rapidly with the shock wave on the leeward of the wing moving forward during the wing folding/unfolding.
In the last chapter, some concluding remarks are given, and the future work is discussed.
引文
[1]. Matthew P. Snyder, Brian Sanders, Franklin E. Eastep, and Geoffrey J. Frank, Vibration and Flutter Characteristics of a Folding Wing, AIAA/ASME/ASCE/ AHS/ASC Structures, Structural Dynamics & Materials Confer,18-21 April 2005, Austin, Texas
[2]. John V., Monish D. T. and Jie R., A reinforcement learning - adaptive control architecture for morphing, AIAA J. 20-22 September 2004-6220
[3]. B. Sanders, F.E. Eastep and E. Foster, Aerodynamic and Aeroelastic Characteristics of Wings with Conformal Control Surfaces for Morphing Aircraft, Journal of Aircraft, 2003, 40(1):94-99,
[4]. Amprikidis M., and Cooper J. E., Development of Smart Spars for Active Aeroelastic Structures, AIAA-2003-1799
[5]. Brian R., Chris P., William A. C., Aircraft sizing with morphing as an independent variable - Motivation, strategies and investigations ,AIAA J. 2002-5840 October 2002
[6].张颖,未来航空器:变形飞机,外国科技动态,2006, (1):43-48
[7]. John B. Davidson, Pawel Chwalowski, Barry S. Lazos, Flight Dynamic Simulation Assessment of a Morphable Hyper-elliptic Cambered Winged Configuration, AIAA-2003-5301
[8]. Leonard D. Wiggins, Matthew D. Stubb, Christopher O. Johnston, Harry H. Robertshaw, Charles F. Reinholtz, and Daniel J. Inman, A Design and Analysis of a Morphing Hyper-Elliptic Cambered Span (HECS) Wing, AIAA-2004-1885
[9]. Cone C.D., The Theory of Induced Lift and Minimum Induced Drag on Non-Planar Lifting Systems, NASA TR-R-139, 1962
[10]. Siochi E. J., Anders J. B., Jr., Cox D. E., Jegley D. C., Fox R. L. and Katzberg S. J., Biomimetics for NASA Langley Research Center, NASA TM-2002-211445, February, 2002
[11]. M.H. Love, P.S. Zink, R.L. Stroud, D.R. Bye, C. Chase,Impact of Actuation Concepts on Morphing Aircraft Structures,AIAA-2004-1724
[12]. D.H. Lee and T.A. Weisshaar. Aeroelastic Studies on a Folding Wing Configuration, AIAA/ASME/ASCE/AHS/ASC Structures, Strucctural Dynamics & Materials Confer18-21.2005.Austin, Texas
[13]. Liu Y.,Kapania R.K., Gern F.H., and Inman D., Equivalent Plate Modeling of Arbitrary Planform Wings, International Conference on Computational Engeineering and Sciences, ICES 2K, Los Angeles, CA August 21-25, 2000
[14]. Gern F.H., Inman D.J., and Kapania R.K., Strctural and Aeroelastic Modeling of General Planform UCAV Wings With Morphing Airfoils, AIAA-2001-1369
[15]. Lee U., Dynamic Continuum Modeling of Beamlike Space Structurs Using Finite-Element Matrices, AIAA J 1990, 28(4):725-731
[16]. Rakesh K. K. and Youhua L., Static and Vibration Analyses of General Wing Structures Using Equivalent Plate Models, AIAA-2000-1434
[17]. Shiv P. Joshi and Zeb Tidwell ,William A. Crossley, and Sekaripuram Ramakrishnan, Comparison of Morphing Wing Strategies Based Upon Aircraft Performance Impacts, 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Confer 19-22 April 2004, Palm Springs, California
[18]. M. Secanell, A. Suleman, and P. Gamboa, Design of a Morphing Airfoil for a Light Unmanned Aerial Vehicle Using High-Fidelity Aerodynamic Shape Optimization, 46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Confer 18-21 April 2005, Austin, Texas
[19].李毅,王生楠,变体飞行器结构振动特性研究,航空计算技术,2006, 6(30)
[20].余晖编译,美国开始试验“变形”飞机,国际航空杂志,2005, (7):50-51
[21]. J.S. Bae, T. M. Seigler, D. J. Inman, Aerodynamic and Aeroelastic Considerations of A Variable-Span Morphing Wing, 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Confer,19-22 April 2004
[22]. Frank S. III a?nd Donald F., Morphing Concept for Quiet Supersonic Jet Boom Mitigation, 43rd AIAA Aerospace Sciences Meeting and Exhibit10-13 January 2005
[23]. K.Hennighausen, G.V.Candler and G.O.Einarsson, Numerical Simulation of Three Dimensional Flow about a Micro Air Vehicle, AIAA-00-3315
[24].张永新,能像鸟一样飞行-NASA研究柔性变形翼飞机,国际航空, 2001,PP60
[25]. David C., Tim S., Frank U.and Matt Mac K., Morphing Inflatable Wing Development for Compact Package Unmanned Aerial Vehicles, 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Confer 19-22 April 2004
[26]. Lian Y.H., Wei Shyy and Raphael H., Shape optimization of a membrane wing for micro air vehicles, 41st Aerospace Sciences Meeting and Exhibit 6-9 January 2003
[27]. Chul H.H., Matthew C., Jou-Young C. and Dimitri N. M., Flexible Multi-Body Design of a Morphing UCAV, AIAA 3rd "Unmanned Unlimited" Technical Conference, Workshop and Exhibit 20-23 September 2004
[28]. V.Venkatakrishnan, A Perspective on Unstructured Grid Flow Solvers, AIAA-95-0067
[29]. Richardson L.F., Weather prediction by numerical process, Cambridge: Cambridge University Press, 1922
[30]. T. J. Baker, Mesh generation:art or science? Progress in Aerospace Sciences,2005, 41: 29-63
[31]. E. Steinthorson, M.S. Liou and L.A. Povinelli, Development of an Explicit Multiblock/ Multigrid Flow Solver for Viscous Flows in Complex Geometries, AIAA-93-2380
[32]. C. Sheng, L. Taylor and D. Whitfield, Multiblock Multigrid Solution of Three Dimensional Incompressible Turbulent Flows about Appended Submarine Configuration, AIAA-95-0203
[33]. J. Flores et.al, Transonic Navier-Stokes solution for a fighter-like configuration, AIAA-87-0032
[34]. R.L. Sorenson, Three-dimensional elliptic grid generation for an F-16. In: Steger JL and Thompson JF, Editors. Three-dimensional grid generation for complex configurations: recent progress(AGAR Dograph No.309, 1988)
[35]. J.L. Steger, F.C. Dougherty and J.A. Benek, A Chimera Grid Scheme, Advances in Grid Generation, Presented at Applied Mechanics, Bioengineering and Fluids Engineering Conference, Houston, Texas, Jun. 20-22, 1983
[36]. D.G. Peace et.al., Development of a Large Chimera Grid System for Space Shuttle, AIAA-93-0533
[37]. Merry M.A. and Shephard M.S., A modified quadtree approach to finite element mesh generation, IEEE Computer Graphics Applications, 1983, 3(1):39-46
[38]. Weatherill N.P., Adaptive mesh generation and flow solution method for Euler equations on unstructured grids,J Comput Phys, 1994, 110:23-38
[39]. Parikh P., Pirzadth S and Lohner R.A package for 3D unstructured grid generation, finite-element flow solution and flow field visualization, NACA CR-182090, 1990
[40]. Kallinderis Y., Khawaja A., and McMorris H. Hybrid prismatic/tetrahedral grid generation for complex geometries, AIAA J 1996, 34:291-298
[41]. Coirier W.J. and Jorgenson P.C.E., A mixed volume grid approach for the Euler and Navier-Stokes equations, AIAA -96-0762
[42]. Karman S.L., SPLITFLOW: a 3D unstructured Cartesian/prismatic grid CFD code for complete geometries, AIAA-95-0343
[43]. Wang Z.J., A Quadtree-based adaptive Cartesian/Quad grid flow solver for Navier-Stokes equations, Computers & Fluids, 1998, 27(4): 529-549
[44]. Zhang L.P., Zhang H.X. and Gao S.C., A Cartesian/unstructured hybrid grid solver and its applications to 2D/3D complex inviscid flow fields, Proceedings of the 7th International Symposium on CFD, Edited by Zhuang FG, International Academic Publishers, Sep. 1997, Beijing, China, PP68-83
[45]. Zhang L.P., Yang Y.J. and Zhang H.X., Numerical simulations of 3D inviscid/viscous flow fields on Cartesian/unstructured/prismatic hybrid grids,Proceedings of the 4th Asian CFD Conference, Edited by Zhang HX, University of Electronic Science and Technology of China Press, Sep. 2000, Mianyang, China, PP93-103
[46]. Joe F. Thompson, Bharat K. Soni and Nigel P. Weatherill, HANDBOOK OF GRID GENERATION,CRC Press,1999
[47].张来平,张涵信,高树椿等,用混合网格数值模拟机动弹复杂流场(大会邀请报告),第十届全国计算流体力学会议论文集,2000年9月,四川,绵阳
[48]. Kazuhiro N., Fumiya T., An Intergrid-Boundary Definition Method For Overset Unstructured Grid Approach, AIAA-99-3304
[49]. Rainald L., Dmitri S., Hong L. and Ravi R., Overlapping Unstructured Grids, AIAA-01-0439
[50]. Grübe B. and Carstens V., Computational of unsteady transonic flow in harmonically oscillating turbine cascades taking into account viscous effects. ASME Journal of Turbomachinery 1998, 120(1):104-111
[51]. Schulze S., Transonic aeroelatic simulation of a flexible wing section, AGARD structures and Materials Panel workshop on Numerical Unsteady Aerodynamics and Aeroelastics Simulation. AGARD R-822, Aalborg, Denmark, 13-17 October 1997, 10:1-20
[52]. Frederic J.B., Consideration on the spring analogy. International Journal for Numerical Methods in Fluids 2000, 32:647-668
[53]. Batina J.T., Unsteady Euler airfoil using unstructured dynamic meshes, AIAA Journal 1990, 28(8):1381-1388
[54]. Slikkeveer P.J., van Loohuizen E.P. and O’Brien SBG, An implicit surface tension algorithm for Picard solvers of surface-tension-dominated free and moving boundary problems. International Journal for Numerical Methods in Fluids 1996, 22:851-865
[55]. Hassan O., Probert E.J. and Morgan K., Unstructured mesh procedures for the simulation of three-dimensional transient compressible inviscid flows with moving boundary components, International Journal for Numerical Methods in Fluids 1998, 27:41-55
[56]. Blom F.J., Leyland P., Analysis of fluid-structure interaction on moving airfoils by means of an improved ALE method, AIAA -97-1770
[57]. Farhat C., Lesoinne M. and Maman N., Mixed explicit/implicit time integration of coupled aeroelastic problems: hree-field formulation, geometric conservation and distributed solution, International Journal for Numerical Methods in Fluids1995, 21:807-835
[58]. Piperno S., Explicit/implicit fluid/structured staggered procedures with a structural prediction and fluid subcycling for 2D inviscid aeroelastic simulation,International Journal for Numerical Methods in Fluids 1997, 25:1207-1226
[59].张来平,王振亚,杨永健,复杂外形的动态混合网格生成方法,空气动力学学报,2004, 22(2):231-236
[60].常兴华,张来平,张涵信,柔性变形体动态网格生成方法,第一届近代气动力(热)会议,四川绵阳,2006
[61].常兴华,段旭鹏,王振亚,张来平,张涵信,基于动态混合网格的鱼体巡游数值模拟,全国第十一届分离流、旋涡和流动控制会议,湖南,张家界,2006
[62]. Nakahashi K., Deiwert G.S., Three-dimensional adaptive grid. AIAA Journal 1986;24(6):948-954
[63]. Liu X.Q., Qin N. and Xia H., Fast dynamic grid deformation based on Delaunay graph mapping. J. Comput. Phys., 2006, 211:405-42
[64].郭正、瞿章华、刘君,用非结构动网格技术计算确定再入飞行器配平攻角,空气动力学学报,2004, 22(3):225-258
[65].杨国伟,钱卫,飞行器跨声速气动弹性数值模拟,力学学报,2005, 37(6):769-776
[66].肖涵山,自适应笛卡尔网格在三维复杂外形流场及机弹分离数值模拟中的应用,博士学位论文,中国空气动力研究与发展中心,2005
[67].张来平,王志坚,张涵信,动态混合网格生成及隐式非定常计算方法,力学学报,2004, 36(6)
[68]. Zhang L. P., Wang Z. Y., Yang Y. J. and Zhang H. X.,3D Dynamic Hybrid Grid Generation and Implicit Solver for Unsteady Flows,The 6th Asian Computational Fluid Dynamics Conference, TAIWAN, OCTOBER 24~ OCTOBER 27, 2005
[69]. Zhang L. P., Chang X. H., Duan X. P., Zhang H. X.,A Block LU-SGS Implicit Dual Time-Stepping Algorithm on Hybrid Dynamic Meshes for Bio-fluid Simulations. Computers & Fluids
[70]. Jameson A., Time Dependent Calculations Using Multigrid with Application to Unsteady Flows past Airfoils and Wings. AIAA-91-1596
[71]. Luo H., Baum J.D., Lohner R., A Fast, Matrix Free Implicit Method for Compressible Flows on Unstructured Grid. J. Comput. Phys., 1998, 146:664~690
[72]. Choi Y.H. and Merkle C.L., Dual-Time Stepping and Preconditioning for Unsteady Computations. AIAA-95-0078
[73]. Van Leer B. and Mulder W.A., Relaxation methods for hyperbolic conservation laws.In:Proceedings of the INRIA Workshop on Numerical Methods for the Euler Equations of Fluid Dynamics, Rocquencourt, France, December, 1983
[74]. Chen R.F. and Wang Z.J.,Fast,block lower-upper Gauss Seidel scheme for arbitrary grids, AIAA J 2000, 38(12):2238-45
[75]. Thomas,P.D. and Lombard C.K.,The Geometric Conservation Law-A Link Between Finite-Difference and Finite-Volume Methods of Flow Computation onMoving Grids,AIAA-78-1208
[76]. Michel L. and Charbel F., Geometric Conservation Laws for Flow Problems with Moving Boundaries and Deformable Meshes, and their Impact on Aeroelastic Computations. Comput. Methods. Appl. Mech. Engrg., 1996, 134:71-90
[77].张来平,非结构网格、矩形/非结构混合网格复杂无粘流场的数值模拟,博士学位论文,中国空气动力研究与发展中心,1996
[78].徐强勋,三维Delaunay四面体网格生成加密技术与应用,上海交通大学硕士学位论文,2002
[79].张来平,杨永健等,三维复杂外形的空间非结构网格的自动生成与应用,计算物理,1999, 16(4)
[80].张来平,张涵信,高树椿,矩形/非结构混合网格技术及在二维/三维复杂无粘流场模拟中的应用,空气动力学学报,1998, 16(1):77-88
[81]. P.R. Spalart and S.R. Allmaras, A one equation turbulence model for aerodynamic flows, AIAA-92-0439
[82]. D. Wilcox, Turbulence modeling for CFD, second edition, DCW Industries, Inc., 1998
[83]. Harris C.D.,Two-Dimensional Aerodynamic Characteristics of the NACA0012 Airfoil in the Langley 8-foot Transonic Pressure Tunnel, NASA TM81927
[84].宋书恒,微型飞行器相关气动力特性的数值模拟,中国空气动力研究与发展中心博士论文,2007
[85]. Christopher L. Rumsery, Thomas B.Gatski, Susan X.Ying and Arild Bertelrud, Prediction of high-lift Flows Using Turbulent Closer Models, AIAA-97-2260
[86]. Mehmet S. and Lakshmi N.S., Dynamic stall alleviation using a deformable leading edge concept-anumerical study.AIAA-2000-0520
[87]. V.J.Modi, On the moving surface boundary-layer control, AIAA-2000-2238
[88]. Razvan F. and Brian E.Wake, Parametric analysis of directed-synthetic jets for improved dynamic-stall performance, AIAA-2003-216
[89]. Mehmet S., Lakshmi N.S., M.S.Chandrasekhara and Chee Tung, Dynamic stall alleviation using a deformable leading edge concept-a numerical study, AIAA-2000-0520
[90]. Bo-sung Lee, Kwanjung Yee, Wandon Joo and Dong-Hu Lee, Passive control of dynamic stall via nose droop with Gurney flap.AIAA-2005-1364
[91].吴介之,马晖扬等,涡动力学引论,高等教育出版社,1993
[92].倪云华,杨爱明,翁培奋,低雷诺数下微型飞行器主动变形翼非定常气动特性数值分析,全国流体力学青年讨论会论文集,四川绵阳,2005,PP125-130
[93].吴锤结,解妍琼,吴介一,“流体滚动轴承”效应及其在流动控制中的应用,全国第九届分离流、漩涡和流动控制会议文集,北京,2002, PP281-290