将压皱成型的波纹板按一定的角度相互交叉堆叠即构成了相互连通的交叉波纹通道。流体经该通道流动的时候，由于受到壁面高低不平的扰动，会使流动方向不断变化，并产生大小不同的漩涡。这一现象，对热质交换具有正面的意义，故其被广泛应用于板式换热器、分馏塔、冷却塔以及直接蒸发冷却填料等热质交换设备。对此交叉波纹通道流动特性的研究，使其传热传质效率更高，阻力更小，具有重要的现实意义和理论价值。

  

图1 直接蒸发冷却填料实物Fig.1 Physical map of direct evaporative cooling packing

  

图2 正弦交叉波纹通道Fig.2 Sinusoidal corrugated channel

一个完整的直接蒸发冷却填料如图1所示，其正弦交叉波纹通道结构如图2所示。这种通道的性能与其材料、波纹高度、波距、波纹倾斜角度、填料厚度、截面尺寸以及运行参数等均有关系。在实验研究方面，Koca等[1,2]对3个不同几何结构，不同厚度的填料在不同风速下的压降进行了实验研究。Franco等[3]对4个不同几何结构的填料在不同水流量和不同风速条件下的压降进行了计算和实验研究。Focke等[4,5]利用水中的邻甲酚酞采用有限扩散电流技术（DNCT）对交叉波纹通道夹角对流体流动形态变化趋势的影响进行了可视化研究。Gaiser等[6]测量了该流道局部压力的变化情况。Stasiek等[7]对由这种通道所构成的板式热交换器其进行了较为详细的实验研究，并给出了雷诺数、波纹波距-波高比以及交叉角度对流动阻力的影响。在数值计算方面，Beshkani等[8]将蒸发冷却用交叉波纹通道简化为一边平直、另一边正弦波纹变化的二维通道。并考虑到结构和流动的对称性，仅对半正弦单元的流动进行研究。其假设水和空气之间流动相互垂直、气体做层流流动、交叉波纹通道表面水膜均匀并忽略其厚度，计算得到流场的速度分布和温度分布，并讨论了进口参数对蒸发效率和压力损失的影响。至于三维流场的计算，目前较多的是针对板式换热器所对应的通道。由于对整个热交换器进行三维计算需要强大的计算机硬件配置，所以对交叉波纹通道热交换器的三维计算大多仅限于对一个具有代表性的单元进行[9-11]。Ciofalo等[10]对一个标准单元格内雷诺数在103～104的流动进行了计算。认为此范围内的流动（尤其是在旋转式空气预热器工作的雷诺数范围1 500～3 000）处于过渡区和弱紊流区。对于如此复杂的通道结构，很难确定用哪一种模型更为合适，故作者分别采用了几个不同的模型，包括层流、标准k-ε模型、低雷诺数模型、直接模拟以及大涡模型。计算结果表明，最好的结果是低雷诺数模型和大涡模型，它能使预测的总沿程阻力系数、平均Nusselt数以及局部Nusselt数比较合理。Mehrabian等[12]对标准单元按照层流假设，进行了微观计算。而对于正弦交叉波纹通道内全流场的计算，Poulter等[13]进行了正交波纹通道宏观的CFD计算，其假设正交波纹通道的流动与两平板间流动相同，忽略了板件复杂结构的影响，给出了大流动分布。Zhao等[14]对这种通道所构成的板式热交换器全流场在不同角度下的性能进行了数值模拟，但该模拟仅给出了二维平面上的压力场和速度场的结果，对流动的细节未作较为详尽的描述。

综上所述，现有的研究大多局限在对一个标准单元格、标准通道或简化结构的研究，其结果与实际情况尚有较大差距。为得到更为贴近实际填料内部气流流动的分布特点，本工作将对由两片相互交叉的正弦波纹板所形成的三维通道，按照实际流动边界条件进行全流场紊流流动的计算，希望得到波纹通道内部的流动规律。

1 物理模型及假设

1.1 物理模型

本工作将以两个填料片之间区域组成的流道作为研究对象，如图3所示，各尺寸如表1所示。两层波纹纸的波纹通道间交叉夹角为90°，且每个波纹纸通道与进口主气流方向夹角为45°。在两波纹片的接触处，考虑到实际蒸发冷却填料采用黏胶剂将两波纹片粘接在一起的特点，如图4所示，以接触点为圆心、将半径为0.75 mm的圆柱面和两个波纹片之间的包罗区域内视作填充黏胶剂，空气在该区域之外的空间内流动，以更加符合实际流动状态。

 

表1 波纹填料尺寸及坐标范围Table 1 Corrugated packing size and coordinate range

  

Packing model Wave height (Hi) /mm Wave distance (P) /mm Ripple angle (θ) / ° X/mm Y/mm Z/mm T7090 7 29.0 90 -37.5-37.5 0-14 -72.5-72.5

  

图3 计算区域结构示意Fig.3 The structure diagram of calculation area

  

图4 接触点附近区域局部放大Fig.4 Partial enlargement of the area near the contact point

1.2 模型假设

为了使模拟趋于简单化，本工作做出如下假设：（1）整个装置处于一种稳定的工况下运行；（2）空气为不可压缩的牛顿流体，忽略质量力；（3）进口速度均匀，且垂直于填料进口断面。

2 控制方程及其求解方法

2.1 控制方程及紊流模型

该流动的控制方程包括连续性方程和运动方程：

虽然各波纹槽交叉断面处的流量有上述不同的变化，但其变化的幅度在3×10-5 kg/s范围内，约为断面流量的12.6%。沿主波纹槽流动仍是主要部分，这与Gaiser等[6]的结论一致。即在两波纹槽夹角为0＜θ＜160°时，除边壁的附近流动方向改变外，大多数气流均沿每一波纹板主波纹槽流动。只有当气流夹角增大后，才会形成Z字形的平行流动。

 

式中，Cµ = 0.09，取l = 0.07 Dh，Dh为填料结构的水力直径。由于沿流动方向上的通道截面积不固定，水力直径参照文献[16]定义为：

 

图8为各波纹槽不同截面上的流量变化曲线，图9为根据图8的数据所整理的各波纹槽流量变化的示意图。从各波纹槽的流量变化来看，由于填料右边壁的阻挡作用，受阻挡作用的各波纹槽在与右边壁接触之前的流量逐渐减少，而且影响的范围较大。如图9所示的波纹槽B，其受影响的截面甚至有4个之多。而在与边壁接触之后，其流量则相应逐渐增加。这与图12各波纹槽不同截面上的压强变化曲线所表现的接近边壁附近的压强较大而中间部分的压强较小的结果相一致。同样，由于在出口断面上，受边壁附近压强较大的影响，接近出口端面的波纹槽C和B，其流量相应也是大大减小，并小于平均值。在出口端面中心位置，由于上下两波纹槽最大流量出现交汇现象，引起出口端面在该位置附近的压强有一定增加，流量相应也有一定的减小，但仍高于平均值。由此，最终造成各波纹槽出口端面的流量相差很大，最大的A波纹槽的流量约比最小的B波纹槽的流量大24.2%。由于计算区域的局限性，本算例所有波纹槽的流量均因左右边壁的影响而改变，没有一个波纹槽在出口端面的流量等于其进口流量。只是在远离左右边壁的中间位置和主流动方向远离出口端面的位置，由于上下波纹槽的对称性，其交叉波纹槽的流量相互影响很小，甚至互不影响。

 

式中，紊流动能（k）和紊流耗散率（ε）分别由式（5）和（6）求出：

保健目标：严格控制呼吸道疾病的发生，预防细菌或病毒性疾病的出现，控制了猪呼吸道疾病，也就减少了其它病原菌的感染机会，减少对其它疾病的药费支出；清除后备母猪体内病原菌及内毒素，抑制体内病毒数量及活性，用药物对猪场内的常见病进行净化；增强后备母猪的体质，提高机体免疫力，促进发情，获得最佳配种率。

 

式中：；ζ为离开壁面的距离；u为平行于壁面的速度。其他各符号的意义以及取值，参见文献[15]。

2.2 边界条件

进口边界条件：空气流速为ux=0，uy=0，uz=u0。

最后，创客教育强调“为了每一个学生的发展”，关注的是人人参与创新。教师需遵循“学生为主体，教师为主导”的基本要求，设计合理的教学模式，进而促进学生的全面发展。

紊流动能：

 

式中，I为紊流强度，在此取5%。

紊流耗散率：

2018年5月,全国生态环境保护大会上,习近平总书记对全面加强生态环境保护,坚决打好污染防治攻坚战,作出了系统部署和具体安排,“习近平生态文明思想”由此确立。

 

式中：是作用在微元体表面上的紊流惯性切应力，采用Boussinesq假设：

 

式中：V代表计算区域的体积；Aw代表计算区域的表面积。

出口边界条件按照局部单向化边界条件处理。

2.3 求解方法

本工作使用Solidwork 2013软件建立了交叉波纹通道填料的三维计算模型。当模型建好后，将其导入到Workbench15.0中并使用ANSYS Meshing对其进行网格划分，采用非结构化四面体网格对物理模型进行网格划分。为有效求解近壁处较大的速度梯度，在近壁面附近布置相对较密的膨胀层，最终划分网格总数为16 829 606，图5显示了整个物理模型的部分网格。

  

图5 网格划分Fig.5 Schematic diagram of meshing

Bering等人(2005)的研究发现，儿童更倾向认为个体死亡后其认知、愿望、情绪功能仍存在，而认为死亡后生物、精神生物、知觉功能会停止。我们的结果重复了儿童来生信念理解的这一特点。

3 模拟结果分析

3.1 结果合理性验证

对所研究对象在不同风速下的沿流动方向上平均压降的数值模拟结果如图6所示，并与文献[17]的产品资料和文献[2]的实验数据进行比较。由图可见，3组数据符合较好，计算结果有效。

  

图6 流动方向压降随风速的变化Fig.6 Flow direction pressure drop with the air velocity

  

图7 下波纹片波纹槽通道图示Fig.7 Lower corrugated groove channel diagram

3.2 各波纹槽流量的变化规律

相互成一定角度排列的交叉通道，实际上就是由两块相同波形、波高和波距的波纹型填料板，以与进口气流方向正负相同波纹倾斜角的方式拼装而成。在实际的计算中，为了便于结果分析和计算合理性检验，将两块波纹板布置成交叉对称形。由于其结构的对称性，并且在计算中不考虑重力的影响，故其流动在各波纹板的波纹槽内理论上亦具有对称性。实际的计算结果显示，上下波纹槽在对应位置处的流量差均小于3.22%，完全可以看作为对称流动。将此两块波纹板水平摆放来看，去掉上面波纹板，可以只研究如图7所示的下波纹板波纹槽内的流动特点。图7中细实线代表下波纹片的波峰线，相邻两细实线之间为下波纹片波纹槽，其分别用大写字母表示。虚线代表上波纹片的波谷线，相邻两虚线之间为上波纹片波纹槽。其中C和D是两个几何上的完整直波纹槽，数据整理过程中，其流量为该槽内的实际流量。而对A，B和E等3个几何上不完整波纹槽，在流体与右边壁冲击之前，其流量为该槽内的实际流量。在与右边壁冲击之后，其流量则应为方向改变90度后上波纹片相应槽内的流量。但考虑到流动的对称性，上波纹片槽内的流量应等于下波纹片对称槽内的流量。这样，在分析各波纹槽内的流量时，按照图7所示的各波纹槽进行分析。为了确保各波纹槽流量计算的准确性，仅计算与上波纹片波谷交叉的截面上的流量，各交叉截面用相应波纹槽编号的下标数字进行区分。其中2～5分别表示波纹槽内部各交叉截面，下标数字1表示进口断面，下标数字6表示出口断面。

免疫接种具有不可替代的作用，然而因重视程度不足，导致很多羊养殖户存在侥幸心理，认为加强饲养管理即可。然而，一部分致病原在不同因素的作用下会导致羊患病[3]。

  

图8 各波纹槽不同截面上的流量变化曲线Fig.8 Flow curve of different sections of each corrugation groove

  

图9 各波纹槽流量变化示意Fig.9 Schematic diagram of flow changes in each corrugation groove

式中，μt为紊动黏度。根据Ciofalo等[10]的结论，采用Jones-Launder的低Re数k-ε模型进行模拟。

在春季，在养殖业的过程中，鸡容易产生多种疾病。在育种过程中，只有仔细观察才能做出更好的判断和分析育种者用它做出更好的判断疾病的处理。否则，它会影响鸡的正常繁殖，并在同一时间，它会影响养鸡业的经济效益。因此，只有通过正确的判断和分析，才能采取更加针对性的预防措施来减少养殖户的损失。

虽然沿波纹槽流动方向上的流量受边界的影响较大，最终使各波纹槽的出流流量也不尽相同。但各波纹槽相互间的流量仍存在着某种规律。以上波纹片的波纹槽为基准，如图10所示，分析其波纹槽I，II，III和IV所对应的交叉单元格中下波纹槽内的流量变化规律。下波纹槽各单元格的流量采用其前后截面流量的平均值计算，其结果如图11所示。图中，单元格位置用其前后两个截面位置代号表示。从图中可以看出，对于上波纹槽I，II和III，其中间的3个单元格上所对应的下波纹槽的流量沿上波纹槽流动方向近似呈线性增加趋势，而且这种线性增加的趋势对于3个通道而言，基本相同。对于上波纹槽道IV，其前两个单元格上所对应的下波纹槽的流量保持着和其他3个波纹槽相同的规律，只是在第3个单元格上其流量才受左边壁的影响而大大减小。由此说明，边壁对各波纹槽流动的影响是整体性的和渐进性的。

  

图10 波纹槽位置示意Fig.10 Corrugated groove position diagram

  

图11 波纹槽流量对比Fig.11 Corrugated groove flow comparison

3.3 各波纹槽压强的变化规律

如前所述，各波纹槽由于受左右边壁的影响，宏观上讲，在沿主流动方向的横截面上，压强分布不均匀，图12为各波纹槽不同截面上的压强变化曲线。由图可看出，在接近左右边壁附近压强较大，

参考文献[14]的方法，在2 mL样液中加入2 mL 0.1 mmol/L DPPH·溶液，混匀室温下避光静置20 min后，在517 nm处测定吸光值，DPPH·清除率按下列公式计算。

而在中心部分压强较小。从沿各波纹槽流动方向上的压强变化来看，在进口的位置，由于流动速度方向急剧发生变化，各波纹槽在此位置处的压力损失均为最大，并近似相同。其后，各波纹槽的压力降则不受进口流动的影响，而主要受与边壁相对位置变化的影响，其压力降各不相同，交错变化。总的来讲，靠近左右边壁附近的压力降小于中间位置的压力降。出口端面各波纹槽平均压强最大值与最小值相差约5.12 Pa，占平均压力降的29.8%。

本文就旅游产品供应角度，提出关于西沙群岛海岛旅游开发的新思路，将海岛旅游和定制服务思路结合研究。定制服务开发对彰显海岛自然环境的休闲性有极大的作用，有利于增强西沙群岛旅游的服务质量，让游客最大限度地参与到产品的设计中，游客体验水平提升，将提高游客满意率，提升海南省旅游实力。同时，开发定制服务旅游方式可减少旅游企业设计产品的成本支出，减少资源浪费，提升运作效率，从而避免市场恶性竞争情况的出现，促进旅游企业不断创新，不断发展。

  

图12 各波纹槽不同截面上的压强变化曲线Fig.12 Pressure curve of different sections of each corrugation

3.4 旋涡的流动规律

对于这种交叉通道，其最大的特点就是各波纹槽除了有主流方向的流动外，还有因槽道交叉影响所产生的如图13所示的二次漩涡流动。在本算例中，各完整单元格内下波纹槽旋涡的流量变化曲线如图14所示。单元格波纹槽旋涡流量定义为各单元格内穿过漩涡中心轴与波纹槽间垂直平面的流量。图中单元格位置同样用其前后两个截面位置代号表示。由图可以看出，各旋涡流量和主流流量之比均在13%～50%，占比相当大。该旋涡的存在，加速了边界层内流动的扰动。Gaiser 等[6,10]均认为，二次流是热质交换的决定性因素，涡强度对传热速率有正面影响，对漩涡的研究具有重要的意义。

  

图13 典型单元格上的二次流漩涡Fig.13 Secondary flow vortex on a typical cell

  

图14 漩涡流量变化曲线Fig.14 Vortex flow curve

从图14可以看出，各完整单元格内下波纹槽旋涡的流量除了靠近边壁附近的单元格外，其他各波纹槽的基本相同，并且沿波纹槽流动方向上呈相同下降衰减的趋势。而此下降衰减的速率沿波纹槽主流动方向上则逐渐减小。在靠近边壁附近，波纹槽E的靠近左边壁的2-3截面间旋涡流量较其他单元波纹槽的要小，而波纹槽B的靠近右边壁的4-5截面间旋涡流量则较其他单元波纹槽的要大。这充分说明，漩涡产生的主要原因是由于各波纹槽内气流方向的改变所致。在进口位置，如图15所示，气流由原来的垂直于进口端面方向突然改变 45°角，变成沿波纹槽道的方向流动，迫使气流在改变方向的同时产生很强的二次流旋涡。该旋涡在随波纹槽主流流动过程中，由于受壁面摩擦的影响而衰减。同时，上波纹片波纹槽内主气流的流动也在一定程度上会使下波纹槽内的流动产生二次流旋涡。所以，上述旋涡在沿波纹槽主流动方向上的衰减会逐渐减小。当波纹槽足够长后，可以预见，其衰减和产生会达到平衡，旋涡流量维持不变。这也正是本工作所计算的旋涡流量之所以比Ciofalo等[10]和Gaiser等[6]所计算的标准单元格充分发展流时旋涡流量大的根本原因。由于在进口位置，各波纹槽流入的流量相同，波纹槽结构相同，这就使得旋涡的流量基本相同。尽管在流动过程中，各波纹槽的主流流量发生变化，但由于该旋涡主要在入口处产生，而各波纹槽的旋涡衰减条件基本相近。故沿流动方向上各不同波纹槽的旋涡流量基本相同。在波纹槽E的靠近左边壁的2-3截面间，由于在2截面前的单元格内左边壁的导流作用，致使该波纹槽内的气流速度方向改变较为缓慢。如图16所示，其二次流漩涡相比其他单元格略微较小。但由于其在2-3截面间由上波纹槽带入的流量较大，其漩涡的流量衰减响应较慢。而波纹槽B在4-5截面漩涡的增大，主要是由于受右边壁的阻滞作用，主流方向的气流被迫旋转90°，流入上波纹槽内。从下波纹槽的角度来看，其就是大量的二次旋涡。

教学包含着教和学两个方面。要想让学生学到并牢固掌握老师所教的知识，必须使他们在学习过程中积极思维，从而发现问题并提出问题，通过老师解答之后增强理解，加深印象。

  

图15 C2截面速度矢量Fig.15 C2 section velocity vector

  

图16 E2截面速度矢量Fig.16 E2 section velocity vector

4 结 论

a）各波纹槽的流量受填料边壁的阻挡作用影响较大。受阻挡作用的各波纹槽在与边壁接触之前的流量逐渐减少。而在与边壁接触之后，其流量则相应逐渐增加。与其相对应的是，接近边壁附近的压强较大而中间部分的压强较小。在出口端面中心位置，由于上下两波纹槽最大流量出现交汇现象，引起出口端面在该位置附近的压强有一定增加，流量相应也有一定的减小。两个因素最终造成各波纹槽出口端面的流量相差很大。本算例中，最大波纹槽的流量约比最小波纹槽的流量大1/3之多。

b）边壁对各波纹槽流动的影响是整体性的和渐进性的。在计算域的中间位置，各波纹槽所对应的交叉波纹槽内的流量之间，近似呈线性变化趋势。而且这种线性变化的趋势对于各波纹槽通道而言，基本相同。

将划分好的网格导入到Ansys Fluent15.0中，并对区域材料以及边界类型进行设定。本工作选用紊流模型中的低Re数k-ε模型，计算方法采用速度和压力耦合的SIMPLE算法，压力项采用线性插值的方法。为了加快迭代收敛的速度，改变动量方程的松弛因子为 0.8，其他保持默认值。再通过求解器对各方程进行迭代求解，各方程参数的残差都默认为10-3，当残差低于这一数值时，即认为收敛，最终得到该流动的数值解。

c）在进口的位置，由于流动速度方向急剧发生变化，各波纹槽在此位置处的压力损失均为最大，并近似相同。其后，各波纹槽的压力降则不受进口流动的影响，而主要受与边壁相对位置变化的影响，其压力降各不相同，交错变化。

d）各波纹槽内存在有较大的二次流漩涡，各旋涡流量和主流流量之比均在 13%～50%。漩涡的产生主要是由于各波纹槽内气流方向的改变所致，各完整单元格内波纹槽旋涡的流量除了靠近边壁附近的单元格外，其他各波纹槽的基本相同，并且沿波纹槽流动方向上呈相同下降衰减的趋势。而此下降衰减的速率沿波纹槽主流动方向上则逐渐减小。

符号说明

Aw——计算区域表面积，m2 p——压强，Pa

C1ε——耗散率方程产生项引入的经验系数 u——速度，m/s

Cµ——k-ε模型引入的紊动黏度经验系数 u0——进口初速度，m/s

C2ε——耗散率方程消失项引入的经验系数 V——计算区域体积，m³

收集该院90例应用微量泵注入胰岛素治疗的ICU糖尿病危重患者，数字表法分组，全面化针对性护理干预组男患者24例，女患者21例。年龄范围35～78岁，平均年龄（54.24±2.78）岁。急性脑血管疾病22例，重症胰腺炎12例，重病肺炎8例，农药中毒3例。血糖异常时间持续1～8 d，平均(3.26±0.35)d。对照组男患者24例，女患者21例。年龄范围35～78岁，平均年龄 （54.01±2.64）岁。急性脑血管疾病21例，重症胰腺炎12例，重病肺炎8例，农药中毒4例。血糖异常时间持续1～8 d，平均(3.22±0.42)d。两组一般资料可比。

Dh——填料结构的水力直径，m μ——流体黏度，Pa·s

f1——低雷诺数模型引入的耗散率方程产生项修正系数 μt——紊动粘度，Pa·s

f2——低雷诺数模型引入的耗散率方程消失项修正系数 ρ——密度，kg/m³

fμ——低雷诺数模型引入的紊动黏度修正系数 ε——紊动耗散率，m2/s3

一起案件，开庭两次，这种情况不多见，但也绝不少见。但是两次庭审之间，时间跨度达一年之久的，着实不多。这约一年的时间里，李凌所做的就是不厌其烦地翻看卷宗和整理证据，李凌说：“那段时间，我每天就是躺在床上了，脑子里都在回想案情，一想到有什么错漏，我就马上爬起来看卷宗。”

Gk——紊动能产生项，Pa/s ζ——离开壁面距离，m

Hi——波纹板波高，m θ——填料交叉角度，°

I——紊流强度 σk——紊动能普朗特数

k——紊动能, m2/s2 σε——紊动耗散率普朗特数

在实际的施工现场，由于人工施工振捣不均匀，会在施工过程中出现蜂窝等一系列问题。混凝土用料速度太快，会增加分离浇灌难度；振捣或者用料配比不科学的话，会导致产生振捣不断向下倾斜、振动出现麻面现象。如果模板卫生清洁工作做不到位或者表面粗糙的话，会附着干燥坚硬的水泥砂砾等杂质，拆开模板的过程中混凝土外表会出现麻面现象。同时，如果模板组装没有做到严密的时候，浇灌建筑混凝土会出现孔洞概率高的现象。

l——填料特征长度，m

参考文献：

[1] Koca R W, Hughes W C, Christianson L L. Evaporative cooling pads: test procedure and evaluation[J]. Applied Engineering in Agriculture,1991, 7(4): 485-490.

[2] Malli A, Seyf H R, Layeghi M, et al. Investigating the performance of cellulosic evaporative cooling pads[J]. Energy Conversion &Management, 2011, 52(7): 2598-2603.

[3] Franco A, Valera D L, Peña A, et al. Aerodynamic analysis and CFD simulation of several cellulose evaporative cooling pads used in Mediterranean greenhouses[J]. Computers & Electronics in Agriculture, 2011, 76(2): 218-230.

[4] Focke W W, Zachariades J, Olivier I. The effect of the corrugation inclination angle on the thermohydraulic performance of plate heat exchangers[J]. International Journal of Heat & Mass Transfer, 1985, 28(8): 1469-1479.

[5] Focke W W, Knibbe P G. Flow visualization in parallel-plate ducts with corrugated walls[J]. Journal of Fluid Mechanics, 1986, 165: 73-77.

[6] Gaiser G, Kottke H V. Flow phenomena and local heat and mass transfer in corrugated passages[J]. Chemical Engineering & Technology,1989, 12(1): 400-405.

[7] Stasiek J, Collins M W, Ciofalo M, et al. Investigation of flow and heat transfer in corrugated passages-I: experimental results[J].International Journal of Heat & Mass Transfer, 1996, 39(1): 149-164.

[8] Beshkani A, Hosseini R. Numerical modeling of rigid media evaporative cooler[J]. Applied Thermal Engineering, 2005, 26(5): 636-643.

[9] Doo J H, Ha M Y, Min J K, et al. Theoretical prediction of longitudinal heat conduction effect in cross-corrugated heat exchanger[J].International Journal of Heat & Mass Transfer, 2012, 55(15/16): 4129-4138.

[10] Ciofalo M, Stasiek J, Collins M W. Investigation of flow and heat transfer in corrugated passages II: numerical simulations[J]. International Journal of Heat & Mass Transfer, 1996, 39(1): 165-192.

[11] Doo J H, Ha M Y, Min J K, et al. An investigation of cross-corrugated heat exchanger primary surfaces for advanced intercooled-cycle aero engines (Part-I: novel geometry of primary surface)[J]. International Journal of Heat & Mass Transfer, 2012, 55(19/20): 5256-5267.

[12] Mehrabian M A, Poulter R, Quarini G L. Hydrodynamic and thermal characteristics of corrugated channels: experimental approach[J].Experimental Heat Transfer, 2000, 13(3): 223-234.

[13] Poulter R, Quarini G L, Mehrabian M. Flow and local heat transfer characteristics of an A.P.V. junior paraflow plate heat exchanger[J].Proceedings of the HTFS Conference, 1994, RS(963): 63-66.

[14] Zhao Y H, Wu Y F, Cheng H J, et al. Numerical Simulation of Corrugated Depth on the Performance of plate heat exchanger[J].Advanced Materials Research, 2014, 860/863: 696-699.

[15] 陶文铨. 数值传热学[M]. 西安: 西安交通大学出版社, 2003: 347-350.

[16] Kim M, Man Y H, Min J K, et al. Numerical study on the cross-corrugated primary surface heat exchanger having asymmetric cross-sectional profiles for advanced intercooled-cycle aero engines[J]. International Journal of Heat & Mass Transfer, 2013, 66(11): 139-153.

[17] He S, Gurgenci H, Guan Z, et al. Pre-cooling with Munters media to improve the performance of natural draft dry cooling towers[J]. Applied Thermal Engineering, 2013, 53(1): 67-77.