Abstract
Many additive manufacturing technologies involve the deposition of particles onto a surface followed by selective, targeted, laser heating. This paper develops a modular computational framework which combines the various steps within this overall process. Specifically, the framework synthesizes the following:
-
particle dynamics, which primarily entails: (a) the movement of the particles induced by contact with the surface, (b) particle-to-particle contact forces and (c) near-field interaction and external electromagnetic fields.
-
laser-input, which primarily entails: (a) absorption of laser energy input and (b) beam interference (attenuation) from particles and
-
particle thermodynamics, which primarily entails: (a) heat transfer between particles in contact by conduction and (b) subsequent thermal softening of the particles.
Numerical examples are provided and extensions are also addressed for two advanced processing scenarios involving solid-liquid-gas phase transformations.
Similar content being viewed by others
Notes
For an early history of the printed electronics field, see Gamota et al. [36]. For reviews of optical coatings and photonics, see Nakanishi et al. [65] and Maier and Atwater [58], for biosensors see Alivisatos [3], for catalysts, see Haruta [41] and for MEMS applications, see Full et al. [35] and Ho et al. [44].
There has been considerable research activity in processing of powders, in particular by compaction, for example, see Akisanya et al. [2], Anand and Gu [5], Brown and Abou-Chedid [8], Domas [20], Fleck [33], Gethin et al. [37], Gu et al. [39], Lewis et al. [55], Ransing et al. [76], Tatzel [87] and Zohdi [101, 102]. The study of “granular” or “particulate” media is wide ranging. Classical examples include the study of natural materials, such as sand and gravel, associated with coastal erosion, landslides and avalanches.
A closely related method, Electron Beam Melting, fully melts the material and produces dense solids that are void free.
\(I\!\!K_{ij}\) can be approximated by an average interfacial value of the \(i-j\) pair, \(I\!\!K_{ij}\approx \frac{I\!\!K_i+I\!\!K_j}{2}\). If the materials are the same, this collapses to simply \(I\!\!K\). As for the mechanical contact, \(A^c_{ij}\) is the contact area associated with the particle pair \((ij)\).
All system parameters can be scaled to describe any specific system of interest. They were selected simply for illustration purposes.
In the idealized limit, the temperature would be constant.
References
Ahmad Z, Rasekh M, Edirisinghe M (2010) Electrohydrodynamic direct writing of biomedical polymers and composites. Macromol Mater Eng 295:315–319
Akisanya AR, Cocks ACF, Fleck NA (1997) The yield behavior of metal powders. Int J Mech Sci 39:1315–1324
Alivisatos P (2004) The use of nanocrystals in biological detection. Nat Biotechnol 22(1):47–52
Allwood J (2005) What is sustainable manufacturing?. Lecture, Cambridge University, Cambridge
Anand L, Gu C (2000) Granular materials: constitutive equations and shear localization. J Mech Phys Solids 48:1701–1733
Avci B, Wriggers P (2012) A DEM-FEM coupling approach for the direct numerical simulation of 3D particulate flows. J Appl Mech 79:010901-1-7
Boyer TH (1988) The force on a magnetic dipole. Am J Phys 56 (8):688–692. Bibcode:1988AmJPh.56.688B. doi:10.1119/1.15501
Brown S, Abou-Chedid G (1994) Yield behavior of metal powder assemblages. J Mech Phys Solids 42:383–398
Carbonell JM, Onate E, Suarez B (2010) Modeling of ground excavation with the particle finite element method. J Eng Mech ASCE 136:455–463
Choi S, Park I, Hao Z, Holman HY, Pisano AP, Zohdi TI (2010) Ultra-fast self-assembly of micro-scale particles by open channel flow. Langmuir 26(7):4661–4667
Choi S, Stassi S, Pisano AP, Zohdi TI (2010) Coffee-ring effect-based three dimensional patterning of micro, nanoparticle assembly with a single droplet. Langmuir 26(14):11690–11698
Choi S, Jamshidi A, Seok TJ, Zohdi TI, Wu. MC, Pisano AP (2012) Fast, high-throughput creation of size-tunable micro, nanoparticle clusters via evaporative self-assembly in picoliter-scale droplets of particle suspension. Langmuir 28(6):3102–3111
Choi S, Pisano AP, Zohdi TI (2013) An analysis of evaporative self-assembly of micro particles in printed picoliter suspension droplets. J Thin Solid Films 537(30):180–189
Cullity BD, Graham CD (2008) Introduction to magnetic Materials, 2 edn. Wiley-IEEE Press, New Jersey, p 103 ISBN 0-471-47741-9
Davis SH (2001) Theory of solidification. Cambridge University Press, Cambridge
Deckard C (1986) Method and apparatus for producing parts by selective sintering U.S. Patent 4,863,538
Demko M, Choi S, Zohdi TI, Pisano AP (2012) High resolution patterning of nanoparticles by evaporative self-assembly enabled by in-situ creation and mechanical lift-off of a polymer template. Appl Phys Lett 99:253102-1–253102-3
Demkowicz (2006) Computing with hp-adaptive finite elements. I. One- and two- dimensional elliptic and Maxwell problems. Chapman & Hall/CRC, Boca Raton
Demkowicz L, Kurtz J, Pardo D, Paszynski M, Rachowicz W, Zdunek A (2007) Computing with Hp-adaptive finite, vol 2. Frontiers: three dimensional elliptic and maxwell problems with applications. CRC Press, Taylor and Francis, Boca Raton
Domas F (1997) Eigenschaft profile und Anwendungsü bersicht von EPE und EPP. Technical report of the BASF Company
Donev A, Cisse I, Sachs D, Variano EA, Stillinger F, Connelly R, Torquato S, Chaikin P (2004) Improving the density of jammed disordered packings using ellipsoids. Science 303:990–993
Donev A, Stillinger FH, Chaikin PM, Torquato S (2004b) Unusually dense crystal ellipsoid packings. Phys Rev Lett 92:255506
Donev A, Torquato S, Stillinger F (2005a) Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles-I Algorithmic details. J Comput Phys 202:737
Donev A, Torquato S, Stillinger F (2005b) Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles-II. Application to ellipses and ellipsoids. J Comput Phys 202:765
Donev A, Torquato S, Stillinger FH (2005c) Pair correlation function characteristics of nearly jammed disordered and ordered hard-sphere packings. Phys Rev E 71:011105
Dornfeld D, Wright P (2007) Technology wedges for implementing green manufacturing. Trans North Am Manuf Res Inst 35:193–200
Duran J (1997) Sands, powders and grains. An introduction to the physics of Granular Matter. Springer, Berlin
Elmore WC, Heald MA (1985) Physics of waves. Dover Publishers, New York
Farhat C, Lesoinne M, Maman N (1995) Mixed explicit/implicit time integration of coupled aeroelastic problems: three-field formulation, geometric conservation and distributed solution. Int J Numer Methods Fluids 21:807–835
Farhat C, Lesoinne M (2000) Two efficient staggered procedures for the serial and parallel solution of three-dimensional nonlinear transient aeroelastic problems. Comput Methods Appl Mech Eng 182:499–516
Farhat C, van der Zee G, Geuzaine P (2006) Provably second-order time-accurate loosely-coupled solution algorithms for transient nonlinear computational aeroelasticity. Comput Methods Appl Mech Eng 195:1973–2001
Feynman RP, Leighton RB, Sands M (2006) The Feynman lectures on physics vol 2. Addison-Wesley, Reading ISBN 0-8053-9045-6
Fleck NA (1995) On the cold compaction of powders. J Mech Phys Solids 43:1409–1431
Frenklach M, Carmer CS (1999) Molecular dynamics using combined quantum & empirical forces: application to surface reactions. Adv Class Trajectory Methods 4:27–63
Fuller SB, Wilhelm EJ, Jacobson JM (2002) Ink-jet printed nanoparticle microelectromechanical systems. J Microelectromech Syst 11:54–60
Gamota D, Brazis P, Kalyanasundaram K, Zhang J (2004) Printed organic and molecular electronics. Kluwer Academic Publishers, New York
Gethin DT, Lewis RW, Ransing RS (2003) A discrete deformable element approach for the compaction of powder systems. Modell Simul Mater Sci Eng 11(1):101–114
Grigoropoulos CP (2009) Transport in laser microfabrication. Cambridge University Press, Cambridge
Gu C, Kim M, Anand L (2001) Constitutive equations for metal powders: application to powder forming processes. Int J Plast 17:147–209
Haile JM (1992) Molecular dynamics simulations: elementary methods. Wiley, New York
Haruta M (2002) Catalysis of gold nanoparticles deposited on metal oxides. Cattech 6(3):102–115
Hase WL (1999) Molecular dynamics of clusters, surfaces, liquids, & interfaces. Advances in classical trajectory methods, vol 4. JAI Press, Stamford
Housholder R (1979) Molding process. U.S. Patent 4,247,508
Ho C, Steingart D, Salminent J, Sin W, Rantala T, Evans J, Wright P (2006) Dispenser printed electrochemical capacitors for power management of millimeter scale lithium ion polymer microbatteries for wireless sensors. 6th International workshop on micro and nanotechnology for power generation and energy conversion applications (PowerMEMS 2006), Berkeley, CA
Huang D, Liao F, Molesa S, Redinger D, Subramanian V (2003) Plastic-compatible low-resistance printable gold nanoparticle conductors for flexible electronics. J Electrochem Soc 150(7):G412–417
Hull C (1984) Apparatus for production of three-dimensional objects by stereolithography, U.S. Patent 4,575,330
Jackson JD (1998) Classical electrodynamics, 3rd edn. Wiley, New York
Johnson K (1985) Contact mechanics. Cambridge University Press, Cambridge
Kansaal A, Torquato S, Stillinger F (2002) Diversity of order and densities in jammed hard-particle packings. Phys Rev E 66:041109
Ko SH, Park I, Pan H, Grigoropoulos CP, Pisano AP, Luscombe CK, Frechet JMJ (2007) Direct nanoimprinting of metal nanoparticles for nanoscale electronics fabrication. Nano Lett 7:1869–1877
Ko SH, Park I, Pan H, Misra N, Rogers MS, Grigoropoulos CP, Pisano AP (2008) ZnO nanowire network transistor fabrication by low temperature, all inorganic nanoparticle solution process. Appl Phys Lett 92:154102
Labra C, Onate E (2009) High-density sphere packing for discrete element method simulations. Commun Numer Method Eng 25(7):837–849
Lesoinne M, Farhat C (1998) Free staggered algorithm for nonlinear transient aeroelastic problems. AIAA J 36(9):1754–1756
Le Tallec P, Mouro J (2001) Fluid structure interaction with large structural displacements. Comput Methods Appl Mech Eng 190(24–25):3039–3067
Lewis RW, Gethin DT, Yang XSS, Rowe RC (2005) A combined finite-discrete element method for simulating pharmaceutical powder tableting. Int J Numer Methods Eng 62:853–869
Lewis RW, Schrefler BA (1998) The finite element method in the static & dynamic deformation & consolidation of porous media, 2nd edn. Wiley Press, New York
Lewis RW, Schrefler BA, Simoni L (1992) Coupling versus uncoupling in soil consolidation. Int J Numer Anal Methods Geomech 15:533–548
Maier SA, Atwater HA (2005) Plasmonics: localization and guiding of electromagnetic energy in metal/dielectric structures. J Appl Phys 2005(98):011101
Macleod HA (2001) Thin flim optcal filters, 3rd edn. I. O. P, Bristol
Martin P (2009) Handbook of deposition technologies for films and coatings, 3rd edn. Elsevier, Boston
Michopoulos G, Farhat C, Fish J (2005) Survey on modeling and simulation of multiphysics systems. J Comput Inf Sci Eng 5(3):198–213
Martin P (2011) Introduction to surface engineering and functionally engineered materials. Scrivener and Elsevier, Amsterdam
Messier R, Giri AP, Roy R (1984) Revised structure zone for thin film physical structure. J Vac Sci Technol A2(2):500
Moelwyn-Hughes EA (1961) Physical chemistry. Pergamon Press, Oxford
Nakanishi H, Bishop KJM, Kowalczyk B, Nitzan A, Weiss EA, Tretiakov KV, Apodaca MM, Klajn R, Stoddart JF, Grzybowski BA (2009) Photoconductance and inverse photoconductance in thin films of functionalized metal nanoparticles. Nature 460:371–375
Niemz MH (2004) Laser tissue interactions-fundamentals and applications. Springer, Berlin
Onate E, Idelsohn SR, Celigueta MA, Rossi R (2008) Advances in the particle finite element method for the analysis of fluid-multibody interaction and bed erosion in free surface flows. Comput Methods Appl Mech Eng 197(19–20):1777–1800
Onate E, Celigueta MA, Idelsohn SR, Salazar F, Surez B (2011) Possibilities of the particle finite element method for fluid-soil-structure interaction problems. Comput Mech 48:307–318
Park J-U, Hardy M, Kang SJ, Barton K, Adair K, Mukhopadhyay DK, Lee CY, Strano MS, Alleyne AG, Georgiadis JG, Ferreira PM, Rogers JA (2007) High-resolution electrohydrodynamic jet printing. Nat Mater 6:782–789
Park I, Ko SH, Pan H, Grigoropoulos CP, Pisano AP, Frechet JMJ, Lee ES, Jeong JH (2008) Nanoscale patterning and electronics on flexible substrate by direct nanoimprinting of metallic nanoparticles. Adv Mater 20:489
Park KC, Felippa CA (1983) Partitioned analysis of coupled systems. In: Belytschko T, Hughes TJR (eds) Computational methods for transient analysis. Elsevier Science Publishers, North-Holland
Piperno S (1997) Explicit/implicit fluid/structure staggered procedures with a structural predictor & fluid subcycling for 2D inviscid aeroelastic simulations. Int J Numer Method Fluids 25:1207–1226
Piperno S, Farhat C, Larrouturou B (1995) Partitioned procedures for the transient solution of coupled aeroelastic problems—Part I: model problem, theory, and two-dimensional application. Comput Method Appl Mech Eng 124(1–2):79–112
Piperno S, Farhat C (2001) Partitioned procedures for the transient solution of coupled aeroelastic problems—Part II: energy transfer analysis and three-dimensional applications. Comput Methods Appl Mech Eng 190:3147–3170
Pö schel T, Schwager T (2004) Computational granular dynamics. Springer, Berlin
Ransing RS, Lewis RW, Gethin DT (2004) Using a deformable discrete-element technique to model the compaction behaviour of mixed ductile and brittle particulate systems. Philos Trans R Soc 362(1822):1867–1884
Rapaport DC (1995) The art of molecular dynamics simulation. Cambridge University Press, Cambridge
Reich-Weiser C, Vijayaraghavan A, Dornfeld DA (2008) Metrics for manufacturing sustainability, Proceedings of IMSEC, ASME, Evanston, IL, Oct 7–10
Ridley BA, Nivi B, Jacobson JM (1999) All-inorganic field effect transistors fabricated by printing. Science 286:746–749
Rojek J, Labra C, Su O, Onate E (2012) Comparative study of different discrete element models and evaluation of equivalent micromechanical parameters. Int J Solids Struct 49:1497–1517. doi:10.1016/j.ijsolstr.2012.02.032
Rosen M, Dincer I, Kanoglu M (2008) Role of exergy in increasing efficiency and sustainability and reducing environmental impact. Energy Policy 36:128–137
Samarasinghe SR, Pastoriza-Santos I, Edirisinghe MJ, Reece MJ, Liz-Marzan LM (2006) Printing gold nanoparticles with an electrohydrodynamic direct write device. Gold Bull 39:48–53
Schlick T (2000) Molecular modeling & simulation. An interdisciplinary guide. Springer, New York
Sirringhaus H, Kawase T, Friend RH, Shimoda T, Inbasekaran M, Wu W, Woo EP (2000) High-resolution inkjet printing of all-polymer transistor circuits. Science 290:2123–2126
Steen WM (1998) Laser material processing. Springer, Berlin
Stillinger FH, Weber TA (1985) Computer simulation of local order in condensed phases of silicon. Phys Rev B 31:5262–5271
Tatzel H (1996) Grundlagen der Verarbeitungstechnik von EPP-Bewä hrte und neue Verfahren. Technical report of the BASF Company
Tersoff J (1988) Empirical interatomic potential for carbon, with applications to amorphous carbon. Phys Rev Lett 61:2879–2882
Torquato S (2002) Random heterogeneous materials: microstructure & macroscopic properties. Springer, New York
Wang JZ, Zheng ZH, Li HW, Huck WTS, Sirringhaus H (2004) Dewetting of conducting polymer inkjet droplets on patterned surfaces. Nat Mater 3:171–176
Wellmann C, Lillie C, Wriggers P (2008) A contact detection algorithm for superellipsoids based on the common-normal concept. Eng Comput 25:432–442
Wellmann C, Lillie C, Wriggers P (2008) Homogenization of granular material modelled by a three-dimensional discrete element method. Comput Geotech 35:394–405
Wellmann C, Lillie C, Wriggers P (2008) Comparison of the macroscopic behavior of granular materials modeled by different constitutive equations on the microscale. Finite Elem Anal Des 44:259–271
Wellmann C, Wriggers P (2011) Homogenization of granular material modeled by a 3D DEM, particle-based methods: fundamentals and applications. In: Onate E, Owen DRJ (eds) pp 211–231
Wellmann C, Wriggers P (2012) A two-scale model of granular materials. Comput Method Appl Mech Eng 205–208:46–58
Widom B (1966) Random sequential addition of hard spheres to a volume. J Chem Phys 44:3888–3894
Wriggers P (2002) Computational contact mechanics. John-Wiley, Chichester
Wriggers P (2008) Nonlinear finite element analysis. Springer, Berlin
Zienkiewicz OC (1984) Coupled problems & their numerical solution. In: Lewis RW, Bettes P, Hinton E (eds) Numerical methods in coupled systems. Wiley, Chichester, pp 35–58
Zienkiewicz OC, Paul DK, Chan AHC (1988) Unconditionally stable staggered solution procedure for soil-pore fluid interaction problems. Int J Numer Methods Eng 26:1039–1055
Zohdi TI (2003) On the compaction of cohesive hyperelastic granules at finite strains. Proc Roy Soc 454(2034):1395–1401
Zohdi TI (2003b) Genetic optimization of statistically uncertain microheterogeneous solids. Philos Trans R Soc Math Phys Eng Sci 361(1806):1021–1043
Zohdi TI (2006) Computation of the coupled thermo-optical scattering properties of random particulate systems. Comput Methods Appl Mech Eng 195:5813–5830
Zohdi TI, Kuypers FA (2006) Modeling and rapid simulation of multiple red blood cell light scattering. Proc R Soc Interface 3(11):823–831
Zohdi TI (2012) Modeling and simulation of the optical response rod-functionalized reflective surfaces. Comput Mech 50(2):257–268
Zohdi TI (2013) Rapid simulation of laser processing of discrete particulate materials. Arch Comput Methods Eng. doi:10.1007/s11831-013-9092-6. pp 1–17
Zohdi TI (2002) An adaptive-recursive staggering strategy for simulating multifield coupled processes in microheterogeneous solids. Int J Numer Methods Eng 53:1511–1532
Zohdi TI (2005) Charge-induced clustering in multifield particulate flow. Int J Numer Methods Eng 62(7):870–898
Zohdi TI (2007) Computation of strongly coupled multifield interaction in particle-fluid systems. Comput Methods Appl Mech Eng 196:3927–3950
Zohdi TI, Wriggers P (2008) Introduction to computational micromechanics. Second Reprinting. Springer, Berlin
Zohdi TI (2010) On the dynamics of charged electromagnetic particulate jets. Arch Comput Methods Eng 17(2):109–135
Zohdi TI (2010) Simulation of coupled microscale multiphysical-fields in particulate-doped dielectrics with staggered adaptive FDTD. Comput Methods Appl Mech Eng 199:79–101
Zohdi TI (2012) Electromagnetic properties of multiphase dielectrics. A primer on modeling, theory and computation. Springer, Berlin
Zohdi TI (2012) Dynamics of charged particulate systems. Modeling, theory and computation. Springer, Berlin
Zohdi TI (2013) Numerical simulation of charged particulate cluster-droplet impact on electrified surfaces. J Comput Phys 233:509–526
Zohdi TI (2014) A direct particle-based computational framework for electrically-enhanced thermo-mechanical sintering of powdered materials. Math Mech Solids. doi:10.1007/s11831-013-9092-6. pp 1–21
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendix 1: Contact area parameter and alternative models
1.1 Contact area parameter
Following Zohdi [116], and referring to Fig. 12, one can solve for an approximation of the common contact radius \(a_{ij}\) (and the contact area, \(A^c_{ij}=\pi a_{ij}^2\)) by solving the following three equations,
and
and
where \(R_i\) is the radius of particle \(i\), \(R_j\) is the radius of particle \(j\), \(L_i\) is the distance from the center of particle \(i\) and the common contact interpenetration line and \(L_j\) is the distance from the center of particle \(j\) and the common contact interpenetration line, where the extent of interpenetration is
The above equations yield an expression \(a_{ij}\), which yields an expression for the contact area parameter
where
1.2 Aternative models
One could easily construct more elaborate relations connecting the relative proximity of the particles and other metrics to the contact force, \(\varvec{\Psi }^{con,n}_{ij}\propto \mathcal{F}({\varvec{r}}_i, {\varvec{r}}_j, {\varvec{n}}_{ij},R_i,R_j, \ldots )\), building on, for example, Hertzian contact models. This poses no difficulty in the direct numerical method developed. For the remainder of the analysis, we shall use the deformation metric in Eq. 2.8. For detailed treatments, see Wellman et al. [91–95] and Avci and Wriggers [6]. We note that with the appropriate definition of parameters, one can recover Hertz, Bradley, Johnson–Kendel–Roberts, Derjaguin–Muller–Toporov contact models. For example, Hertzian contact is widely used, with the assumptions being
-
frictionless, continuous, surfaces,
-
each of contacting bodies are elastic half-spaces, whereby the contact area dimensions are smaller radii of the bodies and,
-
the bodies remain elastic (infinitesimal strains),
results in the following contact force:
which has the general form of \(\varvec{\Psi }^{con,n}=K^*_{ij}\delta ^p_{ij}\), where
-
\(R^*=\left( \frac{1}{R_i}+\frac{1}{R_j}\right) ^{-1}\) and
-
\(E^*=\left( \frac{1-\nu _i^2}{E_i}+\frac{1-\nu _j^2}{E_j}\right) ^{-1}\),
where \(E\) is the Young’s modulus and \(\nu \) is the Poisson ratio, The contact area with such a model has already been incorporated in the relation above, and is equal to \(A^c_{ij}=\pi a^2\) where \(a=\sqrt{R^*\delta _{ij}}\). For more details, we refer the reader to Johnson [48]. Furthermore, we remark that the normal contact between a particle and a wall, with a Hertzian model is given by
where \(R_j=R_w=\infty \) (see Eq. 8.7)
-
\(R^*=R_i\) and
-
\(E^*=\left( \frac{1-\nu _i^2}{E_i}+\frac{1-\nu _j^2}{E_j}\right) ^{-1}\).
It is obvious that for a deeper understanding of the deformation within a particle, it must be treated as a deformable continuum, which would require a highly-resolved spatial discretization, for example using the finite element method for the contacting bodies. This requires a large computational effort. For the state of the art in finite element methods and Contact Mechanics, see the books of Wriggers [97, 98]. For work specifically focusing on the continuum mechanics of particles, see Zohdi and Wriggers [111].
Appendix 2: Phase transformations
To include phase transformations, we consider seven cases, which are implemented in a predictor-corrector manner by first solving
to obtain predicted temperature, and then checking the following:
-
Solid \(\rightarrow \) solid-no melting with \(C_i=C_S\): If \(\theta (t)<\theta _m\) and \(\theta (t+\Delta t)<\theta _m\) then retain Eq. 9.1 with \(C(\theta )=C_S\),
-
Solid \(\rightarrow \) liquid-melting with \(C_i=C_S\): If \(\theta (t)<\theta _m\) and \(\theta (t+\Delta t)\ge \theta _m\) then re-solve Eq. 9.1 with \(C(\theta )=C_S+\frac{\delta \mathcal{P}^{S\rightarrow L}}{\delta \theta }\),
-
Liquid \(\rightarrow \) liquid-melted with \(C_i=C_L\): If \(\theta (t)\ge \theta _m\) and \(\theta (t+\Delta t)\ge \theta _m\) then retain Eq. 9.1 with \(C(\theta )=C_L\),
-
Liquid \(\rightarrow \) solid-solidification with \(C_i=C_L\): If \(\theta (t)\ge \theta _m\) and \(\theta (t+\Delta t)<\theta _m\) then re-solve Eq. 9.1 with \(C(\theta )=C_L+\frac{\delta \mathcal{P}^{L\rightarrow S}}{\delta \theta }\),
-
Liquid \(\rightarrow \) vapor-vaporization with \(C_i=C_L\): If \(\theta (t)<\theta _v\) and \(\theta (t+\Delta t)\ge \theta _v\) then re-solve Eq. 9.1 with \(C(\theta )=C_L+\frac{\delta \mathcal{P}^{L\rightarrow V}}{\delta \theta }\),
-
Vapor \(\rightarrow \) vapor-remains a vapor with \(C_i=C_V\): If \(\theta (t)\ge \theta _v\) and \(\theta (t+\Delta t)\ge \theta _v\) then retain Eq. 9.1 with \(C(\theta )=C_V\),
-
Vapor \(\rightarrow \) liquid-condensation with \(C_i=C_v\): If \(\theta (t)\ge \theta _v\) and \(\theta (t+\Delta t)<\theta _v\) then re-solve Eq. 9.1 with \(C(\theta )=C_V+\frac{\delta \mathcal{P}^{V\rightarrow L}}{\delta \theta }\),
where \(C_S\) is the heat capacity of the solid, \(C_L\) is the heat capacity of the liquid and \(C_V\) is the heat capacity of the vapor and
-
\(0<\delta \mathcal{P}^{S\rightarrow L}\) is the latent heat of melting,
-
\(0<\delta \mathcal{P}^{L\rightarrow S}\) is the latent heat of solidification,
-
\(0<\delta \mathcal{P}^{L\rightarrow V}\) is the latent heat of vaporization,
-
\(0<\delta \mathcal{P}^{V\rightarrow L}\) is the latent heat of condensation and
-
\(0<\delta \theta \) is small and can be thought of as a “bandwidth” for a phase transformation. For more details on melting processes, see Davis [15].
We note that latent heats have a tendency to resist the phase transformations, achieved by adding the positive terms in the denominator, thus enforcing reduced temperature (during the phase transformation).Footnote 6 This approach is relatively straightforward to include within the staggering framework. As a consequence, the number of particles in the system and their heat capacities will also change in the algorithm, and Eq. 4.10 becomes
The subsequent convergence of the thermal calculation is rather quick, since the time steps are extremely small. For more details on convergence on iterative time-stepping schemes, see Zohdi [107–116].
Rights and permissions
About this article
Cite this article
Zohdi, T.I. Additive particle deposition and selective laser processing-a computational manufacturing framework. Comput Mech 54, 171–191 (2014). https://doi.org/10.1007/s00466-014-1012-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00466-014-1012-6