数值模拟方法
数值模拟的概念与方法
他对中科院计算数学所的研究生们说：“从采矿、水库大坝到地下隧道 工程等，世界各国的工程师面临太多的危险。在这些方面，数学是非常有用 的，我们周围的人都需要数学。我希望下一代的数学家们，特别是你们，站 在计算数学与工程之间，最重要的是用发明出的一些数学方法和工具，写出 很好的教科书，把数学交给工程师，追上这个时代。”
且计算精度、计算效率高， 更适用于均质材料和线性性 态情况。
➢软件：Examine2D、Examine3D
C 离散单元法（DEM）
岩体往往为众多的节理或结构面所切割，在某些情况下， 岩体不能视为连续介质，具有明显的不连续性，很难用连 续介质力学方法如有限单元法来处理。
离散单元法是处理非连续介质力学的数值方法，特别适用 于节理岩体的应力分析，在土木工程方面应用广泛，尤其 在边坡稳定分析方面。
.点 (质量)
. . 线(弹簧，梁，杆，间隙)
面 (薄壳, 二维实体,
.. .. 轴对称实体)
.. . .... 体(三维实体)
节点自由度是随连接该节点 单元类型 变化的。
J
三维杆单元 (铰接)
UX, UY, UZ
I
I
L
K
二维或轴对称实体单元
L
UX, UY
I
I
J
P M
L
I
O 三维实体结构单元
N
UX, UY, UZ
载荷 约束
节点 单元
UY ROTY
ROTZ UZ
UX ROTX
结构 DOFs
数值模拟方法在物理仿真中的应用研究
数值模拟方法在物理仿真中的应用研究数值模拟方法是一种通过计算机模拟现实世界物理过程的方法，已经广泛应用于物理仿真领域。
本文将介绍数值模拟方法在物理仿真中的应用研究，并分析其优势和挑战。
一、数值模拟方法概述数值模拟方法是将物理过程转化为计算机可处理的数值问题，通过数值计算的方法对其进行模拟和仿真的过程。
数值模拟方法的核心是将物理过程转化为数学方程，并通过计算机求解这些方程来模拟和预测物理过程的行为和性质。
二、数值模拟方法在物理仿真中的应用1. 流体力学仿真流体力学仿真是数值模拟方法在物理仿真中的一个重要应用领域。
通过数值模拟，可以研究各种流体的运动和行为，了解流体的流动规律和性质。
例如，可以模拟空气动力学，研究飞机在不同飞行状态下的气动性能，优化飞机的设计；还可以模拟液体的流动，研究水流在河道中的流速分布，指导河道治理工程等。
流体力学仿真的过程通常包括建立流体模型、选择适当的数值方法和求解算法、设置边界条件和初值条件，进行数值计算和结果分析。
2. 结构力学仿真结构力学仿真是另一个数值模拟方法在物理仿真中的重要应用领域。
通过数值模拟，可以研究各种物体的力学行为和性能。
例如，可以模拟建筑物的抗震性能，分析结构在地震作用下的应力和变形情况，指导建筑物设计和抗震加固；还可以模拟材料的力学性能，研究材料的应力分布和断裂行为，指导材料的选择和加工工艺等。
结构力学仿真的过程通常包括建立物体模型、选择适当的模拟方法和求解算法、设置边界条件和初值条件，进行数值计算和结果分析。
3. 电磁场仿真电磁场仿真是数值模拟方法在物理仿真中的另一个重要应用领域。
通过数值模拟，可以研究各种电磁现象和电磁器件的性能。
例如，可以模拟电磁波的传播和辐射特性，研究无线通信和天线设计；还可以模拟电磁场的分布和场强分布，研究电磁屏蔽和电磁兼容等问题。
电磁场仿真的过程通常包括建立电磁模型、选择适当的数值方法和求解算法、设置边界条件和初值条件，进行数值计算和结果分析。
数值模拟方法
精品文档数值模拟方法科学研究与解决工程问题的基础在于物理实验与实物观测，例如对金属材料的凝固过程进行物理实验、对天体运行进行观测。
现代科学研究方法的核心是通过实验或观测建立研究对象的数学模型，基于数学模型进行研究与分析。
这种研究方法可以追溯到伽利略的工作，成熟于牛顿的三大定律与微积分。
采用实物模型进行物理实验的研究周期长、投入大，有时甚至无法在实物上进行，如天体物理的研究。
在数学模型上进行的数值模拟研究具有研究周期短、安全、投入少，已经成为不可或缺的工具。
数值模拟方法的应用对象分为三个层次1) 宏观层次：常见的工程建筑、制造设备、零件等;2) 界观层次：材料的微观组织与性能，如金属材料的晶粒度影响其屈服强度;3) 微观层次：基本物理现象与机理，如金属材料凝固时的结晶与晶粒生长过程。
宏观与界观层次的数值模拟方法包括：1) 有限差分方法(Finite Difference Method, FDM) ：微分方程的直接离散方法;2) 有限元单法(Finite Element Method, FEM)用有限尺度的单元的集合来代替连续体，分为Lagrange方法，Euler方法，ALE方法；3) 边界单元方法(Boundary Element Method, BEM) ：一种半解析方法;4) 有限体积方法(Finite Volume Method, FVM) ：把空间划分成有限尺度的体积单元，连续体通过这些在空间上固定的体积单元，单元的空间位置不变；5) 无网格方法(Meshless Method) ：只布置结点，不需要划分单元网格，有权函数。
微观层次的数值模拟方法包括：1) 第一原理法(First Principle Simulation) ：量子力学方法，直接计算原子的电子结构；2) 元胞自动机方法(Cellular Automata) ：把空间用元胞演化、元胞的局部相互作用来描述复杂的、全局的系统。
3) 蒙特卡洛方法(Monte Carlo Method ) ：把颗粒运动定义为随机过程，用势能的变化来判断颗粒运动能否被接受。
随机微分方程的数值模拟方法
随机微分方程的数值模拟方法随机微分方程（Stochastic Differential Equations，简称SDEs）是描述包含随机项的微分方程。
它们在金融学、物理学和生物学等领域中广泛应用，尤其在随机模型建立和数值模拟方面有着重要的作用。
为了模拟和解决随机微分方程，研究者们开发了各种数值模拟方法。
这些方法的目标是通过离散化时间和空间来近似SDE的解，以获得数值解。
在本文中，我将介绍几种常用的数值模拟方法，包括欧拉方法、米尔斯坦方法和龙格-库塔方法。
我们将从简单的欧拉方法开始，逐渐深入探讨这些方法的优点和局限性。
1. 欧拉方法（Euler Method）欧拉方法是最简单和最直接的数值模拟方法之一。
它将区间分成若干小的子区间，然后使用差分逼近来计算每个子区间内的解。
欧拉方法的基本思想是将微分方程中的导数用差分代替，从而将微分方程转化为差分方程。
欧拉方法的数值格式如下：然而，欧拉方法的缺点在于其精度较低，特别是当时间步长较大时。
它也不能很好地处理某些随机微分方程的特殊情况。
2. 米尔斯坦方法（Milstein Method）米尔斯坦方法是对欧拉方法的改进，目的是提高精度。
它通过在欧拉方法的基础上添加额外的项来纠正误差，从而提高数值解的准确性。
米尔斯坦方法的数值格式如下：相比于欧拉方法，米尔斯坦方法在同样的时间步长下通常能够提供更准确的数值解。
然而，对于某些特殊的随机微分方程，米尔斯坦方法也可能存在一些问题。
3. 龙格-库塔方法（Runge-Kutta Method）龙格-库塔方法是一类更为复杂但精度更高的数值模拟方法。
它基于对SDE进行多次逼近来得到数值解，通常可以达到较高的准确性。
龙格-库塔方法的基本思想与常规微分方程的龙格-库塔方法类似，但在计算过程中需要额外考虑随机项的贡献。
相比于欧拉方法和米尔斯坦方法，龙格-库塔方法的数值格式更为复杂，但其准确性和稳定性更高。
总结和回顾：通过本文的介绍，我们对随机微分方程的数值模拟方法有了初步的了解。
数学中的数值模拟方法
数学中的数值模拟方法数学作为一门科学，其应用范围越来越广泛，而数值模拟方法是数学在实际应用中的一个重要组成部分。
数值模拟是利用数学计算机方法，通过对数学模型的数值求解，得到与实际物理过程相对应的数值结果的一种方法。
本文将介绍数学中的数值模拟方法。
一、数值模拟方法的应用数值模拟方法在物理学、化学、生物学、工程学、地球科学等领域均有广泛的应用。
例如，在工程学中，数值模拟可用于模拟过程中的流体力学、热传递、材料力学等。
在物理学中，可用于模拟天体力学中的行星运动、物理量的计算等。
二、有限差分法有限差分法是计算微分方程的一种数值方法。
通过将微分方程中的函数在有限个点上展开，将微分项用差分近似表示，从而将微分方程变为代数方程组。
这种方法可用于求解一维、二维或三维的偏微分方程。
在计算中，有限差分法一般采用迭代方法进行求解。
三、有限元法有限元法主要应用于计算结构力学和固体力学中的问题。
这种方法将结构分解为有限数量的小单元，每个小单元内部的材料和力学特性相同时，对每个小单元进行力学计算，通过将小单元的结果组合成大体系的结果，得到整个结构的受力状态或变形。
四、蒙特卡罗方法蒙特卡罗方法是一种随机数学模拟方法。
它应用随机数的方式解决了一些复杂问题，包括点和粒子运动、概率模型、射线传输等。
利用蒙特卡罗方法，可以在减少计算机运算量的同时，还能得到很好的模拟效果。
五、数值优化方法数值优化方法是一种用于解决优化问题的计算机模拟方法。
在优化问题中，通常需要确定目标函数在一组给定条件下的最大值或最小值。
数值优化方法可以通过迭代计算过程，逐渐接近最优解。
常用的数值优化方法包括模拟退火、遗传算法和粒子群优化等。
六、求解微分方程的方法微分方程是物理学和工程学中常见的数学方法。
可以通过数值模拟方法来求解微分方程。
其中较为常用的有：欧拉法、改进欧拉法、龙格-库塔法、泰勒展开法等。
七、总结数值模拟方法在科学领域中有着广泛的应用。
不同的数值模拟方法适用于不同的问题。
数值模拟方法
数值模拟方法数值模拟方法是一种通过计算机对物理、化学、工程等领域中的现象进行模拟和分析的方法。
它通过建立数学模型，利用数值计算方法对模型进行求解，从而得到所研究系统的一些重要信息。
数值模拟方法已经成为科学研究和工程技术领域中不可或缺的工具之一。
在科学研究中，数值模拟方法可以帮助研究人员更好地理解复杂的物理现象。
例如，在天文学中，科学家们可以利用数值模拟方法来模拟宇宙中恒星的形成和演化过程；在地球科学领域，数值模拟方法可以用来模拟地震波的传播规律。
而在工程技术领域，数值模拟方法则可以帮助工程师们设计更安全、更高效的产品和工艺。
数值模拟方法的核心是建立数学模型。
数学模型是对真实系统的抽象和简化，它可以是基于物理定律的微分方程模型，也可以是基于统计规律的随机模型。
建立好数学模型之后，就需要选择合适的数值计算方法对模型进行求解。
常用的数值计算方法包括有限差分法、有限元法、谱方法等。
这些方法各有特点，适用于不同类型的问题。
在进行数值模拟时，我们需要关注模拟结果的准确性和可靠性。
准确性是指模拟结果与真实系统的符合程度，而可靠性则是指模拟结果的稳定性和可信度。
为了提高模拟结果的准确性和可靠性，我们需要不断改进数学模型和数值计算方法，同时也需要考虑计算机的计算精度和稳定性。
除了关注模拟结果的准确性和可靠性，我们还需要关注模拟的效率。
随着计算机计算能力的不断提高，我们可以利用并行计算、高性能计算等技术来加速数值模拟的过程。
这样可以大大缩短模拟的时间，提高工作效率。
总的来说，数值模拟方法是一种强大的工具，它在科学研究和工程技术中发挥着重要作用。
通过建立数学模型和选择合适的数值计算方法，我们可以更好地理解复杂的现象，设计创新的产品，解决实际的工程问题。
随着计算机技术的不断发展，数值模拟方法也将不断完善，为人类的发展进步提供强大的支持。
数值模拟方法
数值模拟方法数值模拟方法是一种通过计算机模拟数学模型来解决实际问题的方法。
它是利用数值计算方法对不同领域的问题进行模拟和分析，是现代科学技术中的重要工具之一。
数值模拟方法在工程、物理、化学、生物等领域都有广泛的应用，可以帮助人们更好地理解和解决复杂的实际问题。
数值模拟方法的基本思想是将实际问题转化为数学模型，然后利用计算机进行数值计算，得到问题的近似解。
在进行数值模拟时，需要考虑到模型的准确性、计算的稳定性和计算的效率。
因此，数值模拟方法需要结合数学、计算机科学和实际问题的专业知识，进行综合分析和研究。
数值模拟方法的核心是数值计算方法，包括差分法、有限元法、谱方法等。
这些方法都是通过离散化连续问题，将其转化为离散的数学问题，然后利用计算机进行数值计算。
在实际应用中，需要根据具体问题的特点选择合适的数值计算方法，并对计算结果进行合理的分析和解释。
数值模拟方法在工程领域有着广泛的应用。
例如，在航空航天领域，数值模拟方法可以用来模拟飞机的气动性能，优化飞机的设计；在汽车工程领域，可以用来模拟汽车的碰撞安全性能，提高汽车的安全性能；在建筑工程领域，可以用来模拟建筑结构的受力情况，提高建筑结构的稳定性。
通过数值模拟方法，工程师可以更好地理解和分析复杂的工程问题，提高工程设计的效率和质量。
在物理学和化学领域，数值模拟方法也有着重要的应用。
例如，可以利用数值模拟方法模拟材料的结构和性能，研究材料的力学性能、热学性能和电学性能；可以利用数值模拟方法模拟化学反应的动力学过程，研究化学反应的速率和产物分布。
通过数值模拟方法，科学家可以更好地理解和预测物质的性质和行为，为新材料和新药物的设计提供理论支持。
在生物学领域，数值模拟方法也有着重要的应用。
例如，可以利用数值模拟方法模拟生物体内的生物力学过程，研究生物体的运动和变形；可以利用数值模拟方法模拟生物体内的生物化学过程，研究生物体的代谢和信号传导。
通过数值模拟方法，生物学家可以更好地理解和研究生物体的结构和功能，为疾病的诊断和治疗提供理论支持。
力学中的数值模拟方法
力学中的数值模拟方法力学是自然科学中研究物体运动和相互作用的学科。
力学的研究对象包括刚体、弹性体、流体等物质，而这些物质的运动和相互作用往往是非常复杂的。
为了更深入地了解这些现象，研究者们常常采用数值模拟方法。
本文将介绍在力学中常用的数值模拟方法和其应用。
1. 有限元法有限元法是解决力学问题的一种常用数值方法。
它将复杂的物体划分成有限个小元素，在每个小元素上进行基本方程的数值求解。
这些小元素可以是输入自然或几何区域的任意形状和大小。
通过将整个物体分解为由许多这样的小元素组成的形式，有限元法可以轻松处理具有复杂边界和几何形状的问题。
有限元法的一个重要优点是可以模拟多种不同的问题，例如，静力学问题，热力学问题和流体力学问题。
在建筑和航空航天科学中有限元法广泛应用，设计和优化桥梁、飞机机翼和汽车车身。
2. 边界元法边界元法是另一种广泛用于力学课题研究的数值模拟方法。
与有限元法相比，它的计算成本和计算时间更低。
其基本思想是借助几何中的经典定理——格林公式，将原方程转换为涉及单独表面积分的一组方程。
这些方程的求解是通过构造矩阵并进行数值求解得到的。
边界元法在流体动力学中的应用非常广泛，例如模拟液体流动和超声波传播等。
3. 分子动力学模拟分子动力学模拟是一种基于牛顿力学构建计算统计物理学的方法。
它通过建模粒子之间的相互作用来模拟分子系统的力学行为。
由于该方法可以与巨分子水平的化学反应联系起来，这使得它可以在化学和材料科学中应用得非常广泛。
通过使用物理特征的数值模拟，研究者们可以了解更多基于分子层面的成分内部运作和物理过程。
4. 自适应Mesh网格算法有些力学问题中变量可能有非常高的梯度，为解决这种问题，自适应Mesh算法应运而生。
自适应Mesh网格将整个求解域划分成相互交叉的奇下网格或三角形网格。
然后，当解的精度要求在较高的局部变化时，通过极小化给定误差级别来改变不同的小视窗大小，以便能够应对快速变化的解。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
式中bi，j为位置(i，j)处的水底高程，i，j为自由面的高度，kmax为垂直方向上网格单元的最 大数目。本手册中其它向量也运用类似的定义。
4.2动量的半隐式格式
速度场的基本半隐式化演进可被离散到n*的解空间，类似TRIM中的方法，即得：
式中G是一显式源项向量，1表示自由水面的“隐性度”（隐性度一般的取值范围为0.5 < 1 < 1.0已在ELCOM版本1中编码，但未经全面测试）。ＥＬＣＯＭ中默认的半隐格式是对从形 式上看为时间上一阶精确的自由水面演进采用后向欧拉离散法（即1＝1）。 经证明（Casulli和Cattani，1994），求解流体静力学方程组使用的后向欧拉法可被扩展为 方程4.2和4.3的一般二层格式，从形式来看为二阶精确（当1 =1时）。然而在粗网格模拟中， 数值离散精度的增加不一定导致模型技巧的提高。通常，在对许多湖泊和河口进行模拟时，正 压模态是通过CFL条件来求解的，CFL取值可能为5到10之间或大于10。在这种情况下，半隐式 离散化或许稳定，但水流物理过程表达“准确程度”的决定因素为现研究的流态类型。截断误 差的性质对于了解该方法性能至关重要。若采用一阶方法，导程误差为二阶，且在自由水面上 产生阻尼波。若采用二阶方法，导程误差扩散，且在自由水面产生在水域上传播的数值波；通 常使得线性正压波演化为陡峭前沿涌潮，使得受地形影响的表面波产生局地高速流。因此，一 阶方法可以较好描述自由水面形态和局部正压流速，但显示出自由水面惯性响应的多余阻尼。 相比之下，二阶方法以最低的数值耗散将能量保持为表面波的形式，但对于波的形态描述较 差。在静水力学求解中，弥散波将导致贯穿水柱的人工局部加压的产生，这对于求解技术是不 利的。相比之下，表面波的多余阻尼在风速减弱的时候将导致与正压响应相关的大规模运动的 减弱（即正压模式的“震荡”是衰减的）。总体上，高强度增压系统采用后向欧拉格式模拟较 好，因为两、三个周期以前的波能通常与一阶的物理过程不相关。 采用二层隐式离散化（Casulli和Cheng， 1992）或其它显式离散化技术，矩阵A可以表示
算子L()表示平流离散化，B()表示正压离散化，D()表示水平湍流扩散离散化。ELCOM的垂 向扩散的计算既可通过采用涡流粘度项，又可使用垂向混合模型。当采用涡流粘度项时，U = Un，而混合模型可用如下算子表示：
混合算子M()的具体描述见< .> 。
4.3水平扩散离散化
在文献（Casulli和Cheng，1992）提出的原TRIM方法中，当在迹线起点处对水平扩散项进 行离散化处理时，发现增加了其复杂性却并不能给现行方法的精度带来显著提高。因此，方程 4.8和4.9中的水平扩散项(Dx，Dy) 的离散化采用二阶模板，得：
.au/services/models/legacy/cwrsoft/doc/elcom_science_chine... 09/03/2010
Chapter 4
Page 3 of 10
4.4斜压离散化
在x方向上斜压项B被离散为：
式中k＝F表示包含自由面的单元。在y方向上也得到类似的扩散项和斜压项的表达方式。忽 略动量和标量迁移方程组中的垂向扩散项，目前的方案则可免去TRIM中对每个网格水柱的每个 水平速度份量和传输标量组成的三对角线矩阵求逆的步骤。
拉格朗日点处的变量的值可用三线迎风插值法获得：
上述方程便为半拉格朗日法构成了一个线性插值八节点模板。若流场近似于一维均匀流迎风方法。实际上，我们更容易把线性半拉格朗 日法看作在CFL高值时会不断地重置的三维线性迎风格式模板，因为线性半拉格朗日法表现出的 不可接受的数值扩散水平是线性迎风方法的特点。若拉格朗日点处的变量值采用二次插值方法 求得，则数值扩散可得到改善。
.au/services/models/legacy/cwrsoft/doc/elcom_science_chine... 09/03/2010
Chapter 4
Page 2 of 10
为：
中的可以由边界条件确定，a，b，c项分别为：
系数2通过对数值离散化技术的选择来确定。对于2=1的情况，垂向粘性项通过后向欧拉 法离散。对于2= 0的情况，离散化是显式的，且A除对角外其它的项均为0。ELCOM在用户选 择垂向湍流扩散模拟时取2=1，选择混合层模拟时，取2 = 0。 公式和中的源项G可被表示为：
式中，(i，j)下标表示被插值的垂线的位置， 和的取值为+/-{0，1，2}，正负号取决于模板的迎风方向 （U的位置下标增量k的正负号也同理得出）。每条线(i，j)的拉格朗日多项式系数通过如下标准拉格朗日系数 公式(例如：Al_Khafaji和 Tooley，1986）计算：
式中，zp为插值点的垂直坐标；正负号的选择必须满足模版的迎风性。垂向插值之后可得到 y方向上的3次水平插值，形式如下：
二次欧拉－拉格朗日方法
二次半拉格朗日法将使用三线插值（TRIM中已用到）的8点迎风格式模板扩展为使用二次拉 格朗日多项式插值的27点迎风格式模板。运用文献（Casulli 和 Cheng，1992）的符号表示 法，则算子对速度场（x平面）进行内插至位置(i+1/2-a，j-b，k-d)，其中a，b，d为实数，表示 以位置(i+1/2，j，k)为起点的位移。此位置为质点运动轨迹线起点的估计值，该质点穿过时间n 的速度场，经过时间t后移动到位置(i+1/2，j，k)。类似的符号表示法可运用于平面y＝j+1/2。 为方便起见，用上标（p）来标志轨迹线起点，则方程组4.2 和4.3中的平流算子L()（拉普拉斯 算子）可表示为：
4.6 水平对流离散化
在地球物理尺度上进行数值模拟时，难点之一是可能出现多种流态，尤其是内波可能在局部 范围内间竭地产生强烈的垂直运动。在内波分辨率要重点考虑的区域，数值模拟方法的选择是 由全流场中某个小区域的精确性和稳定性需求来决定的。在CFL < O(1)条件下，许多显式空间 离散方法是稳定的，但是当水流方向与网格边线不重合时，在CFL > O(0.5)条件下，这些方法 的三维计算精度通常很差。 二次方法的缺点是其计算量很大。然而在低CFL值（CFL<0.1）的区域，二次型半拉格朗日方 法的求解同样受控于使用二次迎风离散化方法所得的七点迎风格式模板的条件，这种相似性可 用于减少低CFL值区域的计算量要求，而又不会使得整个求解方法的精度有显著的降低。在不同 的区域使用不同方法的概念可被归纳为一种“混合”数值方法的概念。通常的混合模型是根据 某种衡量标准把流场划分为不同的水流区域，分别使用不同的求解方法。根据我们的目的，该 标准即是CFL值，在CFL值低时使用一种离散技术，而CFL值高时使用另一种离散技术。 当前的混合方法已测试到二层格式，在0 < CFL < 2的区域内使用二次半拉格朗日法进行离 散，该CFL值上限为可通过二次SL法计算得出并且在调整模板时无需额外计算量的最大CFL值。 在实际运用中，当CFL>2时，计算精度是不可靠的，因此该上限是一个非常合理的必要条件。在 CFL>2的区域，模型使用线性半拉格朗日离散法（即：TRIM中使用的方法）以使重构模板时的计 算量最小。
最后，x方向上的单次插值如下：
方程组4.22和4.23中的拉格朗日系数通过公式4.21计算，分别用x或y来替代z。拉格朗日差 值法采用二次模板是有利的，因为它降低了八节点线性模板中内波出现的人工阻尼，因而改进 了该方法求解分层流流域的自由运动的能力。这种改进在模拟分层流湖泊时是有必要的，而在 模拟水流物理过程以强制运动为主的河口时可能是次要的。如果不必重新定位模板就能计算出 在取值范围1 < CFL < 2内的流量源项，则二次模板计算的附加必要条件可稍加改进。
Page 5 of 10
维二次拉格朗日插值图解中也对该计算过程加以说明。对于不均匀网格的拉格朗日插值，考虑 采用与计算指数(i，j，k）相关的物理空间坐标(x， y， z）较为方便。如图4.2所示，在y方向 上的三条插值线中运用三维拉格朗日插值法计算得到任意点(x(p)(p) ， y(p)(p) ， z(p)(p) )上的 U(p)值，其中三条插值线是通过9次如下形式的垂向插值得到的：
文献（Casulli，1990）中讨论了有关双线性插值的二维轨迹线起点计算过程，图4.1 的二
.au/services/models/legacy/cwrsoft/doc/elcom_science_chine... 09/03/2010
Chapter 4
.au/services/models/legacy/cwrsoft/doc/elcom_science_chine... 09/03/2010
Chapter 4
Page 4 of 10
4.6.1 线性半拉格朗日方法
对于分层流，半拉格朗日方法效果很理想，因为在0<CFL<1的区域里，该法既精确又稳定。 此外，半拉格朗日法能在CFL值更高的区域里保持稳定，虽然当网格不能很好地处理水流流线 时，该方法在CFL值高的区域的精度会降低。然而，对于分层流，是否能在水平方向上的CFL>1 时提供稳定解对于数值方法的选择是次要的。最大时间步长的选取一般受限于以下两个条件之 一：（1）分层作用最强烈的区域中斜压波的传播速度；（2）在标量迁移过程中可接受的最大 数值扩散。前者要求波的CFL值小于1以维持稳定性，而后者的限制取决于模型结果所需要的网 格分辨率和时段长。 水平对流的半拉格朗日形式，是通过寻找连续空间中的近似点（拉格朗日点）来获得的，这 些拉格朗日点以△t为时间步长随速度场（U，V，W）平流输送至离散点(i，j，k)。用数值表 达，即质点的位置(i，j，k)沿着用速度场U，V，W表示的流线逆向移动。U，V，W流场可由单 时间层或多时间层获得，这取决于预期的精度和计算复杂程度。半拉格朗日离散方法的评述详 见文献（Staniforth和Cote ，1991）。在线性单时间层的半拉格朗日方法里，拉格朗日点通过 以下公式搜索：
.au/services/models/legacy/cwrsoft/doc/elcom_science_chine... 09/03/2010