Engineering Papers⌕ Search

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.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗
DOE OSTI2022

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

We report molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. The SLIC/CDC multiscale model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). The SLIC/CDC model predicts Gibbs energies of solvation for a database of 500 solutes in water with average accuracy better than 1 kcal/mol both for experimental measurements and for explicit-solvent molecular dynamics simulations. The separate SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent accessible surface area (SASA) nonpolar energy mode. The SLIC/SASA model predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems and better than 1.6 kcal/mol average accuracy in ionic liquids. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.

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

Modeling the Environment-Dependent Kinetics of Oxygen Reduction Reaction – a Continuum Model for Electric Double Layer

Here, for proton-exchange-membrane fuel cells (PEMFCs) to achieve broad commercialization, improved energy-conversion efficiency with minimal Pt-based electrocatalyst is required. Because the sluggish rate of oxygen reduction reaction (ORR) limits the efficiency of PEMFCs, the efficiency improvement requires a better understanding of ORR kinetics and mechanism to design better catalyst. To understand the ORR mechanism, theoretical and experimental analyses have been conducted. While previous studies reasonably explained the catalyst-dependent activity on single crystal catalysts in 0.1 M perchloric acid solution, the explicit effect of electrolyte and related microenvironments is not thoroughly understood. The change in the electrolyte alters the electric-double-layer (EDL) structure and thus the local microenvironment at the electrode/electrolyte interface. Thus, the structure of the EDL should be carefully analyzed to uncover the electrolyte-dependent reaction kinetics. In this talk, we propose a multiscale continuum model to predict the EDL structure and examine the effect of perchloric acid concentration on ORR activity on Pt (111). The model includes Density Potential Functional Theory (DPFT) for electron density and Modified Poisson Boltzmann equation for species’ density and electric potential. Also, the interaction between adsorbents and electric field is taken into account by minimizing the grand potential. After model validation with experimentally measured double-layer capacity data as a function of applied potential and concentration, the effect of the perchloric acid concentration (0.02 M – 0.2 M) on ORR activity is analyzed and discussed. It is shown that the model reproduces the specific activity obtained in the experiments when assuming the oxygen adsorption is limiting the rate, which can be attributed to the large energetic barrier for solvent reorganization. Then, extension of the model to PEMFC ionomer electrolytes will be introduced. Overall, the model framework and findings provide insights into the ORR mechanism and guidance on how to tailor catalyst materials for increased PEMFC performance.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Progress on continuum modeling of lithium–sulfur batteries

While lithium–sulfur batteries are a promising next-generation chemistry devices due to their high theoretical energy density, commercialization has been slow due to low coulombic efficiency and poor cycle life. This review explores the ways in which continuum modeling contributes to the understanding of lithium–sulfur (LiS) battery mechanisms and cell-level performance through the lens of micro- and macroscale phenomena. We examine different approaches to modeling important physical phenomena such as reaction mechanisms, cathode microstructure, shuttling, nucleation and precipitation, and transport limitations. This paper also emphasizes the significance and challenge of connecting typical modeling parameters and assumptions to systems-level metrics of a standard state-of-art high performing lithium–sulfur cell. Particularly important, the considerations for high energy density cells and the areas where continuum models can facilitate better collaboration are discussed. Here, we also summarize a few selected works to highlight experimentally-driven modeling, use of electroanalytical techniques, and parameter identification approaches to enable model-based design and advanced battery management systems.

25 ENERGY STORAGE↗
DOE OSTI2023

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development of a concurrent coupled Atomistic - Continuum model to predict the defect and grain structure for Additive Manufacturing process [Slides]

Development of a coupled atomistic-continuum model for metal AM process. Two-way coupling between atomistic (heat flux) and continuum (temperature) domain. Identification of grain and dislocation structure at the atomistic domain. Conversion of dislocation data from discrete atomistic to density form at the continuum. Obtain the spatial distribution of grain structure and dislocation density in additively manufactured metal polycrystal material.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
DOE OSTI2020

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗
DOE OSTI2024

Physics-Based Continuum Modeling for an Aqueous Rechargeable Zn/MnO 2 Battery

This study introduces a framework for modeling the aqueous Zn/MnO 2 rechargeable battery. A reaction system and a physics-based continuum model are proposed based on two reaction types, one involving insertion and the second related to dissolution and deposition of a solid reaction product. The model, fitted to empirical data, predicts voltage behavior and capacity limitations during cycling, identifying electrolytic zinc depletion as a limiting mechanism, depending on the original cell construction. The research suggests the need for further material characterization and reaction analysis, which will advance our understanding and facilitate the development of grid-scale energy storage solutions.

Electrochemistry↗
DOE OSTI2023

Continuum Model Development for Electrochemical Systems (Statement of Work: Ilenia Battiato Subcontract)

The subcontractor shall provide the services of qualified multiscale modelers to perform tasks that contribute to reaching the objectives of the LDRD project entitled Automated and Accelerated Continuum Model Development for Electrochemical Systems (Tracking number 24-ERD-051). These tasks relate to the development, deployment, and validation of multiscale models relevant to electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Continuum Model Development for Flow and Transport in Electrochemical Systems

Stanford University (Subcontractor) shall provide the services of qualified multiscale modelers to perform tasks that contribute to reaching the objectives of the LDRD project entitled Automated and Accelerated Continuum Model Development for Electrochemical Systems (Tracking number 24-ERD-051). These tasks relate to the development, deployment, and validation of multiscale models relevant to flow and transport in electrochemical systems.

42 ENGINEERING↗
DOE OSTI2021

Lattice collective modes from a continuum model of magic-angle twisted bilayer graphene

Here, we show that the insulating states of magic-angle twisted bilayer graphene support a series of collective modes corresponding to local particle-hole excitations on triangular lattice sites. Our theory is based on a continuum model of the magic angle flat bands. When the system is insulating at moiré band filling ν = -3, our calculations show that the ground state supports seven low-energy modes that lie well below the charge gap throughout the moiré Brillouin zone, one of which couples strongly to THz photons. The low-energy collective modes are faithfully described by a model with a local SU(8) degree of freedom in each moiré unit cell that we identify as the direct product of spin, valley, and an orbital pseudospin. Apart from spin and valley-wave modes, the collective mode spectrum includes a low-energy intraflavor exciton mode associated with transitions between flat valence and conduction band orbitals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Continuum Modeling of Porous Electrodes for Electrochemical Synthesis

Electrochemical synthesis possesses substantial promise to utilize renewable energy sources to power the conversion of abundant feedstocks to value-added commodity chemicals and fuels. Of the potential system architectures for these processes, only systems employing 3-D structured porous electrodes have the capacity to achieve the high rates of conversion necessary for industrial scale. However, the phenomena and environments in these systems are not well understood and are challenging to probe experimentally. Fortunately, continuum modeling is well-suited to rationalize the observed behavior in electrochemical synthesis, as well as to ultimately provide recommendations for guiding the design of next-generation devices and components. In this review, we begin by presenting an historical review of modeling of porous electrode systems, with the aim of showing how past knowledge of macroscale modeling can contribute to the rising challenge of electrochemical synthesis. We then present a detailed overview of the governing physics and assumptions required to simulate porous electrode systems for electrochemical synthesis. Here, leveraging the developed understanding of porous-electrode theory, we survey and discuss the present literature reports on simulating multiscale phenomena in porous electrodes in order to demonstrate their relevance to understanding and improving the performance of devices for electrochemical synthesis. Lastly, we provide our perspectives regarding future directions in the development of models that can most accurately describe and predict the performance of such devices and discuss the best potential applications of future models.

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

Eulerian finite element implementations of a dislocation density-based continuum model

In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A data-driven peridynamic continuum model for upscaling molecular dynamics

Nonlocal models, including peridynamics, often use integral operators that embed lengthscales 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 Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This 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. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.

homogenization↗
DOE OSTI2023

Eulerian finite element simulations of the drop weight impact test with a dislocation Density-based continuum model

During a drop weight impact test, the kinetic energy of the falling weight is transferred into the sample resting on the anvil. The plastic deformation in the sample is an important mechanism for the dissipation of the input kinetic energy. We use Eulerian finite element analysis to simulate the deformation and temperature evolution in a 1,3,5-trinitro-1,3,5-triazine (RDX) sample consisting of multiple crystals. In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, i.e., advection. In an effort to reduce the advection errors, we use a rate form of a dislocation density-based continuum model by Luscher et al. Here, the simulations predict localization of plastic deformation, and plastic dissipation as a significant source of heat generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

Long range-corrected (LRC) or range-separated hybrid (RSH) functionals where the long-range (LR) limit of electronic interactions is set to the exact exchange have been shown to correct the tendency of traditional density functional theory (DFT) to underestimate the frontier orbital gap. Consequently, the use of such functionals in calculating electronic excited states using linear response based time-dependent DFT (TDDFT) has been successful in correcting the tendency for underestimating the energies of charge transfer states by DFT-based calculations. More recently formulations of functionals that attenuate the LR limit to address condensed-phase effects to polarize the electronic density have been reported. In particular screened RSH (SRSH) combined with polarizable continuum model (PCM) was benchmarked successfully in reproducing the fundamental gap and charge transfer state energies of molecular systems in the condensed phase. Here we use SRSH-PCM to address triplet excited states, and show its success in obtaining correspondence of the low-lying triplet states to the singlet–triplet gap in a similar way that the fundamental orbital gap corresponds to electron removal and addition energies. Importantly, the accuracy of the SRSH-PCM in calculating triplet excitations stands on the polarization consistent framework in addressing the scalar dielectric constant and without affecting the optimal tuning by triplet energies. Finally, the prospect of even further improving the SRSH-PCM accuracy in calculating triplet states can be achieved by optimal tuning on the basis of the spin multiplicity gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗