多孔介质中化学反应对非等粘流体混合过程影响的格子Boltzmann研究

利用格子Boltzmann方法和GPU计算技术,在孔隙尺度上模拟多孔介质中包含界面化学反应的粘性指进现象,定量分析化学反应对流体混合的影响.采用单浓度变量的双稳态模型来描述界面反应,而各向同性的多孔介质则通过四参数法生成.研究发现化学反应能减小指进界面厚度,抑制流体的混合,甚至会出现反混合现象,并且随着反应速率的增加,影响越明显.     

微孔中反应流体黏性指进混合过程的数值研究

基于国家数值风洞LBM模块和自编程序，本文在孔隙尺度下对微孔中两种可混溶溶液的黏性指进反应混合过程进行了数值研究。针对两种溶液之间的化学反应以及产物的析出，本文重点研究化学反应速率、固态生成物的成核速率以及聚集速率对黏性指进的影响。研究发现：反应会明显改变黏性指进的界面浓度分布，一方面随着反应速率的增加，指进界面处各组分的变化越来越明显；另一方面由于受到扩散作用的限制，当反应速率增大到一定值后，各组分的平均浓度不再发生变化；产物的生成量与反应速率成正相关，与成核速率和聚集速率成负相关；析出的固相组分的浓度与成核速率和聚集速率成正相关；混合长度主要和反应速率相关，成核速率以及聚集速度的改变对混合长度基本没有影响。     

三维溶质枝晶生长数值模拟

建立了在低Péclet数条件下三维溶质枝晶生长的数值模拟模型.该模型采用Zhu和Stefanescu 提出的溶质平衡方法,即根据固/液界面的平衡浓度和实际浓度之差计算固/液界面演化的驱动力.界面的平衡浓度由界面温度和曲率所确定,实际浓度通过采用有限差分法对溶质扩散控制方程进行数值求解而获得.该方法能够合理定量地描述枝晶从初始的非稳态到稳态的生长过程,并且具有较高的计算效率.为了描述具有不同晶体学取向的三维枝晶生长,提出了一种权值平均曲率算法用于计算固/液界面的曲率,在权值平均曲率的算法中耦合了界面能各向异性的因素.该算法简单易实现,并易于从二维推广到三维系统.为了对模型进行验证,将模拟的枝晶尖端稳态生长数据和理论模型的预测结果进行了比较.结果表明,模拟的Al-2wt%Cu合金枝晶尖端稳态生长速率和半径随过冷度的变化接近于Lipton-Glicksman-Kurz解析模型的预测结果.模拟分析了稳态枝晶尖端的形貌,发现三维枝晶尖端是非轴对称的,以四次对称的方式偏离旋转抛物面.最后,应用所建立的模型模拟出具有发达分枝和不同晶体学取向的三维等轴多枝晶生长形貌. 关键词： 微观组织模拟 溶质枝晶生长 权值平均曲率 三维     

界面能各向异性对定向凝固枝晶生长的影响

采用元胞自动机 (cellular automaton, CA) 模型研究了界面能各向异性对二维定向凝固枝晶生长的影响. 模拟结果显示当晶体的择优生长方向与热流方向一致时, 随着界面能各向异性强度的增大, 凝固组织形态由弱界面能各向异性时的海藻晶转变为强界面能各向异性时的树枝晶. 同时, 界面能各向异性强度会影响稳态枝晶尖端状态的选择, 界面能各向异性越强, 定向凝固稳态枝晶尖端半径越小, 尖端界面前沿的液相浓度和过冷度越小. 稳态枝晶生长的尖端状态选择参数与界面能各向异性强度也存在标度律的指数关系, 而枝晶一次间距则受界面能各向异性强度影响较弱. 当晶体的择优生长方向与热流方向呈-40°夹角时随着界面能各向异性强度的增大, 凝固组织形态由海藻晶逐渐转变为退化枝晶, 后又逐渐转变为倾斜枝晶. 关键词： 元胞自动机 枝晶 界面能各向异性     

非预混燃烧大涡模拟的代数二阶矩模型

一种改进的代数二阶矩燃烧模型被应用在Flame-D的大涡模拟中。代数二阶矩模型能够将化学反应速率的脉动项进行模拟,更准确地计算化学反应速率。在实际燃烧过程中,化学反应速率不仅与当地的温度以及组分浓度等相关,更与反应物的混合速率有密切的关系。研究考虑了亚网格尺度内混合速率对化学反应速率的影响,提出了改进的代数二阶矩模型。通过与相同模拟条件下的涡耗散概念模型的模拟结果以及实验结果进行对比,改进的代数二阶矩模型的准确性得到了验证,同时亚网格化学反应速率的模拟对总化学反应速率的影响也得到了研究。     

考虑复杂化学反应的扩散火焰的数值分析

近十年来,由于对燃烧产生污染和燃烧稳定性的进一步关注,必须认识化学反应的详细机理及各种中间产物在燃烧过程中的作用.为此,人们希望找到对复杂化学反应系统进行分析的理论模型和计算方法.它应当能定量地考虑每个基元反应的速率随温度、组分浓度和压力的变化以及化学反应与流动过程的相互作用.     

自定义标量法模拟丙烷燃烧过程

采用自定义标量法模拟丙烷扩散燃烧,该方法通过把反应组分定义为Fluent程序的自定义标量、化学反应速率作为源项求解质量、动量、能量和组分守恒方程,并用化学反应引起的能量变化修正能量方程.考虑了详细的化学反应机理,整个燃烧反应机理包括27种化学物质(不含N2)和83个基元反应.合理地模拟出了丙烷的燃烧过程,并将火焰的长度、温度、丙烷、氧气以及中间产物的分布与实验数据进行比较.     

固体氧化物燃料电池分布参数模拟研究

固体氧化物燃料电池(SOFC)是一种将储存在燃料和氧化剂中的化学能通过电化学反应的方式直接转换为电能和中高温热能的全固态发电装置,被认为是最接近商业化的先进清洁高效发电技术之一。为了考察SOFC在变负荷条件下的稳态和非稳态特性,在MATLAB模拟环境下,建立了板式SOFC的一维稳态和非稳态分布参数模型,该模型充分考虑了SOFC中传热传质、流体流动、化学与电化学反应等因素相互耦合的影响,模拟分析了SOFC在不同条件下的工作参数变化对系统性能的影响。     

不同前驱体模型对航空煤油扩散火焰碳烟颗粒生成的影响研究

采用不同的航空煤油化学反应机理和碳烟成核模型对气态航空煤油扩散火焰中碳烟颗粒的质量浓度和数量浓度进行预测.分别采用航空煤油详细化学反应机理和简化化学反应机理,结合非预混稳态扩散火焰面模型模拟燃烧反应.分别采用C2H2成核模型(基于乙炔浓度)和PAH成核模型(基于多环芳香烃浓度)预测碳烟颗粒浓度分布.研究结果表明,采用详细化学反应机理和PAH成核模型对碳烟体积分数的预测值与试验值吻合很好.相比于C2H2成核模型,采用PAH成核模型对碳烟体积分数的预测精度显著提升.     

多孔介质内混溶流体间粘性指进现象的格子Boltzmann模型

提出一个模拟多孔介质内混溶流体间粘性指进现象的格子Boltzmann模型.采用双分布函数分别求解压力场和浓度场.在浓度场平衡态分布函数中引入与浓度扩散相关的参数.通过调节参数,使粒子碰撞过程中的松弛时间保持恒定.模拟了粘度相同的流体间的混相驱替问题.不同网格下的模拟结果均与解析解吻合良好,验证了模型的可行性.进一步研究粘度比和贝克莱数(Pe)对粘性指进现象的影响.结果显示,增大粘度比会促进"手指"的增长.当粘度比不变时,存在Pe的临界值.当Pe超过临界值时,"手指"前缘会出现分裂现象.对横向平均浓度场的研究显示,混合区域的长度随时间的变化分为两个阶段,它首先随着t^1/2成线性增长,然后随着t成线性增长.     

Dynamics of immiscible radial viscous fingering: A numerical study

In this paper we undertake a numerical investigation of the dynamics of the interface in the problem of immiscible radial viscous fingering in a Hele-Shaw cell when the areal flow rate is maintained constant. Comparison is made with experimental results to check if there is a need to introduce velocity-dependent boundary conditions and to incorporate the effect of thickness of the film left behind by the moving interface. Some new scaling results are suggested by the numerical data. These data, along with those available from laboratory experiments, provide support for a mean field theory for radial immiscible viscous fingering published recently [Phys. Rev. Lett. 65, 2680 (1990)].     

Fingering of chemical fronts in porous media

The influence of chemical reactions on the hydrodynamical fingering instability is analyzed for miscible systems in porous media. Using a realistic reaction scheme, it is shown that the stability of chemical fronts towards density fingering crucially depends on the width and the speed of the front which are functions of chemical parameters. The major difference between the pure and chemically driven fingering is that, in the presence of chemical reactions, the dispersion curves do not vary in time which has important practical experimental consequences. Good agreement with recent experimental data is found.     

Adhesion and fingering in the lifting Hele-Shaw cell: Role of the substrate

The lifting Hele-Shaw cell (LHSC) is used to study adhesion as well as viscous fingering. In the present paper we report a series of observations of development of the interface for different viscous fluids, both Newtonian and non-Newtonian, in a LHSC operated at a constant lifting force. Glass and perspex are used as the plates in two different sets of experiments. The objectives are 1) to measure the time required to separate the plates as a function of the lifting force and 2) to note the force above which viscous fingering appears. We find that for the Newtonian fluids, the plate separation time follows a universal power law with the lifting force, irrespective of fluid and substrate. The non-Newtonian fluids too, with proper scaling obey the same power law. The appearance of fingering, however, depends on the properties of the fluid as well as the substrate. We suggest a modified form of the capillary number which controls the onset of fingering; this new quantity, termed the “fingering parameter” involves the dielectric constants of the substrate and fluid in addition to the viscosity and surface tension.     

Viscous Fingering Patterns in Ferrofluids

Viscous fingering occurs in the flow of two immiscible, viscous fluids between the plates of a Hele–Shaw cell. Due to pressure gradients or gravity, the initially planar interface separating the two fluids undergoes a Saffman–Taylor instability and develops fingerlike structures. When one of the fluids is a ferrofluid and a perpendicular magnetic field is applied, the labyrinthine instability supplements the usual viscous fingering instability, resulting in visually striking, complex patterns. We consider this problem in a rectangular flow geometry using a perturbative mode-coupling analysis. We deduce two general results: viscosity contrast between the fluids drives interface asymmetry, with no contribution from magnetic forces; magnetic repulsion within the ferrofluid generates finger tip-splitting, which is absent in the rectangular geometry for ordinary fluids.     

Influence of double diffusive effects on miscible viscous fingering

Miscible viscous fingering classically occurs when a less viscous fluid displaces a miscible more viscous one in a porous medium. We analyze here how double diffusive effects between a slow diffusing S and a fast diffusing F component, both influencing the viscosity of the fluids at hand, affect such fingering, and, most importantly, can destabilize the classically stable situation of a more viscous fluid displacing a less viscous one. Various instability scenarios are classified in a parameter space spanned by the log-mobility ratios R(s) and R(f) of the slow and fast component, respectively, and parametrized by the ratio of diffusion coefficients δ. Numerical simulations of the full nonlinear problem confirm the existence of the predicted instability scenarios and highlight the influence of differential diffusion effects on the nonlinear fingering dynamics.     

Fluid-fluid interaction during miscible and immiscible displacement under ultrasonic waves

This paper aims at identifying and analyzing the influence of high-frequency, high-intensity ultrasonic radiation at the interface between immiscible (different types of oils and aqueous solutions) and miscible (different types of oil and solvent) fluids. An extensive set of Hele-Shaw type experiments were performed for several viscosity ratios, and interfacial tension. Fractal analysis techniques were applied to quantify the degree of fingering and branching. This provided a rough assessment of the degree of perturbation generated at the interface when the capillary forces along with the viscous forces are effective. Miscible Hele-Shaw experiments were also presented to isolate the effect of viscous forces. We found that ultrasound acts to stabilize the interfacial front, and that such effect is most pronounced at low viscosity ratios. An erratum to this article is available at .     

Fluid mixing from viscous fingering

Mixing efficiency at low Reynolds numbers can be enhanced by exploiting hydrodynamic instabilities that induce heterogeneity and disorder in the flow. The unstable displacement of fluids with different viscosities, or viscous fingering, provides a powerful mechanism to increase fluid-fluid interfacial area and enhance mixing. Here we describe the dissipative structure of miscible viscous fingering, and propose a two-equation model for the scalar variance and its dissipation rate. Our analysis predicts the optimum range of viscosity contrasts that, for a given Péclet number, maximizes interfacial area and minimizes mixing time. In the spirit of turbulence modeling, the proposed two-equation model permits upscaling dissipation due to fingering at unresolved scales.     

Advanced Study of Unsteady Heat and Chemical Reaction with Ramped Wall and Slip Effect on a Viscous Fluid

This paper investigate the effect of slip boundary condition, thermal radiation, heat source, Dufour number,chemical reaction and viscous dissipation on heat and mass transfer of unsteady free convective MHD flow of a viscous fluid past through a vertical plate embedded in a porous media. Numerical results are obtained for solving the nonlinear governing momentum, energy and concentration equations with slip boundary condition, ramped wall temperature and ramped wall concentration on the surface of the vertical plate. The influence of emerging parameters on velocity,temperature and concentration fields are shown graphically.     

Computer simulation of viscous fingering in a lifting Hele-Shaw cell with grooved plates

We simulate viscous fingering generated by separating two plates with a constant force, in a lifting Hele-Shaw cell. Variation in the patterns for different fluid viscosity and lifting force is studied. Viscous fingering is strongly affected by anisotropy. We report a computer simulation study of fingering patterns, where circular or square grooves are etched on to the lower plate. Results are compared with experiments.