留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

加速器用铁氧体-碳化硅混合型高次模阻尼器的研制

陈欣 李晨 赵伟 黄刚 向军 李天涛 杨洁 刘平 秦臻

陈金宝, 相广彪, 王小林, 等. 低辐照剂量率下光纤激光器输出功率的自漂白-辐照平衡实验研究[J]. 强激光与粒子束, 2024, 36: 121001. doi: 10.11884/HPLPB202436.240384
引用本文: 陈欣, 李晨, 赵伟, 等. 加速器用铁氧体-碳化硅混合型高次模阻尼器的研制[J]. 强激光与粒子束, 2025, 37: 024001. doi: 10.11884/HPLPB202537.240154
Chen Jinbao, Xiang Guangbiao, Wang Xiaolin, et al. Experimental study of self-bleaching and radiation equilibrium in output power of fiber lasers at low dose rates[J]. High Power Laser and Particle Beams, 2024, 36: 121001. doi: 10.11884/HPLPB202436.240384
Citation: Chen Xin, Li Chen, Zhao Wei, et al. Development of a ferrite-silicon carbide hybrid high-order mode damper for accelerators[J]. High Power Laser and Particle Beams, 2025, 37: 024001. doi: 10.11884/HPLPB202537.240154

加速器用铁氧体-碳化硅混合型高次模阻尼器的研制

doi: 10.11884/HPLPB202537.240154
基金项目: 国家自然科学青年基金项目(12205280)
详细信息
    作者简介:

    陈 欣,chenxin10698@126.com

    通讯作者:

    李 晨,st.eddie@163.com

  • 中图分类号: TL503.2

Development of a ferrite-silicon carbide hybrid high-order mode damper for accelerators

  • 摘要: 大电流加速器束管中,当带电粒子流通过束管时,会在束管中激励起高频场,为了降低对束流的影响,束管中产生的高次模需要利用阻尼器将高频场能量转换成热量并通过冷却装置导走。介绍了某混合型高次模阻尼器的研制及主要性能。阻尼器采用的吸收材料为铁氧体和碳化硅,吸收材料通过金属化和钎焊实现与金属基板的焊接。通过CST和COMSOL软件分别开展了微波性能仿真和热仿真,对阻尼器的结构进行了优化设计。阻尼器的测试结果表明:该混合型阻尼器的吸收效率与计算结果在1.7 GHz以下频段相接近,在1.7 GHz以上高频段后,仿真吸收效率高于实测结果,相差较大;真空漏率、极限真空、水路耐压均满足超导高频腔设计需求。
  • 在惯性约束聚变(ICF)实验中,中子和包括质子与α粒子等在内的带电粒子是重要的核聚变反应产物[1]。质子产额与质子能谱能够用于诊断燃料-壳层面密度、热斑离子温度、内爆压缩不对称性、燃料-壳层混合效应等重要物理信息[2-4]。目前ICF质子诊断通常使用CR-39固体径迹探测器作为记录介质[5]。尽管CR-39同时具有较好的能量和空间分辨率以及对X射线和电子不灵敏等优势,但是照射后的CR-39探测器需经过一系列离线化学蚀刻处理才能获得质子能谱,并且其质子能谱会受蚀刻的温度和时间影响。这种离线获取能谱的方式,从照射到读取信息通常会有一天以上的间隔,影响了测量的时效性。而且,固体径迹探测器具有较大个体差异,进一步增加了能谱提取过程的不确定性。因此,在进行ICF聚变等离子体相关物理过程的精细研究时,提高测量的准确性与时效性是十分有必要的。Timepix探测器的在线数据处理能力能提供带电粒子类型、能谱、线性能量传递(LET)和方向追踪等方面的信息,提供了提高质子能谱测量时效性与一致性的解决方案[6-7]。D-3He聚变反应产生的质子初始能量为14.7 MeV。实际在以D-3He为燃料的内爆实验中,其质子能谱会被改造,通常为12 MeV与15 MeV[8]。此外,除聚变产物外,ICF内爆实验过程还会产生爆炸碎片[1],这些碎片可能会导致探测器损坏。因此,ICF的诊断系统通常会使用一定厚度的滤片对探测器进行防护。然而质子在通过防护滤片后会有能量损失,且由于能量损失歧离和射程歧离现象,会导致质子能量有所下降以及质子能谱的展宽,同时也会偏离入射方向[9]。单能质子束在通过滤片到达探测器时的能量与初始能量并不一致,在实际实验中也并不能确保质子束始终垂直入射探测器。因此,开展关于Timepix探测器对不同入射角度质子能量响应的研究,对于实现将Timepix探测器用于ICF实验中质子能谱的探测有非常重要的意义。本文采用蒙特卡罗模拟方法,依托于Allpix2硅像素模拟框架开展Timepix探测器对不同能量和不同角度质子响应的研究。在此基础上,给出了利用Timepix探测器在线获取质子能谱的最佳能段,利用适当的滤片使到达探测器的质子能量降低至最佳能段内有利于通过Timepix探测器实现聚变质子能谱的在线获取。

    Timepix探测器提供了高密度的能量敏感像素矩阵,像素间距为55 μm,具有256 × 256个像素。该探测器采用辐射敏感半导体器件通过凸点键合到像素化读出ASIC Timepix芯片。每个像素均配备有高度集成的信号电子器件(放大器、幅度甄别器和数字计数器)。这些像素可单独编程为命中计数模式、到达时间模式和超阈值时间模式[6-7, 10]图1所示为搭载Timepix芯片的Minipix探测器,其体型小巧,可利用高度集成的电源和包括USB接口的读取设备[11- 12]构成紧凑的便携式系统。其在线响应、数据读出和数据预处理功能由集成软件工具PIXet实现[7]

    图  1  搭载Timepix芯片的Minipix探测器[14]
    Figure  1.  The MiniPix detector equipped with a Timepix chip[14]

    Timepix能够高效识别多种粒子种类,包括光子(如X射线)、轻带电粒子(如电子、介子)以及重带电粒子(如质子、α粒子和重离子)。此外,Timepix可探测能量范围广,能覆盖从几MeV/u至几百MeV/u的能量区间,并且在粒子方向上有良好的分辨能力。由于Timepix能识别多种粒子以及其宽能量范围的响应特征,其被越来越多地应用于宇宙空间与粒子放疗等混合辐射场的探测和表征[13]

    Allpix2是一个用户友好的通用开源软件框架,可以用于模拟半导体像素探测器性能[15]。该框架采用模块化结构,将框架核心与实际物理模拟过程分离。如图2所示,框架核心提供了模拟中所有模块通用的四个基本组件:模块运行逻辑、探测器几何结构、信息传递系统和用户界面。Allpix2使用主配置文件与探测器配置文件来设置和执行模拟。主配置文件定义框架参数以及每个单独模块的配置,而探测器配置文件描述了模拟设置中所有探测器的位置与类型。由框架核心调用电荷沉积、粒子输运、电场转换模块完成实例模拟。再通过数字化和输出模块以实现核心与各模块间的信息传递以及输出和存储数据。用户则可以使用ROOT读取存储的信息。

    图  2  Allpix2框架的结构
    Figure  2.  Structure of the Allpix2

    为探究探测器对质子响应的最佳工作能段,在本工作中,模拟分为能量组与角度组进行,探测器与粒子类型设置为Timepix与质子。对于能量组,入射角度固定为垂直于探测平面入射,模拟能量范围在1~16 MeV的探测器响应。对于角度组,分别模拟5.5 MeV与14.7 MeV质子在不同入射角度的探测器响应。具体模拟参数如表1所示,θ为质子入射方向与探测器平面法线的夹角(图3),θ step为角度设置步长;energy为质子束能量,energy step为能量设置步长。

    表  1  模拟参数
    Table  1.  Parameters of simulation
    group detector type particle type θ/(°) θ step/(°) energy of proton/MeV energy step of proton/MeV
    energy Timepix proton 0 / 1-8 1
    0 / 10-16 2
    0 / 6.1-6.5 0.1
    0 / 6.15 /
    angle Timepix proton 0-75 15 5.5 /
    14.7
    下载: 导出CSV 
    | 显示表格
    图  3  Timepix原理示意图
    Figure  3.  Principle schematic of Timepix

    Timepix探测器的工作原理如图3所示。带电粒子(如质子)穿过传感器灵敏区域时,会在其路径上电离产生电子和空穴。这些电子和空穴在耗尽偏置电压的电场中向两极漂移,并在垂直于电场方向水平扩散。单次电离产生的扩散电荷通常会波及多个像素,形成与单个事件相关的像素化簇轨迹(图3)。通过对簇的各像素计数分布分析可得粒子的入射角度,再结合探测器厚度可以得到路径长度等信息[16]。簇可以根据形态被分为3个主要组成部分:电荷水平最高的核心区域、紧邻核心的低电荷水平区域(裙边)以及低电荷水平的直线分支区域。核心与裙边区域组成簇的主体部分[17]。除电荷水平差异外,高能质子与低能质子产生的簇在形态上也存在着显著的区别。图4给出了5.5 MeV和14.7 MeV质子垂直入射产生的簇,相比5.5 MeV的质子,14.7 MeV的质子产生的簇在主体部分之外还出现了分支(图4(b))。这些直线分支来源于δ射线,主要由反冲电子组成[17]。δ射线从主体中心开始向外延伸,能量足够高的δ射线能够穿过裙带边缘形成能被肉眼识别的分支,其能量包含在该部分低电荷水平的信息中。δ射线的分支尾部对角度计算无重要作用,通常在形态处理时被去除。

    图  4  不同能量质子束垂直入射产生的像素化簇
    Figure  4.  Pixelated cluster produced by proton beams in vertical incidence with different energy

    探测器收集到的电荷信号与质子束的能量密切相关,不同能段的质子束产生的簇具有不同的特征。当质子束能量足够高时,所有质子都能穿过传感器灵敏区域,只有部分能量沉积在传感器中,从而导致能量岐离的产生[9]。沉积在传感器灵敏区域中的质子能量在一定范围内变化,使得簇电荷分布展宽。对于能量低的质子束,质子不足以穿过传感器灵敏区域,其能量全部沉积在传感器灵敏区域中,簇电荷的分布不会有显著变化。

    图5展示了不同能量质子束在探测器中产生的簇电荷分布。通过模拟,我们发现质子束在3种不同的能量区间时,簇电荷分布表现出不同的特征:质子束能量范围在1~6 MeV时,全部质子沉积在传感器灵敏区域中,簇电荷分布形成尖峰,簇电荷随质子能量升高而增加(图5(a));质子束能量范围在6.5~16 MeV时,全部质子穿过传感器灵敏区域,电荷分布形成宽峰,簇电荷随质子能量升高而减少且分布范围缩小(图5(b));当质子束能量范围在6.1~6.2 MeV内时,部分质子穿过传感器灵敏区域的同时,另一部分质子完全沉积在其中,簇电荷分布表现为尖峰与宽峰同时存在(图5(c))。在这一特殊情况下,随着质子束能量升高,穿过传感器灵敏区域的质子占据主导,宽峰左移尖峰右移,同时伴随着尖峰计数急剧减少。

    图  5  在不同能量范围的质子束产生的簇电荷分布
    Figure  5.  Charge cluster distribution produced by proton beams in different energy ranges

    图6展示了在1~16 MeV范围内的质子束的簇电荷随能量的变化曲线,以半高宽作为误差棒。能量在1~6 MeV范围的质子,其产生的簇电荷与能量成正比关系,对应簇电荷分布的尖峰部分;能量在6.3~16 MeV范围的质子,其产生的簇电荷与能量成反比关系,对应簇电荷分布的宽峰部分;在质子能量于6.1~6.2 MeV内时,两组曲线的能量重叠,对应簇电荷分布宽峰与尖峰同时存在部分(图6(b))。值得注意的是,不同能量的质子束入射时,可能会在探测器中产生相同的簇电荷。通过对簇电荷分布的分析,可以进一步确定入射质子束的能量区间。

    图  6  簇电荷随能量的变化分布,选取半高宽为误差棒
    Figure  6.  Energy distribution of cluster charge, with the full width at half maximum (FWHM) selected as the error bars

    因此,在将Timepix探测器用于以D-3He为聚变燃料的ICF内爆实验时,配合适当的滤片使质子能量降低至6 MeV以下能区。此时簇电荷展宽较小,质子能量全部沉积在传感器灵敏区域中,簇电荷是质子能量的线性函数,探测效率高。

    像素化簇轨迹的形态和电荷分布与质子的入射角度相关,簇尺寸以像素化簇轨迹主体部分的像素数衡量。在这里对像素化簇轨迹的尺寸进行简化分析,将极角α设为0,并关注于x方向上束线入射方向与法线间的夹角即顶角θ图2)。这使得仅在x方向上的轨迹投影会随顶角θ的变化而发生改变。根据质子束能量将Timepix探测器响应随角度的变化分为两种情况:质子束能量为14.7 MeV,模拟质子能量足够穿过传感器灵敏区域的情况;质子束能量为5.5 MeV,模拟质子能量不足以穿过传感器灵敏区域并最终沉积在其中的情况。当质子束能量为14.7 MeV时,质子拥有足够能量穿过传感器灵敏区域,质子穿过传感器的轨迹的路径深度l随着顶角θ增大而增大。路径深度l与其轨迹投影x1可分别表示为

    ldsecθ
    (1)
    x1dtanθ
    (2)

    式中:d为传感器灵敏区域厚度。当质子束能量为5.5 MeV时,质子的能量不足以穿过传感器灵敏区域,质子沉积在传感器中的轨迹的路径深度l0围绕质子的平均射程涨落。则其轨迹投影x2可以表示为

    x2l0sin(θ)
    (3)

    图7(a) 展示了质子束能量在5.5 MeV以及14.7 MeV时的簇电荷随角度的变化分布。当质子束能量为14.7 MeV时,随顶角θ的增大,l也随之增大,质子损失并沉积在传感器灵敏区域中的能量也就越多,从而导致探测器收集到的电荷增多。然而当顶角θ增加到75°时,除了部分能量沉积导致的电荷宽峰外,还伴随有尖峰的产生。这说明在顶角为75°时,14.7 MeV的质子有部分会沉积在传感器灵敏区域中,无法穿过。当质子束能量仅为5.5 MeV时,质子全部沉积在传感器灵敏区域中,探测器收集到的电荷不会随顶角θ的变化而显著改变。

    图  7  簇电荷与簇尺寸随角度的变化分布
    Figure  7.  Angular distribution of cluster charge and cluster size

    图7(b) 展示了质子能量在5.5 MeV以及14.7 MeV时的簇尺寸随角度的变化分布。当质子能量为14.7 MeV时质子能穿过传感器区域,随着顶角θ逐渐增大,质子轨迹延长,电离产生电荷的区域随之增大,受其轨迹影响的像素也越多,其轨迹投影x1也随之增加。使用正切函数对簇尺寸的角度分布进行拟合,拟合函数与模拟数据符合较好,因此能够穿过传感器的质子的簇尺寸随顶角θ的变化大致为正切分布。然而,质子轨迹的路径深度l与其投影x1并不会随顶角θ的增大而无限增大。在单个Timepix探测器中,受尺寸轨迹的限制,路径深度l最大为传感器灵敏区域的对角线长度,轨迹投影x1最大为像素平面对角线长度。类似的,当质子束能量仅为5.5 MeV时质子全部沉积在传感器灵敏区域中,随着顶角θ逐渐增大,电离产生电荷的区域随之增大,受其轨迹影响的像素也越多,其轨迹投影x2也随之增加。然而,路径深度l0不会随之变化,受其限制x2的增加非常有限。使用正弦函数对沉积簇尺寸x2的角度分布进行拟合,拟合函数与模拟数据符合较好,因此无法穿过传感器的质子的簇尺寸随顶角θ的变化大致为正弦分布。

    若质子能量足够高,能够穿过传感器灵敏区域,其簇尺寸是角度的正切函数,其簇电荷随角度增加而显著增加。相反,当质子能量在6 MeV以下能区时,质子能量全部沉积在传感器灵敏区域内,其簇尺寸是角度的正弦函数,其簇电荷不会随角度改变发生显著变化。此时,质子入射角度不会对Timepix探测器的能量响应产生显著影响。

    本文采用蒙特卡罗方法,在Allpix2框架下模拟了Timepix探测器对质子能量和角度的响应。得到了Timepix探测器对能量在1~16 MeV范围质子的响应,以及5.5 MeV和14.7 MeV质子束产生的簇尺寸与簇电荷在0~75°范围的角度分布。对于无法穿过传感器灵敏区域的质子,簇电荷随质子能量增加线性增加且对角度变化不敏感,簇尺寸与角度关系符合正弦分布;对于能够穿过传感器灵敏区域的质子,簇电荷与能量呈反比关系,簇尺寸随角度的变化符合正切分布;对于可能沉积在传感器灵敏区域边缘的质子,能量响应则同时具备两种特征。但随着能量变化,两种特征之一迅速消失。根据模拟结果,在将Timepix探测器用于ICF聚变质子的在线测量时,应使到达探测器的质子能量降低到6 MeV以下,此时其探测效率最好且对角度的变化不敏感,有利于ICF实验中聚变质子能谱的获取。

  • 图  1  高次模阻尼器结构示意图

    Figure  1.  Schematic of the hybrid high-order mode damper

    图  2  混合型高次模阻尼器,内径为320 mm,高度为445 mm

    Figure  2.  Hybrid high-order mode damper with an inner diameter of 320 mm and a height of 445 mm

    图  3  混合型高次模阻尼器仿真与实测吸收效率(无短路活塞)

    Figure  3.  Simulation and measurement of absorption efficiency of hybrid high-order mode damper (no short piston)

    图  4  混合型高次模阻尼器在吸收功率为10 kW时的水路温度分布仿真结果

    Figure  4.  Simulation results of water temperature distribution for the hybrid high-order mode damper with an absorbed power of 10 kW

    图  5  混合型高次模阻尼器在吸收功率为10kW时的吸波材料温度分布仿真结果

    Figure  5.  Simulation results of absorbing materials temperature distribution for the hybrid high-order mode damper with an absorbed power of 10kW

    图  6  不同吸收功率下进出口水的温升(约0.4 MPa水压)

    Figure  6.  Temperature rise of inlet and outlet water under different absorbed power (about 0.4 MPa water pressure)

    图  7  不同输入功率下吸收功率及吸收效率

    Figure  7.  Absorption power and absorption efficiency under different input power

    表  1  混合型高次模阻尼器在吸收功率为10 kW时的温度分布仿真及测试结果

    Table  1.   Simulation and measurement results of temperature distribution for the hybrid high-order mode damper with an absorbed power of 10 kW

    temperature difference between the inlet and outlet cooling water/℃
    simulation result 2.3
    measurement result 1.1
    下载: 导出CSV

    表  2  混合型高次模阻尼器的性能测试结果

    Table  2.   Performance test results of the hybrid high-order mode damper

    vacuum leak rate/(Pa·L·s−1) ultimate vacuum/Pa water-resistant/MPa
    5×10−10 4.6×10−8 0.92
    下载: 导出CSV

    表  3  混合型高次模阻尼器的设计要求及性能测试结果

    Table  3.   Design requirements and performance test results of the hybrid high-order mode damper

    absorbed
    power/kW
    absorption
    efficiency
    vacuum leak
    rate/
    (Pa·L·s−1)
    ultimate
    vacuum/
    Pa
    water-
    resistant/
    MPa
    temperature difference
    between the outlet and
    inlet cooling water/℃
    design
    requirements
    ≥10 operating frequency band(0.6~3.0 GHz)≥30%,
    critical frequency bands(0.8~1.5 GHz)≥50%
    ≤1×10−7 ≤6.5×10−8 ≥0.9 ≤5
    simulation
    results
    10 operating frequency band(0.6~3.0 GHz)≥53%,
    critical frequency bands(0.8~1.5 GHz)≥61%
    — — — — — — 2.3
    measurement
    results
    10.2 operating frequency band(0.6~3.0 GHz)≥38%,
    critical frequency bands(0.8~1.5 GHz)≥60%
    5×10−10 4.6×10−8 0.92 1.1
    下载: 导出CSV
  • [1] 周文中, 潘卫民, 葛锐, 等. 中国散裂中子源二期双spoke超导腔设计[J]. 强激光与粒子束, 2023, 35:034004 doi: 10.11884/HPLPB202335.220266

    Zhou Wenzhong, Pan Weimin, Ge Rui, et al. Design of the China Spallation Neutron Source phase II double spoke resonator[J]. High Power Laser and Particle Beams, 2023, 35: 034004 doi: 10.11884/HPLPB202335.220266
    [2] 蒲小云, 侯洪涛, 马震宇, 等. 上海光源500MHz超导腔水平测试[J]. 强激光与粒子束, 2019, 31:115104 doi: 10.11884/HPLPB201931.190163

    Pu Xiaoyun, Hou Hongtao, Ma Zhenyu, et al. Horizontal test of 500 MHz superconducting cavity for SSRF[J]. High Power Laser and Particle Beams, 2019, 31: 115104 doi: 10.11884/HPLPB201931.190163
    [3] 米正辉, 沙鹏, 孙毅, 等. BEPCⅡ国产500 MHz超导腔运行综述[J]. 强激光与粒子束, 2018, 30:085103 doi: 10.11884/HPLPB201830.170485

    Mi Zhenghui, Sha Peng, Sun Yi, et al. Operation of domestic 500 MHz superconducting cavity for BEPCⅡ[J]. High Power Laser and Particle Beams, 2018, 30: 085103 doi: 10.11884/HPLPB201830.170485
    [4] Tajima T, Asano K, Furuya T, et al. Bonding of a microwave-absorbing ferrite, TDK IB-004 with copper for the HOM damper of the KEK B-factory SC cavities[C]//Proceedings of the 6th Workshop on RF Superconductivity. 1993: 1914-1916.
    [5] Terui S, Ishibashi T, Shirai M, et al. Development of ferrite higher order mode damper for SuperKEKB vacuum system[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2023, 1053: 168371. doi: 10.1016/j.nima.2023.168371
    [6] Tajima T, Akai K, Asano K, et al. Beam test results on HOM absorber of superconducting cavity for KEKB[C]//Proceedings of 1997 Particle Accelerator Conference. 1997: 3090-3092.
    [7] Tajima T, Asano K, Furuya T, et al. HOM absorbers of superconducting cavities for KEKB[C]//Proceedings of the 5th European Particle Accelerator Conference. 1996: 2127-2129.
    [8] Tajima T, Asano K, Furuya T, et al. Recent development of HOM absorbers for KEKB superconducting cavities[C]//Proceedings of 1997 Workshop on RF Superconductivity. 1997: 709-724.
    [9] Nishiwaki M, Akai K, Furuya T, et al. Developments of HOM dampers for SuperKEKB superconducting cavity[C]//Proceedings of SRF2013. 2014: 1058-1060.
    [10] Moffat D, Barnes P, Kirchgessner J, et al. Design and fabrication of a ferrite-lined HOM load for CESR-B[C]//Proceedings of International Conference on Particle Accelerators. 1993: 977-979.
    [11] Huang Tongming, Pan Weimin, Wang Guangwei, et al. The development of the 499.8 MHz superconducting cavity system for BEPCII[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2021, 1013: 165649. doi: 10.1016/j.nima.2021.165649
    [12] Valles N, Eichhorn R, Hoffstaetter G, et al. HOM studies of the Cornell ERL main linac cavity: HTC-1 through HTC-3[C]//Proceedings of IPAC2013. 2013: 2461-2463.
    [13] Nishiwaki M, Akai K, Furuya T, et al. Developments of SiC damper for superKEKB superconducting cavity[C]// Proceedings of SRF2015. 2015: 1289-1293.
    [14] Xu Wencan, Conway Z A, Daly E, et al. High-power test results for a cylindrical-shell silicon carbide higher-order-mode damper[J]. Physical Review Accelerators and Beams, 2024, 27: 031601. doi: 10.1103/PhysRevAccelBeams.27.031601
    [15] Hao Xuerui, Li Zhongquan, Ye Kuangkuang, et al. 500 MHz higher order mode damped cavity designed for 4th generation synchrotron radiation sources[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2022, 1040: 167273. doi: 10.1016/j.nima.2022.167273
  • 加载中
图(7) / 表(3)
计量
  • 文章访问数:  184
  • HTML全文浏览量:  78
  • PDF下载量:  26
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-05-11
  • 修回日期:  2024-11-27
  • 录用日期:  2024-11-27
  • 网络出版日期:  2024-12-07
  • 刊出日期:  2025-02-15

目录

/

返回文章
返回