[32] Schou J. Physical aspects of the pulsed laser deposition technique: The stoichiometric transfer of material from target to film. Applied Surface Science. 2009;255:5191-5198 [33] Canulescu S, Papadopoulou E, Anglos D, Lippert T, Montenegro MJ, Georgiou S, Döbeli M, Wokaun A. Nanosecond and femtosecond ablation of La0.6Ca0.4CoO3: A comparison between plume dynamics and composition of the films. Applied Physics A. 2011;105:

[34] Bulgakova NM, Bulgakov AV, Bobrenok OF. Double layer effects in laser-ablation

[35] Cleyssens F, Cheesman A, Henley SJ, Ashfold MNR. Studies of the plume accompanying pulsed ultraviolet laser ablation of zinc oxide. Journal of Applied Physics. 2002;92:

[36] Hansen TN, Schou J, Lunney JG. Angle-resolved energy distributions of laser ablated

[37] Amoruso S, Armenante M, Berardi V, Bruzzese R, Velotta R, Wang X. High fluence visibleand ultraviolet laser ablation of metallic targets. Applied Surface Science. 1998;

silver ions in vacuum. Applied Physics Letters. 1998;72:1829-1831

ultraviolet laser ablation of metals. Physical Review B. 2003;67:104101

ultrashort-pulse laser processing of surfaces. Micromachines. 2014;5:1344-1372

Review E. 2004;69:026413

Applied Physics. 2015;117:083108

gold. Physical Review B. 2007;76:035435

plasma plumes. Physical Review E. 2000;62:5624-5634

2016;120:185901

167-176

6886-6894

127-129:1017-1022

[24] Bulgakova NM, Panchenko AN, Zhukov VP, Kudryashov SI, Pereira A, Marine W, Mocek T, Bulgakov AV. Impacts of ambient and ablation plasmas on short- and ultrashort-pulse laser processing of surfaces. Micromachines. 2014;5:1344-1372

[11] Leitz KH, Redlingshöfer B, Reg Y, Otto A, Schmidt M. Metal ablation with short and

[12] Bulgakova NM, Stoian R, Rosenfeld A, Hertel IV, Marine W, Campbell EEB. A general continuum approach to describe fast electronic transport in pulsed laser irradiated materials: The problem of Coulomb explosion. Applied Physics A: Materials Science &

[13] Harilal SS, Freeman JR, Diwakar PK, Hassanein A. Femtosecond Laser Ablation: Fundamentals and Applications. In: Musazzi S, Perini U, editors. Laser-Induced Breakdown

[14] Geohegan DB, Puretzky AA. Dynamics of laser ablation plume penetration through low

[15] Geohegan DB, Puretzky AA. Laser ablation plume thermalization dynamics in background gases: Combined imaging, optical absorption and emission spectroscopy, and

[16] Wood RF, Chen KR, Leboeuf JN, Puretzky AA, Geohegan DB. Dynamics of plume propagation and splitting during pulsed-laser ablation. Physical Review Letters.

[17] Wood RF, Leboeuf JN, Geohegan DB, Puretzky AA, Chen KR. Dynamics of plume propagation and splitting during pulsed-laser ablation of Si in He and Ar. Physical

[18] Harilal SS, Issac RC, Bindhu CV, Nampoori VPN, Vallabhan CPG. Temporal and spatial evolution of C2 in laser induced plasma from graphite target. Journal of Applied Phys-

[19] Harilal SS, Issac RC, Bindhu CV, Nampoori VPN, Vallabhan CPG. Emission characteristics and dynamics of C2 from laser produced graphite plasma. Journal of Applied

[20] Harilal SS, Bindhu CV, Tillack MS, Najmabadi F, Gaeris AC. Plume splitting and sharpening in laser-produced aluminium plasma. Journal of Physics D: Applied Physics.

[21] Harilal SS, Bindhu CV, Tillack MS, Najmabadi F, Gaeris AC. Internal structure and expansion dynamics of laser ablation plumes into ambient gases. Journal of Applied

[22] Diwakar PK, Harilal SS, Hassanein A, Phillips MC. Expansion dynamics of ultrafast laser produced plasmas in the presence of ambient argon. Journal of Applied Physics.

[23] Bulgakov AV, Bulgakova NM. Gas-dynamic effects of the interaction between a pulsed laser ablation plume and the ambient gas: Analogy with an underexpanded jet. Journal

Spectroscopy Theory and Applications. New York: Springer; 2014. p. 143-166

pressure background gases. Applied Physics Letters. 1995;67:197-199

ion probe measurements. Applied Surface Science. 1996;96–98:131-138

ultrashort laser pulses. Physics Procedia. 2011;12:230-238

Processing. 2005;81:345-356

24 Laser Ablation - From Fundamentals to Applications

1997;79:1571-1574

Review B. 1998;58:1533-1543

ics. 1996;80:3561-3565

2002;35:2935-2938

2014;116:133301

Physics. 1997;81:3637-3643

Physics. 2003;93:2380-2388

of Physics D: Applied Physics. 1998;31:693-703

[38] Amoruso S, Bruzzese R, Spinelli N, Velotta R. Characterization of laser-ablation plasmas. Journal of Physics B: Atomic, Molecular and Optical Physics. 1999;32:R131-R172

[52] Gurlui S, Sanduloviciu M, Strat M, Strat G, Mihesan C, Ziskind M, Focsa C. Dynamic space charge structures in high fluence laser ablation plumes. Journal of Optoelectronics

Experimental and Theoretical Studies on the Dynamics of Transient Plasmas Generated by Laser Ablation…

http://dx.doi.org/10.5772/intechopen.70759

27

[53] Gurlui S, Sanduloviciu M, Mihesan C, Ziskind M, Focsa C. Periodic phenomena in laserablation plasma plumes: A self-organization scenario. AIP Conference Proceedings.

[54] Focsa C, Ziskind M, Ursu C, Gurlui S, Pagnon D, Pellerin S, Pellerin N, Dudeck M. Laser-BNSiO2 ceramics interaction: Simulation of the energy deposition on dielectric wall surfaces in Hall thrusters. Journal of Optoelectronics and Advanced Materials. 2008;10:2380-2385

[55] Gurlui S, Agop M, Nica P, Ziskind M, Focsa C. Experimental and theoretical investigations of a laser-produced aluminum plasma. Physical Review E. 2008;78:026405

[56] Ursu C, Gurlui S, Focsa C, Popa G. Space- and time-resolved optical diagnosis for the study of laser ablation plasma dynamics. Nuclear Instruments and Methods in Physics

[57] Focsa C, Nemec P, Ziskind M, Ursu C, Gurlui S, Nazabal V. Laser ablation of AsxSe100x chalcogenide glasses: Plume investigations. Applied Surface Science. 2009;255:5307-5311

[58] Nica P, Vizureanu P, Agop M, Gurlui S, Focsa C, Forna N, Ioannou P, Borsos Z. Experimental and theoretical aspects of Aluminum expanding laser plasma. Japanese Journal

[59] Agop M, Nica P, Gurlui S, Focsa C. Fractal hydrodynamic model of high fluence laser ablation plasma expansion. AIP Conference Proceedings. 2010;1278:612-622

[60] Nica P, Agop M, Gurlui S, Focsa C. Oscillatory Langmuir probe ion current in laser

[61] Cimpoesu RH, Pompilian GO, Baciu C, Cimpoesu N, Nejneru C, Agop M, Gurlui S, Focsa C. Pulsed laser deposition of poly(l-Lactide) acid on nitinol substrate. Optoelec-

[62] Ursu C, Pompilian OG, Gurlui S, Nica P, Agop M, Dudeck M, Focsa C. Al2O3 ceramics under high-fluence irradiation: Plasma plume dynamics through space-and timeresolved optical emission spectroscopy. Applied Physics A: Materials Science &

[63] Agop M, Nica PE, Gurlui S, Focsa C, Paun VP, Colotin M. Implications of an extended fractal hydrodynamic model. European Physical Journal D: Atomic, Molecular, Optical

[64] Focsa C. Laser ablation transient plasma structures expansion in vacuum. IEEE Trans-

[65] Nica P, Agop M, Gurlui S, Bejinariu C, Focsa C. Characterization of aluminum laser produced plasma by target current measurements. Japanese Journal of Applied Physics.

tronics and Advanced Materials, Rapid Communications. 2010;4:2148-2153

produced plasma expansion. Europhysics Letters. 2010;89:65001

actions on Plasma Sciences. 2011;39:Gurlui S, 2820-Gurl2821

and Advanced Materials. 2006;8:148-151

Research Section B. 2009;267:446-450

of Applied Physics. 2009;48:066001

Processing. 2010;101:153-159

2012;51:106102

and Plasma Physics. 2010;56:405-419

2006;812:279-282

[52] Gurlui S, Sanduloviciu M, Strat M, Strat G, Mihesan C, Ziskind M, Focsa C. Dynamic space charge structures in high fluence laser ablation plumes. Journal of Optoelectronics and Advanced Materials. 2006;8:148-151

[38] Amoruso S, Bruzzese R, Spinelli N, Velotta R. Characterization of laser-ablation plasmas. Journal of Physics B: Atomic, Molecular and Optical Physics. 1999;32:R131-R172

[39] Schou J, Amoruso S, Lunney JG. Plume dynamics. In: Phipps C, editor. Laser Ablation

[40] Doggett B, Lunney JG. Langmuir probe characterization of laser ablation plasmas.

[41] Esposito M, Lippert T, Schneider CW, Wokaun A, Donnelly T, Lunney JG, Tellez H, Vadillo JM, Laserna JJ. Pulsed laser ablation of silver: Ion dynamics in the plasma plume.

[42] Toftmann B, Doggett B, Budtz-Jorgensen C, Schou J, Lunney JG. Femtosecond ultraviolet laser ablation of silver and comparison with nanosecond ablation. Journal of Applied

[43] Chen J, Lunney JG, Lippert T, Ojeda-G-P A, Stender D, Schneider CW, Wokaun A. Langmuir probe measurements and mass spectrometry of plasma plumes generated by

[44] Amoruso S, Wang X, Altucci C, de Lisio C, Armenante M, Bruzzese R, Spinelli N, Velotta R. Double-peak distribution of electron and ion emission profile during femto-

[45] Amoruso S, Armenante M, Bruzzese R, Spinelli N, Velotta R, Wang X. Emission of prompt electrons during excimer laser ablation of aluminum targets. Applied Physics

[46] Sunil S, Kumar A, Singh RK, Subramanian KP. Measurements of electron temperature and density of multi-component plasma plume formed by laser-blow-off of LiF-C film.

[47] Kumar A, Singh RK, Thomas J, Sunil S. Parametric study of expanding plasma plume formed by laser-blow-off of thin film using triple Langmuir probe. Journal of Applied

[48] Singh SC, Fallon C, Hayden P, Mujawar M, Yeates P, Costello JT. Ion flux enhancements and oscillations in spatially confined laser produced aluminum plasmas. Physics of

[49] Hutchinson IH. Principles of Plasma Diagnostics. 2nd revised ed. Cambridge: Cam-

[50] Griem HR. Principles of Plasma Spectroscopy. Cambridge: Cambridge University Press;

[51] Kirkwood RK, Moody JD, Kline J, Dewald E, Glenzer S, Divol L, Michel P, Hinkel D, Berger R, Williams E, Milovich J, Yin L, Rose H, MacGowan B, Landen O, Rosen M, Lindl J. A review of laser–plasma interaction physics of indirect-drive fusion. Plasma

laser ablation of La0.4Ca0.6MnO3. Journal of Applied Physics. 2014;116:073303

second laser ablation of metals. Applied Surface Science. 2002;186:358-363

Journal of Physics D: Applied Physics. 2008;41:085211

Journal of Optoelectronics and Advanced Materials. 2010;12:677-680

and Its Applications. New York: Springer; 2007. p. 67-95

Journal of Applied Physics. 2009;105:033306

Physics. 2013;113:083304

26 Laser Ablation - From Fundamentals to Applications

Letters. 1999;75:7-9

Physics. 2009;106:043306

Plasmas. 2014;21:093113

1997

bridge University Press; 2005

Physics and Controlled Fusion. 2013;55:103001

[66] Balika L, Focsa C, Gurlui S, Pellerin S, Pellerin N, Pagnon D, Dudeck M. Laser induced breakdown spectroscopy in a running Hall effect Thruster for space propulsion. Spectrochimica Acta Part B. 2012;74–75:184-189

[79] Nica P, Gurlui S, Osiac M, Agop M, Ziskind M, Focsa C. Investigation of Femtosecond Laser-Produced Plasma from Various Metallic Targets using Langmuir Probe Charac-

Experimental and Theoretical Studies on the Dynamics of Transient Plasmas Generated by Laser Ablation…

http://dx.doi.org/10.5772/intechopen.70759

29

[80] Puretzky AA, Geohegan DB, Haufler RE, Hettich RL, Zheng XY, Compton RN. Laser ablation of graphite in different buffer gases. AIP Conference Proceedings. 1993;288:

[81] Harilal SS, Farid N, Freeman JR, Diwakar PK, LaHaye NL, Hassanein A. Background gas collisional effects on expanding fs and ns laser ablation plumes. Applied Physics A:

[82] Ojeda-G-P A, Schneider CW, Döbeli M, Lippert T, Wokaun A. Plasma plume dynamics, rebound, and recoating of the ablation target in pulsed laser deposition. Journal of

[83] Krása J, Jungwirth K, Gammino S, Krouský E, Láska L, Lorusso A, Nassisi V, Pfeifer M, Rohlena K, Torrisi L, Ullschmied J, Velyhan A. Partial currents of ion species in an

[84] Kelly R, Dreyfus RW. On the effect of Knudsen-layer formation on studies of vaporiza-

[85] Kabashin AV, Nikitin PI, Marine W, Sentis M. Experimental study of spontaneous electric field generated by a laser plasma. Applied Physics Letters. 1998;73:25-27 [86] Kabashin AV, Konov VI, Nikitin PI, Prokhorov AM, Konjević N, Vikor L. Laser plasma generation of currents along a conductive target. Journal of Applied Physics. 1990;68:

[87] Kabashin AV, Nikitin PI, Marine W, Sentis M. Electric fields of a laser plasma formed by optical breakdown of air near various targets. Quantum Electronics. 1998;28:24-28

[88] Wu J, Li X, Wei W, Jia S, Qiu A. Understanding plume splitting of laser ablated plasma: A view from ion distribution dynamics. Physics of Plasmas. 2013;20:113512

[89] Gornushkin IB, Panne U. Radiative models of laser-induced plasma and pump-probe diagnostics relevant to laser-induced breakdown spectroscopy. Spectrochimica Acta B.

[90] Elsied AM, Termini NC, Diwakar PK, Hassanein A. Characteristics of ions emission

[91] Weaver I, Martin GW, Graham WG, Morrow T, Lewis CLS. The Langmuir probe as a diagnostic of the electron component within low temperature laser ablated plasma

[92] Hansen TN, Schou J, Lunney JG. Langmuir probe study of plasma expansion in pulsed laser ablation. Applied Physics A: Materials Science & Processing. 1999;69:S601-S604

from Ultrashort laser produced plasma. Scientific Reports. 2016;6:38256

plume. The Review of Scientific Instruments. 1999;70:1801-1805

teristic. Phys. Plasmas. 2017; in press

Applied Physics. 2017;121:135306

Materials Science & Processing. 2014;117:319-326

expanding laser-created plasma. Vacuum. 2008;83:180-184

tion, sputtering, and desorption. Surface Science. 1988;198:263-276

365-374

3140-3146

2010;65:345-359

[79] Nica P, Gurlui S, Osiac M, Agop M, Ziskind M, Focsa C. Investigation of Femtosecond Laser-Produced Plasma from Various Metallic Targets using Langmuir Probe Characteristic. Phys. Plasmas. 2017; in press

[66] Balika L, Focsa C, Gurlui S, Pellerin S, Pellerin N, Pagnon D, Dudeck M. Laser induced breakdown spectroscopy in a running Hall effect Thruster for space propulsion.

[67] Pompilian OG, Gurlui S, Nemec P, Nazabal V, Ziskind M, Focsa C. Plasma diagnostics in pulsed laser deposition of GaLaS chalcogenides. Applied Surface Science. 2013;278:

[68] Balika L, Focsa C, Gurlui S, Pellerin S, Pellerin N, Pagnon D, Dudeck M. Laser ablation in a running Hall effect Thruster for space propulsion. Applied Physics A: Materials

[69] Dascalu G, Pompilian G, Chazallon B, Nica V, Caltun O, Gurlui S, Focsa C. Rare earth doped cobalt ferrite thin films deposited by PLD. Applied Physics A: Materials Science

[70] Dascalu G, Pompilian G, Chazallon B, Caltun O, Gurlui S, Focsa C. Femtosecond pulsed laser deposition of cobalt ferrite thin films. Applied Surface Science. 2013;278:38-42 [71] Irimiciuc SA, Mihaila I, Agop M. Experimental and theoretical aspects of a laser pro-

[72] Irimiciuc S, Agop M, Nica P, Gurlui S, Mihaileanu D, Toma S, Focsa C. Dispersive effects in laser ablation plasmas. Japanese Journal of Applied Physics. 2014;53:116202

[73] Bulai G, Gurlui S, Caltun OF, Focsa C. Pure and rare earth doped cobalt ferrite laser ablation: Space and time resolved optical emission spectroscopy. Digest Journal of

[74] Bulai G, Dumitru I, Pinteala M, Focsa C, Gurlui S. Magnetic nanoparticles generated by laser ablation in liquid. Digest Journal of Nanomaterials and Biostructures. 2016;

[75] Irimiciuc S, Boidin R, Bulai G, Gurlui S, Nemec P, Nazabal V, Focsa C. Laser ablation of (GeSe2)100x(Sb2Se3)x chalcogenide glasses: Influence of the target composition on the

[76] Irimiciuc SA, Gurlui S, Nica P, Focsa C, Agop MA. Compact non-differential approach for modelling laser ablation plasma dynamics. Journal of Applied Physics. 2017;121:

[77] Irimiciuc SA, Gurlui S, Bulai G, Nica P, Agop M, Focsa C. Langmuir probe investigation of transient plasmas generated by femtosecond laser ablation of several metals: Influence of the target physical properties on the plume dynamics. Applied Surface Science.

[78] Focsa C, Gurlui S, Nica P, Agop M, Ziskind M. Plume splitting and oscillatory behavior in transient plasmas generated by high-fluence laser ablation in vacuum. Appl. Surf. Sci.

plasma plume dynamics. Applied Surface Science. 2017;418B:594-600

Spectrochimica Acta Part B. 2012;74–75:184-189

Science & Processing. 2013;112:123-127

duced plasma. Physics of Plasmas. 2014;21:093509

Nanomaterials and Biostructures. 2015;10:1043-1053

& Processing. 2013;110:915-922

28 Laser Ablation - From Fundamentals to Applications

352-356

11:283-293

083301

2017;417:108-118

2017;424P3:299-309

[93] Dogar AH, Ilyas B, Ullah S, Nadeem A, Qayyum A. Langmuir probe measurements of Nd-YAG laser-produced copper plasmas. IEEE Transactions on Plasma Sciences. 2011;39:897-900

**Chapter 2**

Provisional chapter

**Nanosecond Laser Ablation: Mathematical Models,**

DOI: 10.5772/intechopen.70773

Nanosecond Laser Ablation: Mathematical Models,

The basic mathematical models, computational algorithms, and results of mathematical modeling of various modes of laser action on metals are considered. It is shown that for mathematical description and analysis of the processes of laser heating, melting, and evaporation of condensed media, various theoretical approaches are used: continuum, kinetic, atomistic, etc. Each of them has its own field of applicability, its advantages, and disadvantages. Mathematical description of ns-laser ablation is usually carried out within the framework of continuum approach in the form of hydrodynamic models that take into account reaction of irradiated material to varying density, pressure, and energy both in the target and in the vapor-gas medium. Within the framework of continuum approach, a multiphase, multifront hydrodynamic model and computational algorithm were constructed that were designed for modeling ns-PLA of metal targets embedded in gaseous media. It is shown that proposed model and computational algorithm allow to carry out the simulation of interrelated mechanisms of heterogeneous and homogeneous evaporation of metals manifested as a series of explosive boiling. Modeling has shown that explosive boiling in metals occurs due to the presence of a near-surface temperature maximum. It has been established that in ns-PLA, exposure regimes can be realized in which a phase explosion is the main mechanism of material removal. The verification of reliability of obtained results was carried out by comparing experimental

Keywords: nanosecond pulse, laser action, hydrodynamic model, mathematical modeling,

The pulsed laser ablation (PLA) of condensed media has been intensively studied over the past few decades [1, 2]. The increased interest in PLA is determined by the increasing possibilities

> © The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and eproduction in any medium, provided the original work is properly cited.

© 2017 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

explosive boiling, phase explosion, subsurface temperature maximum

**Computational Algorithms, Modeling**

Computational Algorithms, Modeling

Additional information is available at the end of the chapter

Additional information is available at the end of the chapter

data and calculations with atomistic models.

http://dx.doi.org/10.5772/intechopen.70773

Vladimir I. Mazhukin

Vladimir I. Mazhukin

Abstract

1. Introduction

## **Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling** Provisional chapter Nanosecond Laser Ablation: Mathematical Models,

DOI: 10.5772/intechopen.70773

Vladimir I. Mazhukin

[93] Dogar AH, Ilyas B, Ullah S, Nadeem A, Qayyum A. Langmuir probe measurements of Nd-YAG laser-produced copper plasmas. IEEE Transactions on Plasma Sciences.

[94] Salle B, Chaleard C, Detalle V, Lacour JL, Mauchien P, Nouvellon C, Semerok A. Laser ablation efficiency of metal samples with UV laser nanosecond pulses. Applied Surface

[95] Schou J, Toftmann B, Amoruso S. Pulsed laser deposition: From basic processes to film

[96] Olivier M, Němec P, Boudebs G, Boidin R, Focsa C, Nazabal V. Photosensitivity of pulsed laser deposited Ge-Sb-se thin films. Optical Materials Express. 2015;5:781-793 [97] Afifi MA, Labib HH, Fadel M. Electrical and thermal properties of chalcogenide glass system Se75Ge25-xSbx. Applied Physics A: Materials Science & Processing. 1992;55:167-169

[98] Merches I, Agop M. Differentiability and Fractality in Dynamics of Physical Systems.

[99] Krása J, Lorusso A, Doria D, Belloni F, Nassisi V, Rohlena K. Time-of-flight profile of multiply-charged ion currents produced by a pulse laser. Plasma Physics and Controlled

[100] Doria D, Lorusso A, Belloni F, Nassisi V. Characterization of a nonequilibrium XeCl laserplasma by a movable Faraday cup. The Review of Scientific Instruments. 2004;75:387 [101] Kelly R, Miotello A. Comments on explosive mechanisms of laser sputtering. Applied

[102] Anisimov SI, Luk'yanchuk BS. Selected problems of laser ablation theory. Physics-

[103] Murakami M, Kang YG, Nishihara K, Fujioka S, Nishimura H. Ion energy spectrum of expanding laser-plasma with limited mass. Physics of Plasmas. 2005;12:62706

2011;39:897-900

30 Laser Ablation - From Fundamentals to Applications

Science. 1999;138-139:302-305

Singapore: World Scientific; 2015

Surface Science. 1996;96–98:205-215

Fusion. 2005;47:1339-1349

Uspekhi. 2002;45:293-324

deposition. Proceedings of SPIE. 2005;5830:1-10

Additional information is available at the end of the chapter Vladimir I. Mazhukin

Computational Algorithms, Modeling

http://dx.doi.org/10.5772/intechopen.70773 Additional information is available at the end of the chapter

#### Abstract

The basic mathematical models, computational algorithms, and results of mathematical modeling of various modes of laser action on metals are considered. It is shown that for mathematical description and analysis of the processes of laser heating, melting, and evaporation of condensed media, various theoretical approaches are used: continuum, kinetic, atomistic, etc. Each of them has its own field of applicability, its advantages, and disadvantages. Mathematical description of ns-laser ablation is usually carried out within the framework of continuum approach in the form of hydrodynamic models that take into account reaction of irradiated material to varying density, pressure, and energy both in the target and in the vapor-gas medium. Within the framework of continuum approach, a multiphase, multifront hydrodynamic model and computational algorithm were constructed that were designed for modeling ns-PLA of metal targets embedded in gaseous media. It is shown that proposed model and computational algorithm allow to carry out the simulation of interrelated mechanisms of heterogeneous and homogeneous evaporation of metals manifested as a series of explosive boiling. Modeling has shown that explosive boiling in metals occurs due to the presence of a near-surface temperature maximum. It has been established that in ns-PLA, exposure regimes can be realized in which a phase explosion is the main mechanism of material removal. The verification of reliability of obtained results was carried out by comparing experimental data and calculations with atomistic models.

Keywords: nanosecond pulse, laser action, hydrodynamic model, mathematical modeling, explosive boiling, phase explosion, subsurface temperature maximum

## 1. Introduction

The pulsed laser ablation (PLA) of condensed media has been intensively studied over the past few decades [1, 2]. The increased interest in PLA is determined by the increasing possibilities

© 2017 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

© The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and eproduction in any medium, provided the original work is properly cited.

of its use in a variety of applications, beginning with the already traditional ones: microprocessing [3, 4], pulsed laser deposition (PLD) [5, 6], laser-induced breakdown spectroscopy (LIBS) [7, 8], and new rapidly developing areas: production of nanomaterials [9, 10], surface nanostructuring [11], chemical, and physical synthesis [12]. The total effect of nanotechnology and photonics led to the emergence of a new direction—laser synthesis of colloids [13], which draws general attention to its extensive applicability, primarily in biomedicine [14, 15].

internal structure of plasma plumes and their spatiotemporal evolution in the process of plasma expansion in the background gas at atmospheric pressure is of decisive importance for many engineering applications. For this reason, nanosecond laser ablation continues to be an area of

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

http://dx.doi.org/10.5772/intechopen.70773

33

Any simulation begins with a choice of mathematical model, construction and development of which is given a paramount importance in the computing experiments. For theoretical description and analysis of PLA process of condensed media, various theoretical approaches are used: continuum, kinetic, atomistic (molecular dynamics, etc.). Each of them has its own field of

Atomistic models allow us to conduct research at the atomic level and obtain fundamental knowledge about structure, thermodynamic, and mechanical properties of crystalline materials [32, 33], about physical mechanisms of various processes [34, 35], including the kinetics of heterogeneous and homogeneous phase transitions [36, 37]. The basic methods of atomistic modeling—molecular dynamics (MD) and Monte Carlo (MC), which use as a rule, semiempirical interaction potentials, operate with tens and hundreds of millions of atoms, and allow

However, even with the use of parallel computer platforms, computational costs are enormous, and the space-time scales inherent for PLA processes are beyond the limits of accessibility for atomistic modeling methods. Therefore, in spite of constant progress in the field of designing interatomic potentials and increasing the power of computing systems, the final overcoming of computational constraints is hardly achievable, and continuum models will

Continuous models based on the equations of continuous medium mechanics are realized, as a rule, in the form of hydrodynamic models [18, 24, 28, 30, 31, 38–42] and use the minimum of information and operate with average values of physical characteristics calculated on infinitesimal volume. The methods for solving them are more compact, they have higher accuracy and a relatively small amount of computation. The main shortcomings of continuum approach are manifested in the absence of the possibility of direct investigation of elementary processes in materials and limited possibilities of mathematical description of homogeneous mechanisms of phase transitions of the first kind and calculation of thermophysical, thermodynamic, optical, and other characteristics of matter in a wide range of parameters. These problems are much easier and more fully solved within the framework of atomistic modeling, the results of

In this paper, the application of continuum approach to modeling in preplasma regime of processes dynamics and of main mechanisms of ns-PLA of metal target (Al) in air is considered. The mechanisms of heterogeneous and homogeneous phase transitions interacting with each other are analyzed in detail. Primary attention is paid to model the dynamics of phase explosion of liquid phase of aluminum and the expansion of its fragments in the air, since explosive boiling is considered to be one of the most effective thermal mechanisms of ns laser ablation of materials. Various aspects of this problem have been studied in a number of

which can be used as input parameters in meso- and macrolevel models.

active research in which mathematical modeling plays an increasing role [28–31].

applicability, its advantages, and disadvantages.

calculations in the time range of nanosecond duration.

always remain relevant.

Numerous applications make PLA an attractive direction for fundamental investigations. Despite extensive studies of fundamental properties of laser ablation performed earlier, a number of important physical phenomena still remain insufficiently well studied and understood. Previous studies have established that the nature of interaction of laser radiation depends both on the modes of action: wavelength [16, 17], duration [18, 19], and laser pulse intensity [20, 21] and on thermophysical and optical properties of the target [22, 23], presence of the surrounding gas [24] and its pressure [25, 26]. The greatest differences in the physical mechanisms of laser ablation of metals are observed between short (ns) and ultrashort (ps, fs) pulsed modes [18, 19].

In ultrashort range (fs, pc) of influence, laser radiation freely reaches surface of the target. Absorption of laser radiation by a degenerate electron gas followed by a slowed-down energy exchange between electron and phonon components leads to strong deviation from locally thermodynamic equilibrium of the system as a whole. As a result, laser ablation and its accompanying processes develop after the end of the pulse.

Laser ablation in the nanosecond range is a more complex phenomenon involving many interrelated processes both during action and after the end of laser pulse. Such processes include laser target heating, heterogeneous, and homogeneous phase transformations, taking place in the evaporated matter, formation and expansion of plasma plume, heat transfer, laser and intrinsic plasma radiation transfer, generation and propagation of shock waves, and contact boundaries in the vaporized matter and the surrounding gas environment. In contrast to ultrashort regime, in nanosecond range, with a certain choice of parameters of the action, two experimentally observed and explored phenomena arise—volumetric boiling (phase explosion) of liquid phase of the target and formation of laser plasma in vaporized matter and surrounding gas.

Putting an irradiated target into an external gas environment, which is typical for most PLA applications, significantly complicates the situation. In these cases, the long-lived processes of laser-plasma plume and associated generation and propagation of interacting fronts of shock waves and contact boundaries in the vapor-gas medium play an important role in overall picture of laser ablation [27]. The presence of such a large number of interconnected physical processes creates additional difficulties in determination and investigation of the basic mechanisms of ablation. The information obtained by the methods of instrumental diagnostics is not sufficient due to the lack of data on fundamental phenomena associated with rapidly changing energy (thermal and laser radiation), hydrodynamic fields, and the kinetics of heterogeneous and homogeneous phase transformations in the solid and liquid phases, ionization of vaporized matter, and gas environment. At the same time, the understanding of fundamental physics of internal structure of plasma plumes and their spatiotemporal evolution in the process of plasma expansion in the background gas at atmospheric pressure is of decisive importance for many engineering applications. For this reason, nanosecond laser ablation continues to be an area of active research in which mathematical modeling plays an increasing role [28–31].

of its use in a variety of applications, beginning with the already traditional ones: microprocessing [3, 4], pulsed laser deposition (PLD) [5, 6], laser-induced breakdown spectroscopy (LIBS) [7, 8], and new rapidly developing areas: production of nanomaterials [9, 10], surface nanostructuring [11], chemical, and physical synthesis [12]. The total effect of nanotechnology and photonics led to the emergence of a new direction—laser synthesis of colloids [13], which

Numerous applications make PLA an attractive direction for fundamental investigations. Despite extensive studies of fundamental properties of laser ablation performed earlier, a number of important physical phenomena still remain insufficiently well studied and understood. Previous studies have established that the nature of interaction of laser radiation depends both on the modes of action: wavelength [16, 17], duration [18, 19], and laser pulse intensity [20, 21] and on thermophysical and optical properties of the target [22, 23], presence of the surrounding gas [24] and its pressure [25, 26]. The greatest differences in the physical mechanisms of laser ablation of metals are observed between short (ns) and ultrashort (ps, fs)

In ultrashort range (fs, pc) of influence, laser radiation freely reaches surface of the target. Absorption of laser radiation by a degenerate electron gas followed by a slowed-down energy exchange between electron and phonon components leads to strong deviation from locally thermodynamic equilibrium of the system as a whole. As a result, laser ablation and its accom-

Laser ablation in the nanosecond range is a more complex phenomenon involving many interrelated processes both during action and after the end of laser pulse. Such processes include laser target heating, heterogeneous, and homogeneous phase transformations, taking place in the evaporated matter, formation and expansion of plasma plume, heat transfer, laser and intrinsic plasma radiation transfer, generation and propagation of shock waves, and contact boundaries in the vaporized matter and the surrounding gas environment. In contrast to ultrashort regime, in nanosecond range, with a certain choice of parameters of the action, two experimentally observed and explored phenomena arise—volumetric boiling (phase explosion) of liquid phase of the target and formation of laser plasma in vaporized matter

Putting an irradiated target into an external gas environment, which is typical for most PLA applications, significantly complicates the situation. In these cases, the long-lived processes of laser-plasma plume and associated generation and propagation of interacting fronts of shock waves and contact boundaries in the vapor-gas medium play an important role in overall picture of laser ablation [27]. The presence of such a large number of interconnected physical processes creates additional difficulties in determination and investigation of the basic mechanisms of ablation. The information obtained by the methods of instrumental diagnostics is not sufficient due to the lack of data on fundamental phenomena associated with rapidly changing energy (thermal and laser radiation), hydrodynamic fields, and the kinetics of heterogeneous and homogeneous phase transformations in the solid and liquid phases, ionization of vaporized matter, and gas environment. At the same time, the understanding of fundamental physics of

draws general attention to its extensive applicability, primarily in biomedicine [14, 15].

pulsed modes [18, 19].

32 Laser Ablation - From Fundamentals to Applications

and surrounding gas.

panying processes develop after the end of the pulse.

Any simulation begins with a choice of mathematical model, construction and development of which is given a paramount importance in the computing experiments. For theoretical description and analysis of PLA process of condensed media, various theoretical approaches are used: continuum, kinetic, atomistic (molecular dynamics, etc.). Each of them has its own field of applicability, its advantages, and disadvantages.

Atomistic models allow us to conduct research at the atomic level and obtain fundamental knowledge about structure, thermodynamic, and mechanical properties of crystalline materials [32, 33], about physical mechanisms of various processes [34, 35], including the kinetics of heterogeneous and homogeneous phase transitions [36, 37]. The basic methods of atomistic modeling—molecular dynamics (MD) and Monte Carlo (MC), which use as a rule, semiempirical interaction potentials, operate with tens and hundreds of millions of atoms, and allow calculations in the time range of nanosecond duration.

However, even with the use of parallel computer platforms, computational costs are enormous, and the space-time scales inherent for PLA processes are beyond the limits of accessibility for atomistic modeling methods. Therefore, in spite of constant progress in the field of designing interatomic potentials and increasing the power of computing systems, the final overcoming of computational constraints is hardly achievable, and continuum models will always remain relevant.

Continuous models based on the equations of continuous medium mechanics are realized, as a rule, in the form of hydrodynamic models [18, 24, 28, 30, 31, 38–42] and use the minimum of information and operate with average values of physical characteristics calculated on infinitesimal volume. The methods for solving them are more compact, they have higher accuracy and a relatively small amount of computation. The main shortcomings of continuum approach are manifested in the absence of the possibility of direct investigation of elementary processes in materials and limited possibilities of mathematical description of homogeneous mechanisms of phase transitions of the first kind and calculation of thermophysical, thermodynamic, optical, and other characteristics of matter in a wide range of parameters. These problems are much easier and more fully solved within the framework of atomistic modeling, the results of which can be used as input parameters in meso- and macrolevel models.

In this paper, the application of continuum approach to modeling in preplasma regime of processes dynamics and of main mechanisms of ns-PLA of metal target (Al) in air is considered. The mechanisms of heterogeneous and homogeneous phase transitions interacting with each other are analyzed in detail. Primary attention is paid to model the dynamics of phase explosion of liquid phase of aluminum and the expansion of its fragments in the air, since explosive boiling is considered to be one of the most effective thermal mechanisms of ns laser ablation of materials. Various aspects of this problem have been studied in a number of theoretical and experimental studies [43–56], but there is still no consensus on the mechanism of the phase explosion in metals. In order to obtain detailed information on interaction of heterogeneous and homogeneous mechanisms and data on laser plume morphology, simulation of laser heating, melting, surface evaporation, and evolution of plume in the vapor-gas medium is performed within the framework of new hydrodynamic model with temperature dependences of material properties of the target and explicit tracking of interphase boundary fronts, contact boundary, and shock wave. The release of liquid phase fragments into the atmosphere as a result of phase explosion is modeled by the procedure for introducing a quasinucleus of a new phase (vapor) in the region of near-surface maximum of temperature reaching the value Tmax ≈ 0.9Tcr at time of reaching the maximum permissible overheating of metastable liquid phase.

## 2. Hydrodynamic model

Laser ablation of Al target placed in the air atmosphere is considered. In preplasma regimes of the action in ablation, the processes in three phase states (solid, liquid, and vapor) and two states (perturbed, unperturbed) of external gaseous medium are taken into account. Laser radiation propagates from right to left. The air for the selected action mode is completely transparent to the laser flow, which is partially reflected from the metal surface and partially absorbed in the near-surface layer of the target.

∂ρ ∂t þ

front of shock wave.

∂ð Þ ρ u ∂t þ

∂ð Þ ρε ∂t þ

∂G

2.1. Boundary conditions

undergo a steplike change.

∂ð Þ ρ u ∂x ¼ 0,

> ∂ ρ u2 ð Þ ∂x þ ∂P

∂ð Þ ρ uε

∂x þ κL ρ; T � �G ¼ 0, WT ¼ �λð Þ T ∂T

∂x ¼ 0,

∂x þ ∂WT

Figure 1. Scheme of spatial position of the phases and the direction of motion of interphase, contact boundaries, and the

t > 0, Γs < x < Γsℓ ð Þt < x < Γℓvð Þt < Γυgð Þt < Γgð Þt

Here, ρ, u, ε, T, and P are the density, gas dynamic velocity, internal energy, temperature, and pressure of the substance, respectively, κL and G are the absorption coefficient and laser radiation fluence, WT is heat flux density, and λ is coefficient of thermal conductivity. Indices s, ℓ, v, g denote the belonging of the quantities to solid, liquid, vapor, and air, respectively. In the condensed phase, the value εk has the meaning of enthalpy of liquid and solid phases Hk.

In the hydrodynamic models that easily combine with the kinetic ones, the heterogeneous mechanisms of phase transitions: melting/crystallization and evaporation/condensation are described naturally. A distinctive feature of heterogeneous phase transitions is the presence of sharp interfaces at which the main thermophysical and optical characteristics: enthalpy H, heat capacity Cp and thermal conductivity λ, density ρ, pressure p, and reflectivity R of the surface

In the case of rapid phase transformations that typical for pulsed nanosecond action, a complete set of equations for mass, momentum, and energy flows is used to describe heterogeneous phase transitions. However, considering that the absence of phase equilibrium at interphase boundary is a necessary condition for phase transition to occur, the expressions for

� �,

∂x þ ∂G ∂x

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

k ¼ s, ℓ, υ, g (1)

http://dx.doi.org/10.5772/intechopen.70773

35

∂x ¼ � P ∂u

∂x P ¼ P ρ; T � �, ε ¼ ε ρ; T � �

The mathematical description of ns laser ablation within the continuum approach is realized in the form of hydrodynamic models that allow one to take into account the reaction of condensed (target) and vapor-gas (vaporized matter, gas) media to varying density, pressure, and energy. The processes in each medium are described by a system of nonstationary equations of gas-hydrodynamics supplemented by equations of energy with thermal conductivity, equation of laser radiation transfer, and corresponding equations of state. The 1-D approximation is used for spatial variables. For laser action regimes under consideration, conditions of locally thermodynamic equilibrium (LTE) are assumed for all processes both in the irradiated target and in the vapor-gas medium. Accordingly, hydrodynamic model is formulated in onetemperature approximation.

Phase states and vapor-gas medium are separated among themselves by moving interphase boundaries solid-liquid Гsℓ(t), liquid-vapor Гℓv(t), and by the fronts of contact boundary Гvg(t) and shock-wave Гsh,g(t). The right-hand boundary Гg(t), running along unperturbed gas, is declared moving in order to improve economic efficiency of the computational algorithm. Schematically, their position and direction of motion are shown in Figure 1.

In order to obtain complete information on the kinetics of phase transformations, the morphology and dynamics of laser plume, all moving fronts and boundaries in the course of the solution are subjected to explicit tracking using appropriate relationships, which are simultaneously the boundary conditions for hydrodynamic and energy equations.

Complete system of equations is represented as:

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling http://dx.doi.org/10.5772/intechopen.70773 35

Figure 1. Scheme of spatial position of the phases and the direction of motion of interphase, contact boundaries, and the front of shock wave.

$$\begin{cases} \frac{\partial\rho}{\partial t} + \frac{\partial\{\rho u\}}{\partial x} = 0, \\ \frac{\partial\{\rho u\}}{\partial t} + \frac{\partial\{\rho u^2\}}{\partial x} + \frac{\partial\{p}}{\partial x} = 0, \\ \frac{\partial\{\rho u\}}{\partial t} + \frac{\partial\{\rho u\}}{\partial x} = -\left(P\frac{\partial u}{\partial x} + \frac{\partial W\_T}{\partial x} + \frac{\partial G}{\partial x}\right), \\ \frac{\partial G}{\partial x} + \kappa\_L\{\rho, T\}G = 0, \\ W\_T = -\lambda(T)\frac{\partial T}{\partial x} \\ P = P(\rho, T), \quad \varepsilon = \varepsilon\{\rho, T\} \end{cases} \quad k = s, \ell, \nu, \varrho \tag{1}$$

$$t > 0, \quad \Gamma\_s < \mathbf{x} < \Gamma\_{\mathcal{S}l}(t) < \mathbf{x} < \Gamma\_{\mathcal{U}l}(t) < \Gamma\_{\mathcal{V}g}(t) < \Gamma\_{\mathcal{S}l}(t)$$

Here, ρ, u, ε, T, and P are the density, gas dynamic velocity, internal energy, temperature, and pressure of the substance, respectively, κL and G are the absorption coefficient and laser radiation fluence, WT is heat flux density, and λ is coefficient of thermal conductivity. Indices s, ℓ, v, g denote the belonging of the quantities to solid, liquid, vapor, and air, respectively. In the condensed phase, the value εk has the meaning of enthalpy of liquid and solid phases Hk.

#### 2.1. Boundary conditions

theoretical and experimental studies [43–56], but there is still no consensus on the mechanism of the phase explosion in metals. In order to obtain detailed information on interaction of heterogeneous and homogeneous mechanisms and data on laser plume morphology, simulation of laser heating, melting, surface evaporation, and evolution of plume in the vapor-gas medium is performed within the framework of new hydrodynamic model with temperature dependences of material properties of the target and explicit tracking of interphase boundary fronts, contact boundary, and shock wave. The release of liquid phase fragments into the atmosphere as a result of phase explosion is modeled by the procedure for introducing a quasinucleus of a new phase (vapor) in the region of near-surface maximum of temperature reaching the value Tmax ≈ 0.9Tcr at time of reaching the maximum permissible overheating of

Laser ablation of Al target placed in the air atmosphere is considered. In preplasma regimes of the action in ablation, the processes in three phase states (solid, liquid, and vapor) and two states (perturbed, unperturbed) of external gaseous medium are taken into account. Laser radiation propagates from right to left. The air for the selected action mode is completely transparent to the laser flow, which is partially reflected from the metal surface and partially

The mathematical description of ns laser ablation within the continuum approach is realized in the form of hydrodynamic models that allow one to take into account the reaction of condensed (target) and vapor-gas (vaporized matter, gas) media to varying density, pressure, and energy. The processes in each medium are described by a system of nonstationary equations of gas-hydrodynamics supplemented by equations of energy with thermal conductivity, equation of laser radiation transfer, and corresponding equations of state. The 1-D approximation is used for spatial variables. For laser action regimes under consideration, conditions of locally thermodynamic equilibrium (LTE) are assumed for all processes both in the irradiated target and in the vapor-gas medium. Accordingly, hydrodynamic model is formulated in one-

Phase states and vapor-gas medium are separated among themselves by moving interphase boundaries solid-liquid Гsℓ(t), liquid-vapor Гℓv(t), and by the fronts of contact boundary Гvg(t) and shock-wave Гsh,g(t). The right-hand boundary Гg(t), running along unperturbed gas, is declared moving in order to improve economic efficiency of the computational algorithm.

In order to obtain complete information on the kinetics of phase transformations, the morphology and dynamics of laser plume, all moving fronts and boundaries in the course of the solution are subjected to explicit tracking using appropriate relationships, which are simulta-

Schematically, their position and direction of motion are shown in Figure 1.

neously the boundary conditions for hydrodynamic and energy equations.

metastable liquid phase.

2. Hydrodynamic model

34 Laser Ablation - From Fundamentals to Applications

temperature approximation.

absorbed in the near-surface layer of the target.

Complete system of equations is represented as:

In the hydrodynamic models that easily combine with the kinetic ones, the heterogeneous mechanisms of phase transitions: melting/crystallization and evaporation/condensation are described naturally. A distinctive feature of heterogeneous phase transitions is the presence of sharp interfaces at which the main thermophysical and optical characteristics: enthalpy H, heat capacity Cp and thermal conductivity λ, density ρ, pressure p, and reflectivity R of the surface undergo a steplike change.

In the case of rapid phase transformations that typical for pulsed nanosecond action, a complete set of equations for mass, momentum, and energy flows is used to describe heterogeneous phase transitions. However, considering that the absence of phase equilibrium at interphase boundary is a necessary condition for phase transition to occur, the expressions for conservation laws must be supplemented by appropriate kinetic relations characterizing the degree of nonequilibrium of phase transition.

#### 2.1.1. Left fixed boundary, x = Гs

The condition that flow of mass and heat be equal to zero is used as boundary conditions on the left fixed boundary:

$$\mathbf{x} = \Gamma\_s: \qquad \boldsymbol{\mu} = \mathbf{0}, \quad \boldsymbol{W}\_T = \mathbf{0}; \tag{2}$$

Heterogeneous evaporation is characterized by a sharp phase boundary. Thin (several mean path lengths) nonequilibrium Knudsen layer (KL) is adjacent to the boundary surface. The nonequilibrium of the KL is determined by the flow of matter through the phase boundary. The total flow consists of flow of particles (atoms, molecules) leaving the condensed phase and flow of particles (atoms of molecules) returning from the evaporated matter as a result of

leads to a strong nonequilibrium. The flow of injected particles satisfies the Maxwellian distri-

is determined by conditions of gas-dynamic expansion outside the KL. The behavior of particles inside a nonequilibrium Knudsen layer is described by Boltzmann equation, by solving of which one can determine the unknown parameters Tυ and ρυ that are boundary conditions for the continuum equations of gaseous medium. In the heterogeneous evaporation model, a special role is played by parameter M (M=uυ/usound is Mach number) on the outer side of KL,

In the phase equilibrium state, when the pressure of saturated vapor psat(Tsur) is equal to the external pressure pυ, the parameter M = 0. In the subsonic evaporation regime, when M < 1, the behavior of the interface depends on gas-dynamic perturbations in the vaporized matter flow. The maximum nonequilibrium is determined by the maximum value of mass flow, which is known to be achieved at M = 1, when flow of matter through the interface is maximal, and the recoil pressure at the heated surface is minimal. In this case, such evaporation regime can be realized when the behavior of condensed medium no longer depend on the external gasdynamic problem, which greatly simplifies the description of evaporation process. However, the implementation area for such regime remains uncertain, because due to the nonlinearity of the gas dynamics equations, the transition line between M < 1 and M ≥ 1 depends on desired solution and cannot be determined in advance [64]. Taking into account strong spatiotemporal diversity of processes in the nonequilibrium layer and continual media, KL is usually represented as a strong discontinuity in the gas-dynamic parameters. In such case, the kinetic processes in KL are not explicitly considered, and various phenomenological approaches [65–72] were used to determine boundary conditions on the outer side of the KL, which makes it possible without solving the kinetic problem, to determine the joining conditions under certain assumptions about the type of the nonequilibrium distribution function inside the discontinuity. One of the first works in which intensive evaporation was analyzed on the basis of phenomenological model is the work of Crout [65]. In this work, a nonequilibrium ellipsoidal Maxwellian particle distribution function was used written in an analytical form with the anisotropic in spatial directions longitudinal and perpendicular temperatures. Later papers used Mott-Smith approach [66], which was applied to the structure of the shock waves. With the help of this approach [67, 68] expressions for gas-dynamic parameters at M = 1 were obtained. Later in [69], this approach was extended to the entire range of evaporation 0 ≤ M ≤ 1. In [70, 71], the approximation of the distribution function over the entire evaporation range 0 ≤ M ≤ 1 was carried out by other more complicated expressions that satisfied the additional requirements that were formulated taking into account basic laws of gas dynamics. Following this laws, the used distribution functions must when M ! 1 to provide extreme values for all three complete flows of mass jm,

(+) with parameters of the surface of the condensed medium. The return flow

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

(+) and f

(�) respectively, and that

37

http://dx.doi.org/10.5772/intechopen.70773

collisions. These flows have different distribution functions, f

which determines the degree of nonequilibrium of the phase transition.

bution function f

#### 2.1.2. Model of heterogeneous (surface) melting: crystallization, õ = Гsℓ(t)

For fast phase transitions, heterogeneous melting model consists of a system of equations expressing three conservation laws: mass, momentum, and energy supplemented by kinetic condition for the velocity of melting front υsℓ [57] obtained from molecular-kinetic theory [58] and which is the main characteristic of the process melting-crystallization. In a stationary (laboratory) coordinate system, surface melting-crystallization model written on moving interface of melting x = Гsℓ(t) can be represented as:

$$\mathbf{x} = \Gamma\_{s\ell}(t): \qquad \rho\_s(\mathbf{u}\_s - \nu\_{s\ell}) = \rho\_\ell(\mathbf{u}\_\ell - \nu\_{s\ell})$$

$$p\_s + \rho\_s(\mathbf{u}\_s - \nu\_{s\ell})^2 = p\_\ell + \rho\_\ell(\mathbf{u}\_\ell - \nu\_{s\ell})^2 \tag{3}$$

$$\left(\lambda(T)\frac{\partial T}{\partial \mathbf{x}}\right)\_s - \left(\lambda(T)\frac{\partial T}{\partial \mathbf{x}}\right)\_\ell = \rho\_s L\_m^{ne} \nu\_{s\ell}$$

$$\nu\_{s\ell}(\Delta T\_{s\ell}) = a(3k\_BT\_{s\ell}/m)^{1/2} \left(\exp\left(\beta \frac{L\_m}{k\_BT\_m(\mathbf{p}\_s)} \frac{\Delta T\_{s\ell}}{T\_{s\ell}}\right) - 1\right) \tag{4}$$

$$L\_m^{ne} = L\_m\left(T\_m\left(p\_s\right)\right) + \Delta\mathsf{C}\_{ps}\Delta T\_{s\ell} + \frac{\rho\_s + \rho\_\ell}{\rho\_s - \rho\_\ell}\frac{\left(u\_s - u\_\ell\right)^2}{2}, \quad L\_m\left(T\_m\left(p\_s\right)\right) = L\_{m,0} + \delta\left(T\_m\left(p\_s\right) - T\_{m,0}\right),$$

$$\Delta\mathsf{C}\_{ps} = \mathsf{C}\_{ps} - \mathsf{C}\_{p\ell}, \quad \Delta T\_{s\ell} = T\_{s\ell} - T\_m\left(p\_s\right), \quad T\_m\left(p\_s\right) = T\_{m,0} + \partial p\_{ss}$$

where Lne m is nonequilibrium heat of melting, Lm,0, Tm,0 are equilibrium heat of melting and melting point, respectively, ps is pressure on a solid surface. α, β, and δ are parameters determined from molecular modeling [59, 60]. For Al, α = 0.21, β = 5.28, δ = 6.37 J mol� 1 K�1 , θ = 6.44 � 10�3 K/atm.

#### 2.1.3. Model of heterogeneous evaporation x = Гlv

Investigation of heterogeneous evaporation process began already in the century before last from experimental work: Hertz [61] and theoretical work: Knudsen [62] and continues at the present time, which is determined by its practical importance and not fully clarified features of nonequilibrium behavior of the substance during its evaporation. The mechanism of heterogeneous laser evaporation is realized in subcritical region of surface temperature Tsur and saturated vapor pressure psat = psat(Tsur), for which the following inequalities hold: Tsur < Tcr and psat < pcr. In the case Tsur > Tcr or psat > pcr, then a supercritical laser evaporation regime is realized, at which the state of the substance varies continuously [63].

Heterogeneous evaporation is characterized by a sharp phase boundary. Thin (several mean path lengths) nonequilibrium Knudsen layer (KL) is adjacent to the boundary surface. The nonequilibrium of the KL is determined by the flow of matter through the phase boundary. The total flow consists of flow of particles (atoms, molecules) leaving the condensed phase and flow of particles (atoms of molecules) returning from the evaporated matter as a result of collisions. These flows have different distribution functions, f (+) and f (�) respectively, and that leads to a strong nonequilibrium. The flow of injected particles satisfies the Maxwellian distribution function f (+) with parameters of the surface of the condensed medium. The return flow is determined by conditions of gas-dynamic expansion outside the KL. The behavior of particles inside a nonequilibrium Knudsen layer is described by Boltzmann equation, by solving of which one can determine the unknown parameters Tυ and ρυ that are boundary conditions for the continuum equations of gaseous medium. In the heterogeneous evaporation model, a special role is played by parameter M (M=uυ/usound is Mach number) on the outer side of KL, which determines the degree of nonequilibrium of the phase transition.

conservation laws must be supplemented by appropriate kinetic relations characterizing the

The condition that flow of mass and heat be equal to zero is used as boundary conditions on

For fast phase transitions, heterogeneous melting model consists of a system of equations expressing three conservation laws: mass, momentum, and energy supplemented by kinetic condition for the velocity of melting front υsℓ [57] obtained from molecular-kinetic theory [58] and which is the main characteristic of the process melting-crystallization. In a stationary (laboratory) coordinate system, surface melting-crystallization model written on moving inter-

x ¼ Γsℓ ð Þt : ρsð Þ¼ us � υsℓ ρℓ ð Þ uℓ � υsℓ

ps þ ρsð Þ us � υsℓ

1=2

s � λð Þ T ∂T

2 , Lm Tm ps

melting point, respectively, ps is pressure on a solid surface. α, β, and δ are parameters determined from molecular modeling [59, 60]. For Al, α = 0.21, β = 5.28, δ = 6.37 J mol� 1 K�1

Investigation of heterogeneous evaporation process began already in the century before last from experimental work: Hertz [61] and theoretical work: Knudsen [62] and continues at the present time, which is determined by its practical importance and not fully clarified features of nonequilibrium behavior of the substance during its evaporation. The mechanism of heterogeneous laser evaporation is realized in subcritical region of surface temperature Tsur and saturated vapor pressure psat = psat(Tsur), for which the following inequalities hold: Tsur < Tcr and psat < pcr. In the case Tsur > Tcr or psat > pcr, then a supercritical laser evaporation regime is

m is nonequilibrium heat of melting, Lm,0, Tm,0 are equilibrium heat of melting and

∂x � �

exp β Lm

λð Þ T ∂T ∂x � �

ð Þ us�uℓ 2

υsℓ ð Þ¼ ΔTsℓ αð Þ 3kBTsℓ=m

ΔCps =Cps � Cpℓ , ΔTsℓ = Tsℓ� Tm(ps) , Tm(ps) = Tm , 0 + θps,

ρs�ρℓ

realized, at which the state of the substance varies continuously [63].

x ¼ Γs : u ¼ 0, WT ¼ 0; (2)

2 ¼ pℓ þ ρℓ ð Þ uℓ � υsℓ

ℓ ¼ ρsLne

!

� � � � ¼ Lm, 0 þ δ Tm ps

!

kBTm ps � �

m υsℓ

ΔTsℓ Tsℓ

2

� 1

� � � Tm,0 � �, (3)

(4)

,

degree of nonequilibrium of phase transition.

36 Laser Ablation - From Fundamentals to Applications

face of melting x = Гsℓ(t) can be represented as:

� � � � þ ΔCpsΔTsℓ þ ρsþρℓ

2.1.3. Model of heterogeneous evaporation x = Гlv

2.1.2. Model of heterogeneous (surface) melting: crystallization, õ = Гsℓ(t)

2.1.1. Left fixed boundary, x = Гs

the left fixed boundary:

Lne

m ¼ Lm Tm ps

θ = 6.44 � 10�3 K/atm.

where Lne

In the phase equilibrium state, when the pressure of saturated vapor psat(Tsur) is equal to the external pressure pυ, the parameter M = 0. In the subsonic evaporation regime, when M < 1, the behavior of the interface depends on gas-dynamic perturbations in the vaporized matter flow. The maximum nonequilibrium is determined by the maximum value of mass flow, which is known to be achieved at M = 1, when flow of matter through the interface is maximal, and the recoil pressure at the heated surface is minimal. In this case, such evaporation regime can be realized when the behavior of condensed medium no longer depend on the external gasdynamic problem, which greatly simplifies the description of evaporation process. However, the implementation area for such regime remains uncertain, because due to the nonlinearity of the gas dynamics equations, the transition line between M < 1 and M ≥ 1 depends on desired solution and cannot be determined in advance [64]. Taking into account strong spatiotemporal diversity of processes in the nonequilibrium layer and continual media, KL is usually represented as a strong discontinuity in the gas-dynamic parameters. In such case, the kinetic processes in KL are not explicitly considered, and various phenomenological approaches [65–72] were used to determine boundary conditions on the outer side of the KL, which makes it possible without solving the kinetic problem, to determine the joining conditions under certain assumptions about the type of the nonequilibrium distribution function inside the discontinuity. One of the first works in which intensive evaporation was analyzed on the basis of phenomenological model is the work of Crout [65]. In this work, a nonequilibrium ellipsoidal Maxwellian particle distribution function was used written in an analytical form with the anisotropic in spatial directions longitudinal and perpendicular temperatures. Later papers used Mott-Smith approach [66], which was applied to the structure of the shock waves. With the help of this approach [67, 68] expressions for gas-dynamic parameters at M = 1 were obtained. Later in [69], this approach was extended to the entire range of evaporation 0 ≤ M ≤ 1. In [70, 71], the approximation of the distribution function over the entire evaporation range 0 ≤ M ≤ 1 was carried out by other more complicated expressions that satisfied the additional requirements that were formulated taking into account basic laws of gas dynamics. Following this laws, the used distribution functions must when M ! 1 to provide extreme values for all three complete flows of mass jm, momentum ji and energy je. The model [65] also satisfies these requirements. At the same time, in widely used model [69], the requirement of extrema of all flows at selected point M = 1 is not fulfilled. The calculations showed that total flows, jm, ji, je depending on M have extrema at M = 0.88, 1.18, 1.22, respectively. Nonfulfillment of the requirement of extremes indicates an unsuccessful choice of distribution function for the reverse flow f (�) . The mathematical model of surface evaporation in the Knudsen layer approximation consists of three conservation laws and two additional relations from which the parameters on outer side of the Knudsen layer are determined: temperature Tυ and density ρυ. Velocity uυ and Mach number M are determined from the solution of the equations of gas dynamics:

$$\mathbf{x} = \Gamma\_{l\upsilon}(t): \qquad \begin{aligned} \mathbf{j}\_{l\upsilon}^{\mathfrak{m}} &= \rho\_{\ell}(\mathbf{u}\_{\ell} - \boldsymbol{\nu}\_{l\upsilon}) = \rho\_{\upsilon}(\mathbf{u}\_{\upsilon} - \boldsymbol{\nu}\_{l\upsilon}) \end{aligned} \tag{5}$$

$$\begin{aligned} \boldsymbol{j}\_{l\upsilon}^{\mathfrak{i}} &= \boldsymbol{p}\_{\ell} + \boldsymbol{j}\_{l\upsilon}^{\mathfrak{m}}(\boldsymbol{u}\_{\ell} - \boldsymbol{\nu}\_{l\upsilon}) = \boldsymbol{p}\_{\upsilon} + \boldsymbol{j}\_{l\upsilon}^{\mathfrak{m}}(\boldsymbol{u}\_{\upsilon} - \boldsymbol{\nu}\_{l\upsilon}) \end{aligned} \tag{6}$$

$$\boldsymbol{j}\_{l\upsilon}^{\mathfrak{e}} = -\boldsymbol{j}\_{\ell}^{\mathfrak{T}} + \boldsymbol{j}\_{l\upsilon}^{\mathfrak{m}} \left[ \boldsymbol{H}\_{\ell} + \frac{(\boldsymbol{u}\_{\ell} - \boldsymbol{\nu}\_{l\upsilon})^{2}}{2} \right] - \boldsymbol{\sigma}T^{4} = -\boldsymbol{j}\_{\upsilon}^{\mathfrak{T}} + \boldsymbol{j}\_{l\upsilon}^{\mathfrak{m}} \left[ \boldsymbol{H}\_{\nu} + \frac{(\boldsymbol{u}\_{\upsilon} - \boldsymbol{\nu}\_{l\upsilon})^{2}}{2} \right]$$

where j T ℓ ¼ Wℓ ¼ �λð Þ Tℓ ∂Tℓ ∂x , jT υ ¼ Wυ ¼ �λð Þ Tυ ∂Tυ ∂x . With these expressions taken into account, energy conservation law can be represented as:

$$\mathcal{W}\_{\ell} - \mathcal{W}\_{\upsilon} = \rho\_{\ell} \, \upsilon\_{\ell \upsilon} L\_{\upsilon}^{n \iota} + \sigma T^4,\tag{6}$$

2.1.4. Model of surface condensation

be approximated by the expression [75]:

of gas medium far from it and are extrapolated.

filled with vapor within a condensed medium.

When the inequality psat < p is reached, the direction of phase transition changes, evaporation is replaced by condensation process on the target surface, to which the value M < 0 on outer side of the KL corresponds. Unlike evaporation, surface condensation can also occur in supersonic regime. According to [73] in subsonic condensation regime, gas dynamic quantities Tυ, M, and pυ are related only by one two-parameter dependence pυ/psat = F(T,M), where T=Tυ/Tsur and M are dimensionless parameters of temperature and velocity determined by the state of gas dynamic flow far from KL and are extrapolated from the gaseous medium to the discontinuity surface. The function F(T,M) is determined in advance from the solution of Boltzmann equation, the results of which are then tabulated [74]. The values of the function F(T,M) depend weakly on the parameter T and much stronger on the parameter M. Table values from [74] can

When passing through the point M = �1 to supersonic regime of surface condensation, the boundary conditions change. In this case, all quantities on outer side of KL depend on the state

Thus, the description of the kinetics of heterogeneous evaporation is carried out by two oneparameter dependencies used as boundary conditions (for a known parameter M), and sub-

The greatest difficulty in the continuum approach is the description of homogeneous mechanisms of phase transformations: melting crystallization and evaporation. Homogeneous mechanisms of phase transformations are characterized by nucleation of a new phase in a certain volume of superheated/supercooled matter. Representing them in continuum hydrodynamic models requires considerable additional efforts [57, 76] associated with formation of a cavity

The simplest scheme for simulating volumetric boiling in one-dimensional approximation can be represented by introducing, when certain criteria are satisfied into superheated liquid phase of artificial quasi-nuclei with thickness hi(t) bounded by moving planes of liquid-vapor interface xi(t) = Γℓυ , i(t) , xi + 1(t) = Γℓυ , i + 1(t) where i = 1, 2, … is quasi-nuclei number. The criterion for nucleation beginning is a moment when superheating limit temperature Tℓ,max of liquid phase is reached and spatial coordinate. First quasi-nucleus of vapor phase with initial width Δxυ 1(t) = Γℓυ , 2(t) � Γℓυ , 1(t) ~ 5 nm is placed at this point. For a new region of vapor bounded by two flat surfaces, the initial conditions Tυ , 1 = Tℓ , max , ρυ , i = ρsat(Tℓ , max) , pυ , 1 = psat(Tℓ , max) were set. Under the influence of pressure difference between the inside quasi-nuclei and on the external irradiated surface of the liquid layer with the temperature Tsur, a rapidly expanding cavity filled with vapor is formed and blowout the liquid layer with thickness Δxℓ 1(t) = Γℓυ , 1(t) � Γℓυ(t) in the direction of gaseous medium. The description of the processes in the expanding cavity and in the near-surface liquid layer separated from main target was carried out using

sonic surface condensation is described by only one two-parameter (M,T) dependence.

2.1.5. Model of volumetric boiling of liquid phase heated by a laser pulse

pυ=psat ¼ F Tð Þ ; M ≈ F Mð Þ¼ 0:95 � exp ð Þ �2:42 M (8)

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

http://dx.doi.org/10.5772/intechopen.70773

39

where Lne υ ¼ Le υð Þþ Tℓ Cpυ ð Þþ Tb � Tℓυ ρℓþρυ ρℓ�ρυ ð Þ uℓ�uυ 2 2 is the nonequilibrium heat of evaporation, σ is Stefan-Boltzmann law constant.

To determine parameters Tυ, ρυ, pυ, a modified Crout model was used [71].

$$T\_v = \alpha\_T(M) T\_{\text{surr}} \qquad \rho\_v = \rho\_{\text{sat}} \alpha\_\rho(M), \tag{7}$$

$$p\_{\text{sat}} = p\_b \exp\left[\frac{L\_v(T)}{k\_\text{\textdegree}} \left(\frac{1}{T\_b} - \frac{1}{T\_{\text{surr}}}\right)\right], \quad \alpha\_T(M) = \frac{2\gamma M^2 \left(m^2 + 0.5\right)^2}{\left(1 + \gamma M^2\right)^2 m^2 t^2},$$

$$\alpha\_\rho(M) = \frac{1}{\exp\left(-m^2\right) + \pi^{1/2} m (1 + \text{erf}(m))} \cdot \frac{\left(1 + \gamma M^2\right) m^2}{\gamma M^2 \left(m^2 + 0.5\right)^2},$$

the value m is determined from the equation F(M)(m2 + 0.5)2 � m2 (m2 + 1.5 + a) = 0, where F Mð Þ¼ 1 þ 3γM2�1 γM2 ð Þ �1 2 , a ¼ 2t 2 � 0:5π1=2mt � 1, t ¼ 2m π1=2 þ 1þerfð Þ m exp �m2 ð Þþπ1=2mð Þ 1þerfð Þ m .

ρsat and psat are the saturated vapor density and pressure, pb is the pressure under normal conditions (1.01325 � 105 Pa), usound = (γRT) 1/2 is the sound velocity, M=uυ/usound is the Mach number, Tb is the boiling temperature, kB is Boltzmann's constant. Taking into account psat, the value of ρsat was determined from the equation p=p(ρ,T). For the laser radiation transfer equation on the target surface, the condition Gð Þ¼ Γlυ ð Þ 1 � R Tð Þ sur G0exp � t τ � �2 � � was used where R(Tsur) is the temperature dependence of reflectivity of target surface.

#### 2.1.4. Model of surface condensation

momentum ji and energy je. The model [65] also satisfies these requirements. At the same time, in widely used model [69], the requirement of extrema of all flows at selected point M = 1 is not fulfilled. The calculations showed that total flows, jm, ji, je depending on M have extrema at M = 0.88, 1.18, 1.22, respectively. Nonfulfillment of the requirement of extremes

ematical model of surface evaporation in the Knudsen layer approximation consists of three conservation laws and two additional relations from which the parameters on outer side of the Knudsen layer are determined: temperature Tυ and density ρυ. Velocity uυ and Mach

ℓυ ¼ ρℓ ð Þ¼ uℓ � υℓυ ρυð Þ uυ � υℓυ

ℓυð Þ¼ uℓ � υℓυ pυ þ j

υ ¼ Wυ ¼ �λð Þ Tυ ∂Tυ

Wℓ � Wυ ¼ ρℓ υℓυ Lne

ð Þ uℓ�uυ 2

ρℓ�ρυ

To determine parameters Tυ, ρυ, pυ, a modified Crout model was used [71].

1 Tb � 1 Tsur

� � � �

� σT4 ¼ �j

m

T υ þ j m

υ þ σT4

, αΤð Þ¼ M

exp �m2 ð Þþ π1=2mð Þ 1 þ erfð Þ m � 1 þ γM2 � �m2

ℓυð Þ uυ � υℓυ

Tυ ¼ αΤð Þ M Tsur, ρυ ¼ ρsatαρð Þ M , (7)

γM2

π1=2 þ 1þerfð Þ m

ℓυ Hυ þ ð Þ uυ � υlυ 2

2 is the nonequilibrium heat of evaporation, σ

2γM2 m2 þ 0:5 � �2 1 þ γM2 � �2

ð Þ m2 þ 0:5 2 ,

exp �m2 ð Þþπ1=2mð Þ 1þerfð Þ m .

1/2 is the sound velocity, M=uυ/usound is the Mach

m2t 2 ,

(m2 + 1.5 + a) = 0, where

τ � �2 � �

was used

2

, (6)

∂x . With these expressions taken into

" # (5)

(�)

. The math-

indicates an unsuccessful choice of distribution function for the reverse flow f

number M are determined from the solution of the equations of gas dynamics:

m

2

m

ℓυ ¼ pℓ þ j

ℓυ Hℓ þ ð Þ uℓ � υℓυ 2

" #

x ¼ ΓlυðtÞ : j

38 Laser Ablation - From Fundamentals to Applications

j e ℓυ ¼ �j T ℓ þ j m

υ ¼ Le

F Mð Þ¼ 1 þ 3γM2�1

ℓ ¼ Wℓ ¼ �λð Þ Tℓ

is Stefan-Boltzmann law constant.

psat ¼ pbexp

γM2 ð Þ �1 2 , a ¼ 2t

conditions (1.01325 � 105 Pa), usound = (γRT)

where j T

where Lne

j i

> ∂Tℓ ∂x , jT

account, energy conservation law can be represented as:

υð Þþ Tℓ Cpυ ð Þþ Tb � Tℓυ ρℓþρυ

Lυð Þ T kB

αρð Þ¼ M 1

the value m is determined from the equation F(M)(m2 + 0.5)2 � m2

2 � 0:5π1=2mt � 1, t ¼ 2m

equation on the target surface, the condition Gð Þ¼ Γlυ ð Þ 1 � R Tð Þ sur G0exp � t

where R(Tsur) is the temperature dependence of reflectivity of target surface.

ρsat and psat are the saturated vapor density and pressure, pb is the pressure under normal

number, Tb is the boiling temperature, kB is Boltzmann's constant. Taking into account psat, the value of ρsat was determined from the equation p=p(ρ,T). For the laser radiation transfer When the inequality psat < p is reached, the direction of phase transition changes, evaporation is replaced by condensation process on the target surface, to which the value M < 0 on outer side of the KL corresponds. Unlike evaporation, surface condensation can also occur in supersonic regime. According to [73] in subsonic condensation regime, gas dynamic quantities Tυ, M, and pυ are related only by one two-parameter dependence pυ/psat = F(T,M), where T=Tυ/Tsur and M are dimensionless parameters of temperature and velocity determined by the state of gas dynamic flow far from KL and are extrapolated from the gaseous medium to the discontinuity surface. The function F(T,M) is determined in advance from the solution of Boltzmann equation, the results of which are then tabulated [74]. The values of the function F(T,M) depend weakly on the parameter T and much stronger on the parameter M. Table values from [74] can be approximated by the expression [75]:

$$p\_v/p\_{sat} = F(T,\ M) \approx F(M) = 0.95 \times \exp\left(-2.42\text{ M}\right) \tag{8}$$

When passing through the point M = �1 to supersonic regime of surface condensation, the boundary conditions change. In this case, all quantities on outer side of KL depend on the state of gas medium far from it and are extrapolated.

Thus, the description of the kinetics of heterogeneous evaporation is carried out by two oneparameter dependencies used as boundary conditions (for a known parameter M), and subsonic surface condensation is described by only one two-parameter (M,T) dependence.

#### 2.1.5. Model of volumetric boiling of liquid phase heated by a laser pulse

The greatest difficulty in the continuum approach is the description of homogeneous mechanisms of phase transformations: melting crystallization and evaporation. Homogeneous mechanisms of phase transformations are characterized by nucleation of a new phase in a certain volume of superheated/supercooled matter. Representing them in continuum hydrodynamic models requires considerable additional efforts [57, 76] associated with formation of a cavity filled with vapor within a condensed medium.

The simplest scheme for simulating volumetric boiling in one-dimensional approximation can be represented by introducing, when certain criteria are satisfied into superheated liquid phase of artificial quasi-nuclei with thickness hi(t) bounded by moving planes of liquid-vapor interface xi(t) = Γℓυ , i(t) , xi + 1(t) = Γℓυ , i + 1(t) where i = 1, 2, … is quasi-nuclei number. The criterion for nucleation beginning is a moment when superheating limit temperature Tℓ,max of liquid phase is reached and spatial coordinate. First quasi-nucleus of vapor phase with initial width Δxυ 1(t) = Γℓυ , 2(t) � Γℓυ , 1(t) ~ 5 nm is placed at this point. For a new region of vapor bounded by two flat surfaces, the initial conditions Tυ , 1 = Tℓ , max , ρυ , i = ρsat(Tℓ , max) , pυ , 1 = psat(Tℓ , max) were set. Under the influence of pressure difference between the inside quasi-nuclei and on the external irradiated surface of the liquid layer with the temperature Tsur, a rapidly expanding cavity filled with vapor is formed and blowout the liquid layer with thickness Δxℓ 1(t) = Γℓυ , 1(t) � Γℓυ(t) in the direction of gaseous medium. The description of the processes in the expanding cavity and in the near-surface liquid layer separated from main target was carried out using gas-hydrodynamic system of Eq. (1). The model of heterogeneous evaporation (5)–(7) was used as boundary conditions on interphase planes Γℓυ , i(t) , Γℓυ , i + 1(t). Also, absorption and reflection of the intensity of laser radiation in the formed layer were taken into account. In case of repeated explosive boiling with the formation of next liquid fragment for each of them, as well as the main target, gaseous medium and formed vapor cavities, the solution algorithm remains unified.

2.2. Computational algorithm

coefficient of κL ≥ 107 m1

region.

The differential model (1)–(8) was approximated by a family of conservative finite-difference schemes [80, 81] written on computational grids with dynamic adaptation [79, 82, 83]. The method is based on the idea of transition to an arbitrary nonstationary coordinate system that allows calculations with an arbitrary number of discontinuous solutions, such as shock waves, propagating phase and temperature fronts, contact boundaries, and spalled fragments.

In the mathematical modeling of the process of laser action on metals with an absorption

model, the maximum of the temperature profile Tmax coincides with surface temperature Tsur from which evaporation occurs. At lower values of κL typical for dielectrics as well as for metals when irradiated with an electron beam or X-ray pulsed radiation or in the case of a metaldielectric transition [84, 85], bulk heating of the substance and surface evaporation leads to the formation of temperature maximum below the surface of the substance. As a result, a volumetric explosive boiling of superheated liquid can occur in the region of the temperature maximum. The explosive boiling of superheated liquid is closely related to the concept of metal-dielectric phase transition. Zeldovich and Landau [86] denoted the possibility of a metal-dielectric phase transition for an expanded metal at subcritical temperatures. To perform the metal-dielectric transition, it is necessary that the energy acquired by metal atoms is higher than their binding energy in the crystal lattice and the distance between atoms is equal to the value causing a violation of their short-range order, which leads to the localization of the electrons on atoms [87]. With the advent of lasers capable of evaporating metals and producing plasma on their surfaces, it has become possible to observe laser-induced phase transitions metal-dielectric with the formation of transparency waves in the massive targets [88] and thin metallic films [89] in the subcritical temperature range. In other experimental and theoretical studies of the

results of the appearance of metal-dielectric transitions accompanied by the formation of

Nevertheless, the physical mechanisms of interaction of laser radiation with metals taking into account metal-insulator transition have not yet been fully studied. This causes great difficulties in determining thermophysical and optical characteristics of metals in the vicinity of the critical

Figure 2 shows temperature dependences of thermophysical λ(T), Lv(T), Cp(T) and optical R(T), κL(T) characteristics of aluminum. The curves were Lv(T) and Cp(T) obtained from molecular dynamics calculations, and λ(T), κL(T), and R(T) were constructed on the basis of theoret-

, the model of surface heating and evaporation is widely used. In this

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

http://dx.doi.org/10.5772/intechopen.70773

41

–108 W/cm2 with metallic targets, the

3. Temperature dependences of aluminum properties

interaction of ns-laser pulses with an intensity of 107

transparency waves are reported [53, 54, 85, 88–91].

ical concepts [92–94] and reference data [95–97].

The values of limit of superheat temperature depend on the rate of energy input and can be determined in advance from the molecular dynamics simulation. For Al, depending on the rate of heating, limit of superheat temperature of liquid phase Tmax is in the range Tmax ~ (0.89 ÷ 0.95)Tcr [77].

It should be noted that the use of gas dynamic equations inside the cavities which dimensions are small in the initial moments of time in comparison with the mean free path is a simplifying approximation that makes it easy to take into account the influence of expansion velocity of the cavity.

#### 2.1.6. Moving contact boundary, x = Гvg(t)

On contact boundary and front of the shock wave, well-known standard relations are used [78]. At vapor-air interface, the boundary conditions were set in the form of equal values of velocity, pressure, and temperature:

$$\mathbf{x} = \Gamma\_{l\upsilon}(t) : \quad \mathbf{u}\_{\upsilon} = \mathbf{u}\_{\mathfrak{F}'} \; \; P\_{\upsilon} = P\_{\mathfrak{g}'} \; \; T\_{\upsilon} = T\_{\mathfrak{g}}$$

#### 2.1.7. Moving shock wave, x = Гsh,g(t)

A shock wave in air Гsh,g(t) is a strong nonstationary discontinuity, on which three conservation laws are written in laboratory coordinate system [78]:

$$\begin{aligned} f\_{slt,g}^{\mathfrak{m}} &= \rho\_1 \left( \mu\_1 - \upsilon\_{slt,g} \right) = \rho\_0 \left( \mu\_0 - \upsilon\_{slt,g} \right), \\ f\_{slt,g}^{\mathfrak{j}} &= p\_1 + \rho\_1 \left( \mu\_1 - \upsilon\_{slt,g} \right)^2 = P\_0 + \rho\_0 \left( \mu\_0 - \upsilon\_{slt,g} \right)^2, \\ f\_{slt,g}^{\mathfrak{m}} &\left( \varepsilon\_1 + \frac{\left( \mu\_1 - \upsilon\_{slt,g} \right)^2}{2} \right) - W\_{T,1} = f\_{slt,g}^{\mathfrak{m}} \left( \varepsilon\_0 + \frac{\left( \mu\_0 - \upsilon\_{slt,g} \right)^2}{2} \right) - W\_{T,0} \end{aligned} \tag{9}$$

The indices 0 and 1 denote the values of the quantities on side of the background and the shock wave, respectively.

#### 2.1.8. Right moving boundary, x = Гg(t)

The right boundary on the side of unperturbed gas is declared moving in order to improve economic effectiveness of the computational algorithm [79]. The speed of its motion is found from the differential equation of momentum.

#### 2.2. Computational algorithm

gas-hydrodynamic system of Eq. (1). The model of heterogeneous evaporation (5)–(7) was used as boundary conditions on interphase planes Γℓυ , i(t) , Γℓυ , i + 1(t). Also, absorption and reflection of the intensity of laser radiation in the formed layer were taken into account. In case of repeated explosive boiling with the formation of next liquid fragment for each of them, as well as the main target, gaseous medium and formed vapor cavities, the solution algorithm

The values of limit of superheat temperature depend on the rate of energy input and can be determined in advance from the molecular dynamics simulation. For Al, depending on the rate of heating, limit of superheat temperature of liquid phase Tmax is in the range Tmax ~

It should be noted that the use of gas dynamic equations inside the cavities which dimensions are small in the initial moments of time in comparison with the mean free path is a simplifying approximation that makes it easy to take into account the influence of expansion velocity of the

On contact boundary and front of the shock wave, well-known standard relations are used [78]. At vapor-air interface, the boundary conditions were set in the form of equal values of

x ¼ ΓlυðtÞ : uυ ¼ ug, Pυ ¼ Pg, Tυ ¼ Tg

A shock wave in air Гsh,g(t) is a strong nonstationary discontinuity, on which three conservation

� �2

,

!

u0 � υsh,g � �2 2

� WT, 0

(9)

� �,

� �2 ¼ P0 þ ρ0 u0 � υsh,g

� WT, 1 ¼ j

m sh,g ε0 þ

The indices 0 and 1 denote the values of the quantities on side of the background and the shock

The right boundary on the side of unperturbed gas is declared moving in order to improve economic effectiveness of the computational algorithm [79]. The speed of its motion is found

remains unified.

(0.89 ÷ 0.95)Tcr [77].

2.1.6. Moving contact boundary, x = Гvg(t)

40 Laser Ablation - From Fundamentals to Applications

velocity, pressure, and temperature:

2.1.7. Moving shock wave, x = Гsh,g(t)

j m

j i

j m sh, g ε1 þ

2.1.8. Right moving boundary, x = Гg(t)

from the differential equation of momentum.

wave, respectively.

laws are written in laboratory coordinate system [78]:

sh, g ¼ p1 þ ρ1 u1 � υsh, g

u1 � υsh, g � �2 2

!

� � ¼ ρ0 u0 � υsh, g

sh, g ¼ ρ1 u1 � υsh,g

cavity.

The differential model (1)–(8) was approximated by a family of conservative finite-difference schemes [80, 81] written on computational grids with dynamic adaptation [79, 82, 83]. The method is based on the idea of transition to an arbitrary nonstationary coordinate system that allows calculations with an arbitrary number of discontinuous solutions, such as shock waves, propagating phase and temperature fronts, contact boundaries, and spalled fragments.

## 3. Temperature dependences of aluminum properties

In the mathematical modeling of the process of laser action on metals with an absorption coefficient of κL ≥ 107 m1 , the model of surface heating and evaporation is widely used. In this model, the maximum of the temperature profile Tmax coincides with surface temperature Tsur from which evaporation occurs. At lower values of κL typical for dielectrics as well as for metals when irradiated with an electron beam or X-ray pulsed radiation or in the case of a metaldielectric transition [84, 85], bulk heating of the substance and surface evaporation leads to the formation of temperature maximum below the surface of the substance. As a result, a volumetric explosive boiling of superheated liquid can occur in the region of the temperature maximum. The explosive boiling of superheated liquid is closely related to the concept of metal-dielectric phase transition. Zeldovich and Landau [86] denoted the possibility of a metal-dielectric phase transition for an expanded metal at subcritical temperatures. To perform the metal-dielectric transition, it is necessary that the energy acquired by metal atoms is higher than their binding energy in the crystal lattice and the distance between atoms is equal to the value causing a violation of their short-range order, which leads to the localization of the electrons on atoms [87].

With the advent of lasers capable of evaporating metals and producing plasma on their surfaces, it has become possible to observe laser-induced phase transitions metal-dielectric with the formation of transparency waves in the massive targets [88] and thin metallic films [89] in the subcritical temperature range. In other experimental and theoretical studies of the interaction of ns-laser pulses with an intensity of 107 –108 W/cm2 with metallic targets, the results of the appearance of metal-dielectric transitions accompanied by the formation of transparency waves are reported [53, 54, 85, 88–91].

Nevertheless, the physical mechanisms of interaction of laser radiation with metals taking into account metal-insulator transition have not yet been fully studied. This causes great difficulties in determining thermophysical and optical characteristics of metals in the vicinity of the critical region.

Figure 2 shows temperature dependences of thermophysical λ(T), Lv(T), Cp(T) and optical R(T), κL(T) characteristics of aluminum. The curves were Lv(T) and Cp(T) obtained from molecular dynamics calculations, and λ(T), κL(T), and R(T) were constructed on the basis of theoretical concepts [92–94] and reference data [95–97].

interpreted as a surface spitting of liquid phase when a critical temperature is reached on the surface [99]. This interpretation is not convincing, since it does not allow to determine even

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

Let us consider ns-PLA process of aluminum target under the mode of exposure close to experimental conditions [53]. Al target with the thickness of 104 m is placed in the air under the normal conditions with pressure 1 bar and room temperature (T = 300 K). Laser pulse of

G0 = 6.1 108 W/cm2 is the peak intensity at t = 0, ∞ < t < ∞, β = 4 ln2 falls on the target surface from right to left. Air for the selected exposure mode is completely transparent for laser radiation, which is partially reflected from metal surface and partially absorbed by target

At the initial stage, time evolution of the processes at the surface of the target, in the target and in the gas medium near the target such as appearance of phase fronts (melting, evaporation), contact boundary, and shock wave and associated with formation of new phase media (liquid, vapor) occurs at the leading edge of laser pulse. Heterogeneous melting begins at the moment

interior of the target is heated, the melting front Гsl(t) runs from its surface forming a new region of liquid phase. Further heating leads to appearance of heterogeneous evaporation front Гlv(t) that runs inside the melt. A flow of the evaporated matter was formed at the surface of the melt pushing out air and creating another new phase—vapor. The new area occupied by vapor is limited on the one hand by the moving interface Гlv(t) (evaporating surface) and on the other hand by the moving contact boundary vapor-air Гvg(t). The surface evaporation process is controlled by the surface temperature Tsur(t) and Mach number M(t) on the outer side of KL (0 < M(t) ≤ 1). Evaporation begins at the moment t = 2.3 109 s when the saturated vapor pressure exceeds the pressure of external gas psat(t) > pg(t). The maximum velocity of evaporation front υℓv = 88 ms1 is approximately 1.5 times smaller than velocity of melting front. The vaporized matter flow acting as a piston pushes out cold air, and performing a certain work warms up at peak G0 intensity up to temperature Tvap = 4.7 103 K. Under the pushing action of vapor flow, the compression of cold dense air which turns into a shock wave occurs, t = 2.1 109 s. The spatial structure of erosion plume is shown in Figure 3 by the main

By the moment t = 0, the shock wave propagates with velocity υsh,g = 2.9 km s1 and temperature Tsh,g = 4.4 103 K, ahead of contact boundary moving with speed υsh,g = 2.5 km s1 toward the laser flow. The temperatures of vapor and air are in this case insufficient for

initiation of ionization, and vapor-gas medium remains transparent for laser radiation.

t = 3.9 109 s. The maximum velocity of melting front reaches a value υsl = 130 ms1

) with full width at half maximum (FWHM)

http://dx.doi.org/10.5772/intechopen.70773

, where

43

. As the

2

τ = 5 109 s, with wavelength λL = 1.06 μm and fluence F = 3.5 J/cm2

1/2 G0 exp(β (t/τ)

material layer. The energy release of laser pulse has volumetric nature.

characteristics T(x), ρ(x), u(x), p(x) after the formation of shock wave.

approximately the parameters of explosive boiling.

4.1. Mode of exposure

Gaussian shape G = (β/π)

4.2. Evolution of processes

Figure 2. Temperature dependences of Al properties (a–d).

#### 4. Results discussion

One of the purposes of this work is a detailed study of mechanism of explosive boiling in metals since explosive boiling is considered to be the most efficient thermal mechanism for the laser ablation of materials. At the same time, the difficulties associated with understanding of the mechanism of homogeneous phase transitions in metals occurring under the action of ns-laser pulses are known. In [98], based on analytical solution of thermal model, it was established that laser removal of material from a solid target, with a certain choice of irradiation parameters and material, is determined by the presence of near-surface temperature maximum. The validity of this position was confirmed by the results of numerical solution of thermal model for low-absorbing liquids irradiated by laser pulses [43] and nonmetallic solid materials [47]. For strongly absorbing media mainly metallic, calculations based on thermal model [44, 45, 47, 56] have shown that the magnitude of near-surface temperature maximum is several degrees. On this basis, overheating was excluded from consideration, up to the statement [56] that in metals, maximum temperature is always on the surface of the target, and sub-surface superheating is impossible. Explosive boiling was

interpreted as a surface spitting of liquid phase when a critical temperature is reached on the surface [99]. This interpretation is not convincing, since it does not allow to determine even approximately the parameters of explosive boiling.

#### 4.1. Mode of exposure

Let us consider ns-PLA process of aluminum target under the mode of exposure close to experimental conditions [53]. Al target with the thickness of 104 m is placed in the air under the normal conditions with pressure 1 bar and room temperature (T = 300 K). Laser pulse of Gaussian shape G = (β/π) 1/2 G0 exp(β (t/τ) 2 ) with full width at half maximum (FWHM) τ = 5 109 s, with wavelength λL = 1.06 μm and fluence F = 3.5 J/cm2 , where G0 = 6.1 108 W/cm2 is the peak intensity at t = 0, ∞ < t < ∞, β = 4 ln2 falls on the target surface from right to left. Air for the selected exposure mode is completely transparent for laser radiation, which is partially reflected from metal surface and partially absorbed by target material layer. The energy release of laser pulse has volumetric nature.

#### 4.2. Evolution of processes

4. Results discussion

Figure 2. Temperature dependences of Al properties (a–d).

42 Laser Ablation - From Fundamentals to Applications

One of the purposes of this work is a detailed study of mechanism of explosive boiling in metals since explosive boiling is considered to be the most efficient thermal mechanism for the laser ablation of materials. At the same time, the difficulties associated with understanding of the mechanism of homogeneous phase transitions in metals occurring under the action of ns-laser pulses are known. In [98], based on analytical solution of thermal model, it was established that laser removal of material from a solid target, with a certain choice of irradiation parameters and material, is determined by the presence of near-surface temperature maximum. The validity of this position was confirmed by the results of numerical solution of thermal model for low-absorbing liquids irradiated by laser pulses [43] and nonmetallic solid materials [47]. For strongly absorbing media mainly metallic, calculations based on thermal model [44, 45, 47, 56] have shown that the magnitude of near-surface temperature maximum is several degrees. On this basis, overheating was excluded from consideration, up to the statement [56] that in metals, maximum temperature is always on the surface of the target, and sub-surface superheating is impossible. Explosive boiling was At the initial stage, time evolution of the processes at the surface of the target, in the target and in the gas medium near the target such as appearance of phase fronts (melting, evaporation), contact boundary, and shock wave and associated with formation of new phase media (liquid, vapor) occurs at the leading edge of laser pulse. Heterogeneous melting begins at the moment t = 3.9 109 s. The maximum velocity of melting front reaches a value υsl = 130 ms1 . As the interior of the target is heated, the melting front Гsl(t) runs from its surface forming a new region of liquid phase. Further heating leads to appearance of heterogeneous evaporation front Гlv(t) that runs inside the melt. A flow of the evaporated matter was formed at the surface of the melt pushing out air and creating another new phase—vapor. The new area occupied by vapor is limited on the one hand by the moving interface Гlv(t) (evaporating surface) and on the other hand by the moving contact boundary vapor-air Гvg(t). The surface evaporation process is controlled by the surface temperature Tsur(t) and Mach number M(t) on the outer side of KL (0 < M(t) ≤ 1). Evaporation begins at the moment t = 2.3 109 s when the saturated vapor pressure exceeds the pressure of external gas psat(t) > pg(t). The maximum velocity of evaporation front υℓv = 88 ms1 is approximately 1.5 times smaller than velocity of melting front. The vaporized matter flow acting as a piston pushes out cold air, and performing a certain work warms up at peak G0 intensity up to temperature Tvap = 4.7 103 K. Under the pushing action of vapor flow, the compression of cold dense air which turns into a shock wave occurs, t = 2.1 109 s. The spatial structure of erosion plume is shown in Figure 3 by the main characteristics T(x), ρ(x), u(x), p(x) after the formation of shock wave.

By the moment t = 0, the shock wave propagates with velocity υsh,g = 2.9 km s1 and temperature Tsh,g = 4.4 103 K, ahead of contact boundary moving with speed υsh,g = 2.5 km s1 toward the laser flow. The temperatures of vapor and air are in this case insufficient for initiation of ionization, and vapor-gas medium remains transparent for laser radiation.

where the deepest passing into the region of superheated metastable liquid critical point of

Figure 4. Spatial profiles of temperature maximum, (fragments) before explosive boiling: (a) continual model and (b)

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

http://dx.doi.org/10.5772/intechopen.70773

45

The process of explosive boiling begins at the backside front of the laser pulse. First boiling occurs at the moment t = +1.5 ns when the maximum of the temperature profile reaches the temperature of limiting superheating of liquid phase equal to Tℓ,max = 0.9Tcr = 6840 K. Figure 5a–d shows the spatial profiles of T(x), ρ(x), u(x), p(x). Arising cavity is the result of homogeneous nucleation in superheated liquid phase and is accompanied by sharp increase in pressure pmax = psat(Tℓ,max) which is more than two times greater than the recoil pressure on irradiated free surface psur = 0.55psat(Tsur). This pressure plays a decisive role in the expansion

Figure 5. The spatial profiles of (a) T(x), (b) ρ(x), (c) u(x), (d) p(x) in the target and the gas medium at the moment

liquid-vapor transition is achieved.

4.4. Explosive boiling

atomistic model.

t = 1.7 ns.

Figure 3. The spatial profiles of (a) T(x), (b) ρ(x), (c) u(x), and (d) p(x) in the target and the gas medium at the moment t = �1.9 ns.

#### 4.3. Formation of near-surface temperature maximum

A consistent study of laser ablation is complicated by the fact that heterogeneous and homogeneous mechanisms of both melting and evaporation prove to be interrelated, and this interaction must be taken into account explicitly. Pulsed laser action on materials (including metals) has volumetric nature of energy release. Thus, when target is heated strongly in the vicinity of the critical point (Tcr = 7600 K, ρcr = 0.47 g cm�3 , pcr = 1.42 kbar), thermophysical and optical properties of liquid phase change abruptly. Heat capacity Cp when approaching a critical point tends to infinity Cp ! ∞, thermal conductivity λl and absorption κL coefficients tend to zero, surface absorptivity A = (1�R) ! 1, Figure 2a–c. Under the joint influence of volumetric energy release of laser pulse and energy transfer by surface flow of evaporating substance, a near-surface maximum of temperature Tℓ,max is formed in the depth of liquid phase, for which the following relations are satisfied: Tℓ,max > Tsur, pℓ,max,sat > psur,sat.

Figure 4a shows a fragment of the spatial temperature distribution in near-surface layer of the target with overheating Tℓ,max � Tsur = 170 K at a depth of 30 nm from the irradiated surface at the moment immediately preceding the explosive boiling. For comparison, Figure 4b shows a fragment of temperature profile, calculated under the same conditions with help of molecular dynamics, which showed close results. The maximum was located at a depth of 70 nm, and the superheat value was 150 K. Thus, simulation results obtained by different methods indicate the presence of temperature inhomogeneity in metallic target (Al) caused by a decrease in the temperature of the irradiated surface by heterogeneous evaporation. Thus, the conditions for phase explosion were created with formation of a cavity in the region of maximum temperature, Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling http://dx.doi.org/10.5772/intechopen.70773 45

Figure 4. Spatial profiles of temperature maximum, (fragments) before explosive boiling: (a) continual model and (b) atomistic model.

where the deepest passing into the region of superheated metastable liquid critical point of liquid-vapor transition is achieved.

#### 4.4. Explosive boiling

4.3. Formation of near-surface temperature maximum

44 Laser Ablation - From Fundamentals to Applications

t = �1.9 ns.

vicinity of the critical point (Tcr = 7600 K, ρcr = 0.47 g cm�3

A consistent study of laser ablation is complicated by the fact that heterogeneous and homogeneous mechanisms of both melting and evaporation prove to be interrelated, and this interaction must be taken into account explicitly. Pulsed laser action on materials (including metals) has volumetric nature of energy release. Thus, when target is heated strongly in the

Figure 3. The spatial profiles of (a) T(x), (b) ρ(x), (c) u(x), and (d) p(x) in the target and the gas medium at the moment

optical properties of liquid phase change abruptly. Heat capacity Cp when approaching a critical point tends to infinity Cp ! ∞, thermal conductivity λl and absorption κL coefficients tend to zero, surface absorptivity A = (1�R) ! 1, Figure 2a–c. Under the joint influence of volumetric energy release of laser pulse and energy transfer by surface flow of evaporating substance, a near-surface maximum of temperature Tℓ,max is formed in the depth of liquid

Figure 4a shows a fragment of the spatial temperature distribution in near-surface layer of the target with overheating Tℓ,max � Tsur = 170 K at a depth of 30 nm from the irradiated surface at the moment immediately preceding the explosive boiling. For comparison, Figure 4b shows a fragment of temperature profile, calculated under the same conditions with help of molecular dynamics, which showed close results. The maximum was located at a depth of 70 nm, and the superheat value was 150 K. Thus, simulation results obtained by different methods indicate the presence of temperature inhomogeneity in metallic target (Al) caused by a decrease in the temperature of the irradiated surface by heterogeneous evaporation. Thus, the conditions for phase explosion were created with formation of a cavity in the region of maximum temperature,

phase, for which the following relations are satisfied: Tℓ,max > Tsur, pℓ,max,sat > psur,sat.

, pcr = 1.42 kbar), thermophysical and

The process of explosive boiling begins at the backside front of the laser pulse. First boiling occurs at the moment t = +1.5 ns when the maximum of the temperature profile reaches the temperature of limiting superheating of liquid phase equal to Tℓ,max = 0.9Tcr = 6840 K. Figure 5a–d shows the spatial profiles of T(x), ρ(x), u(x), p(x). Arising cavity is the result of homogeneous nucleation in superheated liquid phase and is accompanied by sharp increase in pressure pmax = psat(Tℓ,max) which is more than two times greater than the recoil pressure on irradiated free surface psur = 0.55psat(Tsur). This pressure plays a decisive role in the expansion

Figure 5. The spatial profiles of (a) T(x), (b) ρ(x), (c) u(x), (d) p(x) in the target and the gas medium at the moment t = 1.7 ns.

1.7 μm, the velocity increased up to 2.3 km/s and density due to cooling changed from 0.65 g/cm3

fragment on the verge of extinction due to surface evaporation fell to hℓ = 1.9 nm, and the

The simulation results are in good qualitative agreement with the experimental data [53]. Quantitatively, the total depth due to release of superheated liquid phase (dΣ = 37.0 nm) differs by about two times in excess. The total removal of matter by the mechanism of explosive boiling considerably exceeds the amount of evaporated matter. The trajectory of shock wave motion almost completely coincides with the experimental dependence [53] up to the moment t = 200 ns. From now, calculated curve is located above the experimental curve. As the fluence

move to leading edge of the laser pulse. The frequency and amount of boiling thus increase to tens and hundreds of ejected liquid fragments, but the overall picture of processes does not

Within the framework of continuum approach, a multiphase, multifront hydrodynamic model and computational algorithm were constructed that were designed for modeling ns-PLA of metal targets embedded in gaseous media. The model contains temperature dependences of optical, thermal, and thermodynamic properties of metal target (Al). The temperature dependences of heat capacity Cp(T) and thermal conductivity λ(T) of aluminum take into account their singularity in the vicinity of critical point. In optical characteristics, the metal-insulator phase transition was taken into account accompanied by the formation of the transparency wave.

The computational algorithm is based on the method of dynamic adaptation, which ideally fit to the study of problems with heterogeneous phase transitions and allows the production of numerical solutions with an explicit allocation of unlimited number of interface boundaries in

The calculations have shown that under the influence of temperature dependences of optical and thermophysical properties of Al the nature of laser heating of the target in hightemperature region changes substantially. The energy release has a volume nature, and irradiated surface is markedly cooled by the process of surface evaporation. An inhomogeneity arises in the spatial temperature profile—a near-surface maximum of temperature. When the limiting superheating temperature is reached, the conditions for explosive boiling were realized with formation of a cavity in the region of the maximum, where the maximum penetration into the region of superheated metastable liquid near the critical point of liquid-vapor transition is achieved. For Al in action modes under the consideration, the typical depths of temperature maximum are d = 20–70 nm with the values ΔTmax = Tℓ,max Tsur = 150–170 K. An approach was proposed to formulation in continuum models of homogeneous nucleation in a superheated liquid phase based on the introduction into hydrodynamic model of quasi-nuclei

condensed medium, contact boundaries, and shock waves in a gaseous medium.

. At the same time, the pressure in the cavity decreased to 200 bar, the thickness of

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

http://dx.doi.org/10.5772/intechopen.70773

47

, the target is heated faster and explosive boiling process begins to

to 1.3 g/cm3

temperature decreased to Tℓ = 5400 K.

increases up to F ~ 10 J/cm2

undergo qualitative changes.

5. Conclusion

Table 1. Number, time of appearance and thickness of spalled fragments.

of cavity filled with vapor and in rapid growth of flyout velocity of thin d1 = 30 nm fragment of liquid metal. Because of small thickness, flying fragment weakly absorbs laser radiation, the main part of which is released in the target increasing its temperature. Both surfaces of spalled fragment and surface of the target are subject to intense surface evaporation. Heterogeneous evaporation leads to decrease in temperature of the fragment and causes strong cooling of target surface, contributing to formation of next near-surface temperature maximum in which a new explosive boiling occurs when the temperature of limit overheating Tℓ,max = 0.9Tcr of metastable liquid is reached. Repeated boiling occurred at a time t = +1.7 ns. The total number of explosive boilings reaches 5. The moments of occurrence and thickness of outgoing fragments are given in Table 1.

Figure 6a–d shows the state of all liquid fragments of explosive boiling at the time t = +3.0 ns. The resulting fragments eventually acquire a sufficiently high flyout velocity and under the influence of laser radiation absorption continue to evaporate until they disappear completely. The general trend in the evolution of all fragments is that over time the cavity size, flyout velocity and density of matter increase, and pressure in the cavity, thickness and temperature of the fragments decrease. So in the initially formed cavity the width increased from 5.0 nm to

Figure 6. The spatial profiles of (a) T(x), (b) ρ(x), (c) u(x), (d) p(x) in the target and the gas medium after five explosive boiling at the moment t = 3.0 ns.

1.7 μm, the velocity increased up to 2.3 km/s and density due to cooling changed from 0.65 g/cm3 to 1.3 g/cm3 . At the same time, the pressure in the cavity decreased to 200 bar, the thickness of fragment on the verge of extinction due to surface evaporation fell to hℓ = 1.9 nm, and the temperature decreased to Tℓ = 5400 K.

The simulation results are in good qualitative agreement with the experimental data [53]. Quantitatively, the total depth due to release of superheated liquid phase (dΣ = 37.0 nm) differs by about two times in excess. The total removal of matter by the mechanism of explosive boiling considerably exceeds the amount of evaporated matter. The trajectory of shock wave motion almost completely coincides with the experimental dependence [53] up to the moment t = 200 ns. From now, calculated curve is located above the experimental curve. As the fluence increases up to F ~ 10 J/cm2 , the target is heated faster and explosive boiling process begins to move to leading edge of the laser pulse. The frequency and amount of boiling thus increase to tens and hundreds of ejected liquid fragments, but the overall picture of processes does not undergo qualitative changes.

## 5. Conclusion

of cavity filled with vapor and in rapid growth of flyout velocity of thin d1 = 30 nm fragment of liquid metal. Because of small thickness, flying fragment weakly absorbs laser radiation, the main part of which is released in the target increasing its temperature. Both surfaces of spalled fragment and surface of the target are subject to intense surface evaporation. Heterogeneous evaporation leads to decrease in temperature of the fragment and causes strong cooling of target surface, contributing to formation of next near-surface temperature maximum in which a new explosive boiling occurs when the temperature of limit overheating Tℓ,max = 0.9Tcr of metastable liquid is reached. Repeated boiling occurred at a time t = +1.7 ns. The total number of explosive boilings reaches 5. The moments of occurrence and thickness of outgoing frag-

Fragment number, n 1 2 3 4 5 Moment of occurrence, ns 1.5 1.7 1.9 2.1 2.4 Fragment thickness d, nm 30 37 26 21 20

Table 1. Number, time of appearance and thickness of spalled fragments.

46 Laser Ablation - From Fundamentals to Applications

Figure 6a–d shows the state of all liquid fragments of explosive boiling at the time t = +3.0 ns. The resulting fragments eventually acquire a sufficiently high flyout velocity and under the influence of laser radiation absorption continue to evaporate until they disappear completely. The general trend in the evolution of all fragments is that over time the cavity size, flyout velocity and density of matter increase, and pressure in the cavity, thickness and temperature of the fragments decrease. So in the initially formed cavity the width increased from 5.0 nm to

Figure 6. The spatial profiles of (a) T(x), (b) ρ(x), (c) u(x), (d) p(x) in the target and the gas medium after five explosive

ments are given in Table 1.

boiling at the moment t = 3.0 ns.

Within the framework of continuum approach, a multiphase, multifront hydrodynamic model and computational algorithm were constructed that were designed for modeling ns-PLA of metal targets embedded in gaseous media. The model contains temperature dependences of optical, thermal, and thermodynamic properties of metal target (Al). The temperature dependences of heat capacity Cp(T) and thermal conductivity λ(T) of aluminum take into account their singularity in the vicinity of critical point. In optical characteristics, the metal-insulator phase transition was taken into account accompanied by the formation of the transparency wave.

The computational algorithm is based on the method of dynamic adaptation, which ideally fit to the study of problems with heterogeneous phase transitions and allows the production of numerical solutions with an explicit allocation of unlimited number of interface boundaries in condensed medium, contact boundaries, and shock waves in a gaseous medium.

The calculations have shown that under the influence of temperature dependences of optical and thermophysical properties of Al the nature of laser heating of the target in hightemperature region changes substantially. The energy release has a volume nature, and irradiated surface is markedly cooled by the process of surface evaporation. An inhomogeneity arises in the spatial temperature profile—a near-surface maximum of temperature. When the limiting superheating temperature is reached, the conditions for explosive boiling were realized with formation of a cavity in the region of the maximum, where the maximum penetration into the region of superheated metastable liquid near the critical point of liquid-vapor transition is achieved. For Al in action modes under the consideration, the typical depths of temperature maximum are d = 20–70 nm with the values ΔTmax = Tℓ,max Tsur = 150–170 K.

An approach was proposed to formulation in continuum models of homogeneous nucleation in a superheated liquid phase based on the introduction into hydrodynamic model of quasi-nuclei with thickness of 1–5 nm. Generation of quasi-nuclei is carried out by the criterion of maximum permissible overheating of initial region. The criterion of limit overheating is determined from molecular dynamics simulation.

[3] Bäuerle D. Laser Processing and Chemistry. Singapore: Springer; 2000

Wiley & Sons; 1994

John Wiley & Sons; 2006

Processing. 2007;87:47-55

ridge: Elsevier; 2016. p. 281-299

berg: Springer-Verlag; 2013

2013;113:143305

Materials Chemistry B. 2016;4:7852-7858

Applied Physics. 2007;101(2):24922-24937

mentals and applications. Chemical Reviews. 2017:1-114

2007

2007

[4] Cheng J, Chang-Sheng Liu S, Shang D, Li W, Perrie G, Dearden K. Watkins. A review of ultrafast laser materials micromachining. Optics and Laser Technology. 2013;46:88-102

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling

http://dx.doi.org/10.5772/intechopen.70773

49

[5] Chrisey DB, Hubler GK, editors. Pulsed Laser Deposition of Thin Films. New York: John

[6] Eason R, editor. Pulsed Laser Deposition of Thin Films. Hoboken: John Wiley & Sons;

[7] Singh JP, Thakur SN. Laser-Induced Breakdown Spectroscopy. Amsterdam: Elsevier;

[8] Cremers DA, Radziemski LJ. Handbook of Laser-Induced Breakdown Spectroscopy.

[9] Barcikowski S, Hahn AVK, Chichkov BN. Properties of nanoparticles generated during femtosecond laser machining in air and water. Applied Physics A: Materials Science &

[10] Nicolodelli G, Lizarelli RFZ, Costa MM, Bagnato VS. Time progression of ultrashort laser ablation in a transparent material. Laser Physics. 2011;21(8):1420-1427

[11] Shukla P, Waugh DG, Lawrence J, Vilar R. Laser surface structuring of ceramics, metals and polymers for biomedical applications: A review. In: Vilar R, editor. Laser Surface Modification of Biomaterials. Techniques and Applications. Amsterdam, Boston, Gamb-

[12] Al-Kattan A, Ryabchikov YV, Baati T, Chirvony V, Sánchez-Royo JF, Sentis M, Braguer D, Timoshenko VY, Anne Estève M, Kabashin AV. Ultrapure laser-synthesized Si nanoparticles with variable oxidation states for biomedical applications. Journal of

[13] Zhang D, Gökce B, Barcikowski S. Laser synthesis and processing of colloids: Funda-

[14] Schmidt V, Belegratis MR, editors. Laser Technology in Biomimetics: Basics and Applications. Series: Biological and Medical Physics, Biomedical Engineering. Berlin, Heidel-

[15] Fadeeva E, Schlie-Wolter S, Chichkov BN, Paasche G, Lenarz T. 5—Structuring of biomaterial surfaces with ultrashort pulsed laser radiation. In: Vilar R, editor. Laser Surface Modification of Biomaterials, Techniques and Applications. No 111. 2016. p. 145-172 [16] Mazhukin VI, Nossov VV, Smurov I. Modeling of plasma-controlled evaporation and surface condensation of Al induced de 1.06 and 0.248 μm laser radiations. Journal of

[17] Hussein AE, Diwakar PK, Harilal SS, Hassanein A. The role of laser wavelength on plasma generation and expansion of ablation plumes in air. Journal of Applied Physics.

Mathematical modeling using the developed technique allowed to obtain a sequence of five explosive boilings for the threshold fluence F = 5.2 J/cm2 and 14 for F = 6 J/cm2 . In the first case, the depth of ablation due to explosive boiling was 37 nm, in the second—70 nm. Due to the surface evaporation, they were 17 nm and 12 nm, respectively. The obtained data make it possible to conclude that phase explosion is the main mechanism of material removal in nanosecond range of laser action with fluence F ≥ 5.2 J/cm2 . The upper values of τ and F are limited by the processes of supercritical expansion and plasma formation and are subject to determination.

In metals, as in nonmetals, the transition from surface evaporation to volume removal of mass occurs in the region of temperature maximum near the critical temperature, where the influence of overheated metastable liquid phase determines the competition between surface evaporation and explosive boiling.

It should be noted that in explosive boiling regimes under investigation, cavity formation takes place in the region of positive values of pressure generated by homogeneous nucleation in superheated liquid phase, in contrast to the case of a series of spalles [57] arising in the region of negative pressures in unloading wave upon action to metal targets of ultrashort picofemtosecond laser pulses.

Validation of obtained results was performed by comparison with experimental data and with calculations with atomistic models [77, 100].

## Acknowledgements

The research was funded by the Russian Scientific Foundation (grant No 15-11-00032).

## Author details

Vladimir I. Mazhukin1,2\*

## References

[3] Bäuerle D. Laser Processing and Chemistry. Singapore: Springer; 2000

with thickness of 1–5 nm. Generation of quasi-nuclei is carried out by the criterion of maximum permissible overheating of initial region. The criterion of limit overheating is determined from

Mathematical modeling using the developed technique allowed to obtain a sequence of five

the depth of ablation due to explosive boiling was 37 nm, in the second—70 nm. Due to the surface evaporation, they were 17 nm and 12 nm, respectively. The obtained data make it possible to conclude that phase explosion is the main mechanism of material removal in nano-

the processes of supercritical expansion and plasma formation and are subject to determination. In metals, as in nonmetals, the transition from surface evaporation to volume removal of mass occurs in the region of temperature maximum near the critical temperature, where the influence of overheated metastable liquid phase determines the competition between surface evap-

It should be noted that in explosive boiling regimes under investigation, cavity formation takes place in the region of positive values of pressure generated by homogeneous nucleation in superheated liquid phase, in contrast to the case of a series of spalles [57] arising in the region of negative pressures in unloading wave upon action to metal targets of ultrashort pico-

Validation of obtained results was performed by comparison with experimental data and with

The research was funded by the Russian Scientific Foundation (grant No 15-11-00032).

1 Keldysh Institute of Applied Mathematics, Russian Academy of Sciences, Moscow, Russia

[1] Miller JC, editor. Laser Ablation. Principles and Applications. Berlin: Springer; 1994 [2] Phipps CR, editor. Laser Ablation and Its Applications. New York: Springer; 2007

. In the first case,

. The upper values of τ and F are limited by

explosive boilings for the threshold fluence F = 5.2 J/cm2 and 14 for F = 6 J/cm2

molecular dynamics simulation.

48 Laser Ablation - From Fundamentals to Applications

oration and explosive boiling.

femtosecond laser pulses.

Acknowledgements

Vladimir I. Mazhukin1,2\*

Author details

References

calculations with atomistic models [77, 100].

\*Address all correspondence to<EMAIL_ADDRESS>
2 National Research Nuclear University "MEPhI", Moscow, Russia

second range of laser action with fluence F ≥ 5.2 J/cm2

[18] Mazhukin VI, Mazhukin AV, Lobok MG. Comparison of nano- and femtosecond laser ablation of aluminium. Laser Physics. 2009;19(5):1169-1178

[33] Meyers MA, Mishra A, Benson DJ. Mechanical properties of nanocrystalline materials.

Nanosecond Laser Ablation: Mathematical Models, Computational Algorithms, Modeling