COMSOL Multiphysics^®6.0 发布亮点

声学模块更新

COMSOL Multiphysics^®6.0 版本为“声学模块”的用户引入了新的压电波，时域显式多物理场接口、用于压力声学的物理场控制网格功能，以及流致噪声。请阅读以下内容，进一步了解这些更新及其他新增功能。

“压电波，时域显式”多物理场接口

在压电波，时域显式多物理场接口中，您可以使用在时域内对压电介质中波传播现象进行建模的新功能，可以对正、逆压电效应进行建模，并可以使用应变-电荷或应力-电荷形式来表示压电耦合。新增的接口利用新增的压电效应，时域显式多物理场耦合，将弹性波，时域显式接口与静电接口相耦合。

新增的接口基于间断伽辽金（dG 或 dG-FEM）法，并采用时域显式求解器。方程组的静电部分在每个时步通过经典有限元法 (FEM) 求解的代数方程组进行求解，这确保了一种计算效率非常高的混合算法，可以求解具有数百万个自由度 (DOF) 的大型模型。这种算法非常适合于集群架构上的分布式计算。您可以在更新的使用压电换能器的超声波流量计和角钢梁无损检测教学案例中查看这一新接口的应用演示。

压力声学的物理场控制网格

物理场控制网格功能现在可用于压力声学，频域和压力声学，瞬态接口。生成的网格遵循有关解析波的最佳实践，确保有足够数量的网格单元。此外，完美匹配层 (PML) 采用结构化网格进行网格划分，周期性条件采用复制网格操作，单个边界层网格用于外场计算。

对于频域，最大频率是自动从研究中选取的。对于时域，最大频率取自物理场设置，确保网格的空间分辨率和求解器时步的时间分辨率一致。对于所有其他声学接口和固体力学接口，仅处理 PML 和周期性条件。所有适用新增的物理场控制网格划分的教学案例都已更新。

新增流致噪声功能

新版本介绍了一种用于流致噪声建模的混合计算气动声学 (CAA) 方法。这种方法基于湍流源与声学方程之间的单向耦合，假设从声场到流场不存在反向耦合，计算方法基于 Lighthill 声学类比（波动方程）的有限元离散化。这种方程形式可以确保已在底层包含固体边界（可以是固定的或振动的）。有两个流致噪声选项可用：Lighthill 类比和较为简单的气动声学波动方程 (AWE) 类比。

新功能依赖于将使用“CFD 模块”求解的大涡模拟 (LES) 流体流动模型与压力声学，频域中的气动声学流动源域特征相耦合，这是通过使用气动声学流动源耦合多物理场耦合和专门的瞬态映射研究来实现。

新增两个高频压力声学物理场接口

两个新增的物理场接口基于高频假设，并基于基尔霍夫-亥姆霍兹积分。第一个是压力声学，渐近散射接口，专门用于散射建模，另一个是压力声学，基尔霍夫-亥姆霍兹接口，主要用于辐射建模。

压力声学，渐近散射

压力声学，渐近散射接口用于模拟高频散射。假设声场为局部平面，因此散射场可以通过解析方式来表示。通过定义面法向阻抗、反射系数或吸声系数，您可以将散射物体表面视为完全反射或具有吸声属性，后两者取决于入射角。该接口可以对球面波和平面波的散射进行建模，其中具有内置功能，可以采用简单的角度考量和更高级的半立方体方法来计算可见性因子。您可以在潜艇高频渐近散射教学案例中查看这一特征的应用演示。

压力声学，基尔霍夫-亥姆霍兹

压力声学，基尔霍夫-亥姆霍兹接口用于对高频辐射进行建模。假设声场为局部平面，这种技术也经常被称为高频边界元法或简称为 HFB 法，通常用于计算高频振动结构的辐射声场，无需对周围流体进行建模。只要流体中的声波波长小于结构和结构模态，这种方法就有效。对于平坦的振动表面，该方法简化为计算瑞利积分。接口中提供用户定义的选项用于访问完整的基尔霍夫-亥姆霍兹积分公式，定义压力及其法向梯度。

新增的压力声学域分解方法

现在您可以使用域分解 (Schwarz)方法求解大规模压力声学（亥姆霍兹）问题，其中使用移位拉普拉斯方法以及与非重叠 Schur 方法相同的用于内部重叠边界的吸收边界条件。这种方法的优点是多重网格可以用作域求解器，并且域分解方法不需要粗化栅格。

集总扬声器边界和内部集总扬声器边界

对于许多电声学应用，我们可以通过将集总 Thiele–Small 表示与有限元法相结合来有效地对扬声器进行建模。电机的电磁部件采用电路接口建模，而声学部分采用压力声学，频域接口求解。该方法假设扬声器膜片振动可以通过活塞运动来描述，并且对微型扬声器特别有意义。内部集总扬声器边界包含空气对膜片两侧的影响。集总扬声器边界可以选择包含通过阻抗的后腔体的柔度。您可以在扬声器驱动器集总模型和戴在仿真人耳上的耳机教学案例中查看这些新增特征的应用演示。

外场计算的扇区对称选项

外场计算特征已扩展，新增了选项，用于处理包含扇区对称的模型。其中两个选项扩展了现有的对称平面功能：扇区对称选项和具有一个对称平面的扇区对称选项。在对放置在无限挡板中的扬声器驱动器进行建模时，后者尤其重要。对于高级扇区对称模型，外场分析还可以通过方位角模数进行扩展。

使用新的外场变量进行优化

在三维中，新的外场变量可用于基于梯度的优化问题，如形状优化或拓扑优化。优化目标现在可以定义为在外场中计算的变量，例如定义辐射方向图或离轴响应的值。该变量仅适用于外场计算特征中的对称平面选项。新变量通过将 _opt 附加到现有外场变量来定义：压力的运算符是pext_opt(x,y,z)，声压级的运算符是Lp_pext_opt(x,y,z)。您可以在矩形扬声器喇叭的形状优化 - 三维教学案例中查看这一更新功能的应用演示。

磁力学多物理场接口

新版本新增了两个物理场接口，用于分析耦合的磁效应和机械效应：磁力学和磁力学，无电流。典型的应用场景是磁场引起固体变形，或者反过来，移动结构改变磁场，例如，这是变压器嗡嗡声的物理解释。新的接口依赖于新的磁机械力多物理场耦合。对某些类型的声学换能器（如平衡电枢换能器）进行建模时，这种耦合特别有意义。这一新功能还需要“AC/DC 模块”，您可以在平衡电枢接收器微型扬声器教学案例中查看相关应用演示。

完美匹配边界 (PMB) 辐射条件

新增的完美匹配边界特征是以辐射条件的形式应用于开放边界的 PML，无需定义域，例如作为几何模型中的一个层。该条件通过 COMSOL Multiphysics^®的额外维度功能自动应用 PML 公式，这也简化了对辐射边界的要求，原因是，它原则上可以具有任何凸面形状。您可以使用不同的选项来控制衰减方向。新增的边界条件可用于所有相关空间维度的压力声学，频域接口。您可以在以下教学案例中查看这一新特征的应用演示：

• piezoacoustic_transducer
• diffuser_schroeder_2d
• acoustic_scattering

“线性欧拉，边界模式”物理场接口

线性欧拉，边界模式接口用于计算和识别存在可以通过理想气体流很好地近似处理的稳态背景平均流的情况下，波导和导管中的传播和非传播模式。该接口对波导的边界、入口或横截面执行特征模态分析。

“线性纳维-斯托克斯，边界模式”物理场接口

线性纳维-斯托克斯，边界模式接口用于计算和识别存在任何稳态等温或非等温背景平均流的情况下，波导和导管中的传播模式和非传播模式。该接口对波导的边界、入口或横截面执行特征模态分析，其中考虑所有热损耗和黏滞损耗效应以及与背景流的相互作用，并根据需要包含声边界层损耗。

线性欧拉和线性纳维-斯托克斯的面外和边界模式分析

在线性欧拉和线性纳维-斯托克斯接口中，您现在可以分别在二维和二维轴对称中添加可选的面外波数和周向波数。在相同的空间维度中，您现在可以使用模式分析研究为这些类型的物理场接口建立 2.5D 仿真。

“弹性波，时域显式”中的计算位移后处理特征

新版本的弹性波，时域显式物理场接口中新增了一个后处理特征，称为计算位移。利用此特征，您可以通过求解一组辅助常微分方程，以最佳方式计算点、沿边、边界或域中的位移。这些新特征作为子特征添加到材料模型中，例如弹性波，时域显式模型或压电材料模型。此特征不影响结果，但仅用于后处理，并生成场变量，用于对位移进行可视化和后处理。由于此特征会添加并求解额外的方程，因此使用它需要额外的计算资源。您可以在各向同性-各向异性试样：弹性波传播教学案例中查看这一新功能的应用演示。

“非线性压力声学，时域显式”中的最小和最大压力后处理特征

非线性压力声学，时域显式接口中新增了一个后处理特征，称为计算最小和最大压力。此特征计算域或边界上随时间和空间变化的最大和最小压力。软件会自动创建两个变量：nate.p_min 和 nate.p_max，可用于后处理，例如，计算聚焦区的大小。您可以在高强度聚焦超声 (HIFU) 在组织体模中的传播教学案例中查看这一特征的应用演示。

用于振动声学仿真的“对，多物理场耦合”

“声学模块”中新增了两个多物理场耦合，用于耦合装配中的声学域和固体域。这是作为“对”条件实现的，支持在装配中两个域之间的界面处使用非共形网格。第一种类型的多物理场耦合，即对，声-结构边界，在装配中用于将压力声学，频域或压力声学，瞬态接口耦合到固体力学接口。第二种类型的多物理场耦合，即对，热黏性声-结构边界，在装配中用于将热黏性声学，频域或热黏性声学，瞬态接口耦合到固体力学接口。这种耦合使用罚函数公式在时域中进行更有效的计算。您可以在模拟同时作为信号发射器和接收器的压电器件教学案例中查看这些新增耦合的应用演示。

热黏性声学的重要更新

对于热声学建模，新版本引入了多个新的和改进的特征。

• 对于端口，新增的平面波选项可用于处理滑移和绝热案例或者非波导模型。先前软件版本中提供的波导模式名称已更新为数值 (0,0) 模式和循环 (0,0) 模式。
• 所有热黏性接口中使用的默认离散化已从拉格朗日单元更改为巧凑边点单元。这将在求解使用结构化网格的模型时带来重要的性能提升。
• 改进了求解器设置，用于求解使用非线性热黏性声学贡献的模型。这些设置对于高度非线性的问题尤为重要。如需在 COMSOL Multiphysics^®5.6 版本的模型中使用这些新设置，您需要将求解器重置为其默认设置。
• 在曲面情况下，热黏性声学与结构力学之间的耦合通过使用“节点”约束类型得到了改进。
• 在边界层网格特征中，新增了厚度明细选项所有层。这样可以更容易地设置跨越整个热黏性声学边界层的网格。
• 在非线性热黏性声学贡献特征中将密度近似设置为二阶的情况下，如果使用热力学系统特征指定空气材料（潮湿或干燥），现在会自动计算材料属性。
• 在热黏性声学，频域接口中，您现在可以为二维轴对称模型添加可选的周向波数，从而可以对复杂传播模式进行高级 2.5D 分析。

您可以在以下模型中查看这些特征的应用演示：

• transfer_impedance_perforate
• nonlinear_slit_resonator
• wax_guard_acoustics

射线声学新增功能

根据方向性释放

在射线声学接口中，新增的带方向性的源节点现在可以用于根据用户定义的空间方向性函数释放具有初始强度或功率的射线分布。在射线声学仿真中定义扬声器源时，这一点尤其重要。

根据外场计算释放

您可以使用新增的根据外场计算释放节点，根据先前研究中的外场计算特征，发射具有强度和相位分布的射线。这有助于多尺度声学仿真，将基于网格的近场解与更远距离的射线追踪仿真相结合。您可以在超声波停车传感器教学案例中查看这一特征的应用演示。

从文件中加载射线坐标时的转变

使用根据数据文件释放节点从文件中加载射线释位置时，您现在可以对初始坐标应用转变。您可以使用膨胀（缩放）、旋转和平移的任意组合。如果初始射线方向也是从文件中加载的，您可以选择对位置和方向应用相同的旋转。

新增的衰减模型

在软件 6.0 版本中，有四种不同的方法来定义衰减系数，控制射线在仿真域中传播时射线强度功率如何降低。您可以指定压力幅值衰减系数或强度幅值衰减系数，单位为奈培/米；也可以输入幅值衰减系数，单位为分贝/单位长度或分贝/波长。衰减系数的所有这些定义也可用于对几何形状之外的空区域中的衰减进行建模。

非局部耦合的简化名称

射线声学接口定义了计算模型中射线表达式的总和、平均值、最大值或最小值的耦合。在 6.0 版本中，这些耦合的名称已简化，以便于使用。

下表列出了新旧耦合名称。

耦合描述	旧名称	新名称
对射线求和	rac.racop1(expr)	rac.sum(expr)
对所有射线求和	rac.racop_all1(expr)	rac.sum_all(expr)
对射线求平均值	rac.racaveop1(expr)	rac.ave(expr)
对所有射线求平均值	rac.racaveop_all1(expr)	rac.ave_all(expr)
对射线求最大值	rac.racmaxop1(expr)	rac.max(expr)
对所有射线求最大值	rac.racmaxop_all1(expr)	rac.max_all(expr)
对射线求最小值	rac.racminop1(expr)	rac.min(expr)
对所有射线求最小值	rac.racminop_all1(expr)	rac.min_all(expr)
在射线的最大值下计算	rac.racmaxop1(expr, evalExpr)	rac.max(expr, evalExpr)
在所有射线的最大值下计算	rac.racmaxop_all1(expr, evalExpr)	rac.max_all(expr, evalExpr)
在射线的最小值下计算	rac.racminop1(expr, evalExpr)	rac.min(expr, evalExpr)
在所有射线的最小值下计算	rac.racminop_all1(expr, evalExpr)	rac.min_all(expr, evalExpr)

旧名称在 6.0 版本中仍然有效，因此无需更新任何现有模型。

带扭转的轴对称

在固体力学接口的二维轴对称中，现在可以包含周向变形，这可以通过在物理场接口的轴对称近似栏中选中包含周向位移复选框来启用。利用该选项，您可以更高效地对轴对称结构的扭转等现象进行建模。

新增的阻尼模型

新版本为机械材料模型添加了新的阻尼模型：

• 波衰减模型本质上是一个黏滞模型，但其参数由材料中弹性波衰减的测量数据给出。它可以在固体力学的线弹性材料中找到。
• 最大损耗因子模型主要用于材料的时域分析，对于这些材料，损耗因子表示在频域中给出了良好的描述。这种阻尼模型适用于支持黏滞阻尼的所有材料模型。
• 在压电材料特征中，除了机械阻尼最大损耗因子之外，还有一种新的介电损耗频域阻尼模型：复介电常数。
• 对于电荷守恒，压电，您现在可以添加两个新的色散模型：德拜和多极德拜。

声学问题后处理新增功能

对于后处理和可视化，新版本包含以下新增功能：

• 射线声学接口的脉冲响应绘图特征中新增了离散傅里叶变换 (DFT)和移动平均选项。
• 当在显示更多选项中启用绘图信息栏后，现在可以使用新增的信息栏进行绘图。新增的栏显示绘图时间，这在射线声学模型或使用新增的基于基尔霍夫-亥姆霍兹的高频接口模型中可能非常重要。
• 在倍频带图中，绘制的物理量现在可以选为连续功率谱密度，频带功率或频带平均功率谱密度。
• 能量衰减子特征中添加了一条警告消息，电平衰减现象偏离正常值时会显示此消息。
• 新增了Thermal Wave颜色表，并将其用作热黏性声学和线性纳维-斯托克斯接口中温度变化图的新默认值。
• 新增了Wave和Wave Light颜色表，并将它们用作表示声压的所有绘图的默认值。
• 更快地生成大型表面图，这一点在地震波在地球内部的传播教学案例中非常明显。
• 接收器数据集的预览速度现在要快得多。

声学模块的其他重要增强功能和更新

• 改进了基于 dG-FEM 的时域显式接口的阻尼公式。新公式更稳定、更精确、性能更高。
• 改进了热黏性声学模型中集总圆近似的处理方法，这适用于狭窄区域声学和内部穿孔板。
• 映射研究结合背景流体流动耦合多物理场特征执行的映射的稳定性得到了改进。对于具有导入几何的模型以及在弯曲边界上具有 CFD 边界层网格的模型，这一点尤其明显。
• 可压缩势流物理场接口现在适用于 PML 域，这简化了背景流的设置，以便在线性势流接口中进行后续分析。

新的和更新的教学案例

COMSOL Multiphysics^®6.0 版本的“声学模块”引入了多个新的和更新的教学案例。