多层介质膜光栅用高反射镜的严格耦合波分析

基于严格耦合波和增强透射矩阵的方法,提出了一种数值稳定的求解多层介质膜光学特性的机理模型.利用该模型计算了多层介质膜光栅用高反射镜的优化设计膜系.使用电子束热蒸发方式制备的多层介质膜光谱特性和理论设计的结果符合得很好,该严格耦合波模型是分析介质膜光学特性有效的稳定的数值求解方法.     

致密多孔介质中气体渗流的格子Boltzmann模拟

为了能用格子Boltzmann方法来正确地刻画致密多孔介质中微尺度流动问题,对单通道模型进行推广,并将其应用于孔隙群里气体渗流的数值模拟.通过对部分具有代表性的多孔介质中的实际流动问题进行模拟,研究渗透率与平均压力和Knudsen数之间的相互关系.基于理论分析及相关文献中的试验结果,验证模拟结果,为用格子Boltzmann方法深入研究气体渗流问题奠定基础.     

格子Boltzmann模型平衡分布边界条件和多棱柱排列渗流流场的数值模拟

介绍了适用于多种流场数值模拟的无滑动格子Boltzmann平衡分布边界条件，这一边界条件是以Bounce-Back方法为基础并满足质量、动量守恒的准则.数值计算结果表明平衡分布边界条件克服了Bounce-Back方法在边界上所产生的滑动速度误差效应.利用平衡分布边界条件数值模拟了由棱柱形充填粒子构成的微尺度渗流流场中的Darcy-Forchheimer方程，通过与Lee 和Yang的数值结果比较，该预测结果是足够可靠的. 关键词： 平衡分布边界条件 渗流介质 Darcy-Forchheimer阻力     

裂缝性油藏水平井产能计算方法

根据裂缝性油藏水平井渗流特征,将渗流区域等效为近井区域和远井区域.远井区域视为连续均匀介质.在近井区域,基于势叠加、镜像反映原理和微元线汇思想,建立考虑天然裂缝和水平井生产段共同渗流的耦合模型.组合两个流动区域得到裂缝性油藏水平井产能评价模型,并进行求解.实例计算结果表明,由于考虑了裂缝之间的干扰,用建立的模型计算的水平井产能小于用理论公式计算的水平井的产能,与实际产能比较接近;随着裂缝密度增加,产量逐渐增加,但增加的幅度降低.     

一种新的模拟渗流运动的数值方法

根据格子Boltzmann方法及相关理论，建立了一个新的模拟渗流运动的数值模型，所得模型没有在边界上采取相应平均措施，同时还避免了一些非物理副产品的出现-实例计算数值结果与精确解符合较好，证明模型可靠- 关键词： 渗流 格子Boltzmann方法 数值模型     

Ghost Fluid方法与双介质可压缩流动计算

应用带有Isobaric修正的GhostFluid方法配合LevelSet方法计算可压缩双介质无粘流动.该方法可以消除计算流体界面时所产生的数值跳动和耗散,且编程上比界面跟踪法简单.应用WENO格式数值求解欧拉方程和LevelSet方程,对由刚性气体状态方程所支配的一二维双介质流动进行数值计算,得到了分辨率较高的计算结果.     

基于等效介质二维光子晶体平面波展开法的理论研究

应用等效介质理论,提出了一种计算二维光子晶体的新平面波展开法,数值模拟了圆介质柱在空气中正方排列结构的能带,模拟结果与常规的平面波计算结果符合的很好.     

金属圆筒内爆轰波对碰效应的欧拉数值模拟

运用基于Youngs技术的多介质欧拉方法及应用该方法研制的二维多介质弹塑性流体动力学欧拉程序meph2d对金属圆筒内爆轰波相互作用过程进行了数值模拟并获得成功，计算图像与实验照片吻合，计算结果与实验测量结果符合较好。说明采用多介质欧拉方法可以对此类爆轰波相互作用及爆轰驱动问题进行有效的数值模拟。     

耦合的van der Pol振子的极限环幅值控制

设计反馈控制器，对一类耦合的van der Pol振子的极限环幅值进行控制. 用近似解析方法求出了控制系统的幅值的控制方程, 得到了控制参数与极限环幅值的函数关系, 使系统的振幅能按需求得到有效地调整. 通过数值计算绘制了在不同控制参数下, 系统响应的时间历程曲线和极限环. 近似解析方法计算得到的结果与数值计算进行了比较，两者是符合的. 这一方法也可以推广应用到其他耦合的van der Pol振子.     

降膜结晶过程的数值模拟

本文采用数值方法对在工业上有着广泛应用前景的降膜结晶分离过程进行了模拟，提出了一套描述降膜结晶过程的降膜流动、传热、传质及结晶相变相互耦合的以控制容积为基础的有限差分的湍流数值计算模型和算法。本模型的计算结果与以前简化模型的计算结果相比，更符合实验值．     

XFEM modeling of hydraulic fracture in porous rocks with natural fractures

Hydraulic fracture (HF) in porous rocks is a complex multi-physics coupling process which involves fluid flow, diffusion and solid deformation. In this paper, the extended finite element method (XFEM) coupling with Biot theory is developed to study the HF in permeable rocks with natural fractures (NFs). In the recent XFEM based computational HF models, the fluid flow in fractures and interstitials of the porous media are mostly solved separately, which brings difficulties in dealing with complex fracture morphology. In our new model the fluid flow is solved in a unified framework by considering the fractures as a kind of special porous media and introducing Poiseuille-type flow inside them instead of Darcy-type flow. The most advantage is that it is very convenient to deal with fluid flow inside the complex fracture network, which is important in shale gas extraction. The weak formulation for the new coupled model is derived based on virtual work principle, which includes the XFEM formulation for multiple fractures and fractures intersection in porous media and finite element formulation for the unified fluid flow. Then the plane strain Kristianovic-Geertsma-de Klerk (KGD) model and the fluid flow inside the fracture network are simulated to validate the accuracy and applicability of this method. The numerical results show that large injection rate, low rock permeability and isotropic in-situ stresses tend to lead to a more uniform and productive fracture network.     

A New Numerical Solution of Fluid Flow in Stratigraphic Porous Media

A new numerical technique based on a lattice-Boltzmann method is presented for analyzing the fluid flow in stratigraphic porous media near the earth's surface. The results obtained for the relations between porosity, pressure, and velocity satisfy well the requirements of stratigraphic statistics and hence are helpful for a further study of the evolution of fluid flow in stratigraphic media.     

Tortuosity for streamlines in porous media

An analysis of tortuosity for streamlines in porous media is presented by coupling the circle and square models. It is assumed that some particles in porous media do not overlap and that fluid in porous media is incompressible. The relationship between tortuosity and porosity is attained with different configurations by using a statistical method. In addition, the tortuosity fractal dimension is expressed as a function of porosity. Those correlations do not include any empirical constant. The percolation threshold and tortuosity fractal dimension threshold of porous media are also presented as: c = 0.32, D T c = 1.07. The predicted correlations of the tortuosity and the porosity agree well with the existing experimental and simulated results.     

三维有序排列多孔介质对流换热的数值研究

利用CFD软件数值研究了颗粒三维有序堆积多孔介质的对流换热问题. 采用颗粒直径分别为14 mm,9.4 mm和7 mm的球形颗粒有序排列构成多孔介质骨架,在多孔骨架的上方有一恒热流密度的铜板. 采用流固耦合的方法研究了槽通道内温度分布和局部对流换热系数的分布以及对流换热的影响因素. 研究结果表明:热渗透的厚度和温度边界层的厚度在流动方向上逐渐增大,并且随流量的增加而减小;当骨架的导热系数比较高时,对流换热随颗粒直径的减小而略有增大;对流换热系数随聚丙烯酰胺溶液浓度的增大而减小,黏性耗散减弱了对流换热. 关键词： 多孔介质 温度场 局部对流换热系数 数值模拟     

A Lattice Boltzmann Model for Fluid-Solid Coupling Heat Transfer in Fractal Porous Media

We report a lattice Boltzmann model that can be used to simulate fluid-solid coupling heat transfer in fractal porous media. A numerical simulation is conducted to investigate the temperature evolution under different ratios of thermal conductivity of solid matrix of porous media to that of fluid. The accordance of our simulation results with the solutions from the conventional CFD method indicates the feasibility and the reliability for the developed lattice Boltzmann model to reveal the phenomena and rules of fluid-solid coupling heat transfer in complex porous structures.     

Numerical simulation for separation of multi-phase immiscible fluids in porous media

Based on a lattice Boltzmann method and general principles of porous flow, a numerical technique is presented for analysing the separation of multi-phase immiscible fluids in porous media. The total body force acting on fluid particles is modified by axiding relative permeability in Nithiarasu＇s expression with an axiditional surface tension term. As a test of this model, we simulate the phase separation for the case of two immiscible fluids. The numerical results show that the two coupling relative permeability coefficients K12 and K21 have the same magnitude, so the linear flux-forcing relationships satisfy Onsager reciprocity. Phase separation phenomenon is shown with the time evolution of density distribution and bears a strong similarity to the results obtained from other numerical models and the flows in sands. At the same time, the dynamical rules in this model are local, therefore it can be run on massively parallel computers with well computational efficiency.     

A new method for simulation and analysis of displacement of fluids in porous media

Using the principle of diffusion-limited aggregation(DLA),a new model is introduced to simulate the displacement of one fluid by another in porous media.The results agree with experiments.apparently they do not leave out film-flow phenomena.Simultaneously,we also present a new numerical method to treat our results by the lattice Boltzmann method(LBM),All these will be helpful for analysing similar subjects.     

随机多孔介质的蒙特卡罗重构与渗流特性模拟

根据多孔介质微观结构的分形尺度标度特征,采用蒙特卡罗方法分别重构随机多孔介质的微观颗粒和孔隙结构,并基于分形毛管束模型研究多尺度多孔介质的气体渗流特性,建立多孔介质微观结构和宏观渗流特性的定量关系。结果表明：分形蒙特卡罗重构的多孔介质微细结构接近真实介质结构,气体渗流特性的计算结果与格子玻尔兹曼模拟数据较为吻合; 多孔介质气体渗透率随着克努森数的增加而增大,孔隙分形维数对于气体渗流的微尺度效应具有显著影响,而迂曲度分形维数对于表观渗透率和固有渗透率的比值影响可以忽略。提出的分形蒙特卡罗方法具有收敛速度快且计算误差与维数无关的优点,有利于深入理解多尺度多孔介质的渗流机理。     

Direct measurement of porous media local hydrodynamical permeability using gas MRI.

The concept of hydraulic permeability is at the core of modeling single phase or multi-phase flow in heterogeneous porous media, as it is the spatial distribution of the permeability that primarily governs the behavior of fluid flow in the medium. To date, the modeling of fluid flow in porous media has been hampered by poor estimates of local permeability. Magnetic Resonance Imaging is well known for its ability to measure non-invasively the local density and flow rate of different fluids saturating porous media [1,2]. In this paper we demonstrate the first non-invasive method for the direct measurement of a single projection of the local permeability tensor of a porous medium using gas-phase MRI. The potential for three-dimensional imaging of the medium permeability is also discussed. The limitations of the method are listed and results are presented in a model porous medium as well as in a real oil reservoir rock.