Research Papers: Heat Transfer in Manufacturing

Numerical Simulation of Transient Three-Dimensional Temperature and Kerf Formation in Laser Fusion Cutting

[+] Author and Article Information
Karim Kheloufi

Laser Material Processing Team,
Centre de Développement des
Technologies Avancées,
P.O. Box 17, Baba-Hassen,
Algiers 16303, Algeria
e-mail: kkheloufi@cdta.dz

El Hachemi Amara

Laser Material Processing Team,
Centre de Développement des
Technologies Avancées,
P.O. Box 17, Baba-Hassen,
Algiers 16303, Algeria
e-mail: amara@cdta.dz

Ahmed Benzaoui

Laboratoire Thermodynamique et
Systèmes Energétiques (LTSE),
Faculté de Physique,
Université des Sciences et de la Technologie,
Houari Boumediene (USTHB),
B.P. 32, El Alia - Bab Ezzouar,
Algiers 16111, Algeria
e-mail: abenzaoui@usthb.dz

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received June 8, 2014; final manuscript received May 7, 2015; published online June 9, 2015. Assoc. Editor: Wilson K. S. Chiu.

J. Heat Transfer 137(11), 112101 (Jun 09, 2015) (9 pages) Paper No: HT-14-1391; doi: 10.1115/1.4030658 History: Received June 08, 2014

In the present study, a three-dimensional transient numerical model was developed to study the temperature field and cutting kerf shape during laser fusion cutting. The finite volume model has been constructed, based on the Navier–Stokes equations and energy conservation equation for the description of momentum and heat transport phenomena, and the volume of fluid (VOF) method for free surface tracking. The Fresnel absorption model is used to handle the absorption of the incident wave by the surface of the liquid metal, and the enthalpy-porosity technique is employed to account for the latent heat during melting and solidification of the material. To model the physical phenomena occurring at the liquid film/gas interface, including momentum/heat transfer, a new approach is proposed which consists of treating friction force, pressure force applied by the gas jet, and the heat absorbed by the cutting front surface as source terms incorporated into the governing equations. All these physics are coupled and solved simultaneously in fluent CFD®. The main objective of using a transient phase change model in the current case is to simulate the dynamics and geometry of a growing laser-cutting generated kerf until it becomes fully developed. The model is used to investigate the effect of some process parameters on temperature fields and the formed kerf geometry.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Schuöcker, D. , 1986, “Dynamic Phenomena in Laser Cutting and Cut Quality,” J. Appl. Phys., B 40(1), pp. 9–14. [CrossRef]
Schuöcker, D. , and Muller, P. , 1987, “High Power Lasers,” Proc. SPIE, 801, pp. 258–264.
Schulz, W. , Simon, G. , Urbassek, H. M. , and Decker, I. , 1987, “On Laser Fusion Cutting of Metals,” J. Phys. D: Appl. Phys., 20, pp. 481–488. [CrossRef]
Schulz, W. , and Becker, D. , 1989, “Mathematical Simulation of Laser Treatment Materials,” European Scientific Laser Workshop on Mathematical Simulation (Lisbon), pp. 178–200.
Vicanek, M. , Simon, G. , Urbassek, H. M. , and Decker, I. , 1987, “Hydrodynamical Instability of Melt Flow in Laser Cutting,” J. Phys. D: Appl. Phys., 20(1), pp.140–145. [CrossRef]
Petring, D. , Abels, P. , and Beyer, E. , 1988,“Absorption Distribution on Idealized Cutting Front Geometries and Its Significance for Laser Beam Cutting,” Proc. SPIE., 1020, pp. 123–131.
Tani, G. , Tomesani, L. , and Campana, G. , 2003, “Prediction of Melt Geometry in Laser Cutting,” Appl. Surf. Sci., 208–209, pp.142–147. [CrossRef]
Yilbas, B. S. , 1997, “A Study Into Laser Cutting Process,” Heat Mass Transfer, 32(3), pp. 175–180. [CrossRef]
Kaplan, A. F. H. , 1996, “An Analytical Model of Metal Cutting With a Laser Beam,” J. Appl. Phys., 79(5), pp. 2198– 2208. [CrossRef]
Gross, M. S. , 2006, “On Gas Dynamic Effects in the Modeling of Laser Cutting Processes,” Appl. Math. Model., 30(4), pp. 307–318. [CrossRef]
Kovalev, O. B. , Orishich, A. M. , Fomin, V. M. , and Shulyat'ev, V. B. , 2001, “Adjoint Problems of Mechanics of Continuous Media in Gas-Laser Cutting of Metals,” J. Appl. Tech. Phys., 42(6), pp. 1014–1022. [CrossRef]
Gross, M. S. , Black, I. , and Muller, W. H. , 2003, “Computer Simulation of the Processing of Engineering Materials With Lasers-Theory and First Applications,” J. Phys. D: Appl. Phys., 36, pp. 307–318. [CrossRef]
Hu, J. , Qiu, X. , Zhang, Z. , and Shen, H. , 2012, “Influence of Oxidation on Flow Structure in Laser-Oxygen Cutting,” J. Appl. Phys., 112, p. 063107. [CrossRef]
Otto, A. , and Schmidt, M. , 2010, “Towards a Universal Numerical Simulation Model for Laser Material Processing,” Phys. Procedia 5(1), pp. 35–46. [CrossRef]
Kohl, S. , and Schmidt, M. , 2013,”Numerical Analysis of the Influence of Beam Characteristics Into the Process Dynamics During Laser Cutting,” Proceedings of ICALEO, pp. 104–109.
Chen, K. , Yao, Y. L. , and Modi, V. , 1999, “Numerical Simulation of Oxidation Effects in the Laser Cutting Process,” Int. J. Adv. Manuf. Tech., 15(11), pp. 835–842. [CrossRef]
Powell, J. , Petring, D. , and Kumar, R. V. , 2009, “Laser–Oxygen Cutting of Mild Steel: The Thermodynamics of the Oxidation Reaction,” J. Phys. D: Appl. Phys. 42, pp. 1–11. [CrossRef]
Modest, M. F. , and Abakians, H. , 1986, “Evaporative Cutting of a Semi-Infinite Body With a Moving CW Laser,” ASME J. Heat Transfer, 108(3), pp. 602–607. [CrossRef]
Chen, K. , and Yao, Y. L. , 1999, “Striation Formation and Melt Removal in the Laser Cutting Process,” J. Manuf. Process. 1.1, pp. 47–53.
Arata, Y. , Maruo, H. , Miyamoto, I. , and Takeuchi, S. , 1979, “Dynamic Behavior in Laser Gas Cutting of Mild Steel,” Trans. JWRI, 8(2), pp. 15–26.
Ivarson, A. , Powell, J. , Kamalu, J. , and Magnusson, C. , 1994, “The Oxidation Dynamics of Laser Cutting of Mild Steel and the Generation of Striations on the Cut Edge,” J. Mater. Process. Technol., 40(3–4), pp. 359–374. [CrossRef]
Chen, K. , Yao,Y. L. , and Modi, V. , 2000, “Gas Jet–Workpiece Interactions in Laser Machining,” ASME J. Manuf. Sci. Eng. 122(3), pp. 429–438. [CrossRef]
Bennon, W. D. , and Incropera, F. P. , 1987, “A Continuum Model for Momentum, Heat and Species Transport in Binary Solid-Liquid Phase Change Systems—I Model Formulation,” Int. J. Heat Mass Transfer, 30(10), pp. 2161–2170. [CrossRef]
Hirt, C. W. , and Nichols, B. D. , 1981, “Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries,” J. Comput. Phys. 39(1). pp. 201–225. [CrossRef]
Butler, K. M. , Ohlemiller, T. J. , and Linteris, G. T. , 2004, “Progress Report on Numerical Modeling of Experimental Polymer Melt Flow Behavior,” International Interflam Conference, Vol. 2, pp. 937–948.
Voller, V. R. , and Prakash, C. , 1987, “A Fixed-Grid Numerical Modeling Methodology for Convection-Diffusion Mushy Region Phase-Change Problems,” Int. J. Heat Mass Transfer, 30(8), pp. 1709–1719. [CrossRef]
Kim, C. H. , Zhang, W. , and DebRoy, T. , 2003, “Modeling of Temperature Field and Solidified Surface Profile During Gas-Metal Arc Fillet Welding,” J. Appl. Phys., 94(4), pp. 2667–2679. [CrossRef]
Kovalev, O. B. , and Orishich, A. M. , 2004, “Theory of Metal Surface Destruction Under the Action of Laser Radiation,” Doklady Phys., 49(3), pp. 175–178. [CrossRef]
Kovalev, O. B. , and Zaitsev, A. V. , 2005, “Modeling of the Free-Surface Shape in Laser Cutting of Metals. 2. Model of Multiple Reflection and Absorption of Radiation,” ASME J. Appl. Mech. Tech. Phys., 46(1), pp. 9–13. [CrossRef]
Ahmadi, B. , Torkamany, M. J. , Jaleh, B. , and Sabaghzadeh, J. , 2009, “Theoretical Comparison of Oxygen Assisted Cutting by CO2 and Yb:YAG Fiber Lasers,” Chin. J. Phys., 47(4), pp. 465–475.
Zaitsev, A. V. , Kovalev, O. B. , Orishich, A. M. , and Fomin, V. M. , 2005, “Numerical Analysis of the Effect of the TEM00 Radiation Mode Polarization on the Cut Shape in Laser Cutting of Thick Metal Sheets,” Quantum Electron., 35(2), pp. 200–204. [CrossRef]
Steen, M. , and Mazumder, J. , 2010, Laser Material Processing, 4th ed. Springer, London.
fluent 6.3, User's Guide, 2006, Fluent, Inc.
Petring, D. , 1994, “Anwendungsorientierte Modellierung desLaserstrahl-schneidens zur rechnergestützten Prozeßoptimierung,” Ph.D. thesis, RWTH Aachen, Germany, Shaker Verlag.
Beyer, E. , 1995, Schweißenmit Laser, Springer, Berlin.
Mahrle, A. , and Beyer, E. , 2009, “Theoretical Aspects of Fibre Laser Cutting,” J. Phys. D, Appl. Phys., 42(17), p. 175507. [CrossRef]


Grahic Jump Location
Fig. 1

Application of gas friction force on the surface melt

Grahic Jump Location
Fig. 2

(a) Computational domain and physical model and (b) Boundary conditions

Grahic Jump Location
Fig. 3

Flowchart for computation process

Grahic Jump Location
Fig. 4

Thermal field in the metal sheet during cutting from the edge; P = 1200 W, Ucut = 30 mm/s, and Ug = 100 m/s

Grahic Jump Location
Fig. 5

Cutting from a point in the top surface

Grahic Jump Location
Fig. 6

Thermal field in the metal during heating stage P = 1200 W, Ucut = 30 mm/s, Ug = 100 m/s (a) t = 0, (b) t = 2 × 10−3, and (c) t = 4.5 × 10−3

Grahic Jump Location
Fig. 7

Thermal field during melting and kerf formation stage P = 1200 W, Ucut = 30 mm/s, Ug = 100 m/s; (a) T = 1.8 × 10−2, (b) t = 2.1 × 10−2, and (c) t = 5.6 × 10−2 s

Grahic Jump Location
Fig. 8

Thermal field in the metal for a fully developed kerf in different views; P = 1200 W, Ucut = 40 mm/s, Ug = 100 m/s

Grahic Jump Location
Fig. 9

Thermal field in the kerf; P = 1200 W, Ug = 100 m/s; (a) Ucut = 10 mm/s and (b) Ucut = 30 mm/s

Grahic Jump Location
Fig. 10

Temperature variation along the melting front (kerf projection on the cutting direction); P = 1200 W, Ug = 100 m/s; (a) Ucut = 10 mm/s and (b) Ucut = 30 mm/s

Grahic Jump Location
Fig. 11

Liquid fraction showing the molten layer in the central plane P = 1200 W, Ug = 100 m/s; (a) Uc=30 mm/s; (b) Uc=20 m/s; (c) Uc=10m/s

Grahic Jump Location
Fig. 12

Velocity magnitude in the liquid film, P = 1200 W, Ug = 100 m/s; (a) Ucut = 10 mm/s and (b) Ucut = 30 mm/s

Grahic Jump Location
Fig. 13

Velocity magnitude variation at the interface along the melting front; (a) Ucut = 10 mm/s and (b) Ucut = 30 mm/s

Grahic Jump Location
Fig. 14

Velocity filed in liquid film

Grahic Jump Location
Fig. 15

Variation of kerf width with cutting speed

Grahic Jump Location
Fig. 16

Effect of assist gas velocity on the kerf width



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In