基于单隙(Single-Clearance)模型估算随机堆积球形颗粒吸声系数的理论方法

《Modelling》：Theoretical Estimation of Sound Absorption Coefficients for Randomly Packed Spherical Granules Using Single-Clearance Model

　　

摘要：本研究旨在建立一种简易的理论方法，用于估算随机堆积颗粒材料(randomly packed granular materials)的吸声特性。研究人员基于已有规则堆积数学模型及利用显微计算机断层扫描(micro Computed Tomography, CT)图像的方法，提出"单隙(single?clearance)模型"—一种可估算吸声系数的理论模型。该模型依据颗粒粒径(particle size, d)与堆积密度(packing density, φ)计算颗粒材料中空隙体积(void volume, Vair)与球体表面积(sphere surface area, Sbeads)，并据此将颗粒材料的堆积结构简化为两平行平面间的隙层(clearance between two planes)。以孔隙率(porosity, 1?φ)沿声波传播方向变化较小的随机堆积结构为对象，考虑该变化对吸声系数影响可忽略不计。在单隙模型中，采用随机紧密堆积(random close packing)堆积密度φ＝0.65计算所得吸声系数，在各粒径下均与基于CT图像推导的理论值基本吻合。因此，仅需粒径与堆积密度即可简易估算吸声系数。

论文解读：《Modelling》期刊所刊此文针对随机堆积球形颗粒吸声性能预测长期依赖CT扫描或复杂数值模拟、难以工程实用的问题，提出可将整体随机堆积等效为单一平行板隙的"单隙(Single?Clearance)模型"，只需粒径(d)与堆积密度(φ)即可估算法向入射吸声系数(normal?incidence sound absorption coefficient, α)，通过与实验值及CT图像理论值对比验证其有效性，为颗粒吸声材料的设计与快速预估提供了简便理论工具。

一、主要关键技术方法概述

研究人员选用粒径d＝4、2、1、0.8、0.6、0.5、0.4、0.3、0.2、0.1、0.05 mm玻璃微珠(glass beads)随机填入内径29 mm、厚20 mm铝合金阻抗管制成随机堆积试样；以Brüel & Kj?r 4206双传声器阻抗管按ISO 10534?2测定法向入射吸声系数作实验基准；理论上先依总体积V_t、堆积密度φ算颗粒体积V_beads＝φV_t及空隙体积V_air＝(1?φ)V_t，再由单球半径r得S_beads＝(3φV_t)/r，等效平行板隙厚度b＝2V_air/S_beads；引入曲折度(tortuosity, α_∞)＝1.45，按Stinson与Allard方法由等效隙内黏性气流Navier?Stokes方程求传播常数γ与特征阻抗Z_c，用一维传递矩阵法(transfer matrix method)算α；另以同批次试样CT分层图像逐层算局部V_air、S_beads后分层传递矩阵乘积累加获CT理论值作参照，单隙模型则以整体平均V_air、S_beads只算一次传递矩阵；分别取φ＝0.65、0.6366(Random Close Packing理论值)、0.60、0.55作参量敏感性分析。

二、研究结果

1. 2.
   Samples and Measurement Apparatus（样品与测量装置）
   列出所用11种粒径玻璃微珠试样规格（表 1）及B&K 4206双传声器阻抗管设备参数（表 2），确认实验按ISO 10534?2标准执行并获可靠法向入射α实测值。
2. 3.
   Theoretical Analysis（理论分析）

3.1. Calculation of Void Volume and Sphere Surface Area（空隙体积与球体表面积的计算）

由样品管总体积V_t与堆积密度φ，导出V_beads＝φV_t、V_air＝(1?φ)V_t；结合单球几何关系S_sphere＝4πr2、V_sphere＝(4/3)πr3推得整体颗粒表面积S_beads＝3φV_t/r。此二全局几何量为后续等效隙建模基础。

3.2. Simplification to Clearance Between Two Planes（简化为两平面间隙层）

将随机堆积多孔结构等效为两无限大平行刚性壁间隙层，隙厚b由等效条件b＝2V_air/S_beads给出，使等效模型同时匹配原堆积的空隙体积与内表面积——圆管或三角管等效模型无法二者兼顾，而平行板隙可满足。

3.3. Normal?Incidence Sound Absorption Coefficient（法向入射吸声系数计算）

在等效平行板隙中视空气为黏性层流，考虑曲折度α_∞＝1.45修正有效密度ρ_s，据Stinson对矩形截面管道及Allard含曲折度处理方法由三维Navier?Stokes方程得等效传播常数γ＝jω√(ρ_sC_s)与特征阻抗Z_c＝√(ρ_s/C_s)，代入一维平面波传递矩阵求得刚性背衬下法向入射吸声系数α。

3.4. Differences Between Analytical Methods（各分析方法差异）

对比三种途径：规则堆积——沿声传播方向细分多层各自算V_air、S_beads再乘各层传递矩阵；随机堆积CT法——数百至上千张垂直切片逐层提取局域V_air、S_beads分层矩阵连乘；单隙模型——全体积取平均V_air、S_beads只建单一等效隙算一次传递矩阵，忽略沿程截面几何波动。

3.5. Verification of Effects of Simplification（简化效果验证）

对d＝4、2、1 mm分析规则与随机堆积沿厚度方向孔隙率与单位体比表面积(specific surface area per m3)分布：规则堆积（尤其是面心立方FCC、六方最密堆积HCP）因周期性排列致孔隙率沿深呈周期性波动，变异系数(coefficient of variation, CV)较大；随机堆积孔隙率与比表面积沿深波动很小，CV远小于规则堆积。理论对比显示：规则堆积取平均几何参数会使α明显偏移逐层计算结果，随机堆积取整体平均值与逐层CT法所得α几乎重合。表明随机堆积沿声传播方向几何参数均一性好，整体简化为单一等效隙对吸声系数计算误差可忽略。

1. 4.
   Comparison of Experimental and Theoretical Values（实验值与理论值比较）
   对全部粒径分别以φ＝0.65、0.6366、0.60、0.55运行单隙模型并与实验值及CT理论值比对：φ＝0.65所得理论曲线在全粒径范围内与CT推导值最接近；大粒径(d≥1 mm)φ＝0.65理论值亦贴近实验值，小粒径(d≤0.1 mm)较低φ(0.60或0.55)拟合更佳——归因于小管壁及装填表层疏松致实测整体φ偏低但中心区近随机紧密堆积，且微细颗粒接触点极窄缝隙(<边界层厚度)占比增大使单隙均匀化低估有效通气隙厚，低φ设定可部分补偿。Knudsen数(Kn＝λ/L, λ为空气分子平均自由程, L为特征长度即隙厚b)分析指出极窄接触区Kn＞0.01偏离纯黏性流假设，属待修正因素。总体而言单隙模型用φ≈0.65可便捷预估随机球形颗粒吸声性能。

三、讨论与结论翻译

研究人员提出单隙(single?clearance)模型用以理论估算随机堆积颗粒材料吸声系数，该法基于粒径与堆积密度求整体空隙体积与球体表面积并将堆积等效为单一平行板隙后借传递矩阵求取法向入射吸声系数。经与需CT图像逐层分割传递矩阵叠乘之既有方法及多粒径实验值对照得出：(1) 随机堆积结构沿声波入射方向孔隙率变动微小，忽略该变动对吸声系数影响可忽略，故可沿声传播方向整体简化为单一平行板隙；(2) 单隙模型中取随机堆积密度φ＝0.65计算之吸声系数于所有测试粒径下均与CT图像理论值基本一致，表明仅凭粒径与堆积密度即可简易估算吸声系数；(3) 该单隙模型预期亦可推广至其它能方便算出孔隙率与表面积且孔隙率沿深波动小的多孔材料。