CN115494639B - 高功率激光光束合成***内通道热效应仿真方法 - Google Patents
高功率激光光束合成***内通道热效应仿真方法 Download PDFInfo
- Publication number
- CN115494639B CN115494639B CN202211373796.6A CN202211373796A CN115494639B CN 115494639 B CN115494639 B CN 115494639B CN 202211373796 A CN202211373796 A CN 202211373796A CN 115494639 B CN115494639 B CN 115494639B
- Authority
- CN
- China
- Prior art keywords
- channel
- distribution
- sub
- field
- optical
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Lenses (AREA)
Abstract
本发明提供了一种高功率激光光束合成***内通道热效应仿真方法,通过先利用流体动力学的方法计算瞬态情况下流场的密度分布,并将流场的密度分布映射为折射率场分布后,以各子光束光路的光路转折处为界划分为多段直通道,对每段直通道引入傍轴近似,在不考虑光线对初始传播方向的偏离的情况下,分别计算各段直通道由于折射率分布不均引起的光程差,并对各段直通道的光程差进行叠加,得到各子光束在出光口处的光程差及波前相位分布,再进一步计算各子光束远场光斑分布,最终计算出合成光束的光束质量因子。通过在***完成光机设计后对其内通道热效应的影响进行精确仿真,获得激光束波前畸变情况,从而可利用上述结果指导光机结构的调整。
Description
技术领域
本发明涉及仿真技术技术领域,尤其涉及一种高功率激光光束合成***内通道热效应仿真方法。
背景技术
随着工业制造和国防安全等领域对激光***的输出功率要求越来越高,单路激光的功率已经不能满足相关应用需求,为实现高功率激光输出,通用的方法是基于各类合束技术的激光光束合成***将多束高功率激光合为一束。高功率激光光束合成***中从激光束入射端到***出口端这段光传输通道称为***内通道,高功率激光在内通道中经过时会加热气体、机械结构件及光学元件,通过热传导及热对流使气体温度升高,造成激光传输通道上气体流场密度不均匀,引起激光波前畸变,降低激光传输质量。
目前,在内通道热效应仿真方面,已有研究人员对内通道介质气体热效应仿真方法进行了研究。现有技术中,期刊《光电子激光》,2004年1月,第15卷第1期,作者柳建等,名称为“封闭充气管道中气体热效应对激光传输的影响”中采用稳相法求解光波传输方程,计算了激光束在封闭管道传输过程中热效应对光束质量的影响。再者,期刊《封闭管道中热耦合效应对激光传输的影响》,2013年7月,第25卷第7期,作者胡鹏等,名称为“封闭管道中热耦合效应对激光传输的影响”中建立了封闭竖直管道热耦合效应物理模型,采用流固耦合方法,数值模拟了封闭管道中介质气体、光学玻璃之间的热相互作用造成流场分布的变化,研究了热耦合作用对近场波前相差和远场光束质量的影响。上述技术方案均是对光通道的热效应仿真,但仅适用于单直通道等简单的光通道,且光束波前的计算算法复杂,因而不适用于复杂光路的仿真。
有鉴于此,有必要设计一种改进的高功率激光光束合成***内通道热效应仿真方法,以解决上述问题。
发明内容
本发明的目的在于提供一种高功率激光光束合成***内通道热效应仿真方法。
为实现上述发明目的,本发明提供了一种高功率激光光束合成***内通道热效应仿真方法,所述高功率激光光束合成***内通道热效应仿真方法包括如下步骤:
S1、依据光束合成***光学设计结果与光机结构模型确定光束合成***光传输内通道的基本三维结构,再将所属基本三维结构模型离散化为有限元模型;在加载激光束等效热源分布和设定边界条件情况下,利用所述有限元模型进行流体动力学仿真,得到光传输内通道有限元模型每个网格节点的温度以及密度场分布;然后,将所述密度场分布转换为折射率场分布;
S2、以各子光束光通道的光通道转折处为界将所述子光束光通道划分为多段直通道,以步骤S1中得到的所述折射率场分布计算各子光束内一定数量的光线通过每段变折射率直通道后引入的光程差,通过插值得到该段直通道引入的光程差分布,对各段直通道的光程差进行叠加,得到各子光束在出光口处的光程差和波前相位分布;
S3、利用步骤S2得到的所述各子光束在出光口处的光程差分布,计算所述各子光束的远场光斑;将所有所述子光束的远场光斑叠加,得到合成光束的远场光斑分布;再根据所述合成光束远场光斑分布计算合成光束的光束质量。
优选的,所述激光光束合成***包括改变子光束的传播方向的光学器件和用于合成光束的衍射光学元件,光束合成过程中所述子光束在所述光学器件的作用下汇聚至所述衍射光学元件的表面,所述衍射光学元件再将所有所述子光束合成为一束。
优选的,在步骤S2中,所述光程差的计算在引入傍轴近似的情况进行。
优选的,在步骤S3中,所述各子光束的远场光斑为根据所述光程差计算近场复振幅分布后,再对所述近场复振幅分布进行二维傅里叶变换计算得到。
本发明的有益效果是:
1、本发明提出的高功率激光光束合成***内通道热效应仿真方法,通过先利用流体动力学的方法计算瞬态情况下流场的密度分布,并将流场的密度分布映射为折射率场分布后,以各子光束光路的光路转折处为界划分为多段直通道,对每段直通道引入傍轴近似,在不考虑光线对初始传播方向的偏离的情况下,分别计算各段直通道由于折射率分布不均引起的光程差,并对各段直通道的波像差进行叠加,得到各子光束在出光口处的波像差及波前相位分布,并进一步计算各子光束远场光斑分布,最终计算出合成光束的光束质量因子;上述过程实现了对高功率激光光束合成***合成光束的远场光斑分布、光束质量的仿真,并能够根据质量因子的数值判断合成光束的质量,为光机结构的设计提供了理论基础。本发明中涉及到的各子光束的波像差仿真具有计算方法简单,仿真通道选取灵活的优点。
2、本发明提出的高功率激光光束合成***内通道热效应仿真方法,通过在***完成光机设计后对其内通道热效应的影响进行精确仿真,获得激光束波前畸变情况,基于上述结果可进一步对光机的结构进行调整,从而实现在光机的设计阶段有效提升光束质量的目的;本发明提出的仿真方法不仅能够满足任意合成技术路径、任意光路结构、任意多子光束路数光束合成***输出光束的光束质量仿真,且适合于简单光通道和光路复杂的高功率激光合成***,具有实用性强的优势。
附图说明
图1为本发明的高功率激光光束合成***内通道热效应仿真方法的仿真计算流程图;
图2为实施例1中的光谱合成***光路示意图;
图3为实施例2中基于DOE的相干合成***光路示意图;
附图标记如下:
1、反射镜;2、衍射光栅;3、透镜;4、衍射光学元件。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面结合附图和具体实施例对本发明进行详细描述。
在此,还需要说明的是,为了避免因不必要的细节而模糊了本发明,在附图中仅仅示出了与本发明的方案密切相关的结构和/或处理步骤,而省略了与本发明关系不大的其他细节。
另外,还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
请参阅图1,本发明提出了一种高功率激光光束合成***内通道热效应仿真方法,激光光束合成***包括改变光束的传播方向的光学器件和用于合成光束的衍射光学元件,光束合成过程中子光束在光学器件的作用下汇聚至衍射光学元件的表面,衍射光学元件再将所有子光束合成为一束。
具体地,对高功率激光光束合成***内通道热效应的仿真方法,包括如下步骤:
S1、合成***内通道三维变折射率场的构建:
首先依据光束合成***光学设计结果与光机结构模型确定光束合成***光传输内通道的基本三维结构,再将三维结构模型离散化为有限元模型;在加载激光束等效热源分布和设定边界条件情况下,利用有限元模型进行流体动力学仿真,得到光传输内通道有限元模型每个网格节点的温度以及密度场分布;然后,将内通道密度场分布转换为折射率场分布;
S2、合成***各子光束在出光口处波像差分布的计算:
在步骤S1得到合成***内通道折射率场分布后,针对每路子光束光通道,以光通道转折处为界划分为多段直通道;接着,引入傍轴近似,在不考虑光线对初始传播方向偏离的情况下,计算各子光束内一定数量的光线通过每段变折射率直通道后引入的光程差,通过插值得到该段直通道引入的光程差分布,对各段直通道的光程差进行叠加后,得到各子光束在出光口处的光程差和波前相位分布;
S3、合成光束的光束质量的计算:
利用傅里叶转换将步骤S2得到的各子光束在出光口处的波像差分布,转化为各子光束的远场光斑;将所有子光束的远场光斑叠加,即得到合成光束的远场光斑分布;再根据合成光束远场光斑分布计算合成光束的光束质量;需要说明的是,上述过程中的为空间坐标位置。
下面结合具体的实施例对本发明的高功率激光光束合成***内通道热效应仿真方法作进一步说明:
实施例1
本实施例对反射式光谱合成***光路进行内通道热效应仿真,反射式光谱合成***光路的示意图如图2所示,光束合成的过程按照如下方式进行:具有一定光焦度的反射镜1将具有不同角度入射和波长的子光束汇聚在衍射光栅2的表面,衍射光栅2将所有子光束合成为一束,针对上述激光合成***内通道的热效应仿真按照如下步骤进行:
S1、合成***内通道三维变折射率场的构建:
首先依据光束合成***光学设计结果与光机结构模型确定光束合成***光传输内通道的基本三维结构,再将三维结构模型离散化为有限元模型;对光传输内通道有限元模型在加载激光束等效热源分布并设定边界条件情况下进行计算流体动力学仿真,得到光传输内通道有限元模型每个网格节点的温度以及密度场分布;然后,将内通道密度场分布转换为折射率场分布;其中,密度场分布与折射率场分布之间的转换按照如下关系式进行:
S2、合成***各子光束在出光口处波像差分布的计算:
在步骤S1得到合成***内通道折射率场分布后,针对每路子光束光通道,以光通道转折处为界划分为多段直通道;接着,引入傍轴近似,在不考虑光线对初始传播方向偏离的情况下,计算各子光束内一定数量的光线通过每段变折射率直通道后引入的光程差,通过插值得到该段直通道引入的光程差分布,对各段直通道的光程差进行叠加,得到各子光束在出光口处的光程差和波前相位分布;其中,光程差的计算式为:,的计算式为:,表示光束在通道中的任意时刻的折射率分布,表示空气密度分布未受扰动时的光程长度;
S3、合成光束光束质量的计算:
利用步骤S2得到的各子光束在出光口处的波像差分布,计算各子光束的远场光斑;将所有子光束的远场光斑叠加,得到合成光束的远场光斑分布;再根据合成光束远场光斑分布计算合成光束光束质量;其中,子光束的远场光斑计算式为:,合成光束的光束质量的具体计算式为:,其中通过计算得到。
实施例2
本实施例对透射式基于衍射光学元件的相干合成***光路进行内通道热效应仿真,相干合成***光路的示意图如图3所示,光束合成的过程按照如下方式进行:利用透镜3将具有不同入射角和波长相同的子光束汇聚在衍射光学元件4的表面,衍射光学元件4将所有子光束合成为一束;本实施例中针对上述激光合成***内通道的热效应仿真方法与实施例1基本一致,在此不再赘述。
综上所述,本发明提出的高功率激光光束合成***内通道热效应仿真方法,通过利用流体动力学的方法实现了高功率激光光束合成***内通道热效应仿真,应用上述方法可对光机内通道热效应的影响进行精确仿真,以获得激光束波前畸变情况,进而采取相应措施对光束质量进行控制。通过上述方式提供了一种计算方法简单,且适用于简单光通道和光路复杂的非单一光通道的高功率激光合成***。
以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。
Claims (8)
1.一种高功率激光光束合成***内通道热效应仿真方法,其特征在于,所述激光光束合成***包括改变子光束的传播方向的光学器件和用于合成光束的衍射光学元件,光束合成过程中所述子光束在所述光学器件的作用下汇聚至所述衍射光学元件的表面,所述衍射光学元件再将所有所述子光束合成为一束;所述高功率激光光束合成***内通道热效应仿真方法包括如下步骤:
S1、依据光束合成***光学设计结果与光机结构模型确定光束合成***光传输内通道的基本三维结构,再将所属基本三维结构模型离散化为有限元模型;在加载激光束等效热源分布和设定边界条件情况下,利用所述有限元模型进行流体动力学仿真,得到光传输内通道有限元模型每个网格节点的温度以及密度场分布;然后,将所述密度场分布转换为折射率场分布;
S2、以各子光束光通道的光通道转折处为界将所述子光束光通道划分为多段直通道,以步骤S1中得到的所述折射率场分布计算各子光束内一定数量的光线通过每段变折射率直通道后引入的光程差,通过插值得到该段直通道引入的光程差分布,对各段直通道的光程差进行叠加,得到各子光束在出光口处的光程差和波前相位分布;其中,在步骤S2中,所述密度场分布与折射率场分布之间的转换按照如下关系式进行:,是Gladstone Dale常数;
S3、利用步骤S2得到的所述各子光束在出光口处的光程差分布,计算所述各子光束的远场光斑;将所有所述子光束的远场光斑叠加,得到合成光束的远场光斑分布;再根据所述合成光束远场光斑分布计算合成光束的光束质量。
6.根据权利要求1所述的高功率激光光束合成***内通道热效应仿真方法,其特征在于,在步骤S2中,所述光程差的计算在引入傍轴近似的情况进行。
8.根据权利要求1所述的高功率激光光束合成***内通道热效应仿真方法,其特征在于,在步骤S3中,所述各子光束的远场光斑为根据所述光程差计算近场复振幅分布后,再对所述近场复振幅分布进行二维傅里叶变换计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211373796.6A CN115494639B (zh) | 2022-11-04 | 2022-11-04 | 高功率激光光束合成***内通道热效应仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211373796.6A CN115494639B (zh) | 2022-11-04 | 2022-11-04 | 高功率激光光束合成***内通道热效应仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115494639A CN115494639A (zh) | 2022-12-20 |
CN115494639B true CN115494639B (zh) | 2023-02-17 |
Family
ID=85115490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211373796.6A Active CN115494639B (zh) | 2022-11-04 | 2022-11-04 | 高功率激光光束合成***内通道热效应仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115494639B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117787021B (zh) * | 2024-02-28 | 2024-05-07 | 中国人民解放军海军工程大学 | 激光远场能量密度估计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382665A (zh) * | 2008-10-24 | 2009-03-11 | 中国科学院上海光学精密机械研究所 | 相干光束同轴合成的方法 |
CN104793339A (zh) * | 2015-05-13 | 2015-07-22 | 中国人民解放军国防科学技术大学 | 一种基于自适应精密光程调节的光纤激光相干合成*** |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5052017A (en) * | 1988-12-01 | 1991-09-24 | Coherent, Inc. | High power laser with focusing mirror sets |
CA2260565A1 (en) * | 1999-01-26 | 2000-07-26 | Cymbolic Sciences International Inc. | Method for combining the output of a laser diode array into a single telecentric stripe |
WO2008089953A1 (en) * | 2007-01-22 | 2008-07-31 | Carl Zeiss Smt Ag | Method for improving imaging properties of an optical system, and optical system |
CN107832532A (zh) * | 2017-11-14 | 2018-03-23 | 哈尔滨工业大学 | 一种气动光传输效应对高速飞行器成像质量数值计算方法 |
US11042025B2 (en) * | 2019-09-20 | 2021-06-22 | Raytheon Company | Optical data communication using micro-electro-mechanical system (MEMS) micro-mirror arrays |
JP2023548820A (ja) * | 2020-10-27 | 2023-11-21 | アイピージー フォトニクス コーポレーション | 体積ブラッググレーティングにより広帯域レーザビームをスペクトル結合するための装置 |
CN113341580B (zh) * | 2021-05-08 | 2022-06-07 | 西安电子科技大学 | 一种相干激光合成*** |
-
2022
- 2022-11-04 CN CN202211373796.6A patent/CN115494639B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382665A (zh) * | 2008-10-24 | 2009-03-11 | 中国科学院上海光学精密机械研究所 | 相干光束同轴合成的方法 |
CN104793339A (zh) * | 2015-05-13 | 2015-07-22 | 中国人民解放军国防科学技术大学 | 一种基于自适应精密光程调节的光纤激光相干合成*** |
Non-Patent Citations (2)
Title |
---|
激光通道传输热特性对远场光束质量的影响;刘顺发等;《强激光与粒子束》;20040630;第703-705页 * |
高功率激光热效应对合束***的影响;孙毅 等;《光学精密工程》;20151130;第3098-3105页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115494639A (zh) | 2022-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115494639B (zh) | 高功率激光光束合成***内通道热效应仿真方法 | |
Lukin | Adaptive optics in the formation of optical beams and images | |
CN113820857B (zh) | 一种产生完美平顶光束/平顶涡旋光束的方法 | |
CN114755837B (zh) | 一种全庞加莱球偏振阵列光束的产生方法及装置 | |
Ahn et al. | Numerical simulation of high-energy laser propagation through the atmosphere and phase correction based on adaptive optics | |
Muslimov et al. | Spectrograph with a composite holographic dispersive element | |
Gavril'eva et al. | Approaches to cross-talk noise reduction in modal holographic wavefront sensors | |
CN115774343A (zh) | 新型数字化艾里类光束的产生与调控装置及方法 | |
Korolkov et al. | New methods of manufacturing high-aperture computer-generated holograms for reference wavefront shaping in interferometry | |
US4850662A (en) | HOE and indirect method of constructing same | |
JP6757307B2 (ja) | 回折素子の設計方法 | |
CN103513425A (zh) | 基于全息光学元件的机载激光雷达对地观测照明成像光照匀化*** | |
CN115993696B (zh) | 一种焦斑控制方法及其*** | |
Akhmetov et al. | Modeling and optimization of optical designs with composite holographic elements | |
Krasin et al. | Measurement of wavefront curvature using computer-generated Fourier holograms | |
Gorbunov et al. | Method of spherical phase screens for modeling the propagation of diverging beams in inhomogeneous media | |
Kashef et al. | Enhancement of thermal blooming effect on free space propagation of high power CW laser beam | |
Shtanko et al. | Compensation of aberrations of off-axis parabolic mirrors by means of digital holography | |
Zhai et al. | High uniformity Bessel beam generated by the axicon with a high-order curved surface | |
Lukin et al. | Four-dimensional computer dynamic model of an atmospheric optical system | |
TW202316094A (zh) | 光學計算方法 | |
Yang et al. | Optical terminal analysis of a multigrating tiled compressor in a PW-class CPA-laser | |
Shih et al. | Metasurface hybrid lens inverse design with diffractive/refractive co-optimization | |
Roggemann et al. | Two deformable mirror concept for correcting scintillation effects in laser beam projection | |
Furieri | Extended Field of View using Multi Conjugated Deformable Lenses |
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 |