Search NASA⌕ Search

SEARCH · Search NASA

Results for “density model”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2018

The Interplay Between Radiation Pressure and the Photoelectric Instability in Optically Thin Disks of Gas and Dust

In optically thin disks, dust grains are photoelectrically stripped of electrons by starlight, heating nearby gas and possibly creating a dust clumping instability-the photoelectric instability (PeI)-that significantly alters global disk structure. In the current work, we use the Pencil Code to perform the first numerical models of the PeI that include stellar radiation pressure on dust grains in order to explore the parameter regime in which the instability operates. In some models with low gas and dust surface densities, we see a variety of dust structures, including sharp concentric rings. In the most gas- and dust-rich models, nonaxisymmetric clumps, arcs, and spiral arms emerge that represent dust surface density enhancements of factors of ∼5-20. In one high gas surface density model, we include a large, low-order gas viscosity and find that it observably smooths the structures that form in the gas and dust, suggesting that resolved images of a given disk may be useful for deriving constraints on the effective viscosity of its gas. Our models show that radiation pressure does not preclude the formation of complex structure from the PeI, but the qualitative manifestation of the PeI depends strongly on the parameters of the system. The PeI may provide an explanation for unusual disk morphologies, such as the moving blobs of the AU Mic disk, the asymmetric dust distribution of the 49Ceti disk, and the rings and arcs found in the HD 141569A disk.

photoelectric instability (PeI)↗
DOE OSTI2018

Investigation of the Electron Density Variation During the 21 August 2017 Solar Eclipse

Abstract This paper presents a comparison of modeled and measured electron densities for the 21 August 2017 solar eclipse across the USA. The location of the instrument was (43.81°N, 247.32°E) where the maximum obscuration of 99.6% occurred at 17.53 hr UT on 21 August. The solar apparent time was 9.96 hr, and the duration of the eclipse was 2.7 hr. It was found that if it is assumed that there are no chromosphere emissions at totality, ~30% coronal emission remaining at totality gave the best fit to the electron density variation at 150 km. The 30% coronal emission estimate has uncertainties associated with respect to uncertainties in the solar spectrum, the measured electron density, and the amount of chromosphere emissions remaining at totality. The agreement between the modeled and measured electron densities is excellent at 150 km with the assumed 30% coronal emission at totality. At other altitudes, the agreement is very good, but the altitude profile would be improved if the model peak electron density ( N m F 2 ) decayed more slowly to better match the data. The minimum N m F 2 in the model occurs ~10 min after totality when it decreases to 0.55 from its noneclipse value. The minimum of the N m F 2 data occurs between 6 and 10 min after totality but is ~15% larger. The total electron content decreases to 0.65 of its preeclipse value. These relative changes agree well with those predicted by others prior to the eclipse.

Reinisch, B. W.↗
DOE OSTI2018

Quantum Chemical Calculation of p K a s of Environmentally Relevant Functional Groups: Carboxylic Acids, Amines, and Thiols in Aqueous Solution

Developing accurate quantum chemical approaches for calculating p K a s is of broad interest. Useful accuracy can be obtained by using density functional theory (DFT) in combination with a polarizable continuum solvent model. However, some classes of molecules present problems for this approach, yielding errors greater than 5 p K units. Various methods have been developed to improve the accuracy of the combined strategy. These methods perform well but either do not generalize or introduce additional degrees of freedom, increasing the computational cost. The Solvation Model based on Density (SMD) has emerged as one of the most commonly used continuum solvent models. Nevertheless, for some classes of organic compounds, e.g., thiols, the p K a s calculated with the original SMD model show errors of 6–10 p K units, and we traced these errors to inaccuracies in the solvation free energies of the anions. To improve the accuracy of p K a s calculated with DFT and the SMD model, we developed a scaled solvent-accessible surface approach for constructing the solute–solvent boundary. By using a “direct” approach, in which all quantities are computed in the presence of the continuum solvent, the use of thermodynamic cycles is avoided. Furthermore, no explicit water molecules are required. Three benchmark data sets of experimentally measured p K a values, including 28 carboxylic acids, 10 aliphatic amines, and 45 thiols, were used to assess the optimized SMD model, which we call SMD with a scaled solvent-accessible surface (SMD sSAS ). Of the methods tested, the M06-2X density functional approximation, 6-31+G(d,p) basis set, and SMD sSAS solvent model provided the most accurate p K a s for each set, yielding mean unsigned errors of 0.9, 0.4, and 0.5 p K units, respectively, for carboxylic acids, aliphatic amines, and thiols. As a result, this approach is therefore useful for efficiently calculating the p K a s of environmentally relevant functional groups.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Accurate modeling of F-region electron densities

In the past year, we have made considerable progress in a number of areas including algorithm development, completion of two major case studies, and the development of a new EUV flux model. As a result, there has been a major improvement in our ability to model global emissions in support of NASA's imaging plans. Activity highlights include the following: developed a new algorithm to allow physical models to reproduce observed NmF2; investigated the relationship between NmF2 and F10.7 at Millstone Hill during 1990; developed a new solar EUV flux model; statistical survey of anomalously high nighttime electron T(sub e) at Millstone Hill; conducted a case study of the March 1990 magnetic storm; and conducted a comparison between theory and data of magnetically quiet behavior of the winter ionosphere at Millstone Hill.

Source record↗
DOE OSTI2024

Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state

Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Multiple-relaxation-time color-gradient lattice Boltzmann model for simulating two-phase flows with high density ratio

In this paper, we propose a color-gradient lattice Boltzmann (LB) model for simulating two-phase flows with high density ratio and high Reynolds number. The model applies a multi-relaxation-time (MRT) collision operator to enhance the stability of the simulation. A source term, which is derived by the Chapman-Enskog analysis, is added into the MRT LB equation so that the Navier-Stokes equations can be exactly recovered. Also, a new form of the equilibrium density distribution function is used to simplify the source term. To validate the proposed model, steady flows of a static droplet and the layered channel flow are first simulated with density ratios up to 1000. Small values of spurious velocities and interfacial tension errors are found in the static droplet test, and improved profiles of velocity are obtained by the present model in simulating channel flows. Then, two cases of unsteady flows, Rayleigh-Taylor instability and droplet splashing on a thin film, are simulated. In the former case, the density ratio of 3 and Reynolds numbers of 256 and 2048 are considered. The interface shapes and spike/bubble positions are in good agreement with the results of previous studies. In the latter case, the droplet spreading radius is found to obey the power law proposed in previous studies for the density ratio of 100 and Reynolds number up to 500.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

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

Modeling dynamics of acute HIV infection incorporating density-dependent cell death and multiplicity of infection

Understanding the dynamics of acute HIV infection can offer valuable insights into the early stages of viral behavior, potentially helping uncover various aspects of HIV pathogenesis. The standard viral dynamics model explains HIV viral dynamics during acute infection reasonably well. However, the model makes simplifying assumptions, neglecting some aspects of HIV infection. For instance, in the standard model, target cells are infected by a single HIV virion. Yet, cellular multiplicity of infection (MOI) may have considerable effects in pathogenesis and viral evolution. Further, when using the standard model, we take constant infected cell death rates, simplifying the dynamic immune responses. Here, we use four models—1) the standard viral dynamics model, 2) an alternate model incorporating cellular MOI, 3) a model assuming density-dependent death rate of infected cells and 4) a model combining (2) and (3)—to investigate acute infection dynamics in 43 people living with HIV very early after HIV exposure. We find that all models qualitatively describe the data, but none of the tested models is by itself the best to capture different kinds of heterogeneity. Instead, different models describe differing features of the dynamics more accurately. For example, while the standard viral dynamics model may be the most parsimonious across study participants by the corrected Akaike Information Criterion (AICc), we find that viral peaks are better explained by a model allowing for cellular MOI, using a linear regression analysis as analyzed by R 2 . These results suggest that heterogeneity in within-host viral dynamics cannot be captured by a single model. Depending on the specific aspect of interest, a corresponding model should be employed.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1981

A model of temporal variations in ozone density in the Martian atmosphere

Temporal variations of the ozone density profile in the Martian atmosphere at high latitudes are calculated for the course of a Martian year, taking into account seasonal and diurnal variations in temperature, water vapor and solar radiation. Calculations are based on a model including 35 neutral photochemical reactions, and vertical eddy diffusion using a time step of 12 min for the region from the surface to 240 km altitude. Results are found to be in better agreement with Mariner 9 observations of the time and magnitude of the seasonal maximum than previous model calculations. The diurnal variation is predicted to be small near the solstices, with the nighttime ozone density greater than the daytime and the magnitude of the difference dependent on season. Opposite temporal variations are predicted for ozone densities above and below about 25 km, and an ozone density maximum at 35-40 km is obtained. It is suggested that the effects of an aerosol layer may not be important in enhancing predicted ozone concentration, and may even decrease it.

Shimazaki, T.↗
DOE OSTI2018

Improved yield and control of spectra from high-intensity laser-generated neutron beams

Abstract Kinetic modeling of laser-ion beam generation from the “break-out afterburner” (BOA) has been modeled for several deuteron-rich solid-density target foils. Modeling the transport of these beams in a beryllium converter shows as much as a fourfold increase in neutron yield over the present state of the art through the use of alternative target materials. Additionally, species-separation dynamics during the BOA can be exploited to control the hardness of the neutron spectra, of interest for, for example, enhancing penetrability in shielded material in active neutron interrogation settings.

laser based neutron sources, active interrogation↗
DOE OSTI2018

Effects of Electrolyte Buffer Capacity on Surface Reactant Species and the Reaction Rate of CO 2 in Electrochemical CO 2 Reduction

In the aqueous electrochemical reduction of CO 2 , the choice of electrolyte is responsible for the catalytic activity and selectivity, although there remains a need for more in-depth understanding of electrolyte effects and mechanisms. In this study, using both experimental and simulation approaches, we report how the buffer capacity of the electrolytes affects the kinetics and equilibrium of surface reactant species and the resulting reaction rate of CO 2 with varying partial CO 2 pressure. Electrolytes investigated include KCl (nonbuffered), KHCO 3 (buffered by bicarbonate), and phosphate-buffered electrolytes. Assuming 100% methane production, the simulation successfully explains the experimental trends in maximum CO 2 flux in KCl and KHCO 3 and also highlights the difference between KHCO 3 and phosphate in terms of pKaas well as the impact of the buffer capacity. To examine the electrolyte impact on selectivity, the model is run with a constant total current density. Using this model, several factors are elucidated, including the importance of local pH, which is not in acid/base equilibrium, the impact of buffer identity and kinetics, and the mass-transport boundary-layer thickness. The gained understanding can help to optimize CO 2 reduction in aqueous environments.

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

Vector interaction enhanced bag model for astrophysical applications

For quark matter studies in astrophysics the thermodynamic bag model (tdBAG) has been widely used. Despite its success it fails to account for various phenomena expected from QCD. We suggest a straightforward extension of tdBAG in order to take the dynamical breaking of chiral symmetry and the influence of vector interactions explicitly into account. As for tdBAG the model mimics confinement in a phenomenological approach. It is based on an analysis of the Nambu–Jona-Lasinio (NJL) model at finite density. Furthermore, we demonstrate how NJL and bag models in this regime follow from the more general and QCD-based framework of the Dyson–Schwinger equations in a medium by assuming simple gluon contact interaction. Based on our simple and novel model, we construct quark hadron hybrid equations of state and systematically study chiral and deconfinement phase transitions, the appearance of s-quarks, and the role of vector interaction. We further study these aspects for matter in β-equilibrium at zero temperature, with particular focus on the current ∼2 M{sub ⊙} maximum mass constraint for neutron stars. Our approach indicates that the currently only theoretical evidence for the hypothesis of stable strange matter is an artifact of tdBAG and results from neglecting the dynamical breaking of chiral symmetry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Microstructure-Sensitive Uncertainty Quantification for Crystal Plasticity Finite Element Constitutive Models Using Stochastic Collocation Methods

Uncertainty quantification (UQ) plays a major role in verification and validation for computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allows numerical predictions, providing the bridge from microstructure to materials properties and performances. Several constitutive models have been proposed in the last few decades to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method, which is mathematically rigorous and has been widely used in the field of UQ, to quantify the uncertainty of three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curve, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Wave spectra of a shoaling wave field: A comparison of experimental and simulated results

Wave profile measurements made from an aircraft crossing the North Carolina continental shelf after passage of Tropical Storm Amy in 1975 are used to compute a series of wave energy spectra for comparison with simulated spectra. Results indicate that the observed wave field experiences refraction and shoaling effects causing statistically significant changes in the spectral density levels. A modeling technique is used to simulate the spectral density levels. Total energy levels of the simulated spectra are within 20 percent of those of the observed wave field. The results represent a successful attempt to theoretically simulate, at oceanic scales, the decay of a wave field which contains significant wave energies from deepwater through shoaling conditions.

Morris, W. D.↗
NASA NTRS1992

Helium in the polar wind

The coupled time-dependent continuity, momentum, and energy equations for four charged species (H(+), He(+), O(+) and electrons) moving through a stationary upper neutral atmosphere (composed of O2, N2, O, He, and H) were solved along open geomagnetic field lines using a modified version of the Gombosi et al. (1985) model. The thermospheric densities and temperatures, obtained from the MSIS-86 model (Hedin, 1987), were used to calculate neutral-ion sources, sinks, and collision frequencies. The various cross sections and conductivities employed by the Gombosi et al. model were also updated. The paper presents the first time-dependent solutions for the He(+) component of the polar wind. Special attention is paid to the solar cycle variation of the ion upflows.

Gombosi, T. I.↗