SinoTechIntel Academic Portal
Open AccessDOI: 10.1016/S1003-6326(26)67070-7Original Research

CFD modeling of gas−liquid mixing and splashing in industrial-scale bottom-blowing furnace: Effects of melt viscosity

School of Energy Science and Engineering, Central South University, Changsha 410083, China

Read Executive PreviewQuick FAQ
CFD modeling of gas−liquid mixing and splashing in industrial-scale bottom-blowing furnace: Effects of melt viscosity
Graphical Abstract / Figure
Published In
Transactions of Nonferrous Metals Society of China (中国有色金属学报)
Published:January 15, 2026Edition:Vol. 32, Issue 1 • pp. 100-112Citation:Cheng-lin LI 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

  • • • Optimal gas−liquid mixing is achieved at melt viscosities of 0.01 and 0.5 Pa·s, whereas minimal splashing occurs within 0.1−0.25 Pa·s (approximately 1076−1100 °C); this non-overlapping window forces operators to prioritize either mixing or splash suppression, directly affecting reaction completeness and metal recovery. • • Shear stress on the furnace lining increases with melt viscosity and rises particularly rapidly in the 0.1−0.25 Pa·s range, making this interval a critical erosion-risk zone that accelerates refractory wear and shortens campaign life. • • The regions surrounding the oxygen lances are the most susceptible to erosion, identifying localized lance-adjacent refractory as the primary failure point and directing maintenance and design efforts to this area. • • A recommended operating viscosity of 0.01−0.1 Pa·s minimizes erosion while maintaining acceptable mixing, providing a concrete control target for industrial bottom-blowing lead smelting that co-processes lead-containing waste materials.
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

The co-processing of lead-containing waste materials in oxygen bottom-blowing lead smelting alters melt viscosity, which in turn governs gas−liquid mixing efficiency, splashing behavior, and furnace lining erosion. A multi-fluid volume of fluid (VOF) model, incorporating experimentally determined high-lead slag viscosity values, was applied to an industrial-scale bottom-blowing furnace to resolve these coupled phenomena. The simulations demonstrate that optimal gas−liquid mixing occurs at viscosity values of 0.01 and 0.5 Pa·s, whereas minimal splashing is confined to the viscosity range of 0.1−0.25 Pa·s, corresponding to temperatures of approximately 1076−1100 °C. The regions surrounding the oxygen lances are identified as the most susceptible to erosion. Shear stress increases with melt viscosity, rising particularly rapidly within the 0.1−0.25 Pa·s interval. To minimize erosion, the melt viscosity should be maintained within 0.01−0.1 Pa·s. These findings provide quantitative thresholds for viscosity control in industrial bottom-blowing furnaces processing lead-containing secondary materials, enabling trade-offs between mixing intensity, splashing suppression, and lining protection to be optimized on an operational basis.

1. Introduction

Oxygen-enriched bottom-blown smelting is widely deployed in non-ferrous metallurgy for its high energy efficiency, intense smelting rate, and capacity to process complex low-grade ores. The introduction of lead-containing waste materials, such as slag and dust, into the feed modifies melt composition and raises viscosity, which can cause incomplete reactions, excessive splashing, and accelerated erosion of the furnace lining. Existing experimental studies have established that viscosity governs bubble size, momentum transfer, and flow regime transitions, but these findings have not been translated into quantitative, furnace-scale design thresholds for industrial bottom-blowing vessels.

This work addresses that gap by coupling experimentally determined high-lead slag viscosity values with a multi-fluid volume of fluid (VOF) model of an industrial-scale bottom-blowing furnace. The simulations resolve gas−liquid mixing efficiency, splashing behavior, and wall shear stress across a viscosity range relevant to lead smelting, thereby identifying the operating windows that balance mixing performance against splashing and lining erosion. The resulting viscosity thresholds provide a directly actionable basis for feed blending and thermal control in furnaces that co-process lead-bearing secondary materials.

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

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

Cite This Research Paper
Cheng-lin LI, You YAN, Ao HAO, Wei-wen HU, Zi-lin YANG, Dong-ling WU, Liu LIU, Hong-jie YAN (2026). CFD modeling of gas−liquid mixing and splashing in industrial-scale bottom-blowing furnace: Effects of melt viscosity. Transactions of Nonferrous Metals Society of China (中国有色金属学报). https://doi.org/10.1016/S1003-6326(26)67070-7
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 dominant failure mechanism of the furnace lining under high-viscosity melts, and where does it initiate?

Shear stress exerted by the gas−liquid flow on the refractory increases with melt viscosity, rising particularly rapidly in the 0.1−0.25 Pa·s range. The regions surrounding the oxygen lances experience the highest shear stress and are the most susceptible to erosion, making lance-adjacent refractory the primary failure location.

Can a single viscosity setpoint simultaneously deliver optimal mixing and minimal splashing?

No. Optimal gas−liquid mixing occurs at 0.01 and 0.5 Pa·s, while minimal splashing is observed only within 0.1−0.25 Pa·s (approximately 1076−1100 °C). These windows do not overlap, so operators must select a viscosity that prioritizes either mixing intensity or splash suppression, or accept a compromise outside both optima.

What viscosity range should be targeted to minimize erosion without sacrificing acceptable mixing?

The recommended operating viscosity is 0.01−0.1 Pa·s. This range avoids the rapid shear-stress escalation that occurs between 0.1 and 0.25 Pa·s while retaining mixing performance at the lower bound of 0.01 Pa·s.

How does the multi-fluid VOF model improve on prior experimental viscosity studies for industrial-scale design?

Prior experimental work identified qualitative trends, such as larger bubbles and enhanced momentum transfer with increasing viscosity, but did not provide furnace-scale thresholds. The multi-fluid VOF model integrates experimentally determined slag viscosity values and resolves local shear stress, splashing, and mixing in an industrial-scale bottom-blowing furnace, yielding quantitative operating windows rather than qualitative trends.

What are the operational implications of co-processing lead-containing waste materials on melt viscosity control?

Incorporating slag and dust into the smelting feed modifies material composition and raises melt viscosity, which can push the melt into the 0.1−0.25 Pa·s erosion-critical range. Maintaining viscosity within 0.01−0.1 Pa·s through feed blending and thermal management is therefore necessary to prevent accelerated lining wear while preserving reaction efficiency and metal recovery.

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