数值模拟方法进展
湍流燃烧数值模拟的研究与进展
湍流燃烧数值模拟的研究与进展湍流燃烧是指在燃烧过程中，燃料与氧化剂在湍流的条件下相遇和反应。
湍流燃烧数值模拟是一种通过计算机模拟湍流燃烧过程的方法，可以提供燃烧器内部的流场和温度分布等信息，对于燃烧器的设计和优化具有重要的意义。
本文将对湍流燃烧数值模拟的研究与进展进行探讨。
首先，湍流模型的选择是湍流燃烧数值模拟的一个关键问题。
湍流现象十分复杂，需要选择适当的湍流模型来模拟湍流流动。
常用的湍流模型有雷诺平均应力模型(RANS)和大涡模拟(LES)。
RANS是一种将湍流场分为均匀部分和涡旋部分的统计方法，适用于模拟湍流较为稳定的情况；而LES则能模拟较为精细的湍流结构，但计算量较大。
根据具体问题的复杂程度和计算资源的限制，选择适当的湍流模型具有重要意义。
其次，化学反应模型的建立是湍流燃烧数值模拟的另一个关键问题。
燃烧过程中涉及到多种化学反应，需要建立合适的化学反应模型来描述燃烧反应。
常见的化学反应模型有简化化学反应模型和详细化学反应模型。
简化化学反应模型基于简化的反应机理，计算速度较快；而详细化学反应模型则基于包含大量反应步骤的反应机理，计算速度较慢但结果更精确。
根据具体问题的要求和计算资源的限制，选择适合的化学反应模型具有重要意义。
此外，边界条件的设定也是湍流燃烧数值模拟的一个关键问题。
边界条件的合理设定可以保证计算结果的准确性。
常用的边界条件有Inflow Boundary Condition、Outflow Boundary Condition、Wall Boundary Condition等。
对于湍流燃烧数值模拟，还需要考虑湍流场的边界条件，例如由湍流脉动引起的湍流输运方程中的涡粘性项的边界条件等。
最后，计算方法的选择也对湍流燃烧数值模拟的结果和计算速度有着重要的影响。
常用的计算方法有有限差分法(FDM)、有限元法(FEM)和有限体积法(FVM)等。
这些方法在计算精度和计算速度方面各有优势，需要根据具体问题的要求选择适当的方法。
中尺度大气数值模拟及其进展
中尺度大气数值模拟及其进展中尺度大气数值模拟及其进展中尺度大气数值模拟是指对中尺度大气运动、湍流、边界层、云微物理、辐射传输等过程进行数值模拟的一种方法。
近年来，随着计算机技术的快速发展和观测技术的不断进步，中尺度大气数值模拟的研究已经取得了许多重要的进展，对于气象预报、气候变化研究和环境污染预测等方面都起到了重要的作用。
中尺度大气数值模拟的目标是通过计算空间和时间上的大量物理量，来模拟和预测中尺度大气运动过程。
中尺度大气运动是指介于大尺度天气系统和小尺度湍流系统之间的系统，其典型特征是空间尺度在几十公里到几百公里之间，时间尺度在几分钟到几小时之间。
中尺度大气运动包括了许多重要的现象，如大气锋面、对流云团、飑线等，对于气象预报和气候变化研究具有重要的意义。
中尺度大气数值模拟的基本原理是通过数值方法将大气方程离散化，并通过数值解算得到大气运动的演化过程。
其中，最常用的模型是基于Navier-Stokes方程的大气动力学模型，通过有限差分、谱方法等数值技术对方程进行求解。
此外，为了更好地模拟大气过程，中尺度大气数值模拟还必须考虑到湍流的影响，湍流参数化是其中的关键技术之一。
近年来，随着计算机技术的不断进步，中尺度大气数值模拟的能力也得到了极大的提高。
传统的数值模拟方法需要通过将整个大气划分成若干个网格，然后分别对每个网格进行计算，这种方法在计算量和存储空间上都有较大的挑战。
为了克服这些问题，新型的数值模拟方法应运而生，如有限元方法、有限体积方法和伪谱法等。
这些方法可以更好地处理复杂的地形、不均匀的边界条件和非线性问题，提高了数值模拟的计算效率和精度。
除了数值方法的发展，观测技术的进步也为中尺度大气数值模拟提供了更多的观测数据，从而提高了数值模拟的准确性和可靠性。
现代大气观测技术，如雷达、卫星和飞机观测等，可以提供高时空分辨率的大气观测数据，在验证和改进数值模拟模型方面发挥重要作用。
此外，数据同化技术的应用也为中尺度大气数值模拟提供了新的思路和方法，通过将观测数据与数值模拟结果进行融合，可以进一步提高数值模拟的准确性和预报能力。
飞机设计优化中流场数值模拟方法的研究及应用创新
飞机设计优化中流场数值模拟方法的研究及应用创新引言：飞机设计优化是现代航空工程中的重要研究领域之一。
在飞机设计阶段，通过模拟流场数值，可以提供对飞机的空气动力学性能进行准确评估的有效工具。
本文将对流场数值模拟方法在飞机设计优化中的研究与应用进行深入探讨，旨在探索创新的方法以提高飞机设计效率和性能。
一、流场数值模拟方法的概述流场数值模拟是一种基于计算流体力学（CFD）的技术，通过离散方程组的求解，得到模拟自由空气中的速度、压力、温度等物理量的数值解。
流场数值模拟方法的基本原理是通过数值计算来模拟真实流体运动的物理现象。
二、流场数值模拟方法在飞机设计优化中的应用现状1. 飞行器气动性能预测流场数值模拟方法可用于预测飞行器在不同飞行状态下的气动性能。
通过改变飞行器的几何形状和工况参数，可以预测其升力、阻力、升阻比等性能指标，为飞机设计提供重要的依据。
2. 空气动力学优化设计在飞机设计的过程中，通过优化飞机的气动外形，可以减少阻力、提高升力、改善飞行稳定性和操纵性。
流场数值模拟方法可以高精度地评估不同设计方案的气动性能，为优化设计提供指导。
3. 结构强度分析除了考虑飞机的气动性能，流场数值模拟方法还可以用于分析飞机在飞行和地面操作时所受到的各种载荷，如空气动力载荷、惯性载荷、操纵系统载荷等。
这对于飞机的结构强度和寿命评估非常重要。
三、流场数值模拟方法的研究进展1. 网格生成技术的改进网格生成是流场数值模拟的基础，良好的网格质量对数值模拟结果的准确性和稳定性至关重要。
近年来，研究人员通过改进传统网格生成算法和开发自适应网格生成技术，提高了数值模拟的效率和准确性。
2. 数值模拟算法的发展为了提高数值模拟的计算效率和准确性，研究人员不断改进传统的数值模拟算法，并提出了一些创新的算法。
例如，基于稳定性的数值模拟方法、并行计算技术等，可以有效地缩短数值模拟的计算时间，同时减小数值模拟误差。
3. 模型与物理效应的改进为了更准确地模拟飞机的流场现象，研究人员通过改进数学模型和物理模型，考虑了更多的气动效应，如湍流、化学反应、燃烧等。
数值模拟技术的最新进展
数值模拟技术的最新进展近年来，随着计算机技术的飞速发展，数值模拟技术也越来越得到人们的关注。
数值模拟技术是指利用数值方法对物理过程进行仿真和计算的技术，它在许多领域中都有应用，如机械工程、航空航天、地质勘探等等。
本文主要介绍数值模拟技术在目前的最新进展。
1. 基于深度学习的数值模拟深度学习技术是近年来飞速发展的人工智能技术之一，它通过模仿人类神经网络的结构和学习方式，从大量数据中学习并识别模式。
在数值模拟方面，深度学习技术可以用于建立高效、准确的模型。
例如，科学家可以将深度学习技术用于流体力学模拟中，这使得模拟能够更快速地进行，并且可以更加准确地预测流体力学现象，例如湍流、流体分离等等。
2. 并行计算技术随着计算机硬件技术的迅速发展，现在的计算机系统已经能够支持大规模并行计算。
并行计算技术是指将计算任务分成多个子任务来同时执行，以提高计算速度。
在数值模拟中，大规模并行计算技术可以极大地提高计算效率，从而使得更加复杂的模拟成为可能。
例如，人们可以用并行计算技术来模拟地震波传播过程，这将有助于更好地理解地震现象的本质。
3. 多物理场耦合模拟多物理场耦合模拟是指通过数值算法将不同物理学领域（如流体力学、电磁学、结构力学等）中的方程同时求解，以模拟多物理场耦合的物理现象。
例如，在航空航天工程中，飞机的设计要求同时考虑结构力学、燃烧、流体力学等多个因素，这时就需要用到多物理场耦合模拟技术。
目前，多物理场耦合模拟技术已经成为数值模拟领域中的一个重要方向。
4. 高性能计算技术高性能计算技术是指计算机系统通过优化计算资源的配置，以提高计算效率和性能。
在数值模拟中，高性能计算技术尤其重要，因为模拟过程中需要进行大量的计算，需要在有限的时间内完成计算任务。
例如，在气象预报领域中，高性能计算技术可以帮助气象预报模型更加准确地预测天气情况，提高天气预报的准确度。
5. 高精度数值算法高精度数值算法是指通过提高数值计算方法的精度，以提高模拟效果的技术。
THMC多场耦合作用下岩石力学实验与数值模拟研究进展
THMC多场耦合作用下岩石力学实验与数值模拟研究进展随着科技的不断发展，THMC（热-水-力-化学）多场耦合作用在岩石力学领域中的应用日益广泛。
THMC多场耦合作用是指热、水、力、化学等多种因素相互作用影响岩石力学性质的现象。
在岩石力学实验与数值模拟研究中，THMC多场耦合作用下的岩石力学性质成为研究的热点之一、本文将介绍THMC多场耦合作用下岩石力学实验与数值模拟研究的最新进展。
一、THMC多场耦合作用下岩石力学实验研究进展1.热-水-力-化学耦合实验平台的建立：近年来，越来越多的研究者开始搭建热-水-力-化学（THMC）耦合实验平台，用于研究岩石在多场耦合作用下的力学性质。
这些实验平台不仅可以控制温度、水分、应力等多个因素，还可以监测岩石的物理化学变化，为研究THMC多场耦合作用下的岩石力学性质提供了良好的实验条件。
2.THMC多场耦合作用下岩石强度实验研究：研究者通过实验方法，探讨了THMC多场耦合作用对岩石强度的影响。
实验结果表明，热-水-力-化学多场耦合作用可以显著影响岩石的强度和破坏模式，对岩石的稳定性和安全性产生重要影响。
3.THMC多场耦合作用下岩石渗透性实验研究：研究者还通过实验方法，研究了THMC多场耦合作用对岩石渗透性的影响。
实验结果表明，热-水-力-化学多场耦合作用会导致岩石渗透性的变化，进而影响地下水流动和岩石的稳定性。
1.基于离散元法的THMC多场耦合作用下岩石模拟：离散元法是一种用于模拟岩石颗粒间相互作用的方法，可以很好地模拟THMC多场耦合作用下岩石的行为。
研究者利用离散元法开展了大量的数值模拟研究，揭示了THMC多场耦合作用对岩石结构、力学性质和破坏模式的影响。
2.基于有限元法的THMC多场耦合作用下岩石模拟：有限元法是一种广泛应用于岩石力学领域的数值模拟方法，能够准确地描述THMC多场耦合作用下岩石的力学行为。
研究者对THMC多场耦合作用下的岩石进行了有限元分析，从而揭示了多场耦合作用对岩石应力、变形和破坏的影响规律。
湍流的数值模拟方法进展
3 大涡模拟（LES ）湍流大涡数值模拟（LES ）是有别于直接数值模拟和雷诺平均模式的一种数值模拟手段.利用次网格尺度模型模拟小尺度湍流运动对大尺度湍流运动的影响即直接数值模拟大尺度湍流运动, 将N —S 方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。
3。
1 基本思想很多尺度不同的旋涡一起组成了湍流运动平均流动主要取决于大漩涡的流动，大尺度运动则受到小旋涡的影响。
流动中的大涡实现了动量、能量质量、热量的交换，耗散主要是由于小涡作用的。
大旋涡中受到流场形状、阻碍物的影响，,使大漩涡的各向异性更加明显。
然而小漩涡之间各项同性，相互没有太大的区别，所以建立统一的模型比较容易一些.综上所述,大涡模拟将湍流瞬时运动量通过滤波将运动分成小尺度和大尺度.大尺度的运动受到小尺度的运动的影响可以通过应力项（类似于雷诺应力项）来表示，即为亚格子雷诺应力，以建立这种模型的方法来模拟。
而大尺度则是求解运动微分方程而计算出来的，也就是说大涡模拟,要先过滤掉小尺度的脉动,然后再推出小尺度的运动封闭方程以及大尺度的运动控制方程。
3。
2 滤波函数正如上面提到，大涡模拟要先将流动变量分解成小尺度量和大尺度量，我们把这个作用叫做滤波.滤波运算就是在一区域内按照一定的条件对函数进行加权平均，作用是将高波数滤掉，使低波数保留，滤波函数的特征尺度决定了截断波数的最大波长,下面三种滤波函数是最为常用的主要有以下三种：盒式、富氏截断以及高斯滤波函数.不可压常粘性系数的湍流运动控制方程为N-S 方程:j ij i j j i i x S x P x u u t u ∂⋅∂+∂∂-=∂∂+∂∂)2(1γρ式中:S 拉伸率张量，表达式为：2/)//(i j j i ij x u x u S ∂∂+∂∂=；γ分子粘性系数；ρ流体密度。
设将变量i u 分解为方程（11）中i u 和次网格变量(模化变量）'i u ,即'+=i i i u u u ，i u 可以采用Leonard 提出的算式表示为:(11）式中)(x x G '-称为过滤函数，显然G(x)满足x d x u x x G x u i i '''-=⎰+∞∞-)()()(⎰+∞∞-=1)(dx x G3.3 控制方程将过滤函数作用与N —S 方程的各项，得到过滤后的湍流控制方程组:由于无法同时求解出变量i u 和j i u u ,所以将j i u u 分解成i j i j ij u u u u τ=⋅+，ij τ即称为次网格剪切应力张量(亦称为亚格子应力)。
地下水数值模拟的研究与应用进展
地下水数值模拟的研究与应用进展【摘要】地下水数值模拟是地下水研究领域的重要工具，随着数值模拟方法的不断发展，其在水资源管理、环境保护和地质勘探等领域的应用也越来越广泛。
本文从数值模拟方法的发展、在水资源管理中的应用、在环境保护中的应用、在地质勘探中的应用以及未来发展方向等方面进行了系统的总结和探讨。
研究表明，地下水数值模拟在提高水资源利用效率、保护地下水资源、指导环境管理和勘探地下资源等方面具有重要意义。
加强地下水数值模拟的研究和应用，将对促进资源有效利用和环境保护具有积极的推动作用。
未来，我们需要进一步完善数值模拟方法，提高模拟精度，探索更广泛的应用领域，推动地下水数值模拟在各领域的发展和应用。
【关键词】地下水数值模拟、研究、应用、发展、水资源管理、环境保护、地质勘探、未来发展方向、重要性、总结、展望1. 引言1.1 地下水数值模拟的研究与应用进展地下水数值模拟是通过数学模型和计算机仿真技术，对地下水系统的水文地质特征进行描述和分析的一种方法。
随着计算机技术的不断进步和地下水问题的日益凸显，地下水数值模拟在水资源管理、环境保护、地质勘探等领域中扮演着重要角色。
在过去的几十年中，地下水数值模拟方法得到了长足发展。
从最初的一维流动模型，到如今的三维多孔介质模型，模拟精度和可靠性不断提高。
各种数值模拟软件的涌现，也为地下水研究提供了便利。
地下水数值模拟在水资源管理中的应用主要包括水资源评价、水资源保护、水资源规划等方面。
通过模拟地下水流动、水质变化等过程，可以更好地指导水资源管理工作，保障人民的饮用水安全。
在环境保护领域，地下水数值模拟被广泛应用于地下水污染源追踪、地下水保护区划定等方面。
通过模拟地下水流动和污染传输，可以及早发现、预防和处理地下水污染事件，减轻环境压力。
地下水数值模拟还在地质勘探领域发挥重要作用。
通过模拟地下水对地下结构的影响，可以为石油、矿产勘探提供重要参考依据。
未来，地下水数值模拟方法将继续发展，模拟精度将进一步提高。
中尺度大气数值模拟及其进展
中尺度大气数值模拟及其进展中尺度大气数值模拟及其进展一、引言大气数值模拟是一种使用数学方程和计算机算法来模拟大气运动和气象现象的方法，它不仅能够帮助预测和研究天气、气候变化等现象，还可为决策提供重要参考。
在气象学研究领域，中尺度大气数值模拟被广泛应用，具有重要的意义。
本文将介绍中尺度大气数值模拟的基础理论和方法，并探讨其在气象学领域中的进展。
二、中尺度大气数值模拟的基础理论和方法中尺度指大气运动的空间尺度在几十到几百公里之间，时间尺度在几小时到几天之间。
中尺度大气数值模拟的基础理论是对大气运动和物理过程的基本方程进行数学化处理，建立相应的模型。
其中，最常用的模型是基于质量守恒、动量守恒、热量守恒和状态方程的Navier-Stokes方程。
为了简化计算，通常还采用了一些物理参数化方案，如湍流参数化、云微物理参数化等。
中尺度大气数值模拟的方法可以分为欧拉法和拉格朗日法。
欧拉法是在空间网格上离散化基本方程，通过数值迭代求解得到大气场的时空分布。
拉格朗日法则是跟踪气体的运动轨迹，通过将大气分成许多气团来模拟大气运动。
三、中尺度大气数值模拟在气象学领域的应用中尺度大气数值模拟在气象学领域有着广泛的应用。
首先，它可以用于天气预报，通过模拟大气运动，结合实时观测数据，可以提供准确的天气预报结果。
其次，中尺度大气数值模拟还可以用于研究气象灾害，如暴雨、台风等的形成机制和前后过程，从而为灾害预防和减轻提供科学依据。
此外，中尺度大气数值模拟还可以用于研究气候变化，如模拟气候系统中的能量和水分交换，探索气候变化的内在机制。
四、中尺度大气数值模拟的进展随着计算机技术的不断发展和模型改进，中尺度大气数值模拟在气象学领域取得了许多重要的进展。
首先，模拟精度显著提高，模型对大气物理过程的描述更加准确。
其次，模拟时间和空间分辨率不断增加，模拟结果更加细致。
此外，数据同化技术的应用使得模拟结果与实况数据更加吻合，提高了模式的可信度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
三、精细建模及新型模拟技术
◆条件模拟技术 ◆流线法模拟技术 ◆多模型统一方程技术
（二）并行求解技术
并行油藏数值模拟实现的并行环境 硬件环境、并行编程环境 ◆共享内存式并行计算机：SGI Power Challenge， Origin 2000等； ◆分布式内存并行计算机：IBM SP2, 曙光2000 ◆局部共享全局分布内存式并行计算机：IBM－SP3， 曙光3000等； ◆工作站网络集群； ◆微机网络集群。
一、微分方程离散化及网格 （三）非结构网格技术
PEBI网格 （Perpendicular Bisection Grids）是指 网格块内的任一点距其 自身网格节点的距离比 到任一其他网格节点的 距离都要近的空间区域， 它是一种局部正交网格。
PEBI网格
二、偏微分方程组的求解技术
（一）预处理共轭梯度法求解技术
二、偏微分方程组的求解技术
（一）预处理共轭梯度法求解技术
3) 预处理的正交极小化方法 研究表明，正交极小化方法的收敛速度与矩阵的A 性质关系很大。 4） 用不完全分解作预处理 最常用且较有效的预处理方法是不完全分解。对 任意可逆矩阵都存在惟一的分解。
二、偏微分方程组的求解技术
（一）预处理共轭梯度法求解技术
油藏数值模拟所面对的是一个大型、稀疏、非对 称的矩阵方程的求解问题 目前，解上述矩阵最有效的方法就是预处理共轭 梯度类解法（Preconditional Gradient Methods）。
二、偏微分方程组的求解技术
（一）预处理共轭梯度法求解技术
1) 共轭梯度法（CG方法） 共轭梯度法是一种迭代法，最多迭代步（是系数 矩阵的阶数）就会收敛到真解, 它又是一种直接解法。 2) 正交极小化方法（OTHOMIN） 由于油藏数值模拟的矩阵是非对称的，有时甚至 非正定的，故对共轭梯度法加以改进，以适应油藏数 值模拟的需要。
一、微分方程离散化及网格 （一）网格加密技术
常规网格加密
静态局部网格加密
一、微分方程离散化及网格 （一）网格加密技术
局部网格加密经历了由早期的静态局部网格 加密（Static Grid Refinement）,到动态局部网 格加密（Dynamic Grid Refinement），由杂交网 格（Hybrid Grid）到多水平组合网格 （Multilevel Composite Grid ），由解决串行 问题到解决并行问题这样一个发展历程。
一、微分方程离散化及网格 （一）网格加密技术
组合网格系统
一、微分方程离散化及网格 （一）杂交网格技术
在近井附近采用正交曲线网格，和矩形网格结合 起来处理水的推进问题，形成了杂交网格技术。
粗网格－杂交网格
细网格－杂交网格
一、微分方程离散化及网格 （三）非结构网格技术
90年代随着三维地质模型的越来越精细化，要求 数值模拟向精细化方向发展，要求数学模型采用张量 形式的渗透率等。非结构网格（包括PEBI和CVFE）就 是针对上述网格技术中存在的实际问题而提出的一种 新的网格技术。
一、微分方程离散化及网格 （一）网格加密技术 静态的局部网格加密的网格系统在整个计算过
程中是固定不变的，常用于断层、尖灭、裂缝及井附 近的加密。
动态的局部网格加密（Heinemann,1983）（如
图9.6.3所示）主要用于追踪水驱前沿。
一、微分方程离散化及网格 （一）网格加密技术 动态的局部网格加密（Heinemann,1983）（如
常规串行油藏模拟的 求解流程
时 间 步 循 环
牛 顿 迭 代
不管模拟的区域有多大， 模拟历史时间有多长，计算 过程都是一个时间步一个时 间步地进行
二、偏微分方程组的求解技术
（二）并行求解技术
3 3 1 4 2 1 2 4
平行算法示意图
采用区域分解的方法将油藏划分为若干子区域 然后利用并行算法并行求解各子区域问题，最终 达到完成整个油藏的数值模拟的计算。
油藏数值模拟的技术进展
●数值模拟软件技术 ●数值模拟应用技术
油藏数值模拟的技术进展
一、微分方程离散化及网格 （一）网格加密技术 早期的局部网格加密方法,粗细网格是一 起编号，统一求解，导致了稀疏矩阵结构不 规则。 由于引入了加密区域，使得系数矩阵的 带宽增加，计算量大大地增加了
一、微分方程离散化及网格 （一）网格加密技术 局部网格加密技术是解决精细描述大型、 复杂的油藏数值模拟问题的有效方法 。 从80年代初引入到油藏数值模拟后，油 藏数值模拟工作者一直在做着不懈的努力。
图9.6.3所示）主要用于追踪水驱前沿。
动态局部网格加密
一、微分方程离散化及网格 （一）网格加密技术
多水平组合网格方法（Multilevel Composite 多水平组合 Grid Scheme 简称ML）以及快速自适应组合网格方 法（Fast Adaptive Composite Grid简称FAC） （McComick ,Thomas ,1986）对局部区域进行逐级 加密，可以避免大网格与小网格相邻,而且FAC非常 适用于开窗技术。所谓组合网格就是由局部粗网格 和各级局部加密的细网格组成，而全局粗网格又叫 基础网格。
（二）并行求解技术
并行油藏数值模拟实现的并行环境 硬件环境、并行编程环境
研制开发微机并行油藏模拟软件是今后发展 的一个主要研究方向。
（三）网格粗化技术
网格粗化技术的关键主要是渗透率的网格粗化。 主要有以下几种方法： 1)按流度加权进行网格粗化。 2)按静态模拟的方法使得粗网格的渗透率产生的流 量与细网格产生的流量基本相当。 3)进行动态的预模拟处理
预处理——嵌套分解（nested）方法 残量约束法（gemres方法）等
预处理正交极小化类的方法在实际应用中是相当 有效的，有时能比一般的迭代法快到上百倍。
二、偏微分方程组的求解技术
（二）并行求解技术
输入静态数据 初始化计算 输入动态数据 岩石性质计算 PVT性质计算 井产量处理 系数计算 线性系统求解 解校正 结束