arxiv: v1 [cs.ce] 11 Feb 2014

Transcription

1 International Journal of Modern Physics C c World Scientific Publishing Copany arxiv: v1 [cs.ce] 11 Feb 2014 VALIDATION EXPERIMENTS FOR LBM SIMULATIONS OF ELECTRON BEAM MELTING REGINA AMMER and ULRICH RÜDE Departent of Coputer Science 10 Syste Siulation, University of Erlangen-Nureberg, Cauerstr. 11, Erlangen, Gerany regina.aer@fau.de ulrich.ruede@fau.de MATTHIAS MARKL, VERA JÜCHTER and CAROLIN KÖRNER Chair of Metals Science and Technology, University of Erlangen-Nureberg, Martensstr. 5, Erlangen, Gerany atthias.arkl@fau.de vera.juechter@ww.uni-erlangen.de carolin.koerner@ww.uni-erlangen.de Received Day Month Year Revised Day Month Year This paper validates 3D siulation results of electron bea elting (EBM) processes coparing experiental and nuerical data. The physical setup is presented which is discretized by a three diensional (3D) theral lattice Boltzann ethod (LBM). An experiental process window is used for the validation depending on the line energy injected into the etal powder bed and the scan velocity of the electron bea. In the process window the EBM products are classified into the categories, porous, good and swelling, depending on the quality of the surface. The sae paraeter sets are used to generate a nuerical process window. A coparison of nuerical and experiental process windows shows a good agreeent. This validates the EBM odel and justifies siulations for future iproveents of EBM processes. In particular nuerical siulations can be used to explain future process window scenarios and find the best paraeter set for a good surface quality and dense products. Keywords: 3D theral lattice Boltzann ethod; free surface; electron bea elting; validation experients. PACS Nos.: Hf, 123.1K 1. Introduction Electron bea elting (EBM) is an additive anufacturing ethod used to produce coplex etallic structures layer by layer fro etal powder 1. EBM opens new opportunities in any industries, ranging fro aircraft anufacturers to producers of edical iplants. However, until now the parts cannot be anufactured at sufficient speed to ake the econoically viable for any but specific very high 1

2 2 Regina Aer, Matthias Markl, Vera Jüchter, Carolin Körner, Ulrich Rüde value applications. In order to accelerate the building process and iprove the accuracy a better understanding of the bea-powder interaction is necessary. This can be gained by 3D siulations of the process. Basedon the workof Körneret al. 2 we use a 3Dtheral LBM for the discretization for the EBM processes. This odel includes hydrodynaic effects, like elt flow, capillarity and wetting, as well as theral effects, like absorption of the bea energy, elting and solidification. Detailed inforation for parallel, optiized 3D absorption algoriths is found in 3 and the description of the odeling aspects, e.g., the generation of the Ti-Al6-V4 powder particles is given in 4. The ipleentation of the odel is ebedded in the walberla fraework (widely applicable lattice Boltzan solver fro Erlangen) which is a lattice Boltzann based fluid flow solver with highly parallelized kernels 5,6. In 4 basic validation exaples are already shown, like the test of the Stefan proble as benchark for the siulation of solid-liquid phase transition, and odest EBM process exaples, e.g., the elting of a spot regarding the interaction of powder particles and electron bea. In this paper we extend the validation exaples focusing on the EBM processes, e.g, hatching scenarios, and copare the nuerical results with experiental data. These nuerical results deonstrate the good potential of the theral LBM approach to understand and predict coplex processes like the EBM depending on therodydnaic as well as on fluid dynaic phenoena. The reainder of the paper is organized as follows: Section 2 describes the 3D theral LBM used for the siulation of the EBM process. The next section defines line energy, scan velocity, porosity and swelling before showing the experiental process window. Subsequently, the generation of the nuerical process window is described, shown and copared with the experiental one. Section 4 concludes the validation experients and outlines future research topics. 2. Nuerical ethods Based on the lattice Boltzann equation the isotheral lattice Boltzann ethod (LBM) is introduced by 7 as an iproveent of its predecessor, the lattice gas autoata (LGA) 8. LBM is an alternative, esoscopic approach for the nuerical solution of the Navier-Stokes equations 9 and can be parallalized easier than other traditional CFD ethods whose coputation is based on the conservation of acroscopic quantities. LBM has also less coputational storage and work requireents than olecular dynaic ethods which consider icroscopic physics. Theral LB ethods are classified into three categories: the ultispeed LBM 10,11, the double-distribution or ultidistribution approach 12,13,14 and hybrid approach where the scalar teperature equation is solved by finite differences or finite volue ethods 15. The ultispeed LBM can be seen as an extension of the isotheral LBM where the equilibriu distribution function depends also on the teperature and has higher nonlinear velocity ters. This ethod suffers fro nuericalinstabilities which haveto be stabilized asit isshown in 16. An additional

3 Validation Experients for LBM Siulations of Electron Bea Melting 3 disadvantage is that only one Prandtl nuber can be siulated by ultispeed LBM which liits the range of applications. The ultidistribution LBM overcoes these drawbacks using a separate distribution function for the teperature field. For the siulation of the EBM process a theral two-distribution LBM is used which is described followed by the treatent of the free surface boundary condition D theral lattice Boltzann ethod The odel for the EBM process is based on the work of Körner et al. in 2 where a 2D theral LBM is discussed. The idea of LBM is solving the Boltzann equation in the hydrodynaic liit for a particle distribution function (pdf) in the physical oentu space. This pdf f(x,v,t) indicates the probability of finding a particle with the acroscopicvelocity v(x,t) at position x and tie t. For the theral LBM a second pdf set h(x,v,t) is used for the nuerical solution ofthe teperature field. The collision operator is approxiated by the BGK-collision operator 17,11, Ω f = 1 τ f (f f eq ), Ω h = 1 τ h (h h eq ), (1) where f eq and h eq denote the Maxwell equilibriu distribution 19. τ f is the relaxation tie related to the viscosity ν of the fluid and given by ν = c 2 s t(τ f 0.5) with the lattice speed of sound c 2 s = x2 3 t. τ 2 h belongs to the second pdf set h i and there exists the relation to the theral diffusivity k = c 2 s t(τ h 0.5). The D3Q19 stencil is used to discretize the icroscopic space, i.e., a finite set of 19discretevelocitiese i andlatticeweightsω i aregiven 19.The Maxwellequilibriu distribution is discretized by a Taylor series expansion 20 and is given by, [ f eq i (x,t) = ω i ρ 1+ (e i u) c 2 + (e i u) 2 ] s 2c 4 u2 s 2c 2, h eq i (x,t) = ω ie s [ 1+ (e i u) c 2 s (2) where u denotes the velocity of the liquid. It is sufficient to use a linearized equilibriu distribution function h eq i for the teperature field. The acroscopic quantities density ρ, oentu ρu and energy density E are coputed in the following way, ρ = i f i, ρu = i e i f i, E = i ] h i. (3) Both pdf sets areone-waycoupled by the velocityof the fluid u which is used for the coputation ofh eq i. The discretizationofthe EBMprocessby a 3Dcoupled-theral LBGK 4 can be written in a strea-collide algorith, streaing { f i (x+e i t,t+ t) = f i (x,t) h i (x+e i t,t+ t) = h i (x,t), (4) collide { fi (x,t+ t) = f i (x,t+ t)+ 1 τ f (f eq i (ρ,u) f i (x,t+ t))+f i h i (x,t+ t) = h i (x,t+ t)+ 1 τ h (h eq i (E,u) h i(x,t+ t))+φ i, (5)

4 4 Regina Aer, Matthias Markl, Vera Jüchter, Carolin Körner, Ulrich Rüde where F i denotes an external force ter, e.g., gravity, and Φ i is the energy injected into the syste by the electron bea, [ (ei u) F i = ω i ρ c 2 + (e ] i u)e i s c 4 g, Φ i (x,t) = ω i E b (x,t), (6) s where E b denotes the energy of the electron bea. Detailed inforation how the bea energy and absorption are coputed is given in 3,4. Next, the free surface boundary condition is outlined Free surface lattice Boltzann ethod Körner et al. introduce in 21 a volue-of-fluid based free surface lattice Boltzann approach. This ethod neglects the gas phase and works copletely on the liquid phase assuing that the liquid phase covers the behavior of the flow copletely. In order to ensure ass and oentu conservation a boundary condition is iposed at the interface which separates liquid and gas phase. More inforation regarding the free surface boundary condition used for this EBM odel can be found in Validation experients The following validation experients do not only cover particular aspects as the solid-liquid phase transition or only concentrate on qualitative accordance but consider also the coplete EBM process in a quantitative way Definitions and experiental setup We use the following experiental setup for the validation experients. A cuboid of size (15x15x10) 3 is generatedby hatching the Ti-Al6-V4 powder particles layer by layer (copare the red arrows in Fig. 1). The hatching differs by line energy and scan velocity of the electron bea. The line energy is defined by, E L = U I v scan = P bea v scan, (7) where U denotes the acceleration voltage in [V], Fig. 1. Experiental setup. I the current in [A] and v scan the scan velocity in [ s ] of the electron bea. The paraeter set (E L,v scan ) defines the electron bea. The quality of EBM products is classified into three categories, naely porous, good and swelling. Good saples have a sooth surface and a relative density higher than 99.5%. If the teperature during the process is too high swelling can occur and the diensional accuracy cannot be guaranteed. On the other side, if the line energy is too low and the relative density is saller than 99.5%, then the saple is porous (cf. Fig. 2). Next, a experiental

5 Validation Experients for LBM Siulations of Electron Bea Melting 5 porous good swelling Increasing line energy Fig. 2. Categories of saples. process window is shown in Fig. 3 which is copared with nuerical siulation results Experiental data Fig. 3 shows the quality of saples with different line energies in [ kj ] and scan velocities in [ s ] of the electron bea. The red circles stand for a porous paraeter Line energy [kj/] good porous swelling Scan velocity [/s] Fig. 3. Process window of experiental data. Built teperature 1000 K, line offset 100 µ, layer thickness 50 µ, focused electron bea. set, the blue squares for an optial, good paraeter set and the yellow rhobus for a saple where swelling occurs because of too high teperatures at the surface. It can be observed that the higher the scan velocity the lower the ebedded line energy has to be to get an optial saple result Nuerical results We validate our 3D nuerical EBM odel against these experiental data of Fig. 3. Because of the high coputational cost of the 3D siulations we only odel the hatching of one powder layer instead of the ultiple layers.

6 6 Regina Aer, Matthias Markl, Vera Jüchter, Carolin Körner, Ulrich Rüde bea offset siulated powder bed hatching lines Fig. 4. Sketch of siulation scenario. One siplified scenario for the exeplary paraeter set of (3.2 s,0.2 J ) does already require 8 copute nodes of LiMa a for 4 wall clock hours. LiMa is the copute cluster where all 3D siulations are done. However, this siplification is in good agreeent with the situation expected for the full setupasitwillbeshowninthissection. We also iniize the siulated powder particle layer, i.e., we focus only a rectangular powder layer and have seven hatching lines (cf. Fig. 4). Thus, a bea offset per layer is defined where the electron bea is on but outside the siulated powder particle bed. In Fig. 5 the quality of the nuerical experients with different line energy and scan velocity is shown. In order to deterine if the nuerical saple is porous the relative density of the saple is easured and if it is less than 99.5% the saple is porous. The definition of nuerical swelling is ore difficult. Our odel does not include evaporation and thus no evaporation pressure. Hence, nuerical swelling effects cannot occur. We have to find another way how we can deterine regions nuerically where swelling exists experientally. Therefore, we use the nuerical coputed teperature which is higher than in real EB elting processes because the cooling effect during the process of evaporation is issing. Furtherore, nuerical artefacts cause outliers in the teperature field. In order to overcoe these nuerical influences we use an averaging process for the teperature values. If these averaged values are still higher than 7500K we assue that swellings occur in the experiental setup. In coparison with the experiental data in Fig. 3 it is observed that all porous saples are consistent in the siulation results (cf. Fig. 5). Furtherore, the nuerical and experiental results for scan velocities of 3.2 s, 4.0 s and 6.4 s agree with each other. Sall differences between nuerical and experiental data are only found for 0.8 s, 1.6 s, 2.4 s, 4.8 s and 5.6 s for higher line energies; for v scan = 0.8 kj s in experiental data swelling already occur for 0.4 and 0.5 but in nuerical data no swelling is easured. The sae effect is seen for scan velocities 1.6 and 2.4 s where swelling occurs for saller line energies experientally than nuerically.for higher scanvelocities of4.8 and 5.6 s the reverseeffect is observed, i.e., oreswelling occurnuerically than experientally. Forv scan = 4.8 s swelling arises nuerically already for a line energy of 0.25 kj while the experiental setup is still good. The sae is seen for v scan = 5.6 kj s ; line energies higher than 0.25 lead to swelling in nuerical siulations, but in experiental data only energies equal 0.3 kj arise swelling. a

7 Validation Experients for LBM Siulations of Electron Bea Melting 7 Line energy [kj/] swelling porous good Scan velocity [/s] Fig. 5. Process window of nuerical data. x = , t = s, siulated powder doain (1.44x0.64x0.24) , bea offset = Both effects can be explained by the focus of the electron bea gun. The focus is ore precise for saller scan velocities iposed by a saller electron bea power P bea, e.g, v scan = 1.6 s or v scan = 2.4 s lead to a precise, sall focus and thus, energy is brought onto a saller area. For these velocities swelling occurs for saller line energies experientally than nuerically because the EB gun focus is constant for different bea power and scan velocities, respectively. For higher scan velocities (larger P bea ) like 4.8 s the focus of the gun spreads and, subsequently, the sae aountofenergyis broughtinto alargerareaofpowderparticlesand the axiu teperature is saller. As a consequence in experiental data less swelling occurs for higher scan velocities (larger P bea ) than in nuerical data. (a) t = s (b) t = s (c) t = s Fig. 6. Hatching of one layer, (0.15 kj, 6.4 s ). Fig. 6 shows hatching of one line with paraeter set (0.15 kj, 6.4 s ) for three different tiesteps. TheinverseGaussiandistributed powderparticles 4 aredeonstrated the and the free surface is visualized by an isosurface. In Fig. 6(a) the electron bea elted one line and powder particles in this first line are not yet

8 8 Regina Aer, Matthias Markl, Vera Jüchter, Carolin Körner, Ulrich Rüde copletely elted. The next Fig. 6(b) shows the second elted line by a coingback electron bea. Here, the particles of the first line are now elted copletely. In the third Fig. 6(c) the electron bea has elted seven hatching lines and the whole powder layer is elted. 4. Conclusion This paper describes the 3D theral LB ethod for the nuerical siulation of EBM processes. Experiental data of hatching a cuboid with different line energies and scan velocities are shown in a process window. The different saples are classified into porous, good and swelling, depending on the paraeter set (line energy, scan velocity). For the validation all paraeter sets are tested nuerically and a nuerical process window is generated by 3D siulations of hatching one powder layer. Nuerical and experiental process window are highly concordant especially in the region of sall line energies because paraeter sets which are porous in experiental process window are that also in the nuerical one. Only in the region between good and swelling differences exist which are reducible to the focus of the EB gun. The precision of the focus depends on the bea power but in our siulations the focus is constant. In order to accelerate the building process of EBM anufacturing it is iportant to extend the process window, i.e., to raise up the question how the process window will look like for a higher bea power and therefore higher scan velocities. After the validation experients in this paper which justify the correct behavior of our fraework and its beneficing to siulate EBM processes we will try to find better paraeter sets regarding quality of EBM products. Acknowledgents Our work is supported by the European Union Seventh Fraework Progra Research for SME s with full title High Productivity Electron Bea Melting Additive Manufacturing Developent for the Part Production Systes Market and grant agreeent nuber We also gratefully thank the Geran Research Foundation (DFG) for funding the experiental work within the Collaborative Research Centre 814, project B2. References 1. P. Heinl et al., Adv. Eng. Mater (2007). 2. C. Körner et al., J. Mater. Process. Tech (2011). 3. M. Markl et al., Procedia Coput. Sci (2013). 4. R. Aer et al, Coput. Math. Appl (2014) 5. C. Feichtinger et al., J. Coput. Sci (2011). 6. H. Köstler et al., IT-Inforation Technology (2013). 7. G.R. McNaara et al., Phys. Rev. Lett (1988). 8. U. Frisch et al., Phys. Rev. Lett (1986).

9 Validation Experients for LBM Siulations of Electron Bea Melting 9 9. F. Higuera et al., Europhys. Lett (1989) 10. F.J. Alexander et al., Phys. Rev. E 47 R2249 (1993) 11. Y. Chen et al., Phys. Rev. E (1994) 12. F. Massaioli et al., Europhys. Lett (1993) 13. X. Shan, Phys. Rev. E (1997) 14. X. He et al., J. Coput. Phys (1998). 15. P. Lalleand et al., Int. J. Mod. Phys B (2003) 16. G.R. McNaara et al., J. Stat. Phys (1995). 17. P. Bhatnagar etal., Phys. Rev (1954). 18. S. Chen et al., Annu. Rev. Fluid Mech (1998). 19. X. He et al., Phys. Rev. E (1997) 20. X. He et al., Phys. Rev. E 55 R6333 (1997) 21. C. Körner et al., J. Stat. Physics (2005).