简述几种常用数值方法的优势及适用性

本文涉及的产品
可观测链路 OpenTelemetry 版,每月50GB免费额度
函数计算FC,每月15万CU 3个月
容器服务 Serverless 版 ACK Serverless,952元额度 多规格
简介: 这篇文章简要介绍了几种常见的数值计算方法在解决偏微分方程中的应用。偏微分方程被用来描述自然界的物理现象,而当它们难以直接求解时,数值方法如有限元法(FEM)、有限差分法(FDM)、有限体积法(FVM)和边界元法(BEM)成为常用工具。有限元法适用于复杂问题,但计算量大;有限差分法简单但对复杂几何形状处理困难;有限体积法对网格质量要求低,适合流体力学;边界元法则仅需处理边界,适合边界主导的问题。各种方法各有优劣,常结合实际问题组合使用。例如,云道智造的伏图仿真平台集成多种数值方法,支持多物理场耦合仿真。

“天地盈虚自有时”,世间万物的变化都有规律可循。寻找万物规律并预测未来是人类一直孜孜以求的事情。

随着科学技术的发展,人们发现许多问题需要使用多个变量的函数来描述。十八世纪中叶,现代数学家们开始用偏微分方程(Partial differential equation)描述自然界物理场的变化规律。

然而,大多数偏微分方程难以有效求解。数值方法是应用最广泛的偏微分方程求解方法之一,包括有限元法、有限差分法、有限体积法、边界元法等。这些方法通过将偏微分方程离散化为有限维的计算问题,然后使用数值方法求解这些问题的解。

本文简要介绍几种常用的数值计算方法。

1. 有限元法(FEM, Finite Element Method)

有限元法是通过将连续空间区域离散为有限个小单元,将偏微分方程转化为每个单元内部的局部方程,然后将这些局部方程组合起来构成整个区域的方程组,通过求解方程组从而得到偏微分方程的近似解。关于有限元法的历史,可以参考《从无限到有限:有限元法的诞生》。

有限元法的优势在于具有高度的适应性、灵活性和计算精度,能够适应各种复杂的几何形状和边界条件(特别是对椭圆型问题有更好的适应性),可以通过调节有限元的数量和自由度的精度来改变求解的近似程度,还可以分析结构的局部细节问题。

但相比于其他方法,有限元法计算量大,需要消耗大量算力和时间;同时由于对于网格的质量敏感,不适合高度非结构化的问题,也难以处理可压缩流体等复杂问题。

有限元法是目前工业软件领域最常用的数值方法,可以广泛应用于结构力学、流体力学、热力学等领域。比如在结构分析中,可以用于分析结构的刚度、应力、振动等特性,常用于飞机、桥梁、建筑物等工程结构的设计和优化。

基于云道智造通用仿真PaaS平台伏图(Simdroid)开发的大型LNG储罐有限元分析软件(点击文字可了解更多),使用有限元法对LNG储罐进行结构和热仿真分析

2. 有限差分法(FDM, Finite Difference Method)

有限差分法的基本思想是把求解域划分为差分网格,用有限的网格节点来代替连续的求解域,并使用Taylor级数展开等方法,把定解问题中的微商换成差商,从而把原问题离散化为差分格式,进而求出数值解。这是一种将微分问题转化为代数问题的近似数值解法。

有限差分法是数值解法中最经典的方法,发展较早且较为成熟。相比于其他方法,有限差分法较为“简单粗暴”,直观易懂、通用性强,适用于简单几何形状和均匀网格的问题,但难以处理复杂几何形状和边界条件,且其精度取决于离散化程度。因此在工业软件领域,有限差分法的应用并不多见。

3. 有限体积法(FVM, Finite Volume Method)

有限体积法又称有限容积法、控制体积法,将求解域划分为有限的离散控制体积,对每个控制体积内部的平衡方程进行积分,从而得到一组离散方程,然后通过求解离散方程组得到近似解。

有限体积法具有良好的收敛性和稳定性,对边界条件的处理相对简单;相比于有限元法,对网格质量要求较低,更容易处理复杂的几何体和非均匀网格。

该方法主要应用于流体力学和热力学等领域。比如在流体力学中,可以用于求解不可压缩流体或可压缩流体的守恒方程,如Navier-Stokes方程等,常用于流体的流动模拟和分析。在进行流固耦合分析时,能够完美和有限元法进行融合。

云道智造伏图电子散热(Simdroid EC,点击文字可申请试用)是针对电子元器件、设备等散热的专用热仿真模块,采用有限体积法求解器,支持流热耦合计算,提供高精度的离散计算方法,同时结合电子散热相关行业经验,提供高保真的仿真模拟

4. 边界元法(BEM, Boundary Element Method)

边界元法是在经典积分方程法和有限元法基础上发展起来的一种数值方法,与有限元法在求解域内划分单元的思想不同,边界元法只在定义域的边界上划分单元,将边界积分方程离散化为线性代数方程组,通过求解这些方程组得到边界节点物理量。

边界元法只需要在边界上进行离散,降低了求解问题的维数,减少了计算的自由度,计算速度快、精度高;对无穷边界或边界条件占主导的问题具有较好的适用性,但较难应用在非均匀介质问题或大规模复杂问题。

作为一种准确而高效的计算方法,边界元法广泛应用于结构分析、电磁场分析、流体力学等工程领域。比如在结构分析中,边界元法在解决弹性、弹塑性、断裂力学等方面具有很大优势;在电磁场分析中,可以用于求解静电场、磁场、电磁场耦合等问题,提供各种电场特性如电势、电场强度、电流分布等信息。

数值方法不一而足,每种方法在求解不同类型问题时各具优缺点。除上述几种方法,还有物质点法(适合模拟涉及材料特大变形和断裂破碎等问题)、谱方法(适合计算流体力学复杂流场问题)、时域有限差分(适合求解电磁波场分布)、格子玻尔兹曼(求解CFD问题)、绝对坐标法(求解多体动力学问题)等。

各种数值方法并非只能单独使用,通常会结合实际问题选择适合的方法组合使用,以获得更精确、稳定、高效的解决方案。

工业软件也是如此。

云道智造自主研发的通用多物理场仿真PaaS平台伏图(Simdroid,点击文字可申请试用具备多种自主可控的通用求解器,支持多物理场耦合仿真。其中隐式结构、显式动力学、热、低频电磁、高频电磁等模块采用有限元法,流体、电子散热模块采用有限体积法,多体动力学模块采用绝对坐标法。下一步,将在高频电磁模块引入边界元法和时域有限差分。此外,采用有限元法和边界元法的声学模块、采用有限元法的电化学模块以及采用物质点法的物质点3D模块都在开发中。


参考资料:

  • 《微分方程数值求解——有限差分法》
  • 《有限元(FEM) 、有限差分(FDM)和有限体积(FVM) 的优势和劣势》
  • 《数值计算的六大方法》
  • 《一篇文章入门“求解器”开发(全篇)》
相关文章
ly~
|
27天前
|
存储 算法 编译器
游戏开发中,C 语言的性能优势体现在哪些方面?
在游戏开发中,C 语言凭借其对硬件的直接访问和内存操作的精准控制,能够显著提升性能。它允许开发者手动管理内存,优化数据存储和读取,充分利用显卡等硬件资源,实现流畅的图形渲染和音效处理。作为一种接近底层的语言,C 语言具有高效的执行速度,适用于物理引擎和碰撞检测等高性能需求模块,并且提供了丰富的运算符和数据类型,便于实现高效的算法。此外,C 语言代码具有良好的可移植性和跨平台性,支持多种操作系统和硬件平台,减少了多平台发布的开发成本。编译器提供的优化选项和手动代码优化的灵活性进一步提升了游戏的整体性能。
ly~
63 5
|
6月前
|
安全 物联网 5G
5g技术的优缺点是什么
5g技术的优缺点是什么
583 0
|
30天前
|
设计模式 监控 安全
Python多线程编程:特性、挑战与最佳实践【1】
Python多线程编程:特性、挑战与最佳实践【1】
28 0
|
30天前
|
设计模式 监控 安全
Python多线程编程:特性、挑战与最佳实践
Python多线程编程:特性、挑战与最佳实践
33 0
|
3月前
|
安全 编译器 C语言
深入分析宏定义的优势、局限性及最佳实践。
深入分析宏定义的优势、局限性及最佳实践。
91 14
|
3月前
|
开发框架 开发者 数据库管理
模块化开发和传统开发的优缺点有哪些
【8月更文挑战第26天】模块化开发和传统开发的优缺点有哪些
83 2
|
5月前
|
区块链 决策智能 UED
目前Layer2 解决方案有什么优缺点
目前Layer2 解决方案有什么优缺点
62 1
|
运维 监控 Java
03-揭秘大厂性能方案的奥秘!(下)
03-揭秘大厂性能方案的奥秘!
85 0
|
运维 监控 架构师
03-揭秘大厂性能方案的奥秘!(上)
03-揭秘大厂性能方案的奥秘!
85 0
|
存储 机器学习/深度学习 自然语言处理
转:算法与程序设计在文档管理系统中的优势与场景应用
算法与程序设计在文档管理系统中具有许多优势,这些优势使得系统更高效、安全、智能,并提供更好的用户体验,使得用户可以更方便地管理和利用文档资源。
81 0