|
关于激光加热的热源加载形式
: _. B# k _- N' \9 R; ^( _/ w7 W) X# z @0 y+ \
看了版上一些讨论关于如何加载激光热源的帖子,觉得好多人做计算时, 都只求获得加载的形式,而不求理解其物理意义。其实如果从物理意义出发,得到加载热源的形式是很容易的。. r' N: K! {7 |6 ^' v9 @, J
2 x0 k! X7 s, `$ A' _例如在平板上加载一个空间上呈高斯分布的激光热源,其功率随时间而变化,
$ h/ ` z( Q$ |2 y4 ^0 M则需要以下参数:
) Z3 b/ l3 r1 [1。 激光功率P[t]。 单位: W , Z( ?2 w% L2 w. [
2。 激光半径r, 单位: m 9 H; V0 ^+ ]8 g+ b i* ~8 w p- |
3。表面反射率R。单位: 无 6 Q' ~( i$ e# _7 w/ K2 V8 T5 T
4。optical absorption coefficient a 单位: m^-1, 或 optical penetration depth:a^-1 单位: m g( x- B# e8 o* f4 | H* R
8 I+ |3 ~5 V2 V7 p9 B激光热源可以用表面热流(W/m^2) 或体热源 (W/m^3)来表示,关键取决与平板厚度d和optical penetration depth a^-1 的比较。( _0 G5 ]$ F2 \7 q) {
: E; v" E7 r9 Z5 u1 D' T Q1。表面热流形式 (d>>a^-1)
) _$ X) k( o3 u- ]/ |) A& H
) L6 s- I; E0 KQ(x,y,t)= P[t]*(1-R)/(Pi*r^2)*exp(-(x^2+y^2)/r^2) / {. u$ s1 |, ~9 \9 d9 m4 x
1 E/ n$ c( a' S0 J* n6 P* y" K
其中P[t]*(1-R)是真正被平板吸收热激光功率。 exp(-(x^2+y^2)/r^2)代表空间高斯分布, 有些文献中会出现exp(-2*(x^2+y^2)/r^2),这是由于高斯半径定义为激光光强1/e和1/e^2的区别。2 K8 H# S& P# h5 f
! G) F+ y$ N5 p4 C/ |- E, ]5 s2。体热源形式 (d~a^-1)
; r- _+ Q$ l8 q2 Y: K1 h
+ t. {: [4 q f; _' vQ(x,y,z,t)= P[t]*(1-R)*a/(Pi*r^2)*exp(-(x^2+y^2)/r^2)*exp(-a*z)
/ _1 z8 H% l% H. J: {" l
! g& w& J: t$ @7 o( l; [1 t y可以看出和表面热流形式相比,多出一项exp(-a*z)代表光在物体中的吸收。/ a8 \+ ]; a: A$ N( P% y5 F- _9 [
另外前面也要多乘个系数a,所以单位是W/m^3.
. U/ q( E; a2 N; w0 d; v- d2 Y
以上的公式其实可以从求解光在物体中传播的Maxwell方程得到,因为光也是电磁波。
- W0 G# Z8 W; G2 f) k/ N9 @. P; r0 S+ I; Y, F" G8 K* y
当 (d>>a^-1)时,表面热流形式和体热源形式对于解传热方程是一样的。
6 @9 T, z8 P8 E( w9 C对于金属,a^-1一般在10~50nm之间。
0 f/ `5 j. t8 M, M% T" b, ^. _9 c; Q1 X# n0 k
对于其他空间分布的激光光源,只用替换后面的exp形式。 |
|