Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Process intensification approach to enhancing heat and mass transfer: Radio frequency (RF) assisted drying of paper and board

Conventional multi-cylinder drying of paper and board typically relies on both conductive drying from steam-heated cylinders and convective drying, where heated air flows over the paper web surface. Conduction primarily contributes to heat transfer, while convection is the main driver of mass transfer. However, conventional drying systems are heavily dependent on steam, usually powered by fossil fuels, and are often energy-inefficient with high levels of waste. Additionally, these surface-driven processes result in a lower percentage of energy absorption compared to the energy supplied, leading to significant energy losses. To improve this long-standing process, an experimental system was developed to investigate a process intensification approach involving the integration of Radio Frequency (RF) heating, a volumetric electromagnetic technology, alongside traditional conduction and convection drying methods. This study also emphasizes the use of sensors to continuously monitor key parameters such as moisture content, supply system temperatures, sample temperatures, air flows, and drying rates in real-time. The effect of RF as an auxiliary energy source in localized environments at varying moisture levels was explored to optimize industrial drying systems, quantify potential improvements, and provide insights for future studies. Experimental results from trials combining RF with convection and with the base case alternating conduction-convection drying processes are presented. It is shown that RF is a viable process intensification approach for paper drying improving the drying rate and energy intensity at higher moisture contents. These findings offer valuable insights for process intensification and contribute to the process development, modeling, and simulation of advanced paper drying techniques.

42 ENGINEERING↗
DOE OSTI2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Deep groundwater↗
DOE OSTI2021

The detrimental effect of elemental contaminants when using B additions to improve the creep properties of a Ni-based superalloy

The effect of Si contamination when using B to improve the creep properties of a Ni-based superalloy was investigated using advanced characterization techniques and first-principles simulations on alloys with high and low B levels with varying Si contents. The positive effect of B segregation along grain boundaries on the creep properties was mitigated by the presence of Si which showed a similar segregation preference. Density functional theory calculations were used for validation by calculating grain boundary cleavage energies. Silicon was shown to decrease grain boundary cohesion which offsets the positive effect of B.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Real-Space Local Dynamics of Molten Inorganic Salts Using Van Hove Correlation Function

Molten inorganic salts are attracting resurgent attention because of their unique physicochemical properties, making them promising media for next-generation concentrating solar power systems and molten salt reactors. The dynamics of these highly disordered ionic media is largely studied by theoretical simulations, while the robust experimental techniques capable of observing local dynamics are not well-developed. To provide fundamental insights into the atomic-scale transport properties of molten salts, we report the real-space dynamics of molten magnesium chloride at high temperatures employing the Van Hove correlation function obtained by inelastic neutron scattering. Overall, our results directly depict the distance-dependent dynamics of a molten salt on the picosecond time scale. This study demonstrates the capability of the developed approach to describe the locally correlated- and self-dynamics in molten salts, significantly improving our understanding of the interplay between microscopic structural parameters and their dynamics that ultimately control physical properties of condensed matter in extreme environments.

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

Defect Annihilation Pathways in Directed Assembly of Lamellar Block Copolymer Thin Films

Defects in highly ordered self-assembled block copolymers represent an important roadblock toward the adoption of these materials in a wide range of applications. This work examines the pathways for annihilation of defects in symmetric diblock copolymers in the context of directed assembly using patterned substrates. Past theoretical and computational studies of such systems have predicted minimum free energy pathways that are characteristic of an activated process. However, they have been limited to adjacent dislocations with opposite Burgers vectors. By relying on a combination of advanced sampling techniques and particle-based simulations, this work considers the long-range interaction between dislocation pairs, both on homogeneous and nanopatterned substrates. As illustrated here, these interactions are central to understanding the defect structures that are most commonly found in applications and in experimental studies of directed self-assembly. More specifically, it is shown that, for dislocation dipoles separated by several lamellae, multiple consecutive free energy barriers lead to effective kinetic barriers that are an order of magnitude larger than those originally reported in the literature for tightly bound dislocation pairs. Furthermore, it is also shown that annihilation pathways depend strongly on both the separation between dislocations and their relative position with respect to the substrate guiding stripes used to direct the assembly.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ring-Opening Dynamics of the Cyclopropyl Radical and Cation: the Transition State Nature of the Cyclopropyl Cation

Here, we provide compelling experimental and theoretical evidence for the transition state nature of the cyclopropyl cation. Synchrotron photoionization spectroscopy employing coincidence techniques together with a novel simulation based on high-accuracy ab initio calculations reveal that the cation is unstable via its allowed disrotatory ring-opening path. The ring strains of the cation and the radical are similar, but both ring opening paths for the radical are forbidden when the full electronic symmetries are considered. These findings are discussed in light of the early predictions by Longuet-Higgins alongside Woodward and Hoffman; we also propose a simple phase space explanation for the appearance of the cyclopropyl photoionization spectrum. The results of this work allow the refinement of the cyclopropane C–H bond dissociation energy, in addition to the cyclopropyl radical and cation cyclization energies, via the Active Thermochemical Tables approach.

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

Purification of Lithium-Ion Battery Black Mass through Tailored Alkaline Corrosion

Obtaining high-purity material outputs is crucial to the viability of novel process aimed at direct recycling of lithium-ion batteries. Metallic impurities in recycled cathodes have been shown to inhibit performance, thereby threatening mainstream acceptance of recycled battery products. Thus, shredded black mass (BM) must be purified to remove metallic contaminants, and specifically Al and Cu originating from the electrode current collectors. We herein explore a process to ionize solid copper and aluminum to ionic form based on tailored alkaline chemistry, without incurring damage to the target cathode material (Li(NixMnyCo1-x-y)O2; NMC). Al and Cu corrosion may be enhanced through the addition of chloride salt, elevated temperatures, and the use of ultrasonication - all of which disrupt the formation of passivating films on the metallic surface, and thereby increase corrosion rate. We demonstrate optimized parameters for Al and Cu corrosion both from a kinetic and overall process cost perspective. Further, we analyze the impact of these conditions on the structural (XRD, SEM), chemical (EDS, ICP), and electrochemical (impedance, cycling, dQ/dV) properties of NMC, and suggest that the present purification method does not significantly disrupt NMC performance. Finally, we present preliminary results from a promising bench-scale demonstration of this purification process applied to a simulated black mass.

ADVANCED PROPULSION SYSTEMS,ENERGY STORAGE↗
DOE OSTI2024

The Role of Cold Pools in Modulating Convective Organization During the MJO

Abstract In this study, we investigate the role of cold pools in modulating convective organization throughout the Madden‐Julian Oscillation (MJO) life cycle using a modeling approach that combines Eulerian and Lagrangian techniques. First, we conduct a simulation using the soundings and forcing of the DYNAMO/AMIE campaign. The simulation shows a lag of several days between the precipitation rate peak time associated with MJO and the highest convective organization time. Second, to analyze the role of cold pools, we consider a series of 2‐day simulations conducted at different stages of the MJO. The simulation results suggest that cold pools are larger and last longer during the mature stages of the MJO, possibly because of decreased environmental surface latent heat fluxes and stronger downdrafts. These lead to the formation of moist rings at the leading edges of cold pools, facilitating the formation of more convective cores and increasing the degree of convective organization.

Geology↗
DOE OSTI2020

Dimer rattling mode induced low thermal conductivity in an excellent acoustic conductor

A solid with larger sound speeds usually exhibits higher lattice thermal conductivity. Here, we report an exception that CuP2 has a quite large mean sound speed of 4155 m s –1 , comparable to GaAs, but single crystals show very low lattice thermal conductivity of about 4 W m –1 K –1 at room temperature, one order of magnitude smaller than GaAs. To understand such a puzzling thermal transport behavior, we have thoroughly investigated the atomic structures and lattice dynamics by combining neutron scattering techniques with first-principles simulations. This compound crystallizes in a layered structure where Cu atoms forming dimers are sandwiched in between P atomic networks. In this work, we reveal that Cu atomic dimers vibrate as a rattling mode with frequency around 11 meV, which is manifested to be remarkably anharmonic and strongly scatters acoustic phonons to achieve the low lattice thermal conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evolution of dislocations during the rapid solidification in additive manufacturing

Materials processed by fusion-based additive manufacturing (AM) typically exhibit relatively high dislocation densities, along with cellular structures and elemental segregation. This representative structural feature significantly influences material performance; however, post-mortem microstructure characterizations of AM materials cannot capture the dynamic evolution of dislocations during the manufacturing process, thereby offering limited mechanism-based guidance for further advancing AM techniques and facilitating the qualification and certification of AM products. In this study, we conduct operando high-energy synchrotron X-ray diffraction experiments on wire-laser directed energy deposition of 316 L stainless steel. Through a unique configuration, our operando synchrotron experiments semi-quantitatively probe the dislocation density in solid phases and their dynamic changes during solidification and subsequent cooling. By integrating this advanced synchrotron technique with multi-physics simulation, in-situ neutron diffraction, and multi-scale electron microscopy characterization, our mechanistic study aims to elucidate the effects of rapid cooling and subsequent thermal cycling on the dislocation generation and evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Photoionization of the iodine 3d, 4s, and 4p orbitals in methyl iodide

Ionization of the I 3d, 4s, and 4p orbitals in methyl iodide (CH 3 I) has been studied by using synchrotron radiation to measure the total ion yield and by recording photoelectron spectra with linearly polarized radiation in two polarization orientations. The complete photoelectron spectrum of CH 3 I has been recorded at several photon energies, and bands due to the C 1s, I 3d, 4s, 4p, and 4d atomic-like orbitals, as well as the molecular orbitals, have been observed and assigned. In the vicinity of the I 3d 5/2 and 3d 3/2 ionization thresholds at 626.8 and 638.3 eV, respectively, the ion yield displays weak structure in the pre-edge region due to transitions into valence or Rydberg states, and, at higher energies, a shoulder and a broad maximum attributed to the I 3d 5/2 → εf and the I 3d 3/2 → εf shape resonances, respectively. The absorption spectrum calculated using time-dependent density functional theory, within the Tamm-Dancoff approximation, has allowed assignments to be proposed for the valence and Rydberg states. The Stieltjes imaging technique has been used to simulate the absorption spectrum above the ionization threshold and indicates that transitions into the f(l = 3) continuum channel dominate. This conclusion has been corroborated by a Continuum Multiple Scattering–Xα (CMS–Xα) calculation. The asymmetric broadening of the photoelectron bands associated with the I 3d orbital, due to post collision interaction, is taken into account in our experimental analysis. Experimentally derived photoelectron anisotropy parameters for the I 3d orbital are in good agreement with the theoretical predictions obtained with the CMS–Xα approach. The I 3d shake-up/shake-off photoelectron spectrum has been recorded, and assignments have been proposed for several of the satellites. The M 4 N 45 N 45 and M 5 N 45 N 45 Auger electron yields have been measured, and that for the M 5 N 45 N 45 decay exhibits a maximum due to interchannel coupling between the 3d5/2 and 3d3/2 continua. The photoelectron band associated with the I 4p orbital has an unusual appearance. Based upon previous theoretical work for the analogous Xe 4p orbital, it appears that the initial I 4p -1 hole state decays rapidly through Coster-Kronig and super-Coster-Kronig transitions. This leads to a redistribution of the spectral intensity associated with the I 4p orbital and results in a photoelectron spectrum containing a single structured band together with an extended continuum. Another continuum is observed on the high binding energy side of the peak due to the 4s orbital, and we assign this to super-Coster-Kronig transitions into the 4p -1 4d -1 continuum.

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

Assessing models of force-dependent unbinding rates via infrequent metadynamics

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent “bonds” being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. Furthermore, this work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future.

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

Universal and interpretable classification of atomistic structural transitions via unsupervised graph learning

Materials processing often occurs under extreme dynamic conditions leading to a multitude of unique structural environments. These structural environments generally occur at high temperatures and/or high pressures, often under non-equilibrium conditions, which results in drastic changes in the material's structure over time. Computational techniques, such as molecular dynamics simulations, can probe the atomic regime under these extreme conditions. However, characterizing the resulting diverse atomistic structures as a material undergoes extreme changes in its structure has proved challenging due to the inherently non-linear relationship between structures as large-scale changes occur. Here, we introduce SODAS++, a universal graph neural network framework, that can accurately and intuitively quantify the atomistic structural evolution corresponding to the transition between any two arbitrary phases. We showcase SODAS++ for both solid–solid and solid–liquid transitions for systems of increasing geometric and chemical complexity, such as colloidal systems, elemental Al, rutile and amorphous TiO 2 , and the non-stoichiometric ternary alloy Ag 26 Au 5 Cu 19 . Finally, we show that SODAS++ can accurately quantify all transitions in a physically interpretable manner, showcasing the power of unsupervised graph neural network encodings for capturing the complex and non-linear pathway, a material's structure takes as it evolves.

36 MATERIALS SCIENCE↗
DOE OSTI2021

PT symmetry, pattern formation, and finite-density QCD

A longstanding issue in the study of quantum chromodynamics (QCD) is its behavior at nonzero baryon density, which has implications for many areas of physics. The path integral has a complex integrand when the quark chemical potential is nonzero and therefore has a sign problem, but it also has a generalized ${\mathscr{PT}}$ symmetry. We review some new approaches to ${\mathscr{PT}}$-symmetric field theories, including both analytical techniques and methods for lattice simulation. We show that ${\mathscr{PT}}$-symmetric field theories with more than one field generally have a much richer phase structure than their Hermitian counterparts, including stable phases with patterning behavior. The case of a ${\mathscr{PT}}$-symmetric extension of a φ 4 model is explained in detail. The relevance of these results to finite density QCD is explained, and we show that a simple model of finite density QCD exhibits a patterned phase in its critical region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Unveiling the internal structure of the Hercules supercluster

ABSTRACT We have investigated the structure of the Hercules supercluster (SCL160) based on data originally extracted from the Sloan Digital Sky Survey (SDSS-DR7). We have traced the mass distribution in the field through the numerical density, weighted by the r′ luminosity of the galaxies, and classified objects based on their spatial position and redshift. This has allowed us to address not only the kinematics of the supercluster as a whole, but also the internal kinematics of each cluster, which had not been explored further before. We have confirmed that the Hercules supercluster is composed of the galaxy clusters A2147, A2151, and A2152. A2151 consists of five subclusters, A2147 two, and A2152 at least two. They form the heart of the Hercules supercluster. We also have found two other gravitationally bound clusters, therefore increasing the known members of the supercluster. We have estimated a total mass of 2.1 ± 0.2 × 1015 M⊙ for the Hercules supercluster. To determine the dynamical masses in this work, we have resorted to the M200−σ scaling relation and the caustic technique. Comparing both methods with simulated data of bimodal merging clusters, we found that the caustic, as well as the σ-based masses, is biased through the merger age, showing a boost just after pericentric passage. This is not in line with the principle of the caustic method, which affirms it does not depend on the cluster dynamical state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Lensing reconstruction in post-Born cosmic microwave background weak lensing

The study of the cosmic microwave background (CMB) lensing potential has established itself by now as a robust way of probing the physics of large-scale structure growth. The most common estimators of the lensing potential are derived under the assumption of Gaussianity of the matter distribution and in the Born approximation of the photon diffusion. In this paper we study the performance of quadratic estimators when applied to realistic sky maps extracted from multiple-lens ray tracing techniques in cosmological $N$-body simulations. These are expected to model accurately the effects due to the non-Gaussianity of the matter distribution induced by its nonlinearity and the deviation from the Born approximation. We show that both these effects on their own lead to reconstruction biases, but these tend to partially cancel each other when both these effects are considered together. We forecast the impact of these biases on the estimation of cosmological parameters for future high-sensitivity CMB experiments like CMB-S4. We find that the cold dark matter density, $Ω_{cdm}$, the optical depth to reionization $τ$, the amplitude of primordial inflationary perturbations, $A_s$, and the sum of neutrino masses, $M_ν$, could be biased at the $1–2σ$ level, if no external data set is used. We also observe a reduction of the bias if external data like baryon acoustic oscillations is included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Semi-inclusive deep-inelastic scattering, parton distributions, and fragmentation functions at a future electron-ion collider

We present a quantitative assessment of the impact a future electron-ion collider would have in the determination of parton distribution functions in the proton and parton-to-hadron fragmentation functions through semi-inclusive deep-inelastic electron-proton scattering data. Specifically, we estimate the kinematic regions for which the forthcoming data are expected to have the most significant impact in the precision of these distributions, computing the respective correlation and sensitivity coefficients. Using a reweighting technique for the sets of simulated data with their realistic uncertainties for two different center-of-mass energies, we analyze the resulting new sets of parton distribution functions and fragmentation functions, which have significantly reduced uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗