Research on Cracked Conditions in Nickel Chrome Alloy Ni50Cr33W4.5Mo2.8TiAlNb, Obtained by Direct Laser Deposition

: Nowadays, additive manufacturing (AM) is a powerful way to make complex-shaped components for airspace engineering from nickel-based superalloys. So, while nickel-based superalloys could easily be we L-DED in sheet-metal thicknesses, they suffered from strain-age cracking and solidiﬁcation during AM or in the post-weld aging treatment. This is attributed to the fact that besides the limitation of γ (cid:48) - phase forming elements (Al and Ti), as they form by AM very rapidly and reduce ductility, the majority of the superalloys contain carbide-forming elements such as Cr, Mo, and W. The precipitation of carbides, which is very effective in strengthening, develops cracks in the heat-affected zone (HAZ) during AM. The difference in isochoric heat capacities and the thermal expansion coef-ﬁcient (TEC) at the phase boundary leads to the appearance of dangerous local destruction energy. If the area of the interfacial interface is sufﬁciently extended, then the accumulation of this energy reaches a level sufﬁcient for a crack formation. We have proposed a crack initiation criterion (CIC) for assessing the dangerous level of fracture energy. The CIC was derived from an estimate of the local energy balance from the heat transfer equation for the two-phase area. Practical approbation of the criterion was carried out after L-DED of samples from Ni50Cr33W4.5Mo2.8TiAlNb (EP648) alloy powder with an increased carbon content based on the study of the chemical composition near the crack formed during solidiﬁcation. Using the proposed criterion provides an opportunity to give the rank to carbide-forming elements according to the degree of their inﬂuence on the fracture energy. Thus, the release of aluminum carbide turned out to be 5.48 times more dangerous than the release of titanium carbide and more than 5 times more dangerous than the release of tungsten carbide and molybdenum.


Introduction
Additive manufacturing (AM) is considered to be a revolutionary technology [1], which allows engineers to create complex shaped components in an end-to-end workflow that quickly moves from CAD files to the final part [2][3][4]. Numerous unique advantages have enabled AM technology to demonstrate significant potential in complex metal fabrication, multi-material/multi-functional design, and integrated manufacturing [2,4] As a result, innovative design may be obtained using AM for nickel-based superalloys [5]. AM can provide a competitive advantage by enhancing heat transfer and reducing turbine Metals 2022, 12,1902 2 of 17 blade pressure loss in complex microchannel cooling systems. [6]. In this regard, electron beam-powder bed fusion (E-PBF) [7][8][9] and laser beam-powder bed fusion (L-PBF) [10,11] have shown the potential for the processing of nickel-based superalloys [12,13]. High-γ nickel-based superalloys are abundantly used in aircraft and power generation turbines, rocket engines, and other high-temperature environments because of their excellent corrosion and oxidation resistance and superior mechanical properties [1]. This phase based on Ni3(Ti, Al) that has an ordered FCC L12 structure [14] constitutes the precipitate used to strengthen the alloy. The γ phase is coherent with the matrix of the superalloy having a lattice parameter that varies by around 0.5%. As particles of γ precipitates aggregate, they decrease their energy states by aligning along the <100> directions, forming cuboidal structures [15].
The physical phenomena occurring in the melting pool with AT are largely similar to the phenomena occurring in welding. Unfortunately, cracks were frequently observed in AM nickel-based superalloys [16][17][18][19]. So, while nickel-based superalloys could easily be we L-DED in sheet-metal thicknesses, they suffered from strain-age cracking during welding or in the post-weld aging treatment. These cracks were found in the heat-affected zone (HAZ) of the weld. This is attributed to the fact that the majority of the superalloys contain 3 to 5% (Al + Ti) and are strengthened γ phase. This precipitation is very effective in strengthening cracks developed in HAZ during welding as it forms very rapidly and reduces ductility [20].
The overall complexity associated with processing advanced hot-section high-γ nickelbased superalloys through additive manufacturing (AM) originates from non-weldability due to the presence of the secondary -γ precipitate phases in high proportions (i.e., more than 60 vol. %) in these alloys, arising from combined additions of aluminum (Al) and titanium (Ti) [21,22]. While these additions provide excellent benefits in improving overall high-temperature creep properties in new castings, they cause different types of cracking such as solidification, grain boundary liquation, strain-age, and hot tearing [1].
The presence of high-γ fraction (which is an obvious peril of Ti + Al addition) increases the likelihood of cracking during post-weld heat treatment known as "strain-age cracking" [23]. The alloys also become susceptible to ductility-dip cracking associated with the formation of grain boundary carbides [24]. Based upon the compositions of 1.0% Al + 0.84% Ti (in wt. %) versus 0.28% Cr + 0.043% Co (in wt. %), the relationship between the alloy chemistry and the weld susceptibility is illustrated in Figure 1 [25]. [2,4] As a result, innovative design may be obtained using AM for nickel-based supera loys [5]. AM can provide a competitive advantage by enhancing heat transfer and r ducing turbine blade pressure loss in complex microchannel cooling systems. [6]. In th regard, electron beam-powder bed fusion (E-PBF) [7][8][9] and laser beam-powder bed fu sion (L-PBF) [10,11] have shown the potential for the processing of nickel-based supe alloys [12,13]. High-γ′ nickel-based superalloys are abundantly used in aircraft an power generation turbines, rocket engines, and other high-temperature environmen because of their excellent corrosion and oxidation resistance and superior mechanic properties [1]. This phase based on Ni3(Ti, Al) that has an ordered FCC L12 structure [14 constitutes the precipitate used to strengthen the alloy. The γ′ phase is coherent with th matrix of the superalloy having a lattice parameter that varies by around 0.5%. As part cles of γ′ precipitates aggregate, they decrease their energy states by aligning along th <100> directions, forming cuboidal structures [15].
The physical phenomena occurring in the melting pool with AT are largely simila to the phenomena occurring in welding. Unfortunately, cracks were frequently observe in AM nickel-based superalloys [16][17][18][19]. So, while nickel-based superalloys could easi be we L-DED in sheet-metal thicknesses, they suffered from strain-age cracking durin welding or in the post-weld aging treatment. These cracks were found in th heat-affected zone (HAZ) of the weld. This is attributed to the fact that the majority of th superalloys contain 3 to 5% (Al + Ti) and are strengthened γ′ phase. This precipitation very effective in strengthening cracks developed in HAZ during welding as it forms ver rapidly and reduces ductility [20].
The overall complexity associated with processing advanced hot-section highnickel-based superalloys through additive manufacturing (AM) originates from non-weldability due to the presence of the secondary -γ′ precipitate phases in high pro portions (i.e., more than 60 vol. %) in these alloys, arising from combined additions o aluminum (Al) and titanium (Ti) [21,22]. While these additions provide excellent benefi in improving overall high-temperature creep properties in new castings, they cause di ferent types of cracking such as solidification, grain boundary liquation, strain-age, an hot tearing [1].
The presence of high-γ′ fraction (which is an obvious peril of Ti + Al addition) in creases the likelihood of cracking during post-weld heat treatment known as "strain-ag cracking" [23]. The alloys also become susceptible to ductility-dip cracking associate with the formation of grain boundary carbides [24]. Based upon the compositions o 1.0% Al + 0.84% Ti (in wt. %) versus 0.28% Cr + 0.043% Co (in wt. %), the relationshi between the alloy chemistry and the weld susceptibility is illustrated in Figure 1 [25].   Researchers have proposed several potential crack initiation mechanisms in nickelbased superalloys: (i) The first mechanism is solid-state cracking. Harrison et al. suggested that residual thermal stresses would lead to solid-state cracking. In this case, the cracks usually showed the characteristics of clean and straight fracture surfaces. To improve the thermal shock resistance of nickel alloys in AM, Harrison et al. proposed a slight increase in the concentration of substituent atoms that strengthen the crystal lattice of the solid solution, which increases the ultimate tensile strength and yield strength at elevated temperature, thereby suppressing the tendency to form cracks [26]. (ii) The second mechanism is liquation cracking. Liquation cracks also occur as a result of the separation of liquid films. Liquid films are formed as a result of the remelting of low-melting precipitates or the eutectic phase [27,28]. This mechanism of crack formation differs from solid-state cracking [29]. Since the introduction of the main elements forming the γ -phase (Al + Ti) is limited by the weldability area, the heat resistance of superalloys for AM is ensured by the introduction of refractory carbideforming elements, in particular Cr, Mo, and W, which form MC, M 23 C 6 , and M 6 C (M = metal) [30]. A large amount of carbides provides dispersion strengthening and grain boundary stabilization. Han et al. [19] hypothesized that the Mo and Cr-rich carbides in Hastelloy X would dissolve and then form a liquid eutectic film with a low melting point. Liquation cracks are also commonly observed in the HAZ (Heat Affected Zone) of a Hastelloy X weld in earlier papers [31][32][33]. More recently, Zhang et al. [34] reported that the continuous formation of γ-M 23 C 6 eutectic at grain boundaries is responsible for liquation cracks in Hastelloy X by directed energy deposition (DED). In addition, they also hypothesized that solidification cracks could be produced from liquation cracks when liquid phase replenishment during solidification is insufficient. (iii) The third mechanism is solidification cracking. Cracks are formed at the final stage of solidification (fS is close to 1, and fS is the fraction of solid matter), when replenishment from the liquid phase becomes difficult [35][36][37][38][39][40], which is the result of the stretching of liquid films in the mushy zone [41]. Montero-Sistiaga et al. [42] further testified to solidification cracking in Hastelloy X, where the surface of the cracks had the appearance of a smooth surface of dendrites. Recent studies on solidification cracking have shown that strain rate, rather than stress or strain, plays a more important role in cracking, and some modern cracking models are based on this hypothesis [38,[43][44][45].
According to the cracking criterion proposed by Kou [38], cracking can start when the separation rate of grains dε local /dt exceeds the sum of the grain growth rate relative to each other (lateral growth rate) and the feed rate from the liquid phase along grain boundaries (GB).
As for solidification cracking, many theories were developed for cracking during solidification. Lippold [29] reviewed the following five theories of hot tearing or solidification cracking: the shrinkage-brittleness theory of Pumphrey and Jennings [46], the strain theory of Apblett and Pellini [47], the generalized theory of Borland [48], the modified generalized theory of Matsuda et al. [49], the technological strength theory of Prokhorov [50], and the vulnerable range theory of Clyne and Davies [51].
The theory of Rappaz, Drezet, and Gremaud (the RDG model) [44] has been the most widely cited theory of solidification cracking since it was developed in 1999, although it was not considered by Lippold [29]. The RDG model considers stress in the transverse direction of the mushy zone and replenishment from the liquid phase through the mushy zone, from the ends of columnar dendrites to their roots. However, the model does not take into account the grain boundary where solidification cracking occurs. Instead, the entire mushy zone is viewed as a single area of a semi-solid structure through which interdendritic fluid enters the dendrite roots. More recently, Kou [35] proposed a criterion for solidification cracking. In contrast to the RDG model [44], Kou [35] took into account the influence of three factors along the grain boundary between adjacent closely spaced columnar dendrites in the mushy zone: (1) Cracking inducing transverse tensile strain rate separating dendrites; (2) Opposing cracking counter speed of lateral growth of dendrites, facilitating their connection; (3) The rate of fluid replenishment to prevent the formation of voids (cracks). He proposed the following criterion for solidification cracking based on volume indicators: cracking can occur if the result of factor 1 exceeds the result of the sum of factors 2 and 3. Based on the criterion, Kou [35,52] further proposed the ratio 1 2 dT d( f S) as a simple indicator of solidification cracking susceptibility near fS = 1. He showed that the characteristic dendrite radius r of columnar dendrites is proportional to 1 2 f S, and thus the lateral growth rate dr dt is proportional to 1 2 d( f S) dt . For a given dT dt , are higher if the ratio dr dt is smaller for grains that grow and merge with each other to resist cracking. A slower increase in the characteristic radius-dr dt contributes to an increase in the length of the channels between the columnar dendrites, which slows down the replenishment from the liquid phase, which is necessary to avoid cracking. Maximum 1 2 dT d( f S) usually occurs near fS = 1, so a convenient index value is the maximum value 1 2 dT d( f S) at 1 2 f S = 0.99 [52], beyond which cracking is unlikely due to extensive interboundary bonding between grains. The RDG model [44] also assumed the absence of cracking beyond fS = 0.98, i.e., 1 2 f S = 0.99. The index has been found to work for Al alloys [53,54], Mg alloys [55], carbon steels [56], and Ni-based alloys [57]. Again, Scheul-based solidification curves can be used to relate this index to the slope solidification curve (Sheila): 1 Regardless of the mechanism of crack formation, the cause of its formation is the difference in chemical composition at points on both sides of the phase boundary. In this case, the difference in thermodynamic potentials between the points leads to crack formation if this difference exceeds the value characteristic of this material.
To quantitatively assess the risk of crack formation depending on the difference in the chemical composition of the phases in the components of superalloy powders, we will propose a method for determining and a criterion for assessing the dangerous level of fracture energy-the crack initiation criterion (CIC). The CIC will be derived from an estimate of the local energy balance from the heat transfer equation for the two-phase area. This criterion will be tested in the analysis of crack formation conditions during laser direct energy deposition (L-DED) of samples from Ni50Cr33W4.5Mo2.8TiAlNb alloy powder with an overestimated carbon content, which leads to excessive precipitation of carbides, which are the cause of crack formation. Based on the microchemical analysis of the γ-matrix in the area of cracking of the alloy samples obtained by L-DED and the composition of the precipitates of primary carbides, the critical values of CIC for characteristic precipitates will be determined.

Theoretical Foundations of the Research Method
In previous works [58,59], the authors of the article showed that the heat balance equation at the interface between two phases can serve as a basis for identifying the energy of discontinuity and identifying the conditions for the initiation of cracks, given the possible displacement of the phases relative to each other. Without taking into account heat transfers to the external environment, the heat balance equation near the phase boundary or an area with a large gradient in the change in chemical composition takes the following form: where f ρ is the phase fraction denoted by the index (w) near the phase boundary; w i , v i is the velocity of phases (w), (v) (points of a continuous medium with a certain chemical composition); c w Ω , c v Ω are the isochoric heat capacities of phases; c w σ , c v σ are the isobaric heat capacities of phases; α w m , α v m are the coefficients of linear expansion of phases; c Ω is the T is the material derivative of temperature; b i is the force constant at the interface.
Consider the condition for a crack formation. Let the speed of the medium in the local region corresponding to the state (1 − f ρ ) be: Then, in the region with the state f ρ the speed of the medium is defined as: where k i is the discontinuity in velocity corresponding to the condition of crack formation, and in the considered local area: In this case, f ρ can be interpreted as a Heaviside function defined by (5) on a section with coordinates Taking into account (3), (4), (5), the heat balance Equation (1) can be written as a differential Equation (6): where: In general, Equation (7) can be represented as: After simplifying (8), we obtain: Suppose that in (2) and (3), v oi = 0 and k i are small by definition, then analysis (9) shows that the function argument tg must be small because the function tg increases sharply with the growth of the argument. For small values of the argument tg(y) y→0 ≈ y. In this case, it follows from (9) The integration constant in (9) is found from the second of the boundary conditions (11), which follow from (2) and (3): From (11) we get: From the first boundary condition (11) follows the dependence for the discontinuity of the flow velocity (13): Consider the condition for the formation of a hypothetical crack with a size δ over time t on the interval (13): In the interval , the material derivative of the temperature change in the interfacial area will take the form: Considering only the rate of temperature change at a material point, we obtain an approximate estimate of the material derivative of temperature: Using (16), from dependence (14) we obtain the average value of the discontinuity interval: Let us assume that: In this case, from dependence (17), taking into account (7) and (18), it follows: where Λ can be interpreted as a condition of discontinuity equal to the width of a hypothetical crack divided by the length of the interval on which the velocity jump is recorded k i . The recalculation of the isochoric heat capacity to the isobaric heat capacity was carried out according to the Magnus-Lindemann equation [60]: where n is the number of atoms in the compound (n = l + m + k); T * m is the melting point of the compound with chemical formula A l B m D k .
Substituting the Magnus-Lindemann Equation (20) into dependence (19) and integrating over temperature, we obtain: For two arbitrarily close points located at a distance 2x 0i , the following condition will be satisfied: i.e., in this case, there will be no discontinuity (crack formation). The condition for maintaining the continuity of the medium Λ(0) > 0 corresponds to the middle of a small interval [−x i0 , x i0 ], x i0 → 0 , while f ρ → 1/2 according to (5).
In this case (21), taking into account (22), takes the form: The final temperature T 0 in (23) corresponds to the temperature of normal conditions, and the initial T m to the melting temperature of the phase that crystallizes last. Let it be a more fusible phase, denoted by the superscript v, i.e., T m = T v m . The second phase, denoted by the index w, is considered a more brittle separation from the first phase.
Taking into account the Magnus-Lindemann Equation (20), the dependence for the continuity condition (23) will take the form: Let us restrict continuity condition (24) from below, for which we introduce the following inequalities into (24): Inequality (25) is based on the fact that for the vast majority of alloys, the more brittle phase (w), which is separated from the phase (v) in the form of precipitation, has a lower coefficient of thermal expansion. An example is the precipitation of brittle carbides along the grain boundaries of the EP648 alloy (Table 1). For a complex substance with the structural formula A m A B m B D m D , where A, B, D are the components that make up the complex substance, and m A , m B , m D are, respectively, the number of atoms in the molecule of this substance, the isochoric heat capacity is determined according to the Neumann-Kopp rule [63]: From the Koref rule for the Debye law for isochoric heat capacity, it follows that the heat capacity of an elementary substance included in a complex structure with a higher Metals 2022, 12,1902 8 of 17 melting point will be greater than the heat capacity of the same substance included in a structure with a lower melting point.
Taking into account (28), it follows from (27) that: Substitution (29) into inequality (26) does not change the sign of inequality: In (30), the expression in square brackets is nothing else than the rule of mixtures for the isochoric heat capacity of the components. Passing to the atomic concentration of components (c at ) L at a point, we can obtain an estimate of the continuity preservation condition in a local neighborhood of this point: In (31), the coefficient 24.304 has the dimension J/K·mol. Dependence (31) characterizes the dimensionless specific internal energy, the excess of which can lead to crack formation with a local change in the molar isochoric heat capacity.
Cracks usually occur either along the phase interface (Figure 2a,b) or near the interface area in places of the greatest inhomogeneity of the chemical composition, for example, along grain boundaries (Figure 2c) [64].
Metals 2022, 12,1902 9 of 18 In (30), the expression in square brackets is nothing else than the rule of mixtures for the isochoric heat capacity of the components. Passing to the atomic concentration of components ( ) at L с at a point, we can obtain an estimate of the continuity preservation condition in a local neighborhood of this point: In (31), the coefficient 24.304 has the dimension J/K·mol. Dependence (31) characterizes the dimensionless specific internal energy, the excess of which can lead to crack formation with a local change in the molar isochoric heat capacity.
Cracks usually occur either along the phase interface (Figure 2a,b) or near the interface area in places of the greatest inhomogeneity of the chemical composition, for example, along grain boundaries (Figure 2с) [64].  . the interface between points (v) and (w) based on the dependences (30), (31) can be found by formula (32), where CIC is the crack initiation criteria (crack formation criterion).
Dependence (32) determines the potential difference of internal energy at points (w) and (v) at distance s, which can lead to the formation of a crack (Figure 3).

Materials and Methods
To determine the appearance of cracks in AM of parts from superalloys and to demonstrate the applicability of the proposed criterion, a study was conducted on L-DED with samples ( Figure 4) with length 40 mm, height 9 mm, and width a multiple of steps (7-8 steps) from Ni50Cr33W4.5Mo2.8TiAlNb (EP648) alloy on the installation manufactured by Institute of Laser and Welding Technologies (ILIST, St. Petersburg, Russia) based on 6-axis robot Fuji Automatic Numeric Control (FANUC, Oshino-mura, Japan) ( Figure 5). Fusing of nickel-based powder with the chemical composition shown in Table  2 was carried out on a pre-cleaned substrate-a plate of low carbon austenitic stainless steel AISI 316L 2 mm thick.

Materials and Methods
To determine the appearance of cracks in AM of parts from superalloys and to demonstrate the applicability of the proposed criterion, a study was conducted on L-DED with samples ( Figure 4) with length 40 mm, height 9 mm, and width a multiple of steps (7-8 steps) from Ni50Cr33W4.5Mo2.8TiAlNb (EP648) alloy on the installation manufactured by Institute of Laser and Welding Technologies (ILIST, St. Petersburg, Russia) based on 6-axis robot Fuji Automatic Numeric Control (FANUC, Oshino-mura, Japan) ( Figure 5). Fusing of nickel-based powder with the chemical composition shown in Table 2 was carried out on a pre-cleaned substrate-a plate of low carbon austenitic stainless steel AISI 316L 2 mm thick.

Materials and Methods
To determine the appearance of cracks in AM of parts from superalloys and to demonstrate the applicability of the proposed criterion, a study was conducted on L-DED with samples ( Figure 4) with length 40 mm, height 9 mm, and width a multiple of steps (7-8 steps) from Ni50Cr33W4.5Mo2.8TiAlNb (EP648) alloy on the installation manufactured by Institute of Laser and Welding Technologies (ILIST, St. Petersburg, Russia) based on 6-axis robot Fuji Automatic Numeric Control (FANUC, Oshino-mura, Japan) ( Figure 5). Fusing of nickel-based powder with the chemical composition shown in Table  2 was carried out on a pre-cleaned substrate-a plate of low carbon austenitic stainless steel AISI 316L 2 mm thick.    All welding tracks were deposited sequentially using technological scanning parameters (Table 3) with a pause of 40 s between each track. These parameters were chosen as the closest to industrial parameters. Two layers of the base material were previously fused onto the substrate to ensure that the substrate material at L-DED does not introduce an error into the results of the experiment. These layers were cut off together with the substrate after the samples were produced.
The obtained samples were examined for the presence of cracks. Analysis of the microstructure was carried out on thin sections of cross-section using a Zeiss Axio Vert (Göttingen, Germany) optical microscope A1 Mat with 200× and 500× magnification for each sample. To improve the visibility of the cracks, a gray filter was used in a bright field. The etching of the samples was carried out in a solution of the following acids: H2SO4-HCl-HNO3-HF in a proportion of 180-180-120-30 mL, respectively, by immersion for 1 h.
The processing of the obtained images of the microstructure was carried out in a specialized software system SIAMS 800 (Yekaterinburg, Russia). Along the boundaries of the identified cracks, a study of the chemical composition was carried out. The elemental composition of the interface was determined by energy dispersive analysis.
To determine the chemical composition, a Phenom ProX (Thermo Scientific Inc., Waltham, MA, USA) electron microscope with an EDS detector for energy dispersive analysis was used. The analysis was carried out on samples before etching with an accelerating voltage of 15 kV. The chemical composition was measured at the boundary of two materials with a step of 18-20 μm at a magnification of 1000×.

Results and Discussion
Data on the measurement of the chemical composition for the key components of the EP648 samples obtained by the L-DED method at points 1-12 ( Figure 6) are shown in Table 4. Typical data of trace element EDS analysis with SEM are presented in Figure 7.  All welding tracks were deposited sequentially using technological scanning parameters (Table 3) with a pause of 40 s between each track. These parameters were chosen as the closest to industrial parameters. Two layers of the base material were previously fused onto the substrate to ensure that the substrate material at L-DED does not introduce an error into the results of the experiment. These layers were cut off together with the substrate after the samples were produced.
The obtained samples were examined for the presence of cracks. Analysis of the microstructure was carried out on thin sections of cross-section using a Zeiss Axio Vert (Göttingen, Germany) optical microscope A1 Mat with 200× and 500× magnification for each sample. To improve the visibility of the cracks, a gray filter was used in a bright field. The etching of the samples was carried out in a solution of the following acids: H 2 SO 4 -HCl-HNO 3 -HF in a proportion of 180-180-120-30 mL, respectively, by immersion for 1 h.
The processing of the obtained images of the microstructure was carried out in a specialized software system SIAMS 800 (Yekaterinburg, Russia). Along the boundaries of the identified cracks, a study of the chemical composition was carried out. The elemental composition of the interface was determined by energy dispersive analysis.
To determine the chemical composition, a Phenom ProX (Thermo Scientific Inc., Waltham, MA, USA) electron microscope with an EDS detector for energy dispersive analysis was used. The analysis was carried out on samples before etching with an accelerating voltage of 15 kV. The chemical composition was measured at the boundary of two materials with a step of 18-20 µm at a magnification of 1000×.

Results and Discussion
Data on the measurement of the chemical composition for the key components of the EP648 samples obtained by the L-DED method at points 1-12 ( Figure 6) are shown in Table 4. Typical data of trace element EDS analysis with SEM are presented in Figure 7.   The isochoric heat capacity of the key elements was calculated by the Debye method using the dependence: where θ is the Debye's temperature, which is defined as     The isochoric heat capacity of the key elements was calculated by the Debye method using the dependence: where θ is the Debye's temperature, which is defined as The isochoric heat capacity of the key elements was calculated by the Debye method using the dependence: where θ is the Debye's temperature, which is defined as In Equations (33) and (34), is Planck's constant; k is Boltzmann's constant; ν is the vibration frequency of atoms; x is the parameter, determined on the basis of the solid-state theory [65]; and T is the temperature (all calculations are made for room temperature, T = 298 K). The Debye model has been quite successful in approximating the thermal energy properties of a variety of solids.
The characteristic Debye temperatures of chemical elements are known from the reference literature, for example, [60]. So, the isochoric heat capacity can be calculated by the method given in [66]. Table 5 contains data on the isochoric heat capacity of the key chemical elements of the EP648 samples obtained by the L-DED method. The calculated isochoric heat capacity at the points of control of the chemical composition is shown in Table 6. The calculation was carried out according to the Neumann-Kopp rule. The results of the continuity preservation condition (CPC) calculation for these points according to dependence (31) are also presented in Table 6. Table 6. Thermodynamic properties and CIC calculations at media points ( Figure 6). Analysis of the data in Table 6 shows that despite the significant difference in chemical composition (Table 4), the values of continuity conditions (31) for points 1-12 ( Figure 4) for the γ-phase matrix on both sides of the crack are almost the same. This suggests that the cause of the crack initiation is the precipitation of carbides, which have a thermodynamic potential different from the background of the rest of the medium and, consequently, the values of the condition for maintaining continuity. Before turning to dependence (32) to quantitatively determine the thermodynamic conditions for crack formation, let us dwell on the analysis of the conditions for carbide precipitation.

Isochoric Heat Capacity Cv at Points 1-12 (J·mol
The complication of alloying leads to a change in the phase composition and microstructure of nickel alloys: the appearance of excess phases during crystallization, which are eutectic (γ-γ ), and M 6 C carbides based on (Ni,Co) 3 (W,Cr) 3 n of alloying solution tungsten [67][68][69]. According to the data of [70], in order to decrease the negative effect on cracking, the elements are arranged in the sequence Al, Ti, W, Mo, Cr. In a similar sequence, these elements are arranged in decreasing order of their γ -forming ability: Al > Ti > Nb > W > Mo > Cr [71].
It is also clear that fast cooling (quenching) of Ni-Cr-Mo based alloys with Ni content more than 50% (in wt. %) from high temperature (e.g., 1250 • C) displays the γ phase with carbide precipitates. Only a long time low-temperature aging regime (e.g., 700-600 • C) leads to long-range ordering (LRO) and the formation of phases based on complex crystalline structures (such as σ, µ, or P phase) [72]. During additive manufacturing of Ni-Cr-Mo-based alloys cooling occurs quite quickly, so the dominant type of structure will be σ. Moreover, it is clear that fast cooling (quenching) of Ni-Cr-Mo-based alloys with Ni content of more than 50% (in wt. %) leads to precipitation of Hf, Ta, etc. carbides near the grain boundaries.
On slow cooling or under the influence of high MC temperature, the carbides are converted into M 6 C-type and M 23 C 6 -type carbides. At high temperature in In625 superalloy [73], the MC-type carbides degenerate to M 6 C-type and M 23 C 6 -type carbides. The precipitation of an M 6 C-type carbide phase within the matrix and the formation of near continuous films, comprising discrete M 6 C/M 23 C6 carbide particles, occur after prolonged (∼70,000 h) service at temperatures close to but less than 600 • C. The shape of MC-type carbides is an irregular lump and those of the M 6 C and M 23 C 6 are bone-like or long-striplike, and most are located near the grain boundaries. A long carbide band along the grain boundary is formed by M 6 C and M 23 C 6 .
Let us consider the values of the crack formation criterion (32) in relation to the precipitation of some Al, Ti, W carbides alongside a solid solution of the γ-phase matrix with continuity preservation condition CPC = 0.0564 (Table 3). The calculation results are shown in Table 7.
As initially follows from (19), continuity preservation condition (CPC) and crack initiation criteria CIC are dimensionless quantities. These are interpreted through the ratio of the infinitely small width of a hypothetical crack to the length of the conditional segment connecting the points at the crack boundaries, in which the mutual change in the chemical composition, and hence the thermodynamic characteristics, is significant from the point of view of crack formation.
It can be seen from the data in Table 7 that aluminum carbide has the highest tendency to form cracks in the γ-phase matrix (CIC = 0.20096) compared to titanium and tungsten carbides at CIC = 0.03630 and CIC = 0.03960, respectively. In this regard, a limited amount of Al and Ti carbide-forming elements is introduced into the alloy. Thus, the main reason for the crack formation during DED is the high concentration of tungsten carbides WC and MoC, which precipitate during alloying mainly along the grain boundaries. It should be noted that WC with a melting point of 3141K has a higher crystallization range compared to MoC, the melting point of which is 2960K. A quantitative indicator of the risk of crack formation (appearance of discontinuities) is the CIC indicator proposed by the authors, which for the presented experimental data at the interface with carbides W turned out to be equal to 0.03960. In addition to CIC, an equally important factor is the concentration of carbide-forming elements (points 2, 4, 5, 6, 10 - Table 4, Figure 6), which, with an increased carbon content, leads to the formation of dangerous carbide precipitates.

Conclusions
Additive manufacturing imposes special requirements on the production of nickelbased superalloys due to their increased cracking susceptibility. This primarily refers to the restrictions on the number of Al and Ti elements forming the γ -phase. Secondly, the restrictions concern carbon, the increased content of which, together with the carbideforming elements W, Mo, etc., contribute to the intensive compact precipitation of primary carbides along the grain boundaries. Due to the heterogeneous concentration distribution of chemical elements, high-speed cooling, and the existence of metastable phases inherent in AM, this segregation is more intense than in traditional casting and can lead to cracking. The root cause of cracks is the difference in the chemical composition of the release materials and the γ-matrix. On the basis of the heat transfer equation at the boundary of a two-phase medium, the condition for the existence of areas with cracks in alloys at the interface was obtained. Such areas include, for example, precipitation of carbides (primary and secondary) along grain boundaries in nickel-based superalloys. The condition takes into account the difference in isochoric-isoentropic potentials, TEC, and crystallization temperatures of two phases. Practical approbation of the condition was carried out in the experimental study of the EP648 alloy with a high carbon content (C < 0.2% wt.) obtained by the L-DED method.
In our study, the next main idea of the method was realized. Based on the results of the chemical composition analysis at points on the boundary with cracks, the value of the continuity condition was determined, which was compared with the same indicator for precipitates of primary carbides of the MC type. The difference between these indicators is a characteristic of the accumulated fracture energy in the area of the phase interface. If the concentration of precipitates of MC-type carbides along any direction increases, and the total length of the phase interface reaches a certain critical level, then a crack is formed. Increased carbon content in the powder of the alloy under study led to the intense formation of carbides of the WC, MoC, and Al 4 C 3 types during the alloying of the powder. If, as was shown by calculations, the criterion for the formation of a CIC crack at the boundary of the γ matrix or γ/γ eutectic phase and MC carbide is significant (CIC Al 4 C 3 = 0.20096), then an excessive concentration of these carbides along the grain boundaries leads to the formation of cracks in the process solidification or strain-age cracking. As recommendations, we can offer: the use of a powder with a low carbon content (C < 0.05% wt.); reduced content of Al having high CIC; providing technological modes and conditions for AM with a high cooling rate to prevent excessive coagulation of carbides.
To comment on the usage of the obtained criteria, some limitations should be mentioned: (i) The obtained criteria in the form of complex characteristics summarize the significant differences in isochoric heat capacities and thermal expansion coefficients of phases in the local areas of precipitations. However, without additional studies near the crack area, it is impossible to accurately determine the dangerous criteria rates for the observed material and the conditions of its L-DED. (ii) Since the criteria are obtained as local characteristics, their application does not take into account the volume fraction and the sizes of precipitations.
Nevertheless, with the help of these criteria, even without additional studies, it is possible to provide a comparative analysis of the crack initiation danger of certain precipitations.