研究和求解器更新

COMSOL Multiphysics® 5.3 版本新增了一个用于 CFD 仿真的求解器,以及一个用于电磁和腐蚀边界元法仿真的求解器。请阅读以下内容,了解 COMSOL Multiphysics® 5.3 版本中有关研究和求解器的所有更新。

CFD 的代数多重网格 (AMG) 求解器

平滑聚集代数多重网格 (SA-AMG) 方法已扩展为与 COMSOL Multiphysics® 中 CFD 的以下专用平滑器配合使用:SCGS、Vanka 和 SOR 线。

使用几何多重网格 (GMG) 求解器时需要考虑三个级别的网格,对于复杂几何和尺寸,在网格剖分和求解时经常会出现问题。而 SA-AMG 求解器只需要一个网格层,对于大型问题和“复杂”几何来说,这使网格剖分过程变得简单许多,使求解过程变得稳定许多。

例如,在太阳能电池板的流-固耦合模型(见图)中,支撑电池板的支杆和梁与周围的空气域相比较小。这个尺寸差异导致很难对空气域以及较小的零件和组件进行网格剖分,如果要创建三个不同尺寸的网格,将更加困难。SA-AMG 求解器只需要一个网格层,更容易获得。

自适应与网格剖分序列和误差估计集成

修订了稳态、参数化和特征值问题的 h 自适应算法,从而可以节省中间解和网格。此外,得到的自适应网格现在以网格剖分特征为基础,可以在需要时从自动求解过程无缝切换到手动自适应。

为了实现这一点,引入了两个网格剖分特征:适应尺寸表达式适应 特征根据解中的错误表达式或所需网格单元尺寸的表达式来细化网格。或者,通过将尺寸表达式 节点添加到“模型开发器”的网格剖分序列,可以使用表达式改变整个建模空间的单元尺寸。请参阅网格更新页面,了解更多详细信息。

适应功能和误差估计也已统一,这样,适应方法使用的误差估计现在可用作对结果进行后处理的因变量。此外,现在也可以对(偏微分方程残差)的 L2 误差估计进行后处理。

功能提升进一步促使网格初始化 方法现在能够对网格执行自适应操作,而不仅仅是三角形和四面体类型。之所以能够实现这一点是因为自适应网格是通过参考 特征构建的,同时保留了原始网格序列以及新的尺寸表达式 网格剖分特征。

如有需要,常规细化最长细化 方法现在还可以自动将网格转换为三角形或四面体。这意味着用户不需要向网格剖分序列添加任何转换 特征。

COMSOL Multiphysics 5.3 版本中使用自适应网格细化功能的基准模型。 这一 Euler Bump 基准模型现在使用 COMSOL Multiphysics® 中新增的自适应网格细化功能。
这一 Euler Bump 基准模型现在使用 COMSOL Multiphysics® 中新增的自适应网格细化功能。


有关使用新的自适应网格细化的 App,请访问以下“案例库”路径:
CFD_Module/High_Mach_Number_Flow/euler_bump

边界元法问题的快速迭代求解器

密集直接求解器现在可用于求解最适用于通过边界元法 (BEM) 求解的问题,这对于不太适用于有限元法 (FEM) 建模的 App 非常有用。

该直接求解器的求解时间大约与问题中自由度 (DOF) 的立方成正比。换句话说,求解时间会随着问题尺寸的增大而显著增加。为了缓解这一状况,根据快速矩阵-矢量乘法为迭代求解器提供了支持,从而能够使用所谓的 ACAACA+ 压缩来压缩矩阵。这些替代方案对应于两种不同版本的自适应交叉近似方法,这是基于远场近似的一种快速矩阵乘法方法。

软件提供了两个预条件器:稀疏近似逆 (SAI)直接预条件器,这两个预条件器都可用于所谓的矩阵的近场部分。矩阵的近场部分是稀疏的,与全矩阵相比,存储和求解所需的内存少很多。SAI 预条件器是一个显式预条件器示例,近似于矩阵逆,而不是矩阵本身。直接预条件器使用更常见的矩阵 LU 分解。

BEM 已在用于求解偏微分方程的常规物理场接口中实现,用于求解 AC/DC 模块中的静电 App,还用于求解腐蚀模块电镀模块中的电化学电流密度 App。

使用 COMSOL Multiphysics 5.3 版本中的 BEM 模拟的石油钻塔,重叠部分显示牺牲阳极的放大图。 使用边界元法 (BEM) 对海水中石油钻塔的静电属性进行建模。诸如尺寸、零件数和几何的常规复杂度以及钻塔所在的无限区域,这些参数都使这样一个示例优先选用 BEM 进行建模。重叠部分是石油钻塔的部分放大,显示了诸如牺牲阳极(大型钻塔结构旁的细杆)等更详细的部分。
使用边界元法 (BEM) 对海水中石油钻塔的静电属性进行建模。诸如尺寸、零件数和几何的常规复杂度以及钻塔所在的无限区域,这些参数都使这样一个示例优先选用 BEM 进行建模。重叠部分是石油钻塔的部分放大,显示了诸如牺牲阳极(大型钻塔结构旁的细杆)等更详细的部分。


有关使用 BEM 的示例,请访问以下“案例库”路径:
ACDC_Module/Capacitive_Devices/capacitor_tunable
ACDC_Module/Applications/touchscreen_simulator
ACDC_Module/Tutorials/capacitive_pressure_sensor_bem

混合 BEM/FEM 问题的求解器支持

有时,多物理场问题可以通过只一种数值方法来求解,然而,针对不同物理场使用不同数值方法(边界元法 (BEM) 和有限元法 (FEM))有可能进行最佳求解。如果矩阵存储是 FEM 部分的最佳稀疏格式,并且是 BEM 部分的密集或无矩阵格式,则可使用混合 BEM/FEM 模型。这样可以对矩阵的各个 FEM 和 BEM 部分使用单独的预条件器/平滑器。

例如,可以使用带有混合预条件器的有效迭代求解器。FEM 部分可以像往常一样自由地进行预处理,而 BEM 部分可以与上述预条件器之一配合用于近场矩阵。迭代方法通过基于混合矩阵的方法/无矩阵方法来计算残差,从而最大限度地利用不同种类的快速矩阵-矢量产品。

精确边界通量变量的灵敏度

现在可以使用“正演”灵敏度方法从边界通量 变量获得灵敏度贡献,这些是可用于某些物理场接口(例如,质量和传热)的精确边界通量变量。对于这些接口,在其各自的离散化 栏下方选中计算边界通量 复选框以访问和使用这些变量。

组合解

可以将两个解对象组合为一个解或数据集。当后处理需要一个解/数据集时或者需要一个解作为新仿真的输入时,这非常有用。可以组合瞬态、参数化和特征频率解,而这些解可以串联起来,也可以求和。


间断 Galerkin 方法的基于网格的性能改进

实现了多项改进,以加速间断 Galerkin (dG) 方法并减少内存占用。其中一项改进是引入一种新的网格度量方法,用于计算所使用的显式时间步进方法的稳定时间步。该度量值是三角形最大内切圆和四面体最大内切球的直径。这个度量方法可以更好地确定稳定的时间积分中需要什么时间步,同时更好地描述手边用于仿真的网格单元的特征。

通过新的网格质量优化程序实现了第二项改进。此操作过程应与 dG 方法配合使用,以进一步增加显式时间步进方法的稳定时间步。此方法可更改网格,避免单元格过小,否则会限制稳定的时间步。在三维中生成四面体网格时,请使用新的避免过小的单元 网格选项(见图)。


例如,请参阅包含 750 万个自由度 (DOF) 的“通用传播时间法超声波流量计”教学模型。在采用 4 核 Intel Core™ i7 处理器(频率为 3.60 GHz)且内存为 32 GB 的台式计算机上进行测试运行时,在 COMSOL Multiphysics® 5.2a 版本中求解声学问题需要 7 小时零 5 分钟,需要 6.0 GB 内存。在 5.3 版本中,使用新的稀疏装配方法,相同的研究求解只需 5 小时零 1 分钟,需要 5.8 GB 内存。这表示速度提升了约 30%,占用的内存也略微减少。

间断 Galerkin 方法基于多核的性能改进

在多核系统上运行模型时,已经减少了内存占用量,其中使用了残余矢量的新稀疏装配方法。所需内存减少,而且不取决于所使用的 CPU 内核数。此外,初始化期间所需的内存已大大减少。这项改进加速了此方法的运行速度,因为避免了不必要的内存分配。

进一步研究使用多核求解时的内存改进时,除了网格剖分和网格度量参数带来的改进,我们可以比较在 COMSOL Multiphysics® 5.2a 和 5.3 版本中都求解过的、所使用的网格和时间步完全相同的模型。这里分析了使用立方基函数的椭球内的简单波动方程仿真,使用 6 个内核在 3.60 GHz 的 Intel® Xeon® CPU E5-1650 v4 上进行比较。即使没有网格改进,CPU 时间也减少约 18%。使用更多核心时,预计内存减少会更加明显。

Size Version 5.2a Version 5.3 Improvement
Medium (6.7 MDOF/t = 0.05) 74 sec./4.1 GB 61 sec./3.2 GB 18%/22%
Large (20 MDOF/t = 0.05) 307 sec./10 GB 250 sec./7.3 GB 19%/27%


Intel、Intel Core 和 Xeon 是 Intel Corporation 或其子公司在美国和/或其他国家/地区的商标。