Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lattice Boltzmann”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Scale translation yields insights into gas adsorption under nanoconfinement

This work describes a scale-translating simulation framework to investigate gas adsorption behavior in nanoconfined pores. The framework combines molecular simulations (MSs), equation of state (EoS), and lattice Boltzmann (LB) simulations. MSs reveal the physics of methane adsorption in nano-sized pores, where input values of fugacity coefficients are optimized based on EoS predictions. Then, an LB free-energy model, which incorporates a viral EoS, upscales intermolecular forces and estimates adsorption behavior via a proposed fluid–wall interaction model. Armed with the values of the LB interaction parameter as a function of pressure, the LB model is used to predict fluid behavior in irregular nanopores, and the results are validated against reference MS data. The LB model is then used to study adsorption behavior at a continuum scale in representative organic shale nanopores based on finely characterized Vaca Muerta shale samples. Furthermore, the results show that methane adsorption could significantly increase contained fluids by 10%–25% in pores smaller than 20 nm. However, in larger pores (40 nm to 90 nm), adsorption's impact diminishes to 2%–3%, suggesting sorption's negligible role beyond a 40 nm pore size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗
DOE OSTI2024

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗
DOE OSTI2022

A coupled DEM-IMB-LBM model for simulating methane hydrate exploitation involving particle dissolution

The coupled discrete element and lattice Boltzmann method using an immersed moving boundary scheme was extended to simulate methane hydrate exploitation involving mass transport and particle dissolution. In this coupled DEM-IMB-LBM model, a new Dirichlet-type thermal boundary condition is extended to simulate moving curved boundaries with constant concentration. A novel periodic boundary including an efficient searching algorithm for particle contact is proposed to reduce the computational cost and boundary effect. So, this model is validated by two numerical examples: a circular particle with concentration convection-diffusion moving in a horizontal channel and mass transport from a cylinder particle in a simple shear flow. The numerical results obtained from the proposed model agree well with previous studies. To further demonstrate the capacity of the proposed model, simulations of methane hydrate exploitation including two formations in marine sediments are carried out. The numerical results indicate that the coupled DEM-IMB-LBM is not only capable of simulating the dissolution of hydrate particles at the grain level, but also recover the sand erosion and migration process in a fundamental perspective during the methane hydrate exploitation process.

42 ENGINEERING↗
DOE OSTI2025

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗
DOE OSTI2022

Accurate hydrodynamic force and torque correlations for prolate spheroids from Stokes regime to high Reynolds numbers

Detailed simulations of flow around various prolate spheroids are performed using multi-relaxation time lattice Boltzmann method (LBM). The simulations are performed in the Reynolds numbers range 0.1 ≤ Re ≤ 2000 at different incident angles 0°≤ φ ≤ 90° for various prolate spheroids of aspect ratios 1 ≤ λ ≤ 8. The Re is based on the volume equivalent sphere diameter. From the simulations, accurate correlations for average drag, lift and torque coefficients (C D , C L , and C T respectively) are proposed. The mean deviations between the correlations and the simulation data for C D , C L , and C T are 2%, 6.5%, and 4.7% respectively. Furthermore, the correlations are fitted with the physics in mind such that they can be extrapolated beyond the regimes simulated. The fits are designed such that correlations reduce to analytical (C D and C L ) solutions available for prolate spheroids in the limit Re ≈ 0. At high Re, the correlations mimic the asymptotic flattening of C D as for spheres and therefore, we expect the correlations to be valid until the critical Re (≈ 10 5 ). Furthermore, extensive comparison of our correlations with other high-fidelity simulations from literature demonstrate the accuracy of our work compared to other correlations from literature. This work enables researchers to perform accurate unresolved Euler-Lagrangian simulations for a wide range of elongated particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Experimental and numerical investigation of fracture conductivity between non-smooth rock surfaces with and without proppant

The enhancement of fracture conductivity is vital for the efficient recovery of subsurface resources, such as geothermal energy and petroleum hydrocarbons. Proppants, granular materials injected into hydraulic fractures to maintain their conductivity, have been studied primarily in the context of smooth fractures (i.e., fractures between smooth rock surfaces). However, non-smooth fractures (i.e., fractures between rough rock surfaces) are common in geoenergy reservoirs and thus require further investigations. In this study, we conducted laboratory measurements of fracture conductivity on shale slabs with non-smooth surfaces and carried out numerical simulation using the lattice Boltzmann (LB) method, which aimed to investigate the conductivity of non-smooth fractures with and without proppants placement. When ceramic proppant with an areal concentration of 2 lb/ft 2 was placed in the fracture, the conductivity was enhanced by roughly 3-8 times compared to fractures without proppant. In fractures with proppant, gas-measured conductivity was higher than that measured with water due to proppant embedment caused by water. The experiments demonstrate the advantages of using proppant in fractures, even if the rock surface roughness can provide certain fracture conductivity via the self-propping mechanism. For fractures without proppants, high rock surface roughness is not necessarily favorable for enhancing fracture conductivity because the self-propping mechanism requires shear slip along the fracture surface. If there is no shear slip, high rock surface roughness can cause a detrimental effect on the fracture conductivity due to the interlocking effect. Utilizing advanced experimental equipment and LB modeling, this research explores the interplays between proppant placement, fracture geometry, and stress conditions to develop a comprehensive understanding of the productivity in non-smooth fractures. Further, the outcomes of this investigation indicate the importance of creating fractures with surface roughness during hydraulic fracturing and will contribute to the development of more efficient stimulation techniques for subsurface energy extraction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

An uncertainty visualization framework for large-scale cardiovascular flow simulations: A case study on aortic stenosis

We present a generalizable uncertainty quantification (UQ) and visualization framework for lattice Boltzmann method simulations of high Reynolds number vascular flows, demonstrated on a patient-specific stenosed aorta. The framework combines EasyVVUQ for parameter sampling with large-eddy simulation turbulence modeling in HemeLB, and executes ensembles on the Frontier exascale supercomputer. Spatially resolved metrics, including entropy and isosurface-crossing probability, are used to map uncertainty in pressure and wall shear stress fields directly onto vascular geometries. Two sources of model variability are examined: inlet peak velocity and the Smagorinsky constant. Inlet velocity variation produces high uncertainty downstream of the stenosis where turbulence develops, while upstream regions remain stable. Smagorinsky constant variation has little effect on the large-scale pressure field but increases WSS uncertainty in localized high-shear regions. In both cases, the stenotic throat manifests low entropy, indicative of robust identification of elevated WSS. By linking quantitative UQ measures to three-dimensional anatomy, the framework improves interpretability over conventional 1D UQ plots and supports clinically relevant decision-making, with broad applicability to vascular flow problems requiring both accuracy and spatial insight.

Hemodynamics↗
DOE OSTI2024

The Impact of Wettability on the Co-moving Velocity of Two-Fluid Flow in Porous Media

Abstract The impact of wettability on the co-moving velocity of two-fluid flow in porous media is analyzed herein. The co-moving velocity, developed by Roy et al. (Front Phys 8:4, 2022), is a novel representation of the flow behavior of two fluids through porous media. Our study aims to better understand the behavior of the co-moving velocity by analyzing simulation data under various wetting conditions. We analyzed 46 relative permeability curves based on the Lattice–Boltzmann color fluid model and two experimentally determined relative permeability curves. The analysis of the relative permeability data followed the methodology proposed by Roy et al. (Front Phys 8:4, 2022) to reconstruct a constitutive equation for the co-moving velocity. Surprisingly, the coefficients of the constitutive equation were found to be nearly the same for all wetting conditions. On the basis of these results, a simple approach was proposed to reconstruct the relative permeability of the oil phase using only the co-moving velocity relationship and the relative permeability of the water phase. This proposed method provides new information on the interdependence of the relative permeability curves, which has implications for the history matching of production data and the solution of the associated inverse problem. The research findings contribute to a better understanding of the impact of wettability on fluid flow in porous media and provide a practical approach for estimating relative permeability based on the co-moving velocity relationship, which has never been shown before.

Engineering↗
DOE OSTI2020

Investigating the Interaction Between Circulating Tumor Cells and Local Hydrodynamics via Experiment and Simulations

Introduction: The biological and mechanical properties of circulating tumor cells (CTCs) in combination with the hemodynamics affect the preference of metastatic sites in the vasculature. Despite the extensive literature on the effects of biological properties on cell adhesion, the effects of hydrodynamic forces on primary attachment remains an active area of research. Here, using simulations in conjunction with experimentation, we provide new insight into the interplay of CTCs dynamics and local hydrodynamics. Methods: A flow experiment of CTC attachment was performed within a bioprinted, double branching endothelialized vessel. Simulations of fluid flow and CTC transport in the reconstructed and idealized bifurcated vessel were respectively performed by HARVEY, our in-house massively parallel computational fluid dynamics solver. HARVEY is based on the lattice Boltzmann and finite element methods to model the fluid and cells dynamics. The immersed boundary method is employed for resolving the fluid–structure interaction. Results: CTC attachment was quantified experimentally at all regions of the complex vessel. The results demonstrate a clear preference for CTCs to attach at the branch points. To elucidate the effect of the vessel topology on the location of attachment, a fluid-only simulation was performed assessing the differences in the hydrodynamics along the vessel. CTC transport in idealized bifurcated vessels was subsequently studied to examine the effects of cell deformability on the local hydrodynamics patterns and, thus, the preference of attachment sites. Conclusions: The current work provides evidence on the correlation of the hydrodynamics forces arising from the vessel topology and CTC properties on the attachment regions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Controlling solute channel formation using magnetic fields

Solute channel formation introduces compositional and microstructural variations in a range of processes, from metallic alloy solidification, to salt fingers in ocean and water reservoir flows. Applying an external magnetic field interacts with thermoelectric currents at solid/liquid interfaces generating additional flow fields. This thermoelectric (TE) magnetohydrodynamic (TEMHD) effect can impact on solute channel formation, via a mechanism recently drawing increasing attention. To investigate this phenomenon, we combined in situ synchrotron X-ray imaging and Parallel-Cellular-Automata-Lattice-Boltzmann based numerical simulations to study the characteristics of flow and solute transport under TEMHD. Observations suggest the macroscopic TEMHD flow appearing ahead of the solidification front, coupled with the microscopic TEMHD flow arising within the mushy zone are the primary mechanisms controlling plume migration and channel bias. Two TE regimes were revealed, each with distinctive mechanisms that dominate the flow. Further, we show that grain orientation modifies solute flow through anisotropic permeability. These insights led to a proposed strategy for producing solute channel-free solidification using a time-modulated magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗
DOE OSTI2021

Reconstructing porous media using generative flow networks

One area of intense scientific interest for the study of sandstones, carbonates, and shale at the pore scale is the use of limited image and petrophysical data to generate multiple realizations of a rock’s pore structure. Such images aid efforts to quantify uncertainty in petrophysical properties, including porosity–permeability transforms. We develop and evaluate a deep learning-based method to synthesize porous media volumes using a so-called generative flow model trained on x-ray microscope images of rock texture and pore structure. These models are optimized on a log-likelihood objective and they synthesize large and realistic three-dimensional images. Further, we demonstrate the rapid generation of sandstone image volumes that display realism as gauged by quantitative comparison of topological features using Minkowski functionals of porosity, specific surface area, and the Euler–Poincaré characteristic (i.e., pore connectivity). We also evaluate the single-phase permeability using Navier–Stokes and lattice Boltzmann methods and show that transport properties of the generated samples match measured trends.

58 GEOSCIENCES↗
DOE OSTI2021

GeoMicro3D: A novel numerical model for simulating the reaction process and microstructure formation of alkali-activated slag

Highlights: • The GeoMicro3D model was innovatively developed for simulating the reaction process and microstructure formation of AAS. • In GeoMicro3D, the dissolution of slag and reactions of aqueous ions are fundamentally described in the local lattice cell. • A novel approach is proposed to speed up the nucleation simulations. • GeoMicro3D was implemented and verified with the relevant experimental data and thermodynamic calculation results. For the first time, this study developed a novel model, named GeoMicro3D, to simulate the reaction process and microstructure formation of alkali-activated slag. The GeoMicro3D model consists of four modules that are designed to simulate, respectively: (i) the initial spatial distribution of real-shape slag particles in alkaline activator, (ii) the dissolution of slag and diffusion of ions via the transition state theory and lattice Boltzmann method, respectively, (iii) the spatial distribution of reaction products using a nucleation probability theory, and (iv) the chemical reactions with thermodynamic modelling. Afterwards the GeoMicro3D model was implemented and verified. The simulation results were discussed and compared with the relevant experimental data and thermodynamic calculation results using GEMS. A good agreement was found in the comparisons, showing the strong simulation capability of GeoMicro3D.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗