Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reasonable accuracy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Computing x-ray absorption spectra from linear-response particles atop optimized holes

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn–Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. Herein, this is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ~0.6 eV rms error for the K-edges of C–F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Viscosity Measurement using Drop Coalescence in Microgravity

We present in here details of a new method, using drop coalescence, for application in microgravity environment for determining the viscosity of highly viscous undercooled liquids. The method has the advantage of eliminating heterogeneous nucleation at container walls caused by crystallization of undercooled liquids during processing. Also, due to the rapidity of the measurement, homogeneous nucleation would be avoided. The technique relies on both a highly accurate solution to the Navier-Stokes equations as well as on data gathered from experiments conducted in near zero gravity environment. The liquid viscosity is determined by allowing the computed free surface shape relaxation time to be adjusted in response to the measured free surface velocity of two coalescing drops. Results are presented from two validation experiments of the method which were conducted recently on board the NASA KC-135 aircraft. In these tests the viscosity of a highly viscous liquid, such as glycerine at different temperatures, was determined to reasonable accuracy using the liquid coalescence method. The experiments measured the free surface velocity of two glycerine drops coalescing under the action of surface tension alone in low gravity environment using high speed photography. The free surface velocity was then compared with the computed values obtained from different viscosity values. The results of these experiments were found to agree reasonably well with the calculated values.

Antar, Basil N.↗
NASA NTRS1993

Accuracy considerations in the computational analysis of jet noise

The application of computational fluid dynamics methods to the analysis of problems in aerodynamic noise has resulted in the extension and adaptation of conventional CFD to the discipline now referred to as computational aeroacoustics (CAA). In the analysis of jet noise accurate resolution of a wide range of spatial and temporal scales in the flow field is essential if the acoustic far field is to be predicted. The numerical simulation of unsteady jet flow has been successfully demonstrated and many flow features have been computed with reasonable accuracy. Grid refinement and increased solution time are discussed as means of improving accuracy of Navier-Stokes solutions of unsteady jet flow. In addition various properties of different numerical procedures which influence accuracy are examined with particular emphasis on dispersion and dissipation characteristics. These properties are investigated by using selected schemes to solve model problems for the propagation of a shock wave and a sinusoidal disturbance. The results are compared for the different schemes.

Scott, James N.↗
NASA NTRS2014

Assessment of the Accuracy of the Conventional Ray-Tracing Technique: Implications in Remote Sensing and Radiative Transfer Involving Ice Clouds.

A fundamental problem in remote sensing and radiative transfer simulations involving ice clouds is the ability to compute accurate optical properties for individual ice particles. While relatively simple and intuitively appealing, the conventional geometric-optics method (CGOM) is used frequently for the solution of light scattering by ice crystals. Due to the approximations in the ray-tracing technique, the CGOM accuracy is not well quantified. The result is that the uncertainties are introduced that can impact many applications. Improvements in the Invariant Imbedding T-matrix method (II-TM) and the Improved Geometric-Optics Method (IGOM) provide a mechanism to assess the aforementioned uncertainties. The results computed by the II-TMþIGOM are considered as a benchmark because the IITM solves Maxwell's equations from first principles and is applicable to particle size parameters ranging into the domain at which the IGOM has reasonable accuracy. To assess the uncertainties with the CGOM in remote sensing and radiative transfer simulations, two independent optical property datasets of hexagonal columns are developed for sensitivity studies by using the CGOM and the II-TMþIGOM, respectively. Ice cloud bulk optical properties obtained from the two datasets are compared and subsequently applied to retrieve the optical thickness and effective diameter from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Additionally, the bulk optical properties are tested in broadband radiative transfer (RT) simulations using the general circulation model (GCM) version of the Rapid Radiative Transfer Model (RRTMG) that is adopted in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM, version 5.1). For MODIS retrievals, the mean bias of uncertainties of applying the CGOM in shortwave bands (0.86 and 2.13 micrometers) can be up to 5% in the optical thickness and as high as 20% in the effective diameter, depending on cloud optical thickness and effective diameter. In the MODIS infrared window bands centered at 8.5, 11, and 12 micrometers biases in the optical thickness and effective diameter are up to 12% and 10%, respectively. The CGOM-based simulation errors in ice cloud radiative forcing calculations are on the order of 10Wm(exp 2).

remote sensing↗
DOE OSTI2024

Multi-fidelity Fourier neural operator for fast modeling of large-scale geological carbon storage

Deep learning-based surrogate models have been widely applied in geological carbon storage (GCS) problems to accelerate the prediction of reservoir pressure and CO2 plume migration. Large amounts of data from physics-based numerical simulators are required to train a model to accurately predict the complex physical behaviors associated with this process. In practice, the available training data are always limited in large-scale 3D problems due to the high computational cost. Therefore, we propose to use a multi-fidelity Fourier neural operator (FNO) to solve large-scale GCS problems with more affordable multi-fidelity training datasets. FNO has a desirable grid-invariant property, which simplifies the transfer learning procedure between datasets with different discretization. Here, we first test the model efficacy on a GCS reservoir model being discretized into 110 k grid cells. The multi-fidelity model can predict with accuracy comparable to a high-fidelity model trained with the same amount of high-fidelity data with 81% less data generation costs. We further test the generalizability of the multi-fidelity model on a same reservoir model with a finer discretization of 1 million grid cells. This case was made more challenging by employing high-fidelity and low-fidelity datasets generated by different geostatistical models and reservoir simulators. We observe that the multi-fidelity FNO model can predict pressure fields with reasonable accuracy even when the high-fidelity data are extremely limited. The findings of this study can help for better understanding of the transferability of multi-fidelity deep learning surrogate models.

58 GEOSCIENCES↗
DOE OSTI2026

Benchmarking of three DWM-based wake models at below-rated wind speeds

Wind turbine wake models are essential tools for predicting power losses and structural loads in wind farms. Among these, the dynamic wake meandering (DWM) model, included as a recommended approach in the International Electrotechnical Commission design standard, is a widely used engineering-fidelity method that balances accuracy and computational cost. This study compares the performance of three DWM-based wake model implementations (from the Technical University of Denmark, the National Renewable Energy Laboratory, and the Institute for Energy Technology) under below-rated wind speed conditions. Model predictions of wake flow, power output, and structural loads for a four-turbine row are evaluated across different ambient turbulence levels and wind-direction misalignments and compared against high-fidelity large-eddy simulation results. All three models captured the overall wake evolution and mean turbine performance with reasonable accuracy; their predicted time-averaged thrust and power were typically within 5 %–10 % of the large-eddy simulation benchmark. However, notable differences emerged in wake structure and unsteady load predictions, with discrepancies increasing for turbines further downstream. These differences highlight the importance of modelling choices such as wake summation and turbulence treatment, which strongly influence power-deficit and fatigue-load predictions. Comparison with large-eddy simulations reveals each approach's strengths and weaknesses, indicating where improvements are needed. Overall, the findings point to specific refinements for DWM models to improve their fidelity, ultimately enabling more robust wake predictions for wind farm design and operation.

17 WIND ENERGY↗
NASA NTRS1989

Compression behavior of delaminated composite plates

The response of delaminated composite plates to compressive in-plane loads was investigated. The delaminated region may be either circular or elliptical, and may be located between any two plies of the laminate. For elliptical delaminations, the axes of the ellipse may be arbitrarily oriented with respect to the applied loads. A model was developed that describes the stresses, strains, and deformation of the sublaminate created by the delamination. The mathematical model is based on a two dimensional nonlinear plate theory that includes the effects of transverse shear deformation. The model takes into account thermal and moisture induced strains, transverse pressures acting on the sublaminate, and contact between the sublaminate and plate. The solution technique used is the Ritz method. A computationally efficient computer implementation of the model was developed. The code can be used to predict the nonlinear-load-strain behavior of the sublaminate including the buckling load, postbuckling behavior, and the onset of delamination growth. The accuracy of the code was evaluated by comparing the model results to benchmark analytical solutions. A series of experiments was conducted on Fiberite T300/976 graphite/epoxy laminates bonded to an aluminum honeycomb core forming a sandwich panel. Either circles or ellipses made from Teflon film were embedded in the laminates, simulating the presence of a delamination. Each specimen was loaded in compression and the strain history of the sublaminate was recorded far into the postbuckling regime. The extent of delamination growth was evaluated by C-scan examination of each specimen. The experimental data were compared to code predictions. The code was found to describe the data with reasonable accuracy. A sensitivity study examined the relative importance of various material properties, the delamination dimensions, the contact model, the transverse pressure differential, the critical strain energy release rate, and the relative growth direction on the buckling load, the postbuckling behavior, and the growth load of the sublaminate.

Peck, Scott O.↗
NASA NTRS1974

Computer experiments in preparation of system identification from transient rotor model tests, part 2

System identification methods which can extract model rotor paramenters with reasonable accuracy from noise polluted blade flapping transient measurements were developed. Usually parameter identification requires data on the state variables, that is on deflections and on rate of deflections. The small size of rotor models makes it, however, difficult to measure more than the blade flapping deflections. For the computer experiments it was, therefore, assumed that only noisy deflection measurements are available. Parameter identifications were performed for one and two unknown parameters. Both rotating coordinates and multiblade coordinates were used. It was found that data processing with a digital filter allowed by numerical differentiation a sufficiently accurate determination of the rates of deflection and of the accelerations to obtain reasonable parameter estimates with a simple linear estimator.

Hohenemser, K. H.↗
DOE OSTI2022

Application of machine learning to evaluating and remediating models for energy and environmental engineering

Machine learning (ML) algorithms have been increasingly successful in their applications to solve energy and environmental engineering problems. ML algorithms have the advantage of being able to solve highly nonlinear issues effectively. Furthermore, considering the limited sample size of data collected in energy and environmental engineering, obtaining a ML model with reasonable accuracy is simple. Unfortunately, the vast majority of the current applications of ML algorithms lack effective screening of dominant factors and comprehensive model validation, which weakens the predictive ability of the models. The present study takes the minimum miscible pressure (MMP) of CO 2 - oil systems as an example. It establishes a systematic and robust predictive model to address this issue. Based on 147 sets of slim tube tests, the predictive models of the MMPs are investigated by application of eight ML algorithms. The paper concludes that most of the published ML models in the field of energy and environmental engineering prediction are not reliable. Furthermore, it addresses the main reasons for the poor performance of some predictive models built by ML and provides guidelines on how to make such models robust. Further, to the best of our knowledge, this is the first study to point out the defects of current ML modeling methods and propose countermeasures for their application in energy and environmental engineering problems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRSDate not supplied

Mixed-Domain Charge Transport in S-Se Alloys as a Li-S Battery Cathode Material

Lithium-sulfur batteries are emerging candidate systems for the next-generation long-range electric aircraft application owing to their potential to deliver high energy density per weight. Their cathodes are to be made primarily with sulfur, but sulfur by itself has much too low electron conductivity to serve as a useful cathode. An idea that has been suggested to overcome this bottleneck is to alloy sulfur with selenium in the hope that electron mobility, which is thought to improve charge transport at some expense of increased weight and reduced energy density. However, understanding of these alloy structures and their transport mechanism is insufficient, and electron mobility at varying degrees of selenium content is not well charted, which are critical for determining the optimum selenium content and testing whether this strategy is generally feasible. The difficulty of computationally characterizing these alloy systems is rooted in the structure. The structures of both sulfur and selenium exhibit eight-atom rings that pack together in various orientations to form their respective crystals. For one, because these crystals (and presumably their alloys as well) feature multiple polymorphs such that their structural coherence is rather unclear. Secondly, these are also relatively large systems (32 atoms per cell) with low-symmetry, which complicate computation both in terms of accuracy and efficiency. Thirdly, such semi-molecular, semi-crystalline structures lead to a combination of band-transport characteristics and hopping-transport characteristics, each of which is a domain with its own physics and set of computational, theoretical challenges. In addition, there is a general dearth of experimental transport data for sulfur-selenium alloys. In this study, we make a comprehensive attempt to tackle this problem using a wide array of first-principles methods. We generate special quasirandom structures to simulate alloy structures, and compute their first-principles Raman spectra for comparison with experimental data to ensure their structural soundness. We then proceed to use recently developed, state-of-the-art tool (AMSET) in order to efficiently compute electronic mobilities under band transport of pure sulfur, pure selenium, as well as their alloys of various compositions at a reasonable accuracy. We then sample numerous dimer configurations of nearest-neighbor eight-atom-ring-pairs and compute electronic hopping rates between them, which yields hopping mobility. A combination of these efforts lead to a general mapping of electron transport behaviors throughout the alloy range. Preliminary results show that introduction of selenium into sulfur initially damages electron mobility due to disorder but eventually improves it beyond that of pure sulfur with additional selenium content, ultimately peaking at the pure-selenium limit. Band transport dominates for holes in sulfur and electrons in selenium, but in all other cases, hopping transport is dominant. Ongoing efforts include determination of charge concentration and conductivity, as well as performing multiphysics modeling to determine the optimum selenium content for best tradeoff between conductivity and energy density for aircraft range.

Junsoo Park↗
DOE OSTI2026

Evaluating the limitations of Bayesian metabolic control analysis

Bayesian Metabolic Control Analysis (BMCA) is a promising framework for inferring metabolic control coefficients in data-limited scenarios, combining Bayesian inference with linear-logarithmic (lin-log) rate laws. These metabolic control coefficients quantify how changes in enzyme activities affect steady-state fluxes and metabolite concentrations across a metabolic network. However, its predictive accuracy and limitations remain underexplored. This study systematically evaluates BMCA’s ability to infer elasticity values, flux control coefficients (FCC), and concentration control coefficients (CCC) under varying data availability conditions using three synthetic metabolic network models. We demonstrate that BMCA predictions are highly dependent on the inclusion of flux and enzyme concentration data, with the omission of these datasets leading to severe inaccuracies. In our synthetic, enzyme-perturbation datasets, external metabolite concentrations had minimal impact and, in some cases, their exclusion improved predictions; when external-nutrient perturbations were introduced and those concentrations were observed, gains were at most modest. Additionally, we find that posterior estimation with both ADVI and HMC can underestimate large-magnitude elasticities in our synthetic settings, with ADVI showing somewhat higher variance under strong up-regulation; thus, recovering |elasticity| ≳ 1.5 remains challenging regardless of the inference engine. ADVI also fails to accurately infer allosteric interactions, even when regulatory effects are strong. While BMCA maintains reasonable accuracy in partially recovering the rankings of the highest FCC values, its estimates of absolute values remain constrained by prior assumptions and data limitations. Our findings reveal the BMCA algorithm’s strengths and weaknesses, providing guidance on its application in metabolic engineering, and highlighting the need for methodological refinements to enhance its predictive capabilities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

CpFe(CO) 2 Radical Generated from Dinuclear [CpFe(CO) 2 ] 2 and Mononuclear (Cp)(CO) 2 Fe(H): Density Functional Theory Is Accurate for One, But Not Both

Density functional theory (DFT) methods remain the most practical approach to calculating properties and reaction mechanisms of transition metal complexes. While the accuracy of DFT methods has been evaluated for some properties of mononuclear organometallic complexes there has been a general lack of evaluation for dinuclear organometallic complexes, in particular bonding changes related to reaction mechanisms. Here, this work evaluated DFT and coupled cluster methods for the accuracy of calculating the CpFe(CO) 2 radical (Fp•) generated from dinuclear [CpFe(CO) 2 ] 2 (Fp 2 ) and mononuclear [(Cp)(CO) 2 Fe(H)] (Fp-H). This transition metal radical fragment was evaluated because dinuclear complexes built with it have recently shown a variety of unique reactions but has proven challenging to accurately calculate with DFT methods. Here we show that DFT methods provide a surprising wide range of fragmentation energies for Fp 2 and lower and mid rung DFT methods as well as DLPNO–CCSD(T) perform well for this dissociation energy. The highest rung double-hybrid methods have a large range in the Fp 2 dissociation energy, and the energy greatly depends on the amount of MP2 correlation energy included. For generating Fp• from Fp-H the lower and mid rung methods that worked well for Fp 2 showed significant error. Double-hybrid methods unfortunately are only accurate for the Fe–H bond if they are very inaccurate for the Fp 2 dissociation energy. While DLPNO–CCSD(T) is not perfect, and not close to chemically accurate for the Fe–H bond, it does provide reasonable accuracy for both Fp 2 and Fp-H dissociation energies.

density functional theory↗
NASA NTRS2012

Digital Architecture for a Trace Gas Sensor Platform

A digital architecture has been implemented for a trace gas sensor platform, as a companion to standard analog control electronics, which accommodates optical absorption whose fractional absorbance equivalent would result in excess error if assumed to be linear. In cases where the absorption (1-transmission) is not equivalent to the fractional absorbance within a few percent error, it is necessary to accommodate the actual measured absorption while reporting the measured concentration of a target analyte with reasonable accuracy. This requires incorporation of programmable intelligence into the sensor platform so that flexible interpretation of the acquired data may be accomplished. Several different digital component architectures were tested and implemented. Commercial off-the-shelf digital electronics including data acquisition cards (DAQs), complex programmable logic devices (CPLDs), field-programmable gate arrays (FPGAs), and microcontrollers have been used to achieve the desired outcome. The most completely integrated architecture achieved during the project used the CPLD along with a microcontroller. The CPLD provides the initial digital demodulation of the raw sensor signal, and then communicates over a parallel communications interface with a microcontroller. The microcontroller analyzes the digital signal from the CPLD, and applies a non-linear correction obtained through extensive data analysis at the various relevant EVA operating pressures. The microcontroller then presents the quantitatively accurate carbon dioxide partial pressure regardless of optical density. This technique could extend the linear dynamic range of typical absorption spectrometers, particularly those whose low end noise equivalent absorbance is below one-part-in-100,000. In the EVA application, it allows introduction of a path-length-enhancing architecture whose optical interference effects are well understood and quantified without sacrificing the dynamic range that allows quantitative detection at the higher carbon dioxide partial pressures. The digital components are compact and allow reasonably complete integration with separately developed analog control electronics without sacrificing size, mass, or power draw.

Gonzales, Paula↗
DOE OSTI2020

Direct Quantification of Heat Generation Due to Inelastic Scattering of Electrons Using a Nanocalorimeter

Transmission electron microscopy (TEM) is arguably the most important tool for atomic-scale material characterization. A significant portion of the energy of transmitted electrons is transferred to the material under study through inelastic scattering, causing inadvertent damage via ionization, radiolysis, and heating. In particular, heat generation complicates TEM observations as the local temperature can affect material properties. Here, the heat generation due to electron irradiation is quantified using both top-down and bottom-up approaches: direct temperature measurements using nanowatt calorimeters as well as the quantification of energy loss due to inelastic scattering events using electron energy loss spectroscopy. Combining both techniques, a microscopic model is developed for beam-induced heating and to identify the primary electron-to-heat conversion mechanism to be associated with valence electrons. Building on these results, the model provides guidelines to estimate temperature rise for general materials with reasonable accuracy. This study extends the ability to quantify thermal impact on materials down to the atomic scale.

42 ENGINEERING↗
DOE OSTI2020

Machine learning to predict biomass sorghum yields under future climate scenarios

Crop yield modeling is critical in the design of national strategies for agricultural production, particularly in the context of a changing climate. Forecasting yields of bioenergy crops at fine spatial resolutions can help to evaluate near-term and long-term pathways for scaling up bio-based fuel and chemical production, and for understanding the impacts of abiotic stressors such as severe droughts and temperature extremes on potential biomass supply. In this work we used a large dataset of 28,364 Sorghum bicolor yield samples (uniquely identified by county and year of observation), environmental variables, and multiple approaches to analyze historical trends in sorghum productivity across the USA. We selected the most accurate machine learning approach (a variation of the random forest approach) to predict future trends in sorghum yields under four greenhouse gas (GHG) emission scenarios and two irrigation regimes. We identified irrigation practices, vapor pressure deficit, and time (a proxy for technological improvement) as the most important predictors of sorghum productivity. Our results showed a decreasing trend of sorghum yields over future years (on average 2.7% from 2018 to 2099), with greater decline under a high GHG emissions scenario (3.8%) and in the absence of irrigation (4.6%). Geographically, we observed the steepest predicted declines in the Great Lakes (8.2%), Upper Midwest (7.5%), and Heartland (6.7%) regions. Our study demonstrates the use of machine learning to identify environmental controllers of sorghum biomass yield and predict yields with reasonable accuracy. These results can inform the development of more realistic biomass supply projections for bioenergy if sorghum production is scaled up. (c) 2020 Society of Chemical Industry and John Wiley & Sons, Ltd

09 BIOMASS FUELS↗
DOE OSTI2021

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of spatial energy distribution-induced porosity on mechanical properties of laser powder bed fusion 316L stainless steel

Laser powder bed fusion (LPBF) additive manufacturing (AM) offers a variety of advantages over traditional manufacturing, however its usefulness for manufacturing of high-performance components is currently hampered by internal defects (porosity) created during the LPBF process that have an unknown impact on global mechanical performance. By inducing porosity distributions through variations in print energy density and inspecting the resulting tensile samples using computed tomography, nearly 50,000 pores across 75 samples were identified. Porosity characteristics were quantitatively extracted from inspection data and compared with mechanical properties to understand the strength of relationships between porosity and global tensile performance. Useful porosity characteristics were identified for prediction of part performance. Results indicate that ductility and strain at ultimate tensile strength are the global tensile properties most significantly impacted by porosity and can be predicted with reasonable accuracy using simple porosity shape descriptors such as volume, diameter, and surface area. Moreover, it was found that the largest pores influenced behavior most significantly. Specifically, pores in excess of 125 µm in diameter were found to be a sufficient threshold for property estimation. These results establish an initial understanding of the complex defect-performance relationship in AM 316L stainless steel and can be leveraged to develop certification standards and improve confidence in part quality and reliability for the broader set of engineering alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lattice Boltzmann simulation of the dissolution of slag in alkaline solution using real-shape particles

Highlights: • A dissolution numerical model was proposed to capture the real dissolution kinetics of slag in alkaline solution. • The log forward dissolution rate of Si was described as a function of NBO/T and solution pH. • A threshold solid volume fraction of 0.688 was found for a voxel in 3D, 63.8% larger than that for a pixel in 2D. • The proposed dissolution numerical model provides a reliable alternative to study the dissolution kinetics of slag. A dissolution numerical model was proposed in this study to capture the real dissolution kinetics of slag in alkaline solution. It consists of three modules, i.e. (i) simulation of the initial particle parking structure of slag in alkaline solution using real-shape particles of slag, (ii) simulation of the chemical reactions between slag and solution based on the transition state theory, and (iii) simulation of the physical transport of aqueous ions using the lattice Boltzmann method. This dissolution numerical model was verified using experimental results, showing reasonable accuracy. After verification, the dissolution numerical model was implemented to study the influences of temperature and particle shape using a proper recipe of slag in alkaline solution. This recipe was designed to avoid solid phase precipitation or gel formation via thermodynamic analysis. The simulation results showed faster dissolution kinetics of slag when using higher temperatures and more irregular particle shapes.

36 MATERIALS SCIENCE↗