SEARCH · Search NASA

Results for “Parallel”

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.

SourceResearch connectionFrom yearThrough year

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

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Physicochemical control of solvation and molecular assembly of charged amphiphilic oligomers at air-aqueous interfaces

Hypothesis : Understanding the rules that control the assembly of nanostructured soft materials at interfaces is central to many applications. We hypothesize that electrolytes can be used to alter the hydration shell of amphiphilic oligomers at the air-aqueous interface of Langmuir films, thereby providing a means to control the formation of emergent nanostructures. Experiments : Three representative salts – (NaF, NaCl, NaSCN) were studied for mediating the self-assembly of oligodimethylsiloxane methylimidazolium (ODMS-MIM + ) amphiphiles in Langmuir films. The effects of the different salts on the nanostructure assembly of these films were probed using vibrational sum frequency generation (SFG) spectroscopy and Langmuir trough techniques. Experimental data were supported by atomistic molecular dynamic simulations. Findings : Langmuir trough surface pressure – area isotherms suggested a surprising effect on oligomer assembly, whereby the presence of anions affects the stability of the interfacial layer irrespective of their surface propensities. In contrast, SFG results implied a strong anion effect that parallels the surface activity of anions. These seemingly contradictory trends are explained by anion driven tail dehydration resulting in increasingly heterogeneous systems with entangled ODMS tails and appreciable anion penetration into the complex interfacial layer comprised of headgroups, tails, and interfacial water molecules. These findings provide physical and chemical insight for tuning a wide range of interfacial assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

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↗

DOE OSTI2024

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ $\ggg$ 1 (with $\mathcal{X}$∥, $ \mathcal{X}$⟂ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. Further, the implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ~ (NΔt) 1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ = 10 10 ) .

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Hardware acceleration for HPS algorithms in two and three dimensions

We provide a flexible, open-source framework for hardware acceleration, namely massively-parallel execution on general-purpose graphics processing units (GPUs), applied to the hierarchical Poincaré–Steklov (HPS) family of algorithms for building fast direct solvers for linear elliptic partial differential equations. To take full advantage of the power of hardware acceleration, we propose two variants of HPS algorithms to improve performance on two- and three-dimensional problems. In the two-dimensional setting, we introduce a novel recomputation strategy that minimizes costly data transfers to and from the GPU; in three dimensions, we modify and extend the adaptive discretization technique of Geldermans and Gillman [1] to greatly reduce peak memory usage. We provide an open-source implementation of these methods written in JAX, a high-level accelerated linear algebra package, which allows for the first integration of a high-order fast direct solver with automatic differentiation tools. We conclude with extensive numerical examples showing our methods are fast and accurate on two- and three-dimensional problems.

Fast direct solvers↗

DOE OSTI2026

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗

DOE OSTI2025

A novel approach to increase accuracy in remotely sensed evapotranspiration through basin water balance and flux tower constraints

Remote sensing-derived evapotranspiration (RSET) products capture the spatiotemporal variations of evapotranspiration (ET) from field to basin scales with unprecedented details. However, their accuracy varies across RSET estimation methods and diverse hydroclimate regions. While ET modeling efforts to account for biophysical processes and controlling parameters have made good progress in recent years, a parallel approach of integrating in-situ ET with RSET could reduce biases in RSET products. Basin water balance ET (WBET) and flux tower ET are widely applied to evaluate RSET accuracy, yet such ET measurements are rarely used for RSET bias corrections, especially for large area applications. To address this issue, we propose a novel approach: the water balance equivalence (WABE) method, which generates spatially continuous WBET for correcting biases in RSET products. The WABE method computes synthetic WBET by integrating observed WBET and flux tower-derived FLUXCOM ET, which fills the spatial gaps of observed WBET and generates a spatially continuous WBET dataset. Synthetic WBET (2002–2015 annual average) of eight-digit hydrologic unit code (HUC8) basins across the conterminous United States (CONUS), constituting 44 % (887 out of 2035 basins) of CONUS basins, was determined within 2.0 % (RMSE = 12 %) of observed WBET at CONUS and between 1–12 % (RMSE = 3–33 %) across 18 regions in CONUS. With WABE-based bias corrections, the overall annual bias of RSET decreased from 10 % (RMSE = 34 %) to 6 % (RMSE = 26 %) across 37 flux tower sites. The WABE method offers a new approach for RSET accuracy improvement and shows great promise for large area implementations with a potential to yield substantial benefits for building accurate basin water budgets and water management decisions.

Khand, Kul↗

DOE OSTI2025

The Elements of Life, Photosynthesis and Genomics

I am a Professor of Biochemistry, Biophysics and Structural Biology and Plant and Microbial Biology at the University of California in Berkeley. I was born and raised in India, emigrated to the United States to attend university, earning a B.S. in Molecular Biology and a Ph.D. in Biochemistry at the University of Wisconsin in Madison. Following post-doctoral studies with Lawrence Bogorad at Harvard University where I became interested in genetic control of trace element quotas, I joined the department of Chemistry and Biochemistry at UCLA. One of the first to appreciate essential trace metals as potential regulators of gene expression, I articulated the details of the nutritional Cu regulon in Chlamydomonas. In parallel, I used genetic approaches to discover the genes governing missing steps in tetrapyrrole metabolism, including the attachment of heme to apocytochromes in the thylakoid lumen and the factors catalyzing the formation of ring V in chlorophyll. After biochemistry and classical genetics, I embraced genomics, taking a leadership role on the Joint Genome Institute’s efforts on the Chlamydomonas genome and more recently, contributing to high quality assemblies of several genomes in the green algal radiation, and large transcriptomic and proteomic datasets — focusing on the diel metabolic cycle in synchronized cultures and acclimation to key environmental and nutritional stressors — that are well-used and appreciated by the community. Finally, a new venture in Berkeley is the promotion of Auxenochlorella protothecoides as the true “green yeast” and as a platform for engineering algae to produce useful bioproducts.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Generalized grain boundary constitutive description implemented in a strain-gradient large-strain FFT-based formulation: Application to nano-metallic laminates

Herein this paper presents a general treatment of grain boundary constitutive behavior in the context of strain-gradient (SG) plasticity, and its numerical implementation in a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. Two novel grain boundary constitutive equations are proposed, allowing for more accurate description of the Burgers vector flow at the grain boundary. The capabilities of the generalized SG-LS-EVPFFT formulation are illustrated for the case of kink-band formation during layer-parallel compression of nano-metallic laminates (NMLs), requiring consideration of the interaction between dislocations and interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Mechanical properties, strain hardening, and fracture behavior of ultrasonic additively manufactured Zircaloy-4 after low-temperature neutron irradiation

Ultrasonic additive manufacturing (UAM) is a solid-state, layer-by-layer advanced manufacturing process that has the potential to create custom spatially controlled composites with embedded wires and sensors for nuclear component manufacture. For this work, to assess the feasibility of using UAM for nuclear-relevant materials research, the technique was used to produce a 3.5-mm-thick Zircaloy-4 plate for irradiation testing. The UAM Zircaloy-4 specimens were irradiated in the High Flux Isotope Reactor at a target irradiation temperature of 117 °C to 2.9 displacements per atom (dpa) to assess differences in irradiation-hardening behavior as a function of alloy processing path. The UAM and reference baseplate (BP) materials increased in yield strength by 372±27 MPa and 346±21 MPa, respectively, and both suffered significant reductions in uniform and total elongation attributed to irradiation hardening at low-temperature. Although the materials had similar nanoscale defect structures, including nanoscale black dot/loop features and strain-induced dislocation channels, the UAM material’s processing-related defects resulted in accelerated strain localization and failure as demonstrated by lower post-irradiation uniform elongation of UAM specimens (0.5 %) compared to BP (1.5 %) material. The UAM material also showed considerable anisotropy in mechanical response due to crack propagation along weld boundaries, resulting in differences in strength & ductility when tested parallel and perpendicular to the prior UAM build orientation. Therefore, although the fundamental irradiation response of UAM-processed Zircaloy-4 was phenomenologically comparable to that of BP reference material, additional optimization of the UAM processing is needed to produce irradiation-resistant and nuclear-relevant materials.

Digital image correlation↗

DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗

DOE OSTI2025

Large language models for batteries

Large Language Models (LLMs) are advanced artificial intelligence systems capable of solving diverse tasks using language, reasoning, and external tools. Despite their growing deployment in academia and industry, their potential remains underexplored in battery research. This review presents a comprehensive overview of existing and emerging applications of LLMs in batterie field, addressing two critical questions: What can LLMs offer to support battery-related tasks, and how to develop more effective models for this purpose. We begin by outlining the principles of LLMs and criteria for selecting appropriate models and tools for battery research and development. We then explore their roles in text-mining, data interpretation, and the development of intelligent battery systems. In parallel, we discuss technical challenges, such as data standardizing and sharing, model evaluation, and tool integration. Lastly, we propose future research directions with short-, medium-, and long-term goals and highlight more broad perspectives for connecting experts and cross-disciplinary collaborations.

SoC↗

DOE OSTI2025

Numerical and experimental analysis of mechanically induced failure in electric vehicle battery modules

Mitigating thermal runaway and cell-to-cell propagation is essential for improving the safety of electric and hybrid vehicles. Enhancing digital twin capabilities to predict battery mechanical abuse is particularly critical for automotive and aerospace applications, where crashworthiness is a key concern. Understanding failure conditions and propagation in battery modules during mechanical abuse is complex due to interactions between structural deformation, heat transfer, electrochemical processes, exothermic reactions and mechanical fracture. While prior studies have focused on modeling cell-level behavior, extending these models to module or pack level is necessary for a system level understating of electric vehicle safety. This study develops coupled large deformation finite element models that simultaneously solve for electrochemistry, material failure, internal short circuit and thermal runaway propagation. The models account for mechanical and thermal interactions between lithium-ion cells and other battery components while the contact interfaces are evolving with time. Model-predicted voltage, temperature and force responses are compared with experimental data for validation. The results demonstrate that the approach captures key failure mechanisms, including thermal propagation through heat transfer, electrical propagation from short circuits in parallel-connected cells, and mechanical propagation via penetration and crack formation. These findings show that computational models are valuable tools for understanding battery module failure and providing insight that can reduce the need for extensive experimental testing.

25 ENERGY STORAGE↗

DOE OSTI2025

Understanding the deformation behavior of the γ rich transformative Fe 38.5 Mn 20 Co 20 Cr 15 Si 5 Cu 1.5 complex concentrated alloy using in situ synchrotron diffraction

In situ tensile testing coupled with synchrotron x-ray diffraction was used to study the deformation behavior of metastability-engineered Fe 38.5 Mn 20 Co 20 Cr 15 Si 5 Cu 1.5 complex concentrated alloy. Monitoring the evolution of phase fraction and strain hardening response allowed the determination of true critical stress for the onset of the transformation-induced plasticity (TRIP) to be ∼375 MPa, preceded by slip starting at ∼255 MPa. In situ EBSD was used to validate the critical stress for transformation at the microstructure level and observe slip traces to confirm prior slip activity before the transformation. Further, a modeling framework based on stacking fault energy (SFE) was developed to predict the critical stress for transformation. Modeling suggested the SFE of the alloy to fall nearly 15 mJ/m 2 , which agrees well with SFE values calculated using synchrotron peak shifting (12 mJ/m 2 ) and thermodynamic calculation (11 mJ/m 2 ). As a result of γ-fcc to ε-hcp phase transformation, new {0002} ε planes emerged parallel to unaligned {111} γ planes with the tensile loading following S-N orientation relationship. Such selective emergence of new diffraction rings corresponding to ε phase is understood based on the reorientation of γ crystals with reference to tensile axis. In conclusion, this approach can be extended to effectively design alloys based on critical stress required for activating different deformation mechanisms to further push the limits of the strength-ductility envelope.

Complex concentrated alloy↗

DOE OSTI2025

Effects of magnetic field assisted heat treatment on the microstructure and mechanical properties of Fe-0.63 %C alloy

This study investigates the influence of an applied magnetic field on the microstructural evolution and mechanical properties of hypoeutectoid steels subjected to heat treatment. Tensile tests and microstructural analysis were performed on samples processed under varying magnetic field strengths (0 T, 5 T, and 9 T) and different austenitization incubation times. The results indicate that the application of a magnetic field alters the fraction of proeutectoid ferrite phase without changing the cooling rates and heat treatment process. Additionally, pearlite microstructural features such as lamellar spacing and misorientation angles exhibit variations under different field strengths. While the pearlite nodule diameter remains largely unaffected, an increase in percentage elongation and strength is observed in the 5 T treated sample, attributed to changes in microstructural features with the magnetic field. Additionally, the percentage elongation is reduced in the samples heat treated with reduced austenitization incubation times. The study further demonstrates that low-angle misorientations increased in the samples taken parallel to the magnetic field direction, influencing the mechanical response. These findings suggest that applying a magnetic field during heat treatment provides an additional driving force for phase transformations, offering a manufacturing process for tailoring microstructures and optimizing mechanical properties. Moreover, integrating magnetic fields in heat treatment processes has potential benefits in energy efficiency.

High magnetic field↗

DOE OSTI2024

Solidification and crystallographic texture modeling of laser powder bed fusion Ti-6Al-4V using finite difference-monte carlo method

Laser powder bed fusion (LPBF) additive manufacturing makes near-net-shaped parts with reduced material cost and time, rising as a promising technology to fabricate Ti-6Al-4V, a widely used titanium alloy in aerospace and medical industries. However, LPBF Ti-6Al-4V parts produced with 67° rotation between layers, a scan strategy commonly used to reduce microstructure and property inhomogeneity, have varying grain morphologies and weak crystallographic textures that change depending on processing parameters. Here, this study predicts LPBF Ti-6Al-4V solidification at three energy levels using a finite difference-Monte Carlo method and validates the simulations with large-area electron backscatter diffraction (EBSD) scans. The developed model accurately shows that a <001> texture forms at low energy and a <111> texture occurs at higher energies parallel to the build direction but with a lower strength than the textures observed from EBSD. A validated and well-established method of combining spatial correlation and general spherical harmonics representation of texture is developed to calculate a difference score between simulations and experiments. The quantitative comparison enables effective fine-tuning of nucleation density (N 0 ) input, which shows a nonlinear relationship with increasing energy level. Future improvements in texture prediction code and a more comprehensive study of N 0 with different energy levels will further advance the optimization of LPBF Ti-6Al-4V components. These developments contribute a novel understanding of crystallographic texture formation in LPBF Ti-6Al-4V, the development of robust model validation and calibration pipeline methodologies, and provide a platform for mechanical property prediction and process parameter optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗