SEARCH · Engineering Papers

Results for “continuum modeling”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Long-range hydrogen-binding effects of carbide interfaces in iron

A micromechanics model was developed to evaluate the elastic binding energy between carbide precipitates and hydrogen interstitials using Eshelby's equivalent inclusion method. In this work, density functional theory (DFT) simulations were performed to obtain the material-specific quantities, e.g., lattice constants and the elastic constants, for the continuum model. Using this model, we find that for coherent carbide precipitates, hydrogen atoms are more likely to bind on the broad surfaces of the disk-like precipitates, which is consistent with experimental observations. For semicoherent and incoherent precipitates, our model suggests that it is possible for semicoherent precipitates to have significant hydrogen binding capability while there is no hydrogen-binding capability of incoherent precipitates, which also agrees with experimental findings. In addition, several factors that influence the binding energies between hydrogen atoms and carbide precipitates were quantitatively analyzed, including the precipitate size, morphology, orientation, and interface. These collective results include both the position and the value of the strongest hydrogen-binding interaction for a wide range of carbide stoichiometries, which contributes to our understanding of hydrogen trapping in steel-based materials.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets—Part I: Mesoscale simulations

Recently a mesoscale model was developed that includes thermomechanically consistent constitutive equations with dissipation due to porous compaction, inelastic distortional deformation rate, and slipping on persistent dry joint surfaces. In this Part I of a two-part paper, we describe enhancements of the mesoscale model which include: normal compliance with hysteresis, non-persistent joints with finite area, joint sets with angular deviation, and the weakening effect of fluid pressure in the joints. Examples study the influences of these effects on average strength. These computationally expensive simulations are used to inspire functional forms in a continuum model, described in Part II, which is needed for large scale simulations of real applications.

42 ENGINEERING↗

DOE OSTI2025

Cohesive instability in elastomers: insights from a crosslinked Van der Waals fluid model

Abstract The resistance to volumetric deformations displayed by polymer networks is largely due to secondary and tertiary interactions between neighboring polymer chains. These interactions are both entropic and enthalpic in nature but are fundamentally different from the entropic forces that resist shearing in these networks. In this paper, we introduce a new depiction of elastomers as a crosslinked Van der Waals fluid. Starting from first principles, we develop constitutive equations that are implemented in a continuum model as well as a discrete network model. Our models predict that the failure of polymer networks may be driven by an instability in the underlying polymer bulk ‘fluid’ or by the breaking of polymer chains, depending on the loading path taken. The results of this study indicate that material failure in elastomers exposed to a purely triaxial state, such as in a poker chip experiment, may be driven by an entirely different mode of instability than those deformed in pure shear, such as in a uniaxial tension experiment.

Lamont, Samuel C.↗

DOE OSTI2020

A mathematical model of solid-state dewetting of barium thin films on W (112)

Solid state dewetting occurs when a thin solid film is heated. The temperature of dewetting depends on the thickness of the film; dewetting can be observed in the range of 1/3 to 1/2 of the bulk melting temperature. While remaining solid, the film behaves in a manner similar to liquids dewetting and agglomerating to forming islands or droplets. One of the possible mechanisms is the conversion of a metastable thin film geometry into a more stable form. Heating the metastable film gives the film atoms higher mobility, and the films retract, dewetting the surface. This atomic motion can be restricted due to surface anisotropy. Herein we present in situ emission microscopy observations of barium thin films deposited onto W (112) by thermal evaporation. From the modeling viewpoint, the evolution of the film in this system could be divided in four stages: (i) the nucleation and growth of the thin film as a simply connected region; (ii) formation of droplets/islands/hillocks; (iii) nucleation of holes; (iv) evolution of the components of the disconnected film to their equilibrium state. We present a continuum model that is qualitatively consistent with the evolution of the film observed at the initial stage of the experiment and discuss the later stages of the evolution of surface structures.

36 MATERIALS SCIENCE↗

DOE OSTI2023

In Situ Raman Mapping of Si Island Electrodes and Stress Modeling as a Function of Lithiation and Size

Si is known for cracking and delamination during electrochemical cycling of a battery due to the large volume change associated with Li insertion and extraction. However, it has been found experimentally that patterned Si island electrodes that are 200 nm thick and less than 7 μm wide can deform in a purely elastic manner. Inspired by this, we performed in situ Raman stress characterization of model poly-crystalline Si island electrodes using an electrochemical cell coupled with an immersion objective lens and designed for a short working distance. A 5 μm wide Si island electrode showed a parabolic stress profile during lithiation, while for a 15 μm Si island electrode, a stress plateau in the center of the electrode was observed. A continuum model with coupled electro-chemo-mechanical (ECM) physics was established to understand the stress measurement. A qualitative agreement was reached between modeling and experimental data, and the critical size effect could be explained by the Li diffusive flux as governed by competition between the Li concentration and hydrostatic stress gradients. Below the critical size, the stress gradient drives Li toward the edges, where the electrode volume is free to expand, while above the critical size, the stress plateau inhibits Li diffusion to the edge and forces destructive stress relief by cracking. Furthermore, this work represents a promising methodology for in situ characterization of ECM coupling in battery electrodes, with suggestions provided for further improvement.

25 ENERGY STORAGE↗

DOE OSTI2020

Data-driven learning of nonlocal models: from high-fidelity simulations to constitutive laws

We show that machine learning can improve the accuracy of simulations of stress waves in one-dimensional composite materials. We propose a data-driven technique to learn nonlocal constitutive laws for stress wave propagation models. The method is an optimization-based technique in which the nonlocal kernel function is approximated via Bernstein polynomials. The kernel, including both its functional form and parameters, is derived so that when used in a nonlocal solver, it generates solutions that closely match high-fidelity data. The optimal kernel therefore acts as a homogenized nonlocal continuum model that accurately reproduces wave motion in a smaller-scale, more detailed model that can include multiple materials. We apply this technique to wave propagation within a heterogeneous bar with a periodic microstructure. Several one-dimensional numerical tests illustrate the accuracy of our algorithm. The optimal kernel is demonstrated to reproduce high-fidelity data for a composite material in applications that are substantially different from the problems used as training data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Modeling the Effect of Material Properties on Liquid-Alkaline Water Electrolysis

Liquid-alkaline water electrolyzers (LAWEs) use electricity to drive the conversion of water to H 2 and O 2 gas. These devices benefit from the use of low-cost nickel electrodes and metal-oxide separators, but suffer from lower current densities and higher cell voltages than proton-exchange-membrane water electrolyzers. Identifying the inefficiencies that result in this poor performance is key to mitigating losses and optimizing LAWEs. Here, we report an experimentally-validated 1-D continuum model of a LAWE that elucidates the gradients within the cell, simulates H 2 crossover, and projects the energy improvements made possible by modulating the properties of the electrodes and separator. The model captures the Nernstian polarization losses and the distribution of gas- and liquid-phases within the electrodes, enabling quantification of energy losses associated with kinetic, ohmic, and bubble-induced (mass-transport) resistances. Simulations demonstrate that LAWE can achieve energy intensities of 50 kWh kg −1 of H 2 at 1 A cm −2 using improved electrode and separator properties.

08 HYDROGEN↗

DOE OSTI2025

QM Investigation of Rare Earth Ion Interactions with First Hydration Shell Waters and Protein-Based Coordination Models

Here, conventional methods for extracting rare earth metals (REMs) from mined mineral ores are inefficient, expensive, and environmentally damaging. Recent discovery of lanmodulin (LanM), a protein that coordinates REMs with high-affinity and selectivity over competing ions, provides inspiration for new REM refinement methods. Here, we used quantum mechanical (QM) methods to investigate trivalent lanthanide cation (Ln 3+ ) interactions with coordination systems representing bulk solvent water and protein binding sites. Energy decomposition analysis (EDA) showed differences in the energetic components of Ln 3+ interaction with representatives of solvent (water, H 2 O) and protein binding sites (acetate, CH 3 COO – ), highlighting the importance of accurate description of electrostatics and polarization in computational modeling of REM interactions with biological and bioinspired molecules. Relative binding free energies were obtained for Ln 3+ with coordination complexes originating from binding sites in PDB structures of a lanthanum binding peptide (PDB entry 7CCO) and LanM, with explicit consideration of the first hydration shell waters, according to quasi-chemical theory (QCT). Beyond the first shell, the bulk solvent environment was represented with an implicit continuum model. Ln 3+ interactions with (H 2 O) 9 and both binding site models became more favorable, moving down the periodic series. This trend was more pronounced with the protein binding site models than with water, resulting in affinity increasing with periodic number, except for the last REM, Lu 3+ , which bound less favorably than the preceding element, Yb 3+ . Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. Conversely, the previously reported experimental data for LanM show a preference for the earlier lanthanides; this is likely due to longer-range interactions and cooperative effects, which are not represented by the reduced models. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. In contrast to the previously reported experimental data for LanM, the peptide preferentially binds the earlier lanthanides. This difference likely arises due to longer-range interactions and cooperative effects not represented by the peptide. Further investigation of Ln 3+ interactions with whole proteins using polarizable molecular mechanics models with explicit solvent is warranted to understand the influence of longer-ranged interactions, cooperativity, and bulk solvent. Nevertheless, the present work provides new insights into Ln 3+ interactions with biomolecules and presents an effective computational platform for designing specific single-site REM binding peptides more efficiently.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Effects of lithium intercalation in twisted bilayer graphene

In this work, we investigate the effects of lithium intercalation in twisted bilayers of graphene, using first-principles electronic structure calculations. To model this system we employ commensurate supercells that correspond to twist angles of 7.34° and 2.45°. From the energetics of lithium absorption we demonstrate that for low Li concentration the intercalants cluster in the AA regions with double the density of a uniform distribution. The charge donated by the Li atoms to the graphene layers results in modifications to the band structure that can be qualitatively captured using a continuum model with modified interlayer couplings in a region of parameter space that has yet to be explored either experimentally or theoretically. Thus, the combination of intercalation and twisted layers simultaneously provides the means for spatial control over material properties and an additional knob with which to tune moiré physics in twisted bilayers of graphene, with potential applications ranging from energy storage and conversion to quantum information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Cristae formation is a mechanical buckling event controlled by the inner mitochondrial membrane lipidome

Cristae are high–curvature structures in the inner mitochondrial membrane (IMM) that are crucial for ATP production. While cristae–shaping proteins have been defined, analogous lipid–based mechanisms have yet to be elucidated. Here, we combine experimental lipidome dissection with multi–scale modeling to investigate how lipid interactions dictate IMM morphology and ATP generation. When modulating phospholipid (PL) saturation in engineered yeast strains, we observed a surprisingly abrupt breakpoint in IMM topology driven by a continuous loss of ATP synthase organization at cristae ridges. We found that cardiolipin (CL) specifically buffers the inner mitochondrial membrane against curvature loss, an effect that is independent of ATP synthase dimerization. To explain this interaction, we developed a continuum model for cristae tubule formation that integrates both lipid and protein–mediated curvatures. This model highlighted a snapthrough instability, which drives IMM collapse upon small changes in membrane properties. We also showed that cardiolipin is essential in low–oxygen conditions that promote PL saturation. These results demonstrate that the mechanical function of cardiolipin is dependent on the surrounding lipid and protein components of the IMM.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Low-energy moiré phonons in twisted bilayer van der Waals heterostructures

We develop a low-energy continuum model for phonons in twisted moiré bilayers, based on a configuration-space approach. In this approach, interatomic force constants are obtained from density functional theory (DFT) calculations of untwisted bilayers with various in-plane shifts. This allows for efficient computation of phonon properties for any small twist angle, while maintaining DFT-level accuracy. Based on this framework, we show how the low-energy phonon modes, including interlayer shearing and layer-breathing modes, vary with the twist angle. As the twist angle decreases, the frequencies of the low-energy modes are reordered and the atomic displacement fields corresponding to phonon eigenmodes break translational symmetry, developing periodicity on the moiré length scale. Here, we demonstrate the capabilities of our model by calculating the phonon properties of three specific structures: bilayer graphene, bilayer molybdenum disulfide (MoS 2 ), and molybdenum diselenide-tungsten diselenide (MoSe 2 -WSe 2 ).

2-dimensional systems↗

DOE OSTI2021

Ab initio calculations of material properties for modeling debris (FY2021 Progress Report)

This report summarizes the work done during the first 6 months of the project “Ab initio calculations of material properties for modeling debris” that is part of the Rapid Response Research (R3) Venture. The main focus for this fiscal year (FY) has been the evaluation of classical interatomic potentials for the description of the condensation of uranium oxides in the regimes relevant to fireball conditions. We also delineate how these functions will be used within the molecular dynamics and Monte Carlo frameworks to compute thermodynamic properties of these materials. These descriptions are to be utilized within continuum modeling formulations to generate a more accurate characterization of the homogeneous nucleation process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Data driven learning of robust nonlocal models

Nonlocal models use integral operators that embed length-scales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials' behavior, we propose a learning framework to extract, from MD data, an optimal nonlocal model as a surrogate for MD displacements. Our framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. The efficacy of this approach is demonstrated with several numerical tests for single layer graphene both in the case of perfect crystal and in the presence of thermal noise.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Mineral-associated organic matter is heterogeneous and structured by hydrophobic, charged, and polar interactions

The formation of mineral-associated organic matter (MAOM) is a key phenomenon that may explain the slow turnover rates of carbon in soil organic matter (SOM). Despite this, important details pertaining to the structure and dynamics of MAOM remain unknown. In the present study, we use replica-exchange molecular dynamics simulations to gain insight into the structure of MAOM on the surface of prototypical phyllosilicate clay and Fe-oxide minerals, montmorillonite and goethite, fine-grained minerals that strongly impact soil carbon dynamics in temperate and tropical regions, respectively. We examine the impact of aqueous chemistry through the presence of either Na + or Ca 2 + charge balancing counterions. Our results are consistent with the hypothesized multilayer sorption (“onion-skin”) model of MAOM and help to explain previous observations regarding the patchy distribution of SOM on mineral surfaces. In particular, the SOM coatings are partial and laterally heterogeneous, and water retains extensive access to mineral surfaces even when significant SOM sorption occurs. Low molecular weight neutral SOM molecules ( < 200 Da) infrequently interact with the mineral surfaces nor their sorbed organic matter coatings and are increasingly labile with decreasing molecular weight. This observation is inconsistent with a central feature of the predominant soil continuum model of SOM and suggests that further iterations of the conceptual model may be required.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

First-principles Modeling and Design of Solid-State Interfaces for the Protection and Use of Lithium Metal Anodes

Li-ion batteries are one of the most advanced energy storage technologies in use today. Li-ion batteries are used in a multitude of applications ranging from consumer electronics, medical devices, sensors and grid storage. However, improving the capacity and energy density delivered by current Li-ion technology requires advanced materials research into novel chemical systems. In this project we have focused particularly in the use of solid-state electrolytes with lithium metal electrodes. Research into all solid-state batteries (ASSB) with Li metal electrodes has significantly expanded in recent years, however most studies reported experimental findings, which left substantial room for theoretical and modeling work as a tool to understand and determine design principles allowing reliable and safe use of ASSBs with Li metal. Among the remaining obstacles preventing reliable use of ASSBs with a Li metal electrode, the stability of the interface between the solid electrolyte and Li metal, and the propagation/dendrite formation of Li metal and resulting mechanical degradation of the electrolyte are key phenomenon that are yet to be fully understood. In the current project we have addressed these two coupled phenomena using first principles calculations and mesoscale continuum modeling. We have obtained chemical and electrochemical stability windows for several solid electrolyte materials. Additionally, from mathematical and numerical modeling of Li protrusion and dendrite initiation during plating and stripping we have determined design criteria in terms of chemical, electrochemical, and mechanical properties and operating conditions for which stable deposition can occur. We also considered the effects of mixed electronic-ionic conduction in solid electrolytes, which has more recently been suggested as another important mechanism involved in ASSB failure. Throughout our work we have successfully addressed important questions necessary for the use of ASSB’s. We have determined guiding principles for materials properties and operating conditions necessary to operate ASSB’s. And have proposed novel solid electrolyte materials with predicted chemical stability an ionic conductivity. Although this represents significant progress in our understanding, open questions remain in order to fully develop reliable and safely operate ASSBs with Li metal. Future work, building on this project will require further experimental, theoretical and simulation efforts to address remaining questions.

25 ENERGY STORAGE↗

DOE OSTI2025

Computationally Guided and Experimentally Validated Design of Custom Chelators for Critical Mineral Recovery

Selective, high throughput separation of target critical metals from complex environments such as fly ash leachates and mining process streams presents a significant challenge for economical production. Custom chelators and sorbents are an attractive technology for selective metal extraction, however it can be difficult to predict their performance, and significant experimental efforts are often required to develop chelating technologies. Here, we present a computational strategy focused on modelling chelator-metal binding interactions and benchmark these results versus experimental data. A computational pipeline combining forcefield, semiempirical, and meta-GGA methods with a thermodynamic framework optimized for error cancellation has been developed to predict binding energies of chelator complexes towards critical mineral recovery applications. This approach, originally validated on [2.2.2] cryptates binding mono- and divalent cations, demonstrated robust predictive capabilities with an R2 of 0.850 against experimental aqueous binding energies. The workflow includes metadynamics for exploring high-dimensional potential energy surfaces and a cluster-continuum model for accurate yet computationally efficient solvation modeling. Error cancellation between solvation energies of free and chelator-coordinated ions enables faster convergence, even with finite cluster sizes. Initial studies on the cryptates revealed consistent metal-ligand coordination patterns, with systematic variations influenced by ion size and charge, highlighting key structural features linked to binding selectivity. Further studies of a proprietary chelator have resulted in identification of previously unreported selectivity towards economically significant metals, which in-house experiments have confirmed, demonstrating the feasibility of this approach. By applying this methodology to new chelators targeting critical minerals such as lithium, cobalt, nickel and other strategic metals, we aim to accelerate the discovery of next-generation chelators for efficient recovery, recycling, and separation processes. This computational framework serves as the backbone of a high-throughput design pipeline tailored for sustainable resource utilization and may be applied to a wide range of systems to meet experimental needs.

computational materials↗

DOE OSTI2024

Analytical Model for Atomic Relaxation in Twisted Moiré Materials

By virtue of being atomically thin, the electronic properties of heterostructures built from two-dimensional materials are strongly influenced by atomic relaxation. The atomic layers behave as flexible membranes rather than rigid crystals. Here we develop an analytical theory of lattice relaxation in twisted moiré materials. We obtain analytical results for the lattice displacements and corresponding pseudo gauge fields, as a function of twist angle. We benchmark our results for twisted bilayer graphene and twisted WSe 2 bilayers using large-scale molecular dynamics simulations. Our single-parameter theory is valid in graphene bilayers for twist angles 𝜃 ≳ 0.7°, and in twisted WSe 2 for 𝜃 ≳ 1.6°. Furthermore, we also investigate how relaxation alters the electronic structure in twisted bilayer graphene, providing a simple extension to the continuum model to account for lattice relaxation.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗