Home Latest Articles
Latest Articles
  • Tao LIANG, Yue LI, Bin YANG, Jianglin LIU, Yaqiang LI, Meijie LIANG
    Journal of the Chinese Ceramic Society. 2026, 54(10): 3499-3517.
    Introduction

    Machine learning has been widely applied to various concrete materials tasks, such as macroscopic performance prediction and mix proportion optimization. The performance of machine learning models depends strongly on the quantity and quality of the available datasets. However, constructing large-scale and high-quality concrete mix proportion datasets remains challenging because concrete experiments are time-consuming, costly, and easily affected by material sources, curing conditions, and human operations. At present, small-sample modeling problems are mainly addressed through two strategies: model optimization and data augmentation. Model optimization methods mainly improve the extraction of information from existing small-sample datasets by adjusting model structures or training strategies. However, they do not fundamentally increase the amount of data, and their improvement in prediction accuracy is still limited by the original sample size. Conventional data augmentation methods mainly rely on generative adversarial networks, which learn the real data distribution through adversarial training between a generator and a discriminator, thereby generating virtual samples to expand the original dataset. Nevertheless, when dealing with complex tabular data, generative adversarial networks may encounter unstable training, insufficient sample diversity, and distribution shifts in generated data, resulting in limited consistency between virtual and real data. To address these issues, this study proposes MixDiff, a denoising diffusion probabilistic model (DDPM) for generating virtual concrete mix proportion data under small-sample conditions. This model aims to generate high-quality virtual concrete mix proportion data.

    Methods

    A small-sample dataset containing 100 sets of concrete mix proportions was first established. The input features included water, cement, fly ash, mineral fines, river sand, manufactured sand, coarse aggregate, water-reducing agent, total binder content, water-to-binder ratio, total fine aggregate content, and sand ratio. The target variables included 7 d compressive strength, 28 d compressive strength, slump, expansion degree, initial setting time, and final setting time. MixDiff was developed using a fully connected neural network and a denoising diffusion probabilistic model. The model consisted of a sinusoidal timestep embedding module and four stacked residual fully connected blocks. Each residual block contained fully connected layers, batch normalization layers, and SiLU activation functions. In the forward diffusion process, Gaussian noise was progressively added to the normalized data according to a cosine variance schedule. In the reverse generation process, a neural network was used to predict the noise and progressively reconstruct concrete mix proportion samples with realistic distribution characteristics from standard Gaussian noise. Subsequently, the quality of virtual data generated by MixDiff was quantitatively evaluated from two perspectives: marginal distribution and feature correlation. Under five-fold cross-validation and grid search, several mainstream machine learning models were further used to compare the prediction performance of the baseline mode and the augmentation mode. Finally, 15 mix proportions were selected from the generated virtual dataset for experimental validation to evaluate the engineering plausibility of the generated data.

    Results and discussion

    The quantitative evaluation results showed that the virtual dataset generated by MixDiff showed high similarity to the real dataset. The average Wasserstein distance, K-S statistic, and J-S divergence between the two datasets were 0.019 5, 0.095 2, and 0.120 7, respectively. The Frobenius distance between their Pearson correlation matrices was 0.332 7, indicating that MixDiff effectively preserved the marginal distribution characteristics of the real data and the correlation structure among variables. The machine learning comparison results showed that the prediction performance of all models improved significantly after the introduction of MixDiff-generated data. In the baseline mode, the average R2 of the seven models was 0.817 1, and increased to 0.953 7 in the data augmentation mode. Meanwhile, the average RMSE decreased from 14.20 to 7.39, and the average MAE decreased from 8.77 to 4.66. Correspondingly, R2 increased by 16.72%, while RMSE and MAE decreased by 47.96% and 46.86%, respectively. Among the machine learning models, the deep neural network achieved the best performance after data augmentation, with R2, RMSE, and MAE values of 0.959 4, 7.00, and 4.52, respectively. In addition, the standard deviations of the five-fold cross-validation results were significantly reduced in the data augmentation mode. This indicates that MixDiff-generated data not only improved the prediction accuracy of machine learning models but also enhanced model stability and generalization ability. Experimental validation of 15 virtual mix proportions showed that the data generated by MixDiff had high plausibility in six key performance indicators, including 7 d compressive strength, 28 d compressive strength, slump, expansion degree, initial setting time, and final setting time. The average MAPE of each performance indicator was lower than 5%, and the overall average MAPE was 3.8%.

    Conclusions

    This study proposed MixDiff, a denoising diffusion probabilistic model for concrete mix proportion data generation. The results show that the virtual data generated by MixDiff were highly consistent with the real data in terms of marginal distribution and feature correlation. The Frobenius distance between the Pearson correlation matrices was only 0.332 7. After introducing MixDiff-generated data, the prediction performance of multiple machine learning models, including RF, SVR, KNN, XGBoost, LightGBM, CatBoost, and DNN, was significantly improved, with R2 values of all models increasing above 0.94. Based on the average performance of the seven machine learning models for the six output performance indicators, R2 increased by 16.72%, while RMSE and MAE decreased by 47.96% and 46.86%, respectively. Among these models, DNN achieved the highest R2 and the lowest RMSE in the data augmentation mode. Experimental validation of 15 virtual mix proportions showed that the generated results were highly consistent with the measured results in terms of 7 d compressive strength, 28 d compressive strength, slump, expansion degree, initial setting time, and final setting time. The MAPE of each performance indicator was below 5%, and the overall average MAPE across the six performance indicators was 3.8%. These results indicate that MixDiff can learn the underlying relationships between concrete mix proportions and macroscopic properties, and it can serve as an effective method for concrete mix proportion data augmentation under small-sample conditions.

  • Xu WU, Suhang LI, Shuai NING, Lin WANG, Changjiao LI, Jiayuan YE, Bo YANG
    Journal of the Chinese Ceramic Society. 2026, 54(10): 3518-3529.
    Introduction

    The chemical composition of cement clinker, particularly the percentage of key oxides, is crucial for evaluating cement quality and is vital for engineering-accident liability tracing. However, traditional methods for chemical composition analysis rely heavily on unhydrated original cement samples, which are often missing, improperly stored, or ambiguously labeled in engineering management, resulting in an insufficient basis for subsequent identification. Once cement hydrates, active minerals undergo irreversible reactions, making it challenging to invert the original clinker composition directly from the hardened paste. In contrast, backscattered electron (BSE) images of the hardened cement paste contain sufficient physical information to distinguish the original chemical compositions of different cement clinkers, thus providing a potential pathway for post-event tracing. Although deep learning methods have been explored for microstructural image mapping to infer macroscopic properties, conventional regression models are typically "black-box" point predictors that cannot quantify prediction uncertainty. The conditional probability density estimation method can not only provide the most probable predicted value, but also output the complete estimated probability-density curve along with its confidence assessment, thereby enabling quantitative characterization of prediction uncertainty. Accordingly, this study proposes a conditional probability density estimation method based on a Deconvolutional Density Network (DDN). Utilizing BSE images of hardened paste together with encoded water-to-cement (W/C) ratios and curing ages as conditional inputs, this model directly predicts the full conditional probability density of the main clinker oxide contents, providing a new objective and quantitative pathway for quality liability tracing.

    Methods

    Cement clinker samples from 47 different plants were used to prepare cement paste specimens with W/C ratios ranging from 0.25 to 0.70 and curing ages from 1 d to 28 d. The hydration was terminated by immersion in anhydrous ethanol, and the samples were polished for subsequent scanning. BSE images were acquired via a JSM-7610F field-emission SEM at 200× magnification and an accelerating voltage of 15 kV. For each sample, two scans were performed, and nine non-overlapping regions were captured per scan, resulting in a dataset of 844 independent BSE images. The proposed model consists of a feature extraction layer based on a residual network (ResNet) with self-supervised learning and cross-layer skip connections, a conditional encoder, and deconvolutional layers. The ResNet extracts multiscale morphological features from local textures to global patterns. The conditional encoder embeds the W/C ratio and curing-age conditions into the image features. The DDN is then employed to output the conditional probability density distributions for the percentages of eight oxides (CaO, SiO2, Al2O3, Fe2O3, SO3, Na2O, MgO, and K2O). The deconvolutional layers map the fused features to a smooth probability distribution over 256 discretized intervals. To evaluate the model, accuracy (ACC) under different confidence levels (α), mean absolute error (MAE), and mean squared error (MSE) were adopted as evaluation metrics.

    Results and discussion

    The estimation results of the proposed model indicate that, at a confidence level of α = 90%, the model achieves an accuracy of over 90% in estimating the percentage contents of five oxides (SiO2, Al2O3, Fe2O3, SO3, and Na2O) on the test set. When the confidence level is relaxed to α = 75%, the estimation accuracy for all eight oxides exceeds 90%. This could be attributed to the DDN's capability of modeling complete probability density distributions instead of simple point estimates, effectively resolving the nonunique-solution problem stemming from the intrinsic information ambiguity of BSE imaging. In terms of MAE and MSE metrics, the model exhibits relatively higher absolute errors for components with higher contents, such as CaO. This is generally attributed to the nature of absolute error metrics, which are positively correlated with the oxide content. Conversely, low-content oxides like Na2O show smaller estimation errors, possibly because their larger fluctuation ranges among different samples allow the model to capture clearer mapping relationships. In addition, ablation experiments further reveal that removing the W/C ratio and age inputs leads to a decrease in estimation accuracy across all confidence levels, demonstrating the effectiveness of the conditional encoder in integrating critical hydration parameters. Compared with the deterministic ResNet regression model, the proposed DDN model exhibits significantly better stability of estimation performance under higher precision requirements. This result demonstrates that the DDN model effectively mitigates the inherent uncertainty in the mapping from BSE images to chemical compositions, a limitation typical of deterministic regression models that lack reliability quantification.

    Conclusions

    The main conclusions of this study are summarized as follows. The proposed probability density estimation method successfully maps BSE images of hardened cement paste to the oxide contents of clinker. Experimental results demonstrate that the model achieves an accuracy of over 90% for five oxides at a confidence level of 90%. The adoption of the Deconvolutional Density Network enables the model not only to output the most probable value of each oxide content but also to provide a complete probability density curve with corresponding confidence assessment, thereby overcoming the limitation of traditional regression models that lack quantifiable reliability for quality tracing applications. Moreover, the ResNet-based feature extraction with cross-layer skip connections effectively captures multiscale features from BSE images, while the integration of W/C ratio and curing age via the conditional encoder significantly enhances the model predictive capability. This study demonstrates that the probabilistic retrospective estimation based on BSE images and DDN offers a robust and objective means to quantitatively infer the original clinker composition from hardened paste. As such, it holds promise as a transformative technique for cement-quality arbitration in engineering-accident liability investigations.

  • Tong LI, Peiyang GUO, Changhao LIU, Huisu CHEN
    Journal of the Chinese Ceramic Society. 2026, 54(10): 3530-3551.

    Magnesia cements have attracted increasing attention as low-carbon and functional binders because of their diverse reaction pathways, rapid hardening characteristics, and potential applications in rapid repair, fire-resistant boards, waste immobilization, carbon sequestration, and specialized engineering environments. Typical magnesia cement systems include magnesium silicate hydrate cement (M-S-H), magnesium phosphate cement (MPC), magnesium oxysulfate cement (MOS), magnesium oxychloride cement (MOC), and reactive magnesia cement (RMC). These systems involve complex reactions among MgO, soluble salts, siliceous or aluminosilicate components, water, and CO2, leading to the formation of various hydration and carbonation products, such as M-S-H/M-(A)-S-H gels, 3Mg(OH)2·MgCl2·8H2O, 5Mg(OH)2·MgCl2·8H2O, 3Mg(OH)2·MgSO4·8H2O, 5Mg(OH)2·MgSO4·7H2O, and KMgPO4·6H2O. Because these products are often highly sensitive to mix proportion, temperature, pore solution chemistry, relative humidity, CO2 concentration, and aggressive ions, it remains difficult to fully elucidate their phase assemblage evolution and long-term durability using experimental characterization alone. Thermodynamic modeling provides an effective approach to predict stable phase assemblages, pore solution composition, saturation indices, and phase stability boundaries under defined boundary conditions and has thus become an important tool for understanding the hydration, modification, degradation, and carbonation mechanisms of magnesia cements.

    This review summarizes recent progress in thermodynamic modeling of phase assemblage evolution in magnesia cements. First, the basic principles of thermodynamic modeling are introduced, including the role of ion activity, saturation index, chemical equilibrium, Gibbs free energy minimization, and mass action law. Particular attention is paid to the selection of activity coefficient models, because magnesia cement systems, especially MOC, MOS, and salt-exposed systems, commonly contain highly concentrated electrolyte solutions. Second, this review collects and summarizes the thermodynamic databases currently available for magnesia cements and further discusses their main limitations. The reviewed databases cover M-S-H, MPC, MOS, MOC, and carbonation-related systems, including poorly crystalline M-S-H/M-(A)-S-H gels, magnesium phosphate hydrates, basic magnesium sulfate phases, magnesium oxychloride hydrate phases, and magnesium carbonate hydrates. These data provide the basis for thermodynamic modeling of hydration, modification, salt attack, leaching, and carbonation processes. Nevertheless, important gaps remain in the existing databases. For M-S-H-based systems, thermodynamic data are still insufficient to describe the continuous variation of Mg/Si and Al/Si ratios in M-S-H and M-(A)-S-H gels. For MPC, the thermodynamic data for intermediate phosphate hydrates and retarder-related complexes remain incomplete. For MOS and MOC, more reliable thermodynamic parameters are needed for magnesium sulfate/chloride phases, especially under variable temperature and high-salinity conditions. In addition, the incorporation of silica fume, fly ash, slag, metakaolin, and other industrial by-products introduces Si, Al, Ca, and Fe into magnesium-based systems, potentially leading to M-S-H, M-(A)-S-H, C-S-H, hydrotalcite-like phases, calcium phosphate phases, and other secondary products; however, thermodynamic data for these multicomponent reaction products are still insufficient. For carbonation thermodynamic modeling, further improvement is needed in the thermodynamic data of metastable magnesium carbonate hydrate phases, as well as in the ion interaction parameters for multicomponent Mg-Cl-SO4-Na-Ca-K-CO2-H2O systems. Therefore, future database development should move from the compilation of individual phase data toward internally consistent multicomponent databases that can better represent the complex reaction environments of magnesia cements.

    Third, the applications of thermodynamic modeling in different magnesia cement systems are reviewed. For M-S-H, modeling studies show that the Mg/Si ratio strongly controls the equilibrium among M-S-H gel, Mg(OH)2, and SiO2. Low Mg/Si systems tend to form silica-rich M-S-H with better chemical stability, whereas high Mg/Si systems are more likely to generate Mg(OH)2, which may reduce resistance to salt attack and carbonation. The introduction of metakaolin, sodium carbonate, or other aluminosilicate components can further induce M-(A)-S-H gel, hydrotalcite-like phases, Al(OH)3, and zeolite-type phases. For MPC, thermodynamic modeling reveals that the Mg/PO4 ratio, pH, phosphate type, and temperature govern the transition from acidic magnesium phosphate phases to intermediate phases and ultimately to more stable struvite-type products. Modeling has also been applied to assess the effects of wollastonite, carbonation, chloride attack, sulfate attack, and seawater exposure on MPC phase assemblages. For MOS, modeling studies indicate that nMgO/nMgSO4 and nH2O/nMgSO4 ratios strongly affect the formation of 5-1-7, which is closely associated with compressive strength and water resistance. Chemical modifiers, especially citric acid and tartrate species, can alter Mg2+ activity and promote the formation of 5-1-7. For MOC, thermodynamic modeling has been widely used to clarify the stability of P3, P5, and Mg(OH)2 as functions of nMgO/nMgCl2, nH2O/nMgCl2, temperature, modifiers, carbonation, chloride attack, and sulfate attack. These studies explain why P5 can provide high early strength but tends to decompose into Mg(OH)2 under leaching or dilution conditions. Finally, carbonation modeling of RMC and other magnesia cements is discussed. The phase composition of magnesium carbonate hydrates is highly sensitive to CO2 concentration, humidity, temperature, and reaction time. Because full-equilibrium calculations tend to overpredict thermodynamically stable phases such as MgCO3 or Mg5(CO3)4(OH)2·4H2O, metastable equilibrium modeling is often performed by suppressing selected stable phases to better reproduce experimentally observed intermediate products. This approach has been applied to RMC, M-S-H, MPC, and MOC carbonation, revealing the importance of kinetic limitations and environmental boundary conditions.

    Overall, thermodynamic modeling has provided important insights into the hydration, modification, salt attack, leaching, carbonation, and durability of magnesia cements. Nevertheless, several limitations remain, including incomplete thermodynamic data for key hydrates and carbonates, insufficient ion interaction parameters for high-salinity multicomponent systems, simplified descriptions of low-crystalline gels and metastable phases, and limited coupling with kinetic processes. Future studies should further improve thermodynamic databases, validate predictions using XRD, TG, NMR, FTIR, SEM/EDS, ICP-OES, and pore solution analysis, and integrate thermodynamic modeling with kinetic, transport, humidity-controlled, and performance models. Such advances will enhance the predictive capability of thermodynamic modeling and support the rational design, durability assessment, and low-carbon applications of magnesia cements.

    Summary and Prospects

    Thermodynamic modeling has become a useful tool for understanding hydration, phase evolution, modification mechanisms, and the durability of magnesia cements. It can predict stable phases, pore solution chemistry, saturation indices, and phase stability boundaries under defined conditions, providing theoretical support for mix proportion design and performance optimization. Existing thermodynamic modeling studies have investigated the effects of raw-material mix proportions, temperature, modifiers, carbonation, and aggressive ions on the phase assemblages of M-S-H, MPC, MOS, MOC, and RMC systems, and considerable progress has been made in understanding the formation and transformation of M-S-H/M-(A)-S-H gels, magnesium phosphate hydrates, basic magnesium sulfate and chloride hydrate phases, and magnesium carbonate hydrates. However, several limitations remain. Thermodynamic data for M-S-H/M-(A)-S-H gels, magnesium carbonate hydrates, metastable intermediate products, and the temperature-dependent parameters of some hydration products are still incomplete. In addition, high-ionic-strength pore solutions require appropriate activity coefficient models, but Pitzer parameters for multicomponent systems remain insufficient. Moreover, thermodynamic modeling mainly reflects equilibrium or quasi-equilibrium states and cannot directly describe MgO dissolution, nucleation, growth, transport, or early hydration kinetics. Future work should focus on improving thermodynamic databases, refining activity models for high-salinity systems, and coupling thermodynamic modeling with experimental characterization, kinetic modeling, and transport simulations. This will enhance its reliability in guiding material design, durability assessment, and low-carbon applications.

  • Taiyu DUAN, Xin ZHANG, Yao JI, Weichao WANG, Qinyuan ZHANG
    Journal of the Chinese Ceramic Society. 2026, 54(10): 3554-3563.
    Introduction

    Fluorophosphate (FP) laser glass offers high rareearth solubility, low nonlinear refractive index, and a large Abbe number, making it a promising gain medium for high-power fiber lasers. The determination of the glass-forming region is fundamental to developing new laser glass systems. Previous studies have relied heavily on trialanderror approaches, yielding only limited data at prohibitive cost in terms of manpower, finances, and material resources. In this work, thermodynamic calculations are employed to investigate the glass-forming region of multicomponent FP glass systems. The results reveal that thermodynamic calculations can effectively predict the glass-forming region of FP glasses. Moreover, differences in cation field strength and interactions with anions are identified as critical factors governing the glassforming ability and structural stability of FP glasses. Based on these theoretical predictions, a ytterbiumdoped, high-gain FP laser glass and fiber are designed and fabricated through compositional optimization, and initial 1.0 μm laser output is successfully demonstrated.

    Methods

    The thermodynamic calculations are based on the free energy balance between liquid and solid phases in binary systems. By solving the equilibrium equation, the melting temperature is obtained as a function of composition, and the eutectic point is determined by locating the minimum of this function. For ternary and higher-component systems, the binary eutectic points are treated as “quasi-compounds”, and the calculation is performed iteratively to construct the theoretical glass-forming region. The predicted regions were validated against experimental data reported in the literature. Raman spectroscopy combined with peak-fitting analysis was used to quantify Qn species, while density and refractive index measurements were carried out for physical property evaluation. Yb3+-doped FP glasses and fibers were then fabricated based on the predicted glass-forming region, and their laser performance was characterized.

    Results and discussion

    Quantitative comparison between predicted and experimental glass-forming regions shows overlap ratios of 74.5%, 100%, and approximately 100% for the three systems examined, confirming that the method effectively locates the main compositional range of glass formation. Raman analysis reveals a transition from Q2-dominated chain structures at high phosphate contents to Q0/Q1 units at low phosphate contents, with the Q2 fraction decreasing from 73.4% to 8.4% as phosphate content drops from 80%(mole fraction) to 40%. Cation field strength is identified as a critical factor: the higher field strength of Zn2+ (0.43) relative to Ba2+ (0.24) leads to stronger network bonds but greater sensitivity to composition variation, resulting in a narrower glass-forming region. Guided by the theoretical predictions, Yb3+-doped FP fibers were fabricated with a peak gain coefficient of 5.56 dB/cm at 1064 nm and stable dual-wavelength laser emission at 1023.3 nm and 1031.5 nm under 974 nm pumping.

    Conclusions

    This work demonstrates that thermodynamic calculations based on eutectic point prediction provide a valid pre-screening tool for locating glass-forming regions in multicomponent FP systems. Cation field strength governs the glass-forming window and structural stability. Guided by theoretical predictions, Yb3+-doped FP fibers were successfully fabricated with stable 1.0 μm laser output, validating the effectiveness of this combined theoretical-experimental approach for developing novel laser glass materials.

  • Naiyong PENG, Zhen XU, Xian ZHANG, Xiaobing MA, Dandan ZHANG, Zhi HONG, Xiao XIAO, Wenxiao ZHAO, Ping LU, Yinsheng XU
    Journal of the Chinese Ceramic Society. 2026, 54(10): 3740-3749.
    Introduction

    Insulating materials play a crucial role in high-voltage power transmission, protecting equipment and ensuring the reliability of the power grid by effectively isolating current and voltage. Conventional ceramic insulators (such as porcelain or aluminium oxide insulators) are widely used in high-voltage power transmission, but they have inherent limitations, including a short service life, high manufacturing costs and many surface defects, which are likely to cause electrical breakdown in the material. Furthermore, defects such as internal porosity in ceramics can lead to deterioration of their mechanical and electrical properties. In contrast, glass insulators, owing to their exceptional mechanical strength and excellent stability, have become one of the materials of choice for long-span transmission lines. Among these, sodium-calcium silicate (NCS) glass offers significant advantages, including high transparency, excellent physicochemical stability, ease of manufacture, and low cost, and thus holds great potential for application in the field of high-voltage power transmission and transformation equipment. However, the mechanical strength and electrical properties of NCS matrix glass are insufficient to meet the demands of modern power equipment for high electrical strength, high mechanical strength and low dielectric loss. Therefore, the glass composition was optimized, and the NCS matrix glass system was modified by doping with B2O3 and K2O. The evolution of the network structure, optical transparency, mechanical strength and electrical properties with varying doping ratios was investigated to develop a material suitable for high-performance, long-life, high-voltage insulators.

    Methods

    The matrix glass composition was an NCS system. Two sets of glass samples were prepared using the high-temperature melting-quenching method. In Group B, B2O3 was added at varying molar fractions (0-20%(mole fraction)) while keeping the other components constant. In Group N, using the mixed alkali effect, K2O was introduced to partially replace Na2O, while maintaining the total content of Na2O and K2O constant. Subsequently, high-purity raw materials (SiO2, CaCO3, Na2CO3, H3BO3, K2CO3) were thoroughly mixed, placed in an alumina crucible and melted for 2 h at 1400 ℃. The melt was cast into preheated copper molds and annealed at 500 ℃ for 2 h to relieve internal stresses. The glass was cut into samples measuring 10 mm × 10 mm × 1 mm and 20 mm × 20 mm × 1 mm, and the surfaces were ground flat using a metallographic grinding and polishing machine. The material structure was characterized using differential scanning calorimetry (DSC), X-ray diffraction (XRD), and Raman spectroscopy. The optical properties were determined by UV-vis spectrophotometry. Mechanical strength was measured using a three-point bending test. The electrical properties included volume resistivity (measured at room temperature using a resistivity meter), dielectric loss, and dielectric breakdown strength (under DC voltage).

    Results and discussion

    The addition of B2O3 significantly alters the glass network structure. When a small amount of B2O3 is introduced, its [BO3] triangular units tend to replace some of the [SiO4] tetrahedra of SiO2, forming an irregular boron-oxygen network structure. This structural transformation reduces the proportion of bridging oxygen atoms in the glass network, leading to a decrease in network connectivity. As the B2O3 content increases further, [BO4] tetrahedra gradually form, the number of bridging oxygen atoms increases, the glass network structure becomes denser, and the structural stability of the glass is enhanced. This structural change enhances the rigidity of the network. When 5%(mole fraction) B2O3 was added, the breakdown strength of the sample reached a maximum of approximately 31.93 kV/mm, with a flexural strength of 108.37 MPa. However, when the B2O3 content was excessive, a large number of [BO3]/[BO4] units appeared in the glass, leading to the breaking of the Si-O-Si chains and a decrease in network polymerization, which severely affected glass performance.

    Compared to Na+ ions, the larger radius of K+ ions means that introducing K2O increases steric hindrance within the glass network, resulting in a denser glass structure and thereby improving flexural strength. When 4% K2O was added, the sample exhibited a flexural strength of 110.55 MPa, and the dielectric loss reached its lowest value (0.0027). The results indicate that an appropriate amount of K2O doping for partial replacement of Na2O can improve the glass's flexural strength and stiffness; however, the introduction of excessive K2O may lead to an increase in defects within the glass network, thereby reducing the material's mechanical properties.

    Conclusions

    In summary, the incorporation of an appropriate amount of B2O3 into NCS glass, combined with the mixed alkali effect resulting from the substitution of Na2O by K2O, is an effective strategy for enhancing the mechanical strength and electrical properties of glass insulating materials. The results indicate that as the B2O3 content increases, the material's flexural strength exhibits a marked nonlinear change; when 5% B2O3 is added, both the flexural strength and electrical breakdown strength of the material are enhanced. The introduction of an appropriate amount of K2O (partially replacing Na2O) alters the polarizability of the glass via the mixed alkali effect, thereby enhancing the material's electrical breakdown strength while maintaining low dielectric loss. When the B2O3 and K2O contents were 5% and 4%, respectively, the prepared sample exhibited high mechanical strength (flexural strength > 110 MPa), excellent insulation properties (room-temperature volume resistivity > 1012 Ω·m, dielectric loss < 0.003 at 50 Hz) and high electrical strength (breakdown strength > 32 kV/mm). These findings provide a theoretical basis and material-design guidance for glass insulating materials used in high-performance electrical equipment.

  • Hui RONG, Han JIA, Dee LIU, Zhihua LIU
    Journal of the Chinese Ceramic Society. 2026, 54(10): 3795-3812.

    Marine biofouling poses a significant threat to the durability of coastal infrastructure, with the attachment of macroorganisms such as oysters inducing complex effects on concrete substrates. This review systematically examines the attachment process of oysters on marine concrete, the multifaceted factors governing this process, and the subsequent influence on material durability, distinguishing between corrosive and protective mechanisms. The analysis further evaluates current strategies for either inhibiting or facilitating oyster colonization to support targeted management of marine structures.

    The review begins by elucidating the biological attachment process of oysters, from larval settlement through metamorphosis, highlighting the critical role of substrate selection. Key determinants of attachment are systematically categorized into environmental conditions, substrate properties, chemical gradients, and biological mediators. Environmental conditions including thermal regime, hydrodynamic forces, and light availability establish broader colonization patterns. Substrate properties such as surface microtopography, material composition, roughness, and color, critically influence larval settlement selectivity and adhesion strength. Chemical gradients, particularly pH value and specific ions, modulate larval behavior and enzymatic activity. Biological mediators, including conspecific cues and interspecific interactions, govern the spatial distribution and density of oyster populations. These factors collectively determine the settlement density and adhesion strength of oysters on concrete surfaces.

    The core of this review addresses the dual impact of oyster attachment on concrete performance. Corrosion can occur, driven by bioacids in larval adhesives that react with cement hydrates to form expansive products, leading to matrix cracking. Post-mortem decomposition of oyster tissue further produces aggressive species that lower surface pH and facilitate biogenic acid attack. Conversely, a protective mechanism operates when a high-density, mature oyster community establishes a durable biofilm-shell composite. This layer physically seals surface pores, refines pore structure, and significantly impedes the ingress of aggressive agents such as chloride ions and oxygen, while also mitigating physical erosion from water flow. The transition from net corrosion to net protection is critically governed by the achievement of a threshold oyster coverage density and biofilm maturity, with low-density attachment providing insufficient protection and potentially exacerbating localized degradation.

    The review further summarizes strategies to manage oyster attachment. Preventive measures include modifying concrete mix designs with pozzolans or slag to enhance inherent durability, physical removal methods such as manual cleaning and ultrasonic techniques, and the development of advanced protective coatings. These coatings encompass natural antifouling agent-based coatings, surface-property-tuning coatings including biomimetic and amphiphilic designs, and stimuli-responsive intelligent coatings that offer targeted, environmentally benign alternatives. However, they face challenges in long-term stability and scalable application. Conversely, inductive measures to promote protective oyster colonization focus on optimizing concrete as a settlement substrate through the incorporation of calcium-based inducers such as bovine bone powder and calcium carbonate, and tailoring surface characteristics including roughness and dark color to significantly enhance larval attachment rates and aim to engineer stable, high-density biofilms for long-term concrete preservation.

    Summary and prospects

    Although significant progress has been achieved in understanding oyster-concrete interactions, the dualistic corrosion-protection mechanisms of oyster attachment remain complex and incompletely understood, with the evolutionary dynamics of oysters as dominant sessile organisms on concrete substrates and their lifecycle impacts on concrete integrity subject to ongoing debate owing to the intricate interplay of diverse marine environmental factors, seasonal variations, and developmental phases that lead to differential growth patterns and distinct effects on concrete durability. Future research should prioritize conducting systematic in situ field trials within authentic marine environments to enable precise mechanistic analysis of adhesion and corrosion processes, facilitate quantitative assessment of oyster life-cycle effects on concrete performance, and establish comprehensive models capable of predicting the transition between corrosive and protective outcomes based on critical thresholds of coverage density, biofilm maturity, and local hydrodynamic regimes. Furthermore, although oyster attachment can be detrimental under conditions of low-density colonization, the formation of stable, dense biofilms may confer significant protective benefits to concrete substrates by impeding aggressive aggressive-agent ingress, refining surface pore structure, and mitigating physical erosion, warranting in-depth investigations into this protective mechanism through integration of biological principles such as selective breeding and behavioral conditioning of target oyster species to promote rapid colonization of intertidal concrete surfaces and engineer cohesive fouling communities capable of providing effective long-term material preservation, ultimately contributing to the development of controllable, large-scale bioprotective strategies for sustainable marine infrastructure management.

  • Zengqi ZHANG, Weijie DU, Yu LIU, Shanliang MA, Xiaoming LIU, Rui JING
    Journal of the Chinese Ceramic Society. 2026, 54(10): 3436-3452.
    Introduction

    Ultra-high-performance concrete (UHPC) is widely applied in bridges, prefabricated structures, and advanced construction owing to its high strength, dense microstructure, and excellent durability performance. Its performance is strongly governed by the type and dosage of supplementary cementitious materials (SCMs), including silica fume (SF), ground granulated blast-furnace slag (GGBS), superfine GGBS (sGGBS), fly ash (FA), and fly ash microspheres (sFA), as well as polycarboxylate superplasticizer (PCE) dosage and curing age (CA). However, the relationships among mixture parameters, water distribution, hydration, and strength development are highly nonlinear, and conventional empirical models often fail to provide accurate predictions.

    Water film thickness (WFT) is a key microscopic parameter reflecting the free water covering the surfaces of solid particles, directly influencing particle lubrication, dispersion, and fresh-state behaviour in high-solid-volume UHPC systems. Compressive strength (CS) is another essential index for mechanical evaluation. Accurate prediction of both WFT and CS is crucial for mixture design and performance optimization. Nevertheless, UHPC experiments are costly and time-consuming, leading to limited datasets. To overcome this small-sample challenge, this study proposes an integrated framework combining conditional tabular generative adversarial network (CTGAN) for data augmentation, Bayesian optimization (BO)-based eXtreme Gradient Boosting (XGBoost) for predictive modelling, and interpretable machine-learning methods (SHAP and PDPs) for mechanism analysis.

    Methods

    A total of 160 experimental data points were obtained from binary, ternary, and quaternary cementitious systems. The input variables included cement content, PCE dosage, SF, GGBS, sGGBS, FA, sFA, and CA, while the output variables were WFT and CS. P·I 42.5 cement was used as the reference cement, and quartz sand with particle sizes of 0.84-2.00 mm was used as the aggregate. In the binary systems, one mineral admixture was used to replace cement by volume at replacement levels of 0%, 5%, 10%, 15%, and 20%. In the ternary systems, the SF content was fixed at 10%. In the quaternary systems, the GGBS content was fixed at 20% based on the ternary system. WFT was calculated using the wet packing method, and CS was measured after curing for 3, 7, 28 d, and 90 d.

    CTGAN was used for data augmentation to expand the original small dataset. The model learned the joint distribution of the real experimental data through adversarial training between the generator and discriminator. After optimization, 1000 synthetic samples were generated for WFT prediction, and 900 synthetic samples were generated for CS prediction. The quality of the synthetic data was evaluated by comparing histograms, cumulative distribution curves, and correlation matrices between real and synthetic datasets.

    The synthetic datasets were used to train XGBoost regression models. BO was applied to optimize hyperparameters, including tree depth, learning rate, number of estimators, subsample ratio, column-sampling ratio, minimum child weight, and regularization. Model performance was evaluated by R2, MAPE, RMSE, and MAE. SHAP and PDPs were employed to interpret feature contributions and nonlinear effects.

    Results and discussion

    The comparison between real and synthetic datasets showed that CTGAN effectively preserved the main distributional characteristics of the original data. The generated samples had value ranges and distribution patterns similar to those of the real samples. The cumulative distribution curves also showed good consistency between real and synthetic datasets. In addition, the correlation matrices indicated that CTGAN captured not only the marginal distribution of individual variables but also the correlation structure among key variables. Therefore, the generated data could provide a reliable data basis for BO-XGBoost model training under small-sample conditions.

    The BO-XGBoost model exhibited excellent prediction performance. For WFT, test R2 reached 0.986, with training and test RMSE values of 0.909 3 nm and 0.387 6 nm, respectively. For CS, test R2 was 0.976, with RMSE values of 5.244 8 MPa and 3.837 9 MPa. Predicted values clustered around the 45° ideal line, with most samples within ±10% error. Residuals were concentrated and symmetric, indicating no systematic bias, and loss curves showed stable convergence without overfitting. These results demonstrate that CTGAN augmentation effectively improved the ability of XGBoost to capture complex nonlinear relationships under small-sample conditions.

    The SHAP and PDP analyses further revealed the influence mechanisms of different variables. For WFT, sFA showed a strong positive contribution, mainly because its fine and spherical particles improved particle dispersion and increased effective water-film formation. SF had a positive effect at low-to-medium dosages, but excessive SF could reduce WFT because of its high specific surface area. GGBS and sGGBS tended to show negative effects at high dosages, which may be related to denser particle packing and compression of effective water-film space. PCE showed a nonmonotonic effect: moderate PCE improved WFT, while excessive PCE could weaken the stability of the water-film structure.

    For CS, CA was the most important factor and showed a strong positive contribution, reflecting continuous hydration and microstructure densification with curing time. sGGBS, SF, and sFA generally promoted strength development. sGGBS contributed to strength through microfilling and potential pozzolanic activity, while SF improved strength through pore refinement and matrix densification. In contrast, GGBS and FA showed negative effects at high dosages because excessive cement replacement reduced clinker content and limited the availability of calcium hydroxide (Ca(OH)2) for secondary reactions. These results demonstrate that different mineral admixtures have dosage-dependent and material-specific effects on UHPC performance.

    Conclusions

    This study establishes a multi-response prediction and interpretation framework integrating CTGAN, BO-XGBoost, SHAP, and PDPs for UHPC under small-sample conditions. CTGAN effectively generated high-quality synthetic data while preserving original distribution and correlation characteristics. The optimized BO-XGBoost model achieved outstanding predictive accuracy and generalization, with test R2 of 0.986 for WFT and 0.976 for CS. The interpretability analysis demonstrated that the effects of mineral admixtures on WFT and CS are highly dependent on material type and dosage: sFA and SF are beneficial for WFT within suitable ranges, whereas sGGBS, SF, and sFA promote CS; excessive GGBS and FA impair CS due to limited reaction activity. The proposed framework offers an effective data-driven approach for UHPC mixture design, performance prediction, and mechanism interpretation when only limited experimental data are available.

  • Zhendong JIANG, Yonggang JIANG, Mengyao ZHAO, Junzong FENG, Liangjun LI, Yijie HU, Jian FENG
    Journal of the Chinese Ceramic Society. 2026, 54(10): 3782-3794.

    Driven by the rapid development of hypersonic vehicles and long-duration spacecraft, increasingly severe aerothermal environments have imposed unprecedented demands on thermal protection systems, requiring novel thermal insulation materials that combine low density, ultralow thermal conductivity, and excellent high-temperature durability. Aerogels, with their three-dimensional nanoporous networks, ultrahigh porosity, and suppressed solid- and gas-phase heat transport, are promising candidates for thermal protection applications. Among them, conventional silica aerogels, despite their excellent insulating properties, are typically limited to service temperatures below 800 ℃ due to structural densification and performance degradation at elevated temperatures. In contrast, alumina aerogels can operate above 1000 ℃ while maintaining low thermal conductivity and demonstrating superior high-temperature stability, making them a focal point in the field of high-temperature thermal insulation.

    However, alumina aerogels still face serious challenges under high-temperature conditions. During heating, alumina undergoes successive phase transitions; in particular, above approximately 1000 ℃, the transformation to thermodynamically stable α-Al2O3 involves reconstruction from a cation-vacancy-rich spinel-like structure to a close-packed corundum structure. This transition causes rapid densification, substantial volumetric shrinkage, and particle sintering. As a result, the specific surface area and porosity collapse rapidly, causing the thermal insulation performance to deteriorate sharply at high temperatures, which severely restricts practical application.

    To address these issues, recent efforts to improve the high-temperature performance of alumina aerogels can be grouped into three main strategies: 1) structural-unit engineering, in which highly crystalline or anisotropic building blocks are designed to enhance the thermal stability of the aerogel skeleton; 2) extrinsic doping, in which elements such as Si, Ti, Zr, and La are introduced to suppress the nucleation and grain growth of α-Al2O3; and 3) isolation phase construction, in which thermally stable secondary phases are formed between alumina particles to retard sintering after phase transformation. This review summarizes recent progress in these three strategies, analyzes of their underlying mechanisms, and discusses future directions for developing alumina aerogels with higher temperature resistance.

    Summary and prospects

    This review systematically summarizes recent advances in enhancing the temperature resistance of alumina aerogels through three key strategies: structural-unit engineering, extrinsic doping, and isolation phase construction. These approaches have effectively delayed alumina phase transformation and suppressed high-temperature sintering, increasing the short-term service temperature from approximately 1000 ℃ to as high as 1600 ℃ while maintaining low thermal conductivity and structural integrity at elevated temperatures. However, for applications requiring higher temperatures and longer exposure times, future research should prioritize multiscale synergistic design, infrared opacification, mechanical reinforcement, and deeper investigation of the microscopic mechanisms of high-temperature phase transformation and sintering. Integrating the integration of emerging technologies such as multiscale simulation and artificial intelligence is expected to accelerate the development of alumina aerogels that simultaneously achieve ultrahigh-temperature structural stability and durable thermal insulation functionality.

  • Wei ZHAO, Linran HUANG, Kaiqi ZOU, Jiamin WU, Qingsong WEI
    Journal of the Chinese Ceramic Society. 2026, 54(9): 3013-3031.

    Binder Jetting (BJ) additive manufacturing emerges as a promising forming process for fabricating geometrically complex ceramic components, including thin walls, internal channels, and topology-optimized features. The process achieves forming through the selective deposition of a liquid binder onto a powder bed, typically without involving a melt pool. Consequently, it is inherently compatible with high-melting-point, brittle ceramic powders and is less constrained by melt-related instabilities. Despite these advantages, achieving robust manufacturing remains challenging. Defects and variability rarely originate from an isolated processing step. They emerge from the interaction of a series of strongly coupled steps like powder spreading, binder deposition and infiltration, drying/curing and depowdering, thermal removal of organics (debinding), and sintering. Heterogeneities introduced early in the process chain are often transmitted and amplified, ultimately manifesting as non-uniform shrinkage, warpage, cracking, residual porosity, and property scatter. In this context, full-process numerical modeling is increasingly required, serving for mechanistic interpretation and for predictive process-window design.

    This review consolidates recent advances on numerical modeling of ceramic Binder Jetting (BJ) and proposes an organizing principle to enhance the coherence and utility of the literature, framing the process chain through transferable state variables that quantify process inheritance. Modeling approaches are evaluated not as isolated stage-specific problems, but according to three dimensions, i.e., 1) the specific physical quantities predicted; 2) the transferability of these quantities as consistent inputs to subsequent stages; and 3) the practicality of calibrating and validating the relevant quantities through accessible experiments. Key transferable state variables highlighted include the spatial distribution of powder-bed porosity and pore-size statistics; permeability and its directional dependence; local binder saturation and organic distribution in the green body; evolving transport coefficients and gas-generation rates during debinding; and the initial relative-density and defect fields that control densification and shape evolution during sintering. Centering on these variables enables an explicit clarification of inter-stage coupling logic and charts a path toward engineering-grade predictability.

    Powder spreading modeling defines the initial mesoscale structure controlling all downstream responses. The discrete element method (DEM) is widely adopted as it enables the direct representation of particle rearrangement, contact interactions, and localized compaction under blades or rollers. The existing studies demonstrate that particle size distribution, morphology, cohesion, layer thickness, and recoater kinematics jointly regulate bed density, surface roughness and structural anisotropy. For full-process simulation, key advantage of DEM converts these attributes into downstream-ready descriptors-including porosity maps, pore-network statistics and anisotropic permeability closures.

    Binder deposition and infiltration determine the saturation field and organic distribution, which simultaneously control green body strength and subsequent gas evolution. At the component level, Darcy-type infiltration models remain an effective tool for rapidly exploring printing conditions. When permeability and capillary effects are expressed through measurable parameters (e.g., viscosity, surface tension, and wetting), these models can be enhanced to incorporate drying, evaporation, and curing/solidification kinetics. When pore-scale interfacial dynamics and droplet impact dominate, higher-resolution approaches are required. Interface-capturing computational fluid dynamics (CFD) and lattice Boltzmann method (LBM) simulations have been used to resolve transient interfaces and flow pathways in complex pore geometries, while coupled DEM-CFD frameworks have been introduced to capture infiltration-induced particle rearrangement in fine powders.

    Debinding is typically the most failure-prone stage since gases generated by binder decomposition must vent through a pore network while the solid skeleton remains weak and thermally stressed. Consequently, predictive modeling necessitates a coupled treatment of reaction kinetics, transport and mechanical response. Decomposition kinetics is typically constrained through thermogravimetric analysis (TGA), commonly implemented with multi-step schemes that represent complex degradation pathways. Gas transport is treated through diffusion, convection, and pressure-driven Darcy flow, with permeability and effective diffusivity evolving with porosity and temperature. Coupled heat and mass transfer in actual geometries can be solved through the finite element method (FEM) implementations, with a failure risk evaluated by relating pore pressure, strength reduction, and constrained deformation.

    Sintering consolidates the debound microstructure and determines densification, grain growth, and macroscopic shape change. Continuum constitutive laws embedded within FEM remain a dominant route for component-scale prediction, providing that temperature-dependent viscosities/creep parameters are calibrated and initial relative-density fields are supplied. These models are essential for predicting anisotropic shrinkage and warpage induced by inherent heterogeneity. In the microstructural scale, the Monte Carlo Potts models and phase-field approaches reveal the mechanisms of neck growth, pore evolution, grain-boundary migration, and coarsening. Despite high computational cost, they provide a basis for reduced-order kinetics and constitutive relations applicable to engineering-scale simulations, linking fundamental mechanisms to manufacturability.

    A consistent thread across all stages is that predictive capability relies on effective parameter identification and validation. Accordingly, this review recommends establishing a validation matrix aligned with transferable variables. Powder spreading should be validated via bed density, surface roughness, and anisotropy metrics; infiltration via penetration depth, saturation mapping, and spatial green strength, debinding via multi-rate TGA/DTG constraints, independent characterization of pore-network evolution, and statistical measures of cracking/distortion as well as sintering via shrinkage anisotropy, warpage, final density, and defect populations. Designing validation around variables required for downstream transfer allows uncertainty to be localized and gradually reduced, thus avoiding implicit redistribution across the chain.

    Summary and Prospects

    Future developments in the field are likely to follow three main trajectories. First, standardized definitions and transfer rules for transferable variables are required to reduce inconsistencies among models and to support reproducible coupled workflows. Second, the integration of uncertainty-aware calibration and sensitivity analysis into process-window design is warranted as fundamental constitutive relations and kinetic parameters are constrained by measurement limitations and condition dependence. Third, integrated workflows combining high-fidelity simulation, reduced-order modeling, and data-driven surrogates are essential for fast component-scale screening and optimization. Driven by these advances, full-process modeling can evolve from retrospective interpretation toward prospective control of process-structure-property linkages in ceramic BJ, thereby improving reliability, dimensional fidelity, and performance consistency across diverse material categories.

  • Peihang ZHAO, Yixi ZHUANG, Rongjun XIE
    Journal of the Chinese Ceramic Society. 2026, 54(9): 3087-3107.

    Persistent luminescence and mechanoluminescence are two distinct luminescent phenomena that do not require continuous in-situ optical excitation. They generate optical signals through the storage and delayed release of charge carriers from traps and mechanical stimulus-induced energy conversion, respectively. Different excitation mechanisms and temporal response characteristics both share fundamental similarities in the processes of carrier generation, trapping, and detrapping. Advances on nanomaterial synthesis and defect engineering enable the extension of persistent luminescent and mechanoluminescent materials from bulk and microscale systems to the nanoscale systems, where significantly enhance size, interface, and trap effects provide new opportunities for performance optimization and application expansion. This review systematically summarizes recent progress on persistent luminescent and mechanoluminescent nanomaterials, with a focus on typical nanoscale preparation strategies, luminescence mechanisms, and a unified understanding of carrier detrapping processes. Representative applications in low-background bioimaging and biomedicine, stress-responsive structural health monitoring, flexible devices, and information display are also outlined. Finally, key challenges and future trends are discussed in terms of performance preservation and structural synergy in nanoscale, quantitative characterization of luminescence properties, and engineering-oriented applications.

    Summary and prospects

    This review summarizes recent progress on the synthesis, mechanisms, and applications of persistent luminescent and mechanoluminescent nanomaterials. We emphasize that advances across both top-down and bottom-up routes substantially improve control over particle size, morphology, and-critically-the defect/trap structures governing carrier trapping and detrapping. We highlight a trap-centric perspective that links delayed emission in persistent luminescence with stress-assisted carrier release in mechanoluminescence via comparing the physical bases of the two phenomena. Translating these materials to the nanoscale can broad their utility, particularly for autofluorescence-free bioimaging and diagnostics, as well as self-powered, visually intuitive stress sensing and flexible optoelectronic devices. Collectively, these studies illustrate a shift from viewing luminescence as an isolated metric toward engineering integrated, stimulus-responsive nano-systems tailored to application demands.

    Some key challenges remain despite a rapid development. A central limitation is the trade-off between achieving precise nanoscale structural control and preserving the deep, stable trap populations needed for strong afterglow or mechanoluminescent output. In addition, a predictive and quantitative understanding of how interfaces, surface states, and localized mechanical fields regulate carrier dynamics in nanomaterials is still lacking. For practical deployment, a long-term stability under environmental and cyclic loading, biocompatibility, and scalable yet reproducible fabrication must be addressed. Future work should involve: 1) advancing in-situ/operando characterization and computation to map trap distributions and carrier pathways, 2) developing hybrid or composite architectures that balance structural integrity with optical functionality, and 3) validating application-driven prototypes in realistic settings from deep-tissue theranostics to structural health monitoring.