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 1,405 records · Page 78

On the TFNS Subgrid Models for Liquid-Fueled Turbulent Combustion

This paper describes the time-filtered Navier-Stokes (TFNS) approach capable of capturing unsteady flow structures important for turbulent mixing in the combustion chamber and two different subgrid models used to emulate the major processes occurring in the turbulence-chemistry interaction. These two subgrid models are termed as LEM-like model and EUPDF-like model (Eulerian probability density function), respectively. Two-phase turbulent combustion in a single-element lean-direct-injection (LDI) combustor is calculated by employing the TFNS/LEM-like approach as well as the TFNS/EUPDF-like approach. Results obtained from the TFNS approach employing these two different subgrid models are compared with each other, along with the experimental data, followed by more detailed comparison between the results of an updated calculation using the TFNS/LEM-like model and the experimental data.

Liu, Nan-Suey↗

Exploring Standing and Reflected Slow-mode Waves in Flaring Coronal Loops: A Parametric Study Using 2.5D MHD Modeling

Abstract Recent observations of reflected propagating and standing slow-mode waves in hot flaring coronal loops have spurred our investigation into their underlying excitation and damping mechanisms. To understand these processes, we conduct 2.5D magnetohydrodynamic (MHD) simulations using an arcade active region model that includes a hot and dense loop. Our simulations allow for in-depth parametric investigations complementing and expanding our previous 3D MHD modeling results. We excite these waves using a large-amplitude, flow pulse applied at one footpoint of the loop in two distinct models as motivated by observations from the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). The first model (Model 1) incorporates classical compressive viscosity coefficient, while the second model (Model 2) adopts a 10-times enhanced viscosity coefficient. We obtain the following major results: (1) Our 2.5D MHD simulations reinforce previous conclusions derived from 1D and 3D MHD models that significantly enhanced viscosity is crucial for the rapid excitation of standing slow waves with damping times consistent with observations by Wang et al. (2015). (2)We uncover that nonlinearity in Model 1 delays the conversion of a reflected propagating wave into a standing wave. In contrast, Model 2 exhibits a much weak influence of nonlinearity on the excitation time of standing waves, thanks to the suppression of these effects by enhanced viscosity. (3) Our results reveal that the transverse temperature structure holds more influence on wave behavior than the density structure. In Model 1, increased loop temperature contrast significantly enhances wave trapping within the structure, mitigating the impact of temperature-dependent viscous damping. Conversely, in Model 2, the impact of temperature structure on wave behavior weakens in comparison to the effect of viscosity. (4) Model 1 displays evident nonlinear coupling to the fast and kink magnetoacoustic waves and pronounced wave leakage into the corona. Model 2 exhibits significantly weaker effects in this regard. Analyzing three observed wave events by SDO/AIA aligns with Model 2 predictions, providing further support for the substantial viscosity increase. Our 2.5D study unravels the complex interplay of wave-flow phenomena and nonlinear processes in coronal loops, extending our previous 1D modeling results to incorporate more realistic loop geometry. This provides insights into scenarios where 3D effects may be neglected, thereby enhancing our understanding of the intricate dynamics of the solar corona.

Tongjiang Wang↗

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The structure and electronic properties of tetrahedrally bonded hydrogenated amorphous carbon

In this work, we have synthesized hydrogenated and deuterated amorphous carbon materials that have a density, 2.7 ± 0.1 g/cm 3 , consistent with almost entirely tetrahedral bonding. In hydrogen-free tetrahedral amorphous carbon, the presence of a minority of sp2 bonded atoms leads to localized states that could be passivated with hydrogen by analogy with hydrogenated amorphous silicon. Neutron diffraction analysis demonstrated that the local bonding environment is consistent with ab initio models of high density hydrogenated tetrahedral amorphous carbon and with the related tetrahedral molecular structure neopentane. The optical bandgap of our material, 4.5 eV, is close to the bandgap in the density of states determined by scanning tunneling spectroscopy (4.3 eV). This bandgap is considerably larger than that of hydrogen-free tetrahedral amorphous carbon, confirming that passivation of sp 2 associated tail-states has occurred. Both the structural and electronic measurements are consistent with a model in which the tetrahedrally bonded carbon regions are terminated by hydrogen, causing hopping conductivity to dominate.

36 MATERIALS SCIENCE↗

Oxygen magnetic polarization, nodes in spin density, and zigzag spin order in oxides

Recent studies on Ba 2 CoO 4 (BCO) and SrRuO 3 (SRO) have unveiled a variety of intriguing phenomena, such as magnetic polarization on oxygens, unexpected nodes in the spin density profile along bonds, and unusual zigzag spin patterns in triangular lattices. Here, using simple model calculations supplemented by density functional theory we explain the presence of nodes based on the antibonding character of the dominant singly occupied molecular orbitals along the transition metal (TM) to oxygen bonds. Our simple model also allows us to explain the net polarization on oxygen as originated from the hybridization between atoms and mobility of the electrons with spins opposite to those of the closest TM atoms. Overall, our results are not limited to BCO and SRO, but they are generic and qualitatively predict the net polarization expected on any ligands, according to the spin order of the closest TM atoms and the number of intermediate ligand atoms. Finally, we propose that a robust easy-axis anisotropy would suppress the competing 120° antiferromagnetic order to stabilize the zigzag pattern order as ground state in a triangular lattice. Our generic predictions should be applicable to any other compound with characteristics similar to those of BCO and SRO.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Community Land Model Version 5: Description of New Features, Benchmarking, and Impact of Forcing Uncertainty

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

58 GEOSCIENCES↗

The last Balmer line and H gamma in model B stars

Profiles of the last few Balmer lines (H12-H27) and the full profiles of the H gamma line were computed for a grid of hydrogen line blanketed model atmospheres. The use of the last hydrogen line and the equivalent width of H gamma provide a quick means of estimating the effective temperature and surface gravity of B stars. Regression curves relating the quantum number of the last hydrogen line to the electron density in the model atmosphere are given. Futhermore, it is shown that the difference in the gravity determination deduced from H gamma equivalent widths using the Edmonds-Schuter-Wells formalism and the Griem wing formulae is, at worst, 0.18 in the common logarithm.

Fischel, D.↗

No need for dark matter in galaxy clusters within Galileon theory

Modified gravity theories with a screening mechanism have acquired much interest recently in the quest for a viable alternative to General Relativity on cosmological scales, given their intrinsic property of being able to pass Solar System scale tests and, at the same time, to possibly drive universe acceleration on much larger scales. Here, we explore the possibility that the same screening mechanism, or its breaking at a certain astrophysical scale, might be responsible of those gravitational effects which, in the context of general relativity, are generally attributed to Dark Matter. We consider a recently proposed extension of covariant Galileon models in the so-called ''beyond Horndeski'' scenario, where a breaking of the Vainshtein mechanism is possible and, thus, some peculiar observational signatures should be detectable and make it distinguishable from general relativity. We apply this model to a sample of clusters of galaxies observed under the CLASH survey, using both new data from gravitational lensing events and archival data from X-ray intra-cluster hot gas observations. In particular, we use the latter to model the gas density, and then use it as the only ingredient in the matter clusters' budget to calculate the expected lensing convergence map. Results show that, in the context of this extended Galileon, the assumption of having only gas and no Dark Matter at all in the clusters is able to match observations. We also obtain narrow and very interesting bounds on the parameters which characterize this model. In particular, we find that, at least for one of them, the general relativity limit is excluded at 2σ confidence level, thus making this model clearly statistically different and competitive with respect to general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗

Incorporating Beam Attenuation Calculations into an Integrated Data Analysis Model for Ion Effective Charge

Here, we have created a forward model for charge-exchange impurity density measurements that incorporates neutral beam attenuation measurements self-consistently for the purpose of determining the ion-effective charge Z eff . The model is constructed within an integrated data analysis framework to include a self-consistent calculation of neutral beam attenuation due to multiple impurity species into the measurement of a single impurity density. The model includes measurements of the beam Doppler-shift spectrum and shine-through particle flux to determine the neutral beam particle density which is attenuated by ion collisions. Synthetic data are generated from the diagnostic forward model using statistical and calibration uncertainties. These “noisy” data are used in the analysis to evaluate how accurately Z eff is determined. Methods of experimental design are employed to calculate the information gained from different diagnostic combinations. The analysis shows that while attenuation measurements alone do not provide a unique impurity density measurement in the case of a multispecies inhomogeneous plasma, they do provide an effective measurement of the Z eff profile and place constraints on the impurity density profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear triggered lightning models for use in finite difference calculations

Two nonlinear triggered lightning models have been developed for use in finite difference calculations. Both are based on three species of air chemistry physics and couple nonlinearly calculated air conductivity to Maxwell's equations. The first model is suitable for use in three-dimensional modeling and has been applied to the analysis of triggered lightning on the NASA F106B Thunderstorm Research Aircraft. The model calculates number densities of positive ions, negative ions, and electrons as a function of time and space through continuity equations, including convective derivative terms. The set of equations is closed by using experimentally determined mobilities, and the mobilities are also used to determine the air conductivity. Results from the model's application to the F106B are shown. The second model is two-dimensional and incorporates an enhanced air chemistry formulation. Momentum conservation equations replace the mobility assumption of the first model. Energy conservation equations for neutrals, heavy ions, and electrons are also used. Energy transfer into molecular vibrational modes is accounted for. The purpose for the enhanced model is to include the effects of temperature into the air breakdown, a necessary step if the model is to simulate more than the very earliest stages of breakdown. Therefore, the model also incorporates a temperature-dependent electron avalanche rate. Results from the model's application to breakdown around a conducting ellipsoid placed in an electric field are shown.

Rudolph, Terence↗

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↗

Assessing Huygens Probe Entry, Descent, and Landing at Titan Simulation using Dragonfly Atmosphere Model

Dragonfly is a New Frontiers Program mission that will deliver a rotorcraft to Saturn's moon, Titan. This mission follows Huygens as the previous mission that successfully landed a vehicle on Titan. A flight mechanics simulation of Dragonfly's Entry, Descent, and Landing sequence has been developed using the Program to Optimize Simulated Trajectories II. The simulation incorporates several subsystem models, including aerodynamics, gravity, and mass properties, to fully capture the multi-body six degree of freedom dynamics. Among all the subsystem models that inform the Entry, Descent, and Landing dynamics, the atmosphere model of Titan is a critical component. The atmosphere model characterizes the density, temperature, pressure, and winds that the entry vehicle experiences during the descent. This impacts several aspects of the descent such as the peak heating, aerodynamics, parachute release conditions, the dynamics of the vehicle and parachutes, and the landing ellipse. In the course of developing Dragonfly, an updated model of the Titan atmosphere has been created corresponding to Dragonfly's arrival in the mid-2030s, approximately one Titan year after the Huygens mission successfully landed a probe on Titan. This work leverages previous work done to investigate Huygens EDL sequence to assess the atmosphere model developed for Dragonfly. This is done by utilizing the updated Titan atmosphere model, the Dragonfly atmosphere model, within the Huygens POST2-based flight simulation with the goal of characterizing the differences between the atmospheric models and assessing how the Dragonfly atmosphere model impacts Huygens entry dynamics.

Huygens↗

Assessing Huygens Probe Entry, Descent, and Landing at Titan Simulation using Dragonfly Atmosphere Model

Dragonfly is a New Frontiers Program mission that will deliver a rotorcraft to Saturn's moon, Titan. This mission follows Huygens as the previous mission that successfully landed a vehicle on Titan. A flight mechanics simulation of Dragonfly's Entry, Descent, and Landing sequence has been developed using the Program to Optimize Simulated Trajectories II. The simulation incorporates several subsystem models, including aerodynamics, gravity, and mass properties, to fully capture the multi-body six degree of freedom dynamics. Among all the subsystem models that inform the Entry, Descent, and Landing dynamics, the atmosphere model of Titan is a critical component. The atmosphere model characterizes the density, temperature, pressure, and winds that the entry vehicle experiences during the descent. This impacts several aspects of the descent such as the peak heating, aerodynamics, parachute release conditions, the dynamics of the vehicle and parachutes, and the landing ellipse. In the course of developing Dragonfly, an updated model of the Titan atmosphere has been created corresponding to Dragonfly's arrival in the mid-2030s, approximately one Titan year after the Huygens mission successfully landed a probe on Titan. This work leverages previous work done to investigate Huygens EDL sequence to assess the atmosphere model developed for Dragonfly. This is done by utilizing the updated Titan atmosphere model, the Dragonfly atmosphere model, within the Huygens POST2-based flight simulation with the goal of characterizing the differences between the atmospheric models and assessing how the Dragonfly atmosphere model impacts Huygens entry dynamics.

Entry Descent Landing↗

Probabilities and statistics for backscatter estimates obtained by a scatterometer

Methods for the recovery of winds near the surface of the ocean from measurements of the normalized radar backscattering cross section must recognize and make use of the statistics (i.e., the sampling variability) of the backscatter measurements. Radar backscatter values from a scatterometer are random variables with expected values given by a model. A model relates backscatter to properties of the waves on the ocean, which are in turn generated by the winds in the atmospheric marine boundary layer. The effective wind speed and direction at a known height for a neutrally stratified atmosphere are the values to be recovered from the model. The probability density function for the backscatter values is a normal probability distribution with the notable feature that the variance is a known function of the expected value. The sources of signal variability, the effects of this variability on the wind speed estimation, and criteria for the acceptance or rejection of models are discussed. A modified maximum likelihood method for estimating wind vectors is described. Ways to make corrections for the kinds of errors found for the Seasat SASS model function are described, and applications to a new scatterometer are given.

Pierson, Willard J., Jr.↗

Modeling of the reactant conversion rate in a turbulent shear flow

Results are presented of direct numerical simulations (DNS) of spatially developing shear flows under the influence of infinitely fast chemical reactions of the type A + B yields Products. The simulation results are used to construct the compositional structure of the scalar field in a statistical manner. The results of this statistical analysis indicate that the use of a Beta density for the probability density function (PDF) of an appropriate Shvab-Zeldovich mixture fraction provides a very good estimate of the limiting bounds of the reactant conversion rate within the shear layer. This provides a strong justification for the implementation of this density in practical modeling of non-homogeneous turbulent reacting flows. However, the validity of the model cannot be generalized for predictions of higher order statistical quantities. A closed form analytical expression is presented for predicting the maximum rate of reactant conversion in non-homogeneous reacting turbulence.

Frankel, S. H.↗

Density and stress distribution in the moon.

A model is presented for the lateral variations of density within the moon. The model gives rise to a gravitational potential which is equal to the observed potential at the lunar surface; moreover, it minimizes the total shear-strain energy of the moon. The model exhibits lateral variations of about plus or minus 0.25 g/cc within 50 km depth. The variations, however, reduce to plus or minus 0.06 and plus or minus 0.008 g/cc within layers at 50 to 135 and 135 to 235 km respectively, and they become negligible below this region. The associated stress differences are found to be about 50 bar within 600 km depth, having their maximum values of about 90 bars at a depth of about 250 km. On the basis of these stress differences a strength of about 100 bar is concluded for the upper 400 km of the lunar interior for the last 3.3 b.y.

Arkani-Hamed, J.↗

A new upper limit to the local population II density.

An upper limit to the local population II density is derived, in terms of the mass-luminosity ratio, on the basis of U, B, and V photometric observations of several thousand stars with V magnitudes between 12 and 18 in a region near the North Galactic Pole. The photographic and photoelectric photometry and the reduction procedures are discussed. Models of the density distribution and luminosity function of the population II stars are used to predict their distribution in color and apparent magnitude. The derived local density of population II is found to be significantly lower than previous estimates. Possible causes for this discrepancy are considered.

Weistrop, D.↗

Effects of angular-momentum-changing collisions on dielectronic satellite spectra

A generalization of the low-density corona-model theory of the intensities of satellite lines of the resonance lines emitted by multiply charged helium-like ions in the X-ray spectra of high-temperature plasmas is developed in order to take account of electron collisional transitions between close autoionizing levels. A formalism for determining the population densities of the autoionizing levels and the satellite line intensities is presented, and the quantum mechanical calculation of autoionization, radiative decay and electron-induced collisional transitions is discussed. Generalized emission rate coefficients corresponding to dielectronic recombination and inner-shell-electron collisional excitation are tabulated as functions of electron density and temperature. The procedure is then used to illustrate the effects of angular-momentum-changing electron collisions on the intensities of the satellites of the 1s2p 1P to 1s2 1S resonance line of the ions Al XII and Si XIII. It is concluded that the collisional transitions have their most significant effect on the intensities of satellite transitions from metastable autoionizing states, which may be useful in determining plasma densities.

Jacobs, V. L.↗