半浸桨不同半径切面入水的水动力特性

半浸桨工作时的特殊工作环境及其水动力特性的研究一直是研究工作中的难点，即便是试验研究，也由于试验精度限制、尺度效应、边界干扰等不利因素的存在，并不能将半浸桨的敞水性能测试得十分准确.Seyyedi等在爱尔兰进行了841-B半浸桨的敞水试验，得到的试验结果与Olofsson试验的推力和转矩系数进行比对，尽管二者随进速变化的趋势有较好的一致性，但是推力和转矩系数在高进速时存在较大的误差.随着近几年计算流体力学技术的发展，国外学者Yari等以及国内学者任振等采用数值模拟的方法，建立了较为完整的数值模拟手段，分别研究了浸深比、斜流角、双桨耦合效应等对半浸桨水动力特性的影响机理以及半浸桨工作过程中的流固耦合机理，为半浸桨在实际工程中的应用提供了技术指导.

式中Y表示输出变量，K 和 L分别代表资本与劳动投入变量(输入变量)，t表示时期，而A(t)表示TFP在时期t的水平。

但是在对半浸桨进行水动力特性预报的过程中，一方面存在计算量较大的问题；另一方面，数值模拟的过程中，预报值的推力系数偏小，转矩值偏大，使得数值模拟方法获得的效率值偏小.因此，近些年Nasrin等采用数值模拟的方法，对半浸桨的杯型切面的入水过程进行水动力数值模拟，验证了随着进速系数的不断增加，半浸桨杯型切面的通气腔会不断向随边方向移动的现象.但是，这些学者在模拟半径为的半浸桨杯型切面入水的过程中，所采用的剖面一般为半浸桨的0.7无因次半径处的有效剖面.然而，随着半浸桨杯型切面的无因次半径位置选择的不同，不仅杯型切面的几何结构形式会产生差异，其垂向的入水速度也会产生较大差异.因此，还需要进一步研究半浸桨杯型切面的切面位置对其入水过程的影响.

亚兰阿姨这番“人、树之喻”，似乎是我所见闻过的评人之论中，最新颖独特的了。也希望当今之人，能对上代知识文化者的经验见地有所借鉴。

本文采用RANS方法，并结合标准-湍流模型、流体体积(VOF)方法以及重叠网格技术，对半径为的半浸桨在06、07及08处的半浸桨杯型切面入水过程进行数值模拟，最终得到不同工况下的水动力特性、自由液面形式、通气腔形式、流场形式以及表面压力分布等，并分析了半浸桨切面位置对其水动力特性的影响机理.

1 计算前处理

1.1 模型参数

选取小底升角楔形体以30 m/s的速度入水的工况进行数值模拟方法和半解析方法之间的验证.在划分好与1.3节相同的网格模型后，监测如图4所示的点和′点压力值的时历曲线，并将数值模拟结果与半解析方法、Oger等的数值模拟方法进行比对，结果如图5所示，图中′为测点所受压强.

1.2 计算工况设置

本文计算的半浸桨杯型切面入水过程中的进速系数、横向受力系数、垂向受力系数以及效率的定义如下：

(1)

(2)

(3)

(4)

式中：为进速；为转速；为直径；、分别为二维杯型切面的向、向受力；为水的密度；为25网格的向计算域厚度；和分别为进速和切向速度；和分别为三维半浸桨受力的推力系数和转矩系数.

根据已有的计算经验可知，半浸桨工作状态的转速大小为 2 254 r/min，进速系数的范围为0.7～1.3.因此，通过换算，可以计算出0.6、07及08无因次半径处的半浸桨杯型切面入水过程中的垂向速度分别为17.69、20.64 及23.59 m/s.在数值模拟的过程中，通过不断改变来流速度的大小，可以改变进速系数的大小.

鄱阳湖区二期防洪工程是在鄱阳湖治理一期工程的基础上，为进一步解决鄱阳湖区防洪问题而进行的一项工程。从1998年开始，先后实施了四个单项工程、第五个单项工程和第六个单项工程，简述如下：

1.3 计算模型设置

..计算域网格划分 本文采用2.5的方式建立计算网格，即建立一种垂向尺度较小的三维六面体形式的网格，由于重叠网格区域在数值模拟的过程中会垂向运动且该物理模型存在自由液面，所以需要对重叠网格区域、重叠网格运动域以及自由液面处的网格尺度进行加密，网格的划分方式如图2所示.

当选取的半浸桨杯型剖面位置不同时，一方面杯型剖面的结构形式不同，尤其是随边位置的差异；另一方面垂直入水速度的大小之间的差异，会使得叶剖面入水时的流场特性产生一定差异，而流场特性的差异一般会影响杯型剖面表面的压力分布，进而影响杯型剖面入水过程的敞水特性以及通气特性.因此，对其流场特性的研究十分重要.不同位置的杯型剖面在不同通气状态时的速度云图如图10～12所示，图中’为速度幅值的大小.

护理人员要密切监测患者病情变化,对患者的疼痛部位、持续时间及缓解方式做好详细记录,必要时24h监测患者血压、心率等生命体征。用药过程注意药物疗效及不良反应,对于冠心病合并心绞痛患者而言,其需要长期服药才能控制病情,而老年患者因生理机能退化,记忆能力下降,服药依从性差,常会产生不按医嘱服药,不配合治疗等情况。因此,护士应该耐心、仔细的将冠心病合并心绞痛的治疗机制讲予患者听,让其了解用药的重要性,提高患者用药依从性,提高患者治疗效果。

1.4 数值模拟的验证

选定的计算模型是半径为的841-B型螺旋桨在0.6、07、08位置处的杯型剖面.图1为半浸桨模型和不同切面对应位置的杯型剖面形式，其中为无因次半径；表1为841-B型半浸桨的具体模型参数.

观察图5中楔形体高速入水的抨击压力曲线可以发现，相较于其他已有的数值模拟方法，现有数值模拟方法能够获得较为精准的结果.在不同测点的压力值比对方面，本文采用的数值模拟手段同半解析法得到的结论十分近似，尤其在<0.001 5 s时，无论是压力还是压力演变趋势，数值模拟的结果均与半解析解相差无几.但是随着楔形体入水深度的不断增加，自由液面附近液滴的飞溅状态以及自由液面的形式会发生较为剧烈的变化，进而使得不同数值模拟方式及半解析解法的精度有所下降.但从受力趋势来看，本文数值模拟所得结果优于其他已有数值模拟方法.

2 计算结果

2.1 自由液面形式及通气形式的变化

为了得到不同切面位置的半浸桨杯型切面入水过程中气腔随进速大小的发展趋势，按照图6的方式提取了通气腔的长度和厚度，得到不同切面位置下通气腔的无因次长度()和无因次厚度()随进速系数的变化趋势，如图7所示.除此之外，将0° 斜流角工况下的不同无因次半径位置处杯型剖面入水过程的自由液面形式进行比对，对比结果如图8所示.

由图3可知，相较于网格1的计算结果，另外两种网格在数值模拟的过程中有非常相近的受力结果.因此为了提高计算效率，本文选择网格数为 824 416 的网格2.

除此之外，由图9所示，随着杯型剖面位置选择的不同，即越接近叶梢位置，横向力系数和效率会越大，垂向力系数相差不大.这是由于半浸桨运行时的转速是一定的，而剖面位置越接近其梢部，线速度就会越大，即杯型剖面垂直入水的速度也越大.而垂直入水的速度越大，会影响到剖面吸力面与压力面两侧的流体流动速度，使得两侧的压力差增大，影响到其横向受力大小，使其横向受力增大.但是由于压力的垂向分量较小，使得剖面位置对垂向受力的影响不大.根据式(4)可知，敞水效率与半浸桨杯型剖面入水过程的横向受力和垂向受力之比成正相关，敞水效率也会随着杯型剖面所选取的位置越接近梢部而增大.但是另一方面，随着进速系数的不断增大，在这几个杯型切面入水均达到部分通气状态时，它们之间横向受力的差距十分微小.

2.2 敞水特性分析

为了比较不同位置的剖面对半浸桨杯型剖面入水过程的敞水特性的影响，对=0.6，0.7，0.8时的横向力系数、垂向力系数以及敞水效率进行比对，如图9所示.

通过图9可知，无论杯型剖面的位置如何选择，其入水过程的横向力系数、垂向力系数都会在过渡状态下产生比较大的波动，并在部分通气状态后随着进速系数的增加而降低；而效率曲线则都在较大的进速系数位置处达到峰值，并在达到峰值后随进速的不断增加而降低.

而根据杯型剖面位置选择的不同，其完全通气状态和部分通气状态之间的过渡状态发生时的进速系数以及效率峰值位置的进速系数会有所差异.具体而言，随着杯形剖面选择的位置越接近叶梢，力的较大波动发生时的进速系数会较大；而0.6杯型剖面的最大效率发生在进速系数为1.1附近，0.7和08杯型剖面的最大效率发生在进速系数为1.2附近.造成该现象的原因如2.1节的分析，即杯型剖面的位置越接近梢部，完全通气状态与部分通气状态之间的过渡状态会更趋向于发生在一个较大的进速系数位置，从而对完全通气状态和部分通气状态的转变形成一个滞后作用，而当杯型剖面达到过渡状态时，其横向力系数及垂向力系数会发生较大的波动.除此之外，不同位置杯型切面入水过程的效率峰值也会随着切面位置接近叶梢而发生一定滞后作用.

HIV属于病毒科慢病毒属中的人类慢病毒组，为直径100～120 nm球形颗粒，由核心和包膜两部分组成。核心由衣壳蛋白（CA，p24）所组成，衣壳内包括两条完全一样的病毒单股正链RNA、核壳蛋白（NC）和病毒复制所必需的酶类，含有反转录酶 （RT，p51/p66）、 整合酶 （IN，p32） 和蛋白酶（PR，p10）。病毒的最外层为包膜，来源于宿主细胞膜的膜质结构，其中嵌有外膜糖蛋白gp120和跨膜糖蛋白gp41；包膜结构之下的是基质蛋白（MA，p17），形成一个病毒内壳。

通过观察图7发现，在进速系数较小时，3种切面位置的杯型切面在入水的情况下，通气腔均有较大的长度和厚度，而随着半浸桨杯型切面入水过程中进速的不断增加，该通气腔的长度和厚度均有所减小，即发生如图8所示的通气腔随着进速的不断提高而向随边方向收缩的现象.总体而言，相对于0.6和08位置的半浸桨的杯型切面而言，0.7位置的半浸桨杯型切面在入水过程中形成的通气腔的厚度以及部分通气状态情况下的通气腔长度均较大，即半浸桨的有效切面通气现象的程度较为严重.从图7(a)可以看出，在完全通气状态过渡到部分通气状态的过程中，通气腔的长度会在某一进速系数下大幅度收缩.除此之外，随着无因次半径的位置不断接近半浸桨的叶梢，该通气腔剧烈收缩发生时的进速系数会有所增大，即通气腔的剧烈收缩过程会有一定程度的滞后.相同的结果可以在图8中显示，如图8(a)所示，当进速系数为0.8，=0.7，0.8时，会达到完全通气状态；而当=0.6时，为部分通气状态.如图8(b)所示，在半浸桨杯型剖面入水过程的进速系数为0.9时，=0.6，0.7处的半浸桨杯型剖面的通气状态已经达到部分通气状态，而=0.8半浸桨杯型剖面处的通气状态则仍处于过渡状态.

2.3 流场特性分析

..网格收敛性分析 为了验证网格的无关性，本文取3种加密形式的网格：网格1，网格数为 634 501；网格2，网格数为 824 416；网格3，网格数为 1 085 275.计算0.7位置的杯型切面在进速系数为0.8时的工况，并比较其横向力和垂向力随时间()变化的曲线，如图3所示.

本研究也有一定的局限性：首先，本研究为回顾性研究，本身存在潜在的偏倚；其次，本研究为单中心、鼻咽癌好发地区收集的病例，这可能影响模型的广泛适用性，因此需要进一步多中心、多地区数据收集验证。

随边压力面附近的流速高速区在完全通气状态下，无论杯型切面的剖面位置如何选取，其面积大小几乎相同，而通过分析完全通气状态下半浸桨杯型切面吸力面附近的通气腔的面积分布可知，覆盖整个吸力面表面的通气腔的面积非常大，流速非常快，可以进一步影响到与其相邻的随边压力面附近的流场分布，增大该区域的面积.因此，在完全通气状态下，无论半浸桨杯型剖面的剖面位置如何选取，其压力面随边附近的高速区面积变化不大.但是随着进速系数的增加，在半浸桨杯型切面的入水过程达到部分通气状态后，该压力面随边的高速区面积随着剖面位置选择的不同会发生相当大的变化.具体而言，在部分通气状态下，随着杯型切面的位置选择越接近叶梢，这个压力面随边附近的高速区会有所缩小.这是由于半浸桨杯型剖面的入水过程达到部分通气状态后，不同切面的通气腔缩小的程度会有所差异，进而影响到通气腔的大小以及自由面的抬升程度.而在过渡状态，气腔在吸力面的分布不稳定，不仅会在压力面随边附近形成高速区，在吸力面的导边处也会形成面积较小的高速区.相较于压力面随边位置的高速区面积变化的差异，吸力面导边附近的低速区面积随杯型剖面的选择位置而产生的差异则较小.

2.4 表面压力分布特性分析

半浸桨不同剖面位置的杯型剖面在入水过程中的流场特性存在相当大的差异，反映到杯型剖面在入水过程中的周围压力分布，进而影响到其压力面和吸力面之间的压力差值.使得不同位置的杯型剖面在入水过程中的横向受力和垂向受力产生了如图9所示的变化趋势.为了进一步研究不同剖面位置的压力分布特性，提取了不同通气状态下的半浸桨杯型切面的压力分布，如图13～15所示.

通过比对图10～12可以发现，尽管杯型切面的剖面位置的选取有所差异，但无论是完全通气状态、部分通气状态还是过渡状态的流场分布特性均具有相似之处，即无论杯型切面的位置如何选择，其导边位置的低速区、随边压力面附近的高速区以及气腔中空气流动的高速区均存在.而在速度场中，杯型切面位置选择的差异主要反映在由于垂直入水时的速度不同而造成的这几个高速区以及低速区的面积大小分布.

另一方面，B公司亦可采取制作经销商评估表的形式，列出可用于风险和绩效评估的项目及每个项目的打分建议，由公司业务人员和区域经理每月针对经销商的表现进行打分上报，同时对于主动暴露风险的业务人员予以一定金额的奖励。通过打分表的指标提醒业务人员，在经销商绩效管理中应该关注哪些点，鼓励业务人员通过趋势分析提前发现经销商在哪些指标上存在潜在风险，进行提前预防和解决。

通过观察图13所示的杯型剖面完全通气状态下的表面压力分布可知，在压力面上，0.6处杯型切面的压力值较小，而在0.7位置，压力面的压力值会有较大程度提升，但是若杯型切面的位置更接近梢部，即0.8的位置，这个压力值的提升程度会大大减小.而在达到了部分通气状态后(见图14)，由于通气腔会有所收缩，所以吸力面的压力分布不会如图13一样均为大气压的值，而是会在导边附近发生不同程度的压力下降现象.对于过渡状态的压力分布，可以通过图15看出，杯型剖面吸力面通气腔的大小有所差异，使得吸力面的压力分布同杯型剖面的位置选取有较大差异.

通过图14可以发现，当杯型切面的入水过程达到部分通气状态后，无论杯型切面的位置如何选择，其导边附近的横向受力会与半浸桨的推力方向相反，进而使得半浸桨的杯型剖面在高进速情况下所产生的充当推力的横向受力减小，也正因如此，在达到部分通气状态后，其横向受力的大小会随着进速的增加而大大减小.尽管不同剖面位置的压力分布有所差异，但是整体来看，充当推力的合力大小和与推力相反的力的大小的差值几乎相同，使得在高进速的情况下，半浸桨杯型切面所受的横向力大小几乎相同.

而从图15所反映出来的过渡状态下半浸桨杯型剖面的压力可以看出，由于不同位置杯型剖面入水过程中的压力面附近的自由液面的抬升程度几乎相同，所以压力面的压力分布在过渡状态下没有较大的差异.但是，过渡状态吸力面上通气腔的不稳定性，一方面使得吸力面的通气腔分布有所差异，导致通气腔的压力值有所差异；另一方面，影响了流场的流动方式，进一步使得水动力特性产生了一定差异，影响吸力面上的压力分布.

3 结论

采用数值模拟的方法，对不同无因次半径切面位置半浸桨杯型切面的入水过程进行了有关水动力特性以及通气特性的分析，通过对比不同无因次半径位置的杯型切面入水过程的敞水特性、流场特性、通气腔分布形式以及表面压力的分布特性，对不同切面位置对半浸桨杯型切面入水过程的影响机理进行了分析，得到如下结论.

(1) 对于841-B型号的半浸桨，其有效切面为0.7位置的杯型切面，相较于0.6位置的切面以及0.8位置的切面而言，有效切面入水过程形成的通气腔的厚度以及部分通气状态情况下的通气腔长度均较大，即半浸桨有效切面通气现象的剧烈程度更加明显.因此在半浸桨剖面设计的过程中，其有效剖面的形式设计优劣会对半浸桨的通气特性有较大影响，设计时应着重考量.

与全国其他地方戏剧如东北二人转、河南豫剧等相比，吕剧文化发展伴有许多困难和问题，发展相对滞后。面对巨大危机和挑战，吕剧在新时期新形势下何去何从？答案只有一个，出路只有创新。在新时期，宣传做大吕剧文化要走多元创新之路，主要应从吕剧宣传推介、人才培养、精品创编、表演水平和舞台设计等多个方面着手。

(2) 随着半浸桨杯型切面位置逐渐接近半浸桨的梢部，完全通气状态与部分通气状态之间的过渡状态会更趋向于发生在一个较大的进速系数位置，即会对完全通气状态和部分通气状态的转变形成一个滞后作用.除此之外，随着无因次半径剖面位置越靠近梢部，横向力系数和效率会越大，垂向力系数相差不大.

(3) 杯型剖面结构形式的不同，尤其是随边位置的差异，且垂直入水的速度大小之间的差异，使得叶剖面入水时的流场特性具有一定差异，进而导致不同切面位置的杯型切面在入水过程中表面压力分布的形式发生改变，进而影响到不同切面位置半浸桨杯型切面入水过程中的横向力系数以及垂向力系数的大小.