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 181 records · Page 10

Hydrogen Dissociation Controls 1-Hexyne Selective Hydrogenation on Dilute Pd-in-Au Catalysts

Increasing the selectivity of the catalytic hydrogenation of alkynes to alkenes is of major importance for the processing of petrochemicals and the production of fine chemicals. Achieving high selectivity for alkene formation at high conversions, however, remains a long-standing challenge in heterogeneous catalysis. Here, the mechanism and origin of the high selectivity of dilute Pd-in-Au catalysts has been studied by a combination of first-principles calculations, microkinetic simulations, and isotopic exchange hydrogenation experiments using Pd 0.04 Au 0.96 nanoparticles embedded in raspberry colloid-templated silica. The Pd is predominantly in the form of isolated atoms, only surrounded by Au atoms, based on prior studies. The simulations indicate that the rate-limiting process for 1-hexyne hydrogenation on Pd monomers in Au(111) is H 2 dissociation, which has a large free energy barrier of 0.86 eV at 363 K and 0.2 bar of H 2 . The C–H bond formation steps, on the other hand, proceed with lower barriers, which contrasts with previous studies of extended Pd catalysts. The microkinetic simulations identify the sizable H 2 dissociation barrier and the small barrier for the hydrogenation of 1-hexyne as key factors that lead to a high selectivity for the production of 1-hexene from 1-hexyne, even at high conversion. The unconventional H 2 dissociation limiting process in combination with the low coverage of weakly bound hydrocarbon intermediates explains the near-zero order of 1-hexyne found experimentally. Furthermore, the partial hydrogenation of 1-hexyne to form 1-hexene is shown to be an irreversible process from our isotopic exchange hydrogenation experiments and is explained by the strongly exothermic nature of the reaction. Diluting active species, like Pd, in a less active host metal, like Au, hence appears promising as a means of tuning the binding energy of reactants and altering reaction profiles, leading to distinct kinetic behavior for an optimal catalytic activity and selectivity. Furthermore, the combination of microkinetic modeling, density functional theory calculations, and isotopic exchange experiments is thus demonstrated to be an effective approach to modeling important catalytic phenomena.

1-hexyne↗

Measurement and modeling of proton-induced reactions on arsenic from 35 to 200 MeV

72 As is a promising positron emitter for diagnostic imaging that can be employed locally using a 72 Se generator. However, current reaction pathways to 72 Se have insufficient nuclear data for efficient production using regional 100–200 MeV high-intensity proton accelerators. In order to address this deficiency, stacked-target irradiations were performed at LBNL, LANL, and BNL to measure the production of the 72 Se/ 72 As positron emission tomography (PET) generator system via 75 As (p, x) between 35 and 200 MeV. This work provides the most well-characterized excitation function for 75 As ( p , 4 n ) 72 Se starting from threshold. Additional focus was given to report the first measurements of 75 As (p, x) 68 Ge and bolster an already robust production capability for the highly valuable 68 Ge/ 68 Ga PET generator. Thick target yield comparisons with prior established formation routes to both generators are made. In total, high-energy proton-induced cross sections are reported for 55 measured residual products from 75 As, nat Cu , and nat Ti targets, where the latter two materials were present as monitor foils. These results were compared with literature data as well as the default theoretical calculations of the nuclear model codes TALYS, COH, EMPIRE, and ALICE. Reaction modeling at these energies is typically unsatisfactory due to few prior published data and many interacting physics models. Therefore, a detailed assessment of the talys code was performed with simultaneous parameter adjustments applied according to a standardized procedure. Particular attention was paid to the formulation of the two-component exciton model in the transition between the compound and preequilibrium regions, with a linked investigation of level density models for nuclei off of stability and their impact on modeling predictive power. This paper merges experimental work and evaluation techniques for high-energy charged-particle isotope production in an extension to an earlier study of this kind.

43 PARTICLE ACCELERATORS↗

Ion plume/S-band carrier interaction study

A study was performed to determine the effects of a mercury ion thruster plume on an S-band telecommunication carrier. Experiments were carried out on a 30 cm thruster in a JPL test chamber. Results from simple analytical models were compared with the above measurements and major discrepancies were discovered. Modifications to the electron density model provided a qualitative explanation, but further work is necessary for a quantitative answer. The results indicate the effects of the plume, on S and X Band telecommunications will be minor, with the possible exception of critical angle blockage.

Stanton, P.↗

Simulation of cryogenic turbopump annular seals

San Andres employed the NBS software package MIPROPS to account for density's dependence on pressure in the simulation of liquid annular seals. His example on a LH2 seal showed a significant change in the mass coefficient compared to a constant density model. San Andres, Yang, and Childs extended this analysis by including the pressure and temperature dependence of density, specific heat, viscosity, volumetric expansion, and thermal conductivity in a coupled solution of the energy, momentum, and continuity equations. Their example showed very significant changes in stiffness and inertia for a high speed (38,000 rpm), large L/D ratio (0.5) LOX seal, as compared to their constant temperature results. The current research rederived the San Andres-Yang-Childs (SYC) analysis and extended it to include not only the Moody friction model of SYC but also the Hir's friction model. The derivation begins with obtaining the local differential equations of continuity, momentum, and energy conservation in the seal. These equations are averaged across the film thickness to obtain the resulting 'bulk flow' differential equations. Shear stress and convective heat loss through the stator (seal) and rotor are related to the Moody and Hir's friction factor model. The Holman analogy is employed to relate heat conduction in or out of the fluid film's boundary layer to the friction induced shear stress.

Palazzolo, Alan B.↗

Partial redistribution effects in the formation of hydrogen lines in quiescent prominences

Departures from complete frequency redistribution (CRD) in hydrogen lines are investigated for solar prominences. Partial redistribution effects (PRD) are found both in the wings (their already known lowering) and in the central part of the L alpha line; a new feature is evidenced here: the partially coherent scattering in the near wings of the line leads to a double-peaked profile mirroring the incident solar radiation. With a low density model, we obtain a good agreement with OSO 8 observed profiles. On the contrary, the PRD computed L beta profile (lower density, no reversal) departs from the observed one, a result which calls for more progress in terms of non-LTE transfer and modelling.

Heinzel, P.↗

Analytical model for the density distribution in the Io plasma torus

An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.

Mei, YI↗

Semianalytical Propagation of Satellite Orbits about an Arbitrary Central Body

Precision mean element (PME) satellite theories play a key role in orbit dynamics analyses. These theories employ: nonsingular orbital elements comprehensive force models Generalized Method of Averaging Numerical interpolation concepts The Draper Semianalytical Satellite Theory (DSST) (Refs. 1 - 6), whose development was led by the author, and the independently-developed Universal Semianalytical Method (USM) (Ref. 7) are examples of such theories. These theories provide the capability to tailor the force modeling to meet the desired computational speed vs. accuracy trade-off. The flexibility of such theories is demonstrated by their ability to include complicated atmosphere density models and spacecraft models in the perturbation theory context. The value of high speed satellite theories, in this era of computational plenty, is that they allow new ways of looking at astrodynamical problems such as orbit design (Refs. 8, 9) and atmosphere density updating (Refs. 10, 11). In the mid to late-1980 s, the geodynamics community led the development of very precise geopotential models such as GEM T2 and GEM T3 (Ref. 12), and with the subsequent analysis of the TOPEX flight data, JGM-2 and JGM-3 (Ref. 13). These were high degree and order geopotentials, at least 50 x 50. In 1993, the DSST implementation in the GTDS program was extended to include the 50 x 50 geopotential models (Ref. 14). The 50 x 50 geopotential, J2000 integration coordinate system, and solid Earth tide capabilities were integrated in GTDS by Scott Carter (Ref. 15). This capability demonstrated 1 m accuracy versus the TOPEX Precise Orbit Ephemerides. Subsequently the DSST Standalone program was also extended to include high degree and order geopotential models (Ref. 5). More recently GTDS has been hosted in the Linux PC environment. However, all of these efforts have been limited to modeling the motion of an artificial Earth satellite. They did not consider the additional complexities associated with lunar, planetary, or other natural satellite orbiters. Such complexities include: additional coordinate systems (associated with the direction of the north pole of rotation and the prime meridian of the new central bodies) (Ref. 16) normalized gravity model coefficients (desirable for high degree and order fields) (Ref. 17) indirect oblateness

Cefola, Paul J.↗

A review of orbital debris environment modeling at NASA/JSC

Orbital debris modeling in three areas is presented in this paper. These are: (1) debris environment evolution, as characterized by a numerical processor using detailed mission model projections, (2) debris cloud evolution for the time immediately after breakup, and (3) two-dimensional flux density modeling to characterize both flux levels and directionality for a debris environment.

Reynolds, Robert C.↗

Exospheric temperatures during solar cycle 22

An analysis is conducted of the orbital decay of the Long Duration Exposure Facility (LDEF) satellite in order to develop an empirical model of the solar radio flux at 10.7 cm and of the exospheric temperature driving current atmospheric-density models at high radio fluxes. The Jacchia-Lineberry (1982) atmospheric model is used to fit a set of orbit semimajor axis observations to a predicted set whose only free variable was the exospheric temperature. A comparison of the exospheric temperatures with those predicted by the Jacchia (1977) model indicates increasing divergences above 1100 K, suggesting needed modifications to the model.

Badhwar, Gautam D.↗

Subsurface Seismic Structure Around the Sanford Underground Research Facility

To address the fundamental challenges of understanding processes associated with enhanced or engineered geothermal systems (EGS), teams from the EGS Collab SIGMA-V project are conducting stimulations at the Sanford Underground Research Facility. Various modalities of geophysical observations are collected and analyzed, and different types of parameters are modeled to help investigate these critical processes. Among these parameters, the three-dimensional stress field is one of the primary objectives. Seismic velocity models have been used to estimate the 3D stress field for other regions. A detailed seismic velocities and density model can help us infer spatial variations of the stress field in the study area. We inverted the 3D seismic structure around the site by using body-wave travel times, surface-wave dispersions and gravity observations simultaneously. The resulting model will be a starting point for more detailed investigations and could also be used to study subsurface geological and stress field variations of the broad area.

Chai, Chengping↗

Subsurface Seismic Structure Around the Sanford Underground Research Facility

To address the fundamental challenges of understanding processes associated with enhanced or engineered geothermal systems (EGS), teams from the EGS Collab SIGMA-V project are conducting stimulations at the Sanford Underground Research Facility. Various modalities of geophysical observations are collected and analyzed, and different types of parameters are modeled to help investigate these critical processes. Among these parameters, the three-dimensional stress field is one of the primary objectives. Seismic velocity models have been used to estimate the 3D stress field for other regions. A detailed seismic velocities and density model can help us infer spatial variations of the stress field in the study area. We inverted the 3D seismic structure around the site by using body-wave travel times, surface-wave dispersions and gravity observations simultaneously. The resulting model will be a starting point for more detailed investigations and could also be used to study subsurface geological and stress field variations of the broad area.

Chai, Chengping↗

Post-glacial relaxation of a viscously stratified compressible mantle

The postglacial relaxation phenomenon of a viscously compressible nonself-gravitating spherical-shell model is investigated. Analytical solutions to an exponentially depth-dependent viscosity are developed with various types of compressible density models such as those with an exponential dependence of the radius and an algebraic-root dependence of the radius. The solutions of a viscously compressible multilayered model with constant thermodynamic properties and viscosity are developed by the propagator matrix method. The results show that inferences of deep mantle viscosity from postglacial rebound would be hampered by mantle compressibility for long-wavelength harmonics because of the smaller excitation of compressible eigenfunctions in the lower mantle. The relaxation times and velocity fields are more sensitive to higher viscosity contrast for the exponentially varying viscosity than for models with discrete jumps in the viscosity structure.

Wu, J.↗

Cellular Automaton Modeling of Dynamic Recrystallization of Nimonic 80A Superalloy Based on Inhomogeneous Distribution of Dislocations Inside Grains

A modified cellular automaton (CA) model in the time and space scale is developed to simulate the dynamic recrystallization (DRX) mechanism of Nimonic 80A nickel-base superalloy during hot deformation. Some necessary thermo-mechanical parameters used as input data in CA approach such as work hardening and softening coefficient are fitted according to the previous hot compression experiments. Based on the experimental results, nucleation kinetics is modified as a function of Zener-Holloman parameter. A dislocation density model at transient recrystallization boundary is proposed, in which the dislocation density is a function of temperature, strain rate and annealed dislocation density. The simulated results are in good agreement with experimental results, proving that the modified CA model is capable of predicting both flow behavior and microstructural evolution of Nimonic 80A during hot deformation. According to the experimental and simulated results, the DRX behavior is discussed in detail. It demonstrates that the flow stress and microstructural evolution of DRX have a relationship with the deformation temperature and strain rate.

36 MATERIALS SCIENCE↗

Collapse of the electron gas from three to two dimensions in Kohn-Sham density functional theory

Here under pressure, a quasi-2D electron gas can collapse toward the true 2D limit. In this limit, the exact exchange-correlation energy per electron has a known finite limit, but general-purpose semilocal approximate density functionals, like the local density approximation (LDA) and Perdew-Burke-Ernzerhof generalized gradient approximation (PBE GGA), are known to diverge to minus infinity. Here we consider a model density for a non-interacting electron gas confined to a thickness $L$ by infinite-barrier walls, with a fixed 2D density $1/(\pi (r_s^{\text{2D}})^2)$, and $r_s^{\text{2D}}=4$ Bohr. We estimate that LDA, PBE, and the strongly constrained and appropriately normed (SCAN) meta-GGA are accurate for the exchange-correlation energy over a wide quasi-2D range, $1.5< L/r_s^{\text{2D}} < 3.85$, but not for smaller $L$. Of these functionals, only SCAN tends to a finite limit when $L$ tends to 0. Since the non-interacting kinetic energy, treated exactly in Kohn-Sham theory, dominates in this limit within a deformable jellium model, all of the general-purpose functionals can estimate the pressure required to achieve any thickness (with SCAN and LDA better than PBE). This pressure vanishes around $L/r_s^{\text{2D}} = 3.85$, where the 3D electron density is roughly that of the valence electrons in metallic potassium, and reaches about 20 GPa at $L/r_s^{\text{2D}} = 1.5$ and 400 GPa at $L/r_s^{\text{2D}} = 0.6$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The role of particle diffusion in the lower transition region: Revised interpretation of emission measures

Our energy-balance models of the lower transit ion region were presented in a previous paper (New Models of the Chromosphere and Transition Region). Here we show the influence of particle diffusion on the calculated hydrogen and helium number densities for a given temperature-density model (model C in the preceding paper). We have also solved the statistical equilibrium and radiative transfer equations for a 13-level He I atom (22 radiative transitions) and a 6-level He II ion (15 radiative transitions) together with He III. The resulting He I and He II level-1 number densities and He III density are shown as functions of temperature. Diffusion substantially increases eta(sub HeI) for T greater than 35,000 K and decreases eta(sub HeI) (while increasing eta(sub(HeII)) between 9,000 and 25,000 K. Including the effects of diffusion also increases eta(sub HeIII) for T less than 60,000 K.

Avrett, Eugene H.↗

Analysis of extreme-ultraviolet observations of a polar coronal hole

Emission gradient curves for extreme-ultraviolet resonance lines of lithiumlike ions have been constructed from spectroheliograms of a northern polar coronal hole observed on August 14, 1973, with the Harvard experiment on Skylab. An emission-measure analysis indicates both reduced density and coronal temperature in the coronal hole. The boundary geometry of the coronal hole is determined, and a temperature-density model that is consistent with the observed intensities is constructed. The model gives a conductive flux of 60,000 erg/sq cm per sec at 1.03 solar radii. The boundary geometry and density distribution are combined with typical solar-wind parameters at the earth to determine an outflow velocity of 15 km/s at 1.08 solar radii.

Mariska, J. T.↗

An information theory approach to the density of the earth

Information theory can develop a technique which takes experimentally determined numbers and produces a uniquely specified best density model satisfying those numbers. A model was generated using five numerical parameters: the mass of the earth, its moment of inertia, three zero-node torsional normal modes (L = 2, 8, 26). In order to determine the stability of the solution, six additional densities were generated, in each of which the period of one of the three normal modes was increased or decreased by one standard deviation. The superposition of the seven models is shown. It indicates that current knowledge of the torsional modes is sufficient to specify the density in the upper mantle but that the lower mantle and core will require smaller standard deviations before they can be accurately specified.

Graber, M. A.↗

A Coupled Plasma and Sheath Model for High Density Reactors

We present a coupled plasma and collisionless; sheath model for the simulation of high density plasma processing reactors. Due to inefficiencies in numerical schemes and the resulting computational burden, a coupled multidimensional plasma and sheath simulation has not been possible model for gas mixtures and high density reactors of practical interest. In this work we demonstrate that with a fully implicit algorithm and a refined computational mesh, a self-consistent plasma and sheath simulation is feasible. We discuss the details of the model equations, the importance of ion inertia, and the resulting sheath profiles for argon and chlorine plasmas. We find that at low operating pressures (10-30 mTorr), the charge separation occurs only within a 0.5 mm layer near the surface in a 300 mm inductively coupled plasma etch reactor. A unified model eliminates the use of off-line or loosely coupled sheath models with simplifying assumptions which generally lead to uncertainties in ion flux and sheath electrical properties.

Deepak, Bose↗