Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATOMIC MODELS”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The design of a photoionization front experiment using the Z-Machine as a driving source and estimated measurements

Radiation-driven heat fronts are present in the early universe during reionization, the circumstellar medium of supernovae, and in high-energy-density physics experiments. Dedicated experiments to observe and diagnose the behavior of these types of heat fronts can improve our understanding of these phenomena. A simulation study of photoionization fronts using the HELIOS-CR radiation hydrodynamics code provides an experimental design for the Z-Machine at Sandia National Laboratory using a measurement-calibrated input radiation flux to drive the photoionization front. The simulations use detailed atomic physics and non-diffusive radiation transport in 1D to determine an optimal gas pressure of 0.75 atm for an experiment in N gas as well as the effects of increasing the thickness of the window that seals the gas cell. Post-processing of these simulations demonstrates that ratios of atomic rate coefficients place the heat front in a physics regime where photoionization dominates the energy deposition. Here, to see the sensitivity of the simulations to changes in the model and spatial grid, this analysis performed resolution, atomic model detail, and radiation transport angular grid studies showing less than 10% deviation from the nominal model for increased complexity, when possible. An effort to emulate 3D geometric effects on the radiation flux using an artificial attenuation scheme has shown that, even for conservative estimates of the flux, simulations still produce a photoionization front. Estimations of a streaked, visible spectroscopy measurement using SPECT3D showed that line emission measurements are present early in time and that later in time thermal emission should become dominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The multiscale coarse-graining method. XI. Accurate interactions based on the centers of charge of coarse-grained sites

It is essential to be able to systematically construct coarse-grained (CG) models that can efficiently and accurately reproduce key properties of higher-resolution models such as all-atom. To fulfill this goal, a mapping operator is needed to transform the higher-resolution configuration to a CG configuration. Certain mapping operators, however, may lose information related to the underlying electrostatic properties. In this paper, a new mapping operator based on the centers of charge of CG sites is proposed to address this issue. Four example systems are chosen to demonstrate this concept. Within the multiscale coarse-graining framework, CG models that use this mapping operator are found to better reproduce the structural correlations of atomistic models. The present work also demonstrates the flexibility of the mapping operator and the robustness of the force matching method. For instance, important functional groups can be isolated and emphasized in the CG model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Nucleon Isovector Pairing in Nuclei: Microscopic Approach, Boson Representation, and Collective Model

Nucleon pair correlations in atomic nuclei are analyzed within a nuclear microscopic model with residual isovector pairing forces. These are formulated in the boson representation of fermion operators whereby the collective mode of pair excitations can be isolated without restricting the size of the one-particle basis. This method allows one to analyze the fluctuations in the nonsuperfluid phase of nuclear matter, its phase transition to the superfluid phase, and strong pair correlations. The performance of the method is exemplified by numerical results for the nuclei in the vicinity of the doubly magic {sup 56}Ni nucleus.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Modeling the influence of nozzle-generated turbulence on diesel sprays

The physical mechanisms governing spray breakup in direct injection engines, such as aerodynamic induced instabilities and nozzle-generated cavitation and turbulence, are not well understood due to the experimental and computational limitations in resolving these processes. Recent x-ray and visible extinction measurements have been con-ducted with a targeted interest in the spray formation region in order to characterize the distribution of droplet sizes throughout the spray. Detailed analysis of these measurements shows promise of yielding insight into likely mechanisms governing atomization, which can inform the improvement of spray models for engine computational fluid dynamic (CFD) codes. In order to investigate potential atomization mechanisms, we employ a joint experimental and computational approach to characterize the structure of the spray formation region using the Engine Combustion Network Spray D injector. X-ray tomography, radiography and ultra-small angle x-ray scattering measurements conducted at the Advanced Photon Source at Argonne National Laboratory quantify the injector geometry, liquid fuel mass and Sauter mean diameter (SMD) distributions under non-vaporizing conditions. Diffused back-illumination imaging measurements, conducted at the Georgia Institute of Technology, characterize the asymmetry of the spray structure. The selected range of injection pressures (50 – 150 MPa) and ambient densities (1.2 – 22.8 kg/m3) allow for the influence of aerodynamic forces on the spray to be studied in a controlled and systematic manner, while isolating the atomization process from the effects of vaporization. In comparison to high ambient density conditions, the spray is observed to be more asymmetric at low ambient density conditions. Although several mechanisms may cause asymmetries in the nozzle exit flow conditions and ultimately the spray distribution, irregularities in the internal nozzle geometry were identified, suggesting an increased sensitivity of the spray structure to internal nozzle surface finish imperfections at such conditions. The presence of these asymmetries may influence the ability to interpret line-of-sight measurements and their derived SMD values and trends from a single viewing angle of the spray. With this consideration in mind, the measured local sensitivities to ambient density suggest that for ambient densities less than 2.4 kg/m3, aerodynamic effects are likely suppressed, allowing the influence of turbulent-induced breakup to be isolated. In concert with the experimental measurements, we utilize three-dimensional, CFD Lagrangian-Eulerian spray simulations in CONVERGE to evaluate the details of the predicted spray structure. In particular, we compare measured and predicted sensitivities of the SMD distribution to changes in injection and ambient conditions from three different atomization models, namely Kelvin Helmholtz (KH), KH Aerodynamics Cavitation Turbulence (KH-ACT), and the newly developed KH-Faeth hybrid model. While none of the existing hybrid spray models were able to replicate the experimentally observed sensitivities, it was found that the scales characterizing the KH-Faeth model show promise of capturing the experimentally observed trends if the effects of secondary droplet breakup are neglected. These results inform recommendations for future experiments and computational studies that can guide the development of an improved spray breakup model.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2014

Density functional theory and evolution algorithm calculations of elastic properties of AlON

Different models for aluminum oxynitride (AlON) were calculated using density functional theory and optimized using an evolutionary algorithm. Evolutionary algorithm and density functional theory (DFT) calculations starting from several models of AlON with different Al or O vacancy locations and different positions for the N atoms relative to the vacancy were carried out. The results show that the constant anion model [McCauley et al., J. Eur. Ceram. Soc. 29(2), 223 (2009)] with a random distribution of N atoms not adjacent to the Al vacancy has the lowest energy configuration. The lowest energy structure is in a reasonable agreement with experimental X-ray diffraction spectra. The optimized structure of a 55 atom unit cell was used to construct 220 and 440 atom models for simulation cells using DFT with a Gaussian basis set. Cubic elastic constant predictions were found to approach the experimentally determined AlON single crystal elastic constants as the model size increased from 55 to 440 atoms. The pressure dependence of the elastic constants found from simulated stress-strain relations were in overall agreement with experimental measurements of polycrystalline and single crystal AlON. Calculated IR intensity and Raman spectra are compared with available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mechanistic insights into carbon–carbon coupling on NiAu and PdAu single-atom alloys

Carbon–carbon coupling is an important step in many catalytic reactions, and performing sp 3 –sp 3 carbon–carbon coupling heterogeneously is particularly challenging. It has been reported that PdAu single-atom alloy (SAA) model catalytic surfaces are able to selectively couple methyl groups, producing ethane from methyl iodide. In this work, we extend this study to NiAu SAAs and find that Ni atoms in Au are active for C–I cleavage and selective sp 3 –sp 3 carbon–carbon coupling to produce ethane. Furthermore, we perform ab initio kinetic Monte Carlo simulations that include the effect of the iodine atom, which was previously considered a bystander species. We find that model NiAu surfaces exhibit a similar chemistry to PdAu, but the reason for the similarity is due to the role the iodine atoms play in terms of blocking the Ni atom active sites. Specifically, on NiAu SAAs, the iodine atoms outcompete the methyl groups for occupancy of the Ni sites leaving the Me groups on Au, while on PdAu SAAs, the binding strengths of methyl groups and iodine atoms at the Pd atom active site are more similar. These simulations shed light on the mechanism of this important sp 3 –sp 3 carbon–carbon coupling chemistry on SAAs. Furthermore, we discuss the effect of the iodine atoms on the reaction energetics and make an analogy between the effect of iodine as an active site blocker on this model heterogeneous catalyst and homogeneous catalysts in which ligands must detach in order for the active site to be accessed by the reactants.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Advancing Turbulent Spray and Combustion Models for Compression Ignition Engine Simulations

This thesis seeks to investigate the turbulent mixing influence on spray atomization and combustion processes encountered in compression ignition diesel engines. Despite greater thermal efficiency of diesel engine than spark ignition engine, the nature of stratified air-fuel mixture and non-premixed flame gives rise to unacceptable levels of nitrogen oxides (NOx) and particulate matter (PM), thus the use of diesel engines has often been limited to heavy-duty vehicle and industrial power sources. However, recent advancement in diesel engine combustion strategies, e.g. low temperature com-bustion (LTC), has demonstrated promising pathways towards improvement in the engine-out pollutants. Therefore, particularly in the effort of computer-aided en-gine design tasks, such a new engine design concept requires more accurate modeling techniques applicable over a broader range of engine operating conditions than those of conventional engine strategies. In the notion of challenges in new engine operat-ing conditions, this thesis aims to present successful implementation of improvement in numerical modeling techniques in high-pressure spray atomization and resulting turbulent spray flame of interest. Three-dimensional Computational Fluid Dynamics (CFD) in in-cylinder turbu-lent combustion is considered an integral part of engine design progress, but rather a cost-prohibitive to apply over a broad range of engine relevant conditions. In spite of successful use of existing spray atomization modeling, prior researchers have pointed out some degree of failure in LTC targeted injection strategies. Furthermore, finite rate and strong nonlinearity of chemistry influenced by local turbulent mixing still re-main in challenges to account for in cost-efficient CFD analysis. In this context, a new attempt of hybrid spray primary breakup modeling is presented and demonstrated in successful application aimed at LTC technique. In addition, the Representative Interactive Flamelets (RIF) model with aid of multi-flamelets approach is extensively assessed in terms of predictive capability against classical combustion model. The combustion model employed in this study are fully examined in the general diesel combustion metric, e.g., ignition delay and flame lift-off length as well as newly sug-gested test metric, combustion recession. The combustion recession has been recently idenfied, but still remain largely unknown. Since the governing physics of this phe-nomenon is characterized by turbulent mixing coupled with finite rate chemistry, this can be considered as a relevant test metric for turbulent combustion models. In ad-dition, very recent experimental studies have introduced a new non-sooting diesel combustion technique by manipulating direct injection method. The ducted fuel in-jection (DFI) has thus been demonstrated with its potential of low soot emissions. Knowing that the duct equipped ahead of injector nozzle was identifed to enhance turbulent mixing, investigations of DFI combustion may prove the effectiveness of turbulece-chemistry interaction modeling. This thesis presents comprehensive under-standings of aformentioned diesel combustion techniques in terms of several important physics keywords, e.g., turbulent mixing and detailed chemistry

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Enhanced analysis of experimental x-ray spectra through deep learning

X-ray spectroscopic data from high-energy-density laser-produced plasmas has long required thorough, time-consuming analysis to extract meaningful source conditions. There are often confounding factors due to rapidly evolving states and finite spatial gradients (e.g., the existence of multi-temperature, multi-density, multi-ionization states, etc.) that make spectral measurements and analysis difficult. Here, in this paper, we demonstrate how deep learning can be applied to enhance x-ray spectral data analysis in both speed and intricacy. Neural networks (NNs) are trained on ensemble atomic physics simulations so that they can subsequently construct a model capable of extracting plasma parameters directly from experimental spectra. Through deep learning, the models can extract temperature distributions as opposed to single or dual temperature/density fits from standard trial-and-error atomic modeling at a significantly reduced computational cost compared to traditional trial-and-error methods. These NNs are envisioned to be deployed with high repetition rate x-ray spectrometers in order to provide detailed real-time analysis of experimental spectra.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Integrated Data Analysis to Improve Te and Zeff Measurement on NSTX-U

The multi-energy soft x-ray (ME-SXR) camera is a new concept for a fusion plasma diagnostic based on a pixelated detector with adjustable photon cutoff thresholds. It is a single, compact, solid-state detector which combines spatial and spectral sensitivity and can be configured remotely at run time, making it ideal for a fusion environment. The detector can be configured for application to a wide range of physics problems including impurity transport, electron temperature profiles, and non-thermal effects. The detector has been installed on the Madison Symmetric Torus (MST) as part of a collaboration between Princeton Plasma Physics Laboratory and the University of Wisconsin–Madison to develop the technology and analysis techniques. A full diagnostic forward model of the ME-SXR diagnostic has been developed to enable analysis of the raw photon count data. This model incorporates the atomic physics underlying the plasma emission, diagnostic geometry, and detector response characteristics. Employing this model, a Bayesian analysis approach is used to extract electron temperature from the raw data. This analysis approach uses all available information, clearly quantifies prior information, and handles uncertainties transparently. The ME-SXR camera is being used to study the 3D (helical, stellarator-like) state that forms in the MST device. This helical core structure has improved energy confinement. The transition to the helical state is accompanied by a rise in electron temperature in the central part of the helical structure, but then this enhancement fluctuates. The characteristics of this fluctuation are being compared to a magneto-hydrodynamics model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

A Method for Combining Experimentation and Molecular Dynamics Simulation to Improve Cohesive Zone Models for Metallic Microstructures

Fracture processes within a material begin at the nanometer length scale at which the formation, propagation, and interaction of fundamental damage mechanisms occur. Physics-based modeling of these atomic processes quickly becomes computationally intractable as the system size increases. Thus, a multiscale modeling method, based on the aggregation of fundamental damage processes occurring at the nanoscale within a cohesive zone model, is under development and will enable computationally feasible and physically meaningful microscale fracture simulation in polycrystalline metals. This method employs atomistic simulation to provide an optimization loop with an initial prediction of a cohesive zone model (CZM). This initial CZM is then applied at the crack front region within a finite element model. The optimization procedure iterates upon the CZM until the finite element model acceptably reproduces the near-crack-front displacement fields obtained from experimental observation. With this approach, a comparison can be made between the original CZM predicted by atomistic simulation and the converged CZM that is based on experimental observation. Comparison of the two CZMs gives insight into how atomistic simulation scales.

Hochhalter, J. D.↗
DOE OSTI2026

Systems-To-Atoms (S2A): enabling hydrogen for climate security

The project addresses a critical gap in hydrogen infrastructure by integrating system-level energy models with atomic-scale material simulations in a unified Systems-to-Atoms (S2A) framework. The motivation stems from the need to develop efficient, cost-effective, and durable hydrogen transport and utilization technologies to support decarbonization of hard-to-electrify sectors such as heavy-duty transportation. Current system models lack awareness of material performance mechanisms, while material-scale models do not account for system-level usage and variability. To bridge this divide, the team developed a co-simulation capability linking techno-economic analyses, reactor/process-flow modeling, and molecular-scale catalysis simulations. Applied to hydrogen delivery in California, the framework enabled comparative evaluations of compressed, cryogenic, and liquid organic hydrogen carrier (LOHC) pathways, highlighting how catalyst operation and unit process efficiency influence overall performance. The results demonstrate that no single material or transport mode is universally optimal; instead, heterogeneous solutions tuned to specific operational contexts deliver better performance. The project delivers a new capability for cross-scale material co-design, advancing hydrogen infrastructure readiness and informing DOE and LLNL missions in climate and energy resilience.

organic↗
NASA NTRS2004

Spatial Distributions of Metal Atoms During Carbon SWNTs Formation: Measurements and Modelling

Experiments and modelling have been undertaken to clarify the role of metal catalysts during single-wall carbon nanotube formation. For instance, we wonder whether the metal catalyst is active as an atom, a cluster, a liquid or solid nanoparticle [1]. A reactor has been developed for synthesis by continuous CO2-laser vaporisation of a carbon-nickel-cobalt target in laminar helium flow. The laser induced fluorescence technique [2] is applied for local probing of gaseous Ni, Co and CZ species throughout the hot carbon flow of the target heated up to 3500 K. A rapid depletion of C2 in contrast to the spatial extent of metal atoms is observed in the plume (Fig. 1). This asserts that C2 condenses earlier than Ni and Co atoms.[3, 4]. The depletion is even faster when catalysts are present. It may indicate that an interaction between metal atoms and carbon dimers takes place in the gas as soon as they are expelled from the target surface. Two methods of modelling are used: a spatially I-D calculation developed originally for the arc process [5], and a zero-D time dependent calculation, solving the chemical kinetics along the streamlines [6]. The latter includes Ni cluster formation. The peak of C2 density is calculated close to the target surface where the temperature is the highest. In the hot region, C; is dominant. As the carbon products move away from the target and mix with the ambient helium, they recombine into larger clusters, as demonstrated by the peak of C5 density around 1 mm. The profile of Ni-atom density compares fairly well with the measured one (Fig. 2). The early increase is due to the drop of temperature, and the final decrease beyond 6 mm results from Ni cluster formation at the eutectic temperature (approx.1600 K).

Cau, M.↗
DOE OSTI2022

Fluid and hybrid simulations of the ionization instabilities in Hall thruster

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. Further, the predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Atomistic Study of the Plastic Deformation of Transition Metals and High Entropy Alloys (Final Report)

The overarching goal of the research performed in the studies supported by the DOE Basic Energy Sciences has been investigation of mechanical properties of materials employed in technologies related to energy production and conversion using atomic level computer modeling. The materials studied are body-centered-cubic (BCC) transition metals, including α-iron, their alloys and most recently also compositionally complex, so called high entropy alloys (HEAs), which are believed to become future materials in high temperature applications. The goal of the research has been to achieve atomic level understanding of the deformation and fracture of these materials starting from the electronic structure, coarse grain the material energetics into potentials describing interaction between atoms and proceed via atomistic modeling and mesoscopic models to investigation of the glide of dislocations and thus advance the theory of plastic deformation of studied materials.

36 MATERIALS SCIENCE↗
NASA NTRS2008

Multiscale Modeling of Damage Processes in fcc Aluminum: From Atoms to Grains

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, current analysis is limited to small domains and increasing the size of the MD domain quickly presents intractable computational demands. A preferred approach to surmount this computational limitation has been to combine continuum mechanics-based modeling procedures, such as the finite element method (FEM), with MD analyses thereby reducing the region of atomic scale refinement. Such multiscale modeling strategies can be divided into two broad classifications: concurrent multiscale methods that directly incorporate an atomistic domain within a continuum domain and sequential multiscale methods that extract an averaged response from the atomistic simulation for later use as a constitutive model in a continuum analysis.

Glaessgen, E. H.↗
DOE OSTI2014

On the temperature dependence of H-U{sub iso} in the riding hydrogen model

The temperature dependence of hydrogen U{sub iso} and parent U{sub eq} in the riding hydrogen model is investigated by neutron diffraction, aspherical-atom refinements and QM/MM and MO/MO cluster calculations. Fixed values of 1.2 or 1.5 appear to be underestimated, especially at temperatures below 100 K. The temperature dependence of H-U{sub iso} in N-acetyl-l-4-hydroxyproline monohydrate is investigated. Imposing a constant temperature-independent multiplier of 1.2 or 1.5 for the riding hydrogen model is found to be inaccurate, and severely underestimates H-U{sub iso} below 100 K. Neutron diffraction data at temperatures of 9, 150, 200 and 250 K provide benchmark results for this study. X-ray diffraction data to high resolution, collected at temperatures of 9, 30, 50, 75, 100, 150, 200 and 250 K (synchrotron and home source), reproduce neutron results only when evaluated by aspherical-atom refinement models, since these take into account bonding and lone-pair electron density; both invariom and Hirshfeld-atom refinement models enable a more precise determination of the magnitude of H-atom displacements than independent-atom model refinements. Experimental efforts are complemented by computing displacement parameters following the TLS+ONIOM approach. A satisfactory agreement between all approaches is found.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Exact solution of two simple non-equilibrium electron-phonon and electron-electron coupled systems: The atomic limit of the Holstein-Hubbard model and the generalized Hatsugai-Komoto model

One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales—for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections, but strongly correlated electrons require a more complex description, and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. Then, we examine the photoemission spectra of this model in and out of equilibrium. Next, we use similar methods to exactly solve an extended version of the Hatsugai-Komoto model and examine its behavior in and out of equilibrium. These calculations lead us to propose using the first moment of the photoemission spectra to signal non-equilibrium changes in electron-electron and electron-phonon couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗