一维非均匀Kerr介质的传输矩阵算法

研究了一维非均匀Kerr介质的子层逆向递推传输矩阵算法,给出了数值实现方案.此方法可以方便而精确地模拟非均匀Kerr介质的传输行为.应用此算法模拟一维Kerr缺陷光子晶体的光学双稳开关过程时, 得出了和非线性时域有限差分方法相一致的结果.     

三维非均匀不稳定渗流方程的自适应网格粗化算法

将渗透率自适应网格技术应用于三维非均匀不稳定渗流方程的网格粗化算法中,在渗透率或孔隙度变化异常区域自动采用精细网格,用直接解法求解渗透率或孔隙度变化异常区域的压强分布,在其它区域采用不均匀网格粗化的方法计算,即在流体流速大的区域采用精细网格.用该方法计算了三维非均匀不稳定渗流场的压降解,结果表明三维非均匀不稳定渗流方程的三维非均匀自适应网格粗化算法的解在渗透率或孔隙度异常区的压强分布规律与采用精细网格的解非常逼近,在其它区域压强分布规律与粗化算法的解非常逼近,计算速度比采用精细网格提高100多倍.     

表征体元尺度渗流的离散统一动理学格式

结合表征体元尺度的通用渗流模型,提出离散统一动理学格式（DUGKS）渗流方法,分别用均匀网格和非均匀网格计算二维Poiseuille、Couette、方腔流等经典渗流问题,检验DUGKS渗流方法的有效性和非均匀网格应用的优势,将DUGKS渗流方法应用到裂缝系统中.     

二维高速碰撞问题欧拉数值模拟的混合网格计算

提出了适用于二维平面或轴对称多介质流体力学两步欧拉数值方法中输运计算的混合网格界面处理.在一个混合网格中,将界面近似看作直线.整个方法分为3步:①用混合网格周围的8个网格的介质面积份额确定界面的法线方向;②用混合网格的介质面积份额或体积份额确定界面的直线方程;③用此直线方程求出通过网格边界的流.给出了用此方法所做的测试、数值计算及与其它算法的比较.     

气体动理学格式模拟多孔介质流动的可行性

将气体动理学格式（GKS）拓展到模拟多孔介质内的低速渗流,并检验在孔隙尺度上模拟不可压缩低速流动的可行性与有效性.结果表明：GKS具有二阶空间精度,能够较精确地计算多孔介质的渗透率;相比于单松弛格子玻尔兹曼方法,GKS能够精确实现壁面无滑移边界条件,从而正确反映渗透率与黏性无关的特性;对于Berea砂岩切片结构中的复杂流动,模拟结果与实验吻合较好,能较精确地计算渗透率.给出GKS模拟达西渗流的马赫数选取准则,为研究多孔介质流动提供新的工具.     

多孔介质中盐指现象的数值模拟

运用基于杂交网格的高精度数值方法研究了多孔介质中的盐指现象.该算法将基于边界拟合坐标下的高精度有限差分法和高精度的泊松方程快速求解器有效地结合在一起,从而达到提高整体的计算精度、计算效率和稳定性的目的.通过比较不同孔隙率的多孔介质对盐指对流的传热传质效应的影响,发现在标准孔隙率较低的多孔介质中,盐度扩散的速度明显比热扩散的速度快,盐指很快触及上下壁面,使得上下层的盐度梯度迅速减小,这是与非多孔介质具有明显差异之处. 关键词： 多孔介质 双扩散对流 盐指     

地层渗透参数反演计算

利用渗流源网格压力(及井底压力)对网格渗透率导数关系(称为渗透率压力场敏感系数)实现由三维渗流方程压强数值解计算井底压力对网格渗透参数导数的计算.基于渗透率压力场敏感系数的计算结果构造地层渗透参数反演方程,由多井系统井压监测数据实现反演地层渗透参数的计算.算例显示:反演方法对均匀(或分块均匀)油藏模型的计算取得了较好的结果,反演迭代10次后计算结果与模型参数的最大相对误差收敛于1%之内.该方法对油藏动态模拟研究及实现油藏参数动态跟踪具有实际意义.     

多网格谱元法及其在高性能计算中的应用*

针对传统谱元法在每个单元内只能存在单一均匀介质,应用在复杂非均匀介质的波传播模拟中可能造成极大计算规模的问题,发展了多网格谱元法。该方法在谱元法单元内引入独立的辅助网格,用于精细描述单元内的介质和外力分布变化,在较稀疏的主网格上进行波场的求解。基于声波和弹性波方程推导了多网格谱元法公式,并对几种典型模型进行了波场的数值模拟。与传统谱元法的对比结果表明,此算法在复杂非均匀介质的弹性波传播模拟中可以利用较少的网格点数达到不低于传统算法的精度。此外,实现了并行化的多网格谱元法,获得了较好的并行效率。     

大孔穴多孔介质内湍动特性研究

建立了用于过滤燃烧数值模拟的大孔穴多孔介质3-D均质网络模型。利用该模型研究了四组不同参数多孔介质内压力梯度与渗流速度之间的变化规律。数值模拟与实验结果吻合得很好。数值计算表明:相对于孔隙率,孔径亦即渗透率对压力梯度的影响更为显著;入口效应仅局限在距入口界面2～3个单元长度内,几何参数和渗流速度对其影响可以忽略;多孔介质的存在抑制了湍动能的发展,在受入口效应影响的区域,湍动能衰减得比较剧烈,之后变化趋缓,最终与耗散率形成动态平衡。     

计算中子增殖率的时间强迫碰撞方法

为获得更多粒子在特定时间步内的碰撞信息,给出3种解决介质非均匀的时间强迫碰撞抽样目的:离散化抽样目的,强截面抽样目的和均匀时间抽样目的.将强迫碰撞目的由只适合在各自均匀的几何区域内分别强迫碰撞,扩展为也适合在特定时间步内中子可能经过任意非均匀介质的强迫碰撞.通过一个简单的增殖系统的中子增殖率计算问题,对时间强迫碰撞的3种目的和直接模拟目的进行了计算比较.结果表明,3种时间强迫碰撞目的的效率相差不大,但明显优于直接模拟目的,当步长从2×10^-5缩至10^-6μs,FOM(计算效率或称优度)约高出2至4个量级.     

裂缝性介质多尺度深度学习模型

结合人工神经网络建立裂缝介质多尺度深度学习流动模型.基于一套粗网格和一套细网格,通过在粗网格上训练数据,多尺度神经网络能够以较少的自由度训练出准确的神经网络.并在粗网格上通过求解局部流动问题获得多尺度基函数,结合神经网络进一步得到精细网格的解.基于离散裂缝的流动方程可视为多层网络,网络层数依赖于求解时间步数.阐述裂缝介质多尺度机器学习数值计算格式的建立,介绍如何使用多尺度算法构建离散裂缝模型的多尺度基函数,并采用超样本技术进一步提高计算准确性.数值结果表明,多尺度有限元算法与机器学习结合是一种有效的流体流动模拟算法.     

Numerical methods for multiscale transport equations and application to two-phase porous media flow

We discuss numerical methods for linear and nonlinear transport equations with multiscale velocity fields. These methods are themselves multiscaled in nature in the sense that they use macro and micro grids, multiscale test functions. We demonstrate the efficiency of these methods and apply them to two-phase flow in heterogeneous porous media.     

A high-resolution gas-kinetic scheme with minimized dispersion and controllable dissipation reconstruction

In order to simulate multiscale problems such as turbulent flows effectively, the high-order accurate reconstruction based on mini- mized dispersion and controllable dissipation (MDCD) is implemented in the second-order accurate gas-kinetic scheme (GKS) to improve the accuracy and resolution. MDCD is firstly extended to non-uniform grids through the modification of dissipation and dispersion coefficients for uniform grids based on the local stretch ratio. Remarkable improvements in accuracy and resolution are achieved on general grids. Then a new scheme, MDCD-GKS is constructed, with the help of MDCD reconstruction, not only for conservative variables, but also for their gradients. MDCD-GKS shows good accuracy and efficiency in typical numerical tests. MDCD-GKS is also coupled with the improved delayed detached-eddy simulation (IDDES) hybrid model and applied in the fine simulation of turbulent flow around a cylinder, and the prediction is in good agreement with experiments when using the relatively coarse grid. The high accuracy and resolution of the developed GKS guarantee its high efficiency in practical applications.     

Advances in multi-domain lattice Boltzmann grid refinement

Grid refinement has been addressed by different authors in the lattice Boltzmann method community. The information communication and reconstruction on grid transitions is of crucial importance from the accuracy and numerical stability point of view. While a decimation is performed when going from the fine to the coarse grid, a reconstruction must performed to pass form the coarse to the fine grid. In this context, we introduce a decimation technique for the copy from the fine to the coarse grid based on a filtering operation. We show this operation to be extremely important, because a simple copy of the information is not sufficient to guarantee the stability of the numerical scheme at high Reynolds numbers. Then we demonstrate that to reconstruct the information, a local cubic interpolation scheme is mandatory in order to get a precision compatible with the order of accuracy of the lattice Boltzmann method.These two fundamental extra-steps are validated on two classical 2D benchmarks, the 2D circular cylinder and the 2D dipole–wall collision. The latter is especially challenging from the numerical point of view since we allow strong gradients to cross the refinement interfaces at a relatively high Reynolds number of 5000. A very good agreement is found between the single grid and the refined grid cases.The proposed grid refinement strategy has been implemented in the parallel open-source library Palabos.     

Nondissipative and energy-stable high-order finite-difference interface schemes for 2-D patch-refined grids

A class of finite-difference interface schemes suitable for two-dimensional cell-centered grids with patch-refinement and step-changes in resolution is presented. Grids of this type are generated by adaptive mesh refinement methods according to resolution needs dictated by the physics of the problem being modeled. For these grids, coarse and fine nodes are not aligned at the mesh interfaces, resulting in hanging nodes. Three distinct geometries are identified at the interfaces of a domain with interior patch-refinement: edges, concave corners and convex corners. Asymptotic stability in time of the numerical scheme is achieved by imposing a summation-by-parts condition on the interface closure, which is thus also nondissipative. Interface stencils corresponding to an explicit fourth-order finite-difference scheme are presented for each geometry. To preserve stability, a reduction in local accuracy is required at the corner geometries. It is also found that no second-order accurate solution exists that satisfies the summation-by-parts condition. Tests using the 2-D scalar advection equation and an inviscid compressible vortex support the stability and accuracy of these stencils for both linear and nonlinear problems.     

Lattice Boltzmann method on quadtree grids for simulating fluid flow through porous media: A new automatic algorithm

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a synthetic porous medium using the LB method and on the quadtree grid structure. To this end, the LB method was used on nonuniform grids coupled with a technique for image reconstruction which resulted in the quadtree grids for simulation of fluid flow through porous media. Accuracy and efficiency of this algorithm is compared against the conventional LB method based on uniform grids. While the decrease in computational time in the proposed LB method on nonuniform grids is found to be significant regarding the size of the initial and reconstructed images, the same level of accuracy is obtained when compared with the conventional LB method on uniform grids.     

求解对流-扩散-反应问题的改进有限积分法

将求解偏微分方程的有限积分法应用于对流-扩散-反应问题,发现对于非对流占优的对流扩散问题,有限积分法的精度比QUICK法高一个数量级,比传统的有限体积法高两个数量级.处理对流占优的对流-扩散-反应问题时,对流项的离散时引进加权参数,通过调节该参数反映输运的方向性.结果表明这种改进的有限积分法的精度比传统的有限体积法至少高四个数量级,同时明显改进了原来的有限积分法的精度和稳定性.对于对流占优的对流-扩散-反应问题,即使采用粗网格,计算结果也未出现非物理振荡现象,表明改进的有限积分法具有很好的稳定性.     

Numerical model of two-dimensional heterogeneous combustion in porous media under natural convection or forced filtration

A novel mathematical model and original numerical method for investigating the two-dimensional waves of heterogeneous combustion in porous media are proposed and described in detail. The mathematical model is constructed within the framework of the model of interacting interpenetrating continua and includes equations of state, continuity, momentum conservation and energy for solid and gas phases. Combustion, considered in the paper, is due to the exothermic reaction between fuel in the porous solid medium and oxidiser contained in the gas flowing through the porous object. The original numerical method is based on a combination of explicit and implicit finite-difference schemes. A distinctive feature of the proposed model is that the gas velocity at the open boundaries (inlet and outlet) of the porous object is unknown and has to be found from the solution of the problem, i.e. the flow rate of the gas regulates itself. This approach allows processes to be modelled not only under forced filtration, but also under free convection, when there is no forced gas input in porous objects, which is typical for many natural or anthropogenic disasters (burning of peatlands, coal dumps, landfills, grain elevators). Some two-dimensional time-dependent problems of heterogeneous combustion in porous objects have been solved using the proposed numerical method. It is shown that two-dimensional waves of heterogeneous combustion in porous media can propagate in two modes with different characteristics, as in the case of one-dimensional combustion, but the combustion front can move in a complex manner, and gas dynamics within the porous objects can be complicated. When natural convection takes place, self-sustaining combustion waves can go through the all parts of the object regardless of where an ignition zone was located, so the all combustible material in each part of the object is burned out, in contrast to forced filtration.