波尔兹曼数字岩芯致密砂岩渗透率研究

致密砂岩渗透率在油气勘探开发、地应力测量及水库地质灾害等领域具有重要意义,但鉴于目前实验技术局限性,温压耦合渗透率测量尚无法通过实验手段实现. 在已有研究工作基础上,首次建立了基于D3Q27数字岩芯模型,并计算了高温压耦合低渗砂岩渗透率问题. 首先,以鄂尔多斯盆地某油田延长组致密砂岩为例,利用X射线CT断层成像技术岩芯获取10μm, 5μm, 2μm分辨率致密砂岩内部结构数据,应用基于量子力学第一性原理的D3Q27格子波尔兹曼数字岩芯模型建立数值模型. 进而,利用自编3DLBM程序分别计算了不同分辨率渗透率随围压(0～200MPa)、孔隙压(0～65MPa)和温度(25℃～180℃)变化规律,通过与Inc AUTOLAB2 000C岩石测试分析系统实验结果对比,验证了程序的可靠性,得到低渗砂岩断层最佳分辨率;最后,在并行CPU-GPU平台上计算了高温高压耦合(0℃～400℃, 0～1.4 GPa)下致密砂岩渗透率值及其各向异性随温压变化规律,并讨论了致密砂岩中水在达到超临界状态后对致密岩石内部结构的影响.     

不同发射深度下导弹水下点火气水流体动力计算

从流体动力角度研究了不同发射深度下，导弹水下点火这一非定常非线性过程。整个系统分为外部水流场、喷管流场和燃气泡流场三个区域加以考虑。水流场采用不可压势流模型，用边界元方法求解；喷管内流场采用非定常一元流动模型，用特征线差分法求解，并设置了激波检测功能；燃气泡采用基于质量和能量守恒的零维计算模型。在时间域中用步进方法实现了三个流场的耦合求解。给出了四种发射深度下的数值计算结果，展示了导弹水下点火的一     

列车与结构动态耦合分析的并行计算方法

在分析结构动态响应时考虑列车与结构的动态耦合作用,采用详细三维有限元模型会带来计算量太大的问题。本文采用并行计算方法,根据列车与结构动态耦合模型的计算特点,设计实现了列车结构耦合均衡的分区算法,并以两个工程应用为例,利用该方法对列车结构三维数值模型进行分区计算,结果表明该分区方法比递归坐标二分法有更好的并行效率。     

带自由表面水流与离散颗粒耦合模型研究

基于考察泥沙运动的细观行为特征,采用离散单元法(DEM)模拟泥沙颗粒运动,结合带自由表面的水动力学计算模型,建立了CFD-DEM耦合数值模型。计算程序开发基于Fortran语言来实现。耦合模型中实现了硬球模型和软球模型两种颗粒碰撞模型,应用范围较广。作为自由表面水流与泥沙颗粒流数值模型的初步研究,在模型建立的基础上,对模型做了基本的验证。分别通过单颗粒静水沉降和混合颗粒群分选两个计算工况,验证了模型的正确性及模拟精度。该耦合模型可进一步丰富带自由表面水流条件下泥沙运动的研究手段。     

多场耦合系统动力学仿真方法研究进展

现代工程系统往往是以复杂结构/机械系统为主体,融合热、流、电磁等若干子系统的多场耦合系统.此类系统动力学建模复杂、计算难度大,给系统动态特性高效精确评估与设计优化带来前所未有的挑战,有关其高效精确动力学仿真方法的研究愈发受到关注.本文详细回顾了复杂工程系统多场耦合动力学仿真方法研究成果和进展,包括:多场耦合动力学建模与数值求解基本策略、网格变形处理方法、耦合数据交换技术、数值计算效率等问题,在此基础上详细讨论了单一和混合不确定性条件下多场耦合系统不确定性分析及可靠性评估方法,以期为相关研究提供有益的借鉴和参考.     

填充床中气体渗流与化学波耦合问题的研究

研究等温和显著气固反应条件下填充床内反应气体浓度(物质)波推进与混合气体渗流的相互作用,指出化学反应对流动的影响包括两个方面,反应过程中混合气体质量的变化和密度的变化。混合气体流动将反过来影响反应进程。分析表明,按耦合模型和非耦合模型得到的速度场完全不同;按耦合模型,反应气体的浓度(物质)波阵面的推进对混合气体的流场有显著影响,因此按耦合模型计算的混合气体流场强烈地依赖于时间;忽略化学反应引起的混合气体密度变化的耦合模型,将导致一个质量消失的汇(或质量生成的源),因此将引起混合气体渗流速度的明显变化,并可能导致物理上不合理的结果;按耦合模型和非耦合模型计算的浓度场也有很大差别;当反应气体与惰性气体摩尔质量相差较大时,不能忽略反应过程中混合气体密度的变化;研究表明对于显著气固反应不能忽略化学反应与气体渗流的相互作用。     

基于龙格库塔法的多输出物理信息神经网络模型

物理信息神经网络(physics-informed neural networks, PINN)由于嵌入了物理先验知识,可以在少量训练数据的情况下获得自动满足物理约束的代理模型,受到了智能科学计算领域的广泛关注.但是, PINN的离散时间模型(PINN-RK)无法同时近似多个物理量相互耦合的偏微分方程系统,限制了其处理复杂多物理场的能力.为了打破这一限制,文章提出了一种基于龙格库塔法的多输出物理信息神经网络(multi-output physics-informed neural networks based on the Runge-Kutta method, MO-PINN-RK), MO-PINN-RK模型在离散时间模型的基础上采用了并行输出的神经网络结构,通过将神经网络划分为多个子网络,建立了多个神经网络输出层.采用不同输出层近似不同物理量的方式, MO-PINN-RK模型不仅可以同时表征多个物理量,而且还能够实现求解偏微分方程系统的目的.另外, MO-PINN-RK克服了PINN离散时间模型仅适用于一维空间的局限性,将其应用范围扩展到了更为普遍的多维空间.为了验证MO-PIN...     

多物理场耦合作用分析的近场动力学理论与方法

广义来说, 近场动力学(peri-dynamics,PD)是假设每个物质点在承受一定范围内的非接触相互作用下,研究整个物理系统演化过程的理论,为涉及非连续和非局部相互作用的问题提供了一个统一的数学框架,具有广泛的适用性.在简要介绍诸多工程对于多物理场模型和数值计算软件的迫切需求后,针对现有商用软件在处理结构非连续演化问题时遇到的瓶颈,引入近场动力学理论和方法. 概述近场动力学固体力学模型,系统阐述近场动力学扩散模型和近场动力学多物理场耦合建模的研究现状和进展,主要涉及电子元器件、电子封装和岩土工程领域的多物理场耦合建模,包括热--力、湿--热--力、热--氧、热--力--氧、力--电、热--电、力--热--电、多孔介质的水--力流固相互作用等非耦合、半耦合与完全耦合模型,强调发展耦合方程数值解法的重要性.最后对扩散问题和多物理场耦合问题的近场动力学理论模型、数值算法和工程应用做进一步展望.      

一种复合材料桨叶气动弹性优化设计方法

考虑到直升机旋翼流场的复杂性,准确的气动力计算需要采用计算流体力学(CFD)方法,而旋翼桨叶由于展弦比较大,几何非线性效应突出,采用计算流体力学和有限元分析(CFD-FEA)方法实现桨叶的单次双向流固耦合分析就需要大量的时间,对优化设计而言,计算量难以承受。针对CFD/FEA耦合计算气动弹性特性的精度和高效性问题,通过PCA提取耦合系统的特征,基于径向基(RBF)神经网络建立气动力降阶模型,代替CFD求解器用于旋翼桨叶的气动弹性分析。将其计算结果与CFD/FEA耦合计算结果进行了对比。研究结果表明,该降阶模型是可行、高效、精确的,可以快速准确地进行复合材料直升机桨叶气动弹性优化设计研究。     

Review of peridynamics for multi-physics coupling modeling

广义来说, 近场动力学(peri-dynamics,PD)是假设每个物质点在承受一定范围内的非接触相互作用下,研究整个物理系统演化过程的理论,为涉及非连续和非局部相互作用的问题提供了一个统一的数学框架,具有广泛的适用性.在简要介绍诸多工程对于多物理场模型和数值计算软件的迫切需求后,针对现有商用软件在处理结构非连续演化问题时遇到的瓶颈,引入近场动力学理论和方法. 概述近场动力学固体力学模型,系统阐述近场动力学扩散模型和近场动力学多物理场耦合建模的研究现状和进展,主要涉及电子元器件、电子封装和岩土工程领域的多物理场耦合建模,包括热--力、湿--热--力、热--氧、热--力--氧、力--电、热--电、力--热--电、多孔介质的水--力流固相互作用等非耦合、半耦合与完全耦合模型,强调发展耦合方程数值解法的重要性.最后对扩散问题和多物理场耦合问题的近场动力学理论模型、数值算法和工程应用做进一步展望.     

无单元伽辽金法的并行计算

对无单元伽辽金法的并行计算进行了详细研究,并将其应用于弹性动力学问题。使用并行桶搜索算法进行节点搜索,使用并行几何搜索算法进行样点搜索,讨论了移动最小二乘MLS(Moving Least Squares)形函数及其导数的并行计算和方程组的并行求解,并利用多层图形划分实现负载平衡。最后给出了并行无单元伽辽金法应用于弹性动力学的计算流程和实例。计算结果表明无单元伽辽金法具有很高的并行性和很好的并行效率,对其进行并行计算具有非常重要的意义。     

A fire field model implemented in a parallel computing environment

This paper describes the implementation of a fire field model in the parallel computing environment offered by multiple transputers. The fire model is built into the general purpose SIMPLE-based CFD code HARWELL-FLOW3D. The technique of domain decomposition has been applied tb convert the conventional serial version of FLOW3D into a code capable of efficiently utilizing an arbitrary number of transputers. Fire simulations consisting of up to 24 000 computational cells are performed on parallel systems with up to 15-processors. The run time for this simulation has been reduced from over 4 days on a single processor to just over 8 h on the 15-processor system. An interactive graphics system has also been developed which runs in parallel with the main computations.     

Quasi-real-time simulation of rotating drum using discrete element method with parallel GPU computing

Real-time simulation of industrial equipment is a huge challenge nowadays. The high performance and fine-grained parallel computing provided by graphics processing units (GPUs) bring us closer to our goals. In this article, an industrial-scale rotating drum is simulated using simplified discrete element method (DEM) without consideration of the tangential components of contact force and particle rotation. A single GPU is used first to simulate a small model system with about 8000 particles in real-time, and the simulation is then scaled up to industrial scale using more than 200 GPUs in a 1D domain-decomposition parallelization mode. The overall speed is about 1/11 of the real-time. Optimization of the communication part of the parallel GPU codes can speed up the simulation further, indicating that such real-time simulations have not only methodological but also industrial implications in the near future.     

A mesh patching method for finite volume modelling of shallow water flow

A new mesh‐patching model is presented for shallow water flow described by the 2D non‐linear shallow water (NLSW) equations. The mesh‐patching model is based on AMAZON, a high‐resolution NLSW engine with an improved HLLC approximate Riemann solver. A new patching algorithm has been developed, which not only provides improved spatial resolution of flow features in particular parts of the mesh, but also simplifies and speeds up the (structured) grid generation process for an area with complicated geometry. The new patching technique is also compatible with increasingly popular parallel computing and adaptive grid techniques. The patching algorithm has been tested with moving bores, and results of test problems are presented and compared to previous work. Copyright © 2005 John Wiley & Sons, Ltd.     

PARALLEL FINITE ELEMENT ANALYSIS OF HIGH FREQUENCY VIBRATIONS OF QUARTZ CRYSTAL RESONATORS ON LINUX CLUSTER

Quartz crystal resonators are typical piezoelectric acoustic wave devices for frequency control applications with mechanical vibration frequency at the radio-frequency （RF） range. Precise analyses of the vibration and deformation are generally required in the resonator design and improvement process. The considerations include the presence of electrodes, mountings, bias fields such as temperature, initial stresses, and acceleration. Naturally, the finite element method is the only effective tool for such a coupled problem with multi-physics nature. The main challenge is the extremely large size of resulted linear equations. For this reason, we have been employing the Mindlin plate equations to reduce the computational difficulty. In addition, we have to utilize the parallel computing techniques on Linux clusters, which are widely available for academic and industrial applications nowadays, to improve the computing efficiency. The general principle of our research is to use open source software components and public domain technology to reduce cost for developers and users on a Linux cluster. We start with a mesh generator specifically for quartz crystal resonators of rectangular and circular types, and the Mindlin plate equations are implemented for the finite element analysis. Computing techniques like parallel processing, sparse matrix handling, and the latest eigenvalue extraction package are integrated into the program. It is clear from our computation that the combination of these algorithms and methods on a cluster can meet the memory requirement and reduce computing time significantly.     

On the mechanical modeling of anisotropic biological soft tissue and iterative parallel solution strategies

Biological soft tissues appearing in arterial walls are characterized by a nearly incompressible, anisotropic, hyperelastic material behavior in the physiological range of deformations. For the representation of such materials we apply a polyconvex strain energy function in order to ensure the existence of minimizers and in order to satisfy the Legendre–Hadamard condition automatically. The 3D discretization results in a large system of equations; therefore, a parallel algorithm is applied to solve the equilibrium problem. Domain decomposition methods like the Dual-Primal Finite Element Tearing and Interconnecting (FETI-DP) method are designed to solve large linear systems of equations, that arise from the discretization of partial differential equations, on parallel computers. Their numerical and parallel scalability, as well as their robustness, also in the incompressible limit, has been shown theoretically and in numerical simulations. We are using a dual-primal FETI method to solve nonlinear, anisotropic elasticity problems for 3D models of arterial walls and present some preliminary numerical results.     

Portable data‐parallel surface reconstruction on a uniform rectilinear grid

With the increasing heterogeneity and on‐node parallelism of high‐performance computing hardware, a major challenge is to develop portable and efficient algorithms and software. In this work, we present our implementation of a portable code to perform surface reconstruction using NVIDIA's Thrust library. Surface reconstruction is a technique commonly used in volume tracking methods for simulations of multimaterial flow with interfaces. We have designed a 3D mesh data structure that is easily mapped to the 1D vectors used by Thrust and at the same time is simple to use and uses familiar data structure terminology (such as cells, faces, vertices, and edges). With this new data structure in place, we have implemented a piecewise linear interface reconstruction algorithm in 3 dimensions that effectively exploits the symmetry present in a uniform rectilinear computational cell. Finally, we report performance results, which show that a single implementation of these algorithms can be compiled to multiple backends (specifically, multi‐core CPUs, NVIDIA GPUs, and Intel Xeon Phi processors), making efficient use of the available parallelism on each. We also compare performance of our implementation to a legacy FORTRAN implementation in Message Passing Interface (MPI) and show performance parity on single and multi‐core CPU and achieved good parallel speed‐ups on GPU. Our research demonstrates the advantage of performance portability of the underlying data‐parallel programming model.     

3-D Fracture Propagation Simulation and Production Prediction in Coalbed

NｏｍｅｎｃｌａｔｕｒｅＡ(ｘ,ｔ)—Ｃｒｏｓｓ＿ｓｅｃｔｉｏｎａｒｅａｏｆｆｒａｃｔｕｒｅａｔｔｉｍｅｔａｎｄｓｉｔｅｘ,ｍ2 ;Ｂｇ,Ｂｗ —Ｃｏａｌｂｅｄｇａｓａｎｄｗａｔｅｒｖｏｌｕｍｅｆａｃｔｏｒ,ｄｉ ｍｅｎｓｉｏｎｌｅｓｓ;ｂ—Ｌａｎｇｍｕｉｒｆａｃｔｏｒ,(ＭＰａ) - 1 ;Ｃｇ,ＣＲ,Ｃｗ —Ｃｏａｌｂｅｄｇａｓ,ｒｏｃｋａｎｄｗａｔｅｒｃｏｍ ｐｒｅｓｓｉｂｉｌｉｔｙｃｏｅｆｆｉｃｉｅｎｔ,(ＭＰａ) - 1 ;Ｃｐ —Ｓａｎｄｒａｔｉｏ ,ｄｉｍｅｎｓｉｏｎｌｅｓｓ;Ｃｔ(ｘ ,ｔ)—Ｔｏｔａｌｆｒａｃｔｕｒｉｎ…     

并行求解周期性三对角方程组的谢-莫方法

针对工程计算中出现的周期性三对角线性方程组,本文利用谢尔曼-莫里森方法修正周期性三对角线性方程组,将其转化为求解两个三对角线性方程组的问题,结合分治法进行并行计算。在对周期性三对角线性方程组并行化分析的基础上,在集群系统中利用可移植消息传递标准MPI进行并行程序设计。结果表明,使用该方法求解周期性三对角线性方程组,既具有数值稳定性,又收到较好的并行效果。     

基于集群的并行有限元分析研究

利用面向对象的方法实现了基于集群的并行有限元分析.首先介绍构建高性能计算集群的系统分析新方法,并针对建成集群系统的特点,对并行消息传递库MPI进行了并行语义分析,进而对其主要的实现函数进行面向对象的重构而建立起OO消息传递库,以此为基础分别利用两种不同的策略实现了面向对象的并行有限元分析(OOParaFEA: Object-Oriented Parallel FEA),一种是通过对传统的基于域分解方法的并行PCG算法进行改造而在集群平台实现,另一种是对已有的基于子结构方法的串行有限元分析程序加以扩展,加入系统方程组并行求解器而达到有限元分析并行化的目的.多个分析算例表明,基于集群进行面向对象的并行有限元分析可以有效提高计算效率,为进一步的网络化CAD/CAE研究奠定良好基础.