Acoustic field modeling and measurement of the parallel graphene based thermo-acoustic actuator
-
摘要: 主动流动控制技术是提升航空航天飞行器气动特性的有效手段，其中激励器是这些控制技术的核心组件，如等离子体激励器、合成射流激励器以及振荡射流激励器等。本文提出了一种基于热声效应的并联石墨烯热声激励器，该激励器具有结构简单、输入功率低、控制频率宽以及结构适应性强等优点，能够适应飞行器各种复杂曲壁面安装环境和变飞行工况，具有良好的应用前景。具体而言，该热声激励器利用石墨烯材料比热容极低和热导率高的特性，通过焦耳加热原理向周围空气辐射周期性的声场以进行声激励控制。首先，通过热声理论对石墨烯热声激励器进行声场建模，建模中加入了组合声源叠加原理，优化了声波相位差和声场指向性计算方法；其次，改进了石墨热声烯激励器的电路连接方式，采用并联石墨烯薄膜有效提升了声场声压幅值；最后，搭建了石墨烯热声激励器的声压测试平台，验证并研究了并联石墨烯热声激励器的声场特征，分析了输入功率、频率、测试距离等因素的影响。Abstract: The active flow control technology is an effective method to improve the aerodynamic characteristics of aerospace vehicles, in which actuators are the core of these control technologies such as plasma actuator, synthetic jet actuator and oscillating jet actuator. Present study proposes a parallel graphene actuator based on the thermo-acoustic effect. The actuator has the advantages of simple structure, low input power, wide controlling frequency and strong structural adaptability. It can adapt to various complex curved wall installations of aircraft environment and variable working conditions, resulting in a good application prospect. Specifically, the thermo-acoustic actuator utilizes the extremely low specific heat capacity and high coefficient of heat conductivity characteristics of graphene materials. Through the Joule heating principle, the graphene films can radiate a periodic sound field to the surrounding air for acoustic excitation control. According to this principle, firstly this paper uses the thermo-acoustic theory to model the sound field of the graphene actuator, adds the principle of combined sound source superposition to the modeling, and optimizes the computation method of acoustic wave phase difference and sound field directivity. Secondly, it improves the circuit connection method of the parallel graphene actuator films so that it effectively increases the sound pressure amplitude of the sound field. Finally, by establishing a sound pressure test platform in a semi-anechoic chamber room and analyzing the influence of input power, frequency, test distance and other factors on output sound pressure, this paper studies and verifies the sound field of the parallel graphene thermo-acoustic actuator.
-
0 引 言
随着“工业4.0”和“中国制造2025”的提出，增材制造技术在航空、航天和汽车等领域的应用不断深入，以铝合金为代表的轻质合金材料需求快速增长。铝合金具有资源丰富、密度低、比强度较高的优点。2019年SmarTech报告指出，2018年增材制造铝合金粉出货量增长43%，且未来10年将以每年近50%的速度增长。
增材制造要求铝合金粉末具有粒径集中度高（15～53 μm和53～150 μm）、球形度高、含氧量低、无卫星粉和空心粉等特性。金属液体雾化是增材制造用金属粉末的主要制备方式[1]，主要基于气雾化、等离子体旋转电极雾化[2-3]和等离子体雾化[4]等原理。采用这些方法制备的金属粉末粒径大、粒径分布宽、球形度差且含有卫星粉和空心粉，无法为增材制造提供高性能粉末。与气雾化和旋转电极雾化相比，采用转盘雾化技术制备的金属粉末具有粒径小、粒径集中度高、球形度高、无卫星粉和空心粉、成本低廉等优势，更加适用于增材制造。
受限于高温热烧蚀，转盘雾化技术目前主要应用于制备低温锡粉和功能材料等，转盘雾化制备铝合金及更高温度合金粉末的相关实验研究仍在开展之中。20世纪90年代，佐藤富雄等[5-7]采用钼盘片制备了铝镁锆合金粉末，后又进一步分析了粉末高温酸化特性[8]。Angers等研究了转盘表面涂层技术[9]并采用倒置转盘雾化技术制备了2024铝合金粉末[10]，研究了AZ91镁合金在转盘雾化中产生的不规则颗粒[11]。Labrecque等[12]也对倒置转盘离心雾化制备镁合金技术进行了研究。陈振华等[13-14]提出了金属熔体多级快速凝固过程成粉理论并制取了铝锌合金粉末。温树德[15]基于转盘雾化技术研究了铅和铝雾化，并对8640钢进行了试探性雾化。Öztürk等[16-18]基于带鳍片转盘，采用喷嘴冷却转盘方法研究了锡、铝等金属的雾化技术。Eslamian等[19]在采用转盘离心雾化6061铝合金的过程中添加碳化硅粉末，形成了金属基复合粉末。刘英杰等[20]采用转盘雾化法制备了AlSi10Mg铝合金粉末，并对粉末各项性能和3D打印结构件性能进行了测试。Plookphol等[21]分析了转盘离心雾化制备SAC305粉末的影响参数。Zhang和Zhao[22-25]研究了离心雾化转杯几何形状和关键工艺参数对锡粉粒径和分布的影响并对金属液体的流动发展进行了分析。Alejandra等[26]研究了向雾化颗粒吹气对离心雾化制备铜粉的影响。王东祥等[27-30]分析了转盘雾化的若干流体力学问题，并采用熔融铝进行了离心雾化高温相似试验研究。从现有研究成果看，转盘雾化制备铝合金粉末技术目前仍处于实验室研究阶段，以制备少量样品和分析粉末性能为主，转盘表面金属液体的流体动力学问题研究和系统的粉末粒径分布规律研究较少。
本文基于独立开发的一套高温转盘雾化制粉实验装置，针对1060纯铝和AlSi10Mg铝合金，研究金属液体的铺展流动特征，分析流量、转速、盘径和转盘表面构型等对金属粉末性能的影响规律，基于线性回归理论分析得到粒径公式，以期为铝合金粉末的转盘离心雾化制备提供理论和技术支撑。
1 实验材料和方法
1.1 实验材料
采用1060纯铝和AlSi10Mg铝合金作为实验材料，化学成分见表1和2，表3和4分别列出了两者的物理性质（表中温度为实验中熔融铝液温度）。
表 1 1060纯铝的化学成分（%）Table 1 Chemical composition of 1060 pure aluminum (%)Al Fe Si Cu Zn V Mn Mg Ti 99.6 0.35 0.25 0.05 0.05 0.05 0.03 0.03 0.03 表 2 AlSi10Mg铝合金的化学成分（%）Table 2 Physical properties of 1060 pure aluminum (%)Al Si Zn Mg Fe Mn Ti Cu Ni Sn Pb else 9.70 0.60 0.50 0.50 0.40 0.10 0.07 0.05 0.05 0.05 注：“else”表示铝合金中除其他化学成分外，其余皆为Al的含量。 表 3 1060纯铝的物理性质Table 3 Physical properties of 1060 pure aluminumTemperature
/KDensity
/（g·cm−3）Viscosity
/（Pa·s）Surface tension
coefficient /（N·m−1）1150 2.310 0.926 × 10−3 0.843 Temperature
/KDensity
/（g·cm−3）Viscosity
/（Pa·s）Surface tension
coefficient/（N·m−1）1150 2719 1.3 × 10−3 0.826 1.2 高温转盘雾化制粉实验装置
本文基于现有研究积累，开发了一套高温转盘雾化制粉实验装置（图1和2）。该装置主要由加热炉、坩埚（内有塞杆）、转盘、高速电机、雾化室、真空泵、循环水冷系统（包括冷却水池、水泵、冷却水管道等）、氮气充装及循环冷却系统（含氮气钢瓶、风机、冷却风嘴及相关管道等）、集粉罐、高清摄像机、计算机等组成。加热炉为电阻式，电压380 V，最高温度1200 ℃，温度误差 ± 1 ℃。坩埚、塞杆都为石墨材质，坩埚容积8.8 L，塞杆头部为锥面，与坩埚下底面内锥孔配合密封。雾化室直径2.5 m，高度3.5 m，上半部为圆柱体，下半部为80°倒锥体。高速电机为标准电主轴，采用变频器控制转速，最高转速6 × 104 r/min。雾化室壁面为夹层冷却结构，以循环冷却水增压冷却。水泵为不锈钢多级离心泵（上海雨洲泵业制造有限公司，32CDLF4–160型），流量4 m3/h。
![]()
图 1 高温转盘雾化制粉实验装置示意图Fig. 1 Schematic diagram of high-temperature rotary disc atomi-zation powder-making experimental device为防止粉末氧化并提高球形度，雾化室内部为惰性气体环境：先以真空泵将雾化室抽真空，再充入高纯氮气至微正压。真空泵为旋片式（上海万经泵业制造有限公司，2X–100型），极限压力6 × 10−2 Pa。采用标准40 L钢瓶供应氮气（体积分数99.9%），减压后充入雾化室并保持微正压100 Pa，防止外部空气进入。为快速冷却飞离转盘边缘的金属液滴，同时防止转盘上表面热量经转轴向下传递至高速电机内部，设计了一路氮气循环冷却系统，对转盘和转轴周围液滴进行冷却。冷却氮气采用高压风机（浙江亚士霸电机有限公司，HG–3000型，最大风压38 kPa，流量318 m3/h）增压。雾化室侧壁开了2个标准视窗（视窗直径约90 mm），视窗外安装2台高清摄像机，视频信号通过网线传输至计算机进行实时显示和存储。雾化室下部安装了收集粉末的集粉罐。
![]()
图 2 高温转盘雾化制粉实验装置照片Fig. 2 Photograph of high temperature disc atomization powder making experimental device1.3 不同构型的转盘设计
转盘是雾化制粉装置的核心部件，需承受高温金属液体的持续冲刷，应具有耐高温、耐熔蚀的特性；为保证金属液体与转盘充分接触，转盘上表面应具有良好的润湿性；为防止金属液体的热量沿转轴向下传递至高速电机内部，转盘还应具有低热导率。本文设计了平面、锥面和弧面等3种构型的转盘（图3），用以研究不同构型转盘的雾化性能。为保证润湿性及耐高温、耐熔蚀特性，转盘以耐高温、耐熔蚀的复合材料制作，并进行了提高润湿性的加工处理。转盘具体参数见后文实验结果介绍。
1.4 转盘雾化实验
实验中，将雾化室抽真空至压力100 Pa以下，然后充入体积分数99.9%的氮气至微正压100 Pa（表压）。向坩埚中投入一块铝锭（直径80 mm，长150 mm），加热炉升温至880 ℃，铝锭完全熔化为铝液。启动高速电机至设定转速，拔出塞杆，铝液垂直流下至转盘中心，迅速沿盘面铺展流动至转盘边缘，高速甩出雾化为液滴。液滴在表面张力作用下收缩为球形并在飞行过程中与氮气快速对流换热，凝固为粉末，最终下落汇集至集粉罐中。图4为转盘雾化制备球形铝粉实验的典型照片。从图中可以看到：亮红色的高温铝液在铺展过程中与转盘发生热传递，使得转盘升温变红；转盘周围布满大量凝固后飘落的铝粉。
![]()
图 4 转盘雾化制备球形铝粉的实验图片Fig. 4 Picture of spherical aluminum powder prepared by centrifugal atomization of rotating disc1.5 粉末取样及筛分处理
受水跃或沉积层的影响，部分铝液并非均匀铺展于盘面，而是斜向飞溅而出，形成大片铝渣或铝丝。研究者普遍认为铝渣或铝丝不属于离心雾化产物，需以筛网将其过滤，收集粒径750 μm以下的粉末进行统计分析。采用二分器从收集的粉末中取出50 g样品进行筛分，将各筛网上的粉末收集至样品瓶并精确称重，基于筛下质量百分比计算粉末累积质量分数。
使用的分析测试仪器包括：扫描电子显微镜（Hitachi SU1510，二次电子分辨率3.0 nm）；粉末二分器（国标3 mm × 32格，浙江欣锴实验仪器设备有限公司）；机械振筛机（河南亚创机械设备有限公司，YC–200型）；标准分样筛（直径200 mm，60、100、150、300、400目，上海言锦筛网加工厂）；高精度电子秤（3000 ± 0.02 g）。
2 实验结果与讨论
2.1 转盘表面铝液铺展流动形貌表征
图5为雾化实验结束后转盘上表面的照片（实验转速分别为18000、24000、30000 r/min）。以18000 r/min转速的转盘为例（图6）：在直径25 mm范围内，铝膜铺覆稀疏且非常薄，命名为“光洁区（clean area）”；在直径25～44 mm范围内，转盘表面铺覆了薄薄的铝层，手触有磨砂感，命名为“磨砂区（frosted area）”。在直径44 ～48 mm范围内，磨砂层上表面又覆了一层非常薄的铝膜，命名为“铺展区（spreading area）”；“铺展区”之外（直径大于48 mm）即为常见的“沉积层区（deposit layer area）”。
从图5和6可以看到：越靠近转盘边缘，沉积层越厚（可达1～2 mm）；在部分实验中，甚至还出现了沉积层上又覆一层沉积层的“二次沉积”现象；随着转速升高，“磨砂区”范围扩大，在部分实验中，沉积层出现了螺旋棱或径向棱，相关研究认为该现象可能是由转盘和氮气的对流换热所造成的：转盘边缘处的金属液膜较薄，加热量小于冷却量，液膜更容易降温凝固并累积。
沉积层容易造成离心雾化过程中出现“水跃”，产生较大的铝渣，降低粉末收集率，因此实验中应尽量减少沉积层。在转盘雾化实验中，由于铝液流动不够充分且在转盘停止后易迅速凝固为尺寸较大的铝渣，难以观察转盘表面各典型区域，导致相关研究较少，因此，转盘表面各典型区域的产生原因还亟待开展进一步的研究分析。
图7为转速30000 r/min的实验结束后的转盘照片。从图中可以看到转盘边缘存在许多“锯齿”（约420个），雾化液滴就是在这些锯齿点位上形成。该实验现象与佐藤富雄等[6]的研究结果一致。
2.2 粉末微观形貌表征
选取粒径在15～53 μm之间的1060纯铝和AlSi10Mg合金铝粉末分析微观形貌。图8为粉末的扫描电镜图片。图9[20]给出了典型的气雾化制备的AlSi10Mg合金铝粉末形貌进行对比。从图中可以看出，在同样的惰性气体环境下，转盘雾化粉末的球形度好、表面光洁、大小均匀、无卫星粉，而气雾化粉末中充满了大小不一的卫星粉。造成这种差异的主要原因在于：在气雾化制粉过程中，液体被气流破碎为大小不一的液滴，在湍流和旋涡作用下，发生多次碰撞和黏附，小颗粒和大颗粒液滴团聚为卫星粉；而在转盘雾化制粉过程中，液体从转盘中心向边缘铺展、分散流动，在转盘边缘雾化为液滴后，沿转盘切线方向飞出，继续在雾化室内分散运动，而雾化室内氮气比较稳定，液滴在飞行过程中不易碰撞和团聚，避免了卫星粉的出现。卫星粉会导致粉末流动性变差，在增材制造中产生孔隙缺陷。因此，转盘雾化制备的粉末球形度和流动性更好，非常适合增材制造。
![]()
图 8 粒径15～53 μm的粉末扫描电镜图片Fig. 8 Scanning electron microscope images of powders between 15-53 μm![]()
2.3 粉末粒径分布分析
从某次典型实验（转速39900 r/min，转盘直径59 mm）制备的1060纯铝粉中取样50 g进行数据分析，得到该样品的粒径分布曲线，如图10所示（左侧纵坐标为差分分布，右侧为累积分布）。为表征粉末的粒径集中度，采用单峰Extreme函数对差分曲线进行拟合分析，函数表达式如下：
![]()
图 10 粉末粒径差分分布、累积分布及Extreme函数拟合Fig. 10 Particle size differential distribution, cumulative distribution, and Extreme fit$$ \begin{split}& y = {y_0} + A{e^{( - {e^{ - z}} - z + 1)}} \\& z = (x - {x_c})/w \end{split} $$ (1) 式中：x、y分别为横、纵坐标；z为计算中间值；xc为中心点横坐标；y0为偏移量；A为幅值；w为宽度。
从图中可以看到，实验数据峰值出现在105.7 μm处。Extreme函数拟合曲线的y0 = 0.08343，A = 0.44052，中心xc = 102.18724，宽度w = 9.27，说明该样品粒径分布非常窄。0～150 μm累积分布为83.04%，对应增材制造中使用的15～53 μm（选择性激光熔化SLM，Selective Laser Melting）和53～150 μm（激光近净成型LENS，Laser Engineered Net Shaping）这2个区间，说明转盘雾化制备的粉末粒径分布非常符合增材制造的要求。
2.4 分裂模式对粒径分布的影响
在转盘离心雾化过程中，由于液体流速和转盘转速不同，液体存在3种分裂模式，即直接滴状分裂（DDF）、丝状分裂（LD）和膜状分裂（FD）[32]，如图11所示。
随着液体流速增大，分裂模式逐渐由DDF模式向LD模式和FD模式变化。Champagne等[33-34]发现，对于各种金属（如Al、Cu和Zn），DDF模式下液滴的最终形状为球形，而LF模式下更易形成椭球形状。从DDF模式到LF模式和LF模式到FD模式的转变流量公式[34]为：
$$ {Q_t} = \frac{{K{\sigma ^{0.88}}{D^{0.68}}}}{{{\mu ^{0.17}}{\omega ^{0.6}}{\rho ^{0.71}}}} $$ (2) 式中：D为转盘直径，单位m；ω为转速，单位rad/s；ρ为密度，单位kg/m3；σ为表面张力，单位N/m；μ为动力学黏性系数，单位Pa·s；K为常数（从DDF向LD转变时，K = 0.07；从LD向FD转变时，K = 1.33）。从式（2）可以看出：当转速ω增大、转盘直径D减小时，转变流量Qt降低，有利于从DDF模式过渡到LF模式，并最后过渡到FD模式。
为研究不同分裂模式对雾化效果的影响，采用不同孔径的坩埚导液管来控制铝液流量。图12为相同条件下（转速18000 r/min，转盘直径59 mm，转盘为弧面盘（弧面半径90 mm），加热温度880 ℃）、金属流量Q分别为5.33和1.23 mL/s的1060纯铝粉的粒径差分分布和累积分布。从粒径分布曲线（红色和黑色曲线）可以看到：与大流量（5.33 ml/）粉末相比，小流量（1.23 mL/s）粉末的粒径分布曲线更窄，即小流量粉末粒径集中度更高，细粉末含量也更高。从粒径累积分布曲线（蓝色和绿色曲线）可以看到：小流量下0～150 μm粉末的质量分数为63.3%，粉末中位径d50 = 137.1 μm；大流量的0～150 μm粉末含量为51.52 %，中位径d50 = 151.3 μm。
![]()
图 12 不同流量下的粉末粒径分布和累积分布Fig. 12 Powder particle size distribution and cumulative percentage at different flow rates结合表3中1060纯铝的物理性质，根据式（2）可以计算得到不同转速下DDF模式向LD模式转变的临界流量Q1和LD模式向DD模式转变的临界流量Q2，如表5所示。从表中可以看到，对于上文转速18000 r/min、转盘直径59 mm的实验，Q1 = 1.28 mL/s。实验中的小流量1.23 mL/s低于临界流量Q1，属于DDF模式；而大流量5.33 mL/s处于Q1和Q2之间，属于LD模式。
表 5 离心分裂模式转变临界流量Table 5 Centrifugal atomization critical flow rate at different speedsω/（r·min−1） D/mm Q1 /（mL·s−1） Q2 /（mL·s−1） 6000 59 2.47 46.93 12000 59 1.63 30.97 18000 59 1.28 24.28 24000 59 1.08 20.43 本节实验得到的数据表明，DDF模式下的液滴中位径小于LD模式（即小流量下的中位径更小），粒径分布也更窄，这与Champagne的分裂模式理论一致。因此，在其他条件相同时，采用较小流量使分裂模式变为DDF模式，能够获得更小的粉末中位径和更高的细粉率。
2.5 中位径影响参数规律
转盘雾化粉末的中位径与液体密度、表面张力、转速、转盘直径等参数相关。Walton等[35]给出了计算离心雾化粉末中位径d50的公式：
$$ {d_{50}} = 3.8{\left( {\frac{\sigma }{{\rho D{\omega ^2}}}} \right)^{0.5}} $$ (3) Labrecque等[12]结合体积流量的影响，给出了AZ91镁合金离心雾化粉末中位径的理论公式：
$$ {d_{50}} = 4.27 \times {10^6} \cdot \frac{1}{{{\omega ^{0.95}}}} \cdot \frac{1}{{{D^{0.61}}}} \cdot {\left( {\frac{\sigma }{\rho }} \right)^{0.42}} \cdot {Q^{0.12}} $$ (4) Champagne等[36]提出了预测离心雾化粉末中位径的公式：
$$ {d_{50}} = \frac{{3.65{Q^{0.06}}}}{{\omega {D^{0.58}}}}{\left( {\frac{\sigma }{\rho }} \right)^{0.46}} $$ (5) Dunkley等[37]给出了球形液滴中位径计算公式：
$$ {d_{50}} = {\left( {\frac{{6\sigma }}{{\rho {\omega ^2}R}}} \right)^{0.5}} $$ (6) Tornberg[38]提出了一个中位径理论公式：
$$ {d_{50}} = \sqrt {\frac{{4.8\sigma }}{{\rho {\omega ^2}R\left[ {1 - {{\left( {\omega R} \right)}^{0.93}}\mu } \right]}}} $$ (7) 国为民等[39]认为，熔化的薄层液态金属雾化为极小的液滴瞬间抛出时，所受的离心力理论上与表面张力相等。离心雾化粉末中位径的理论公式[39]为：
$$ {d_{50}} = \sqrt {\frac{{12\sigma }}{{\rho {\omega ^2}D}}} $$ (8) 以上各式中：R为转盘半径，单位m；Q为金属液体流量，单位m3/s。
以上理论模型虽然不尽相同，但液滴中位径的变化规律基本一致，在液体物性参数固定的情况下，都是随着转速和转盘直径的增大而减小；对于不同的液体，密度越大、表面张力越小，中位径越小。
2.5.1 不同转速对中位径的影响
在其他参数固定（加热温度880 ℃，转盘直径60 mm，弧面半径90 mm，金属流量5.33 mL/s）的情况下，分别进行了转速为6000、12000、18000、24000 r/min的离心雾化实验，4个转速对应的中位径d50分别为266.8、179.0、147.7、116.0 μm。图13为实验值与理论公式（3）～（8）的计算值对比。从图中可以看到，各曲线的整体趋势都是中位径随转速增大而减小：当转速在10000 r/min以下时，理论计算值明显大于实验值，Champagne公式的计算值更接近实验值；当转速在10000 r/min以上时，实验值位于各计算值之间，国为民（Guo）和Labrecque公式的计算值与实验值吻合较好。
![]()
图 13 不同转速下的粉末中位径实验值与计算值对比Fig. 13 Comparison of experimental data and theoretical values of d50 at different rotational speeds实验刚开始时，转盘初始温度较低，铝液会迅速在转盘表面凝结，使转盘动平衡变差，若此时转盘转速较高，极易导致转盘剧烈振动甚至损坏。为保护转盘和电极，在实验中采用了逐渐增大转盘转速直至实验所需转速的方式。因此，高转速实验条件下制备的粉末中会掺混一部分较低转速下产生的相对较粗的粉末，降低了高转速实验条件下的细粉率，导致通过实验得到的d50值偏大，而实际高转速下的d50相对较小，与理论计算值更为接近。针对这一问题，需要进一步改进实验方法和分析手段来进行验证。
2.5.2 不同转盘直径对粉末中位径的影响
在其他参数固定（加热温度880℃，转盘弧面半径90 mm，金属流量2.94 mL/s，转速24000 r/min）的情况下，针对直径40、50、60 mm的转盘进行实验，将实验值与理论计算值进行对比，结果如图14所示。可以看到，粉末中位径的实验值和计算值都呈现出随转盘直径增大而减小的规律：直径为40 mm时，中位径为123.4 μm，与Tornberg公式的计算值最为接近；直径为50和60 mm时，中位径分别为96.4和86.5 μm，与Champagne公式的计算值最为接近。
![]()
图 14 不同转盘直径下的粉末中位径实验值与计算值对比Fig. 14 Experimental data and theoretical values of median diameter of powders with varying disc diameters2.5.3 不同转盘构型对粉体中位径的影响
在其他参数固定（加热炉温度880 ℃，转盘直径60 mm，转速24000 r/min，流量2.94 mL/s）的情况下，基于图3中的3种转盘构型，分别针对平面盘、12°锥面盘、弧面盘A（弧面半径90 mm）和弧面盘B（弧面半径80 mm）进行实验。对实验得到的粉末进行取样和筛分，分别统计d10、d50和d90（粒径累计分布为10%、50%、90%时对应的粒径），结果如图15所示。从图中可以看到：在4种转盘构型下，d10都约为54.0 μm，而d50、d90的差异较大；弧面盘B的d50最小（85.5 μm），平面盘的d50最大（99.8 μm）；12°锥面盘的d90最小（197.8 μm），弧面盘A的d90最大（257.6 μm）。因此，若需要粒径非常小的粉末，可以采用d50最小的弧面盘B（弧面半径80 mm）；若需绝大部分粉末都较细，可以采用d90最小的12°锥面盘。
![]()
图 15 不同转盘构型下的粉末粒径统计Fig. 15 Particle size statistics for different disc configurations2.5.4 中位径线性拟合公式
通过量纲分析可知，粉末中位径d50与雾化参数[σ/（ρDω2）]0.5成线性关系。图16给出了1060纯铝粉的中位径d50实验值随雾化参数[σ/（ρDω2）]0.5的变化情况及线性拟合曲线。中位径线性回归拟合公式为：
![]()
图 16 1060纯铝粉中位径实验值及线性拟合曲线Fig. 16 Data of median diameter of 1060 pure aluminum with atomization parameters and linear fitting curve$$ {d_{50}} = 3.22{\left( {\frac{\sigma }{{\rho D{\omega ^2}}}} \right)^{0.5}} $$ (9) 拟合曲线的截距为0，拟合斜率为3.22，标准误差为0.102，表明中位径与雾化参数线性关系良好。
3 结 论
本文开展了转盘雾化制备铝合金粉末的实验研究，得到如下结论：
1）铝液在转盘表面铺展流动的过程中，存在“光洁区”“磨砂区”“铺展区”和“沉积层区”等4个典型区域。
2）液体流量降低，会使转盘雾化分裂模式由丝状分裂向直接滴状分裂转变，能够降低粉末中位径并提高细粉率。
3）在本文实验的平面、锥面和弧面转盘中，弧面盘的d50最小，锥面盘的d90最小。
4）1060纯铝粉的中位径d50与雾化参数(σ/ρDω2)0.5的线性关系良好，通过线性回归拟合，提出了一个新的中位径理论公式。
致谢：本文实验操作得到了技术人员樊志勇的帮助，扫描电镜实验得到了雷现奇老师的帮助，特此致谢。
-
![]()
图 1 石墨烯热声激励器实物图与热声效应原理图
Fig. 1 Picture of graphene thermo-acoustic actuator and schematic diagram of thermoacoustic effect
![]()
图 2 并联石墨烯热声激励器示意图、实物图和测量系统示意图
Fig. 2 Schematic diagram, physical diagram and measurement system diagram of parallel graphene thermoacoustic exciter
![]()
图 4 输入频率3 kHz、测试距离1 cm下不同结构的并联石墨烯热声激励器输出声压理论值与测量值对比
Fig. 4 Comparison of the theoretical and measured output sound pressure values of parallel graphene thermo-acoustic actuator with different structures at fin=3 kHz, r=1 cm
![]()
图 5 并联石墨烯热声激励器结构及其对输出声压和幅频特性的影响
Fig. 5 The structure of the parallel graphene thermo-acoustic actuator and its influence on the output sound pressure and amplitude-frequency characteristic
![]()
图 6 输入电频率3 kHz、输入电压25 V下，并联石墨烯热声激励器结构对输出声压和幅频特性的影响
Fig. 6 The influence of the parallel graphene thermo-acoustic actuator’s structure on the output sound pressure and amplitude-frequency characteristic (fin=3 kHz, Uin=25 V)
-
[1] 中国科学院. 中国学科发展战略-新型飞行器中的关键力学问题[M]. 北京: 科学出版社, 2018. [2] LIN J C. Review of research on low-profile vortex generators to control boundary-layer separation[J]. Progress in Aerospace Sciences, 2002, 38(4-5): 389–420. doi: 10.1016/S0376-0421(02)00010-6
[3] JIRASEK A. Vortex-generator model and its application to flow control[J]. Journal of Aircraft, 2005, 42(6): 1486–1491. doi: 10.2514/1.12220
[4] ANUAR M M, NORDIN N, SHARIFF M Z F, et al. Effect of Employing Vortex Generator on Curve Diffuser Performance[J]. Journal of Advanced Mechanical Engineering Applications, 2021, 2(2): 53–60.
[5] WANG J J, LI Y C, CHOI K S. Gurney flap—Lift enhancement, mechanisms and applications[J]. Progress in Aerospace Sciences, 2008, 44(1): 22–47. doi: 10.1016/j.paerosci.2007.10.001
[6] MOUSAVI M, MASDARI M, TAHANI M. Power performance enhancement of vertical axis wind turbines by a novel gurney flap design[J]. Aircraft Engineering and Aerospace Technology, 2021. doi: 10.1108/aeat-02-2021-0052
[7] KUMAR V, ALVI F. Efficient control of separation using microjets[C]//Proc of the 35th AIAA Fluid Dynamics Conference and Exhibit. 2005: 4879. doi: 10.2514/6.2005-4879
[8] 王林, 罗振兵, 夏智勋, 等. 高速流场主动流动控制激励器研究进展[J]. 中国科学(技术科学), 2012, 42(10): 1103–1119. DOI: 10.1360/ze2012-42-10-1103 WANG L, LUO Z B, XIA Z X, et al. Review of actuators for high speed active flow control[J]. SCIENTIA SINICA Technologica, 2012, 42(10): 1103–1119. doi: 10.1360/ze2012-42-10-1103
[9] MOREAU E. Airflow control by non-thermal plasma actuators[J]. Journal of Physics D: Applied Physics, 2007, 40(3): 605–636. doi: 10.1088/0022-3727/40/3/s01
[10] ZHU Z H, FRADERA-SOLER P, JO W, et al. Numerical simulation of the flow around a square cylinder under plasma actuator control[J]. Physics of Fluids, 2021, 33(12): 123611. doi: 10.1063/5.0072081
[11] 罗振兵, 夏智勋. 合成射流技术及其在流动控制中应用的进展[J]. 力学进展, 2005, 35(2): 221–234. DOI: 10.3321/j.issn:1000-0992.2005.02.009 LUO Z B, XIA Z X. Advances in synthetic jet technology and applications in flow control[J]. Advances in Mechanics, 2005, 35(2): 221–234. doi: 10.3321/j.issn:1000-0992.2005.02.009
[12] PALUMBO A, DE LUCA L. Experimental and CFD characterization of a double-orifice synthetic jet actuator for flow control[J]. Actuators, 2021, 10(12): 326. doi: 10.3390/act10120326
[13] GREGORY J W, GNANAMANICKAM E P, SULLIVAN J P, et al. Variable-frequency fluidic oscillator driven by a piezoelectric bender[J]. AIAA Journal, 2009, 47(11): 2717–2725. doi: 10.2514/1.44078
[14] KARA K, KIM D, MORRIS P J. Flow-separation control using sweeping jet actuator[J]. AIAA Journal, 2018, 56(11): 4604–4613. doi: 10.2514/1.J056715
[15] ZAMAN K B M Q, BAR-SEVER A, MANGALAM S M. Effect of acoustic excitation on the flow over a low- Re airfoil[J]. Journal of Fluid Mechanics, 1987, 182: 127. doi: 10.1017/s0022112087002271
[16] NISHIOKA M, ASAI M, YOSHIDA S. Control of flow separation by acoustic excitation[J]. AIAA Journal, 1990, 28(11): 1909–1915. doi: 10.2514/3.10498
[17] KURELEK J W, KOTSONIS M, YARUSEVYCH S. Transition in a separation bubble under tonal and broadband acoustic excitation[J]. Journal of Fluid Mechanics, 2018, 853: 1–36. doi: 10.1017/jfm.2018.546
[18] ANDAN A D, LEE D J. Effect of external acoustic excitation on NACA0015 discrete tonal noise[J]. Applied Acoustics, 2018, 141: 374–381. doi: 10.1016/j.apacoust.2018.07.030
[19] BERNARDINI C, BENTON S I, BONS J P. Understanding leading edge stall physics by acoustic excitation[C]//Proc of the 52nd Aerospace Sciences Meeting. 2014: 1125. doi: 10.2514/6.2014-1125
[20] 朱奇亮, 高永卫, 叶正寅. 不同声激励方式对多段翼型升力特性的影响[J]. 航空工程进展, 2013, 4(2): 232–236. DOI: 10.3969/j.issn.1674-8190.2013.02.016 ZHU Q L, GAO Y W, YE Z Y. Influence of the different acoustic excitations on lift characteristic of a multi-element airfoil[J]. Advances in Aeronautical Science and Engineering, 2013, 4(2): 232–236. doi: 10.3969/j.issn.1674-8190.2013.02.016
[21] WANG Y J, HE Z, LIU D W. Effect of acoustic excitation on flow asymmetry over slender body at high angles of attack[J]. Procedia Engineering, 2012, 31: 182–186. doi: 10.1016/j.proeng.2012.01.1010
[22] TIAN H, XIE D, YANG Y, et al. Static behavior of a graphene-based sound-emitting device[J]. Nanoscale, 2012, 4(11): 3345–3349. doi: 10.1039/c2nr30417a
[23] TIAN H, REN T L, XIE D, et al. Graphene-on-paper sound source devices[J]. ACS Nano, 2011, 5(6): 4878–4885. doi: 10.1021/nn2009535
[24] 卞安华. 石墨烯薄膜声学特性及主动控制研究[D]. 苏州: 苏州大学, 2017. BIAN A H. Study on acoustic characteristics of graphene film and its application in active noise control[D]. Suzhou: Soochow University, 2017.
[25] WANG Z G, CHEN Y F, LI P J, et al. Flexible graphene-based electroluminescent devices[J]. ACS Nano, 2011, 5(9): 7149–7154. doi: 10.1021/nn2018649
[26] 邢倩荷. 石墨烯薄膜扬声器热声特性分析及应用[D]. 苏州: 苏州大学, 2018. XING Q H. Analysis of thermoacoustic characteristics of graphene film loudspeakers and its application[D]. Suzhou: Soochow University, 2018.
[27] DASCHEWSKI M, BOEHM R, PRAGER J, et al. Physics of thermo-acoustic sound generation[J]. Journal of Applied Physics, 2013, 114(11): 114903. doi: 10.1063/1.4821121
[28] 杜功焕, 朱哲民, 龚秀芬. 声学基础[M]. 3版. 南京: 南京大学出版社, 2012. [29] HU H P, ZHU T, XU J. Model for thermoacoustic emission from solids[J]. Applied Physics Letters, 2010, 96(21): 214101. doi: 10.1063/1.3435429
[30] 黄伟. 基于热声效应的石墨烯扬声器研究[D]. 杭州: 浙江大学, 2020. HUANG W. Research on the graphene loudspeakers based on thermoacoustic effect[D]. Hangzhou: Zhejiang University, 2020.
[31] SUK J W, KIRK K, HAO Y F, et al. Thermoacoustic sound generation from monolayer graphene for transparent and flexible sound sources[J]. Advanced Materials, 2012, 24(47): 6342–6347. doi: 10.1002/adma.201201782
-
期刊类型引用(1)
1. 蒋丹丹，张浩然. 试飞中的扰流板安全分析研究. 科技与创新. 2024(05): 163-165 . 
百度学术
其他类型引用(2)
下载:
计量
- 文章访问数: 201
- HTML全文浏览量: 91
- PDF下载量: 19
- 被引次数: 3
