一种基于完美匹配层的无限声场仿真方法技术

本发明专利技术涉及一种基于完美匹配层的无限声场仿真方法，包括以下步骤：在不考虑流动情况下，构建无源声场的亥姆霍兹方程；在计算域周围赋予完美匹配层区域，并使用复数坐标形式表示增加完美匹配层区域后的亥姆霍兹方程；对所述计算域采用三角形非结构网格进行分割，对所述完美匹配层区域采用四边形结构网格进行分割；对分割后的所述计算域和完美匹配层区域采用有限元方法进行离散，得到弱形式的亥姆霍兹方程，并进行求解得到声压分布曲线。本发明专利技术能够有效地解决虚假反射，计算成本高，难以实现等问题。等问题。等问题。

【技术实现步骤摘要】
一种基于完美匹配层的无限声场仿真方法


[0001]本专利技术涉及工业仿真
，特别是涉及一种基于完美匹配层的无限声场仿真方法。

技术介绍

[0002]对于无限声场问题，如水下结构的声辐射、声散射等问题，以及无限域中声传播等问题，通常采用声学边界元，声学无限元，声学有限元等方法，以上数值方法均需要对计算域进行截断处理。
[0003]现有比较主流的截断处理方法为完美匹配层(PML)方法，它是在计算域的周围增加一个吸收层，模拟各向异性材料，吸收辐射到该区域的散射场，计算域边界不会产生虚假反射。但在实际应用中，由于PML必须在感兴趣区域的有限距离上被截断，因此它的外部边界产生了人工反射。此时，对问题进行离散化求解通常会导致误差变大，虽然通过增加PML厚度等方法可以降低计算误差，但是这种方式会大大增加计算成本，导致有时候并不可行。

技术实现思路

[0004]本专利技术所要解决的技术问题是提供一种基于完美匹配层的无限声场仿真方法，能够有效地解决虚假反射，计算成本高，难以实现等问题。
[0005]本专利技术解决其技术问题所采用的技术方案是：提供一种基于完美匹配层的无限声场仿真方法，包括以下步骤：
[0006](1)在不考虑流动情况下，构建无源声场的亥姆霍兹方程；
[0007](2)在计算域周围赋予完美匹配层区域，并使用复数坐标形式表示增加完美匹配层区域后的亥姆霍兹方程；
[0008](3)对所述计算域采用三角形非结构网格进行分割，对所述完美匹配层区域采用四边形结构网格进行分割；
[0009](4)对分割后的所述计算域和完美匹配层区域采用有限元方法进行离散，得到弱形式的亥姆霍兹方程，并进行求解得到声压分布曲线。
[0010]所述步骤(1)中构建的无源声场的亥姆霍兹方程表示为：Δp+k2p＝0，其中，Δ为拉普拉斯算子，p为声场中的声压，k为波数。
[0011]所述步骤(2)中使用复数坐标形式表示增加完美匹配层区域后的亥姆霍兹方程具体为：其中，表示x方向映射系数，表示y方向映射系数，σ
x
表示完美匹配层区域在x方向的吸收函数，σ
y
表示完美匹配层区域在y方向的吸收函数，v
x
表示网格点在x方向到边
界上的法向距离，v
y
表示网格点在y方向到边界上的法向距离，a表示完美匹配层区域在x方向的最小坐标值，b
*
表示完美匹配层区域在y方向的最大坐标值，b表示完美匹配层区域在y方向的最小坐标值。
[0012]所述弱形式的亥姆霍兹方程表示为：其中，Ω
F
为计算域，Ω
A
为完美匹配层区域，表示计算域内网格点的声压，表示有限元中型函数，表示PML区域内网格点的声压，Γ表示计算域边界，g表示物理域的声波振幅对平面外法向量的求导，s表示边界长度离散单元。
[0013]所述完美匹配层区域在x方向的吸收函数所述完美匹配层区域在y方向的吸收函数c为声速。
[0014]有益效果
[0015]由于采用了上述的技术方案，本专利技术与现有技术相比，具有以下的优点和积极效果：本专利技术只需要选定边界，简单设置后可自动划分网格，外推出PML区域，构造的PML区域能够实现对声波的完全吸收，无虚假反射，且在离散化求解是不会令误差变大。
附图说明
[0016]图1是本专利技术实施方式的流程图；
[0017]图2是本专利技术实施方式中PML区域示意图；
[0018]图3是本专利技术实施方式中计算域和PML混合网格示意图；
[0019]图4是本专利技术实施方式的PML方法和特性阻抗方法计算结果对比图；
[0020]图5
‑
图8是本专利技术实施方式的PML方法与声学商业软件Comsol的结果对比图。
具体实施方式
[0021]下面结合具体实施例，进一步阐述本专利技术。应理解，这些实施例仅用于说明本专利技术而不用于限制本专利技术的范围。此外应理解，在阅读了本专利技术讲授的内容之后，本领域技术人员可以对本专利技术作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。
[0022]本专利技术的实施方式涉及一种基于完美匹配层的无限声场仿真方法，该方法应用PML方法求解线性声学中的波传播谐散射问题，即标量亥姆霍兹方程的求解。如图1所示，具体包括以下步骤：
[0023]步骤1，在不考虑流动情况下，构建无源声场的亥姆霍兹方程。该亥姆霍兹方程表示为：Δp+k2p＝0，其中，Δ为拉普拉斯算子，p为声场中的声压，k为波数。
[0024]步骤2，在计算域周围赋予完美匹配层区域，并使用复数坐标形式表示增加完美匹配层区域后的亥姆霍兹方程。具体地说：
[0025]采用复数坐标形式，
[0026]对于二维声场问题，如图2所示，在计算域Ω
F
周围赋予厚度为δ
PML
的PML区域Ω
A
，则亥姆霍兹方程可表示为：其中，表示x方向映射系数，，表示y方向映射系数，σ
x
表示完美匹配层区域在x方向的吸收函数，σ
y
表示完美匹配层区域在y方向的吸收函数，v
x
表示网格点在x方向到边界上的法向距离，v
y
表示网格点在y方向到边界上的法向距离，a
*
表示完美匹配层区域在x方向的最大坐标值，a表示完美匹配层区域在x方向的最小坐标值，b
*
表示完美匹配层区域在y方向的最大坐标值，b表示完美匹配层区域在y方向的最小坐标值。
[0027]综上可知，只要吸收函数为正值，任何传播波进入PML区域后都不会在边界出现反射，同时会迅速衰减。当PML区域足够厚或吸收函数值足够大时，传播波衰减非常快，以至于它们在PML区域外边界的反射非常弱；此外，反射波在通过PML区域返回计算域时进一步衰减，导致从PML区域外边界到计算域外边界的人工反射最小。
[0028]但是，当采用离散数值方法对问题进行求解时，这些优点没有一个被完全保留，因此必须对PML方法进行优化，以实现给定计算成本最佳可能解决方案。
[0029]步骤3，对所述计算域采用三角形非结构网格进行分割，对PML区域采用四边形结构网格进行分割。
[0030]本步骤中将计算域和PML区域采用不同的分割方式，形成了混合网格(见图3)。对于混合网格，计算域采用三角形非结构网格，PML区域采用四边形结构网格，计算精度取决于以下因素：1)单元函数的阶数；2)计算域内的单元大小；3)PML的单元层数；4)PML区域的厚度；5)吸收函数的形式。
[0031]步骤4，对分割后的所述计算域和完美匹配层区域采用有限元方法进行离散，得到弱形式的亥姆霍兹方程，并进行求解得到声压分布曲线。
[0032]本步骤中，对计算域和PML域采用有限元方法进行离散，亥姆霍兹方程方程可写为弱形式：
[0033][0034]其中，表示计算域内网本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于完美匹配层的无限声场仿真方法，其特征在于，包括以下步骤：(1)在不考虑流动情况下，构建无源声场的亥姆霍兹方程；(2)在计算域周围赋予完美匹配层区域，并使用复数坐标形式表示增加完美匹配层区域后的亥姆霍兹方程；(3)对所述计算域采用三角形非结构网格进行分割，对所述完美匹配层区域采用四边形结构网格进行分割；(4)对分割后的所述计算域和完美匹配层区域采用有限元方法进行离散，得到弱形式的亥姆霍兹方程，并进行求解得到声压分布曲线。2.根据权利要求1所述的基于完美匹配层的无限声场仿真方法，其特征在于，所述步骤(1)中构建的无源声场的亥姆霍兹方程表示为：Δp+k2p＝0，其中，Δ为拉普拉斯算子，p为声场中的声压，k为波数。3.根据权利要求2所述的基于完美匹配层的无限声场仿真方法，其特征在于，所述步骤(2)中使用复数坐标形式表示增加完美匹配层区域后的亥姆霍兹方程具体为：其中，表示x方向映射系数，表示y方向映射系数，σ
x
表示完美匹配层区域在x方向的吸收函数，σ...

【专利技术属性】
技术研发人员：钟策，陈赫，
申请(专利权)人：上海索辰信息科技股份有限公司，
类型：发明
国别省市：