空间科学

紧凑型扩束压散角束线设计及相空间演化分析

  • 王建军 ,
  • 李奕霖 ,
  • 铁维昊 ,
  • 扈晓程 ,
  • 方进勇
展开
  • 中国空间技术研究院西安分院, 西安 710100

(1985年—),男,高级工程师。主要研究方向为空间辐照效应、束流动力学。通信地址:陕西省西安市长安区东长安街 504 号空间电子信息技术研究院(710100)电子邮箱:

网络出版日期: 2026-07-22

Design of a Compact Beam-Expansion and Divergence-Reduction Beamline and Analysis of Transverse Phase-Space Evolution

  • Jianjun WANG ,
  • Yilin LI ,
  • Weihao TIE ,
  • Xiaocheng HU ,
  • Jinyong FANG
Expand
  • China Academy of Space Technology (Xi’an), Xi’an 710100, China

Online published: 2026-07-22

摘要

为满足地面电子束辐照实验对大面积、低发散角辐照场的需求,设计一种由两组Halbach四极磁铁、一组电四极磁铁和一组螺线管磁铁构成的紧凑型扩束压散角束线,并采用束流动力学仿真方法对束流包络及横向相空间的演化规律进行了分析。结果表明,束斑90%半径由0.76 mm扩展至25 mm,发散角小于0.05 mrad的粒子比例由3.60%提升至52.31%,束流扩展与发散角压缩效果较为显著。研究验证了永磁与电磁混合构型在约1.9 m紧凑尺度下协同实现扩束与发散角压缩的可行性,为受限空间条件下电子束大面积辐照平台的束流光学设计提供参考。

本文引用格式

王建军 , 李奕霖 , 铁维昊 , 扈晓程 , 方进勇 . 紧凑型扩束压散角束线设计及相空间演化分析[J]. 空间科学与试验学报, 2026 , 3(3) : 61 -67 . DOI: 10.19963/j.cnki.2097-4302.2026.03.005

Abstract

To meet the requirements of ground-based electron irradiation experiments for a large-area irradiation field with low beam divergence, a compact beamline capable of simultaneous beam expansion and divergence compression was developed. The beamline comprises two sets of Halbach quadrupole magnets, one electromagnetic quadrupole magnet, and one solenoid, spanning a total length of approximately 1.9 m. Beam dynamics simulations were carried out to investigate the evolution of the beam envelope and transverse phase space throughout the beamline. The results demonstrate that the 90% beam radius increases from 0.76 mm to 25 mm, while the fraction of particles with divergence angles below 0.05 mrad rises from 3.60% to 52.31%, indicating notable effects on beam expansion and divergence compression. These findings confirm the viability of a hybrid permanent-magnet and electromagnetic lattice for achieving coordinated beam expansion and divergence reduction within a compact footprint, providing a reference for beam-optical design of large-area electron irradiation platforms under space-constrained conditions.

0 引 言

航天器电子器件的抗辐射性能评估通常依赖地面电子辐照试验,电子加速器因能够稳定提供可控电子束,成为开展空间辐射效应地面模拟试验与器件抗辐射性能评估的重要平台[1-4]。然而,加速器引出口电子束束斑通常仅为毫米量级,直接辐照大尺寸样品易造成局部功率密度过高[5],难以满足大面积辐照需求;同时,空间电荷效应会使束流在漂移传输过程中包络扩展并伴随横向发散增强[6],限制有效传输距离,增加样品区束流调控与辐照难度[7]。受实验终端布局、屏蔽结构及系统集成空间约束,如何在紧凑尺度内协同实现束斑扩展、包络匹配与发散角控制,已成为大面积低发散电子束辐照束线设计的关键问题。
为实现大面积辐照,常用方法主要包括小束斑扫描、散射扩束和磁光学扩束。扫描磁铁可通过时变偏转磁场实现束流横向覆盖,已在多类加速器束流调控和辐照装置中得到应用[8-11],但其束流传播方向随扫描过程变化,难以获得静态、低发散的辐照场,且可能引入时间结构相关的均匀性问题。散射扩束通常依靠靶材料增加横向角散以扩大束斑,但会导致能量损失、能散增大及附加辐射增强,不利于辐照场品质保持[12-13]。相比之下,基于四极磁铁的磁光学扩束方法可通过横向相空间调控实现束斑静态整形,并在扩束过程中约束束流发散角[14],更适用于紧凑型大面积低发散辐照束线设计。近年来,Halbach永磁结构[15-16]因具有高梯度、紧凑化和低功耗等优势[17],已在康奈尔-布鲁克海文能量回收直线加速器测试装置等系统中得到应用验证[18-19],但其磁场调节能力有限,对制造误差、安装偏差及入射束流参数波动较为敏感[20-21]。因此,结合永磁元件紧凑优势与电磁元件可调能力的混合磁光构型成为一种可行方案。Anania等[22]在ALPHA-X紧凑型FEL束线中采用永磁固定聚焦与电磁动态修正相结合的方案,验证了该类构型的工程可行性。然而,现有研究多侧重于磁体结构设计或束线整体性能验证,对紧凑尺度下多类型磁元件协同作用引起的束流包络演化、发散角压缩机制及横向相空间重构特性仍缺乏系统分析。
针对上述问题,本文面向空间辐射效应地面模拟试验中的大面积电子束辐照需求,提出一种紧凑型扩束压散角束线方案,将Halbach永磁四极磁铁、电四极磁铁与螺线管磁铁集成于有限空间内,以实现束流扩束与发散角压缩的协同调控。通过磁场分布计算和粒子追踪仿真,分析束流包络、发散角及横向相空间的演化规律,阐明不同磁元件在扩束、匹配和压散角过程中的作用机制,并验证该方案在紧凑空间约束下的可行性,研究结果可为受限空间条件下电子束大面积辐照平台的束流光学设计提供参考。

1 束流参数与束线总体设计

1.1 初始束流参数表征

本研究采用包含10 000个宏粒子的相对论电子束作为入射束流。本文所述发散角定义为粒子运动方向相对于束线轴线的夹角,统计时采用横向合成发散角表征。对第$ i $个粒子,令$ x_{i}^{'}={v}_{x,i}/{v}_{z,i} $$ y_{i}^{'}={v}_{y,i}/{v}_{z,i} $,其中,$ {v}_{x,i} $$ {v}_{y,i} $$ {v}_{z,i} $分别为粒子速度在xyz方向的分量,则第$ i $个粒子的横向合成发散角可表示为:
$ {\theta }_{i}=\arctan \left(\sqrt{x_{i}^{'2}+y_{i}^{'2}}\right)\approx \sqrt{x_{i}^{'2}+y_{i}^{'2}} $
本文中“发散角小于0.05 mrad的电子占比”为满足θi < 0.05 mrad的粒子数占总粒子数的比例。根据上述定义统计得到的初始束流主要参数见表1
表 1 初始电子束参数

Table 1 Initial electron beam parameters

参数指标
90%粒子半径/mm0.76
电子束能量/MeV10.07~10.10
90%相对能散/%1.49
x方向均方根发射度/mm·mrad3.878 9
y方向均方根发射度/mm·mrad3.953 6
小于0.05 mrad电子占比/%3.6
图1给出了其初始横向分布及投影粒子数分布。
图 1 初始电子束横向分布及投影粒子数分布

Fig.1 Initial transverse distribution of the electron beam and projected particle-count distributions

图1可见,初始束流横截面近似呈圆对称高斯分布,90%粒子包络半径为0.76 mm。x方向和y方向归一化均方根发射度分别为3.878 9 mm·mrad和3.953 6 mm·mrad,二者数值接近,表明束流在两个横向方向上的相空间分布具有较好的对称性。在mm·mrad单位制下,束流横向相空间的协方差矩阵为:
$ {\Sigma }=\left(\begin{matrix}0.127\;1 & 0.080\;7 & -0.000\;4 & -0.000\;1\\0.080\;7 & 0.327\;3 & -0.000\;1 & -0.002\;3\\-0.000\;4 & -0.000\;1 & 0.129\;6 & 0.085\;9\\-0.000\;1 & -0.002\;3 & 0.085\;9 & 0.338\;1\end{matrix}\right) $
由式(2)可见,∑11与∑33数值接近,说明束流在xy方向的尺寸分布具有良好对称性;此外,表征横向耦合的非对角项数量级远小于主对角项,表明初始束流处于较好的横向解耦状态。

1.2 束线总体设计

本文提出了如图2所示的紧凑型束线结构。束线扩束段由两组Halbach永磁四极磁铁和一组电四极磁铁组成。永磁四极采用钐钴永磁材料和Halbach环形结构,可以在有限空间内提供较高磁场梯度,完成入射束流的初始光学整形;电四极磁铁的磁极和回轭采用电工纯铁,铁轭为正八边形对称布置。通过调节励磁电流对永磁四极磁铁的场强误差、装配误差及束流初始参数引起的光学失配进行补偿,从而提高束线对工程误差和工况变化的适应能力。束线末端设置一组螺线管磁铁,其产生的近轴对称轴向磁场可对束流施加横向聚焦,并在保持束斑尺寸基本稳定的同时压缩发散角。各磁铁单元间的漂移距离作为重要优化变量,与磁场梯度共同决定束流横向尺寸及相空间分布的演化路径。根据Li等[23]提出的两阶段优化方法对束线参数进行优化设计,优化结果如表2所示,表中漂移距离为该磁铁入口距上一磁铁出口的间距。该束线在总长度约1.9 m内实现了束流横向展宽与发散角压缩的协同调控,为空间受限条件下的扩束压散角束线设计提供了可行方案。
图 2 紧凑型扩束压散角束线结构示意

Fig.2 Schematic layout of the compact beam-expansion and divergence-compression beamline

表 2 优化后束线关键参数

Table 2 Key parameters of the optimized beamline

磁铁序号 磁铁类型 漂移距离/m 内半径/m 外半径/m 极面长度/m 电流/A
1 永磁四极磁铁 0.100 0.040 9 0.050 9 0.105 0
2 永磁四极磁铁 0.648 0.072 7 0.086 2 0.027 2
3 电四极磁铁 0.228 0.100 0 0.277 6 0.050 0 18.787
4 螺线管磁铁 0.700 0.060 0 0.075 0 0.150 0 8.415

2 磁场特性与束流动力学分析

2.1 磁场分布

束线各磁性元件的横截面磁场分布如图3所示,束线中心轴上的轴向磁场分布如图4所示。
图 3 各个磁铁横截面磁场分布

Fig.3 Magnetic field distribution across the cross-sections of each magnet

图 4 束线中心轴上的轴向磁场分布

Fig.4 Axial magnetic-field distribution along the beamline centerline

图3可见,两组Halbach永磁四极磁铁均由16块永磁体拼接形成,其孔径内部呈现典型四极场分布,磁场方向沿相邻象限交替变化,其中第一、第二组永磁四极的磁场方向相反。图4中黑色实线表示四极磁铁沿轴向的磁场梯度分布,红色虚线表示螺线管轴向磁通密度分布。可以看出,第一组永磁四极磁铁具有较高磁场梯度7.922 6 T/m,主要用于在入射阶段建立束流扩展所需的初始相空间条件。第二组永磁四极磁铁梯度为−1.144 1 T/m,其符号与第一组相反,与前置0.648 m漂移段共同作用,使束流在自然演化过程中进一步扩展。后段电四极磁铁和螺线管磁铁的梯度幅值分别为0.419 0 T/m和0.086 0 T(峰值轴向磁通密度),主要承担束流横向匹配与压缩发散角功能。

2.2 束流动力学分析

在上述最优磁场参数条件下,采用粒子追踪方法对包含10 000个宏粒子的电子束进行数值模拟,分析束流横向包络及相空间分布随传播距离的演化规律。图5展示了束流在XZYZ平面内的粒子轨迹。在四极磁铁组与漂移段的共同作用下,束流横向包络沿传播方向呈现明显的交替聚焦与发散特征,两个横向平面的包络尺寸逐步趋于一致。在XZ平面内,第一组四极磁铁对X方向产生去聚焦作用,使束流包络开始展开,进入第二组四极磁铁后,束流受到相反方向的聚焦调制,包络演化趋势发生调整,第三组四极磁铁进一步对束流进行整形与补偿,使其在达到目标横向尺寸后满足螺线管入口匹配要求。相应地,在YZ平面内,束流演化规律与XZ平面相反,即依次经历聚焦、发散和再次聚焦。经过三组四极磁铁交替调控后,束流横截面在两个横向平面内的包络趋于一致,横向分布对称性得到改善。
图 5 束流在XZYZ平面内的粒子轨迹演化

Fig.5 Evolution of particle trajectories in the XZ and YZ planes

图6展示了电四极磁铁出口处的横向相空间分布,xy方向相图均呈近线性分布,最大横向尺寸基本一致,且两方向相图在横向位置小于0.02 m范围内高度重合,表明经前三组磁元件协同调控后,束流在两个横向方向上已获得良好的光学匹配状态。
图 6 电四极磁铁出口处横向相空间分布

Fig.6 Transverse phase-space distributions at the exit of the electromagnetic quadrupole

当束流进入螺线管磁场区域后,其横向相空间发生明显旋转,如图7所示。相椭圆与X轴平行,粒子在Y轴方向的分布范围明显收缩,表明螺线管主要通过相空间取向重构实现发散角压缩。
图 7 螺线管出口处xy方向横向相空间分布

Fig.7 Transverse phase-space distributions in the x and y directions at the exit of the solenoid

为进一步定量评估发散角压缩效果,图8对比了初始束流与螺线管出口处束流的发散角分布。统计结果表明,经过束线调控后,发散角小于0.05 mrad的电子比例由初始的3.60%提高至52.31%,表明束线对发散角的压缩效果显著。
图 8 初始与螺线管出口束流发散角分布对比

Fig.8 Comparison of beam divergence distributions at the initial state and solenoid exit

图9给出了束线调控前后的横向空间分布变化,束流90%粒子半径由0.76 mm扩展至约25 mm,束斑横截面积显著增大,从而降低了单位面积内的电子数密度,有助于减弱空间电荷效应对束流包络的影响,改善束线传输及样品区辐照场的稳定性。
图 9 初始与螺线管出口束流横向空间分布对比

Fig.9 Comparison of transverse spatial distributions at the initial state and solenoid exit

上述结果表明,前级四极磁铁组通过交替聚焦—散焦及漂移传输改变束流横向相空间取向,并逐步形成束斑扩展所需的相空间分布;末级螺线管在轴对称磁场作用下进一步驱动横向相空间椭圆旋转,使束流在保持较大横向尺度的同时实现发散角压缩。从元件功能分配来看,Halbach永磁四极磁铁主要用于前级相空间整形和包络扩展。由于该位置处束流横向尺寸较小,小孔径Halbach结构能够在较短轴向长度内提供较高磁场梯度,有利于提高前级整形段的紧凑性并降低功耗需求。后段电四极磁铁则用于束流光学状态的可调匹配,为螺线管段的相空间旋转提供合适入口条件。相较于纯电磁构型,该方案减少了前级强整形段对励磁电源、散热和安装空间的依赖;相较于纯永磁构型,后段电磁元件保留了对入射束流波动和磁体误差的匹配补偿能力。因此,该混合构型在紧凑性、低功耗和可调匹配能力之间取得了较好的平衡。

3 结 语

本文设计了一种面向地面电子束辐照实验的紧凑型扩束压散角束线,通过磁场计算与粒子追踪仿真,分析了束流横向相空间的演化规律,主要结论如下:
(1)本文设计的束线由2组Halbach永磁四极磁铁、1组电四极磁铁和1组螺线管磁铁组成,可在约1.9 m长度内实现电子束扩束与发散角压缩的协同调控。前级四极磁铁组主要完成束流扩束及横向匹配,末级螺线管通过相空间旋转实现发散角压缩。
(2)仿真结果表明,经束线调控后,90%粒子半径由0.76 mm扩展至约25 mm,发散角小于0.05 mrad的电子比例由3.60%提高至52.31%,束流扩展与发散角压缩效果较为显著。
上述结果验证了永磁与电磁混合构型在紧凑尺度下协同实现扩束与发散角压缩的可行性。该构型兼顾了Halbach永磁四极磁铁的紧凑低功耗特性和电磁元件的可调匹配能力,可为紧凑型电子束辐照装置的束流光学设计提供技术支撑。
1
沈自才, 闫德葵. 空间辐射环境工程的现状及发展趋势[J]. 航天器环境工程, 2014, 31 (3): 229- 240.

SHEN Z C, YAN D K. Present status and prospects of space radiation environmental engineering[J]. Spacecraft Environment Engineering, 2014, 31 (3): 229- 240.

2
王会斌, 呼延奇, 郑悦, 等. 航天器空间辐射效应分析技术现状与思考[J]. 航天器环境工程, 2022, 39 (4): 427- 435.

WANG H B, HU Y Q, ZHENG R, et al. The present situation of technologies for analysis of space radiation effects to spacecraft and some retrospection[J]. Spacecraft Environment Engineering, 2022, 39 (4): 427- 435.

3
陶孟泽, 李宏伟, 张龙龙, 等. 5 MeV电子直线加速器建设与空间辐射效应应用[J]. 航天器环境工程, 2025, 42 (3): 263- 269.

DOI

TAO M Z, LI H W, ZHANG L L, et al. Construction of the 5 MeV electron linear accelerator(LINAC) and its application in space radiation effects[J]. Spacecraft Environment Engineering, 2025, 42 (3): 263- 269.

DOI

4
韩虹, 王子玉, 钟兴旺. 空间环境对飞行器的影响[J]. 空间电子技术, 2011, 8 (1): 1- 5.

HAN H, WANG Z Y, ZHONG X W. Influence of space environment on spacecrafts[J]. Space Electronic Technology, 2011, 8 (1): 1- 5.

5
WANG C Q, XIANG Y H, LI J H, et al. Design of beam homogenization for ground simulation of space electron environment[J]. Nuclear Techniques, 2020, 43 (12): 120201.

DOI

王常强, 向益淮, 李金海, 等. 地面模拟空间电子环境束流均匀化设计[J]. 核技术, 2020, 43 (12): 120201.

DOI

6
姚海波. 相对论电子束长程稳定传输性能研究[D]. 长沙: 国防科技大学, 2021.

YAO H B. The Long-distance and stable transport of relativistic electron beams[D]. Changsha:National University of Defense Technology,2021.

7
KIM D, ANDEREWS H L, CHOI B K, et al. Divergence study and emittance measurements for the electron beam emitted from a diamond pyramid[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2020, 953, 163055.

DOI

8
李金海, 张立锋, 曾自强. 电子辐照加速器束流扫描均匀化研究[J]. 原子能科学技术, 2015, 49 (2): 564- 568.

DOI

LI J H, ZHANG L F, ZENG Z Q. Research of homogenizing scanned beam of electron irradiation accelerator[J]. Atomic Energy Science and Technology, 2015, 49 (2): 564- 568.

DOI

9
DENHAM P, YANG Y, GUO V, et al. High energy electron diffraction instrument with tunable camera length[J]. Structural Dynamics, 2024, 11 (2): 024302.

DOI

10
THIEBERGER P, STASSINOPULOS E G, Van GUNTEN O, et al. Heavy-ion beams for single-event research at brookhaven—present and future[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 1991, 56, 1251- 1255.

11
陈泉. 串列加速器R20支线重离子辐照装置的束流光学计算[R]. 北京: 中国原子能科学研究院年报, 2009(1): 201-203.

12
杨新宇, 郭刚, 彭朝华. 基于八极磁铁的束流均匀化设计方法[J]. 现代应用物理, 2019, 68 (5): 054104- 1-054104-8.

DOI

HAN J H, GUO G, LIU J C, et al. Design of 100-MeV proton beam spreading scheme with double-ring double scattering method[J]. Acta Physica Sinica, 2019, 68 (5): 054104- 1-054104-8.

DOI

13
韩金华, 郭刚, 刘建成, 等. 100 MeV质子双环双散射体扩束方案设计[J]. 物理学报, 2019, 68 (5): 054104- 1-054104-8.

DOI

HAN J H, GUO G, LIU J C, et al. Design of 100-MeV proton beam spreading scheme with double-ring double scattering method[J]. Acta Physica Sinica, 2019, 68 (5): 054104- 1-054104-8.

DOI

14
WEINGARTNER R, FUCHS M, POPP A, et al. Imaging laser-wakefield-accelerated electrons using miniature magnetic quadrupole lenses[J]. Physical Review Special Topics—Accelerators and Beams, 2011, 14 (5): 052801.

DOI

15
BJORK R, SMITH A, BAHL C R H. Analysis of the magnetic field, force, and torque for two-dimensional Halbach cylinders[J]. Journal of Magnetism and Magnetic Materials, 2010, 322 (1): 133- 141.

DOI

16
BLUMLER P, SOLTNER H. Practical concepts for design, construction and application of Halbach magnets in magnetic resonance[J]. Applied Magnetic Resonance, 2023, 54 (11): 1701- 1739.

17
HU X, TIE W, SHAO Y, et al. Design and optimization of the Halbach cylinder for engineering application scenarios[J]. Journal of Magnetism and Magnetic Materials, 2025, 616, 172845.

DOI

18
TSOUPAS N, BERG J S, BROOKS S, et al. Design of a modified Halbach magnet for the CBETA project[R]. Upton, NY: Brookhaven National Laboratory, 2018.

19
BROOKS S. CBETA permanent magnet production run[R]. Upton, NY: Brookhaven National Laboratory, 2019.

20
SCHILLACI F, MAGGIORE M, RIFUGGIATO D, et al. Errors and optics study of a permanent magnet quadrupole system[J]. Journal of Instrumentation, 2015, 10 (5): T05001.

DOI

21
LI Y, HU X, TIE W, et al. Design and performance validation of a compact hybrid-magnet beam expansion and collimation system[J]. AIP Advances, 2026, 16 (1): 015206.

DOI

22
ANANIA M P, CLARK D, ISSAC R C, et al. The ALPHA-X beam line: Toward a compact FEL[C]. Proceedings of the 31st International Free Electron Laser Conference, Liverpool, UK, August 23—28, 2009. 172-175.

23
LI Y, HU X, TIE W, et al. Two-stage optimization method for beam optical systems based on multi-objective intelligent optimization algorithms[J]. IEEE Transactions on Nuclear Science, 2025, 72 (12): 3624- 3635.

DOI

文章导航

/