CN107627024B - 一种激光加热悬浮小球的方法及能量分布计算方法 - Google Patents
一种激光加热悬浮小球的方法及能量分布计算方法 Download PDFInfo
- Publication number
- CN107627024B CN107627024B CN201610605215.5A CN201610605215A CN107627024B CN 107627024 B CN107627024 B CN 107627024B CN 201610605215 A CN201610605215 A CN 201610605215A CN 107627024 B CN107627024 B CN 107627024B
- Authority
- CN
- China
- Prior art keywords
- laser
- energy
- bead
- annulus
- follows
- 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
Abstract
本发明涉及一种激光加热悬浮小球的方法及能量分布计算方法,其特征在于:将悬浮小球置于正多面体的中心,多个激光源布置于正多面体的角点,多束激光通过发射相同功率的激光对小球进行加热,直至小球融化。
Description
技术领域
本发明涉及一种利用激光加热的方法,尤其涉及一种利用激光加热悬浮样品以获得温度均匀分布的一种方法以及一种激光加热悬浮小球的能量分布计算方法。
背景技术
空间材料科学,作为空间科学与应用领域中的重要分支,是传统的材料科学向空间环境下的延伸,是发展材料科学新理论、探索材料制备工艺和拓展材料应用新领域中最活跃的前沿性交叉学科之一。空间具有微重力、超真空、无容器和强辐射等特殊效应,是研究材料熔融、凝固等的理想试验条件。但是,人们迄今所能利用的空间资源仍然十分有限。因此,模拟空间环境中各种效应的地面方法应用而生。悬浮技术就是其中之一,它可以模拟空间环境中的无容器状态。在悬浮技术中,通过对悬浮的材料进行激光加热是使材料获得熔融状态的最有效的途径。
研究发现,在对悬浮材料进行熔融加工的过程中,使材料获得均匀的温度分布是进行精确的材料研究的前提条件之一。然后现有技术中并没有很好的方法能够使得材料获得均匀温度。本发明目的在于通过一种新的激光加热方法加热悬浮材料以使小球表面获得均匀的温度分布。
发明内容
本发明是鉴于上述的情况,提供一种能解决上述问题的激光加热方法,具体而言本发明提供一种一种激光加热悬浮小球的方法,其特征在于:将悬浮小球置于正多面体的中心,多个激光源布置于正多面体的角点。
进一步地,其特征在于:所述激光源的输出功率相同。
进一步地,其特征在于:所述激光源的直径相同,且大于等于小球的直径。
进一步地,其特征在于:所述的小球的形状为圆球形。
进一步地,其特征在于:其特征在于所述的正多面体为正四面体。
进一步地,其特征在于:其特征在于所述的正多面体为正六面、正八面、正十二面或正二十面体。
本发明还提供一种一种激光加热悬浮小球的能量分布计算方法,其采用上面任一方案所述激光加热悬浮小球的方法加热,其特征在于,所述计算步骤为:将每束激光在小球表面的分布进行网格划分,即按等距分为多个圆环,例如分为10个圆环,每一圆环对应的弦心角90°/环数,然后根据激光束能量、激光束直径、小球直径、圆环的面积和吸收率可计算出任一圆环吸收的能量,进而可以算得小球表面的能量分布。
进一步地,其特征在于:对于每一束激光都有一个通过球心的圆面与其垂直,且激光能量在该圆面上是均匀分布的,其能量密度为:
其中I为单束激光的能量,R为激光束的半径;
将照射到球面上的激光进行切向和法向的分解,则法向能量密度:n=msinθ
其中θ为激光束与球面接触点切线的夹角;
球面圆环对应的面积:
ds=2πr2cosθdθ
其中r为小球的半径
假设小球材料的吸收率为α,则每一微圆环表面吸收的能量为:
n×α×ds
这样对于弦心角从θ1变化到θ2的球面圆环来说,其面积为:
s=2πr2(sinθ2-sinθ1)
其吸收的能量为:
然后将每一圆环上吸收的单束激光的能量在圆环面积上进行平均分配,为:
然后将球面上同一区域不同激光束在此的能量密度进行加和,即可得到所有激光束在此的总的能量密度分布。
进一步地,其特征在于:对于每一束激光都有一个通过球心的圆面与其垂直,且激光能量在该圆面上是均匀分布的,其能量密度为:
其中I为单束激光的能量,r为小球的半径,r+0.1为激光束的半径。
将照射到球面上的激光进行切向和法向的分解,则法向能量密度:n=msinθ。
球面圆环对应的面积:
ds=2πr2cosθdθ
假设小球材料的吸收率为α,则每一微圆环表面吸收的能量为:
n×α×ds
这样对于弦心角从θ1变化到θ2的球面圆环来说,其面积为:
s=2πr2(sinθ2-sinθ1)
其吸收的能量为:
进一步地,其特征在于:将每一圆环上吸收的单束激光的能量在圆环面积上进行平均分配,为:
然后只需要将球面上同一区域不同激光束在此的能量密度进行加和即可得到所有激光束在此的总的能量密度分布。
效果:
采用本发明的方案,能够简单方便实现悬浮小球的均匀加热,直至小球融化。
附图说明
图1是激光加热悬浮材料的物理结构示意图。
图2是激光束在小球表面分布的网格划分示意图。
图3是激光束在小球表面的分解示意图。
图4a是直径2mm镍基小球在不同激光束照射下温度随时间的变化关系曲线图。
图4b是直径5mm镍基小球在不同激光束照射下温度随时间的变化关系曲线图。
图5a是直径2mm Al2O3小球在不同激光束照射下温度随时间的变化关系曲线图。
图5b是直径5mm Al2O3小球在不同激光束照射下温度随时间的变化关系曲线图。
图6是直径2mm钨小球在不同激光束照射下温度随时间的变化关系曲线图。
图7是铝基小球吸收系数随温度的变化关系。
具体实施方式
为了使本技术领域人员更好的理解本发明,下面结合附图和实施方法对本发明作进一步的详细描述。
本发明的激光加热球状悬浮材料(小球)的物理结构为,将小球置于一空间正多面体结构的中心,并将激光器放置在空间正多面体的各个角点,参见图1,以正四面体为例,即四个激光源占据正四面体的顶角位置,对发出激光束对位于正四面体中心位置的小球进行垂直加热。
其中,激光源发出的激光束直径可根据需要调整,本发明中激光照射小球表面时,其光束直径设置为大于等于小球直径,优选地,为略大于小球的直径,例如大于小球直径0.2mm,这样一束激光照射小球的最大面积为小球的半球区域。由于在小球表面不同的区域会有不同的光束照射,有些区域会有激光束重叠。
参见图2,为了精确的计算小球表面各点激光束的总功率密度,将每束激光在小球表面的分布进行了网格的划分,可以按等距分为多个圆环,例如分为10个圆环,每一圆环对应的弦心角90°/环数,例如为9°
对于每一束激光来说,都有一个通过球心的圆面与其垂直,且激光能量在该圆面上是均匀分布的,其密度为:
其中I为单束激光的能量,r为小球的半径,r+0.1为激光束的半径(这里激光束直径比小球直径大0.2mm为示例,所以激光束半径以r+0.1表示,但是可以知道,激光束直径可以是其他值)。
将照射到球面上的激光进行切向和法向的分解,如图3所示,则法向能量密度:n=m sinθ。
θ为激光束与球面接触点切线的夹角。
球面微圆环对应的面积:
ds=2πr2cosθdθ
假设小球材料的吸收率为α,则每一微圆环表面吸收的能量为:
n×α×ds
这样对于弦心角从θ1变化到θ2的球面圆环来说,其面积为:
s=2πr2(sinθ2-sinθ1)
其吸收的能量为:
这样就可以求得球面上任一微圆环吸收的能量。对于不同的小球只需知道小球的直径和吸收率以及激光束的能量,就可以获得小球表面的能量分布。
为了计算出小球表面激光重叠区域的能量密度,将每一圆环上吸收的单束激光的能量在圆环面积上进行平均分配,为:
这样只需要将球面上同一区域不同激光束在此的能量密度进行加和即可得到所有激光束在此的总的能量密度分布。
对于图1中的物理模型来说,小球表面激光束的能量分布可以分为三大区域,分别为只有单束激光照射的区域,两束激光重叠的区域,三束激光重叠的区域,如图2所示,每一区域都有相互对称的部分,通过仿真软件可以方便的计算出被圆环分割的每一小块的面积。通过面积与该小块面积上激光能量的分布密度即可求得该块区域上的激光能量。
下面通过实施例和比较例对本发明的方法和效果进行说明:
基于上述模型和计算原理,利用ansys软件对激光对小球的加热效果进行了仿真,并与单束激光和两束激光对称加热的情况进行了对比,其结果如表1所示。其中,加热环境为真空,激光总功率为550W,电光转换效率为0.46,这样对于一束激光的情况,其光功率为253W;对于两束激光的情况,每束激光的能量为126.5W;对于四束激光的情况,每束激光的能量为63.25W。
表1不同激光束照射情况下不同材料不同直径小球的表面最大最小温度及温差
对于吸收系数为0.38、发射系数为0.2、直径为2mm和5mm,熔点在1300℃~1450℃内的镍基合金小球来说,当只有1束、2束和4束激光照射时,其表面的温度分布与时间的关系如图4所示,其表面的最大温差如表1所示。
对于吸收系数为0.44、熔点为2030℃、直径为2mm和5mm的Al2O3小球来说,当只有1束、2束和4束激光照射时,其表面的温度分布与时间的关系如图5所示,其表面的最大温差如表1所示。
对于吸收系数为0.13、发射系数为0.3,熔点为3380℃、直径为2mm的钨小球来说,当只有1束、2束和4束激光照射时,其表面的温度分布与时间的关系如图6所示,其表面的最大温差如表1所示。
由图4、5、6和表1可以看出,随着激光束的数量增加,小球表面最高温度和最低温度的差值逐渐减小,说明温度分布愈加均匀。可见,当采用四束激光以正四面体顶角的方位对小球进行激光加热时可获得均匀的温度分布。
上面方案中是以正四面体为例,其最节省激光源。但是可以知道,正六面、正八面、正十二面、正二十面体就可以,正多面体角点越多则激光综合密度越均匀。并且,上面方案中各激光器的功率、直径相同,可以知道,也可以采用不同功率和/或不同直径的激光源,例如高功率和低功率的激光源交替分布等,只要在整体上保证小球上获得的能量基本均匀即可。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。此外,尽管本说明书中使用了一些特定的术语,但这些术语仅仅是为了方便说明,并不对本发明构成任何限制。
Claims (9)
1.一种激光加热悬浮小球的能量分布计算方法,将悬浮小球置于正多面体的中心,多个激光源布置于正多面体的角点进行加热,将每束激光在小球表面的分布进行网格划分,即按等距分为多个圆环,每一圆环对应的弦心角90°/环数,然后根据激光束能量、激光束直径、小球直径、圆环的面积和吸收率可计算出任一圆环吸收的能量,进而可以算得小球表面的能量分布,其特征在于:对于每一束激光都有一个通过球心的圆面与其垂直,且激光能量在该圆面上是均匀分布的,其能量密度为:
其中I为单束激光的能量,R为激光束的半径;
将照射到球面上的激光进行切向和法向的分解,则法向能量密度:
n=m sinθ
其中θ为激光束与球面接触点切线的夹角;
球面圆环对应的面积:
ds=2πr2cosθdθ
其中r为小球的半径,
假设小球材料的吸收率为α,则每一微圆环表面吸收的能量为:
n×α×ds
这样对于弦心角从θ1变化到θ2的球面圆环来说,其面积为:
s=2πr2(sinθ2-sinθ1)
其吸收的能量为:
将每一圆环上吸收的单束激光的能量在圆环面积上进行平均分配,为:
然后将球面上同一区域不同激光束在此的能量密度进行加和,即可得到所有激光束在此的总的能量密度分布。
2.根据权利要求1所述的能量分布计算方法,其特征在于:所述激光源的直径相同,且大于等于小球的直径。
3.根据权利要求2所述的能量分布计算方法,其特征在于:激光直径大于小球直径0.2mm。
4.一种激光加热悬浮小球的能量分布计算方法,将悬浮小球置于正多面体的中心,多个激光源布置于正多面体的角点进行加热,将每束激光在小球表面的分布进行网格划分,即按等距分为多个圆环,每一圆环对应的弦心角90°/环数,然后根据激光束能量、激光束直径、小球直径、圆环的面积和吸收率可计算出任一圆环吸收的能量,进而可以算得小球表面的能量分布,其特征在于:对于每一束激光都有一个通过球心的圆面与其垂直,且激光能量在该圆面上是均匀分布的,其能量密度为:
其中I为单束激光的能量,r为小球的半径,r+0.1为激光束的半径;
将照射到球面上的激光进行切向和法向的分解,则法向能量密度:
n=m sinθ
其中θ为激光束与球面接触点切线的夹角;
球面圆环对应的面积:
ds=2πr2cosθdθ
假设小球材料的吸收率为α,则每一微圆环表面吸收的能量为:
n×α×ds
这样对于弦心角从θ1变化到θ2的球面圆环来说,其面积为:
s=2πr2(sinθ2-sinθ1)
其吸收的能量为:
将每一圆环上吸收的单束激光的能量在圆环面积上进行平均分配,为:
然后将球面上同一区域不同激光束在此的能量密度进行加和,即可得到所有激光束在此的总的能量密度分布。
5.根据权利要求4所述的能量分布计算方法,其特征在于:所述激光源的直径相同。
6.根据权利要求1或4所述的能量分布计算方法,其特征在于:所述圆环为10个。
7.根据权利要求1或4所述的能量分布计算方法,其特征在于:所述激光源的输出功率相同。
8.根据权利要求1或4所述的能量分布计算方法,其特征在于:所述的小球的形状为圆球形。
9.根据权利要求1或4所述的能量分布计算方法,其特征在于:所述的正多面体为正四面、正六面、正八面、正十二面或正二十面体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610605215.5A CN107627024B (zh) | 2016-07-28 | 2016-07-28 | 一种激光加热悬浮小球的方法及能量分布计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610605215.5A CN107627024B (zh) | 2016-07-28 | 2016-07-28 | 一种激光加热悬浮小球的方法及能量分布计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107627024A CN107627024A (zh) | 2018-01-26 |
CN107627024B true CN107627024B (zh) | 2019-07-26 |
Family
ID=61113394
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610605215.5A Active CN107627024B (zh) | 2016-07-28 | 2016-07-28 | 一种激光加热悬浮小球的方法及能量分布计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107627024B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111521517B (zh) * | 2020-04-10 | 2022-05-10 | 中国科学院上海硅酸盐研究所 | 一种基于双相机视觉的熔融态悬浮椭球液滴图像处理算法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5043131A (en) * | 1989-12-18 | 1991-08-27 | Kms Fusion, Inc. | Ignition of deuterium-trtium fuel targets |
JP2003334677A (ja) * | 2002-03-14 | 2003-11-25 | Onizuka Glass:Kk | レーザ加工装置およびこれを用いた加工方法 |
JP2004103261A (ja) * | 2002-09-04 | 2004-04-02 | Japan Science & Technology Corp | 球状超伝導体及びその製造方法 |
JP2011115819A (ja) * | 2009-12-02 | 2011-06-16 | Hitachi Global Storage Technologies Netherlands Bv | ヘッド・ジンバル・アセンブリの製造方法及びヘッド・ジンバル・アセンブリの接続パッドを相互接続する装置 |
CN105575444A (zh) * | 2016-02-06 | 2016-05-11 | 北京应用物理与计算数学研究所 | 热核聚变发生方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060159394A1 (en) * | 2005-01-19 | 2006-07-20 | Victor Grubsky | Method of and apparatus for manufacturing fiber grating devices |
-
2016
- 2016-07-28 CN CN201610605215.5A patent/CN107627024B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5043131A (en) * | 1989-12-18 | 1991-08-27 | Kms Fusion, Inc. | Ignition of deuterium-trtium fuel targets |
JP2003334677A (ja) * | 2002-03-14 | 2003-11-25 | Onizuka Glass:Kk | レーザ加工装置およびこれを用いた加工方法 |
JP2004103261A (ja) * | 2002-09-04 | 2004-04-02 | Japan Science & Technology Corp | 球状超伝導体及びその製造方法 |
JP2011115819A (ja) * | 2009-12-02 | 2011-06-16 | Hitachi Global Storage Technologies Netherlands Bv | ヘッド・ジンバル・アセンブリの製造方法及びヘッド・ジンバル・アセンブリの接続パッドを相互接続する装置 |
CN105575444A (zh) * | 2016-02-06 | 2016-05-11 | 北京应用物理与计算数学研究所 | 热核聚变发生方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107627024A (zh) | 2018-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yan et al. | Optimization of a discrete dish concentrator for uniform flux distribution on the cavity receiver of solar concentrator system | |
Yorke et al. | The formation of protostellar disks. 2: Disks around intermediate-mass stars | |
Mooney et al. | Star formation rates and the far-infrared luminosity of Galactic molecular clouds | |
Harding et al. | Pulsar pair cascades in a distorted magnetic dipole field | |
Roettiger et al. | Relic radio emission in 3C 388 | |
CN104864278B (zh) | Led自由曲面照明系统 | |
CN107627024B (zh) | 一种激光加热悬浮小球的方法及能量分布计算方法 | |
Gao et al. | Mega-gauss plasma jet creation using a ring of laser beams | |
Witte et al. | Beyond the margin recipe: the probability of correct target dosage and tumor control in the presence of a dose limiting structure | |
Song et al. | A tool for fast flux distribution calculation of parabolic trough solar concentrators | |
Anwar et al. | A thermal battery mimicking a concentrated volumetric solar receiver | |
Grinin et al. | Disk wind in the radiation of young intermediate-mass stars | |
CN103530697A (zh) | 辐射型塔式太阳能热电系统的镜场优化设计方法 | |
Bannikova et al. | Dynamical model of an obscuring clumpy torus in AGNs–I. Velocity and velocity dispersion maps for interpretation of ALMA observations | |
Fukumura et al. | Soft X-Ray Excess from Shocked Accreting Plasma in Active Galactic Nuclei | |
CN111180092B (zh) | 激光排布确定方法、激光装置使用方法、六孔球腔和系统 | |
Kocifaj et al. | Designing of light-pipe diffuser through its computed optical properties: A novel solution technique and some consequences | |
Shcherbakov et al. | Modeling hot gas flow in the low-luminosity active galactic nucleus of NGC 3115 | |
Pavlović et al. | Optical design of a solar parabolic thermal concentrator based on trapezoidal reflective petals | |
US10038107B2 (en) | Enhanced photo-thermal energy conversion | |
Strittmatter et al. | X-Ray Emission from White Dwarfs | |
Law | The effect of differential rotation on the internal structure of low mass stars | |
CN207833104U (zh) | 实现旋转对称三维光强分布自由曲面透镜 | |
Pavlović et al. | Optical design of a solar parabolic concentrating collector based on trapezoidal reflective petals | |
Gus’ kov et al. | Analysis of direct-drive capsule compression experiments on the Iskra-5 laser facility |
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 |