强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

板壳式换热器波纹倾角对换热和阻力性能影响

研究了波纹倾角（α=45°，60°，75°）对板壳式换热器单流道内单相流动与换热过程的影响，分析了波纹流道内的速度分布、湍动能分布、压力分布和温度分布；基于范宁摩擦阻力因子f和Nu与Re的关系，提出了板壳式换热器不同波纹倾角下换热特性和流动阻力特性的预测关联式；通过计算不同波纹倾角下的综合性能评价指标（performance evaluation criteria,PEC）和面积质量因子（j/f）综合评价了流动与换热性能。结果表明：波纹倾角是影响圆形板间流体流动形态的因素之一，随着波纹倾角增大，会出现“十字交叉流”向“曲折流”转变；所提出的关联式能够很好地预测板壳式换热器内阻力与换热性能，阻力特性偏差在±14%范围内，换热特性偏差在±7%范围内；α为45°时，j/f较大，流道内阻力相对小；α为60°时，PEC较大，流道内换热性能相对强。     

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

基于流固耦合方法的水轮机蜗壳结构振动及损伤机理研究

蜗壳流道内的内水压力是引起外围混凝土发生损伤的主要原因.基于流固耦合理论,并引入混凝土弹塑性损伤模型,建立了流体与蜗壳结构耦合振动分析的理论框架,提出了一套水轮机流道内水体流动诱发蜗壳外围混凝土振动损伤的数值计算方法.首先基于有限体积法建立水轮机蜗壳流道流动的数值模型,同时采用有限元方法建立蜗壳结构固体区域的三维有限元模型;进而将流体区域边界上动水压力作为外荷载实时传递给固体区域边界进行三维有限元非线性损伤瞬态分析,实现了大型水轮机蜗壳结构中流体流场到固体应力、位移场的单向耦合三维数值分析.通过计算分析得到了水体流动诱发蜗壳外围混凝土振动的响应规律以及混凝土损伤的发展规律.     

超音速混合化学激光的小信号增益

超音速扩散式化学激光器谐振腔内混合流场的流动情况对于光腔的光学性能有重要影响。本文在考虑流场情况下利用矩形型线与Gauss-Lorentz 形型线模型作了零功率吋的增益计算,同时进行了比较分析,本文方法的计算结果与用Navier-Stokes 方程的计算结果定性比较小信号增益值分布规律相同。     

一种改进的颗粒移动床μ(I)拟流体模型及应用

颗粒移动床在工业领域应用广泛, 发展实用可靠的颗粒移动床模型具有理论和应用价值. 本文基于颗粒流μ(I)模型, 补充局部颗粒体积分数与颗粒局部压力和局部颗粒流密度的关系式, 将移动床内密集颗粒处理成可压缩拟流体, 建立了颗粒流单相可压缩流μ(I)模型, 并建立了颗粒流?壁面摩擦条件, 在计算中对颗粒流拟黏度和拟压力项进行正则化处理. 采用上述模型与方法对3种典型散料在移动床缩口料仓内的流动进行模拟, 与实验对比, 得到了玻璃珠、刚玉球和粗沙的μ(I)模型参数, 分析了3种不同散料在料仓内的颗粒速度、体积分数等分布特性, 模拟结果较好地揭示了料仓内不同物料的整体流和漏斗流特性; 进而以玻璃珠为例, 对移动床颗粒单管绕流流动进行了模拟, 所得结果合理揭示了管流附近的流动特性. 计算结果表明, 对于本文的计算工况, 颗粒体积分数变化最大范围为0.510 ~ 0.461, 绝大部分区域流动惯性数小于0.1, 改进的单相μ(I)模型能合理预测出密集颗粒流移动床内的流动特性, 方法可行且较多相流算法能明显减小计算量.      

小管径管翅式换热器气压胀接成形机理研究

管翅式换热器是制冷行业中最常用的换热器形式,其换热管的胀接性能决定了换热器的传热性能.本文提出了管翅式换热器的三维流-固耦合模型,采用单向流固耦合瞬态数值模拟方法,对小管径管翅式换热器的流体和固体域的流动和变形特征开展了数值研究.计算结果表明:根据换热管和翅片的胀接成形要求和胀后管径要求,气压胀接压力的合理范围为P=12.5 MPa,与理论公式推导值一致.根据管翅应力随时间变化的规律可知,换热管接头处应力远大于其屈服极限66 MPa,翅片接头处应力刚好略大于其屈服极限132 MPa,满足胀接成形要求.胀后的换热管直径随着压力的增加其管径增大,换热管的径向位移在水平方向较小,垂直方向较大,其最大和最小位移差约为0.03 mm.探究了管翅间残余接触压力随胀接压力的变化,残余接触压力随胀接压力的变化可分为三个阶段.结果表明当胀接压力使得翅片内孔发生屈服后,继续增大胀接压力会导致胀接不完全.最后研究了保压时间的影响,结果表明保压时间的增加对胀接效果并没有明显影响.相关结果可为工程实际中小管径管翅式换热器气压胀接工艺提供理论指导.     

爆破抛体加速运动的渗流模型计算

本文对有限球体爆破时渗流损失对抛体动能的影响作了计算。固体介质作为松散体处理,服从库仑内摩擦定律,破坏后的流动状态服从剪胀条件;爆炸生成的高压气体存松散介质内按渗流二次式规律流动,考虑线性摩阻与二次摩阻。 本文对一定块度,不同规模,不同比药量,不同初始松散率情况下的渗流影响作了计算。计算结果有助于进一步阐明工程爆破中抛体质心速度随爆破规模的变化规律,指出除了重力影响之外,还存在渗流影响。     

对二维分离流涡黏性系数非线性分布的新认识

以弱非线性涡黏性模型为出发点,对Delery分离流动实验结果进行分析并获得了非平衡态分离区涡黏性系数与形状因子J之间的非线性关系. 该非线性关系显示在分离起始阶段,涡黏性系数较平衡态先减小,后增大;再附阶段,涡黏性系数较平衡态数值逐渐增大,并在再附点位置接近最大,而后又逐渐减小,恢复到平衡态水平. 总结涡黏性系数的这种非线性发展数学关系式,并将它应用于BL模型,在不添加微分方程的情况下发展出一种适用于分离流动的改进代数湍流模型. 对低速平板流动,跨声速,超声速以及高超声速分离流动的计算结果表明,该改进湍流模型可以较准确地模拟各类复杂分离流动,计算精度明显优于传统代数模型以及一些两方程模型,而计算工作量仍与BL模型相当. 这表明所提出的涡黏性系数非线性发展规律是正确的,且应用在二维分离流动中具有一定的普适性.      

欧拉坐标系下具有锐利相界面的可压缩多介质流动数值方法研究

可压缩多介质流动问题的数值模拟在国防和工业领域内均具有重要的研究价值,诸如武器设计、爆炸安全防护等,通常具有大变形、高度非线性等特点,是一项极具挑战性的研究课题. 本文提出了一种基于 Euler 坐标系的非结构网格、具有锐利相界面的二维和三维守恒型多介质流动数值方法,可用于模拟可压缩流体和弹塑性固体在极端物理条件下的大变形动力学行为. 利用分片线性的水平集函数重构出单纯形网格内分段线性的相界面,并在混合网格内构建出具有多种介质的相界面几何结构,理论上可以处理全局任意种介质、局部 3 种介质的多介质流动问题. 利用传统的有限体积格式来计算单元边界上同种介质间的数值通量,并通过在相界面法向上求解局部一维多介质 Riemann 问题的精确解来计算不同介质间的数值通量,保证了相界面上的通量守恒. 提出了一种非结构网格上的单元聚合算法,消除了由于网格被相界面分割成较小碎片、违反 CFL 条件,进而可能带来数值不稳定的问题. 针对一维多介质 Riemann 问题、激波与气泡相互作用问题、浅埋爆炸问题、空中强爆炸冲击波和典型坑道内冲击波传播问题开展了数值模拟研究,将计算结果与相关的理论、实验结果进行比对,验证了数值方法的正确性和可靠性.      

A new numerical scheme for convective dominated heat transfer in a porous medium with strong temperature gradients

A new numerical scheme, theimplicit correction scheme, has been developed for heat transfer in a porous medium with strong temperature gradients. The scheme includes diffusion, convection and transverse heat transfer processes. By using correction coefficients which are based on transverse heat transfer, the effects of convection enthalpy flow and diffusion are modified. Under suitable limiting conditions, the implicit correction scheme can be reduced to the central-difference, upwind, or power-law scheme. The correction scheme is shown to be especially useful in calculations of the thermal effectiveness of the regenerator in Stirling cycle refrigeration.     

Numerical study of steady/unsteady flow and heat transfer in porous media using a characteristics-based matrix-free implicit FV method on unstructured grids

In this study, a high-resolution characteristic-based finite-volume (FV) method on unstructured grids [Int. J. Numer. Method Eng. 50 (2001) 11; Int. J. Heat Fluid Flow 21 (2000) 432] is extended by a matrix-free implicit dual-time stepping scheme for the numerical simulation of steady and unsteady flow and heat transfer with porous media. The method has been used to study the characteristics of a complex problem: flow and heat transfer in a channel with multiple discrete porous blocks, which was originally proposed by Huang and Vafai [J. Thermophys. Heat Transfer 8 (3) (1994) 563]. In addition, flow and heat transfer in a channel partially or fully filled with porous layers and containing solid protruding blocks with constant heat flux on its lower surface are also investigated in details. Hydrodynamic and heat transfer results are reported for both steady and transient flow cases. In particular, the effects of Darcy and Reynolds numbers on heat transfer augmentation and pressure loss are studied. An in-depth discussion of the formation and variation of recirculation is presented and the existence of optimum porous insert is demonstrated. At high Reynolds numbers the flow in the porous channel exhibits a cyclic characteristics although unlike the non-porous channel flow, the cyclic vortex development is only restricted to a small area behind the last solid block, while temperature changes more slowly and does not exhibit cyclic variations over a long period of time. It is shown that for all the cases studied altering some parametric values can have significant and interesting effects on both flow pattern as well as heat transfer characteristics.     

Buoyancy-driven flow in an enclosure containing time periodic internal sources

A computational study is carried out to investigate the effect of the sinusoidally driven heat source on the fluid flow and heat transfer within a two-dimensional square cavity. The cavity, which has solid walls of constant temperature, is filled with a fluid including uniformly distributed internal heat source. In addition, the effects of the different periods of the sinusoidally driving heat source on heat transfer are investigated and presented as figures.     

Steady and Unsteady Heat Transfer in a Channel Partially Filled with Porous Media Under Thermal Non-Equilibrium Condition

Steady and pulsatile flow and heat transfer in a channel lined with two porous layers subject to constant wall heat flux under local thermal non-equilibrium (LTNE) condition is numerically investigated. To do this, a physical boundary condition in the interface of porous media and clear region of the channel is derived. The objective of this work is, first, to assess the effects of local solid-to-fluid heat transfer (a criterion indicating on departure from local thermal equilibrium (LTE) condition), solid-to-fluid thermal conductivity ratio and porous layer thickness on convective heat transfer in steady condition inside a channel partially filled with porous media; second, to examine the impact of pulsatile flow on heat transfer in the same channel. The effects of LTNE condition and thermal conductivity ratio in pulsatile flow are also briefly discussed. It is observed that Nusselt number inside the channel increases when the problem is tending to LTE condition. Therefore, careless consideration of LTE may lead to overestimation of heat transfer. Solid-to-fluid thermal conductivity ratio is also shown to enhance heat transfer in constant porous media thickness. It is also revealed that an increase in the amplitude of pulsation may result in enhancement of Nusselt number, while Nusselt number has a minimum in a certain frequency for each value of amplitude.     

The Use of Porous Fins for Heat Transfer Augmentation in Parallel-Plate Channels

In this study, a steady, fully developed laminar forced convection heat augmentation via porous fins in isothermal parallel-plate duct is numerically investigated. High-thermal conductivity porous fins are attached to the inner walls of two parallel-plate channels to enhance the heat transfer characteristics of the flow under consideration. The Darcy–Brinkman–Forchheimer model is used to model the flow inside the porous fins. This study reports the effect of several operating parameters on the flow hydrodynamics and thermal characteristics. This study demonstrates, mainly, the effects of porous fin thickness, Darcy number, thermal conductivity ratio, Reynolds number, and microscopic inertial coefficient on the thermal performance of the present flow. It is found that the highest Nusselt number is achieved at fully filled porous duct which requires the highest pumping pressure. The results show that using porous fins requires less pumping pressure with comparable high heat augmentation weight against fully filled porous duct. It is found that higher Nusselt numbers are achieved by increasing the microscopic inertial coefficient (A), the Reynolds number (Re), and the thermal conductivity of the porous substrate k ₂. The results show that heat transfer can be enhanced (1) with the use of high thermal conductivity fins, (2) by decreasing the Darcy number, and (3) by increasing microscopic inertial coefficient.     

Fully-developed conjugate heat transfer in porous media with uniform heating

We propose a computational method for approximating the heat transfer coefficient of fully-developed flow in porous media. For a representative elementary volume of the porous medium we develop a transport model subject to periodic boundary conditions that describes incompressible fluid flow through a uniformly heated porous solid. The transport model uses a pair of pore-scale energy equations to describe conjugate heat transfer. With this approach, the effect of solid and fluid material properties, such as volumetric heat capacity and thermal conductivity, on the overall heat transfer coefficient can be investigated. To cope with geometrically complex domains we develop a numerical method for solving the transport equations on a Cartesian grid. The computational method provides a means for approximating the heat transfer coefficient of porous media where the heat generated in the solid varies “slowly” with respect to the space and time scales of the developing fluid. We validate the proposed method by computing the Nusselt number for fully developed laminar flow in tubes of rectangular cross section with uniform wall heat flux. Detailed results on the variation of the Nusselt number with system parameters are presented for two structured models of porous media: an inline and a staggered arrangement of square rods. For these configurations a comparison is made with literature on fully-developed flows with isothermal walls.     

Heat Transfer and Particle Behaviours in Dispersed Two-Phase Flow with Different Heat Conductivities for Liquid and Solid

Liquid-solid two-phase flow with heat transfer is directly simulated, to investigate the effects of the ratios of heat conductivities (solid to liquid) and bulk solid volume fraction from dense to dilute situations. The interaction between fluid and particles is solved by our original immersed solid approach on a rectangular grid system. A discrete element method with a soft-sphere collision model is applied for particle-particle and particle-wall interactions. Governing equation of temperature is time-updated with an implicit treatment for the diffusion term, which enables robust simulation with particles of very high/low ratios of heat conductivities (from 1/1000 to 1000) to the fluid. The local heat flux at the fluid-solid interface is modelled by a new flux decomposition technique, and incorporated into the implicit scheme of the temperature. The method is applied to a 2-D particulate flow in a natural convection in a square domain at a relatively low Rayleigh number. In the dense condition, for the cases with high ratios of heat conductivity, the heat transfer is promoted by strong convection, while the particles of low ratios of heat conductivity tend to hinder the development of the temperature rise in the flow field, causing a weak convection and low Nusselt number. Under a condition of relatively low solid volume fraction, fixed particles only depress the heat convection as the number of particles and heat conductivity ratio increase. For the cases with freely-moving particles, on the other hand, heat conductivity of particles has a stronger influence on the heat transfer of the system than the number of particles. The above simulation results highlight the effect of temperature distributions within the particles and liquid.     

An updated three-zone heat transfer model for slug flow boiling in microchannels

This work proposes a novel physics-based model for the fluid mechanics and heat transfer associated with slug flow boiling in horizontal circular microchannels to update the widely used three-zone model of Thome et al. (2004). The heat transfer model has a convective boiling nature and predicts the time-dependent variation of the local heat transfer coefficient during the cyclic passage of a liquid slug, an evaporating elongated bubble and a vapor plug. The capillary flow theory, extended to incorporate evaporation effects, is applied to estimate the bubble velocity along the channel. A liquid film thickness prediction method also considering bubble proximity effects, which may limit the radial extension of the film, is included. The minimum liquid film thickness at dryout is set to the channel wall roughness. Theoretical heat transfer models accounting for the thermal inertia of the liquid film and for the recirculating flow within the liquid slug are utilized. The heat transfer model is compared to experimental data taken from three independent studies. The 833 slug flow boiling data points cover the fluids R134a, R245fa and R236fa, and channel diameters below 1 mm. The proposed evaporation model predicts more than 80% of the database to within ±30%. It demonstrates a stronger contribution to heat transfer by the liquid slugs and correspondingly less by the thin film evaporation process compared to the original three-zone model. This model represents a new step towards a complete physics-based modelling of the bubble dynamics and heat transfer within microchannels under evaporating flow conditions.     

Analytical considerations of flow boiling heat transfer in metal-foam filled tubes

Flow boiling in metal-foam filled tube was analytically investigated based on a modified microstructure model, an original boiling heat transfer model and fin analysis for metal foams. Microstructure model of metal foams was established, by which fiber diameter and surface area density were precisely predicted. The heat transfer model for flow boiling in metal foams was based on annular pattern, in which two phase fluid was composed by vapor region in the center of the tube and liquid region near the wall. However, it was assumed that nucleate boiling performed only in the liquid region. Fin analysis and heat transfer network for metal foams were integrated to obtain the convective heat transfer coefficient at interface. The analytical solution was verified by its good agreement with experimental data. The parametric study on heat transfer coefficient and boiling mechanism was also carried out.