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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Orientational anisotropy in simulated vapor-deposited molecular glasses

Enhanced kinetic stability of vapor-deposited glasses has been established for a variety of glass organic formers. Several recent reports indicate that vapor-deposited glasses can be orientationally anisotropic. In this work, we present results of extensive molecular simulations that mimic a number of features of the experimental vapor deposition process. The simulations are performed on a generic coarse-grained model and an all-atom representation of N,N′-bis(3-methylphenyl)-N,N′-diphenylbenzidine (TPD), a small organic molecule whose vapor-deposited glasses exhibit considerable orientational anisotropy. The coarse-grained model adopted here is found to reproduce several key aspects reported in experiments. In particular, the molecular orientation of vapor-deposited glasses is observed to depend on substrate temperature during deposition. For a fixed deposition rate, the molecular orientation in the glasses changes from isotropic, at the glass transition temperature, T{sub g}, to slightly normal to the substrate at temperatures just below T{sub g}. Well below T{sub g}, molecular orientation becomes predominantly parallel to the substrate. The all-atom model is used to confirm some of the equilibrium structural features of TPD interfaces that arise above the glass transition temperature. We discuss a mechanism based on distinct orientations observed at equilibrium near the surface of the film, which get trapped within the film during the non-equilibrium process of vapor deposition.

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

Numerical modelling of emission of a two-level atom near a metal nanoparticle with account for tunnelling of an electron from an atom into a particle

We report the results of numerical modelling of emission of a two-level atom near a metal nanoparticle under resonant interaction of light with plasmon modes of the particle. Calculations have been performed for different polarisations of light by a dipole approximation method and a complex multipole method. Depending on the distance between a particle and an atom, the contribution of the nonradiative process of electron tunnelling from a two-level atom into a particle, which is calculated using the quasi-classical approximation, has been taken into account and assessed. We have studied spherical gold and silver particles of different diameters (10 – 100 nm). The rates of electron tunnelling and of spontaneous decay of the excited atomic state are found. The results can be used to develop nanoscale plasmonic emitters, lasers and photodetectors. (nanooptics)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Semi-grand canonical Monte Carlo simulation of the acrolein induced surface segregation and aggregation of AgPd with machine learning surrogate models

The single atom alloy of AgPd has been found to be a promising catalyst for the selective hydrogenation of acrolein. It is also known that the formation of Pd islands on the surface will greatly reduce the selectivity of the reaction. As a result, the surface segregation and aggregation of Pd on the AgPd surface under reaction conditions of selective hydrogenation of acrolein are of great interest. In this work, we lay out a workflow that can predict the surface segregation and aggregation of Pd on a FCC(111) AgPd surface with and without the presence of acrolein. We use machine learning surrogate models to predict the AgPd bulk energy, AgPd slab energy, and acrolein adsorption energy on AgPd slabs. Then, we use the semi-grand canonical Monte Carlo simulation to predict the surface segregation and aggregation under different bulk Pd concentrations. Under vacuum conditions, our method predicts that only trace amount of Pd will exist on the surface at Pd bulk concentrations less than 20%. However, with the presence of acrolein, Pd will start to aggregate as dimers on the surface at Pd bulk concentrations as low as 6.5%.

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

Atomic density functional and diagram of structures in the phase field crystal model

The phase field crystal model provides a continual description of the atomic density over the diffusion time of reactions. We consider a homogeneous structure (liquid) and a perfect periodic crystal, which are constructed from the one-mode approximation of the phase field crystal model. A diagram of 2D structures is constructed from the analytic solutions of the model using atomic density functionals. The diagram predicts equilibrium atomic configurations for transitions from the metastable state and includes the domains of existence of homogeneous, triangular, and striped structures corresponding to a liquid, a body-centered cubic crystal, and a longitudinal cross section of cylindrical tubes. The method developed here is employed for constructing the diagram for the homogeneous liquid phase and the body-centered iron lattice. The expression for the free energy is derived analytically from density functional theory. The specific features of approximating the phase field crystal model are compared with the approximations and conclusions of the weak crystallization and 2D melting theories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Empirical prediction of saline water atomization pressure loss and spray phase change using local flow pressure analysis

A pressure analysis technique was developed to obtain spray spatial evaporation profile in a thermal desalination process. The technique can replace temperature-based evaporation measurement methods that are challenged by liquid-phase interference. It also provides an alternative to optical methods with high cost and complexity. Fundamental analysis was developed to determine local humidity ratios from changes in local static and dynamic pressures. The analysis was applied to data obtained from an external mixing air-assist atomizer. An empirical model was developed for the pressure loss in atomization process and was used in the analysis of humidity ratio. The results were used to develop empirical correlations for local humidity ratio in spray mediums using inlet conditions. Utilizing inlet conditions to predict evaporation profiles is an important contribution as it eliminates the need to take local spray measurements. Obtaining local data is often difficult and expensive; Our method and correlation circumvent the need for that. Finally, the correlations presented in this paper apply to saline water sprays with total dissolved salt of 0-10 wt%. The salinity-specific correlations are accurate to within ±10% of experimental data. The individual models were aggregated into a unified correlation for salinities within 0-10 wt%. The unified correlation is accurate to within ±18% of experimental data.

42 ENGINEERING↗
DOE OSTI2017

Exploration of Turbulent Atomization Mechanisms for Diesel Spray Simulations

The atomization and initial spray formation processes in direct injection engines are not well understood due to the experimental and computational challenges associated with resolving these processes. Although different physical mechanisms, such as aerodynamic-induced instabilities and nozzle-generated turbulence and cavitation, have been proposed in the literature to describe these processes, direct validation of the theoretical basis of these models under engine-relevant conditions has not been possible to date. Recent developments in droplet sizing measurement techniques offer a new opportunity to evaluate droplet size distributions formed in the central and peripheral regions of the spray. There is therefore a need to understand how these measurements might be utilized to validate unobservable physics in the near nozzle-region. To address this need, we conduct a computational study using 3D CFD simulations in CONVERGE to explore the relationship between the selected primary atomization model and droplet sizes formed in the central and peripheral regions of the spray. Two existing primary atomization models from the literature are studied to characterize the influence of competing aerodynamics and turbulence mechanisms on the spray formation process. We develop and implement a new hybrid primary atomization model to evaluate the influence of the assumed turbulent scaling on the predicted spray structure. Local sensitivity analysis is performed over a wide range of ambient densities, injection pressures and nozzle diameters to compare the response of predicted droplet sizes in different regions of the spray to changes in injection and ambient conditions. Comparison of the predicted spray structure among the three spray models and against available measurements helps identify a set of experimental conditions and measurements that are needed to inform the development of improved atomization and spray breakup models.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Electronic Structure and Mechanical Properties of Solvated Montmorillonite Clay Using Large-Scale DFT Method

Montmorillonite clay (MMT) has been widely used in engineering and environmental applications as a landfill barrier and toxic waste repository due to its unique property as an expandable clay mineral that can absorb water easily. This absorption process rendered MMT to be highly exothermic due to electrostatic interactions among molecules and hydrogen bonds between surface atoms. A detailed study of a large supercell model of structural clay enables us to predict long-term nuclear waste storage. Herein, a large solvent MMT model with 4071 atoms is studied using ab initio density functional theory. The DFT calculation and analysis clarify the important issues, such as bond strength, solvation effect, elasticity, and seismic wave velocities. These results are compared to our previous study on crystalline MMT (dry). The solvated MMT has reduced shear modulus (G), bulk modulus (K), and Young’s modulus (E). We observe that the conduction band (CB) in the density of states (DOS) of solvated MMT model has a single, conspicuous peak at –8.5 eV. Moreover, the atom-resolved partial density of states (PDOS) summarizes the roles played by each atom in the DOS. These findings illuminate numerous potential sophisticated applications of MMT clay.

36 MATERIALS SCIENCE↗
NASA NTRS2003

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗
DOE OSTI2015

Effects of the local structure dependence of evaporation fields on field evaporation behavior

Accurate three dimensional reconstructions of atomic positions and full quantification of the information contained in atom probe microscopy data rely on understanding the physical processes taking place during field evaporation of atoms from needle-shaped specimens. However, the modeling framework for atom probe microscopy has only limited quantitative justification. Building on the continuum field models previously developed, we introduce a more physical approach with the selection of evaporation events based on density functional theory calculations. This model reproduces key features observed experimentally in terms of sequence of evaporation, evaporation maps, and depth resolution, and provides insights into the physical limit for spatial resolution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Modeling the Kondo effect of a magnetic atom adsorbed on graphene

The Kondo effect due to a magnetic atom adsorbed on graphene is investigated using two theoretical approaches, including the non-equilibrium Green's function method within the slave-boson mean-field approximation and a newly developed tight-binding (TB) model with an effectively infinite Hubbard U term. Both methods reveal the presence of a Kondo peak in the local density of state (LDOS) of graphene near the Fermi level. A sharp Kondo peak is only observed in the odd-neighboring C sites of the C atom directly below the magnetic atom placed in top position. The peak intensity is found to decay quickly with respect to the distance between the C site and the magnetic atom. A theoretical model of scanning tunneling microscopy (STM) of graphene is presented as a means to identify the manifestation of the Kondo effect via direct topographic STM measurements. STM simulations show that in the proximity of the magnetic atom, only the C atoms of one sublattice are visible at low bias potential and thus a triangular lattice can be seen, in agreement with recent experiments. Further away from the magnetic atom, the Kondo effect dies down and eventually both sublattices are visible, leading to the recovery of the hexagonal lattice. Furthermore, this work not only provides a new TB method to study the Kondo effect in graphene, but also presents both scanning tunneling spectroscopy and microscopy fingerprints of the Kondo effect to facilitate its experimental verification.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tuning the Product Selectivity of Single-Atom Alloys by Active Site Modification

There is widespread interest in developing catalysts with uniform active sites that consist of single atoms, thereby simplifying the reaction mechanism and improving product selectivity. Here, we examine experimentally how CO can be used to modify the active sites on a strong binding single-atom alloy and examine how this in turn impacts product selectivity for a reaction that has two different pathways. Specifically, we find that CO can be used to selectively block isolated Rh atom active sites in a RhCu(111) model single-atom alloy catalyst surface and promote the dehydrogenation pathway for adsorbed ethyl groups by suppressing the hydrogenation pathway.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine Learning Predictions of Simulated Self-Diffusion Coefficients for Bulk and Confined Pure Liquids

Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. Here in this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.

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

A NEW GENERATION OF LOS ALAMOS OPACITY TABLES

We present a new, publicly available set of Los Alamos OPLIB opacity tables for the elements hydrogen through zinc. Our tables are computed using the Los Alamos ATOMIC opacity and plasma modeling code, and make use of atomic structure calculations that use fine-structure detail for all the elements considered. Our equation of state model, known as ChemEOS, is based on the minimization of free energy in a chemical picture and appears to be a reasonable and robust approach to determining atomic state populations over a wide range of temperatures and densities. In this paper we discuss in detail the calculations that we have performed for the 30 elements considered, and present some comparisons of our monochromatic opacities with measurements and other opacity codes. We also use our new opacity tables in solar modeling calculations and compare and contrast such modeling with previous work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Electromagnetic model for near-field microwave microscope with atomic resolution: Determination of tunnel junction impedance

An electrodynamic model is proposed for the tunneling microwave microscope with subnanometer space resolution as developed by Lee et al. [Appl. Phys. Lett. 97, 183111 (2010)]. Tip-sample impedance Z{sub a} was introduced and studied in the tunneling and non-tunneling regimes. At tunneling breakdown, the microwave current between probe and sample flows along two parallel channels characterized by impedances Z{sub p} and Z{sub t} that add up to form overall impedance Z{sub a}. Quantity Z{sub p} is the capacitive impedance determined by the near field of the probe and Z{sub t} is the impedance of the tunnel junction. By taking into account the distance dependences of effective tip radius r{sub 0}(z) and tunnel resistance R{sub t}(z) = Re[Z{sub t}(z)], we were able to explain the experimentally observed dependences of resonance frequency f{sub r}(z) and quality factor Q{sub L}(z) of the microscope. The obtained microwave resistance R{sub t}(z) and direct current tunnel resistance R{sub t}{sup dc}(z) exhibit qualitatively similar behavior, although being largely different in both magnitude and the characteristic scale of height dependence. Interpretation of the microwave images of the atomic structure of test samples proved possible by taking into account the inductive component of tunnel impedance ImZ{sub t} = ωL{sub t}. Relation ωL{sub t}/R{sub t} ≈ 0.235 was obtained.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Variational formulation of classical and quantum models for intense laser pulse propagation

Highlights: •We consider the theoretical description of intense laser pulses propagating through gases. •We derived reduced models using Lagrangian and Hamiltonian formulations. •We obtain models where the atomic electrons are treated classically or quantum-mechanically. •We derive a number of conserved quantities of the reduced models. -- Abstract: We consider the theoretical description of intense laser pulses propagating through gases. Starting from a first-principles description of both the electromagnetic field and the electron motion within the gas atoms, we derive a hierarchy of reduced models. We obtain a parallel set of models, where the atomic electrons are treated classically on the one hand, and quantum-mechanically on the other. By working consistently in either a Lagrangian formulation or a Hamiltonian formulation, we ensure that our reduced models preserve the variational structure of the parent models. Taking advantage of the Hamiltonian formulation, we deduce a number of conserved quantities of the reduced models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A new method for obtaining model-free viscoelastic material properties from atomic force microscopy experiments using discrete integral transform techniques

Viscoelastic characterization of materials at the micro- and the nanoscale is commonly performed with the aid of force–distance relationships acquired using atomic force microscopy (AFM). The general strategy for existing methods is to fit the observed material behavior to specific viscoelastic models, such as generalized viscoelastic models or power-law rheology models, among others. Here we propose a new method to invert and obtain the viscoelastic properties of a material through the use of the Z-transform, without using a model. We present the rheological viscoelastic relations in their classical derivation and their z-domain correspondence. We illustrate the proposed technique on a model experiment involving a traditional ramp-shaped force–distance AFM curve, demonstrating good agreement between the viscoelastic characteristics extracted from the simulated experiment and the theoretical expectations. We also provide a path for calculating standard viscoelastic responses from the extracted material characteristics. The new technique based on the Z-transform is complementary to previous model-based viscoelastic analyses and can be advantageous with respect to Fourier techniques due to its generality. Additionally, it can handle the unbounded inputs traditionally used to acquire force–distance relationships in AFM, such as ramp functions, in which the cantilever position is displaced linearly with time for a finite period of time.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Proton scattering by a hydrogen atom in an effectively two-body model

It is assumed that the total potential of proton interaction with a hydrogen atom is the sum of the short-range nuclear soft-core Reid potential and the long-range Thomas-Fermi potential. A quantum mechanical analysis of low-energy features of the phase shift and cross section for elastic proton scattering on a hydrogen atom is given for the case of zero total angular momentum. The calculations performed in the present study within a nonlinear version of the variable-phase approach ultimately revealed that, because of a long-range character of the asymptotic behavior of the Thomas-Fermi potential, the respective cross section at low energies oscillates but has a finite number of zeros.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: III. From atomic nuclei to neutron stars

Here, we present BSkG3, the latest entry in the Brussels-Skyrme-on-a-grid series of large-scale models of nuclear structure based on an energy density functional. Compared to its predecessors, the new model offers a more realistic description of nucleonic matter at the extreme densities relevant to neutron stars. This achievement is made possible by incorporating a constraint on the infinite nuclear matter properties at high densities in the parameter adjustment, ensuring in this way that the predictions of BSkG3 for the nuclear Equation of State are compatible with the observational evidence for heavy pulsars with M > 2M ⊙ . Instead of the usual phenomenological pairing terms, we also employ a more microscopically founded treatment of nucleon pairing, resulting in extrapolations to high densities that are in line with the predictions of advanced many-body methods and are hence more suited to the study of superfluidity in neutron stars. By adopting an extended form of the Skyrme functional, we are able to reconcile the description of matter at high densities and at saturation density: the new model further refines the description of atomic nuclei offered by its predecessors. A qualitative improvement is our inclusion of ground state reflection asymmetry, in addition to the spontaneous breaking of rotational, axial, and time-reversal symmetry. Quantitatively, the model offers lowered root-mean-square deviations on 2457 masses (0.631 MeV), 810 charge radii (0.0237 fm) and an unmatched accuracy with respect to 45 primary fission barriers of actinide nuclei (0.33 MeV). Reconciling the complexity of neutron stars with those of atomic nuclei establishes BSkG3 as a tool of choice for applications to nuclear structure, the nuclear equation of state and nuclear astrophysics in general.

Physics↗