CN113849953A - 一种用于空间x射线通信的微焦斑装置设计优化方法 - Google Patents

一种用于空间x射线通信的微焦斑装置设计优化方法 Download PDF

Info

Publication number
CN113849953A
CN113849953A CN202110728154.2A CN202110728154A CN113849953A CN 113849953 A CN113849953 A CN 113849953A CN 202110728154 A CN202110728154 A CN 202110728154A CN 113849953 A CN113849953 A CN 113849953A
Authority
CN
China
Prior art keywords
ray
photon
focal spot
incident
lens
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110728154.2A
Other languages
English (en)
Other versions
CN113849953B (zh
Inventor
李连升
梅志武
王立
石永强
左富昌
陈建武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Control Engineering
Original Assignee
Beijing Institute of Control Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Control Engineering filed Critical Beijing Institute of Control Engineering
Priority to CN202110728154.2A priority Critical patent/CN113849953B/zh
Publication of CN113849953A publication Critical patent/CN113849953A/zh
Application granted granted Critical
Publication of CN113849953B publication Critical patent/CN113849953B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/40—Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明提供了一种用于空间X射线通信的微焦斑装置设计优化方法，该方法以光学系统的微焦斑均方根半径为优化目标，以视场、离焦量和镜片倾斜量为设计变量，充分考虑X射线光子能量和反射率的特征信息，构建了空间X射线通信的微焦斑装置设计优化模型，基于蒙特卡洛法实现空间X射线通信用的微焦斑装置的设计优化。本发明一方面避免传统X射线光学仿真分析方法仅考虑单一能量的X射线光子，且不考虑X光线光子反射率的缺陷；另一方面采用蒙特卡洛法实现不同参数组合条件下X射线聚焦光学系统的微焦斑优化，获得高增益的X射线通信用微焦斑装置的设计优化，为设计人员提供更多合理的选择。

Description

一种用于空间X射线通信的微焦斑装置设计优化方法
技术领域
本发明涉及航天器产品设计技术领域，特别涉及一种用于空间X射线通信的微焦斑装置设计优化方法。
背景技术
空间聚焦型X射线探测器作为深空探测、脉冲星计时与导航、X射线通信等领域的共性基础产品，需求日益迫切。深空探测对高通量实时通信的迫切需求亟需人类开展基于短波光学通信的新体制研究。X射线频率高达 1018Hz，可以实现高通量、低色散、低损耗通信，是一种极具发展前景的通信技术，目前美国已在国际空间站、ESA在水星探测计划中开展X射线通信技术验证。
此外，X射线通信以其高通量、高穿透性、低损耗、无色散等特性在载人返回舱、高超声速飞行器等领域的黑障区通信方面极具发展潜力。飞行器以高超音速再入大气层，由于粘性流和激波加热的作用，飞行器表面附近的空气被电离，形成包裹飞行器的“等离子鞘套”。等离子体对于电磁波有着强烈的反射和吸收的作用，出现通信“黑障”现象。X射线可几乎以无衰减地穿过飞行器外包裹的等离子鞘套，可解决超高速飞行器或载人返回舱再入大气层遇到的“黑障区”通信问题，具有重要军事意义。
其中，X射线通信光学系统受ATP指向精度与稳定度的影响，聚焦式X 射线光学既需要高准直发射，也需要微焦斑聚焦。传统的掠入射聚焦X射线光学系统(Wolter-I)具有良好的轴上聚焦性能，但随着离轴角度增大，其聚焦性能急剧下降，严重影响了X射线天线的光学增益。传统的掠入射聚焦X 射线光学系统已无法满足空间X射线通信等领域的应用需求。
发明内容
本发明所要解决的技术问题是：克服现有Wolter-I聚焦光学系统的不足，提供了一种空间X射线通信微焦斑装置的设计优化方法，该方法基于所构建的微焦斑最小化优化模型，基于光线矢量方法建立X射线光线传输方程，采用基于蒙特卡洛和X射线反射理论进行多种工况下的聚焦性能仿真分析，通过多变量组合分析获取三种典型工况下的最小化焦斑均方根半径，获得优化设计变量值。
本发明所采用的技术解决方案：
一种用于空间X射线通信的微焦斑装置设计优化方法，包括以下步骤：
(1)确定空间X射线通信微焦斑装置设计优化参数，所述设计优化参数包括X射线镜片倾斜量r、视场F和离焦量Dec为设计变量；
(2)根据步骤(1)确定的空间X射线通信微焦斑装置的设计优化参数，构建空间X射线通信微焦斑装置优化模型；
(3)初始参数化，设置P个X射线光子在光学镜头内表面上的入射位置、光子能量和掠入射角；
(4)确定设计变量取值范围；
(5)利用蒙特卡洛光学追踪法和X射线全反射理论，根据步骤(3)设置的X射线光子入射位置坐标，X光子沿直线掠入射至光学镜片的内表面，通过计算入射点的曲率半径以及入射点到光学镜头中心轴线的距离，得到每个X射线光子的实际掠入射角；
(6)根据每个X射线光子的光子能量计算所述光子的临界入射角，计算得到第p个X射线光子的临界入射角为φp；
(7)通过比较每个X射线光子的临界入射角φp与所述光子的实际掠入射角
Figure RE-GDA0003343705270000021
判定该入射的X射线光子是否发生全反射，如果
Figure RE-GDA0003343705270000022
则认定该X射线光子发生全反射，并统计全反射的光子数Ntotal；
(8)计算所述全反射的X射线光子Ntotal在光学镜片内表面的反射角；
(9)根据步骤(8)计算获得所述光子在镜头内表面上的反射角，根据三角函数计算每个光子到达探测器焦平面的位置坐标；
(10)X射线微焦斑聚焦性能分析，开展空间X射线通信微焦斑装置设计。
进一步的，所述步骤(2)中，空间X射线通信的微焦斑装置设计优化模型建立过程如下：
(2a)以50％视场～100％视场范围焦斑均方根半径最小为优化目标；
(2b)设定典型工作场景，F＝50％视场、80％视场、100％视场，同时在仿真分析过程中镜片倾斜量和离焦量组合变化，设每个X光子在焦平面上的位置坐标(X3i、Y3i、Z3i)，坐标系原点位于焦平面中心；i＝1,2,3,…,Nline；计算每个X光子的位置与探测器中心的距离：
Figure RE-GDA0003343705270000031
(2c)根据步骤(2b)构建空间X射线通信的微焦斑装置设计优化模型：
Figure RE-GDA0003343705270000032
进一步的，所述步骤(5)中，蒙特卡洛光学追踪法和X射线全反射计算，具体过程如下：
(5a)在X射线光学系统前端距离d处建立坐标系，模拟脉冲星辐射的X 射线光子从无穷远处以平行光入射，所构建的坐标系以光轴方向为Z轴，XY 轴满足右手定则；
(5b)在d处垂直光轴的平面内随机生成N个样本点，与光轴Z以夹角θ入射X射线光学镜头内表面，经单次反射后聚焦于焦平面；
(5c)构建入射光线、反射光线矢量方程：
设定入射光子三维坐标(X1i、Y1i、Z1i)，入射光线与Z轴夹角为Fov，因此入射光线的直线方程可表示为：
Figure RE-GDA0003343705270000041
(5d)构建X射线光学系统面形方程：
X射线掠入射反射镜的面形公式为：
Figure RE-GDA0003343705270000042
其中，c表示面形顶点曲率，k是二次曲面系数，当k＝-1时，公式表示的面形是抛物面；r是坐标系径向半径，z是面形矢高；
(5e)计算入射光线与X射线光学系统面形交点：
联立空间直线与X射线光学系统面形方程，即可求出X光子入射到X射线光学系统面形的交点(X2i、Y2i、Z2i)；
(5f)对面形公式计算偏导，获得每个入射点的法向量：
Figure RE-GDA0003343705270000043
(5g)构建法线方程
Figure RE-GDA0003343705270000044
(5h)计算反射光线
设定反射光线的方向向量为(ni,nj,nk)，反射光线对应的直线方程为：
Figure RE-GDA0003343705270000045
根据入射角等于反射角，向量积法则，|α×β|＝|α|·|β|sinθ,θ＝(α，β)，根据该式求出反射光线的方向向量，入射向量与切平面法向量的差积，该差积与反射光线向量垂直，点积为0，据此可构建X射线反射光线的直线方程；
(5i)通过上述过程计算反射光线与焦平面的交点(X3i、Y3i、Z3i)
根据所述步骤(3)以及焦平面所在平面方程，即可计算出每个聚焦X射线光子在焦平面上的三维坐标(X3i、Y3i、Z3i)。
进一步的，所述步骤(10)中，对X射线微焦斑聚焦性能分析，具体评价过程如下：
(9a)根据所述(5j)计算出的反射光线与焦平面的焦点(X3i、Y3i、Z3i)，计算该焦点距离焦平面几何中心位置的绝对距离，计算公式如下：
Figure RE-GDA0003343705270000052
(9b)根据镜片倾斜量和离焦量的矩阵组合，综合分析多组合下每次仿真分析所对应的di：
(9c)计算三种典型工作场景下的焦斑均方根最小化值，并确定所对应的设计参数值。
进一步的，步骤(2)中，优化模型以微焦斑均方根半径为优化目标，建立优化模型：
Figure RE-GDA0003343705270000051
其中，k＝1、2、3，分别对应50％视场、80％视场和100％视场三种典型工况，优化目标函数表示为Min(RMS)，即最小化X射线光学聚焦微焦斑均方根半径RMS；di为第i个发生全反射的X射线光子到达探测器焦平面后的位置与探测器中心之间的距离。
进一步的，步骤(3)中，第p个光子的入射位置坐标分别为xp、yp、zp，所述坐标系的原点设定为探测器中心，Z轴为光学镜头的中心轴线；第p个X 射线光子的能量为Ep，Ep的取值范围为Emin～Emax，且服从随机分布；第p个光子的掠入射角为θp，θp在设定的角度范围θmin～θmax内随机分布；p＝1、2、…、 P，P为设定的X射线光子样本量。
进一步的，步骤(5)中，第p个X射线光子的实际掠入射角
Figure RE-GDA0003343705270000065
的计算公式如下：
Figure RE-GDA0003343705270000061
其中，其中，Rp为第p个光子入射点处的镜头曲率半径；dp为第p个光子在光学镜头内表面入射点到光学镜头中心轴线的距离。
进一步的，步骤(4)中，设定X射线聚焦光学的视场F分别为50％、80％和100％三种典型应用场景，设定微焦斑的前后移动范围为-10mm≤Dec≤ 10mm，镜片倾斜量Tilt取值范围为[-5’,5’]。
进一步的，步骤(6)中，
Figure RE-GDA0003343705270000062
其中，f1为光学镜头材料的散射因子。
进一步的，步骤(8)中，
Figure RE-GDA0003343705270000063
其中，αq为第q个发生全反射的X射线光子在镜头内表面上的反射角，
Figure RE-GDA0003343705270000064
为第q个全反射X射线光子的实际掠入射角，θq′为第q个全反射X射线光子的掠入射角；其中，q＝1，2，…，Ntotal。
进一步的，步骤(10)中，综合变化镜片倾斜量和离焦量并计算每种仿真分析结果所对应的聚焦焦斑RMS均方根半径，记录每轮的RMS最小值，然后综合计算出所有最小值RMS中的最小值，及其对应的设计参数，开展空间X 射线通信微焦斑装置设计。
本发明与现有技术相比的优点在于：
(1)在本发明的X射线聚焦光学系统中，首次构建了以三种典型工况下的微焦斑均方根最小化为目标的优化模型，确定了以镜片倾斜量和离焦量多变量组合的总体仿真分析方案，为设计与决策人员提供了更多选择；
(2)本发明基于光线矢量方程构建了任意面形X射线掠入射光学系统仿真方程，解决了传统光学仿真软件无法同时模拟X射线光子能量和反射率等特征信息的难题，可实现高效率、多维度信息的一体化仿真分析，更接近于工程实际，为工程设计人员提供高可信的技术支撑；
(3)本发明采用蒙特卡洛法和X射线全反射理论实现不同参数组合条件下 X射线聚焦光学系统的微焦斑仿真分析，可实现高增益的X射线通信用微焦斑装置的设计优化结果，提高了航天产品开发效率，降低了研制成本。
附图说明
图1为本发明用于空间X射线通信的微焦斑装置设计优化方法的流程图；
图2为本发明中X射线光学系统光线追迹模型示意图；
图3为本发明中X射线掠入射光线追迹仿真分析图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述：
如图1-3所述，本发明对空间X射线通信的微焦斑装置进行设计优化，以光学系统的微焦斑均方根半径为优化目标，以镜片倾斜量和离焦量为设计变量，充分考虑X射线光子能量和反射率的特征信息，针对50％、80％和100％视场三种典型工况，构建了空间X射线通信的微焦斑装置设计优化模型，采用蒙特卡洛法和X射线全反射理论进行光线追迹，进行多组合参数仿真分析与优化。本发明避免传统X射线光学仿真分析方法仅考虑单一能量的X射线光子，且不考虑X光线光子反射率的缺陷；采用蒙特卡洛法实现不同参数组合条件下X射线聚焦光学系统的微焦斑优化，获得高增益的X射线通信用微焦斑装置的设计优化，为设计人员提供更多合理的选择。
如图1所示的流程图，本发明的一种用于空间X射线通信的微焦斑装置设计优化方法，包括以下步骤：
一种用于空间X射线通信的微焦斑装置设计优化方法，包括以下步骤：
(1)确定空间X射线通信微焦斑装置设计优化参数，所述设计优化参数包括X射线镜片倾斜量r、视场F和离焦量Dec为设计变量；
(2)根据步骤(1)确定的空间X射线通信微焦斑装置的设计优化参数，构建空间X射线通信微焦斑装置优化模型，优化模型以微焦斑均方根半径为优化目标，建立优化模型：
Figure RE-GDA0003343705270000081
其中，k＝1、2、3，分别对应50％视场、80％视场和100％视场三种典型工况，优化目标函数表示为Min(RMS)，即最小化X射线光学聚焦微焦斑均方根半径(RMS)；di为第i个发生全反射的X射线光子到达探测器焦平面后的位置与探测器中心之间的距离。
(3)初始参数化，设置P个X射线光子在光学镜头内表面上的入射位置、光子能量和掠入射角，其中，第p个光子的入射位置坐标分别为xp、yp、zp，所述坐标系的原点设定为探测器中心，Z轴为光学镜头的中心轴线；第p个X 射线光子的能量为Ep，Ep的取值范围为Emin～Emax，且服从随机分布；第p个光子的掠入射角为θp，θp在设定的角度范围θmin～θmax内随机分布；p＝1、2、…、 P，P为设定的X射线光子样本量；
(4)确定设计变量取值范围，设定X射线聚焦光学的视场F分别为50％、 80％和100％三种典型应用场景，设定微焦斑的前后移动范围为-10mm≤Dec≤ 10mm，镜片倾斜量Tilt取值范围为[-5’,5’]。
(5)利用蒙特卡洛光学追踪法和X射线全反射理论，根据步骤(3)设置的X射线光子入射位置坐标，X光子沿直线掠入射至光学镜片的内表面，通过计算入射点的曲率半径以及入射点到光学镜头中心轴线的距离，得到每个X射线光子的实际掠入射角。第p个X射线光子的实际掠入射角
Figure RE-GDA0003343705270000091
的计算公式如下：
Figure RE-GDA0003343705270000092
其中，其中，Rp为第p个光子入射点处的镜头曲率半径；dp为第p个光子在光学镜头内表面入射点到光学镜头中心轴线的距离。
(6)根据每个X射线光子的光子能量计算所述光子的临界入射角，其中，计算得到第p个X射线光子的临界入射角为φp，p＝1、2、…、P；
Figure RE-GDA0003343705270000093
其中，f1为光学镜头材料的散射因子。
(7)通过比较每个X射线光子的临界入射角φp与所述光子的实际掠入射角
Figure RE-GDA0003343705270000094
判定该入射的X射线光子是否发生全反射，如果
Figure RE-GDA0003343705270000095
则认定该X射线光子发生全反射，并统计全反射的光子数Ntotal；
(8)计算所述全反射的X射线光子Ntotal在光学镜片内表面的反射角：
Figure RE-GDA0003343705270000096
其中，αq为第q个发生全反射的X射线光子在镜头内表面上的反射角，
Figure RE-GDA0003343705270000097
为第q个全反射X射线光子的实际掠入射角，θq′为第q个全反射X射线光子的掠入射角；其中，q＝1，2，…，Ntotal；
(9)根据步骤(8)计算获得所述光子在镜头内表面上的反射角，根据三角函数计算每个光子到达探测器焦平面的位置坐标。
(10)X射线微焦斑聚焦性能分析，综合变化镜片倾斜量和离焦量并计算每种仿真分析结果所对应的聚焦焦斑(RMS)均方根半径。记录每轮的RMS最小值，然后综合计算出所有最小值RMS中的最小值，及其对应的设计参数，开展空间X射线通信微焦斑装置设计。
步骤(2)中，空间X射线通信的微焦斑装置设计优化模型建立过程如下：
(2a)以50％视场～100％视场范围焦斑均方根半径最小为优化目标；
(2b)设定典型工作场景，F＝50％视场、80％视场、100％视场，同时在仿真分析过程中镜片倾斜量和离焦量组合变化，设每个X光子在焦平面上的位置坐标(X3i、Y3i、Z3i)，坐标系原点位于焦平面中心；i＝1,2,3,…,Nline；计算每个X光子的位置与探测器中心的距离：
Figure RE-GDA0003343705270000101
(2c)根据步骤(2b)构建空间X射线通信的微焦斑装置设计优化模型：
Figure RE-GDA0003343705270000102
步骤(5)中，蒙特卡洛光学追踪法和X射线全反射计算，具体过程如下：
(5a)在X射线光学系统前端距离d处建立坐标系(d＝100mm)，模拟脉冲星辐射的X射线光子从无穷远处以平行光入射，所构建的坐标系以光轴方向为Z轴，XY轴满足右手定则；
(5b)在d处垂直光轴的平面内随机生成N个样本点，与光轴Z以夹角θ入射X射线光学镜头内表面，经单次反射后聚焦于焦平面；
(5c)构建入射光线、反射光线矢量方程：
设定入射光子三维坐标(X1i、Y1i、Z1i)，入射光线与Z轴夹角为Fov，因此入射光线的直线方程可表示为：
Figure RE-GDA0003343705270000111
(5d)构建X射线光学系统面形方程：
X射线掠入射反射镜的面形公式为：
Figure RE-GDA0003343705270000112
其中，c表示面形顶点曲率，k是二次曲面系数，当k＝-1时，公式表示的面形是抛物面；r是坐标系径向半径，z是面形矢高；
(5e)计算入射光线与X射线光学系统面形交点：
联立空间直线与X射线光学系统面形方程，即可求出X光子入射到X射线光学系统面形的交点(X2i、Y2i、Z2i)；
(5f)计算每个入射点的法向量：
Figure RE-GDA0003343705270000113
(5g)构建法线方程
Figure RE-GDA0003343705270000114
(5h)计算反射光线
设定反射光线的方向向量为(ni,nj,nk)，反射光线对应的直线方程为：
Figure RE-GDA0003343705270000115
根据入射角等于反射角，向量积法则，|α×β|＝|α|·|β|sinθ,θ＝(α，β)，根据该式求出反射光线的方向向量，入射向量与切平面法向量的差积：(该差积与反射光线向量垂直，点积为0)，据此可构建X射线反射光线的直线方程；
(5i)通过上述过程计算反射光线与焦平面的交点(X3i、Y3i、Z3i)
根据所属(3)以及焦平面所在平面方程，即可计算出每个聚焦X射线光子在焦平面上的三维坐标(X3i、Y3i、Z3i)。
步骤(10)中，对X射线微焦斑聚焦性能分析，具体评价过程如下：
(9a)根据所属(5j)计算出的反射光线与焦平面的焦点(X3i、Y3i、Z3i)，计算该焦点距离焦平面几何中心位置的绝对距离，计算公式如下：
Figure RE-GDA0003343705270000121
(9b)根据镜片倾斜量和离焦量的矩阵组合，综合分析多组合下每次仿真分析所对应的di：
(9b)计算三种典型工作场景下的焦斑均方根最小化值，并确定所对应的设计参数值。
图2为采用X射线全反射理论和矢量法进行的蒙特卡洛法光线追迹过程示意图。其中，X1Y1Z1坐标系是X射线光子的起始位置，坐标系X2Y2Z2是X射线光子击中X光学系统内表面，并产生全反射，坐标系X3Y3Z3是焦平面所在的位置，也是X射线光子聚焦的位置，可以通过统计每个光子到达焦平面的坐标与其中心O3的距离计算焦斑均方根半径(RMS)。
图3是基于本专利所提方法进行的X射线掠入射光线追迹仿真分析图。其中，X和Y轴分别表示光路横截面坐标轴，Z轴为X射线聚焦光学系统的光轴。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。

Claims (11)

1.一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于，包括以下步骤：
(1)确定空间X射线通信微焦斑装置设计优化参数，所述设计优化参数包括X射线镜片倾斜量r、视场F和离焦量Dec为设计变量；
(2)根据步骤(1)确定的空间X射线通信微焦斑装置的设计优化参数，构建空间X射线通信微焦斑装置优化模型；
(3)初始参数化，设置P个X射线光子在光学镜头内表面上的入射位置、光子能量和掠入射角；
(4)确定设计变量取值范围；
(5)利用蒙特卡洛光学追踪法和X射线全反射理论，根据步骤(3)设置的X射线光子入射位置坐标，X光子沿直线掠入射至光学镜片的内表面，通过计算入射点的曲率半径以及入射点到光学镜头中心轴线的距离，得到每个X射线光子的实际掠入射角；
(6)根据每个X射线光子的光子能量计算所述光子的临界入射角，计算得到第p个X射线光子的临界入射角为φp；
(7)通过比较每个X射线光子的临界入射角φp与所述光子的实际掠入射角
Figure RE-FDA0003343705260000011
判定该入射的X射线光子是否发生全反射，如果
Figure RE-FDA0003343705260000012
则认定该X射线光子发生全反射，并统计全反射的光子数Ntotal；
(8)计算所述全反射的X射线光子Ntotal在光学镜片内表面的反射角；
(9)根据步骤(8)计算获得所述光子在镜头内表面上的反射角，根据三角函数计算每个光子到达探测器焦平面的位置坐标；
(10)X射线微焦斑聚焦性能分析，开展空间X射线通信微焦斑装置设计。
2.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：所述步骤(2)中，空间X射线通信的微焦斑装置设计优化模型建立过程如下：
(2a)以50％视场～100％视场范围焦斑均方根半径最小为优化目标；
(2b)设定典型工作场景，F＝50％视场、80％视场、100％视场，同时在仿真分析过程中镜片倾斜量和离焦量组合变化，设每个X光子在焦平面上的位置坐标(X3i、Y3i、Z3i)，坐标系原点位于焦平面中心；i＝1,2,3,…,Nline；计算每个X光子的位置与探测器中心的距离：
Figure RE-FDA0003343705260000021
(2c)根据步骤(2b)构建空间X射线通信的微焦斑装置设计优化模型：
Figure RE-FDA0003343705260000022
3.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：
所述步骤(5)中，蒙特卡洛光学追踪法和X射线全反射计算，具体过程如下：
(5a)在X射线光学系统前端距离d处建立坐标系，模拟脉冲星辐射的X射线光子从无穷远处以平行光入射，所构建的坐标系以光轴方向为Z轴，XY轴满足右手定则；
(5b)在d处垂直光轴的平面内随机生成N个样本点，与光轴Z以夹角θ入射X射线光学镜头内表面，经单次反射后聚焦于焦平面；
(5c)构建入射光线、反射光线矢量方程：
设定入射光子三维坐标(X1i、Y1i、Z1i)，入射光线与Z轴夹角为Fov，因此入射光线的直线方程可表示为：
Figure RE-FDA0003343705260000031
(5d)构建X射线光学系统面形方程：
X射线掠入射反射镜的面形公式为：
Figure RE-FDA0003343705260000032
其中，c表示面形顶点曲率，k是二次曲面系数，当k＝-1时，公式表示的面形是抛物面；r是坐标系径向半径，z是面形矢高；
(5e)计算入射光线与X射线光学系统面形交点：
联立空间直线与X射线光学系统面形方程，即可求出X光子入射到X射线光学系统面形的交点(X2i、Y2i、Z2i)；
(5f)对面形公式计算偏导，获得每个入射点的法向量：
Figure RE-FDA0003343705260000033
(5g)构建法线方程
Figure RE-FDA0003343705260000034
(5h)计算反射光线
设定反射光线的方向向量为(ni,nj,nk)，反射光线对应的直线方程为：
Figure RE-FDA0003343705260000035
根据入射角等于反射角，向量积法则，|α×β|＝|α|·|β|sinθ,θ＝(α，β)，根据该式求出反射光线的方向向量，入射向量与切平面法向量的差积，该差积与反射光线向量垂直，点积为0，据此可构建X射线反射光线的直线方程；
(5i)通过上述过程计算反射光线与焦平面的交点(X3i、Y3i、Z3i)
根据所述步骤(3)以及焦平面所在平面方程，即可计算出每个聚焦X射线光子在焦平面上的三维坐标(X3i、Y3i、Z3i)。
4.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：
所述步骤(10)中，对X射线微焦斑聚焦性能分析，具体评价过程如下：
(9a)根据所述(5j)计算出的反射光线与焦平面的焦点(X3i、Y3i、Z3i)，计算该焦点距离焦平面几何中心位置的绝对距离，计算公式如下：
Figure RE-FDA0003343705260000041
(9b)根据镜片倾斜量和离焦量的矩阵组合，综合分析多组合下每次仿真分析所对应的di：
(9c)计算三种典型工作场景下的焦斑均方根最小化值，并确定所对应的设计参数值。
5.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：步骤(2)中，优化模型以微焦斑均方根半径为优化目标，建立优化模型：
Figure RE-FDA0003343705260000042
其中，k＝1、2、3，分别对应50％视场、80％视场和100％视场三种典型工况，优化目标函数表示为Min(RMS)，即最小化X射线光学聚焦微焦斑均方根半径RMS；di为第i个发生全反射的X射线光子到达探测器焦平面后的位置与探测器中心之间的距离。
6.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：步骤(3)中，第p个光子的入射位置坐标分别为xp、yp、zp，所述坐标系的原点设定为探测器中心，Z轴为光学镜头的中心轴线；第p 个X射线光子的能量为Ep，Ep的取值范围为Emin～Emax，且服从随机分布；第p个光子的掠入射角为θp，θp在设定的角度范围θmin～θmax内随机分布；p＝1、2、…、P，P为设定的X射线光子样本量。
7.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：步骤(5)中，第p个X射线光子的实际掠入射角
Figure RE-FDA0003343705260000055
的计算公式如下：
Figure RE-FDA0003343705260000051
其中，其中，Rp为第p个光子入射点处的镜头曲率半径；dp为第p个光子在光学镜头内表面入射点到光学镜头中心轴线的距离。
8.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：步骤(4)中，设定X射线聚焦光学的视场F分别为50％、80％和100％三种典型应用场景，设定微焦斑的前后移动范围为-10mm≤Dec≤10mm，镜片倾斜量Tilt取值范围为[-5’,5’]。
9.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：步骤(6)中，
Figure RE-FDA0003343705260000052
其中，f1为光学镜头材料的散射因子。
10.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：步骤(8)中，
Figure RE-FDA0003343705260000053
其中，αq为第q个发生全反射的X射线光子在镜头内表面上的反射角，
Figure RE-FDA0003343705260000054
为第q个全反射X射线光子的实际掠入射角，θq′为第q个全反射X射线光子的掠入射角；其中，q＝1，2，…，Ntotal。
11.根据权利要求1所述的一种用于空间X射线通信的微焦斑装置设计优化方法，其特征在于：步骤(10)中，综合变化镜片倾斜量和离焦量并计算每种仿真分析结果所对应的聚焦焦斑RMS均方根半径，记录每轮的RMS最小值，然后综合计算出所有最小值RMS中的最小值，及其对应的设计参数，开展空间X射线通信微焦斑装置设计。
CN202110728154.2A 2021-06-29 2021-06-29 一种用于空间x射线通信的微焦斑装置设计优化方法 Active CN113849953B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110728154.2A CN113849953B (zh) 2021-06-29 2021-06-29 一种用于空间x射线通信的微焦斑装置设计优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110728154.2A CN113849953B (zh) 2021-06-29 2021-06-29 一种用于空间x射线通信的微焦斑装置设计优化方法

Publications (2)

Publication Number Publication Date
CN113849953A true CN113849953A (zh) 2021-12-28
CN113849953B CN113849953B (zh) 2024-06-25

Family

ID=78975087

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110728154.2A Active CN113849953B (zh) 2021-06-29 2021-06-29 一种用于空间x射线通信的微焦斑装置设计优化方法

Country Status (1)

Country Link
CN (1) CN113849953B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116149053A (zh) * 2023-03-16 2023-05-23 哈尔滨工业大学 一种极紫外波段光同时实现单点聚焦及波前分割的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008076399A2 (en) * 2006-12-15 2008-06-26 Hand Held Products, Inc. Apparatus and method comprising deformable lens element
US20100091947A1 (en) * 2008-10-10 2010-04-15 Niu han-ben Differential Interference Phase Contrast X-ray Imaging System
CN104502376A (zh) * 2014-12-22 2015-04-08 北京师范大学 X射线纳米成像设备及成像分析系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008076399A2 (en) * 2006-12-15 2008-06-26 Hand Held Products, Inc. Apparatus and method comprising deformable lens element
US20100091947A1 (en) * 2008-10-10 2010-04-15 Niu han-ben Differential Interference Phase Contrast X-ray Imaging System
CN104502376A (zh) * 2014-12-22 2015-04-08 北京师范大学 X射线纳米成像设备及成像分析系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LI L: "Monte Carlo-based multiphysics coupling analysis of x-ray pulsar telescope", SENSORS, SYSTEMS, AND NEXT-GENERATION SATELLITES XIX. SPIE, 30 December 2015 (2015-12-30) *
王丽娟: ".高功率密度微焦斑X射线源聚焦系统的设计", 电子显微学报, vol. 38, no. 2, 15 February 2019 (2019-02-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116149053A (zh) * 2023-03-16 2023-05-23 哈尔滨工业大学 一种极紫外波段光同时实现单点聚焦及波前分割的方法

Also Published As

Publication number Publication date
CN113849953B (zh) 2024-06-25

Similar Documents

Publication Publication Date Title
Wendelin et al. SolTrace: a ray-tracing code for complex solar optical systems
CN111060899B (zh) 星地一体化激光雷达回波波形仿真方法及系统
CN109557028B (zh) 一种具有密集光斑图案的多通池
EP2887092A1 (en) Computing radar cross section
Andrews et al. Beam wander effects on the scintillation index of a focused beam
CN105425394A (zh) 高能量高准直角太阳模拟器光学系统
CN109407310B (zh) 一种多通池的设计方法
CN113849953B (zh) 一种用于空间x射线通信的微焦斑装置设计优化方法
CN108051908A (zh) 一种成像光学系统
Peng et al. Simulation of transmitted X-rays in a polycapillary X-ray lens
Pompea et al. Stray light analysis of the Apache Point Observatory 3.5-m telescope system
Landman Optical performance of the reflective surface profile of a heliostat
CN111625957B (zh) 塔式太阳能镜场接收器平面的辐射能密度分布模拟方法
CN115600375A (zh) 一种高精度定日镜能流密度计算方法
Yushkin et al. Mathematical model of orbital and ground-based cross-dispersion spectrographs
CN112380643B (zh) 一种等离子体包覆目标的近场电磁散射建模方法
Gao et al. A low-time complexity semi-analytic Monte Carlo radiative transfer model: Application to optical characteristics of complex spatial targets
Xia [Retracted] Adjustment Method of “FAST” Active Reflector Based on Optimal Fitting Strategy
Early et al. A methodology for reflected laser beam hazard analyses
RU2461929C1 (ru) Способ оптимального размещения и ориентации приемного/передающего излучателя в виде коаксиально расположенных диэлектриков цилиндрической формы в фокальной области используемых коллимирующих поверхностей
Steadman et al. GUERAP III simulation of stray light phenomena
CN116430581B (zh) 部分相干光经过自由空间及光学元件的传播的计算方法
Ho et al. Analytical methods to evaluate flux distributions from point-focus collectors for solar furnace and dish engine applications
Prokhorov et al. Modeling of integrating spheres for photometric and radiometric applications
CN109541560B (zh) 一种基于通讯信号塔的辐射定标方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant