固定床反应器内流动的研究

固定床反应器是化学工业中常用的反应设备。详细研究反应器内的流动,对提高反应器的性能很有意义。本文分别对轴向及径向固定床反应器内的流动作了分析。对轴向床反应器,它的主要问题是流经催化剂床层的流速分布不均匀。作者分析了其原因,提出了用组合式均布装置来解决流速均布的问题,並用冷模实验进行验证。实验结果表明,本装置可使催化剂床层上方截面上的流速分布,相对标准偏差<9％,均布装置的阻力系数ζ≈1.35。对径向床反应器,作者导得了反应器内流体流动的基本方程,並进行实验研究。研究结果表明,反应器结构型式、催化剂床层阻力系数、中心通道与环形通道的面积比以及反应器高径比是影响流经床层流速分布的主要因素。     

高频等离子体尾焰接地熄弧的原因及解决办法

高频等离子体在化工、冶金中有着广泛的应用.国内先后用于制造钛白、光导纤维及氮化硼等新材料.但在试验或生产中常遇到下述问题:与高频等离子体发生器相连的化学反应器一旦接地,系统就无法正常运转.经过反复试验发现,由于高频等离子体尾焰很长,时而会接触连在后面的反应器器壁,而一般反应器又是导电的,这时反应器的(电)接地就会引起熄 ...     

轴向流固定反应器流场计算研究——数学模型与计算方法

本文提出了轴向流固定反应器的流动模型和基本方程。运用有限差会计了这种反应器内部的流场。     

轴向流固定床反应器流场计算研究——数学模型与计算方法

本文提出了轴向流固定床反应器的流动模型和基本方程。运用有限差分法计算了这种反应器内部的流场。     

脉冲流在Y型和倒Y型微反应器中强化传质的格子Boltzmann模拟

连续流微反应器的迅速发展为化学合成技术提供了一条可精准控制的路径。微反应器中流体的流动、混合和传质是反应的物理基础,因此强化传递就能够合理改善混和效果,增加相接触面积,减小微通道尺寸,以缩短分子扩散距离。由于微反应器中对流传递数量级非常低,不容易控制。通过改进微混合器构型和引入脉冲流动等主被动强化措施,可以有效改善对流和传质,影响化合反应。如何量化分析这些影响,是微反应器强化传质和优化控制反应过程的基础。本文基于有效的虚拟串行竞争反应格子Boltzmann模型,通过对Y型和倒Y型的微反应器的流场结构、混合传质和化合反应进行数值研究,定量分析了脉冲流动在不同流场构型下传质和化学反应的影响,所得结论可为连续流微反应器设计以及微反应精准控制提供有意义的参考。     

爆震室壁面条件对脉冲爆震发动机爆震性能的影响

在爆震室内快速形成稳定传播的爆轰波是脉冲爆震发动机的关键．本文利用有限速率化学反应模型,考虑粘性、热对流,基于N-S方程对氢气与空气/氧气为反应混合物的爆震发动机爆震室内流场进行计算．从流场压力、速度、涡量、湍流动能等方面研究爆震室壁面条件对燃烧爆轰性能的影响,分析流场爆轰波压力与流场湍动能的关系,讨论可燃气体燃烧转爆轰的机理．结果表明：爆震室内燃烧爆轰机理受到化学反应能量释放、壁面摩擦效应、壁面与外界热交换的影响．在文中讨论的范围内,相比于半圆形和三角形的爆震室装置,矩形的爆震室增强装置能在更短的时间内得到较高的爆轰波压力和湍动能峰值．壁面粗糙层高度（粗糙度）影响爆震室的燃烧爆轰性质．当壁面粗糙度为0.15mm时,粗糙度对爆轰的激励作用大于抑制作用,能较快形成稳定的爆轰波,且推力为35.5N;随着壁面对流换热系数的增大,爆震室壁面的散热加剧．当壁面对流换热系数大于临界值2.6W/(m2·K)时,爆震室内不能形成稳定的爆震波．     

超音速/高超音速三维边界层的层流控制基金项目

根据可压缩黏性稳定性理论研究了壁面冷却和抽吸对超音速、高超音速三维边界层的层流控制作用.数值结果证明壁面冷却对第一模式起稳定作用,对第二模式有不稳定作用;壁面抽吸对第一、二模式都起稳定作用;直到Me＝7,导致绝热壁边界层转捩的始终是第一模式,Me≥6的冷却壁边界层则是第二模式对转捩起主导作用.壁面冷却能够推迟边界层转捩,但是和二维边界层相比壁面冷却对高速三维边界层的层流控制作用是很有限的.     

高速双锥绕流中热化学与输运模型影响研究

激波与边界层之间相互作用是高超声速飞行中的常见现象,对飞行器气动性能与飞行安全至关重要.对于高焓来流,流场中通常存在复杂的物理化学现象,此时准确模拟流场中激波边界层相互作用的难度大,相关物理化学建模仍有待进一步考察和研究.本文针对最近文献中纯净空气高超声速双锥绕流实验开展数值研究,分别研究了不同热化学模型与输运模型对壁面压力与热流的影响.热力学模型包括完全气体、热力学平衡和非平衡模型,化学模型包括冻结和非平衡化学模型,输运模型包括经典的Wilke/Blottner/Eucken模型与更加复杂的Gupta/SCEBD模型,以及考虑壁面催化/非催化影响的模型.计算了6个不同算例,涵盖了低焓至高焓来流等不同工况.壁面压力与热流的数值计算结果与实验结果符合较好;对于低焓来流,计算结果主要受到分子内能分布的影响,输运模型对计算结果的影响不大;对于高焓来流,一方面计算结果受到化学反应与壁面催化的影响较大,另一方面不同输运模型对计算结果的影响也更加明显.      

弯管内爆轰波传播的流场显示和数值模拟

采用激光纹影系统拍摄了爆轰波在不同位置的流场照片. 用二阶附加半隐的龙格- 库塔法和五阶WENO格式 分别离散欧拉方程时间和空间导数项，用基元反应来描述爆轰化学反应过程，获得了压力、 温度、典型组元质量分数分布及数值胞格结构和爆轰波平均速度. 结果表明：受壁面稀疏波 和压缩波影响，爆轰波阵面发生畸变. 但由于弯管曲率半径较大，未出现爆轰波熄灭. 靠近 凹壁面的激波强度大于凸壁面侧，且凹壁面侧的反应区宽度较凸壁面侧要窄. 弯管出口处的 三波点数目较入口处减少，爆轰波衰减. 在出口直段，受扰动的爆轰波可恢复为自持爆轰波. 爆轰波流场、胞格结构、平均爆轰波速度的计算和实验结果定性一致.     

高速双锥绕流中热化学与输运模型影响研究

激波与边界层之间相互作用是高超声速飞行中的常见现象,对飞行器气动性能与飞行安全至关重要.对于高焓来流,流场中通常存在复杂的物理化学现象,此时准确模拟流场中激波边界层相互作用的难度大,相关物理化学建模仍有待进一步考察和研究.本文针对最近文献中纯净空气高超声速双锥绕流实验开展数值研究,分别研究了不同热化学模型与输运模型对壁面压力与热流的影响.热力学模型包括完全气体、热力学平衡和非平衡模型,化学模型包括冻结和非平衡化学模型,输运模型包括经典的Wilke/Blottner/Eucken模型与更加复杂的Gupta/SCEBD模型,以及考虑壁面催化/非催化影响的模型.计算了6个不同算例,涵盖了低焓至高焓来流等不同工况.壁面压力与热流的数值计算结果与实验结果符合较好;对于低焓来流,计算结果主要受到分子内能分布的影响,输运模型对计算结果的影响不大;对于高焓来流,一方面计算结果受到化学反应与壁面催化的影响较大,另一方面不同输运模型对计算结果的影响也更加明显.     

Numerical simulation of the flow field in a reactor for titanium dioxide production

A physicochemical and fluid dynamic model is formulated for the numerical simulation of the flow field in a reactor for titanium dioxide production, the turbulence motion is described by theKε equation, the governing equations are solved by the SIMPLER algorithm devised by Partankar and Spalding. The velocity, tmperature and concentration fields are obtained for three cases: A) with chemical reaction and thermal insulation on the walls; B) with chemical reaction and wall temperature is 450K; C) without chemical reaction and thermal insulation on the walls, and the physicochemical numerical simulation for the titanium dioxide production has been done. The results of the paper can be used as a theoretical guide for the engineers in the design of such reactors.     

Numerical investigation on cold flow dynamics of supercritical water fluidized bed reactor with inclined distributor: Design and scale up

Supercritical water fluidized bed reactor (SCWFBR) is a novel concept for the gasification of coal and biomass to produce hydrogen. In this work, to enhance the mixing in the axial direction, an inclined distributor is introduced to optimize the flow dynamics in SCWFBR with partitioned fluid supply. Through numerical simulations based on the two fluid model (TFM), the effects of the inclined distributor structure and operating parameters on the solid distribution and the residence time are evaluated with the optimal values determined. Numerical results show that, area ratio = 2:1, SCW velocity ratio = 3:1, flow ratio = 3.36:1 and inclination angle = 20° are the optimal design in this paper. A predictive correlation of the minimum fluidization velocity for the improved SCWFBR is also proposed based on the numerical data. The average error between the correlation and numerical simulation results is approximately 1.4% which strongly demonstrates its capability. Finally, based on the optimal design, the lab-scale reactor is further scaled up and the studies about two scale-up rules are carried out. Only the cold flow is simulated in this study without considering chemical reaction which would be involved in future work.     

CFD modeling using heterogeneous reaction kinetics for catalytic dehydrogenation syngas reactions in a fixed-bed reactor

A comprehensive 2D computational fluid dynamics （CFD） model was developed to simulate the flow behavior and catalytic dehydrogenation reaction of syngas in a heterogenous fixed-bed reactor （FBR）. The model combined the porous medium CFD model with a reaction kinetics model. To acquire an accu- rate reaction kinetics model, a comprehensive reaction mechanism was studied for the heterogeneous catalytic dehydrogenation reaction ofsyngas over a supported metal catalyst. Based on the reaction mech- anism and a statistical test, a reliable kinetics model was proposed. The CFD model combined with the above kinetics model was validated with one set of experimental data. The CFD model was also used to predict key reaction variable distributions such as the temperature and the component concentrations in the reactor.     

THE DYNAMIC BEHAVIORS OF VISCOELASTIC PIPE CONVEYING FLUID WITH THE KELVIN MODEL

Based on the differential constitutive relationship of linear viscoelastic, material, a solid-liquid coupling vibration equation for viscoelastic pipe conveying fluid is derived by the D'Alembert's principle. The critical flow velocities and natural frequencies of the cantilever pipe conveying fluid with the Kelvin model (flutter instability) are calculated with the modified finite difference method in the form of the recurrence formula. The curves between the complex frequencies of the first, second and third mode and flow velocity of the pipe are plotted. On the basis of the numerical, calculation results, the dynamic behaviors and stability of the pipe are discussed. It should be pointed out that the delay time of viscoelastic material with the Kelvin model has a remarkable effect on the dynamic characteristics and stability behaviors of the cantilevered pipe conveying fluid, which is a gyroscopic non-conservative system.     

布居数反转的非平衡流模拟计算——对流动激光器性能的研究

Ⅰ.引言40年代不可逆过程热力学已形成一个新的领域。60年代初对扩散、传热、导电、弛豫以及热电、电磁等不可逆过程的各种唯象处理已被综合为统一的体系——非平衡态热力学。由于热机、气体放电、热核聚变以及天体物理、地球物理等方面的研究和发展,流体力学与非平衡态物理学相结合而形成的非平衡流问题的研究已扩展到广泛的领域,流体力学与非平      

Transient analysis of diffusive chemical reactive species for couple stress fluid flow over vertical cylinder

The unsteady natural convective couple stress fluid flow over a semi-infinite vertical cylinder is analyzed for the homogeneous first-order chemical reaction effect. The couple stress fluid flow model introduces the length dependent effect based on the material constant and dynamic viscosity. Also, it introduces the biharmonic operator in the Navier-Stokes equations, which is absent in the case of Newtonian fluids. The solution to the time-dependent non-linear and coupled governing equations is carried out with an unconditionally stable Crank-Nicolson type of numerical schemes. Numerical results for the transient flow variables, the average wall shear stress, the Nusselt number, and the Sherwood number are shown graphically for both generative and destructive reactions. The time to reach the temporal maximum increases as the reaction constant K increases. The average values of the wall shear stress and the heat transfer rate decrease as K increases, while increase with the increase in the Sherwood number.     

CFD SIMULATION OF FLUID CATALYTIC CRACKING IN DOWNER REACTORS

A mathematical model has been developed for the simulation of gas-particle flow and fluid catalytic cracking in downer reactors. The model takes into account both cracking reaction and flow behavior through a four-lump reaction kinetics coupled with two-phase turbulent flow. The prediction results show that the relatively large change of gas velocity affects directly the axial distribution of solids velocity and void fraction, which significantly interact with the chemical reaction. Furthermore, model simulations are carried out to determine the effects of such parameters on product yields, as bed diameter, reaction temperature and the ratio of catalyst to oil, which are helpful for optimizing the yields of desired products. The model equations are coded and solved on CFX4.4.     

CFD SIMULATION OF FLUID CATALYTIC CRACKING IN DOWNER REACTORS

A mathematical model has been developed for the simulation of gas-particle flow and fluid catalytic cracking in downer reactors. The model takes into account both cracking reaction and flow behavior through a four-lump reaction kinetics coupled with two-phase turbulent flow. The prediction results show that the relatively large change of gas velocity affects directly the axial distribution of solids velocity and void fraction, which significantly interact with the chemical reaction. Furthermore, model simulations are carried out to determine the effects of such parameters on product yields, as bed diameter, reaction temperature and the ratio of catalyst to oil, which are helpful for optimizing the yields of desired products. The model equations are coded and solved on CFX4.4.     

渗流力学的新发展

渗流力学是研究流体在多孔介质内运动规律的科学。自1856年法国工程师达西提出线性渗流定律以来,渗流力学一直在向前发展,但最近10多年来的发展更为迅速,其主要动向是:渗流力学应用范围比以前更广;渗流力学理论以较快的速度不断深化;渗流力学研究手段不断实现现代化。本文主要对几项新型渗流——非等温渗流,物理化学渗流,非牛顿流体渗流,生物流体渗流,细观渗流的发展作了综述。      

基于势流理论的内孤立波与立管作用数值研究

内孤立波是一种发生在水面以下的在世界各个海域广泛存在的大幅波浪, 其剧烈的波面起伏所携带的巨大能量对以海洋立管为代表的海洋结构物产生严重威胁, 分析其传播演化过程的流场特征及立管在内孤立波作用下的动力响应规律对于海洋立管的设计具有重要意义. 本文基于分层流体的非线性势流理论, 采用高效率的多域边界单元法, 建立了内孤立波流场分析计算的数值模型, 可以实时获得内孤立波的流场特征. 根据获得的流场信息, 采用莫里森方程计算内孤立波对海洋立管作用的载荷分布. 将内孤立波流场非线性势流计算模型与动力学有限元模型结合来求解内孤立波作用下海洋立管的动力响应特征, 讨论了内孤立波参数、顶张力大小以及内部流体密度对立管动力响应的影响. 发现随着内孤立波波幅的增大, 海洋立管的流向位移和应力明显增大. 由于上层流体速度明显大于下层, 且在所研究问题中拖曳力远大于惯性力, 因此管道顺流向的最大位移发生在上层区域. 顶张力通过改变几何刚度阵的值进而对立管的响应产生明显影响. 对于弱约束立管, 内部流体的密度对管道的流向位移影响较小.