Search NASA⌕ Search

SEARCH · Search NASA

Results for “agreement dynamics”

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

Detailed study of the eikonal reaction theory for the breakup of one-neutron halo nuclei

Background: One-neutron removal reactions are used to study the single-particle structure of unstable nuclei, and in particular the exotic halo nuclei. The eikonal reaction theory (ERT) has been developed by Yahiro, Ogata, and Minomo [Prog. Theor. Phys. 126, 167 (2011)] to include dynamical effects, which are missing in the usual eikonal description of these reactions. Encouraging results have been obtained for total breakup cross sections in comparison to more elaborate reaction models. Purpose: Here, we extend these comparisons to more differential breakup cross sections expressed as functions of the relative energy or parallel momentum between the core and halo neutron. Method: ERT predictions of these cross sections are compared to state-of-the-art calculations. Results: The hypotheses upon which the ERT is based are confirmed and their range of validity is made clearer. The actual ordering of the evolution operators affects ERT differential cross sections and a specific choice leads to excellent agreement with the reference calculation. Dynamical effects in the treatment of the neutron-target interaction can be significant in the parallel-momentum observable. Conclusions: The role of the different interactions in the dynamics of breakup reactions of one-neutron halo nuclei are better understood and improvements to the ERT are suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

Wintertime nitric acid chemistry - Implications from three-dimensional model calculations

A three-dimensional simulation of the evolution of HNO3 has been run for the winter of 1979. Winds and temperatures are taken from a stratospheric data assimilation analysis, and the chemistry is based on Limb Infrared Monitor of the Stratosphere (LIMS) observations. The model is compared to LIMS observations to investigate the problem of 'missing' nitric acid chemistry in the winter hemisphere. Both the model and observations support the contention that a nitric acid source is needed outside of the polar vortex and north of the subtropics. Observations suggest that HNO3 is not dynamically controlled in middle latitudes. The model shows that given the time scales of conventional chemistry, dynamical control is expected. Therefore, an error exists in the conventional chemistry or additional processes are needed to bring the model and data into agreement. Since the polar vortex is dynamically isolated from the middle latitudes, and since the highest HNO3 values are observed in October and November, a source associated solely with polar stratospheric clouds cannot explain the deficiencies in the chemistry. The role of heterogeneous processes on background aerosols is reviewed in light of these results.

Rood, Richard B.↗
DOE OSTI2020

Nature of the nonequilibrium phase transition in the non-Markovian driven Dicke model

The Dicke model famously exhibits a phase transition to a superradiant phase with a macroscopic population of photons and is realized in multiple settings in open quantum systems. In this paper, we study a variant of the Dicke model where the cavity mode is lossy due to the coupling to a Markovian environment while the atomic mode is coupled to a colored bath. We analytically investigate this model by inspecting its low-frequency behavior via the Schwinger-Keldysh field theory and carefully examine the nature of the corresponding superradiant phase transition. Integrating out the fast modes, we can identify a simple effective theory allowing us to derive analytical expressions for various critical exponents including the dynamical exponent. We find excellent agreement with previous numerical results when the non-Markovian bath is at zero temperature; however, contrary to these studies, our low-frequency approach reveals that the same exponents govern the critical behavior when the colored bath is at finite temperature unless the chemical potential is zero. Furthermore, we show that the superradiant phase transition is classical in nature, while it is genuinely nonequilibrium. In this work, we derive a fractional Langevin equation and conjecture the associated fractional Fokker-Planck equation that captures the system's long-time memory as well as its nonequilibrium behavior. Finally, we consider finite-size effects at the phase transition and identify the finite-size scaling exponents, unlocking a rich behavior in both statics and dynamics of the photonic and atomic observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dynamics and Structure of Molten Neodymium Titanate

Binary rare earth titanates can be vitrified using levitation melt quenching, and their glasses are promising for optical applications. To better understand their thermophysical properties and dynamics, we have used molecular dynamics (MD) and empirical potential structure refinement (EPSR) to simulate molten neodymium titanate, 83TiO 2 –17Nd 2 O 3 . The models are compared with high-energy synchrotron X-ray diffraction of the melt structure over 1060–1960 K and thermophysical property measurements of electrostatically levitated droplets over 1500–2200 K. As temperature increases, the network periodicity increases and the correlation length decreases, both in monotonic trends that suggest continuous structural evolution through the melt and supercooled liquid regimes. At 1900 K, the mean Ti–O and Nd–O coordination numbers are 5.8(5) and 6.7(5), compared to the glass, 5.72(6) and 7.70(26). The equilibrium melting point is measured as 1860(20) K, over 100 K higher than previous estimates. Viscosity data from MD are in good agreement with experimental measurements. The dynamic propensity of Nd is similar to that of the network-forming Ti, suggesting that Nd does not participate as a conventional modifier. The viscosity of this fragile liquid is likely controlled by the breaking of cation–oxygen bonds and rearrangements of Ti–O polyhedra.

molecular dynamics↗
DOE OSTI2023

Extension of Argonne’s Plant Dynamics Code (PDC) Capabilities to TerraPower Direct-Cycle Supercritical CO2 (“Pascal”) Nuclear Reactor Concept. Final CRADA report

Argonne National Laboratory (the Contractor) and TerraPower, LLC, (the Participant), entered into a Cooperative Research and Development Agreement (CRADA) to extend the Plant Dynamics Code (PDC) to model the Pascal reactor plant. PDC is the system-level computer code developed at Argonne, primarily for design and transient analysis of supercritical carbon dioxide (sCO2) energy conversion systems for various power plants. Pascal is a supercritical carbon dioxide (sCO2)-cooled direct-cycle nuclear power reactor under development by TerraPower. New features were added to the code to model components and systems specific to the Pascal plant and to make the code readily available as a tool for design, control, and safety analysis for this and other direct-cycle power plants. TerraPower provided design information of the Pascal reactor and the requirements for the development of the new PDC features. Argonne modified the PDC capabilities to introduce the new features.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2010

Simulation and 'TWINS Observations of the 22 July 2009 Storm

TWINS is the first mission to perform stereo imaging of the Earth's ring current. The magnetic storm on 22 July 2009 is the largest storm observed since TWINS began routine stereo imaging in June 2008. On 22 July 2009, the Dst dropped to nearly -80nT at 7:00 and 10:00 UT. During the main phase and at the peak of the storm, TWINS 1 and 2 were near apogee and moving from pre-dawn to post-dawn local time. The energetic neutral atom (ENA) imagers on the 2 spacecraft captured the storm intensification and the formation of the partial ring current. The peak of the ENA emissions was seen in the midnight-to-dawn local-time sector. The development of this storm has been simulated using the Comprehensive Ring Current Model (CRCM) to understand and interpret the observed signatures. We perform CRCM runs with constant and time-varying magnetic field. The model calculations are validated by comparing the simulated ENA and ion flux intensities with TWINS ENA images and in-situ ion data from THEMIS satellites. Simulation with static magnetic field produces a strong shielding electric field that skews the ion drift trajectories toward dawn. The model's corresponding peak ENA emissions are always eastward than those in the observed TWINS images. On the other hand, simulation with a dynamic magnetic field gives better spatial agreements with both ENA and insitu particle data, suggesting that temporal variations of the geomagnetic field exert a significant influence upon global ring current ion dynamics.

Fok, Mei-Ching↗
DOE OSTI2020

Atomic fluctuations in electronic materials revealed by dephasing

The microscopic origin and timescale of the fluctuations of the energies of electronic states has a significant impact on the properties of interest of electronic materials, with implication in fields ranging from photovoltaic devices to quantum information processing. Spectroscopic investigations of coherent dynamics provide a direct measurement of electronic fluctuations. Modern multidimensional spectroscopy techniques allow the mapping of coherent processes along multiple time or frequency axes and thus allow unprecedented discrimination between different sources of electronic dephasing. Exploiting modern abilities in coherence mapping in both amplitude and phase, we unravel dissipative processes of electronic coherences in the model system of CdSe quantum dots (QDs). The method allows the assignment of the nature of the observed coherence as vibrational or electronic. The expected coherence maps are obtained for the coherent longitudinal optical (LO) phonon, which serves as an internal standard and confirms the sensitivity of the technique. Fast dephasing is observed between the first two exciton states, despite their shared electron state and common environment. This result is contrary to predictions of the standard effective mass model for these materials, in which the exciton levels are strongly correlated through a common size dependence. In contrast, the experiment is in agreement with ab initio molecular dynamics of a single QD. Electronic dephasing in these materials is thus dominated by the realistic electronic structure arising from fluctuations at the atomic level rather than static size distribution. The analysis of electronic dephasing thereby uniquely enables the study of electronic fluctuations in complex materials.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Time dependent effects at the foot of the chromosphere

A dynamical model was proposed which is in agreement with the observations that suggest the solar chromosphere is heated preferentially in magnetic regions. The dynamical model requires an energy flux that is large enough to lead to amplitudes for which the apparent temperatures match the empirical line temperatures in the UTV lines. Consistency with the observations of chromospheric emission then requires that the emission take place not uniformly over the solar surface but in smaller spatial regions.

Kalkofen, W.↗
NASA NTRS2011

Simulation and Twins Observations of the 22 July 2009 Storm

TWINS is the first mission to perform stereo imaging of the Earth's ring current. The magnetic storm on 22 July 2009 was at the time the largest storm observed since TWINS began routine stereo imaging in June 2008. On 22 July 2009, the Dst dropped to nearly .80 nT at 0700 and 1000 UT. During the main phase, and at the peak of the storm, TWINS 1 and 2 were near apogee and moving between predawn and postdawn local time. The energetic neutral atom (ENA) imagers on the two spacecraft captured the storm intensification and the formation of the partial ring current. The peak of the high-altitude ENA emissions was seen in the midnight-to-dawn local time sector. The development of this storm has been simulated using the comprehensive ring current model (CRCM) to understand and interpret the observed signatures. We perform CRCM runs with constant and time-varying magnetic field. The model calculations are validated by comparing the simulated ENA and ion flux intensities with TWINS ENA images and in situ ion data from a THEMIS satellite. Simulation with a static magnetic field produces a strong shielding electric field that skews the ion drift trajectories toward dawn. The model's corresponding peak ENA emissions are always more eastward than those in the observed TWINS images. On the other hand, the simulation with a dynamic magnetic field gives better spatial agreement with both ENA and in situ particle data, suggesting that temporal variations of the geomagnetic field exert a significant influence upon global ring current ion dynamics.

Fok, M.-C.↗
NASA NTRS1992

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗
DOE OSTI2022

Extracting Inelastic Scattering Cross Sections for Finite and Aperiodic Materials from Electronic Dynamics Simulations

Explicit time-dependent electronic structure theory methods are increasingly prevalent in the areas of condensed matter physics and quantum chemistry, with the broad-band optical absorptivity of molecular and small condensed-phase systems nowadays routinely studied with such approaches. Here, in this paper, it is demonstrated that electronic dynamics simulations can similarly be employed to study cross sections for the scattering-induced electronic excitations probed in nonresonant inelastic X-ray scattering and momentum-resolved electron energy loss spectroscopies. A method is put forth for evaluating the electronic dynamic structure factor, which involves the application of a momentum boost-type perturbation and transformation of the resulting reciprocal space density fluctuations into the frequency domain. Good agreement is first demonstrated between the dynamic structure factor extracted from these electronic dynamics simulations and the corresponding transition matrix elements from linear response theory. The method is then applied to some extended (quasi)one-dimensional systems, for which the wave vector becomes a good quantum number in the thermodynamic limit. Finally, the dispersion of many-body excitations in a series of hydrogen-terminated graphene flakes (and twisted bilayers thereof) is investigated to highlight the utility of the presented approach for capturing morphology-dependent effects in the inelastic scattering cross sections of nanostructured and/or noncrystalline materials.

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

High-temperature active oxidation of nanocrystalline silicon-carbide: A reactive force-field molecular dynamics study

Flexible woven SiC ceramics are prone to accelerated fiber embrittlement under high temperature oxidation in dynamic oxygen environments. The nanocrystalline structure of the constituent fibers impacts the reaction kinetics and phase transformations during active oxidation. However, fundamental understanding and quantification of grain boundary effects on oxidation behavior in nanocrystalline SiC remain elusive when temperatures exceed 1500 K. This study deploys large-scale molecular dynamics simulations with a reactive force-field to elucidate the complex roles of atomic oxygen reservoir conditions and grain size on oxidation kinetics and the nature of oxides produced in both monocrystalline and nanocrystalline 3C-SiC between 1100 K and 2000 K. The simulations with dynamically replenished oxygen provide good agreement with oxidation kinetics and activation energies for the monocrystalline Si(100) and C(100) orientations published in the available literature. This study reveals that, by contrast, nanocrystalline SiC samples exhibit two distinct oxidation kinetics with a transition point at 1500 K due to surface melting, which is supported by experimental evidence. The introduction of a grain-boundary network produces a two-fold decrease in oxidation activation energies compared to monocrystalline SiC below 1500 K. Above 1500 K, however, the activation energies rise substantially due to the formation of a liquid Si phase at the SiC/Si oxide interface. Further, it is shown that the stability of the interfacial liquid phase is promoted by incoherent grain boundaries in the crystalline SiC. These findings are important for the deployment of nanocrystalline SiC fibers in advanced thermal protection systems for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

QUIC-fire: A fast-running simulation tool for prescribed fire planning

Coupled fire-atmospheric modeling tools are increasingly used to understand the complex and dynamic behavior of wildland fires. Multiple research tools linking combustion to fluid flow use Navier-Stokes numerical solutions coupled to a thermodynamic model to understand fire-atmospheric feedbacks, but these computational fluid dynamics approaches require high-performance computing resources. We present a new simulation tool called QUIC-Fire to rapidly solve these feedbacks by coupling the mature 3-D rapid wind solver QUIC-URB to a physics-based cellular automata fire spread model Fire-CA. QUIC-Fire uses 3-D fuels inputs similar to the CFD model FIRETEC, allowing this tool to simulate effects of fuel structure on local winds and fire behavior. Results comparing fire behavior metrics to the computational fluid dynamic model FIRETEC show strong agreement. As a result, QUIC-Fire is the first tool intended to provide an opportunity for prescribed fire planners to compare, evaluate and design burn plans, including complex ignition patterns and coupled fire atmospheric feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Load-Displacement Relation and Gap Distribution Between Rough Surfaces: Partial Differential Equations Approach

We develop a theoretical model to predict the load-displacement relation and probability density function for gaps between contacting rough surfaces. We derive partial differential equations from the previous model (Joe et al., 2018), and extend them to the non-adhesive contact problem. The predictions of the present theory are compared with numerical results using the Green's function molecular dynamics algorithm, and a good agreement is obtained.

adhesion↗
DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗
DOE OSTI2024

Spinodal enhancement of fluctuations in nucleus-nucleus collisions

Subensemble Acceptance Method (SAM) [1, 2] is an essential link between measured event-by-event fluctuations and their grand canonical theoretical predictions such as lattice QCD. The method allows quantifying the global conservation law effects in fluctuations. In its basic formulation, SAM requires a sufficiently large system such as created in central nucleus-nucleus collisions and sufficient space-momentum correlations. Directly in the spinodal region of the First Order Phase Transition (FOPT) different approximations should be used that account for finite size effects. Thus, we present the generalization of SAM applicable in both the pure phases, metastable and unstable regions of the phase diagram [3]. Obtained analytic formulas indicate the enhancement of fluctuations due to crossing the spinodal region of FOPT and are tested using molecular dynamics simulations. A rather good agreement is observed. Using transport model calculations with interaction potential we show that the spinodal enhancement of fluctuations survives till the later stages of collision via the memory effect [4]. However, at low collision energies the space-momentum correlation is not strong enough for this signal to be transferred to second and third order cumulants measured in momentum subspace. This result agrees well with recent HADES data on proton number fluctuations at $\sqrt{S_{NN}}$ = 2.4 GeV which are found to be consistent with the binomial momentum space acceptance [5].

Poberezhnyuk, Roman↗
DOE OSTI2024

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗