<figure><img src="image_1.png" /></figure>

HAL
open science

# Design and Assessment of an Austenitic Stainless Alloy for Laser Powder Bed Additive Manufacturing

Mariam Assi, Julien Favre, Marcin Brykala, Franck Tancret, Anna Fraczkiewicz

► To cite this version:

Mariam Assi, Julien Favre, Marcin Brykala, Franck Tancret, Anna Fraczkiewicz. Design and Assessment of an Austenitic Stainless Alloy for Laser Powder Bed Additive Manufacturing. Applied Sciences, 2024, 14 (19), pp.8649. 10.3390/app14198649. emse-04724656

HAL Id: emse-04724656

https://hal-emse.ccsd.cnrs.fr/emse-04724656v1

Submitted on 7 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Design and Assessment of an Austenitic Stainless Alloy for Laser Powder Bed Additive Manufacturing

Mariam Assi<fnref n="1" />, Julien Favre <fnref n="1" />,*, Marcin Brykala <fnref n="2" />, Franck Tancret <fnref n="3" /> and Anna Fraczkiewicz <fnref n="1" />,∗

**Abstract:** Recent developments in metallic additive manufacturing (AM) processes for the production of high-performance industrial pieces have been hampered by the limited availability of reliably processable or printable alloys. To date, most of the alloys used in AM are commercial grades that have been previously optimized for different manufacturing techniques. This study aims to design new alloys specifically tailored for AM processes, to minimize defects in the final products and to optimize their properties. A computational approach is proposed to design novel and optimized austenitic alloy compositions. This method integrates a suite of predictive tools, including machine learning, calculation of phase diagrams (CALPHAD) and physical models, all piloted by a multi-objective genetic algorithm. Within this framework, several material-dependent criteria are examined and their impact on properties and on the occurrence of defects is identified. To validate our approach, experimental tests are performed on a selected alloy composition: powder is produced by gas atomization and samples are fabricated by laser powder bed fusion. The microstructure and mechanical properties of the alloys are evaluated and its printability is compared with a commercial 316L stainless steel taken as a reference. The optimized alloy performs similarly to 316L in terms of coefficient of thermal expansion, hardness and elongation, but has a 17% lower yield strength and ultimate tensile strength (UTS), indicating that further optimization is required.

**Keywords:** alloy design; computational thermodynamics; laser powder bed fusion (LPBF); austenitic stainless alloys

## 1. Introduction

Laser powder bed fusion (LPBF) is becoming increasingly prominent in the aerospace and automotive industries, with a rising demand for materials specifically designed for this process. In the first two decades of metal additive manufacturing (AM) development, the focus on materials primarily involved using conventional grades [1]. Most of the literature has focused on the investigation of existing alloys that show promising results when fabricated using AM techniques. Examples include 304L and 316L in the stainless steel category, AlSi10Mg and AlSi12 in aluminum-based alloys, Ti-6Al-4V in titanium-based alloys, and In718 in nickel-based alloys, among others [2]. Despite the promising results obtained with these alloys in AM, their production requires optimizing a wide range of processing parameters. Indeed, the additive nature of these processes and the associated complex thermal histories introduce many difficulties, making it challenging to produce components that are both dense and structurally robust. Porosity, cracks, residual stresses and microstructure in the as-solidified state can significantly degrade the mechanical properties of AM-produced parts compared with conventionally manufactured pieces [3]. Many alloys are associated with a specific or preferred processing route; for instance, there are nickel base superalloys specifically optimized to be produced by forging (e.g., 718 grade),

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Assi, M.; Favre, J.; Brykala, M.; Tancret, F.; Fraczkiewicz, A. Design and Assessment of an Austenitic Stainless Alloy for Laser Powder Bed Additive Manufacturing. *Appl. Sci.* **2024**, *14*, 8649. https://doi.org/10.3390/app14198649

Academic Editors: Jana Petru, Marek Pagáč, Tomasz Kozior, Jiří Hajnyš and Jakub Mesicek

Received: 27 August 2024

Revised: 18 September 2024

Accepted: 19 September 2024

Published: 25 September 2024

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Mines Saint-Etienne, CNRS, UMR 5307 LGF, Centre SMS, 42023 Saint-Etienne, France</footnote>

<footnote n="2">² National Centre for Nuclear Research, NCBJ, 05-400 Świerk-Otwock, Poland; marcin.brykala@ncbj.gov.pl</footnote>

<footnote n="3">³ Institut des Matériaux de Nantes Jean Rouxel (IMN), Nantes Université, CNRS, 44000 Nantes, France; 
franck.tancret@univ-nantes.fr</footnote>

<footnote>* Correspondence: julien.favre@emse.fr (J.F.); fraczkie@emse.fr (A.F.)</footnote>
casting (e.g., 738), directional solidification (e.g., CMSX-4) or powder metallurgy (e.g., N18), so similar efforts should be undertaken towards the design of AM-specific alloys [4]. Nevertheless, few studies have focused on optimizing alloy composition to improve the overall properties of AM products and most alloy design efforts have been devoted primarily to aluminum and its alloys [5,6]. Improving process robustness, minimizing defects and achieving superior mechanical properties can be undertaken by designing new alloys tailored for AM. This, in turn, has the potential to significantly improve the performance of 3D printed parts. In recent years, several authors have presented successful approaches to design novel compositions using a combination of different computational tools [7–10]. This approach includes the use of computational thermodynamics to predict phases [11] and other physical models to predict material properties [12,13]. However, depending on the specific property of interest, physical models may not always be readily available. In such cases, the use of machine learning (ML) algorithms to build and train models using existing data has been widely adopted [14,15].

Most studies on alloy design and additive manufacturing focus primarily on the acceleration of alloy development by additive manufacturing. This is achieved by producing small quantities of materials for rapid screening and characterization of a wide range of compositions, as in the case of steels [16], or by producing out-of-equilibrium microstructures with unique properties, as in the case of titanium [17] or aluminum alloys [6] for additive manufacturing. Research dedicated to the development of alloys for additive manufacturing is less common and less well documented, with ongoing exploratory work on the precise methodologies and design criteria to be employed. Previous studies have often used a general approach to alloy selection and sorting [18,19]. However, these studies usually do not rely on a comprehensive multi-criteria exploration. In a recent development, Ackers et al. designed titanium alloys using thermodynamic calculations [20], considering constraints to prevent solidification cracking. They specifically tested compositions suitable for LPBF while maintaining associated process conditions and while considering surface defects, phase structure and cracking susceptibility. A similar approach was developed by Mishra and Thapliyal for the design of crack-free aluminum alloys, by reducing the solidification range [18,21]. Similarly, Tang developed a computational method to design new nickel superalloys that are resistant to hot cracking by minimizing the solidification range using the Scheil model [22]. In addition, they incorporated criteria to reduce strain-age cracking and improve mechanical strength and creep resistance. Their approach is multi-objective, focusing on the selection of alloy compositions by setting bounds on various property tables. However, the global objective space has not been fully explored using an optimization algorithm, suggesting the potential for further improvements through such exploration.

It was only with the work of Sabzi et al. [23,24] that a comprehensive multi-objective exploration of alloy compositions for additive manufacturing processes was established and successfully applied. They focused on the printability of stainless steel, proposing a defect prevention method that integrates material-dependent properties. Their approach, validated on 316L stainless steel, enhanced printability when compared with other alloys, but was limited to compositions near the nominal composition of 316L. Despite promising results in existing studies, a method is needed that focuses on material composition while addressing various material-dependent phenomena that are related to defects in current AM parts, including porosity formation, cracking, surface roughness and element loss by selective evaporation. A comprehensive alloy design methodology is still lacking, particularly in three key areas. First, an exhaustive list of optimization criteria to reduce defects, better control microstructure, and improve melt pool stability is lacking. Second, no global optimization has been performed over a wide range of compositions, making it difficult to identify optimal alloys that may differ significantly from commercial alloys. This work addresses these two issues by proposing a general methodology and validating it experimentally. The discussions will provide valuable insights into potential improvements in AM design techniques.
This work aims to contribute to the current progresses in alloy design approaches for the discovery of new austenitic alloys optimized for additive manufacturing processes. These alloys may belong to the family of austenitic stainless steels or austenitic high entropy alloys. In Section 1, we establish optimization criteria to guide the alloy design process. These criteria fall into two categories. The first group consists of criteria specific to the LPBF process, aimed at minimizing the risk of defects. A second group of criteria is related to the microstructure and the desired mechanical properties, with a focus on maintaining properties suitable for mechanical engineering applications. These criteria can be determined through the combined use of computational thermodynamics, physical models and machine learning algorithms for predicting relevant alloy properties. In the second part, we integrate these criteria into a multi-objective genetic algorithm (GA) to automate the optimization process. This algorithm ultimately proposes a set of optimized alloy compositions. Finally, we validate the models through experimental evaluation. An alloy composition selected from the optimized set is produced in powder form, processed using an AM technique (LPBF) and evaluated experimentally for defects, microstructure, as well as physical and mechanical properties. The resulting microstructure and properties are compared with those of 316L, selected as a reference alloy, to determine the potential improvements in printability and performance resulting from the design process.

## 2. Materials and Methods

The powder of the selected alloy was manufactured externally by LERMPS UTBM using gas atomization under argon atmosphere. Particle size was assessed using a laser granulometer Malvern Mastersizer 3000 (Malvern Panalytical, Malvern, UK), with a gas pressure of 3.5 bars. The particle size ranges from $d_{10} = 11$ µm to $d_{90} = 40$ µm, with a median value of $d_{50} = 22$ µm. Using a ProX200 LPBF machine (3D Systems, Rock Hill, SC, USA), cubic (10 × 10 × 10 mm) and rectangular specimens (42 × 14 × 4 mm) for tensile tests were fabricated on a 316L stainless steel substrate. A near-infrared laser source ($\lambda = 1064$ nm) was used. The hatching distance was set to 60 µm and the powder layer thickness was 30 µm. The scanning strategy was 45°/225° relative to the orthogonal axis of the substrate. The focal distance was -6.0 mm. Specimens were arranged in staggered rows on the building platform, with a building sequence oriented in the opposite direction of the Ar gas flow to prevent spatter deposition. All processing conditions were based on ranges optimized for 316L [25] to evaluate the robustness of the new alloy. Cubic test pieces were produced using two scan speeds (v) of 1000 mm/s and 1200 mm/s, with the laser power (P) gradually increased from 120 W to 250 W in increments of 6 W (+2%). As a first indicator of material quality, the surface roughness of the top surface was assessed using optical microscopy. Surface roughness is typically correlated with porosity, as both are usually associated with lack of fusion or unstable melt flow [26]. This technique has proven successful in a previous study for the rapid evaluation of new alloys via LPBF. Using this method, the optimal surface roughness was achieved with a laser power of 144 W and a scan speed of 1200 mm/s, and these parameters were used for the remainder of the study. The chemical composition of the powder and samples produced by LPBF was measured by inductively coupled plasma optical emission spectroscopy ICP-OES. The carbon content was measured by combustion infrared absorption (CIR). The values are listed in Table 1, showing some difference in the powder composition from the nominal values, especially with slightly lower amounts of some alloying elements, like Co and Al. The Co content further decreases slightly after LPBF processing.

<table><caption>**Table 1.** Composition of the selected optimized alloy.</caption><thead><tr><th>wt. %</th><th>Fe</th><th>Cr</th><th>Ni</th><th>Mn</th><th>Co</th><th>Al</th><th>C</th></tr></thead><tbody><tr><td>Nominal</td><td>42</td><td>15.2</td><td>13.7</td><td>11.7</td><td>15</td><td>2.3</td><td>0.02</td></tr><tr><td>Powder</td><td>45.9</td><td>14.3</td><td>13.1</td><td>11</td><td>13.5</td><td>2.1</td><td>0.0064</td></tr><tr><td>LPBF build</td><td>42.5</td><td>15.2</td><td>14.2</td><td>13.9</td><td>12</td><td>2.1</td><td>0.0064</td></tr></tbody></table>
Microstructure observations were conducted on a scanning electron microscope (SEM) Zeiss SUPRA55VP (Zeiss, Oberkochen, Germany) operated at 20 kV in backscattered electrons (BSE) mode and for Electron Back-Scattered Diffraction (EBSD) mapping. Variations in chemical composition were assessed using Energy Dispersive X-ray Spectrometry (EDX) mapping, with an Oxford SDD detector XmaxN 80 (Oxford Instruments, Abingdon-on-Thames, UK). No significant variations in chemical composition were detected within the limits of the EDX technique, indicating the absence of major compositional variations in the microstructure. However, it is possible that minor variations not detectable by this technique may still be present. Phase identification was accomplished using an X-ray Diffractometer (XRD) X'Pert Pro MPD Panalytical (Philips/Panalytical, Malvern, UK) with Cu Kα radiation operated at 40 mA and 45 kV. Divergence slits at 0.5° were set for the incident beam and a graphite monochromator with a Miniprop point detector were used for the diffracted beam. The analyzed angle was between 20° and 120° with a step size of 0.03° and a counting time of 25 s/step.

Thermal expansion tests were performed using a DIL 402 Expedis Select Netzsch (Netzsch, Selb, Germany) dilatometer in a horizontal configuration. Specimens of size 3 × 3 × 12 mm were placed in the alumina sample holder and tested in a temperature range from 25° to 850 °C, under a protective argon gas atmosphere (flow rate of 40 mL/min). Measurements were performed at a heating rate of 3 °C/min. For each material, three measurements were performed, which demonstrated a high repeatability, with a variation of approximately 1% across experiments. An alumina reference sample was used for calibration. The measurements were performed according to the standard recommendations [27].

Tensile tests were performed on flat specimens at room temperature using an Instron 1186 machine at a strain rate of 10⁻³ s⁻¹. Specimens were cut from 42 × 4 × 15 mm³ plates in two orientations: one along the scanning direction (SD) and the other along the building direction (BD). Specimens had a reduced cross-sectional area of 16 mm² and a length of 12 mm. They were tested in the as-machined condition and strain was measured using a clip-on extensometer with an initial gauge length of 10 mm (Instron, Norwood, MA, USA). Hardness tests were performed using a Wolpert Testwell (Testwell, Gennevilliers, France) machine, applying a 10 kg load with a dwell time of 10 s.

## 2.1. Criteria for Alloy Design

The calculated criteria should, whenever possible, accurately reflect the behavior of the material or specific trends, ideally independent of processing conditions. If complete independence is not possible, they should at least be applicable over a wide range of processing conditions. Then, it is desirable to gain a physical understanding of the impact of each feature on the overall outcome. The design criteria can be categorized into four groups. First, a criterion is set to manage the microstructure and ensure the stability of the desired phases. Secondly, a set of criteria is established to prevent the occurrence of solidification cracking. The third group focuses on maintaining melt pool stability, aiming to minimize melt instabilities that can lead to gas entrapment and subsequent pore formation. Lastly, several thermomechanical properties must be controlled, including yield stress and the coefficient of thermal expansion, to ensure a good strength of the material as well as to control the level of distortions and residual stresses.

### 2.1.1. Phase Stability: Targeting an Austenitic Structure

To illustrate and demonstrate the method to design alloys that are specific to AM, we chose to propose austenitic stainless alloys, with the primary objective being the predominant presence of austenite (γ, i.e., a face-centered cubic (FCC) solid solution) in the microstructure. Computational thermodynamics, with the calculation of phase diagrams (CALPHAD) method, can be employed to predict phase formation under equilibrium conditions. However, additive manufacturing generates out-of-equilibrium microstructures and CALPHAD may become inaccurate for such processes, somehow similar to welding. One possible solution is to use the Scheil–Gulliver model [28,29] to predict the solidifica-
tion sequence. This model, based on the CALPHAD method, assumes non-equilibrium conditions in which the solid phase forms instantaneously without diffusion, while the liquid phase stays fully homogeneous and in equilibrium with the last solid formed. It provides predictions of the composition and amount of each phase as the material cools and solidifies. However, this model is not well-suited for predicting the phase fractions after complete cooling. For instance, the Scheil-Gulliver model could predict the formation of δ ferrite at very high temperature, whereas that phase could transform into austenite in solid state during subsequent cooling; similarly, austenite formed at high temperature may transform into α ferrite or martensite upon cooling. In the present work, an empirical approach was preferred.

Therefore, the prediction of phases resulting from AM was approached using an empirical diagram usually applied to welding [30]. This relies on the quantification of γ-stabilizing elements and of elements promoting body-centered cubic (BCC) ferritic phases (α, δ), through quantities respectively named the “nickel equivalent” (Ni<sub>eq</sub>) and the “chromium equivalent” (Cr<sub>eq</sub>). One of the most frequently employed graphic representations for this purpose is the Schaeffler-Delong diagram [30,31]. Yet, in our specific case, in order to account for the influence of as many alloying elements as possible, we consider the Hull equations and diagram which are an extension of the former [32], where contents are expressed in wt.%, as follows:

$$ \begin{aligned} \text{Cr}_{\text{eq}} &= \text{Cr} + 1.21\text{Mo} + 0.48\text{Si} + 2.27\text{V} + 0.72\text{W} + 2.2\text{Ti} + 0.14\text{Nb} + 0.21\text{Ta} + 2.48\text{Al} \\ \text{Ni}_{\text{eq}} &= \text{Ni} + 0.11\text{Mn} - 0.0086\text{Mn}^2 + 0.14\text{Co} + 0.44\text{Cu} + 18.4\text{N} + 24.5\text{C} \end{aligned} $$

This approach can serve as a straightforward tool for estimating the phases that will be present after solidification and cooling to room temperature. It aids in a quick preliminary classification of compositions into two categories: potentially acceptable alloys (referred to as “feasible”) and unacceptable ones (“unfeasible”).

In Figure 1, the diagram illustrates the feasible space, defined with green dashed lines. The objective is to narrow down the search space to alloys that are characterized by a fully austenitic microstructure after solidification and cooling to room temperature. The feasible space is defined by two boundaries: one setting a limit on ferrite content and another on martensite content. Unlike martensite, which is to be completely avoided due to its brittleness, there is a small margin of around 5% allowed for the boundary with ferrite. Indeed, it is generally agreed that alloys containing less than 5% ferrite exhibit properties that are undistinguishable from those of fully austenitic ones. Figure 1 displays the distribution of several existing alloys on the Hull diagram, with ferrite content ranging from 0% to 5%. These alloys include 304, 316, SSW2 [33], 316L-5Mo-5Ni [34] and the so-called Cantor alloy, an equimolar CoCrFeMnNi “high entropy alloy” (HEA) [35]. Furthermore, to encourage the search for alloys that preferably reside within the austenitic region, the distances of the alloy composition from both the ferrite and martensite boundaries (illustrated as d<sub>m</sub> and d<sub>f</sub> from a random alloy X in Figure 1) is defined. An objective for alloy design is then set to maximize the smallest of the two distances. In conclusion, this criterion regarding phase stability combines a constraint on composition to ensure the presence of austenite and an objective to maximize its stability.

To further validate the austenitic nature of the alloy, the output prediction of the Scheil model is also considered. The simulation is conducted using Thermo-Calc 2021b with the TCFE10 database and TC-Python API. Although this model is primarily used to predict hot cracking, as discussed in the following section, it also provides insight into the phase fractions after solidification. To address the risk of the excessive formation of undesirable brittle phases and intermetallics, which could potentially affect mechanical properties, a second constraint has been implemented. This constraint ensures that δ-ferrite and austenite together constitute at least 99% of the phases after solidification, effectively rejecting alloy compositions containing more than 1% of other phases.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Hull diagram with a definition of the feasible search space defined in our study. A stands for austenite, F for ferrite and M for martensite.</figcaption></figure>

### 2.1.2. Solidification Cracking

The issue of solidification cracking has been extensively investigated in conventional manufacturing techniques such as casting and welding, but it also occurs in the specific case of LPBF process [36]. One of the key factors influencing the susceptibility of an alloy to cracking is its solidification temperature range. Therefore, calculating this value serves as a means to quantify the cracking susceptibility of an alloy [18]. More specifically, it is crucial to consider the critical temperature range (CTR) [37,38]. Lowering this range in the mushy zone during solidification minimizes shrinkage, reduces liquid suction in the interdendritic spacing and ultimately lowers the risk of void and crack formation. Definitions vary across sources, but, generally, the CTR is defined as the difference between temperatures at which certain fractions of solids have formed. This study adopts the definition corresponding to fractions of solids between 95% ($T_1$) and 99% ($T_2$), as these values are commonly reported [13,23,24]. To predict the solidification range and associated phases, including their order and the evolution of their fractions, the classical Scheil model from Thermo-Calc is employed. This computational tool enables the calculation of the complete solidification range. This model is only applied to alloys predicted by the Hull diagram to be austenitic, whereas alloys judged unsatisfactory by the Hull criterion are excluded from the Scheil calculation in order to limit computation time.

When considering processes with high solidification rates, various studies have demonstrated that the primary formation of ferrite from the liquid can significantly reduce the susceptibility to solidification cracking compared with the primary formation of austenite. Such a beneficial effect is observed when solidification produces ferrite only (F mode) or ferrite then austenite (FA mode) compared with the production of austenite only (A mode) or austenite then ferrite (AF mode) [39–42]. An important factor contributing to this effect is the higher solubility of harmful impurities (S, P...) in ferrite compared with that in austenite. This results in less segregation during the primary ferrite solidification. The cracking behavior of several austenitic stainless steels is influenced by the impurity content of S and P, as well as the $Cr_{eq}/Ni_{eq}$ ratio [43]. Bollinghaus et al. [44] noted that, as the $Cr_{eq}/Ni_{eq}$ ratio approaches 1.48, the impurity content of S and P, which leads to cracking, becomes significantly higher. Beyond this value (~1.48), cracking becomes rare, even with somewhat elevated impurity contents (>0.2 wt.%). This value of $Cr_{eq}/Ni_{eq}$ corresponds to the range where the solidification mode changes from A or AF to FA. To ensure that the solidification
mode is FA or F, it would be possible to use a minimum value of the Cr<sub>eq</sub>/Ni<sub>eq</sub> ratio as a design constraint. Instead, our choice was to impose the first phase formed during solidification as predicted by Thermo-Calc, which is set to δ-ferrite. Moreover, using Thermo-Calc with the Scheil model makes it possible to estimate the fraction of ferrite formed, unlike the Cr<sub>eq</sub>/Ni<sub>eq</sub> ratio alone. As a precaution, a minimum of 10% δ-ferrite formation is specified to reduce the risk of cracking [45]. However, with the ultimate desired microstructure being austenitic, it is also useful to establish an upper limit on the fraction of δ-ferrite formed during solidification, in addition to the Hull diagram criterion, to avoid an incomplete transformation of ferrite into austenite during cooling. The assessment of several existing alloys indicates that, beyond the 50% of ferrite predicted in a Scheil simulation, a risk exists that the microstructure will not be fully austenitic after welding or AM.

### 2.1.3. Melt Pool Composition and Stability

AM processes involve the use of liquid metal, making it imperative to examine various properties of materials in this state. Among these properties, surface tension (ST) plays a critical role and significantly impacts the wetting conditions of the previously deposited solid phase by the liquid phase. The value of ST influences the ability of the melt pool to form a continuous and stable seam. At high scanning speeds and with low deposited energy, the length of the melt pool increases relative to its diameter, leading to the phenomenon known as the “balling effect” [2]. This happens when the melt zone, which can be roughly seen as a “cylinder”, becomes unstable and splits into beads for a given ratio of its diameter over its length. The melt cylinder is attached to the solid substrate by a contact band; however, the physical phenomenon can be compared with some extent to a free-falling liquid jet and the resulting stability criterion differs only of a constant geometric factor [46,47]. A material-dependent criterion can be determined by considering the stability of a free-falling jet over time. Equation (1) [48] provides the expression of the characteristic time $\tau$ for the breakup of a fluid jet into drops:

$$ \tau \propto \sqrt{\frac{\rho R^3}{\gamma}} \qquad (1) $$

where R denotes the jet radius, $\rho$ signifies fluid density and $\gamma$ stands for the surface tension value. It can be inferred that, reducing the ST value would result in an extended breakup time, potentially longer than the solidification time, thus mitigating the balling effect. This decrease in ST would, overall, enhance the geometric stability of the seam. Due to the absence of a reliable physical model with which to predict the surface tension of metallic alloys as a function of composition, a general machine learning (ML) approach was employed. The model relies on a Bayesian algorithm, specifically Gaussian process regression (GPR); the method and the assessment of the model are detailed elsewhere [49].

Another phenomenon contributing to melt pool instability is vaporization. As reported by King et al. [50], the vaporization of the melt can lead to the occurrence of keyholing. This may result in the formation of pores if gas bubbles become trapped [51]. Reducing the risks associated with such instabilities can be achieved by lowering the vaporization flux of the elements. Additionally, selective vaporization may induce significant composition changes as shown by Mukherjee et al. [52]. To assess the risks related to instabilities arising from vaporization and composition changes, the vaporization flux of each element is estimated. Equation (2) shows the Langmuir expression for the vaporization flux ($J_i$) of element i [53]:

$$ J_i = \lambda_c \frac{P_i}{\sqrt{2\pi M_i T}} \qquad (2) $$

where $P_i$ and $M_i$ are the equilibrium vapor pressure and molecular weight of element i, respectively, T is the temperature and $\lambda_c$ is a positive constant that accounts for the condensation, typically assigned a value of 1 under vacuum conditions. Estimating $J_i$ for individual elements requires the determination of their respective vapor pressure values.
One approach to achieve this is by estimating it at a specific temperature using Alcock equations [54] in the form of a polynomial, as follows:

$$ \log P_i = A_i + \frac{B_i}{T} + C_i \log T + \frac{D_i}{T^3} \quad (3) $$

where $A_i$, $B_i$, $C_i$ and $D_i$ are coefficients specific to the i-th element, with commonly used values sourced from [54]. Equation (3) is applicable within the indicated temperature range and the fitted coefficients depend on the state of the material (solid or liquid). Yet, coefficients are often available over only a limited temperature range, and sometimes for the solid or liquid state only, the latter being the one needed. To overcome this problem, we will assume that the equations can be extended beyond the specified temperature ranges, and that they hold for both solid and liquid states. To verify these hypotheses, the vapor pressure calculated from Alcock equation at temperatures above the provided temperature range were compared in some cases for which the necessary data exists. For instance, aluminum has a vapor pressure of 2 atm at 2610 °C [55], compared with 1.97 atm calculated using the Alcock equation. Additionally, the equations for the liquid state are available for several metals, but for other ones (W, Mo, Cr, Nb and Mn), only equations for the solid can be found. However, when equations are found for both the liquid and solid metal, they seem to give similar values over the range of temperatures of interest. Consequently, it is here estimated that the use of equations established for solid state, when equations for the liquid are not found for the concerned metal, will still give a good estimation. In our study, we set both the average vaporization flux of the elements and the standard deviation of their corresponding values as objectives to be minimized at the boiling temperature of the alloy (estimated by a rule of mixtures), in order to avoid porosity through keyholing and to avoid compositional changes due to differential vaporization, respectively.

### 2.1.4. Physical and Mechanical Properties

* Maximization of the solid solution hardening effect

For crystalline materials, yield stress can be increased by mechanisms such as the grain size effect (Hall–Petch), precipitation hardening, strain hardening or solid solution hardening. In this particular study, we have ignored the effect of precipitation because we focus on the design of single-phase alloys. The influence of grain size was also omitted from the optimization criteria, as it is assumed to be mostly process dependent. Strain hardening is also ignored as we primarily focus on AM—i.e., non-deformed alloys. A remaining possible level of action to increase alloy strength is to maximize solid solution hardening. Labusch has suggested that the solid solution hardening can be represented as follows [56]:

$$ \Delta\sigma_{SSH} = \sum_{i} B_{i}X_{i}^{2} \quad (4) $$

where $X_i$ is the atomic fraction of component i and $B_i$ is a constant defined by the following [12,57]:

* $B_i = 3\mu\epsilon_i^{(4/3)}Z$, with $\mu$ the shear modulus, $\epsilon_i$ the mismatch parameter and $Z$ a fitting constant;

* $\epsilon_i = (\eta_i^2 + \alpha^2 \delta_i^2)^{0.5}$, with $\eta_i$ the term of the elastic misfit.

Using Vegard's law [58], for Fe-based alloys, one can estimate the atomic size misfit $\delta_i$ by the following:

$$ \delta_i = \frac{dr}{rdX_i} = \frac{r_i - r_{Fe}}{(1 - X_i)r_{Fe} + X_i \cdot r_i} \quad (5) $$

with $r_i$ representing the atomic radius of the solute and $r_{Fe}$ the atomic radius of the iron solvent. As for the term of the elastic misfit, in general, the variation of shear modulus with composition is assumed to be linear [59] and so the parameter $\eta_i$ can be estimated as follows:
$$ \eta_i = 2 \frac{\mu_i - \mu_{Fe}}{\mu_i + \mu_{Fe}} \qquad (6) $$

where $\mu_{Fe}$ and $\mu_i$ are the shear moduli of solvent Fe and solute i, respectively. Finally, by considering a value of $\alpha = 16$ which is typically used for FCC alloys to account for the interaction of solute atoms with edge dislocations [60], it is possible to calculate a solid solution strengthening index according to Equations (4)–(6). For alloy design, the solid solution index $\Delta\sigma_{SSH}$ is an objective to be maximized.

### Minimization of residual stress and distortions

The addition and melting of layers in LPBF directly affect residual stresses and distortions in the built component. These issues in AM components result from the spatial temperature gradient due to localized heating and cooling, leading to uneven strain distribution [61,62]. Assessing the risk of distortion and residual stresses is a very complex task but, for a given thermal history (as in thermal shock [63]), the extent of thermal strains scales, overall, with the coefficient of thermal expansion, $\alpha_{CTE}$. The stresses are then proportional to a product of strains and elastic constants. For different alloy elements i, with element atomic fraction $X_i$, $\alpha_{CTE}$ can be estimated using a simple rule of mixtures [64], where $\alpha_{CTE_i}$ represents the coefficient of thermal expansion for pure element i, as follows:

$$ \alpha_{CTE} = \sum_{i} X_{i} \alpha_{CTE_{i}} \qquad (7) $$

Elastic coefficients usually exhibit a limited variation, especially for austenitic stainless alloys (e.g., steels). Therefore, for alloy design it is here proposed to minimize $\alpha_{CTE}$ and to ignore any possible variation in elastic properties.

## 2.2. Optimization of Alloys with a Genetic Algorithm

### 2.2.1. Compositional Space Explored

Austenitic steels contain iron and primarily rely on Ni, Mn and Co as key alloying elements to stabilize the FCC solid solution, along with Cr to ensure corrosion resistance. In this study, we focus on iron-rich alloys, requiring a minimum of 30 wt.% Fe. To ensure good corrosion resistance and intergranular corrosion resistance, we set a minimum Cr content at 15 wt.% and a C content of 200 wt. ppm, similar to low carbon stainless steels. Indeed, increasing chromium and limiting carbon content prevents intergranular precipitation of Cr carbides. Nickel plays a crucial role in stabilizing the austenitic phase, therefore the upper limit for this element is set to a high value of 35% wt. Manganese and cobalt, both serving as austenite stabilizers, have an upper limit of 25%. Adjustments to the maximum concentrations of the remaining elements were made based on data from commercial alloy compositions [65]. The compositional search space is indicated in Table 2.

<table><caption>**Table 2.** The compositional search space considered for the optimization (in wt.%).</caption><thead><tr><th>Element</th><th>Fe</th><th>Cr</th><th>Ni</th><th>Co</th><th>Mn</th><th>Al</th><th>Ti</th><th>Nb</th><th>Mo</th><th>W</th><th>Si</th><th>C</th></tr></thead><tbody><tr><td>Min</td><td>30</td><td>15</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.02</td></tr><tr><td>Max</td><td>100</td><td>35</td><td>35</td><td>25</td><td>25</td><td>5</td><td>2</td><td>2</td><td>5</td><td>2</td><td>5</td><td>0.02</td></tr></tbody></table>

### 2.2.2. Genetic Algorithm for the Global Optimization

The multi-objective optimization was conducted using the NSGA-II algorithm [66], which was implemented using the pymoo Python package [67]. The optimization was performed on a population of 100 alloys, with the same number of offspring, to ensure a constant population size while keeping computation times reasonable. At each optimization step, new individuals were generated through simulated binary crossover with a distribution index set to $\eta = 20$. The algorithm also used polynomial mutation with a distribution index set to $\eta = 20$.
The objective functions—i.e., quantities that must be either maximized or minimized—and constraints—i.e., binary conditions that must be fulfilled—are summarized in Table 3 and were employed to evaluate and sort the various alloys within the population. The optimization algorithm followed a three-step process. In the first step, the Hull diagram was used to rapidly assess the potential of an alloy for being austenitic (constraint *a* in Table 3). If this criterion was not met, the alloy was penalized by setting low fitness scores for all the other criteria and non-feasible constraints, thereby saving computational time. In the second step, the Scheil model was applied using Thermo-Calc, and its outputs were used to estimate constraints *b* to *e* and objective *g* in Table 3. In the third step, additional objectives were calculated using GPR-ML for objective *h* and analytical formulas were used for objectives *f* and *i* to *l*. The convergence of the algorithm is measured by the replacement rate: beyond 40 generations it drops to a rather low value of 1% of the population per step, indicating a good convergence. The optimization loop was terminated after 50 generations, for a replacement rate of 0.5% per step.

<table><caption>**Table 3.** Summary of the different objective functions and constraints.</caption><thead><tr><th>Criteria</th><th>Constraint</th><th>Objective</th><th>Method of Calculation</th></tr></thead><tbody><tr><td rowspan="3">Austenitic structure</td><td>(a) Belongs to feasible space on Hull diagram</td><td>(f) Maximize distance from boundaries</td><td>Hull diagram</td></tr><tr><td>(b) δ &lt; 50%</td><td></td><td></td></tr><tr><td>(c) δ + γ &gt; 99% at the end of solidification</td><td></td><td>Scheil model</td></tr><tr><td rowspan="2">Solidification cracking</td><td>(d) Solidification mode = FA</td><td>(g) Minimize CTR</td><td rowspan="2">Scheil model</td></tr><tr><td>(e) δ &gt; 10%</td><td></td></tr><tr><td>Melt pool stability and balling</td><td></td><td>(h) Minimize ST</td><td>GPR-ML model</td></tr><tr><td>Thermal strains and residual stresses</td><td></td><td>(i) Minimize α<sub>CTE</sub></td><td>Rules of mixtures</td></tr><tr><td>Strength</td><td></td><td>(j) Maximize SSH</td><td>Labusch SSH model</td></tr><tr><td>Porosity</td><td></td><td>(k) Minimize total vaporization flux</td><td>Vaporization flux</td></tr><tr><td>Control of composition</td><td></td><td>(l) Minimize variance of vaporization flux</td><td>Vaporization flux</td></tr></tbody></table>

## 3. Results

### 3.1. *Selection of a Specific Alloy*

The distribution of compositions for the optimized population of alloys is illustrated in Figure 2 using a violin plot. Cr content ranges between 15% and 31 wt.%, while Ni falls within the range of 11% to 22 wt.%. Mn exhibits a broader range, spanning from 0% to 24 wt.%. Notably, a significant amount of Co, ranging from 5% to 20 wt.%, is prevalent in most of the alloys. This elevated cobalt content is likely a result of the criterion for enhancing solid solution hardening and stabilizing the austenitic structure. The content of other elements is generally limited to a few percent: Mo and Al do not exceed 3% and 4 wt.%, respectively. On the other hand, elements like Nb, Ti and W have average values well below a fraction of a percent and are considered negligible, despite their potentially high solid solution hardening capability. It seems possible that the strong tendency of these three elements towards primary segregation conflicts with the objective of minimizing the CTR value. Furthermore, these three elements act as ferrite stabilizers, ultimately raising the Cr<sub>eq</sub> value and rendering their inclusion unfavorable.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Violin plots of the optimized alloy population after 50 generations, with the median point (black diamond) and the selected composition (red dot).</figcaption></figure>

After optimization, all the alloys in the population belong to a Pareto front. Selecting a single optimal alloy is challenging as each performs better in some specific criteria. To determine the best overall option, we ranked the alloys by their average performance across all objectives. The preference ranking organization method for enrichment evaluation (PROMETHEE) method [68], a multi-criteria decision tool, was used to select the best alloys based on expressed preferences. The detailed methodology is explained elsewhere [69] but, to summarize, the user expresses the relative importance of each objective (or property), using a so-called preference function whose parameters are specific to each objective. This allows, based on said preferences, to establish a global ranking of alloys that are, by design in the present case, incomparable as they are all Pareto optimal, i.e., each alloy outperforms all the others on at least one objective. This results in a shortlist of alloys. At this stage, selecting an alloy remains challenging because all materials of the shortlist display rather close characteristics. To aid in selection, we introduced an additional criterion: the temperature range for austenite stability at equilibrium, predicted by Thermo-Calc. Although not initially included in the multi-objective optimization due to computational costs, it is essential for avoiding unwanted phase transformations during post-AM heat treatments. This criterion helped differentiate between otherwise equivalent alloys based on other factors. From the best performing alloys on the shortlist, one composition was selected based on this additional criterion. This composition is given in Table 1 and marked with red dots in Figure 2.

Table 4 presents a comparison of the values of our defined objective functions for the selected material, calculated based on the effective powder composition specified in Table 1, and for 316L. Overall, the optimized alloy outperforms 316L, except for the thermal expansion coefficient, which is slightly higher than that of this reference alloy.

<table><caption>**Table 4.** Calculated objective functions values for the optimized alloy and for 316L stainless steel.</caption><thead><tr><th></th><th>Hull min.<br>Distance (a.u.)</th><th>CTR<br>(K)</th><th>ST<br>[N.m<sup>-1</sup>]</th><th>SSH (a.u.)</th><th>α<sub>CTE</sub><br>[10<sup>-6</sup>/K]</th><th>Vap. Flux<br>Average</th><th>Vap. Flux<br>std.</th></tr></thead><tbody><tr><td>Optimized alloy</td><td>3.5</td><td>77</td><td>1.0</td><td>44</td><td>13</td><td>0.12</td><td>0.25</td></tr><tr><td>316L</td><td>1.3</td><td>92</td><td>1.3</td><td>31</td><td>11</td><td>0.13</td><td>0.34</td></tr></tbody></table>

The Scheil simulations of solidification and phase fractions with temperature for the selected alloy are depicted in Figure 3. The Scheil diagram in Figure 3a reveals that the alloy initially solidifies into δ-ferrite (BCC_A2) for a solid fraction up to 20%, then austenite is formed (FCC_A1). By the end of the solidification, these two phases collectively make up 99% of the material. The predicted CTR is relatively low, with a value of 77 K, compared with the value of 92 K for the 316L alloy. Consequently, the risk of solidification cracking is minimal. According to the equilibrium calculation of Figure 3b, the alloy is
expected to maintain a single FCC structure within the temperature range of 950 to 1230 °C. Below 950 °C, an ordered BCC phase of the B2 type is predicted to form at equilibrium. Nevertheless, the Cr<sub>eq</sub> and Ni<sub>eq</sub> values, in accordance with the Hull diagram, suggest the alloy is anticipated to remain fully austenitic after solidification and cooling in a welding-type process. A negligible amount of chromium-rich M<sub>23</sub>C<sub>6</sub> carbides (<0.005%) is predicted below approximately 800 °C.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** (a) Scheil model and (b) equilibrium Thermo-Calc predictions for the selected alloy.</figcaption></figure>

## 3.2. *Experimental Validation*

### 3.2.1. Evaluation of the Coefficient of Linear Thermal Expansion

Dilatometry measurements were performed for three materials: 316L stainless steel, the optimized alloy after casting, and the same optimized alloy produced by LPBF. Figure 4 shows the variation of the expansion coefficient for the different alloys as a function of temperature. Each point was determined by averaging the values within a range of -50 °C to +50 °C around the experimental point. The coefficient α systematically increases by about 30% with temperature and then stabilizes above 500 °C.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Variation of the thermal expansion coefficient α<sub>CTE</sub> with temperature.</figcaption></figure>

The average coefficient of expansion over the three measurements was found to be nearly identical for both the cast material (α<sub>CTE</sub> = 20.5 × 10<sup>-6</sup> K<sup>-1</sup>) and the LPBF material
(20.4 × 10⁻⁶ K⁻¹). The average value for 316L was slightly lower (19.7 × 10⁻⁶ K⁻¹). While the absolute measured values differ from the predictions and are slightly higher for both alloys, the relative variations between them remain consistent with the predictions (Table 4). In all cases, the experimental values exceeded the predicted ones, revealing the limitations of a linear mixture model in accurately predicting α<sub>CTE</sub>: while it can quickly estimate the relative variation in α<sub>CTE</sub> between different alloys, the accuracy of the predicted absolute values remains poor and they should not be used directly. Nevertheless, as trends seem correctly captured, such a simple tool can be used for alloy design if the property is to be optimized (e.g., minimized), which is the case here.

3.2.2. Microstructure after AM Processing

* Analysis of porosity

To determine the porosity, the alloy produced through LPBF is characterized by SEM
observations: a set of images is taken over a total surface of 4 mm². The images are then
thresholded to capture the pore fraction and size distribution using ImageJ software 1.54f
(Figure 5). An example is shown for the material built at P = 144 W and v = 1200 mm/s,
displaying an average porosity of 0.05%, with a standard deviation of 0.02%. This rate is
notably low compared with those observed in similar alloys, especially considering that
standard processing parameters were used without any specific adjustment. Typically, the
porosity rate reaches values of around 0.8–1% for common alloys like 316L [70,71]. Figure 6
illustrates the pore size distribution, which is roughly a lognormal distribution with an
average size of 4.5 µm and a median value of 3.4 µm. The pore size remains less than or
equal to 10 µm. This is below the critical defect size that can initiate crack propagation,
even under fatigue conditions. Indeed, assuming a very low crack propagation threshold
of 2 MPa·m¹/² to be conservative and a penny-shaped crack of radius 5 µm associated with
a conventional geometrical factor of 2/√π, cracking would be likely to occur only for a
stress amplitude superior to 560 MPa, which is higher than the yield stress. Given that a
pore is normally less critical than a crack and that the actual crack propagation threshold
more likely lies between 5 and 10 MPa·m¹/² for this category of alloys, it is almost certain
that the observed pores cannot become critical defects and initiate any damage.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** (a) Observation of the pores by SEM (red arrows), (b) thresholding of the image to detect pores, and (c) detection using ImageJ software.</figcaption></figure>

* Microstructure analysis

X-ray diffraction analysis was performed to identify the phases present in the powder
and in the LPBF built material (Figure 7). The observed peaks indicate the presence of
two phases in the powder. The primary FCC phase corresponds with the austenite phase
(γ), with a theoretical lattice parameter a = 3.606 Å, according to JCPDS card number
01-081-8775. Additionally, a minor BCC phase, estimated to constitute 3% of the mate-
rial, corresponds with the ferrite phase, characterized by a theoretical lattice parameter
a = 2.882 Å (JCPDS card number 04-003-2920). The LPBF built material characterized in the
same conditions shows the presence of only the austenite phase. The ferrite phase observed
in the powder is therefore not detected in the LPBF built samples. It is possible that the
ferrite phase was first formed during solidification, as predicted, and then transformed into austenite.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Pore size distribution obtained by SEM image analysis for a material built at P = 144 W, v = 1200 mm/s.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 7.** X-ray diffraction phase identification in as-atomized powder and LPBF bulk material. Green stars represent austenite (JCPDS 01-081-8775) and purple stars represent ferrite (JCPDS 04-003-2920).</figcaption></figure>

The crystalline orientation in the melt pools is characterized using EBSD. The inverse pole figure (IPF) map is shown in Figure 8a, illustrating the crystalline orientation along the building direction (BD) with the laser scanning direction (SD) along the X and Y axes. Large grains with melt pool geometry are observed and high angle boundaries (>15°) are represented by black lines. Small variations in the IPF color code within grains suggest the presence of a post-solidification substructure, indicating that a significant portion of the FCC phase formed directly during the solidification step. This observation is consistent with predictions from the Scheil model (Figure 3a), suggesting a likely liquid → γ solidification for much of the observed austenite. It is noteworthy that the microstructure does not consist of long columnar grains, as is commonly observed in the 316L alloy [72]. The crystalline orientation of each successive layer is different from the previous one, indicating that solidification does not occur by epitaxial growth between successive layers. Figure 8b shows an IPF map along BD, with BD normal to the figure and SD along the X and Y
axes. This map confirms the presence of large grains in the melt pool, resulting in a microstructure that follows a lattice pattern. The pole figure corresponding to Figure 8b is shown in Figure 8c for {100} planes, showing that this crystalline orientation is preferentially aligned along the BD and SD axes and is close to a cube texture.

<figure><img src="image_12.png" /><figcaption>**Figure 8.** IPF map indicating the crystalline orientation along BD for (a) the cross-section of melt pools along SD direction, (b) the cross-section along BD and (c) pole figure of {100} planes.</figcaption></figure>

### 3.2.3. Mechanical Properties

In its as-built condition, the average hardness of the alloy is of $216 \pm 6$ HV10, which falls within the upper range of typical hardness values for wrought 316L, usually between 170 and 220 HV [65], and is close to the hardness of 316L produced by LPBF, which has been reported at 231 HV [73]. Tensile tests were conducted at room temperature along two distinct orientations, with the tensile axis aligned either parallel to the building direction (BD) or to the scanning direction (SD) to evaluate a potential anisotropic behavior. Three tests were performed for each orientation. The test results (Figure 9) exhibit good reproducibility, showing yield stress differences below 5 MPa between specimens for both SD and BD orientations. The material tested along SD exhibits the highest strength, with an average yield stress of 520 MPa, an ultimate tensile stress (UTS) of 600 MPa and an elongation to fracture of 40%. In contrast, the material tested along BD displays a yield stress 6% lower (490 MPa) and a UTS 14% lower (525 MPa). These differences may be due to the slight anisotropy of the cellular dendrite microstructure of the samples, as observed in [74], along with the fact that metallurgical defects are more easily introduced into the bonding area between two consecutive layers [75]. While stress-strain curves in horizontally cut samples are quite reproducible, vertically cut samples exhibit notable discrepancies in elongation to fracture, suggesting potential microstructure variations. The mechanical properties of the optimized alloy can be compared with the reference 316L alloy. For example, in the Huang et al. work, a 316L produced by LPBF had a yield strength of 623 MPa, a UTS of 726 MPa, and an elongation of 47% [73]. Our optimized alloy offers a yield strength and UTS that are both 17% lower than 316L for the same elongation. Therefore, future development could focus on further increasing strength, for example by increasing the maximum allowable concentration of alloying elements during optimization.
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Engineering stress–strain curves for the tensile tests on horizontal (SD) and vertical (BD) as-built LPBF samples.</figcaption></figure>

The fracture surfaces of the tensile specimens were observed by SEM (Figure 10a). In all samples, ductile fracture is observed, with microvoids and dimples characteristic of ductile fracture. However, occasional small spherical inclusions were noted, as illustrated in Figure 10 with a red arrow. The EDX analysis of these inclusions (Figure 10b,c) show the presence of Al and O, indicating that they are probably alumina. Interestingly, these inclusions were not initially detected through SEM imaging of the powder. As aluminum exhibits a significantly higher affinity for oxygen compared with other alloying elements, the presence of aluminum oxides may result from a small oxygen fraction within the powder or the atmosphere of the laser chamber. Consequently, these inclusions are likely due to contamination during the process. To mitigate their occurrence, a better control over powder processing and storage conditions, coupled with monitoring and regulating the oxygen content within the LPBF chamber, may prove effective in reducing the incidence of such inclusions. However, the latter are always a few µm in diameter, so that, for the reasons already invoked, pores cannot be at the origin points of cracking under static or even cyclic loading.

<figure><img src="image_14.png" /><figcaption>**Figure 10.** (a) SEM image of the fracture surface and EDX maps for (b) Al and (c) O of the spherical inclusions referred to with red arrows.</figcaption></figure>

The challenge of powder oxidation could nevertheless be addressed in future alloy design research. The introduction of additional criteria may help to penalize elements that are particularly prone to oxidation. For example, one possible method could be to develop a penalty function based on the standard free enthalpies of formation of their respective oxides. This strategy would aim to minimize the risk of oxidation during laser melting and powder storage, thereby avoiding the formation of potentially embrittling (i.e., larger) oxide inclusions.
## 4. Conclusions

In this study, a genetic algorithm was used to optimize a single-phased, iron-rich austenitic alloy for additive manufacturing. The Hull diagram provided feasibility constraints and the Scheil model addressed solidification conditions. The Labusch model was used to evaluate and to increase solid solution strengthening, to promote improved mechanical resistance. Minimizing the coefficient of thermal expansion ($\alpha_{CTE}$), vaporization fluxes and reducing surface tension were key objectives to reduce porosity or distortions. The risk of cracking was also minimized via several criteria. These objectives and constraints were integrated into a multi-objective genetic algorithm, resulting in the identification and selection of a new optimized composition with improved printability. The optimized alloy powder was produced by gas atomization and then used in the LPBF process to produce specimens for microstructural and mechanical characterization.

The study leads to the following conclusions:

* The overall optimization method efficiently converges on a promising set of optimized alloys with expected attractive properties as predicted by the models.

* The selection of a single alloy proves to be challenging and requires the introduction of an additional criterion regarding the stability range of austenite to aid in the decision-making process.

* The selected alloy produced crack-free parts with a low porosity rate, even when printed using standard processing parameters. This confirms the effectiveness of the design process in eliminating critical defects and ensuring high ductility.

* The LPBF material consists of a single FCC phase. As predicted, its coefficient of thermal expansion is slightly higher than that of 316L stainless steel. Tensile tests showed good ductility, with strength values comparable to other LPBF-printed austenitic steels such as 316L. These experimental results confirm the accuracy of the CALPHAD model and other analytical models, and increase confidence in the alloy property predictions made during the optimization stage.

* The overall performance of the alloy was found to be quite comparable to 316L. Expansion coefficient measurements showed only a 3% difference from 316L. Mechanical properties were also comparable, with similar levels of hardness and elongation. However, the yield strength and UTS of the optimized alloy were 17% lower. Therefore, further improvements should be considered in the optimization process to surpass existing austenitic alloys such as 316L.

Future steps could consider additional criteria such as material cost, density, and spatter formation. New criteria could also focus on reducing sensitivity of the powder to oxidation; however, these criteria relate more to the powder itself than to the alloy composition alone, making it difficult to address from an alloy-only perspective. While the study demonstrates feasibility, targeting specific applications has promising industrial potential.

**Author Contributions:** Conceptualization, F.T., A.F. and J.F.; methodology, M.A., J.F., F.T. and A.F.; formal analysis, M.A., J.F., A.F. and F.T.; investigation, M.A. and M.B.; resources, A.F. and F.T.; writing—original draft preparation, M.A. and J.F.; writing—review and editing, J.F., F.T. and A.F.; supervision, J.F., F.T. and A.F.; funding acquisition: J.F. and A.F. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The raw data supporting the conclusions of this article will be made available by the authors on request.

**Conflicts of Interest:** All the authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
## References

1. Clare, A.T.; Mishra, R.S.; Merklein, M.; Tan, H.; Todd, I.; Chechik, L.; Li, J.; Bambach, M. Alloy Design and Adaptation for Additive Manufacture. *J. Mater. Process. Technol.* **2022**, *299*, 117358. [CrossRef]

2. DebRoy, T.; Wei, H.L.; Zuback, J.S.; Mukherjee, T.; Elmer, J.W.; Milewski, J.O.; Beese, A.M.; Wilson-Heid, A.; De, A.; Zhang, W. Additive Manufacturing of Metallic Components—Process, Structure and Properties. *Prog. Mater. Sci.* **2018**, *92*, 112–224. [CrossRef]

3. Sanaei, N.; Fatemi, A. Defects in Additive Manufactured Metals and Their Effect on Fatigue Performance: A State-of-the-Art Review. *Prog. Mater. Sci.* **2021**, *117*, 100724. [CrossRef]

4. Bandyopadhyay, A.; Traxel, K.D.; Lang, M.; Juhasz, M.; Eliaz, N.; Bose, S. Alloy Design via Additive Manufacturing: Advantages, Challenges, Applications and Perspectives. *Mater. Today* **2022**, *52*, 207–224. [CrossRef]

5. Pauzon, C.; Buttard, M.; Després, A.; Charlot, F.; Fivel, M.; Chehab, B.; Blandin, J.-J.; Martin, G. Direct Ageing of LPBF Al-1Fe-1Zr for High Conductivity and Mechanical Performance. *Acta Mater.* **2023**, *258*, 119199. [CrossRef]

6. Buttard, M.; Freixes, M.L.; Josserond, C.; Donnadieu, P.; Chéhab, B.; Blandin, J.-J.; Gault, B.; De Geuser, F.; Martin, G. Ageing Response and Strengthening Mechanisms in a New Al-Mn-Ni-Cu-Zr Alloy Designed for Laser Powder Bed Fusion. *Acta Mater.* **2023**, *259*, 119271. [CrossRef]

7. Tancret, F. Computational Thermodynamics, Gaussian Processes and Genetic Algorithms: Combined Tools to Design New Alloys. *Model. Simul. Mater. Sci. Eng.* **2013**, *21*, 045013. [CrossRef]

8. Tancret, F.; Toda-Caraballo, I.; Menou, E.; Rivera Díaz-Del-Castillo, P.E.J. Designing High Entropy Alloys Employing Thermodynamics and Gaussian Process Statistical Analysis. *Mater. Des.* **2017**, *115*, 486–497. [CrossRef]

9. Deschamps, A.; Tancret, F.; Benrabah, I.-E.; De Geuser, F.; Van Landeghem, H.P. Combinatorial Approaches for the Design of Metallic Alloys. *Comptes Rendus Phys.* **2018**, *19*, 737–754. [CrossRef]

10. Menou, E.; Toda-Caraballo, I.; Rivera-Díaz-del-Castillo, P.E.J.; Pineau, C.; Bertrand, E.; Ramstein, G.; Tancret, F. Evolutionary Design of Strong and Stable High Entropy Alloys Using Multi-Objective Optimisation Based on Physical Models, Statistics and Thermodynamics. *Mater. Des.* **2018**, *143*, 185–195. [CrossRef]

11. Tancret, F. Computational Thermodynamics and Genetic Algorithms to Design Affordable Γ'-Strengthened Nickel-Iron Based Superalloys. *Model. Simul. Mater. Sci. Eng.* **2012**, *20*, 045012. [CrossRef]

12. Sieurin, H.; Zander, J.; Sandström, R. Modelling Solid Solution Hardening in Stainless Steels. *Mater. Sci. Eng.* **2006**, *415*, 66–71. [CrossRef]

13. Menou, E.; Ramstein, G.; Bertrand, E.; Tancret, F. Multi-Objective Constrained Design of Nickel-Base Superalloys Using Data Mining- and Thermodynamics-Driven Genetic Algorithms. *Model. Simul. Mater. Sci. Eng.* **2016**, *24*, 055001. [CrossRef]

14. Guo, Z.; Sha, W. Modelling the Correlation between Processing Parameters and Properties of Maraging Steels Using Artificial Neural Network. *Comput. Mater. Sci.* **2004**, *29*, 12–28. [CrossRef]

15. Kaufmann, K.; Vecchio, K.S. Searching for High Entropy Alloys: A Machine Learning Approach. *Acta Mater.* **2020**, *198*, 178–222. [CrossRef]

16. Knoll, H.; Ocylok, S.; Weisheit, A.; Springer, H.; Jägle, E.; Raabe, D. Combinatorial Alloy Design by Laser Additive Manufacturing. *Steel Res. Int.* **2017**, *88*, 1600416. [CrossRef]

17. Zhang, T.; Liu, C.-T. Design of Titanium Alloys by Additive Manufacturing: A Critical Review. *Adv. Powder Mater.* **2022**, *1*, 100014. [CrossRef]

18. Mishra, R.S.; Thapliyal, S. Design Approaches for Printability-Performance Synergy in Al Alloys for Laser-Powder Bed Additive Manufacturing. *Mater. Des.* **2021**, *204*, 109640. [CrossRef]

19. Hunt, J.; Derguti, F.; Todd, I. Selection of Steels Suitable for Additive Layer Manufacturing. *Ironmak. Steelmak.* **2014**, *41*, 254–256. [CrossRef]

20. Ackers, M.A.; Messé, O.M.D.M.; Hecht, U. Novel Approach of Alloy Design and Selection for Additive Manufacturing towards Targeted Applications. *J. Alloys Compd.* **2021**, *866*, 158965. [CrossRef]

21. Thapliyal, S.; Komarasamy, M.; Shukla, S.; Zhou, L.; Hyer, H.; Park, S.; Sohn, Y.; Mishra, R.S. An Integrated Computational Materials Engineering-Anchored Closed-Loop Method for Design of Aluminum Alloys for Additive Manufacturing. *Materialia* **2020**, *9*, 100574. [CrossRef]

22. Tang, Y.T.; Panwisawas, C.; Ghoussoub, J.N.; Gong, Y.; Clark, J.W.G.; Németh, A.A.N.; McCartney, D.G.; Reed, R.C. Alloys-by-Design: Application to New Superalloys for Additive Manufacturing. *Acta Mater.* **2021**, *202*, 417–436. [CrossRef]

23. Sabzi, H.E.; Rivera-Díaz-del-Castillo, P.E.J. Defect Prevention in Selective Laser Melting Components: Compositional and Process Effects. *Materials* **2019**, *12*, 3791. [CrossRef] [PubMed]

24. Sabzi, H.E.; Maeng, S.; Liang, X.; Simonelli, M.; Aboulkhair, N.T.; Rivera-Díaz-del-Castillo, P.E.J. Controlling Crack Formation and Porosity in Laser Powder Bed Fusion: Alloy Design and Process Optimisation. *Addit. Manuf.* **2020**, *34*, 101360. [CrossRef]

25. Tucho, W.M.; Lysne, V.H.; Austbø, H.; Sjolyst-Kverneland, A.; Hansen, V. Investigation of Effects of Process Parameters on Microstructure and Hardness of SLM Manufactured SS316L. *J. Alloys Compd.* **2018**, *740*, 910–925. [CrossRef]

26. Qiu, C.; Panwisawas, C.; Ward, M.; Basoalto, H.C.; Brooks, J.W.; Attallah, M.M. On the Role of Melt Flow into the Surface Structure and Porosity Development during Selective Laser Melting. *Acta Mater.* **2015**, *96*, 72–79. [CrossRef]

27. ISO 17562; Fine Ceramics (Advanced Ceramics, Advanced Technical Ceramics)—Test Method for Linear Thermal Expansion of Monolithic Ceramics by Push-Rod Technique. ISO: Geneva, Switzerland, 2016.
28. Scheil, E. Bemerkungen Zur Schichtkristallbildung. *Int. J. Mater. Res.* **1942**, *34*, 70–72. [CrossRef]

29. Keller, T.; Lindwall, G.; Ghosh, S.; Ma, L.; Lane, B.M.; Zhang, F.; Kattner, U.R.; Lass, E.A.; Heigel, J.C.; Idell, Y.; et al. Application of Finite Element, Phase-Field, and CALPHAD-Based Methods to Additive Manufacturing of Ni-Based Superalloys. *Acta Mater.* **2017**, *139*, 244–253. [CrossRef]

30. Schaeffler, A.L. Constitution Diagram for Stainless Steel Weld Metal. *Met. Prog.* **1949**, *56*, 680.

31. Delong, W.T. Ferrite in Austenitic Stainless Steel Weld Metal. *Weld. J.* **1974**, *53*, 273.

32. Hull, F.C. Delta Ferrite and Martensite Formation in Stainless Steels. *Weld. J.* **1973**, *52*, 193.

33. Salama, E.; Eissa, M.M.; Tageldin, A.S. Distinct Properties of Tungsten Austenitic Stainless Alloy as a Potential Nuclear Engineering Material. *Nucl. Eng. Technol.* **2019**, *51*, 784–791. [CrossRef]

34. Indacochea, J.E.; Gattu, V.K.; Chen, X.; Rahman, T. *Performance of a Steel/Oxide Composite Waste Form for Combined Waste Steams from Advanced Electrochemical Processes*; University of Illinois: Chicago, IL, USA, 2017.

35. Cantor, B. Multicomponent and High Entropy Alloys. *Entropy* **2014**, *16*, 4749–4768. [CrossRef]

36. Sun, Z.; Tan, X.P.; Descosins, M.; Mangelinck, D.; Tor, S.B.; Lim, C.S. Revealing Hot Tearing Mechanism for an Additively Manufactured High-Entropy Alloy via Selective Laser Melting. *Scr. Mater.* **2019**, *168*, 129–133. [CrossRef]

37. Rappaz, M.; Drezet, J.-M.; Gremaud, M. A New Hot-Tearing Criterion. *Met. Mater. Trans. A* **1999**, *30*, 449–455. [CrossRef]

38. Galy, C.; Le Guen, E.; Lacoste, E.; Arvieu, C. Main Defects Observed in Aluminum Alloy Parts Produced by SLM: From Causes to Consequences. *Addit. Manuf.* **2018**, *22*, 165–175. [CrossRef]

39. Suutala, N. Effect of Manganese and Nitrogen on the Solidification Mode in Austenitic Stainless Steel Welds. *Met. Trans. A* **1982**, *13*, 2121–2130. [CrossRef]

40. Suutala, N. Effect of Solidification Conditions on the Solidification Mode in Austenitic Stainless Steels. *Met. Trans. A* **1983**, *14*, 191–197. [CrossRef]

41. Radhakrishnan, V.M. Hot Cracking in Austenitic Stainless Steel Weld Metals. *Sci. Technol. Weld. Join.* **2000**, *5*, 40–44. [CrossRef]

42. Nishimoto, K.; Mori, H. Hot Cracking Susceptibility in Laser Weld Metal of High Nitrogen Stainless Steels. *Sci. Technol. Adv. Mater.* **2004**, *5*, 231–240. [CrossRef]

43. Shankar, V.; Gill, T.P.S.; Mannan, S.L.; Sundaresan, S. Solidification Cracking in Austenitic Stainless Steel Welds. *Sadhana* **2003**, *28*, 359–382. [CrossRef]

44. Böllinghaus, T.; Herold, H.; Cross, C.E.; Lippold, J.C. *Hot Cracking Phenomena in Welds II*, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2008; ISBN 978-3-540-78627-6.

45. Saluja, R.; Moeed, K. Experimental Investigation Of Solidification-Mode And Response Surface Modeling Of Ferrite-Content In Grade 304L Pulse GMA Welded Plates. *Mater. Today Proc.* **2019**, *18*, 3876–3890. [CrossRef]

46. Gusarov, A.V.; Smurov, I. Modeling the Interaction of Laser Radiation with Powder Bed at Selective Laser Melting. *Phys. Procedia* **2010**, *5*, 381–394. [CrossRef]

47. Chandrasekhar, S. *Hydrodynamic and Hydromagnetic Stability*; Dover Publications: Mineola, NY, USA, 1981; ISBN 978-0-486-64071-6.

48. De Gennes, P.-G.; Brochard-Wyart, F.; Quéré, D. *Capillarity and Wetting Phenomena*; Springer: New York, NY, USA, 2004; ISBN 978-1-4419-1833-8.

49. Assi, M.; Favre, J.; Fraczkiewicz, A.; Tancret, F. Machine Learning-Based Model of Surface Tension of Liquid Metals: A Step in Designing Multicomponent Alloys for Additive Manufacturing. *J. Mater. Sci.* **2022**, *57*, 13466. [CrossRef]

50. King, W.E.; Anderson, A.T.; Ferencz, R.M.; Hodge, N.E.; Kamath, C.; Khairallah, S.A.; Rubenchik, A.M. Laser Powder Bed Fusion Additive Manufacturing of Metals; Physics, Computational, and Materials Challenges. *Appl. Phys. Rev.* **2015**, *2*, 041304. [CrossRef]

51. Martin, A.A.; Calta, N.P.; Hammons, J.A.; Khairallah, S.A.; Nielsen, M.H.; Shuttlesworth, R.M.; Sinclair, N.; Matthews, M.J.; Jeffries, J.R.; Willey, T.M.; et al. Ultrafast Dynamics of Laser-Metal Interactions in Additive Manufacturing Alloys Captured by in Situ X-ray Imaging. *Mater. Today Adv.* **2019**, *1*, 100002. [CrossRef]

52. Mukherjee, T.; Zuback, J.S.; De, A.; DebRoy, T. Printability of Alloys for Additive Manufacturing. *Sci. Rep.* **2016**, *6*, 19717. [CrossRef]

53. He, X.; DebRoy, T.; Fuerschbach, P.W. Alloying Element Vaporization during Laser Spot Welding of Stainless Steel. *J. Phys. D Appl. Phys.* **2003**, *36*, 3079–3088. [CrossRef]

54. Alcock, C.B.; Itkin, V.P.; Horrigan, M.K. Vapour Pressure Equations for the Metallic Elements: 298–2500K. *Can. Metall.* **Q.** **1984**, *23*, 309–313. [CrossRef]

55. Haynes, W.M. *CRC Handbook of Chemistry and Physics*; CRC Press: Boca Raton, FL, USA, 2014; ISBN 978-0-429-17019-5.

56. Labusch, R. A Statistical Theory of Solid Solution Hardening. *Phys. Status Solidi (B)* **1970**, *41*, 659–669. [CrossRef]

57. Fleischer, R.L. Substitutional Solution Hardening. *Acta Metall.* **1963**, *11*, 203–209. [CrossRef]

58. Denton, A.R.; Ashcroft, N.W. *Vegard's Law*. *Phys. Rev. A* **1991**, *43*, 3161–3164. [CrossRef]

59. Ledbetter, H.M.; Reed, R.P. Elastic Properties of Metals and Alloys, I. Iron, Nickel, and Iron-Nickel Alloys. *J. Phys. Chem. Ref. Data* **1973**, *2*, 531–618. [CrossRef]

60. Zander, J.; Sandström, R.; Vitos, L. Modelling Mechanical Properties for Non-Hardenable Aluminium Alloys. *Comput. Mater. Sci.* **2007**, *41*, 86–95. [CrossRef]

61. Kruth, J.-P.; Deckers, J.; Yasa, E.; Wauthlé, R. Assessing and Comparing Influencing Factors of Residual Stresses in Selective Laser Melting Using a Novel Analysis Method. *Proc. Inst. Mech. Eng. Part B J. Eng. Manuf.* **2012**, *226*, 980–991. [CrossRef]
62. Levkulich, N.C.; Semiatin, S.L.; Gockel, J.E.; Middendorf, J.R.; DeWald, A.T.; Klingbeil, N.W. The Effect of Process Parameters on Residual Stress Evolution and Distortion in the Laser Powder Bed Fusion of Ti-6Al-4V. *Addit. Manuf.* **2019**, *28*, 475-484. [CrossRef]

63. Lu, T.J.; Fleck, N.A. The Thermal Shock Resistance of Solids. *Acta Mater.* **1998**, *46*, 4755-4768. [CrossRef]

64. Toda-Caraballo, I.; Galindo-Nava, E.I.; Rivera-Díaz-del-Castillo, P.E.J. Unravelling the Materials Genome: Symmetry Relationships in Alloy Properties. *J. Alloys Compd.* **2013**, *566*, 217-228. [CrossRef]

65. *Ansys GRANTA EduPack Software*, ANSYS, Inc.: Cambridge, UK, 2019.

66. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A Fast and Elitist Multiobjective Genetic Algorithm: NSGA-II. *IEEE Trans. Evol. Computat.* **2002**, *6*, 182-197. [CrossRef]

67. Blank, J.; Deb, K. Pymoo: Multi-Objective Optimization in Python. *IEEE Access* **2020**, *8*, 89497-89509. [CrossRef]

68. Brans, J.P.; Vincke, P. Note—A Preference Ranking Organisation Method: (The PROMETHEE Method for Multiple Criteria Decision-Making). *Manag. Sci.* **1985**, *31*, 647-656. [CrossRef]

69. Figueira, J.; Greco, S.; Ehrogott, M. *Multiple Criteria Decision Analysis: State of the Art* *Surveys*; International Series in Operations Research & Management Science; Springer: New York, NY, USA, 2005; Volume 78, ISBN 978-0-387-23067-2.

70. Yadroitsev, I.; Smurov, I. Selective Laser Melting Technology: From the Single Laser Melted Track Stability to 3D Parts of Complex Shape. *Phys. Procedia* **2010**, *5*, 551-560. [CrossRef]

71. Yusuf, S.; Chen, Y.; Boardman, R.; Yang, S.; Gao, N. Investigation on Porosity and Microhardness of 316L Stainless Steel Fabricated by Selective Laser Melting. *Metals* **2017**, *7*, 64. [CrossRef]

72. Sun, Z.; Tan, X.; Tor, S.B.; Chua, C.K. Simultaneously Enhanced Strength and Ductility for 3D-Printed Stainless Steel 316L by Selective Laser Melting. *NPG Asia Mater.* **2018**, *10*, 127-136. [CrossRef]

73. Huang, G.; Wei, K.; Deng, J.; Liu, M.; Zeng, X. High-Power Laser Powder Bed Fusion of 316L Stainless Steel: Defects, Microstructure, and Mechanical Properties. *J. Manuf. Process.* **2022**, *83*, 235-245. [CrossRef]

74. Amato, K.N.; Gaytan, S.M.; Murr, L.E.; Martinez, E.; Shindo, P.W.; Hernandez, J.; Collins, S.; Medina, F. Microstructures and Mechanical Behavior of Inconel 718 Fabricated by Selective Laser Melting. *Acta Mater.* **2012**, *60*, 2229-2239. [CrossRef]

75. Ma, M.; Wang, Z.; Gao, M.; Zeng, X. Layer Thickness Dependence of Performance in High-Power Selective Laser Melting of 1Cr18Ni9Ti Stainless Steel. *J. Mater. Process. Technol.* **2015**, *215*, 142-150. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.