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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Theoretical Calculation and Validation of the Water Vapor Continuum Absorption

The primary objective of this investigation is the development of an improved parameterization of the water vapor continuum absorption through the refinement and validation of our existing theoretical formalism. The chief advantage of our approach is the self-consistent, first principles, basis of the formalism which allows us to predict the frequency, temperature and pressure dependence of the continuum absorption as well as provide insights into the physical mechanisms responsible for the continuum absorption. Moreover, our approach is such that the calculated continuum absorption can be easily incorporated into satellite retrieval algorithms and climate models. Accurate determination of the water vapor continuum is essential for the next generation of retrieval algorithms which propose to use the combined constraints of multi-spectral measurements such as those under development for EOS data analysis (e.g., retrieval algorithms based on MODIS and AIRS measurements); current Pathfinder activities which seek to use the combined constraints of infrared and microwave (e.g., HIRS and MSU) measurements to improve temperature and water profile retrievals, and field campaigns which seek to reconcile spectrally-resolved and broad-band measurements such as those obtained as part of FIRE. Current widely used continuum treatments have been shown to produce spectrally dependent errors, with the magnitude of the error dependent on temperature and abundance which produces errors with a seasonal and latitude dependence. Translated into flux, current water vapor continuum parameterizations produce flux errors of order 10 W/ml, which compared to the 4 W/m' magnitude of the greenhouse gas forcing and the 1-2 W/m' estimated aerosol forcing is certainly climatologically significant and unacceptably large. While it is possible to tune the empirical formalisms, the paucity of laboratory measurements, especially at temperatures of interest for atmospheric applications, preclude tuning the empirical continuum models over the full spectral range of interest for remote sensing and climate applications. Thus, we propose to further develop and refine our existing far-wing formalism to provide an improved treatment applicable from the near-infrared through the microwave. Based on the results of this investigation, we will provide to the remote sensing/climate modeling community a practical and accurate tabulation of the continuum absorption covering the near-infrared through the microwave region of the spectrum for the range of temperatures and pressures of interest for atmospheric applications.

Ma, Qiancheng↗
DOE OSTI2023

Analysis of bipolar membranes for electrochemical CO 2 capture from air and oceanwater

Carbon dioxide (CO 2 ) must be removed from the atmosphere to mitigate the negative effects of climate change. However, the most scalable methods for removing CO 2 from the air require heat from fossil-fuel combustion to produce pure CO 2 and continuously regenerate the sorbent. Bipolar-membrane electrodialysis (BPM-ED) is a promising technology that uses renewable electricity to dissociate water into acid and base to regenerate bicarbonate-based CO 2 capture solutions, such as those used in chemical loops of direct-air-capture (DAC) processes, and in direct-ocean capture (DOC) to promote atmospheric CO 2 drawdown via decarbonization of the shallow ocean. In this study, we develop an experimentally validated 1D model for the electrochemical regeneration of CO 2 from bicarbonate-based carbon capture solutions and seawater using BPM-ED. For DAC, our experimental and computational results demonstrate that pH swings induced by BPM water dissociation drive the formation of CO 2 at the cation-exchange layer|catholyte interface with energy-intensities of less than 150 kJ mol -1 . However, high rates of bubble formation increase energy intensity at current densities >100 mA cm -2 . Correspondingly, accelerating water dissociation catalysis and enacting bubble removal could enable CO 2 recovery at energy intensities <100 kJ mol -1 and current densities >100 mA cm -2 . For DOC, mass transport limitations associated with low carbon concentrations in oceanwater suggest that DOC is best suited for clean production of acid and base usable in downstream processes. Further, these results provide design principles for industrial-scale CO 2 recovery using BPM-ED.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Enabling 6C Fast Charging of Li–Ion Batteries with Graphite/Hard Carbon Hybrid Anodes

Li-ion batteries that can simultaneously achieve high-energy density and fast charging are essential for electric vehicles. Graphite anodes enable a high-energy density, but suffer from an inhomogeneous reaction current and irreversible Li plating during fast charging. In contrast, hard carbon exhibits superior rate performance but lower energy density owing to its lower initial coulombic efficiency and higher average voltage. In this work, these tradeoffs are overcome by fabricating hybrid anodes with uniform mixtures of graphite and hard carbon, using industrially-relevant multi-layer pouch cells (>1 Ah) and electrode loadings (3 mAh cm –2 ). By controlling the graphite/hard carbon ratio, this study shows that battery performance can be systematically tuned to achieve both high-energy density and efficient fast charging. Pouch cells with optimized hybrid anodes retain 87% and 82% of their initial specific energy after 500 cycles of 4C and 6C fast-charge cycling, respectively. This is significantly higher than the 61% and 48% specific energy retention with graphite anodes under the same conditions. The enhanced performance is attributed to improved homogeneity of the reaction current throughout the hybrid anode, which is supported by continuum-scale modeling. Furthermore, this process is directly compatible with existing roll-to-roll battery manufacturing, representing a scalable pathway to fast charging.

25 ENERGY STORAGE↗
DOE OSTI2023

Composite Lithium Metal Structure to Mitigate Pulverization and Enable Long‐Life Batteries

In lithium metal batteries, non‐uniform stripping of lithium results in pit formation, which promotes subsequent non‐uniform, dendritic deposition. This viscous cycle leads to pulverization of lithium which promotes cell shorting or capacity degradation, symptoms further exaggerated by high electrode areal loading and lean electrolytes. Here, to address this challenge, a composite lithium metal anode is engineered that contains uniformly distributed, nanometer‐sized carbon particles. This composite lithium is shown to strip more uniformly since the growth of non‐uniform pits is intercepted by the carbon particles. This mechanism is corroborated by a continuum electrochemical model. Subsequent lithium deposition on carbon particles is also found to be more uniform than on the surface with irregular pits. Notably, the pulverization rate of composite lithium is 26 times slower than that of commercial lithium. Moreover, in a Li‐S battery with sulfurized polyacrylonitrile cathode, the use of the composite anode extends the cycle life by three times when the areal capacity is 8 mAh cm −2 . The approach of using an engineered lithium composite structure to address challenges during both stripping and plating can inform future designs of lithium metal anodes for high areal capacity operations.

high areal capacity↗
DOE OSTI2023

Optimizing the Composite Cathode Microstructure in All‐Solid‐State Batteries by Structure‐Resolved Simulations

Abstract All‐solid‐state batteries are considered as an enabler for applications requiring high energy and power density. However, they still fall short of their theoretical potential due to various limitations. One issue is poor charge transport kinetics resulting from both material inherit limitations and non‐optimized design. Therefore, a better understanding of the relevant properties of the cathode microstructure is necessary to improve cell performance. In this article, we identify optimization potentials of the composite cathode by structure‐resolved electrochemical 3D‐simulations. In our simulation study, we investigate the influence of cathode active material fraction, density, particle size, and active material properties on cell performance. Special focus is set on the impact of grain boundaries on the cathode design. Based on our simulation results, we can predict target values for cell manufacturing and reveal promising optimization strategies for an improved cathode design.

25 ENERGY STORAGE↗
DOE OSTI2025

Thermal Gradient Effects on Local Hotspot Ignition in 1,3,5,7‐Tetranitro‐1,3,5,7‐tetrazocane (HMX)

Understanding hotspot ignition, growth, and criticality, as well as the timescales of each, is crucial for parameterizing mesoscale and continuum‐level models that rely on a statistical understanding of hotspots. However, these models often consider hotspots to have uniform temperatures or for hotspots of a given temperature and size to always behave the same. Therefore, using molecular dynamics simulations, we assess the influence of thermal distribution effects on hotspot local ignition and time to ignition in 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane (HMX) nanoscale hotspots. Finally, by assessing hotspots with a gradient driven from an initial two‐temperature core–shell setup, we show that small increases in the core temperature under a constant average temperature can lead to order‐of‐magnitude effects on reaction and local ignition timescales.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Design of a Simple Flyer Plate Launcher

The response of materials to shock loading is important to understand for a variety of applications. When shock physics emerged during and after WWII, direct explosive loading or explosively driven plate impact was the primary tool for these studies. Subsequent decades have seen the widespread use of large caliber guns for plate impact studies, laser-shock facilities and pulsed power facilities. INL currently lacks a gun suitable for plate impact or explosives casting and machining facilities; however, it does possess explosives use and handling capabilities. An option for performing plate impact experiments was needed, therefore continuum scale models were utilized to explore a few simple donor-acceptor explosive plane wave lens designs, one of which could be hand packed with plastic explosives to launch flyer plates. 2D simulations were performed to study different geometries in an effort to minimize the difference in shock arrival across the central portion of a small copper flyer plate. A shock wave arrival time difference under 50 ns across 50 percent of the center of the flyer was achieved with a few designs. This work summarizes the computational models and results.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2021

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ameloblastin binding to biomimetic models of cell membranes – A continuum of intrinsic disorder

A 37-residue amino acid sequence corresponding to the segment encoded by exon-5 of murine ameloblastin (Ambn), AB2 (Y67-Q103), has been implicated with membrane association, ameloblastin self-assembly, and amelogenin-binding. Here, our aim was to characterize, at the residue level, the structural behavior of AB2 bound to chemical mimics of biological membranes using NMR spectroscopy. To better define the structure of AB2 using NMR-based methods, recombinant 13 C- and 15 N-labelled AB2 (*AB2) was prepared and data collected free in solution and with deuterated dodecylphosphocholine (dPC) micelles, deuterated bicelles, and both small and large unilamellar vesicles. Amide chemical shift and intensity perturbations observed in 1 H- 15 N HSQC spectra of *AB2 in the presence of bicelles and dPC micelles suggest that a region of *AB2, S6-E36 (murine Ambn S68 – E98), associates with the membrane biomimetics. A CSI-3 analysis of the NMR chemical shift assignments for *AB2 free in solution and bound to dPC micelles indicated the peptide remains disordered except for the adoption of a short, 12-residue α-helix, F10-G21 (murine Ambn F72-G83). In dPC micelles, the NOE NMR data was void of patterns characteristic of long-lived helical structure indicating this helix was transient in nature. A continuum of intrinsic disorder in the membrane-bound state may be responsible for ameloblastin’s ability to dynamically interact with multiple partners at the same site during amelogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Recent progress in the phase-field dislocation dynamics method

The phase-field dislocation dynamics (PFDD) method, originated in 2002, is a continuum dislocation model that uses order parameters to describe dislocation slips in crystalline materials. In the past two decades, and especially since it was last reviewed in 2016, PFDD was advanced significantly in terms of the mathematical formulation, numerical implementation, and applicability. The main purpose of this short review is to summarize recent progress made to improve the energy functional formulation and numerical techniques of PFDD as well as its recent applications. Additionally, some recommendations for future work to further extend the PFDD method are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding metal propagation in solid electrolytes due to mixed ionic-electronic conduction

Metal penetration into a solid electrolyte (SE) is one of the critical problems impeding the practical application of solid-state batteries. In this study, we investigate the conditions under which electronic conductivity of the SE can lead to metal deposition and fracture within the SE. Three different stages for void filling (metal plating initiation, metal growth, and metal compression) in the SE are identified. We show that a micron-size isolated void in the SE near the anode can be quickly filled in by metal and fractured when the developed pressure in the void grows larger than the maximum pressure the SE material can sustain. We find that the anode voltage and applied current density play a significant role in determining the vulnerability to metal deposition. We discuss several strategies to prevent electronic conductivity-driven metal propagation in electrolytes that are not fully dense, including the densified layers between the anode and SE.

25 ENERGY STORAGE↗
DOE OSTI2023

DFT-based investigation of solvation of Nd(III)/Yb(III) cations in room-temperature ionic-liquid

We report a DFT study on the solvation of bare rare-earth (RE) cations [Nd(III) and Yb(III)] in the room-temperature ionic-liquid (RTIL) based on 1,3-dimethyl-imidazolium + , PF 6 - , [MMI][PF 6 ]. In the present article, we aim to investigate whether the chosen RTIL can preferentially solvate any of the RE cations under study. First, the RTIL binding to the RE cation is studied in the gas phase, and then explored with continuum solvation model employing RTIL parameters. First, we considered the formation of the first solvation shell, where the PF 6 - anions directly bond with the Ln(III) cation. We then further surrounded the first solvation shell with [MMI]+ cations to form the second solvation shell. The stability of PF 6 - coordinated complexes is studied by calculating the binding energy values. The number of PF 6 - anions for both Ln(III) cations starts with three ligands and goes until the maximum binding energy in the RTIL medium is reached. The coordination number and binding modes of PF 6 - depend on the size of the RE cation. Comparing the binding energy values between the RE cations for the first solvation shell complexes, it is noticed that Yb-complexes have relatively greater binding energy than Nd-complexes. We have further performed a Boltzmann population analysis to find out the most dominant conformer among each category of PF 6 - coordinated complexes. A study using Natural Population Analysis shows that after complexation with PF 6 - the charge on the Yb-centers lower to a greater extent than Nd-centers, implying stronger bonding between Yb and PF 6 - , which is in harmony with the binding energy analysis. Overall, the computational investigation provides an in-depth understanding of the complexation phenomena of RE cations in RTIL, which can be useful for the practical separation of RE cations using RTIL.

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

Ab initio calculations of scattering and fusion reactions

We review recent advancements achieved within the ab initio no-core shell model with continuum (NCSMC), a unified first-principle (or, ab initio) approach to nuclear bound and continuum states. We highlight its application to reactions of astrophysical interests, including the 7 Be(p, γ) 8 B solar fusion and α-α scattering, and to the structure of the exotic 12 Be nucleus. Furthermore, our findings provide a solid foundation for integrating ab initio calculations with experimental measurements, allowing for more precise evaluations of key thermonuclear reaction rates at astrophysical energies, and offer new insights into the factors that contribute to the emergence of clustering.

Quaglioni, Sofia [Lawrence Livermore National Labo↗
DOE OSTI2025

Ab initio informed 20 Ne(p, pα) 16 O reaction elucidates the emergence of alpha clustering from chiral potentials

We report on the first ab initio informed α knock-out reaction in the intermediate-mass region, with the aim to probe the underlying chiral potential and its impact on the emergence of alpha clustering in this mass region. The theoretical predictions of the α+ 16 O clustering in the 20 Ne ground state, based on the ab initio symmetry-adapted no-core shell model with continuum, yield a triple differential cross section for 20 Ne(p, pα) 16 O that is in a remarkable agreement with the data. This allows us to examine predictions of surface and in-medium α-cluster features that emerge from the underlying realistic nucleon-nucleon interaction with no parameters fitted to nuclear data beyond the two-body system, and to compare these to the successful antisymmetrized molecular dynamics approach.

Sargsyan, G. H. [Michigan State Univ., East Lansin↗
DOE OSTI2020

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

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

G 0 W 0 Ionization Potentials of First-Row Transition Metal Aqua Ions

Here, we report computations of the vertical ionization potentials within the GW approximation of the near-complete series of first-row transition metal (V–Cu) aqua ions in their most common oxidation states, i.e., V 3+ , Cr 3+ , Cr 2+ , Mn 2+ , Fe 3+ , Fe 2+ , Co 2+ , Ni 2+ , and Cu 2+ . The d-orbital occupancy of these systems spans a broad range from d 2 to d 9 . All of the structures were first optimized at the density functional theory level using a large cluster of explicit water molecules that are embedded in a continuum solvation model. Vertical ionization potentials were computed with the one-shot G 0 W 0 approach on a range of transition metal ion clusters (6, 18, 40, and 60 explicit water molecules), wherein the convergence with respect to the basis set size was evaluated using the systems with 40 water molecules. We assess the results using three different density functional approximations as starting points for the vertical ionization potential calculations, namely, G 0 W 0 @PBE, G 0 W 0 @PBE0, and G 0 W 0 @r 2 SCAN. While the predicted ground-state structures are similar to all three exchange-correlation functionals, the vertical ionization potentials were in closer agreement with experiment when using the G 0 W 0 @PBE0 and G 0 W 0 @r 2 SCAN approaches, with the r 2 SCAN-based calculations being significantly less expensive. Computed bond distances and vertical ionization potentials for all structures are in good agreement with available experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗