0 引言

为提高飞机的作战能力，对航空发动机的推重比要求越来越高，矢量喷管作为航空发动机的主要组成部分也面临着高效、轻质的要求[1-2]。传统机械调节矢量喷管结构复杂、质量轻，无法满足高推重比的要求，固定几何气动矢量喷管采用二次流对主流的干扰形成矢量偏转，结构简单、质量轻，可以满足未来发动机的需求[3-4]。

党的十八大以来，以习近平同志为核心的党中央坚持把解决好“三农”问题作为全党工作重中之重，坚持农业农村优先发展，两个“坚持”掷地有声、落地见效，“三农”发展开创出崭新局面。综观十八大以来的6个中央“一号文件”，逻辑上一脉相承、层层递进、浑然一体，实践中加速了城乡融合发展的历史进程。乡村振兴战略的提出，更是奏响了加快实现农业强、农村美、农民富的气势磅礴的主旋律，书写了我国城乡社会向现代化转型过程中致力于协调发展、共同繁荣的浓墨重彩的新篇章。

国外的研究机构开展了各种流体推力矢量控制方式的研究，综合起来主要集中在3种控制方法上，即激波矢量控制技术、喷管喉部偏移技术、反流控制技术。Giuliano等在NASA兰利喷流排气试验装置上专门做了基于激波矢量控制的2元球面收敛/扩散调节片俯仰推力矢量喷管试验[5]；Deere等在流体偏航矢量喷管的基础上，开展了喉部位置偏移法实现推力矢量的数值计算[6-7]；Leavit等对多轴推力矢量喷管的吼道倾斜方法进行了稳态研究[8]；Kenrick等对流体推力矢量喷管的多喷射点方法进行了试验和数值计算[9]；NASA兰利研究中心进行了大尺寸的反流法实现推力矢量的试验研究，Hunter等采用PAB3D计算软件进行了反流法计算[10]。国内在流体推力矢量技术方面也开展了一定的研究工作。乔渭阳等采用试验和数值计算相结合的方法，研究了二次流喷射对流体推力矢量的影响[11]；罗静等采用数值模拟的方法，计算分析了流体喷射对喷管气动矢量角的影响[12]；王占学等完成了基于二次流喷射的流体推力矢量数值计算和实验方案[13-15]。

目前固定几何气动矢量喷管二次流效率低，在现有的二次流流量下，矢量偏转角较小；在非设计点状态下的低状态（如主流落压比3）时喷管推力性能较低；对二次流控制喉道面积的研究很少。为此，本文主要研究了阀门结构的二次流注射方式（包括喉道和扩张段二次流）、二次流气动参数、主流气动参数和引射外界气体对流场结构与喷管性能（推力系数、矢量角和喉道面积控制率）的影响。

1 物理模型

固定几何气动矢量喷管的物理模型如图1所示。主要包括收敛段、扩张段、二次流通道和引射通道。喉道二次流和扩张段二次流为阀门结构的注射方式如图1（a）所示，二次流出口位于喷管内，根据不同状态的需求，通过转动阀门调节二次流的关闭或开启角度。喉道二次流的角度β为60°和74.5°2种情况，扩张段二次流的角度α为30°和60°2种情况,二次流角度指二次流和扩张段的夹角；β＝0°时喉道二次流阀门关闭，α＝0°时扩张段二次流阀门关闭。二次流为常规的注射方式如图1（b）所示，非阀门结构，出口位于扩张段壁面位置，并垂直于扩张段。

本刊主要栏目有：【1】专家论坛，专家就自身研究方向的热点问题发表相应的看法、观点和态度；【2】专家述评，组织报道介绍分析本学科领域科学技术研究现况及发展水平的专家述评；【3】专题报道，集中报道本学科领域的热点问题的专题研究论文；【4】论著，报道本学科领域前沿性基础和临床研究的最新研究成果；【5】临床研究，主要报道本学科领域临床诊疗的新经验、新技术；【6】病例报告，报道本学科领域的典型性临床病例，包括罕见病例；【7】文献综述，报道反映本学科领域研究进展和动态的国内外文献综述。

 

  

图1 固定几何气动矢量喷管物理模型

2 计算方法

2.1 流场计算

运用商业软件进行全3维流场计算，热燃气假设为理想可压缩流体，使用基于密度耦合算法，湍流模型采用SST k-ω模型，主流边界条件设置为压力入口，总温为1100 K，总压为0.3、0.4 MPa；二次流设置为压力入口，总温为288 K，总压为0.3～0.6 MPa；压力出口温度为288 K，总压为0.1 MPa。

涂料生产使用的易挥发有机溶剂引起的VOCs排放的比例最高，占全部排放量的80%以上,其中二甲苯的含量最高。通过采用密闭型设备、提高自动化水平、优化生产工艺，可以有效地从源头上控制VOCs的产生量。对于生产过程中排放出的VOCs可以采用有效的净化处理工艺进行处理。

针对图1（b）所示的几何模型，研究了主流落压比（N PR=3、4）和二次流落压比（S PR＝3、4、5、6）对矢量角的影响，计算结果如图10所示。从图中可见，扩张段二次流注射方式非阀门结构时产生的最大矢量角约为9°；当主流压力不变时，随着二次流压力增加喷管矢量角增大，最大增加约为3°；当二次流压力不变时，随着主流压力增加矢量角约减小3°～4°。

  

图2 对称面网格分布

  

图3 喷管区域网格分布

2.2 轴向推力系数和矢量角计算

（3）政府人员定位偏差引发执行偏差。一是认知偏差。即政府工作人员对政府购买公共服务的政策及法律缺乏正确认识。对购买公共服务的错误定位导致政府部门对其缺乏应有的积极性和主动性，甚至演变为政绩工程。二是理念滞后。政府工作人员在政府由管理走向治理的机构改革与职能变革中，对社会组织缺乏信任。治理理念的滞后不利于合作型政府的建设，严重挫伤社会组织人员积极性，引发政社矛盾。三是角色错位。即认为政府购买公共服务在服务外包的同时，责任和困难也随之推卸给了社会组织，在公共服务供给时出现不作为现象。

N PR和S PR对流场马赫数的影响如图6所示。从图中可见，在N PR一定的情况下，随着S PR增大，二次流对主流造成的斜激波位置会向前移动，说明S PR增大对主流具有更强的扰动作用，使得喷管轴向推力系数减小；在S PR一定的情况下，随着N PR增大，二次流对主流造成的斜激波位置会向后移动，并且扩张段后段低速流动区域减小，说明N PR增大会减弱二次流对主流的扰动作用，并且更接近设计状态，减少过膨胀区域，使得喷管轴向推力系数增大。

 

扩张段二次流角度对流场马赫数的影响如图12所示。从图中可见，随着扩张段二次流角度增大，激波位置变化不太明显，尾喷流偏转角度略有减小，说明扩张段二次流角度对矢量角影响不大。另外，根据图11、12对比可知，注射方式为阀门结构的扩张段二次流产生的矢量偏转效应明显优于注射方式为非阀门结构的扩张段二次流。

橡胶止水带及接缝填料按设计位置安装，并用钢筋固定牢固，预防止水带变形，跑偏和撕裂，保证止水带周围特别是止水带下面的混凝土振捣密实，避免止水带周围形成空穴。

 

式中：U i为等熵完全膨胀时喷管出口速度

 

矢量角δ的大小反映了纵向推力Fy与轴向推力Fx比值的大小，定义为

 

喉道控制率R TAC为喉道面积最大值与最小值之差与最小面积的比值

 

3 计算方法验证

根据数值计算模型，进行了冷态吹风试验，其中N PR=4时，中间非偏转状态的静压数值计算结果与试验结果的对比如图4所示。从图中可见计算结果与试验结果基本吻合，验证了数值计算方法的可行性。

  

图4 壁面中心线静压分布

4 计算结果与分析

4.1 推力系数的影响因素

针对β＝0°、α=60°的几何模型，研究了主流落压比（N PR=3、4）和扩张段二次流落压比（S PR＝3、4、5、6）对推力系数的影响；针对β＝0°的几何模型，在N PR＝4的情况下，研究了α对轴向推力系数的影响；针对 α=30°、β＝74.5°的几何模型，在 N PR＝3、S PR=6的情况下，研究了引射方式对喷管性能的影响；计算结果如图5所示。从图中可见，N PR不变，随着S PR增大喷管轴向推力系数减小，S PR从3增加到6时，轴向推力系数减小约3%～4%；S PR不变，随着N PR增大轴向推力系数增大，N PR从3增加到4，轴向推力系数减小约5%；二次流压力不变，二次流角度α增加，轴向推力系数最大减小1.5%；引射开启后的轴向推力系数比引射关闭时最大提高5%。

  

图5 N PR、S PR、α和injection对推力系数的影响

The paper by Libertino et al is important for readers to appreciate how to identify those patients that are likely to have a good response to revascularization.

 

  

图6 N PR和S PR对流场马赫数的影响

引射对流场马赫数及流场静压的影响分别如图7、8所示。从图中可见，在N PR=3时喷管处于过膨胀状态，扩张段出现较大面积的负压区和低速区，在引射关闭时喷管末端出现外界大气被引到喷管内部的现象，出现反流，降低喷管推力性能；在引射开启后，引射流与主流混合提高喷管流量，增加喷管动量，减少主流区域的负压区，同时减少末端气体引入导致的反流现象。

喷管的推力为出口气流动量与压差在出口面上积分的和，X方向的实际推力F x为式中：m˙为实际质量流量；Ux为出口截面的X方向的速度；P e为出口截面的静压；A e为出口截面的面积。

  

图7 引射对流场马赫数的影响

扩张段二次流角度对流场马赫数的影响如图9所示。从图中可见，随着二次流角度α增大，二次流对主流造成的斜激波位置会向前移动，说明α增大会加强二次流对主流的扰动作用，使得喷管轴向推力系数减小。

  

图8 引射对流场静压的影响

  

图9 扩张段二次流角度对流场马赫数的影响

4.2 矢量角的影响因素

采用结构化网格，总数控制在190万左右。整个计算域中心对称面网格分布如图2所示，喷管区域对称面网格如图3所示。在内流道对网格进行加密，外场网格逐渐变稀。

  

图10 N PR、S PR和α对矢量角的影响

主流压力和二次流压力对流场马赫数的影响图14所示。从图中可见，当N PR不变时，随着S PR增大，二次流出口马赫数增大，二次流后面的低速区增加，主流高速区向中心线方向偏移，喷管有效流通面积减小，喷管喉道控制率提高；当S PR不变时，随着N PR增大，二次流出口马赫数减小，二次流后面的低速区减小，主流高速区向壁面方向偏移，喷管有效流通面积增大，喷管喉道控制率降低。

N PR和S PR对流场马赫数的影响如图11所示。从图中可见，当主流压力不变时，随着二次流压力增加，激波位置向下移动，矢量角增大，说明二次流压力增加对主流影响增大；二次流压力不变，随着主流压力增加，激波位置向上移动，矢量角减小，说明主流压力增加需要相应增加二次流压力才能满足矢量角的需求。

  

图11 N PR和S PR对流场马赫数的影响

定义轴向推力系数C f x为X方向实际推力Fx与总的理想推力F i（主流与二次流理想推力之和）的比值

  

图12 扩张段二次流角度对流场马赫数的影响

4.3 喉道控制率的影响因素

针对α＝0、β=74.5°的几何模型，研究了主流落压比（N PR=3、4）和喉道二次流落压比（S PR＝3、4、5、6）对喉道控制率的影响，计算结果如图13所示。从图中可见，当主流压力不变时，随着喉道二次流压力增加，喉道控制率增大，最大增加约18%；当二次流压力不变时，随着主流压力增加，喉道控制率减小约2%～11%。

  

图13 N PR、S PR和β对喉道控制率的影响

针对β＝0°的几何模型，在主流落压比N PR＝3的情况下，研究了扩张段二次流角度（α=30°、60°）对矢量角的影响，计算结果如图10所示。从图中可见，扩张段二次流注射方式为阀门结构时产生的最大矢量角为17°，明显高于扩张段二次流注射方式为非阀门结构产生的最大矢量角；随着扩张段二次流角度增加，矢量角略有减小约为0.5°。

针对α=60°的几何模型，在N PR＝3、扩张段S PR=3的情况下，研究了喉道两侧二次流角度（β=60°、74.5°）对喉道控制率的影响，计算结果如图14所示。从图中可见，随着喉道二次流角度增大，喉道控制率增加约为5%～20%，最大为0.59。

题名“志归”，所记辞官归里，亦以示归田之志。望着“平畴远风，绿畦如浪”的家园景色，祁彪佳亦觞亦咏，“乃此是志吾之归也，亦曰归固吾志也”[1]160。志归斋有一种平淡冲远、古拙浑厚之美，是一种典型的诗意的栖居地。这种诗意，是以陶渊明为代表的田园诗意，它以牧歌般的旋律缭绕在祁彪佳的心胸，催生着他挥之不去的“归去来”绵绵情思。

 

  

图14 N PR和S PR对流场马赫数的影响

喉道二次流角度对流场马赫数的影响如图15所示。从图中可见，随着喉道二次流角度增大，二次流出口马赫数增大，二次流后面的低速区增加，主流高速区减小，喷管有效流通面积减小，喷管喉道控制率提高。

  

图15 扩张段二次流角度对流场马赫数的影响

5 结论

本文通过对主流、喉道二次流、扩张段二次流、引射等各因素对固定几何气动矢量喷管推力系数、推力矢量角和喉道控制率等气动性能的影响分析，得到以下结论。

（1）在本文的研究范围内，主流落压比由3增加到4时轴向推力系数增加约为5%；扩张段二次流落压比由3增加到6时轴向推力系数约减小3%～4%；主流落压比为3时，引射开启后的轴向推力系数比引射关闭时最大提高5%；扩张段二次流角度由30°增加到60°时轴向推力系数最大减小1.5%；

（2）主流落压比由3增加到4时矢量角减小约3°～4°；扩张段二次流落压比由3增加到6时喷管矢量角增大，最大增加约3°；随着扩张段二次流角度增大，矢量角略有减小；注射方式为阀门结构的扩张段二次流产生的最大矢量角比注射方式为非阀门结构的扩张段二次流约大8°；

（3）喉道二次流落压比由3增加到6时喉道控制率最大增大约18%；主流落压比由3增加到4时喉道控制率减小约2%～11%；喉道二次流角度由60°增大到74.5°时，喉道控制率提高约5%～20%。

引射是1种有效提高低状态（如N PR=3）下固定几何气动矢量喷管推力系数的方式；扩张段二次流注射方式为阀门结构时能有效提高喷管矢量偏转角；喉道二次流注射方式为阀门结构时能有效控制喷管喉道面积；引射和阀门结构的二次流注射方式的提出为固定几何气动矢量喷管的工程化应用提供了新方向。

因此，基于以上分析可以认为，在稳定增长的目标下，如果耐用品的寿命足够长，货币政策当局应选择耐用品部门通胀作为货币政策盯住目标，无论是基于缓和产出波动的角度，还是基于降低政策引致的社会福利损失的角度，均应如此。

参考文献：

[1]贾东兵,周吉利,邓洪伟.固定几何气动矢量喷管综述[J].航空发动机,2012,38（6）:29-33.JIA Dongbing,ZHOU Jili,DENG Hongwei.Summary of fluidic control fixed geometry nozzle technology[J].Aeroengine,2012,38（6）:29-33.（in Chinese）

[2]靳宝林，郑永成.一种有前途的推力矢量技术-流体推力矢量控制喷管[J].航空发动机，2000（4）：52-57.JIN Baolin,ZHENGYongcheng.One promisingthrust vectoring technology-fluidic control vectoring nozzle[J].Aeroengine,2000（4）:52-57.（in Chinese）

[3]Reginald G W.Fluidic thrust vectoring and throat control exhaust nozzle[R].AIAA-2002-4060.

[4]Abeyounis W K，Bennett B D.Static internal performance of an over-expanded,fixed-geometry,nonaxi symmetric nozzle with fluidic pitch-thrust-vectoring capability[R].NASA-TP-3645,1987.

[5]Wing V J.Static investigation of a fixed-aperture exhaust nozzle employing fluidic injection for multiaxis thrust vector control[R].AIAA-97-3149.

[6]Deere K A.Computational investigation of the aerodynamic effects on fluidic thrust vectoring[R].AIAA-2000-3598.

[7]Deere K A，Berrier B L.Computational study of fluidic thrust vectoring using separation control in a nozzle[R].AIAA-2003-3803.

[8]Leavitt L D，Lamb M.Static investigation of a multiaxis thrust-vectoring nozzle concept featuring diagonal throat hinge [R].NASA-TP-3492，1994.

[9]Kenrick A.Experimental and computational investigation of multiple injection ports in a convergent-divergent nozzle for fluidic thrust vectoring[R].AIAA-2003-3802.

[10]Hunter CA,Deere K A.Computational investigation of fluidic counter flow thrust vectoring[R].AIAA-99-2669.

[11]乔渭阳,蔡元虎.基于次流喷射控制推力矢量喷管的实验及数值研究[J].航空动力学报，2001（3）：273-278.QIAO Weiyang,CAI Yuanhu.A study on the two-dimensional thrust vectoring nozzle with secondary flow injection[J].Journal of Aerospace Power,2001，15（3）：273-278.（in Chinese）

[12]罗静，王强，额日其太.两种流体控制方案矢量喷管内流场计算及分析[J].北京航空航天大学学报，2004，30（7）：597-601.LUO Jing,WANG Qiang,Eriqitai.Computational analysis of two fluidic thrust vectoring concepts on nozzle flow field[J].Journal of Beijing University of Aeronautics and Astronautics,2004，30（7）：597-601.（in Chinese）

[13]李志杰，王占学，蔡元虎.二次流喷射位置对流体推力矢量喷管气动性能影响的数值模拟 [J].航空动力学报，2008，23（9）:1603-1608.LI Zhjie,WANG Zhanxue,CAI Yuanhu.Numerical simulation of aerodynamic effects of secondary injection positions on fluidic thrust vectoring nozzle [J].Journal of Aerospace Power,2008，23（9）:1603-1608.（in Chinese）

[14]王占学，王玉男，李志杰，等.基于激波控制的流体推力矢量喷管试验[J].推进技术，2010，31（6）:751-756.WANG Zhanxue，WANG Yunan，LI Zhijie，et al.Experiment on fluidic thrust-vectoring nozzle based on shock control concept[J].Journal of Propulsion Technology,2010，31（6）:751-756.（in Chinese）

[15]刘爱华，王占学.二次流喷射对喷管流场性能的影响[J].推进技术，2007，28（2）:144-147.LIU Aihua,WANG Zhanxue.Effect of secondary flow injection on flow field and performance of nozzle[J].Journal of Propulsion Technology,2007，28（2）:144-147.（in Chinese）