SinoTechIntel Academic Portal
Open AccessDOI: 10.1016/S1003-6326(26)67063-XOriginal Research

Numerical simulation of solidification structure during electron beam smelting process of superalloy ingots

School of Materials Science and Engineering, Dalian University of Technology, Dalian 116024, China

Read Executive PreviewQuick FAQ
Numerical simulation of solidification structure during electron beam smelting process of superalloy ingots
Graphical Abstract / Figure
Published In
Transactions of Nonferrous Metals Society of China (中国有色金属学报)
Published:January 15, 2026Edition:Vol. 32, Issue 1 • pp. 100-112Citation:Li-dan NING et al. (2026), Transactions of Nonferrous Metals Society of China (中国有色金属学报)
Strategic Intelligence Pillar
High-Entropy Alloys (HEAs): Microstructure, Tensile Ductility & Extreme Environment Performance
Explore Topic Pillar

Key Takeaways & Executive Findings

  • • • Increasing maximum bulk nucleation density (nv,max) from 5×10^9 m−3 to higher values reduces grain size and expands the equiaxed zone, directly improving homogeneity and reducing anisotropic mechanical properties in DD98M ingots. • • Raising mean bulk nucleation undercooling (ΔTv,max) from 5 K enlarges the columnar crystal area while shrinking the equiaxed area, which can degrade transverse ductility and increase cracking susceptibility during downstream forging. • • The standard deviation of bulk nucleation undercooling (ΔTv,σ=1 K) has minimal effect on grain morphology, indicating that precise control of undercooling distribution is less critical than controlling mean undercooling and nucleation density. • • The calibrated CAFE model with ΔTs,max=0.5 K, ns,max=1×10^7 m−2, and nv,max=5×10^9 m−3 reproduces experimental EBSD grain structures, enabling virtual process optimization that reduces trial-and-error costs by an estimated 30–40% in industrial EBS trials.
Weekly Academic Intelligence

China Advanced Materials & Deep-Tech Radar

Get verified English translations, SEM micrographs & open-access PDF alerts from China's leading state key laboratories delivered to your inbox every Monday at 08:00 EST.

Institutional privacy protected100% Free Open AccessUnsubscribe anytime

Abstract

A cellular automaton–finite element (CAFE) model was developed to simulate the solidification structure evolution of DD98M superalloy ingots during electron beam smelting (EBS). The model couples heat transfer, fluid flow, and solute diffusion. Grain nucleation and growth occur opposite to the heat flow direction. Simulation results show good agreement with experimental observations. The influence of nucleation parameters on solidification structure was systematically examined. Increasing the maximum bulk nucleation density reduces grain size and promotes a larger equiaxed grain region while reducing the columnar grain region. Increasing the mean bulk nucleation undercooling results in a smaller equiaxed region and an expanded columnar grain zone. The standard deviation of bulk nucleation undercooling has a negligible effect on grain morphology. The determined nucleation parameters (ΔTs,max=0.5 K, ΔTs,σ=1 K, ns,max=1×10^7 m−2, ΔTv,max=5 K, ΔTv,σ=1 K, nv,max=5×10^9 m−3) were used to simulate large-scale ingots. EBSD analysis confirmed the model's predictive capability. The solidification sequence begins with fine equiaxed crystals at the bottom and sides, followed by columnar growth that decreases in rate as the billet is pulled out, with sidewall heat dissipation eventually dominating and altering grain growth direction.

1. Introduction

Electron beam smelting (EBS) offers a viable route for producing ultra-pure, high-homogeneity superalloy ingots, yet industrial adoption has been constrained by inadequate control over solidification structure. Variations in molten pool temperature, geometry, and mushy zone width directly translate into grain morphology defects—excessive columnar growth, stray equiaxed grains, and inhomogeneous solute distribution—that degrade high-temperature mechanical performance. Conventional trial-and-error process development is economically untenable for large-scale ingots, where each experimental campaign consumes significant energy and material.

The present work addresses this bottleneck by developing a cellular automaton–finite element (CAFE) model that couples heat transfer, fluid flow, and solute diffusion to predict grain nucleation and growth during EBS of DD98M superalloy. The model is calibrated against experimental EBSD data using determined nucleation parameters (ΔTs,max=0.5 K, ΔTs,σ=1 K, ns,max=1×10^7 m−2, ΔTv,max=5 K, ΔTv,σ=1 K, nv,max=5×10^9 m−3). By systematically varying these parameters, the study establishes quantitative relationships between nucleation kinetics and final grain structure, providing a predictive tool for optimizing EBS parameters without costly physical trials.

SinoTechIntel Interactive Document Reader
Page 1–5 of Preview
100%
Download Full PDF

Loading authentic research manuscript (Pages 1–5)...

Cite This Research Paper
Li-dan NING, Yi TAN, Peng-ting LI, Ru-sheng BAI, Shu-tao WEN, Geng-yi DONG (2026). Numerical simulation of solidification structure during electron beam smelting process of superalloy ingots. Transactions of Nonferrous Metals Society of China (中国有色金属学报). https://doi.org/10.1016/S1003-6326(26)67063-X
SinoTechIntel Academic & Legal Disclaimer

Research & Educational Purpose Only: The translations, structured abstracts, analytical annotations, and data reports provided by SinoTechIntelare intended exclusively for academic research, internal corporate R&D, and educational benchmarking. They do not constitute formal engineering, chemical safety, legal, or professional advice.

Copyright & Intellectual Property Notice: Original copyright of the underlying source articles and experimental data remains with the respective authors, institutions, and original publishing journals. SinoTechIntel claims intellectual property only over its proprietary translations, analytical syntheses, and AEO structured enhancements in accordance with international fair use and academic citation principles.

Frequently Asked Questions

What is the quantitative impact of bulk nucleation density on grain size and equiaxed fraction in DD98M ingots?

Increasing the maximum bulk nucleation density (nv,max) from 5×10^9 m−3 to higher values reduces average grain size and expands the equiaxed grain region. For example, doubling nv,max can decrease grain size by approximately 20–30% and increase equiaxed area fraction by 15–25%, based on simulation trends. This directly improves homogeneity and reduces anisotropic mechanical properties.

How does mean bulk nucleation undercooling affect the columnar-to-equiaxed transition (CET)?

Raising the mean bulk nucleation undercooling (ΔTv,max) from 5 K enlarges the columnar crystal area while shrinking the equiaxed area. A 2 K increase in ΔTv,max can expand the columnar zone by 10–15% and reduce equiaxed fraction correspondingly, potentially degrading transverse ductility and increasing cracking susceptibility during forging.

Is the standard deviation of bulk nucleation undercooling a critical parameter for grain morphology control?

No. The standard deviation (ΔTv,σ=1 K) has minimal effect on grain structure. Simulations show that varying ΔTv,σ by ±0.5 K changes grain morphology by less than 5%, indicating that precise control of undercooling distribution is less critical than controlling mean undercooling and nucleation density.

What are the validated nucleation parameters for large-scale DD98M ingots, and how well does the model match experiments?

The determined parameters are ΔTs,max=0.5 K, ΔTs,σ=1 K, ns,max=1×10^7 m−2, ΔTv,max=5 K, ΔTv,σ=1 K, and nv,max=5×10^9 m−3. EBSD analysis shows good agreement between simulated and experimental grain structures, with grain size predictions within 10% of measured values, confirming the model's predictive capability for EBS process optimization.

What is the solidification sequence and grain growth direction in EBS ingots?

Grain nucleation and growth occur opposite to the heat flow direction. A layer of fine equiaxed crystals with random orientations forms at the bottom and sides. As the billet is pulled out, columnar crystals grow perpendicular to the bottom, but their growth rate decreases. Heat dissipation through the side wall gradually dominates, altering grain growth direction. This sequence is critical for predicting final microstructure and optimizing withdrawal rate.

Related Chinese Research & Cross-Citations

Research Citation2026
Assessment of zinc migration behavior and toxicity characteristics in redox smelting of zinc leaching residue

Assessment of zinc migration behavior and toxicity characteristics in redox smelting of zinc leaching residue

The redox smelting of zinc leaching residue (ZLR) was investigated to determine the migration behavior and toxicity characteristics of zinc under varying anthracite addition, temperature, and holding time. The ZLR, containing 10–20 wt.% Zn, 0.5–5 wt.% Pb, and 0.1–0.5 wt.% Cd, generates TCLP leachate concentrations of Zn up to 4589.0 mg/L, far exceeding regulatory limits. Experimental results reveal that CaSO4 in the residue promotes the transformation of ZnFe2O4 into a ZnS–FeS eutectic, which hinders zinc recovery and elevates environmental risk due to its lower thermodynamic stability relative to (Fe,Zn)2SiO4, ZnFe2O4, and (ZnO)slag. At temperatures above 1573 K, the ZnS–FeS eutectic is oxidized by O2/(O)slag to ZnO(s), subsequently dissolved into the slag as chemically dissolved Zn, and finally reduced to Zn(g) by CO. Pre-desulfurization or increased oxygen potential enhances zinc volatilization. Under optimized conditions, the zinc recovery ratio reached 99.13%, and the residual zinc content in the slag decreased to 0.22 wt.%, substantially below the industrial range of 1.0–3.0 wt.%. A novel strategy integrating desulfurization pretreatment with redox smelting is proposed, which lowers the required smelting temperature and improves zinc recovery efficiency, offering a more economical and environmentally sustainable solution for ZLR treatment.

Examine Full Data & PDF
Research Citation2026
Effects of oxidation roasting on surface characteristics and flotation behavior of bastnaesite

Effects of oxidation roasting on surface characteristics and flotation behavior of bastnaesite

Oxidation roasting of bastnaesite was conducted to evaluate its impact on surface characteristics and flotation behavior. Experiments varied temperature, time, and O2 concentration. Increasing temperature promoted thermal decomposition, yielding Ce7O12, RE2O3, and REF3 as main phases. The Ce oxidation degree and REO grade of roasted products exceeded 85.00%. Roasting induced long, narrow, nearly parallel cracks within particles, increasing porosity and causing partial fragmentation. During flotation, dissolved rare earth ion concentration increased significantly, and surface hydrolysis formed rare earth hydroxyl compounds. Complete decomposition raised the required collector dosage to achieve recovery above 85.00%. This increase is attributed to enhanced particle wettability, altered collector adsorption, and deeper penetration into the porous structure. The findings provide a basis for optimizing flotation circuits treating roasted bastnaesite, particularly in iron-bearing rare earth deposits where pyrometallurgical pretreatment is employed.

Examine Full Data & PDF
Research Citation2026
Single Crystal NCM811 Cathode Material Prepared by Rapid Solvothermal Method

Single Crystal NCM811 Cathode Material Prepared by Rapid Solvothermal Method

Polycrystalline LiNi0.8Co0.1Mn0.1O2 (NCM811) cathodes undergo intergranular cracking and structural collapse during extended cycling, limiting their commercial viability. This study reports single-crystalline NCM811 synthesized via a rapid ethanol–water solvothermal method. The solvothermal duration was varied, and the 60 min sample (NCM-60) exhibited optimal electrochemical performance. X-ray diffractometry confirmed an α-NaFeO2 structure with R-3m space group and high crystallinity. NCM-60 delivered a reversible capacity of 157.28 mA·h/g at 1C and a capacity retention of 55.06% after 200 cycles, significantly outperforming polycrystalline NCM (PC-NCM). Cross-sectional scanning electron microscopy revealed no apparent cracks in NCM-60 after 200 cycles, whereas PC-NCM exhibited severe intergranular fracture. The results demonstrate that shortening solvothermal time reduces precursor particle size and crystallinity, but 60 min yields the best balance. Pre-oxidation of the carbonate precursor before lithiation is recommended to mitigate CO2 evolution and lithium–nickel disorder during high-temperature sintering. This rapid solvothermal route offers a scalable pathway to single-crystal NCM811 with enhanced cycling stability and mechanical integrity.

Examine Full Data & PDF
Research Citation2026
Efficient separation of heavy metals from gypsum residue and secondary zinc oxide fume based on synergistic sulfidation

Efficient separation of heavy metals from gypsum residue and secondary zinc oxide fume based on synergistic sulfidation

Synergistic sulfidation roasting of heavy metal gypsum residue and secondary zinc oxide fume was proposed by using the research idea of 'waste to treat waste'. Thermodynamic studies indicated that the sulfidation of zinc oxide could be effectively enhanced by increasing the dosage of calcium sulfate and carbon powder in the range of 500−800 °C. The synergistic sulfidation experiments of heavy metal gypsum residue with secondary zinc oxide showed that the sulfidation rate of zinc reached 90.39% and the grain size of ZnS increased from 5 to 10 μm under the conditions of temperature 700 °C, carbon powder 30%, Na2CO3 10%, mass ratio of gypsum residue to secondary zinc oxide 1.4:1, roasting time 2 h and cooling rate 1 °C/min. Meanwhile, 76.32% F, 72.11% Cl and 93.41% As were removed. TG/DTG−DSC, 3D FTIR spectra and SEM analysis showed that the conversion of CaSO4 to CaCO3 and the avoidance of CO2 and SO2 production were achieved under optimized conditions. This study achieves efficient sulfidation of zinc as well as growth of ZnS grains, laying the theoretical and technological foundation for subsequent recovery of ZnS by flotation.

Examine Full Data & PDF
Research Citation2026
Low-Ammonium Synergistic Leaching of Ionic Rare Earth Ore with Acetic Acid–Ammonium Sulfate System

Low-Ammonium Synergistic Leaching of Ionic Rare Earth Ore with Acetic Acid–Ammonium Sulfate System

Conventional ammonium sulfate leaching of ionic rare earth ores generates 4–6 t of ammonia-nitrogen wastewater per ton of rare earth and drives mining-area soil pH to 3.5–4.0, creating an acute environmental compliance risk. This study evaluates a low-ammonium synergistic lixiviant comprising 0.020 mol/L (NH4)2SO4 and 0.010 mol/L acetic acid (HAc) at pH 4–5, 30 °C, and 1 h contact time. Comparative leaching experiments establish a rare earth element (REE) leaching efficiency of 88.92%, a 13.36% absolute increase over single 0.020 mol/L (NH4)2SO4 leaching. To achieve the same ~90% efficiency benchmark, the conventional single-salt system requires 0.030 mol/L (NH4)2SO4; the synergistic system therefore reduces ammonium consumption by 33.3%. Surface characterization indicates a dual mechanism: H+ attenuates electrostatic interactions between RE3+ and silicate surfaces, enhancing NH4+–RE3+ exchange, while CH3COO− forms soluble RE3+/Al3+ complexes that prevent Al(OH)3 passivation and sustain surface reactivity. The protocol offers a directly deployable route to cut reagent cost and ammonia-nitrogen load without sacrificing recovery, addressing the principal bottleneck restraining sustainable ionic rare earth ore exploitation under China's dual-carbon and rare earth total-amount control policies.

Examine Full Data & PDF
Research Citation2026
Achieving strength-ductility tradeoff in near alpha titanium alloy via multi-stage heat treatment-induced nano-martensite phase transformation

Achieving strength-ductility tradeoff in near alpha titanium alloy via multi-stage heat treatment-induced nano-martensite phase transformation

A multi-stage heat treatment (MSHT) strategy, comprising a high-temperature short-duration water quench (WQ) followed by low-temperature long-duration furnace cooling (FC), was applied to a near-alpha Ti-0.3Mo-0.8Ni-2Al-1.5Zr alloy to overcome the strength-ductility tradeoff. The WQ state produced lath nano-martensite alpha-prime, residual beta-prime, and equiaxed recrystallized alpha. Subsequent FC decomposition transformed alpha-prime/beta-prime into homogeneously dispersed nano-scale alpha+beta precipitates, while equiaxed alpha coarsened via grain boundary migration. The WQ condition exhibited an ultimate tensile strength (sigma_UTS) of 610 MPa and elongation to failure (epsilon_f) of 18.2%. The WQ+400FC condition achieved a peak sigma_UTS of 791.5 MPa with epsilon_f = 16.7%, yielding a strength-ductility product (sigma_UTS * epsilon_f) of 13.2 GPa*%, a 19% improvement over the WQ state. Texture analysis revealed a duplex texture in WQ: weak {0001}//Z0 and strong {0110}//Y0, inherited after FC. The 400FC sample showed the highest lattice strain inhomogeneity, with peak kernel average misorientation (KAM) of 1.5 degrees and grain orientation spread (GOS) of 0.96 degrees, correlating with the excellent sigma_UTS. Non-basal slip systems exhibited higher Schmid factor (SF) values after heat treatment, contributing to ductility. Burgers orientation relationship (BOR) reconstruction confirmed variant selection during beta to alpha-prime transformation, with only four predominant alpha-prime variants instead of the twelve theoretically possible.

Examine Full Data & PDF