Search NASA⌕ Search

SEARCH · Search NASA

Results for “stack”

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 433 records · Page 24

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

A computational investigation of high-flux, plate-and-frame membrane modules for industrial carbon capture

In this work, we study the application of membrane-based separation systems for carbon capture, considering plate-and-frame membrane modules. The successful deployment of membrane CO 2 capture system relies on high-performing membranes as well as effective membrane modules that can fully exploit the developed membranes. A plate-and-frame membrane module is especially attractive for CO 2 capture from industrial flue gas due to its lower pressure drop compared to its counterparts such as spiral wound modules and hollow fiber modules. To design better plate-and-frame modules, we investigate their basic unit - a single membrane stack through a combination of computational modeling and experimental investigations. The modeling approach is based on Computational Fluid Dynamics (CFD) to represent a multiphysics problem, including the fluid flow and diffusion processes within a membrane module. We use experimental data collected under different operating conditions to validate the CFD model. Numerical results suggest a good agreement between experiments and model outputs for the CO 2 recovery, CO 2 mole fraction in the retentate and permeate, and stage-cut. The CFD model is able to predict accurately the flow behavior, providing valuable insights on the effects of fluid dynamics on mass transfer of CO 2 . We also carry out a sensitivity analysis to identify the effect of key parameters on the CO 2 recovery and the CO 2 purity of the outlet streams.

CFD simulation↗

Bioelectrochemical crossbar architecture screening platform for extracellular electron transfer

Electroactive microbes can serve as living components in bioelectronic devices, where their unique ability to transfer electrons enables applications in sensing, energy conversion, and synthesis, but they remain challenging to engineer because the bioelectrochemical systems (BESs) used for characterization are low throughput. Here, we present a bioelectrochemical crossbar architecture screening platform (BiCASP) that uses stacked and orthogonally arrayed electrodes to enable individual sample selection for characterization in arrayed formats. This device reports on the current generated by electroactive bacteria on the minute timescale, decreasing the time for data acquisition by several orders of magnitude compared to conventional BESs. This device increases the throughput of screening engineered biological components in cells, identifying mutants of the membrane protein wire MtrA in Shewanella oneidensis that retain the ability to support extracellular electron transfer (EET). BiCASP may be integrated with bioelectronics that need directed evolution of electroactive proteins.

Shewanella↗

Graphite crystals in catalytically-graphitized glass-like carbon

Catalytic graphitization of glass-like carbon leads to enhanced growth of micro-sized graphitic crystals with unusual shapes of wires, filaments, tubes, rods, whiskers, and spirals. Similar particles with axial symmetry are also found in pure glass-like carbon heat-treated at high temperatures. Nonetheless, the presence of a graphitization catalyst, Si in this case, in the heat-treatment process supports the transformation of porous, disordered carbon structure towards the graphitic atomic order and the formation of manifold peculiar polyhedral wires and particles of geometry distinct from the plate-like shape typical for conventional graphite. In contrast to conventional carbon nanotubes and fibers, the graphene layers are stacked perpendicular to the tube axis, while the size of the most common tube fibers can reach up to 10 μm in diameter and 100 μm in length. X-ray diffraction, Raman spectroscopy, scanning and transmission electron microscopy, small-angle X-ray scattering combined with complementary techniques have been used to characterize the structure of the glass-like carbon derived from furfuryl alcohol catalytically-graphitized using Si particles at 3000 °C. Finally, since control of graphite shape is vital to achieving the level of performance required in contemporary applications, the obtained results demonstrate that the catalytic graphitization method may be employed to produce filamentous graphite crystals.

36 MATERIALS SCIENCE↗

TPCpp-10M: Simulated proton-proton collisions in a time projection chamber for AI foundation models

Scientific foundation models hold great promise for advancing nuclear and particle physics by improving analysis precision and accelerating discovery. Yet, progress in this field is often limited by the lack of openly available large scale datasets, as well as standardized evaluation tasks and metrics. Furthermore, the specialized knowledge and software typically required to process particle physics data pose significant barriers to interdisciplinary collaboration with the broader machine learning community. This work introduces a large, openly accessible dataset of 10 million simulated proton-proton collisions, designed to support self-supervised training of foundation models. To facilitate ease of use, the dataset is provided in a common NumPy format. In addition, it includes 70,000 labeled examples spanning three well defined downstream tasks: track finding, particle identification, and noise tagging, to enable systematic evaluation of the foundation model's adaptability. The simulated data are generated using the Pythia Monte Carlo event generator at a center of mass energy of $\sqrt{s}$ = 200 GeV and processed with Geant4 to include realistic detector conditions and signal emulation in the sPHENIX Time Projection Chamber at the Relativistic Heavy Ion Collider, located at Brookhaven National Laboratory. This dataset resource establishes a common ground for interdisciplinary research, enabling machine learning scientists and physicists alike to explore scaling behaviors, assess transferability, and accelerate progress toward foundation models in nuclear and high energy physics. The complete simulation and reconstruction chain is reproducible with the sPHENIX software stack. All data and code locations are provided under Data Accessibility.

Data Analysis, Statistics and Probability (physics↗

Intramolecular redox-site interplay effect on organic electrode for fast-charging and wide-temperature-range sodium-ion batteries

Organic electrode materials (OEMs) hold great promise for sodium-ion batteries (SIBs) due to their exceptional structural tunability and sustainability. However, the development of OEMs with fast redox kinetics and robust structural integrity remains challenging, especially over a wide temperature range. Herein, we propose an effective strategy to address both sluggish redox kinetics and insufficient structural stability in OEMs by constructing an intramolecular redox-site interplay effect. This effect is demonstrated by two hexaazatrinaphthylene-carboxylate isomers, namely HATN-m-COONa and HATN-o-COONa. Systematic experimental and computational results jointly reveal the intramolecular redox-site interplay effect in HATN-o-COONa decreases the rigid π-π stacking interactions and minimizes the skeleton structural distortion, offering faster redox kinetics and enhanced structural integrity in HATN-o-COONa compared to HATN-m-COONa (without intramolecular redox-site interplay effect). Consequently, HATN-o-COONa exhibits superior rate performance (258 mA h g−1 at 10 A g−1) and enhanced cycle stability (93% after 1000 cycles at 5 A g−1) compared to HATN-m-COONa. More importantly, HATN-o-COONa demonstrates exceptional wide-temperature adaptability, ranging from -40 °C (315 mA h g−1 at 0.1 A g−1) to 60 oC (343 mAh g-1 at 5 A g-1). This work establishes a promising design rationale for developing fast-charging and wide-temperature adaptable OEMs for energy storage systems.

Gao, Yawei [ORNL] (ORCID:0000000225672853)↗

Analyzing inference workloads for spatiotemporal modeling

Ensuring power grid resiliency, forecasting climate conditions, and optimization of transportation infrastructure are some of the many application areas where data is collected in both space and time. Spatiotemporal modeling is about modeling those patterns for forecasting future trends and carrying out critical decision-making by leveraging machine learning/deep learning. Once trained offline, field deployment of trained models for near real-time inference could be challenging because performance can vary significantly depending on the environment, available compute resources and tolerance to ambiguity in results. Users deploying spatiotemporal models for solving complex problems can benefit from analytical studies considering a plethora of system adaptations to understand the associated performance-quality trade-offs. To facilitate the co-design of next-generation hardware architectures for field deployment of trained models, it is critical to characterize the workloads of these deep learning (DL) applications during inference and assess their computational patterns at different levels of the execution stack. In this paper, we develop several variants of deep learning applications that use spatiotemporal data from dynamical systems. We study the associated computational patterns for inference workloads at different levels, considering relevant models (Long short-term Memory, Convolutional Neural Network and Spatio-Temporal Graph Convolution Network), DL frameworks (Tensorflow and PyTorch), precision (FP16, FP32, AMP, INT16 and INT8), inference runtime (ONNX and AI Template), post-training quantization (TensorRT) and platforms (Nvidia DGX A100 and Sambanova SN10 RDU). Overall, our findings indicate that although there is potential in mixed-precision models and post-training quantization for spatiotemporal modeling, extracting efficiency from contemporary GPU systems might be challenging. Instead, co-designing custom accelerators by leveraging optimized High Level Synthesis frameworks (such as SODA High-Level Synthesizer for customized FPGA/ASIC targets) can make workload-specific adjustments to enhance the efficiency.

97 MATHEMATICS AND COMPUTING↗

Bridging paradigms: Designing for HPC-Quantum convergence

Here, this paper presents a comprehensive software stack architecture for integrating quantum computing (QC) capabilities with High-Performance Computing (HPC) environments. While quantum computers show promise as specialized accelerators for scientific computing, their effective integration with classical HPC systems presents significant technical challenges. We propose a hardware-agnostic software framework that supports both current noisy intermediate-scale quantum devices and future fault-tolerant quantum computers, while maintaining compatibility with existing HPC workflows. The architecture includes a quantum gateway interface, standardized APIs for resource management, and robust scheduling mechanisms to handle both simultaneous and interleaved quantum–classical workloads. Key innovations include: (1) a unified resource management system that efficiently coordinates quantum and classical resources, (2) a flexible quantum programming interface that abstracts hardware-specific details, (3) A Quantum Platform Manager API that simplifies the integration of various quantum hardware systems, and (4) a comprehensive tool chain for quantum circuit optimization and execution. We demonstrate our architecture through implementation of quantum–classical algorithms, including the variational quantum linear solver, showcasing the framework’s ability to handle complex hybrid workflows while maximizing resource utilization. This work provides a foundational blueprint for integrating QC capabilities into existing HPC infrastructures, addressing critical challenges in resource management, job scheduling, and efficient data movement between classical and quantum resources.

97 MATHEMATICS AND COMPUTING↗

Oak Ridge Computing Academy: An HPC cluster deployment and management pilot

The High Performance Computing Technologies (HPCT) course is a hands-on High Performance Computing (HPC) cluster deployment and management training program offered as part of the International School for Advanced Studies (SISSA) and the International Center for Theoretical Physics (ICTP) Master in High Performance Computing (MHPC) specialization. Here, this training program introduces students to key concepts in cluster configuration. which include networking, software stack provisioning, job scheduling, and monitoring. The publicly available course materials feature several examples and underlying methods that are broadly applicable to cluster deployment and management. This paper discusses the design of a new workforce development program at the Oak Ridge National Laboratory that is based on HPCT, the Oak Ridge Computing Academy (ORCA). The ORCA pilot program was hosted by the Oak Ridge Leadership Computing Facility (OLCF) in Summer 2025. As a part of this discussion, HPCT and ORCA course contents and infrastructure are outlined, ORCA participant experiences are detailed, and potential opportunities for improvement are discussed.

Education↗

A two-stage optical fusion framework for wildfire severity mapping across the conterminous United States

Accurate wildfire severity mapping (WSM) is essential for post-fire recovery planning, erosion risk assessment, ecosystem monitoring, and disaster risk reduction. Although Landsat and Sentinel optical imagery have been widely used for burn severity assessment, the added value of fusing multiple optical sensors has not been sufficiently quantified across diverse fire events, particularly since the launch of Landsat-9. This study evaluates whether multisensor optical fusion improves wildfire severity mapping relative to single-sensor baselines using Sentinel-2, Landsat-8, and Landsat-9 imagery across 40 wildfire events in the conterminous United States. We tested a two-stage fusion framework that combines feature-level fusion with pixel-level dimensionality reduction. First, feature-level fused datasets were created through early fusion by combining standardized post-fire bands from each sensor into a single predictor stack. Both raw reflectance bands and pairwise spectral transforms were retained to capture within- and cross-sensor spectral interactions. Second, Linear Discriminant Analysis was applied to both single-sensor and fused datasets to produce comparable low-dimensional feature spaces. Six machine-learning classifiers were then used to benchmark model performance with repeated spatially buffered train–test splits. Results show that Landsat-9 was the strongest single-sensor baseline. Among the fusion strategies, Sentinel-2 + Landsat-9 produced the most consistent improvement and reduced performance variability. Landscape-condition analysis further showed that this fusion was most beneficial in shrubland-dominated and high-terrain fires, where it achieved the highest overall mean accuracy and the fewest failures. In contrast, its benefits were less reliable in evergreen forests, mixed vegetation, and low- to moderate-elevation terrain. In operational settings, the Sentinel-2 + Landsat-9 configuration offers a practical solution for post-fire recovery planning, erosion-risk assessment, watershed management, and ecological monitoring when field observations are available and timely satellite-based information is needed.

Landsat↗

In-situ observations of cyclic deformation in an extruded Mg-2Nd-1Y-0.1Zr-0.1Ca alloy

In this study, the evolution of deformation mechanisms during cyclic loading in an extruded, solution-treated Mg–2Nd–1Y–0.1Zr–0.1Ca alloy was investigated using a combination of in-situ loading, scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), and focused ion beam (FIB) nanofabrication. The initial microstructure exhibited a random crystallographic texture with no preferred grain orientation. Flat, rectangular dog-bone specimens were subjected to load-controlled, fully reversed fatigue for 50 cycles, during which the same region was sequentially mapped to track microstructural changes. After 10 cycles of loading deformation twins were observed. During tensile reloading detwinning or narrowing of those twinned regions occurred. After 20 cycles, detwinning ceased and residual twins remained in the material. SEM imaging revealed numerous surface slip traces after cyclic loading. EBSD-assisted slip trace analysis identified the activation of prismatic and pyramidal < c+a> slip systems during low-cycle fatigue. Site-specific scanning transmission electron microscopy (STEM) further revealed that deformation was also accommodated by basal < a> slip and the dissociation of < c+a> dislocations. Center-of-symmetry (COS) analysis confirmed that the dissociation of < c+a> dislocations resulted in the formation of I₁ intrinsic stacking faults after cyclic loading. These findings provide new insights into the complex interplay of dislocation mechanisms governing fatigue deformation in rare-earth-containing Mg alloys.

Cyclic deformation↗

Insights from FEED studies for retrofitting existing fossil power plants with carbon capture technology

Recent United States Department of Energy (DOE) sponsored front-end engineering design (FEED) studies for retrofitting existing fossil-fueled power plants with state-of-the-art carbon capture technology contain previously overlooked real-world design considerations for near-term deployment of carbon capture. Insights from examining seven recently published FEED study reports are summarized in this paper. This includes a discussion of the design, performance, and cost implications associated with (1) location-specific considerations such as water availability, land availability, and accessibility; (2) host-plant-specific factors such as flue gas specifications, allowable degree of integration between the capture system and host plant, and operational mode; and (3) miscellaneous factors such as market conditions, permitting requirements, and business case incentives. In conclusion, this manuscript highlights (1) water availability as a key design and cost driver, with host plant steam extraction increasing capture system cooling water availability, (2) modularization and constructability impacts on the number of capture trains, (3) the impacts of host plant operational mode and capacity factor on the business case for installing capture, and (4) the merit of continued research, development, and demonstration efforts addressing steam extraction, host plant tie-in at the stack, solvent reclamation and air emissions control.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗

Grid-responsive hydrogen production: Capital utilization and current density vs. efficiency in variable electricity markets

To achieve low-cost hydrogen production from water electrolyzers, grid tied electrolysis may need to operate dynamically to minimize the cost of supplying energy to the electrolyzer stack and produce hydrogen during low-cost hours and turn off/down during high-cost hours. Operating systems in this way can decrease capital utilization (capacity factor) and electricity costs. This strategy would shift the dominant cost drivers away from electricity (and thus efficiency) to the capital costs of the system, due to the underutilized capital when operating at low-capacity factors. Increasing the operational current density of the system could, in effect, reduce the capital cost of the system while producing hydrogen at a lower efficiency on a per unit energy basis. In the variable electricity cost profiles analyzed in this paper, increasing the current density for liquid alkaline from 0.5 A/cm2 to 1.5 Ac/m2 and proton exchange membrane electrolyzers from 2 A/cm2 to 4 A/cm2 resulted in substantial reductions in the levelized cost of hydrogen. Additionally, as capacity factors and electricity costs decrease, the optimal operating current density of the electrolyzer systems analyzed increases. These findings suggest R&D efforts should focus on increasing the operational current densities, reducing the turn down ratios, and understanding the durability implications of those strategies on low-temperature liquid alkaline and proton exchange membrane electrolyzers.

08 HYDROGEN↗

High-Throughput Uniformity and Defect Monitoring in Low-Temperature Electrolysis Porous Transport Layers Using X-Ray Radiography

Effective quality control (QC) for manufacturing proton exchange membrane water electrolysis (PEMWE) components is critical to enabling widespread adoption of the technology for hydrogen generation. This study investigates X-ray radiography as a novel, high-throughput, potentially in-line QC technique for detecting defects and assessing material property distributions in titanium-based porous transport layers (PTLs) which constitute a crucial component of low temperature PEMWE stacks. We obtain radiographs of a set of fifteen PTLs and model their absorbance of the broadband radiation as a second-order polynomial to account for the non-monoenergetic radiation source used in this study. The resulting model serves as a basis for predicting the areal density and porosity distributions of the PTLs. We find radiography successful in detecting multiple instances of defects, including holes/depressions, cracks, and excess material on the surface or in the pores of the material, demonstrating its potential as a robust in-line QC tool for PTL manufacturing.

08 HYDROGEN↗

In planta production of the nylon precursor beta-ketoadipate

Beta-ketoadipate (βKA) is an intermediate of the βKA pathway involved in the degradation of aromatic compounds in several bacteria and fungi. Beta-ketoadipate also represents a promising chemical for the manufacturing of performance-advantaged nylons. We established a strategy for the in planta synthesis of βKA via manipulation of the shikimate pathway and the expression of bacterial enzymes from the βKA pathway. Using Nicotiana benthamiana as a transient expression system, we demonstrated the efficient conversion of protocatechuate (PCA) to βKA when plastid-targeted bacterial-derived PCA 3,4-dioxygenase (PcaHG) and 3-carboxy-cis,cis-muconate cycloisomerase (PcaB) were co-expressed with 3-deoxy-D-arabinoheptulosonate 7-phosphate synthase (AroG) and 3-dehydroshikimate dehydratase (QsuB). This metabolic pathway was reconstituted in Arabidopsis by introducing a construct (pAtβKA) with stacked pcaG, pcaH, and pcaB genes into a PCA-overproducing genetic background that expresses AroG and QsuB (referred as QsuB-2). The resulting QsuB-2 x pAtβKA stable lines displayed βKA titers as high as 0.25 % on a dry weight basis in stems, along with a drastic reduction in lignin content and improvement of biomass saccharification efficiency compared to wild-type controls, and without any significant reduction in biomass yields. Using biomass sorghum as a potential crop for large-scale βKA production, techno-economic analysis indicated that βKA accumulated at titers of 0.25 % and 4 % on a dry weight basis could be competitively priced in the range of $2.04-34.49/kg and $0.47-2.12/kg, respectively, depending on the selling price of the residual biomass recovered after βKA extraction. This study lays the foundation for a more environmentally-friendly synthesis of βKA using plants as production hosts.

Kazaz, Sami↗

Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels

Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

42 ENGINEERING↗