Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
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 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.↗
DOE OSTI2020

Investigating effects of soil chemicals on density of small mammal bioindicators using spatial capture-recapture models

Monitoring the ecological impacts of environmental pollution and the effectiveness of remediation efforts requires identifying relationships between contaminants and the disruption of biological processes in populations, communities, or ecosystems. Wildlife are useful bioindicators, but traditional comparative experimental approaches rely on a staunch and typically unverifiable assumption that, in the absence of contaminants, reference and contaminated sites would support the same densities of bioindicators, thereby inferring direct causation from indirect data. We demonstrate the utility of spatial capture-recapture (SCR) models for overcoming these issues, testing if community density of common small mammal bioindicators was directly influenced by soil chemical concentrations. By modeling density as an inhomogeneous Poisson point process, we found evidence for an inverse spatial relationship between Peromyscus density and soil mercury concentrations, but not other chemicals, such as polychlorinated biphenyls, at a site formerly occupied by a nuclear reactor. Although the coefficient point estimate supported Peromyscus density being lower where mercury concentrations were higher (β = –0.44), the 95% confidence interval overlapped zero, suggesting no effect was also compatible with our data. Estimated density from the most parsimonious model (2.88 mice/ha; 95% CI = 1.63–5.08), which did not support a density-chemical relationship, was within the range of reported densities for Peromyscus that did not inhabit contaminated sites elsewhere. Environmental pollution remains a global threat to biodiversity and ecosystem and human health, and our study provides an illustrative example of the utility of SCR models for investigating the effects that chemicals may have on wildlife bioindicator populations and communities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Feedback density and causal complexity of simulation model structure

Measures of simulation model complexity generally focus on outputs; we propose measuring the complexity of a model’s causal structure to gain insight into its fundamental character. This article introduces tools for measuring causal complexity. First, we introduce a method for developing a model’s causal structure diagram, which characterises the causal interactions present in the code. Causal structure diagrams facilitate comparison of simulation models, including those from different paradigms. Next, we develop metrics for evaluating a model’s causal complexity using its causal structure diagram. We discuss cyclomatic complexity as a measure of the intricacy of causal structure and introduce two new metrics that incorporate the concept of feedback, a fundamental component of causal structure. The first new metric introduced here is feedback density, a measure of the cycle-based interconnectedness of causal structure. The second metric combines cyclomatic complexity and feedback density into a comprehensive causal complexity measure. Finally, we demonstrate these complexity metrics on simulation models from multiple paradigms and discuss potential uses and interpretations. These tools enable direct comparison of models across paradigms and provide a mechanism for measuring and discussing complexity based on a model’s fundamental assumptions and design.

97 MATHEMATICS AND COMPUTING↗