激光加載過程中金屬熔化的數(shù)值計算——論文開題答辯ppt

1、激光加載過程中金屬熔化的數(shù)值計算,年級：2010級姓名：嚴向峰學號：201005089指導教師：徐讓書,報告大綱,課題背景國內外科研現(xiàn)狀課題目標建模分析算例結果分析結論,一.課題背景,激光加載過程中壁面溫度場以及相變過程取決于壁面加載條件，激光輻照所引起的質量遷移也同入射表面的環(huán)境條件、附面層流場狀況等因素密切相關。亞音速空中飛行物近壁區(qū)速度附面層和溫度附面層結構復雜，附面層中輸運現(xiàn)象旺盛，與激光輻照作用中的傳熱、傳質

2、過程以及多相流強烈耦合，從而可能產(chǎn)生氣動加熱，氣動冷卻以及固體靶材熔化后所產(chǎn)生的流固耦合等諸多現(xiàn)象，對激光對壁面材料和結構的加熱、熔化、氣化、力學效應以及等離子現(xiàn)象產(chǎn)生很大影響。 課題研究亞音速空中目標附面層流場，以及近壁區(qū)中動量、能量和質量輸運對壁面激光加載特性的影響，揭示在有氣流條件下激光輻照效應機制，著重研究壁面已存在金屬熔池狀態(tài)下的兩相流流動，對于評估目標抗激光防護技術在實際飛行環(huán)境中的效能有重要意義。,

3、二.國內外科研現(xiàn)狀,在激光對目標物壁面加載特性方面，國內開展了激光作用下物體的溫度場、熔化和燒蝕的仿真計算和實驗研究。流場對加載特性的影響：,李海燕等人以美國ABL系統(tǒng)激光武器為例，分析了高超聲速飛行器附面層流場氣體影響激光武器毀傷效應的一些主要因素。試件附面層流場中的氣體具備大氣層氣體的所有特性，激光在流場中傳輸時必然經(jīng)歷與在大氣中傳輸時的相同過程。近壁面氣流流動狀況與溫度條件決定了附面層流場中動量、能量和質量的輸運過程，對于靶材

4、發(fā)生表面熔化后與氣流相互作用所形成的熔池流動、熔融物質運動、靶材材料氣化擴散、固壁內的溫度場分布等都有較大影響。,,王偉平等研究了鋼片在空氣自然對流及低速（Ma = 0.1）切向氣流這兩類不同試驗環(huán)境下使用重復頻率YAG激光輻照下的溫度響應。通過理論分析激光作用下的氣流與靶面的換熱通量，發(fā)現(xiàn)其與激光參數(shù)（強度、作用時間）、氣流參數(shù)（物理性質，流速）及靶參數(shù)（物理性質，幾何尺寸）都有關系。實驗結果表明，在激光作用下，靶材的溫升與靶前表面的

5、對流換熱狀況有較大關系，當靶面存在低速附面層流場時，靶材的溫升明顯低于其處于空氣自然對流環(huán)境下所采集數(shù)據(jù)，證實了氣動冷卻效應對激光加載特性的影響。,右圖為YAG激光加載時固壁溫度場隨時間的變化，能明顯觀測到氣動冷卻效果,,杜太焦等研究了高速目標在激光輻照下的溫升效應，采用來自風洞實驗的準則方程式作為邊界上的對流換熱系數(shù)，并在計算中考慮到了材料參數(shù)隨溫度的變化。陳敏孫等采用有限差分法對激光輻照下復合材料壁面非穩(wěn)態(tài)一維傳熱模型進行數(shù)值求解，

6、考慮到了附面層氣流流場對于能量加載過程的影響。但由于算法及計算側重點的限制，切向氣流對流換熱以及熱分解氣體對表面熱交換覆蓋效應的影響僅以邊界條件的形式代入計算。鄭艷麗等采用有限元方法求解鋼板在第二類和第三類邊界條件下的導熱非穩(wěn)態(tài)泊松方程，并考慮材料氧化反應熱，得到在吹空氣、吹氮氣以及無吹氣條件下靶材的溫度場時空分布。 張健的數(shù)值計算研究結果比以上研究更接近于考慮流動與輻照效應的耦合，他利用有限元方法分析了亞聲速條件下運動目

7、標在激光輻照全過程的溫度場和熱應力場的分布與演化規(guī)律。由于高速流場的存在，計算結果顯示出了明顯的氣動冷卻效果，并且輻照區(qū)域的兩側由于受熱程度不同而出現(xiàn)了熱應力分布不均的現(xiàn)象。但是其采用的計算方案將流場處理為簡化的二維不可壓縮流動模型，因而未能反映流場特別是附面層結構對壁面激光加載特性影響的規(guī)律和機理。,,在激光作用下金屬熔化和燒蝕問題中，理論研究、實驗研究和數(shù)值計算仿真的研究成果加深了對相關現(xiàn)象和機制的理解。激光加載下金屬的熔化與燒

8、蝕：,楊立新等采用有限元方法研究了激光加工過程中熔池熔化和凝固的過程，與之前的研究有所不同的是，通過在計算中引入液相分數(shù)標量，以熱焓源項導入熔化潛熱，應用“Darcy”源項導入動量方程，從而建立起統(tǒng)一的液體和固體連續(xù)方程。規(guī)避了研究熔池特性過程中視熔化邊界為移動邊界條件，將流體區(qū)與固體區(qū)分開單獨求解這種非最優(yōu)方案。通過比較不同機理主導作用下的熔池形狀，分析了熱傳導、Marangoni和熱浮力三種不同的換熱機理對于計算結果的影響。模擬的結

9、果顯示，在多重換熱機理作用下，熔池內部換熱更加活躍，熔池形狀將變得更寬更淺，與只存在熱傳導機制下的熔池形狀有明顯差別。表面張力在整個熔池形成，發(fā)展以至最后停止加載后的凝固階段都起著非常重要的作用。,,綜合換熱機理作用下熔池內部流動 綜合換熱機理作用下熔池形狀差異,當存在金屬熔池時，熱傳導、熱浮力與熱表面張力的綜合作用使得熔池內部出現(xiàn)旋流運動，這將使激光加載能量更快的傳遞至熔

10、池邊緣，使得熔池形狀更寬更淺,,汪任憑等建立了激光深熔焊接過程中氣、液、固三相統(tǒng)一的數(shù)學模型，并采用了VOF方法追蹤自由液面，模型中考慮了表面張力、熱浮力和反沖壓力的作用，通過數(shù)值計算得到了匙孔變化的動態(tài)過程與相應的溫度分布。在加熱初期，激光能量主要用于加熱金屬基體，熱傳導是這一階段熱量傳遞的主要方式。隨著加載過程的進行，靶材表面經(jīng)歷熔化至汽化的過程，在較大的蒸汽反沖作用下，在熔化的金屬內部形成匙孔。匙孔的存在有利于金屬材料對激光能量的

11、吸收。,,目前，國外公開發(fā)表的相關文獻中鮮見關于空中飛行器附面層流場對激光武器毀傷效應的直接研究，對流場模擬裝置及其流場的研究報道也較少。,激光加載致使金屬熔化這一領域，多位研究者均使用數(shù)值模型的方式來模擬激光加工過程中的熔池特性。Basu B等利用CFD方法計算了穩(wěn)態(tài)情況下激光加載熔化問題，在仿真過程中將熔化界面處理成為移動邊界，對此邊界兩側的固/液相區(qū)域分開求解，網(wǎng)格重新劃分的工作量較大。 且這樣的處理手段還將導致仿真結果對初始條件

12、的巨大依賴，Lee S Y等人在隨后的數(shù)值仿真模擬中發(fā)現(xiàn)，流動區(qū)初始熔池形狀的不同選擇將帶來明顯差異的仿真結果。,,分析上述國內外研究成果，不難發(fā)現(xiàn)以下不足：大多數(shù)的研究中，附面層流場對于壁面激光加載特性的影響僅進行至壁面材料熔化之前。而涉及到熔化之后的研究通常為激光加工及工藝范疇，不考慮切向氣流的影響，罕有同時考慮兩者共同作用時的研究成果。模型建立過程中，對激光輻照能量的加載方式進行了簡化處理，大多采用邊界換熱條件引入計算，或

13、采用更為簡單的固定空間體熱源方式，不能夠準確的反映出激光加載的特點。在壁面金屬熔池形成與發(fā)展的計算過程中，采用的簡化手段使得計算結果嚴重依賴于初始條件。計算時大多將熔化邊界的發(fā)展簡單的看做由金屬內部熱量傳導所引起，不考慮內部沸騰及平衡汽化等現(xiàn)象。,三.課題目標,實現(xiàn)有附面層氣流條件下飛行目標物的激光輻照熱力效應的耦合數(shù)值計算，著重研究靶材表面已存在金屬熔池情況下激光加載時的兩相流流場以及溫度場，并在后續(xù)計算中考慮相變及化學反應等因素

14、。以下是本次課題所研究的主要內容： 1.亞音速附面層流場中激光加載金屬壁面熔化/凝固相變過程計算； 2.激光輻照能量加載過程中亞音速氣流掠過金屬熔池時的多相流計算； 3.計算過程中將考慮平衡汽化效應的質量輸運現(xiàn)象對仿真結果的影響。,四.建模分析,右圖為本次課題所采用的初始計算域，包含近壁流場，噴嘴以及“無限遠”流場等三部分的大計算域。大計算域的網(wǎng)格數(shù)量巨大，對計算資源的消耗極其可觀，無法用于激光能量加載的

15、直接計算。但可以在大計算域中進行附面層流場的計算，進而獲取精確的流場參數(shù)作為小計算域的流動邊界條件,,,模型的精簡,在對模型進行能量加載計算之前，必須對三維模型進行簡化。簡化后的計算域將只包含靶材平面以及其上方部分氣相（Air）空間。小計算域模型如右圖。為了盡可能的降低模型對計算資源的消耗，在試算之后對固相域的網(wǎng)格數(shù)量進行了進一步的削減。這一操作的根據(jù)及目的將在其后說明。,,滑移壁面,壓力進口,壓力出口,氣相域,固相域,輻照表面,物性

16、參數(shù)設置,模型中選擇的靶材為鋼鐵，在選擇其物性參數(shù)時，部分參數(shù)根據(jù)鋼鐵中含碳量的不同略有差異，不過這并不影響對于激光在存在附面層條件下壁面加載機理的研究。靶面材料物性參數(shù)如下表：,模型中流場風速為0.3Ma，在大計算域中采集某一截面處的流場數(shù)據(jù)，使用profile方式將其導入精簡后的計算域的壓力進口。 在精簡模型中，僅保留了激光光斑所在區(qū)域的16×16×10（mm）固相區(qū)域，其區(qū)域邊界距離光斑邊緣最近距

17、離為3mm，氣相區(qū)域仍保持30×30×20（mm）原有體積不變。,,,根據(jù)計算過程中的溫度變化范圍，空氣及靶材材料隨溫度變化的物性參數(shù)如下,,,在表中可以看到，金屬類材料隨著自身溫度的上升，其導熱性將不斷下降。從前人的研究中得知，利用Q600SDT差式掃描量熱儀對HT62銅鋅合金的比熱容進行測定，能夠顯著觀察到材料導熱系數(shù)隨溫度上升而下降趨勢，并在材料熔化后能夠觀測到階躍式下降。 本次課題中在靶材物性

18、參數(shù)的設置中能夠體現(xiàn)出靶材導熱系數(shù)隨溫度的上升而下降，暫時沒有考慮熔化后導熱性能突降的特性。在FLUENT設置上述隨溫度變化的物性參數(shù)時，均采用分段線性插值方式導入計算模型。,激光加載的計算與分析,,,非傅立葉效應,,,UDF（用戶自定義函數(shù)）分析,,激光能量具有指向性強這一突出特點，但FLUENT的輻射模型中對于輻射的離散處理無法滿足這一要求，因此需要采用其他手段將此類定向能（激光，電子束）加載于模型之上。在之前的研究中，大多數(shù)采用將

19、輻射能轉化為熱流邊界條件的方式，加載于計算模型之中。這樣的處理方式，在出現(xiàn)移動邊界之后顯然是不適用的。　　在課題研究中，采用UDF方法來實現(xiàn)激光能量的加載。由于激光在金屬中的趨膚深度很小，對紅外波段吸收層厚度也只有幾十微米，因此可以認定：激光能量在氣液/氣固相界面處即完成了吸收與反射的全過程。根據(jù)這一物理特點編寫UDF，將激光能量以內熱源形式添加于兩相界面所在的網(wǎng)格之中?！　≡谝淹瓿傻腢DF中，將實現(xiàn)激光能量在相界面處的動態(tài)加載，滿

20、足能量載荷的高斯分布要求，能夠模擬平衡氣化過程中出現(xiàn)的質量輸運過程。,,,,能量加載示意圖,UDF設計中的注意事項,,由于在UDF中采用的相界面判斷準則為相體積分數(shù)0 < αq ≤1，因此在出現(xiàn)液滴飛散脫離靶材表面，但又處于加載區(qū)域的情況下，會對液滴下方區(qū)域造成遮蔽效果，但UDF在液滴所在的平面法線方向上掃描過程中會出現(xiàn)多個相界面。顯然，在多個相界面中只有一個是計算中的受熱表面，在程序設計中已考慮此種情況對計算精度的影響。,多個符

21、合判斷條件的網(wǎng)格,左圖存在三個單元網(wǎng)格滿足相界面判斷條件，但激光能量將只加載于擁有最大Z軸坐標的網(wǎng)格中。 因此在程序編寫過程中，通過首次的線程循環(huán)得到符合相體積分數(shù)條件的網(wǎng)格后，還將通過二次循環(huán)掃描來確定當前網(wǎng)格是否是沿靶材平面法線方向上擁有相同（X，Y）坐標且符合相體積分數(shù)條件的網(wǎng)格中Z軸坐標值最大的一個。,仿真過程的進一步簡化,,在課題研究的初期，采用LES方法計算大計算域的模型，約50min完成一次迭代步長（10-

22、6s量級），預計完成整個計算需用時400多天，這顯然無法實現(xiàn)。因此要謀求仿真過程的進一步簡化。 所采用的簡化方式如下： 1.UDF中涉及到網(wǎng)格的二次掃描，因此減少網(wǎng)格數(shù)量可以明顯提升計算效率； 2.使用其他湍流模型完成計算（k-ε）； 3.利用相似準則處理模型，增大迭代時間步長。,相似準則在模型中的應用,,當模型計算至靶面材料熔化后產(chǎn)生的多相流階段時，更有流動特征時

23、間尺度（不大于10-5s）與傳熱特征時間尺度（約10-2s）差異過大的問題，嚴重影響計算的效率。 在前人的研究中，證實了相似準則在激光壁面加載模型中的適用性。但之前的研究中均不涉及壁面受熱熔化后的過程，因此僅采用幾何相似準則來建立縮比計算模型。 本文利用量綱理論分析多項流流動的相似準則，再結合兩相間作用的定性分析建立起一套完整的多相流對流換熱的相似準則體系。并通過改變物性參數(shù)比例尺來調整非穩(wěn)態(tài)模型的計算特征時

24、間，縮小流動特征時間尺度與傳熱特征時間尺度之間的差異，從而能夠選用更大的迭代時間步長，達到提升計算效率的目的。,,,若要使得簡化后的計算模型保持與原型相似的流動傳熱現(xiàn)象，則相似模型中流體歐拉數(shù)Eu、雷諾數(shù)Re、馬赫數(shù)Ma、普朗特數(shù)Pr、弗勞德數(shù)Fr、韋伯數(shù)We、傅立葉數(shù)Fo、流場邊界條件及初始條件應與原始模型相同。相似準則的推導過程不再贅述。 使用如下參數(shù)的對比模型進行試算：,1:1,1:2,由左側兩圖可知，兩個模型的相界

25、面沿氣流流動方向的推移情況幾乎一致，模擬效果良好。而1：1模型實際計算時長10412s， 1：2模型實際計算時長5981s，相似處理后的模型縮短計算時間約42.6%，顯著提升計算效率。,相似準則適用性及誤差分析,,在論文中，對上述相似處理方法進行的適用性的驗證與誤差分析，此處不在贅言。這套多相流對流換熱的相似準則在兩相間存在一相高粘度流體，且重力對流動過程影響較低的情況下，選擇恰當?shù)目s放比例（原則上不大于4，本次課題選擇的比例為2）進行

26、相似變換后的模型計算出的結果與原始模型相符，誤差約為6.5%。 經(jīng)過模型試算，存在附面層影響下的激光壁面加載模型滿足以上限定條件，由于降低流場風速，故在非穩(wěn)態(tài)計算中使用的時間步長可增大2-4倍，計算效率顯著提升。其計算誤差將主要體現(xiàn)在人為降低50%對流換熱效率（抑制氣動冷卻效果），但考慮到氣流冷卻效果相對于激光能量的加熱效果來說是個小量，綜合來看，這套相似準則體系適用于本次課題的計算，在后續(xù)的計算中，將默認

27、采用比例系數(shù)為2的相似模型計算流場結果。,五.算例結果分析,由于靶面材料在熔化之前的計算模擬過程已經(jīng)由林銳在《高速近壁流場壁面加載特性CFD研究》一文中詳細分析，故此處不再贅述，僅給出以下重要仿真結果的分析說明：由于壁面溫度的升高，對流換熱換熱量增大，在熱損失增大而加載能量總量不變的情況下，壁面材料溫度上升的趨勢逐漸變緩，附面層氣流剪應力和Nu數(shù)均與壁面溫度變化成反比例關系。在亞聲速范圍內，激光加熱溫升遠大于氣流恢復溫度，可以不

28、考慮氣動加熱。切向附面層氣流增加了激光加熱材料所需的能量或功率，延緩靶面材料達到熔點的時間約12.06%。對流損失的功率密度不可忽略，影響激光損傷閾值。,無附面層流場影響下的激光壁面加載機理,,,下圖為靶材上表面首層位于輻照中心點的網(wǎng)格溫升曲線。計算結果顯示，加載進行至0.248s時網(wǎng)格溫度到達靶面金屬材料的凝固點1630K，繼續(xù)加載能量直至0.304s后模型認定網(wǎng)格中所有固相金屬材料均已加熱至熔化狀態(tài)，溫度繼續(xù)上升。在整個熔化過程中，

29、此溫升平臺期占據(jù)時長約18.4%，對靶材熔化初期的影響不能忽略。 隨著激光輻照能量的增大，此平臺期持續(xù)時間將縮短，但仍能夠在材料熔點附近觀察到溫度波動。,,,,上圖展示了熔化過程進行伊始，靶材金屬材料中上表面的溫度分布?？梢缘弥?，在激光輻照區(qū)域具有相當大的溫度梯度，而在此區(qū)域外的壁面溫度幾乎未受影響。盡管金屬有著很大的導熱系數(shù)，但是在激光加載這類型的瞬時強熱作用下，其加熱作用僅限于光斑附近的表面。,輻照區(qū)域,,,,由于激光

30、能量加載的高斯分布，輻照光斑邊緣處所加載的能量尚不及輻照中心處的1/10，因此越遠離輻照中心處的靶面區(qū)域溫度達到熔化溫度所需的能量加載時間就越長。如上圖所示，熔池半徑在熔化開始后不到0.03s的時間中便已增大至2mm，隨后熔化邊界的推移速度便出現(xiàn)了明顯的下降，至0.56s時熔池半徑才增大至4mm。,,,,,如右圖所示，計算于0.58s截止時，熔池中最高溫度2538K，緊靠靶材上表面輻照區(qū)域的空氣受熱明顯，氣液兩相之間換熱旺盛。與此相對的

31、，由熔化邊界所隔開的固液兩相間存在較大的溫度梯度，金屬物質間導熱效率不及氣流的換熱效率，其熔池正下方首層網(wǎng)格中溫度尚不足750K，驗證了熔融金屬較低的導熱性。,當計算模型中沒有引入附面層氣流影響時，熔池的發(fā)展及溫度場分布都呈現(xiàn)出以輻照中心對稱的特點，在計算時長中，熔池邊界尚未擴展到激光光斑邊緣處，熔池深度也沒有進一步的發(fā)展。,存在附面層流場影響下的激光壁面加載機理,,,無風狀態(tài)下的壁面加載特性在整體上呈現(xiàn)中心對稱，當引入附面層氣流的影響

32、之后，這種對稱性將被打破，能量加載特性將在氣流流向方向上呈現(xiàn)不對稱的特點。 模型中引入附面層流場后，分析其速度及溫度梯度分布：,附面層速度梯度,附面層溫度梯度,,,,從速度梯度示意圖中可以看到：氣流進口處的氣流數(shù)據(jù)由于采用profile文件的方式導入，已屬于穩(wěn)定流動，因此看不到采用氣流噴嘴時噴嘴處速度梯度的劇烈振蕩。氣流在由進口流至輻照區(qū)域低壓區(qū)之前約3mm處時，速度梯度出現(xiàn)激增，近壁區(qū)域氣流速度降低，與主流之間的速度差增

33、大。而在輻照區(qū)域及其后的流動過程中速度梯度保持一個相對小值，形成穩(wěn)定的附面邊界層流場。隨著邊界層流場的形成與擴展，降低了氣流與壁面之間由于對流換熱形成的熱損失，對于激光能量加載過程是有利的。 由溫度梯度曲線可以知道：氣流在達到輻照區(qū)域之前幾乎沒有溫度的變化。流經(jīng)輻照區(qū)域時，溫度梯度急劇減小。隨著靠近輻照中心處，溫度梯度的絕對值達到最大值。當氣流通過輻照中心后，壁面溫度與氣流溫度間的差距逐漸減小，溫度梯度開始增大。

34、 在靶面材料熔化初期，附面層氣流對激光加載特性的影響在輻照區(qū)前端主要表現(xiàn)為氣動冷卻效果，而在輻照區(qū)后端，在邊界層厚度增加及氣流溫度高于壁溫（熔化開始初期）的雙重作用下，附面層氣流在距輻照中心3mm之后的輻照區(qū)域甚至出現(xiàn)了對壁面的加熱效果。,熔池半徑發(fā)展的計算與分析,,,,三張圖片分別表示當熔池深度為0.5mm/1.0mm/1.5mm時，熔池半徑沿附面層氣流流向上的變化趨勢,,,,通過對熔池實際發(fā)展過程的數(shù)據(jù)分析，發(fā)現(xiàn)在熔池形成的初期

35、，能夠直觀的看到，熔池邊界向輻照區(qū)域前緣推移的速度（黑線）低于向后緣推移的速度（紅線）。 造成這一現(xiàn)象的主要原因是：輻照區(qū)域前緣及中心處的氣流速度邊界層厚度在氣體升溫及壁面不光滑兩種機理的作用下迅速增大。隨著邊界層厚度的增大，近壁面的空氣速度下降，其通過表面對流換熱效應帶走的熱量也開始下降，且越接近后緣處氣流速度越低，帶走的熱量也就越少。且氣流在流經(jīng)輻照區(qū)域后緣時，空氣溫度已高于壁面溫度，熱量開始由氣流傳至壁面。輻照區(qū)域前

36、后兩端的熱環(huán)境存在差異，造成了熔池發(fā)展的不對稱。,,,氣流方向,熔池不對稱發(fā)展示意圖（0.21325s）,,,,當熔化開始深入靶材金屬材料內部后，熔池的發(fā)展趨勢將和初期發(fā)展呈現(xiàn)明顯的區(qū)別，此時熔化邊界向輻照區(qū)域前緣推移的推移速度（黑線）反而高于向后緣推移的速度（紅線），這一現(xiàn)象在熔化深入至第三層網(wǎng)格（1.5mm）后尤為明顯。 隨著熔池深度的增加，熔融金屬液將越來越難以被吹離熔池，最終部分熔融金屬液將滯留于熔池后半端，避免了

37、激光能量直接輻照在熔池底部的金屬壁面上，變相起到了隔離作用。這一現(xiàn)象帶來的保護作用遠遠高于氣流對壁面的加熱效果，因此熔化邊界向輻照區(qū)域后緣推移的速度大大降低，而輻照前緣部分直接暴露于激光輻照之下，熔化速度快于后緣。,熔池不對稱發(fā)展示意圖（0.36825s）,,,氣流方向,熔池深度發(fā)展分析,,當附面層氣流不存在時，經(jīng)0.56s的輻照時間，熔化仍停留于激光加載區(qū)域，未出現(xiàn)向靶材深處發(fā)展的現(xiàn)象。 當模型中加入附面層氣流的影響后，

38、從數(shù)據(jù)分析中可以得知，當靶材表面熔化深入至0.5mm時用時0.248s，熔化深入至1mm用時0.20825s，深入至1.5mm時用時0.065s，熔化邊界向靶材內部推移的速度呈明顯加速趨勢。 出現(xiàn)這一現(xiàn)象的主要原因，是因為熔池形成后熔融金屬脫離輻照區(qū)域后形成的凹坑，使得激光加載能量能夠直接深入靶材內部，對于能量加載是有利的。若引入化學反應模型后，隨著熔融金屬不斷脫離熔池或與空氣發(fā)生氧化反應，熔池深度將不斷增大，靶材最終將

39、被熔穿。,,六.結論,（1）課題中使用基于商業(yè)CFD軟件建立的用戶自定義函數(shù)程序在模型計算全程中表現(xiàn)良好，能夠實現(xiàn)其預定功能。（2）激光加載具有極高的過熱度，為高度不平衡過程，熔化過程進行迅速。（3）在能量加載過程中，輻照區(qū)域處壁面升溫明顯，金屬材料有著大的熱導率，但激光加載的熱作用可以認為僅限于輻照光斑所在的表面。（4）附面層氣流對激光壁面加載特性的影響不容忽略。但在熔化過程的不同階段，流場對激光壁面加載特性的影響方式也存在差異

40、（5）當存在附面層氣流作用時，能量加載初期形成的淺表熔池中的熔融金屬被氣流帶離輻照區(qū)域，熔池向靶材內部發(fā)展的速度大大高于無氣流作用下的計算模型。,對今后工作的建議,,（1）采用LES方法實現(xiàn)對三維環(huán)境下開啟VOF模型的物質熔化過程的計算研究，獲取更高的計算精度。（2）計算不同速度附面層氣流作用下對激光能量加載特性影響的研究數(shù)據(jù)，并在保證計算穩(wěn)定性的前提下延長計算時間。（3）本文計算分析了激光輻照下靶材表面金屬熔化的初期過程，當熔化