SEARCH · Engineering Papers

Results for “Probability density function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Empirical radius formulas for canonical neutron stars from bidirectionally selecting features of equations of state in extended Bayesian analyses of observational data

Significant advancement in Bayesian inference of nuclear equation of state (EOS) from gravitational wave and x-ray observations of neutron stars (NSs) has been made by the nuclear astrophysics community especially since GW170817. By extending the traditional Bayesian analysis which normally ends at presenting the marginalized posterior probability distribution functions (PDFs) of individual EOS parameters and their correlations (or sometimes only the Pearson correlation coefficients which are only reliably useful when the variables are linearly correlated while they are actually often not), we search for a data-driven and robust empirical formula for the radius 𝑅 1.4 of canonical NSs in terms of the characteristic EOS parameters (features). We also identify the single most important but currently poorly known EOS parameter for determining the 𝑅 1.4 . Using three regression-model-building methodologies: bidirectional stepwise feature selection, least absolute shrinkage selection operator (LASSO) regression, and neural network regression on a large set of posterior EOSs and the corresponding 𝑅 1.4 values inferred from earlier comprehensive Bayesian analyses of NS observational data, we systematically and rigorously develop the most probable 𝑅 1.4 formulas with varying statistical accuracy and technical complexity. Here, the most important EOS parameters for determining 𝑅 1.4 are found consistently in each of the feature selection processes to be (in order of decreasing importance): curvature 𝐾 sym , slope 𝐿, skewness 𝐽 sym of nuclear symmetry energy, skewness 𝐽 0 , incompressibility 𝐾 0 of symmetric nuclear matter, and the magnitude 𝐸 sym ⁡(𝜌 0 ) of symmetry energy at the saturation density 𝜌 0 of nuclear matter.

Bayesian methods↗

NASA NTRS2013

NASA Tech Briefs, May 2013

Topics include: Test Waveform Applications for JPL STRS Operating Environment; Pneumatic Proboscis Heat-Flow Probe; Method to Measure Total Noise Temperature of a Wireless Receiver During Operation; Cursor Control Device Test Battery; Functional Near-Infrared Spectroscopy Signals Measure Neuronal Activity in the Cortex; ESD Test Apparatus for Soldering Irons; FPGA-Based X-Ray Detection and Measurement for an X-Ray Polarimeter; Sequential Probability Ratio Test for Spacecraft Collision Avoidance Maneuver Decisions; Silicon/Carbon Nanotube Photocathode for Splitting Water; Advanced Materials and Fabrication Techniques for the Orion Attitude Control Motor; Flight Hardware Packaging Design for Stringent EMC Radiated Emission Requirements; RF Reference Switch for Spaceflight Radiometer Calibration; An Offload NIC for NASA, NLR, and Grid Computing; Multi-Scale CNT-Based Reinforcing Polymer Matrix Composites for Lightweight Structures; Ceramic Adhesive and Methods for On-Orbit Repair of Re-Entry Vehicles; Self-Healing Nanocomposites for Reusable Composite Cryotanks; Pt-Ni and Pt-Co Catalyst Synthesis Route for Fuel Cell Applications; Aerogel-Based Multilayer Insulation with Micrometeoroid Protection; Manufacturing of Nanocomposite Carbon Fibers and Composite Cylinders; Optimized Radiator Geometries for Hot Lunar Thermal Environments; A Mission Concept: Re-Entry Hopper-Aero-Space-Craft System on-Mars (REARM-Mars); New Class of Flow Batteries for Terrestrial and Aerospace Energy Storage Applications; Reliability of CCGA 1152 and CCGA 1272 Interconnect Packages for Extreme Thermal Environments; Using a Blender to Assess the Microbial Density of Encapsulated Organisms; Mixed Integer Programming and Heuristic Scheduling for Space Communication; Video Altimeter and Obstruction Detector for an Aircraft; Control Software for Piezo Stepping Actuators; Galactic Cosmic Ray Event-Based Risk Model (GERM) Code; Sasquatch Footprint Tool; and Multi-User Space Link Extension (SLE) System.

Source record↗

NASA NTRS1994

Photoionization of the diffuse interstellar medium and galactic halo by OB associtations

Assuming smoothly varying H I distributions in te Galactic disk, we have calculated the geometry of diffuse II regions due to OB associations in the Galactic plane. Near the solar circle, OB associations with a Lyman continuum (Lyc) photon luminosity Psi(sub Lyc) = 3.3 x 10(exp 7) cm(exp -2) s(exp -1), produce H II regions that are density bounded in the vertical direction (H II chimneys) allowing Lyc to escape the gaseous disk and penetrate into the Galactic halo. We provide analytic formulae for the Lyc escape fraction as functions of S(sub 0) O-star catalog of Garmany and a new Lyc stellar stellar Lyc stellar flux calibration, we find a production rate of Lyc photons by OB associations within 2.5 kpc of Psi(sub Lyc) = 3.3 x 10(exp 7) cm(exp -2) s(exp -1). Integrating the fraction of Lyc photons that escape the disk over our adopted luminosity function of OB associations, we estimate that approximately 7% of the ionizing photons, or Phi(sub Lyc) = 2.3 x 10(exp 6) cm(exp -2) s(exp -1), escape each side of the H I disk layer and penetrate the diffuse ionized medium ('Reynolds layer'). This flux is sufficient to explain the potoionization of this, although we have not constructed a model for the observed H-alpha emission and pulsar dispersion measures that is fully consistent with the absorption rate of Lyc in the H II layer. Since our quiescent model does not account for the effects of dynamic chimneys and superbubbles, which should enhance Lyc escape, we conclude the O stars are the probable source of ionizing radiation for the Reynolds layer. For a random distribution of OB associations throughout the disk, the Lyc flux is nearly uniform for heights Z is greater than approximately 0.8 kpc above the midplane.

Dove, James B.↗

NASA NTRS1998

HST Spatially Resolved Spectra of the Accretion Disc and Gas Stream of the Nova-Like Variable UX Ursae Majoris

Time-resolved eclipse spectroscopy of the nova-like variable UX UMa obtained with the Hubble Space Telescope/Faint Object Spectrograph (HST/FOS) on 1994 August and November is analysed with eclipse mapping techniques to produce spatially resolved spectra of its accretion disk and gas stream as a function of distance from the disk centre. The inner accretion disk is characterized by a blue continuum filled with absorption bands and lines, which cross over to emission with increasing disk radius, similar to that reported at optical wavelengths. The comparison of spatially resolved spectra at different azimuths reveals a significant asymmetry in the disk emission at ultraviolet (UV) wavelengths, with the disk side closest to the secondary star showing pronounced absorption by an 'iron curtain' and a Balmer jump in absorption. These results suggest the existence of an absorbing ring of cold gas whose density and/or vertical scale increase with disk radius. The spectrum of the infalling gas stream is noticeably different from the disc spectrum at the same radius suggesting that gas overflows through the impact point at the disk rim and continues along the stream trajectory, producing distinct emission down to 0.1 R(sub LI). The spectrum of the uneclipsed light shows prominent emission lines of Lyalpha, N v lambda1241, SiIV Lambda 1400, C IV Lambda 1550, HeII Lambda 1640, and MgII Lambda 2800, and a UV continuum rising towards longer wavelengths. The Balmer jump appears clearly in emission indicating that the uneclipsed light has an important contribution from optically thin gas. The lines and optically thin continuum emission are most probably emitted in a vertically extended disk chromosphere + wind. The radial temperature profiles of the continuum maps are well described by a steady-state disc model in the inner and intermediate disk regions (R greater than or equal to 0.3R(sub LI) ). There is evidence of an increase in the mass accretion rate from August to November (from V = 10 (exp -8.3 +/-0.1) to 10(exp -8.1 +/- 0.1 solar mass yr(exp -1)), in accordance with the observed increase in brightness. Since the UX UMA disc seems to be in a high mass accretion, high-viscosity regime in both epochs, this result suggests that the mass transfer rate of UX UMA varies substantially (approximately equal to 50 per cent) on time-scales of a few months. It is suggested that the reason for the discrepancies between the prediction of the standard disk model and observations is not an inadequate treatment of radiative transfer in the disc atmosphere, but rather the presence of addition important sources of light in the system besides the accretion disk (e.g., optically thin contiuum emission from the disk wind and possible absorption by circumstellar cool gas).

Baptista, Raymundo↗

DOE OSTI2025

Illuminating the mechanistic impacts of an Fe-quaterpyridine functionalized crystalline poly(triazine imide) semiconductor for photocatalytic CO 2 reduction

The strategy of incorporating earth-abundant catalytic centers into light-absorbing architectures is desirable from the viewpoint of low cost, low toxicity, and versatility for activating small molecules to produce solar-based fuels. Herein, we show that an Fe-quaterpyridine molecular catalyst can be anchored to a light-absorbing, crystalline, carbon nitride (PTI), to yield a molecular-catalyst/material hybrid, Fe-qpy-PTI, capable of facilitating CO 2 reduction to CO selectively (up to ∼97–98%) in aqueous solution under low-intensity light irradiation. This hybrid material leverages the ability of the Fe-qpy catalyst to bind CO 2 upon a one-electron reduction, as achieved by transfer of excited electrons from the carbon-nitride semiconductor. At a low incident power density of only 50 mW cm −2 , the catalytic activity of the hybrid material was measured across a range of catalyst loadings from 0.1–3.8 wt%, yielding CO rates of up to 596 μmol g −1 h −1 for a 3.8 wt% loading during a 3 h experiment. Over the course of 8 h, the hybrid material attained a CO evolution rate of 608 μmol g −1 h −1 and 305 turnovers for a TOF of ∼38 h −1 and an apparent quantum yield of 2.6%. Higher light intensities provided an initial increase in activity but negatively impacted photocatalytic rates with time, with an AQY of 0.6% at 150 mW cm −2 and 0.4% at 250 mW cm −2 . Transient absorption spectroscopy results showed electron survival probabilities consistent with the trends in observed product rates. Computational modeling was also used to evaluate and understand the mechanistic pathway of the high product selectivity for CO versus H 2 . These results thus help unveil key factors for leveraging the mechanistic understanding of molecular catalysts for CO 2 reduction for pairing with light absorbing semiconductors and establishing optimal conditions to attain maximal rates in aqueous solution.

McGuigan, Scott [Baylor University, Waco, TX (Unit↗

DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗

NASA NTRS2004

Risk Assessment and Control through Countermeasure System Iplementation for Long-term Crew Exposure to Microgravity

Experience with the International Space Station (ISS) program demonstrates the degree to which engineering design and operational solutions must protect crewmembers from health risks due to long-term exposure to the microgravity environment. Risks to safety and health due to degradation in the microgravity environment include crew inability to complete emergency or nominal activities, increased risk of injury, and inability to complete safe return to the ground due to reduced strength or embrittled bones. These risks without controls slowly increase in probability for the length of the mission and become more significant for increasing mission durations. Countermeasures to microgravity include hardware systems that place a crewmember s body under elevated stress to produce an effect similar to daily exposure to gravity. The ISS countermeasure system is predominately composed of customized exercise machines. Historical treatment of microgravity countermeasure systems as medical research experiments unintentionally reduced the foreseen importance and therefore the capability of the systems to function in a long-term operational role. Long-term hazardous effects and steadily increasing operational risks due to non-functional countermeasure equipment require a more rigorous design approach and incorporation of redundancy into seemingly non- mission-critical hardware systems. Variations in the rate of health degradation and responsiveness to countermeasures among the crew population drastically increase the challenge for design requirements development and verification of the appropriate risk control strategy. The long-term nature of the hazards and severe limits on logistical re-supply mass, volume and frequency complicates assessment of hardware availability and verification of an adequate maintenance and sparing plan. Design achievement of medically defined performance requirements by microgravity countermeasure systems and incorporation of adequate failure tolerance significantly reduces these risks. Future implementation of on-site monitoring hardware for critical health parameters such as bone mineral density would allow greater responsiveness, efficiency, and optimized design of the countermeasures system.

Gernand, Jeremy M.↗

DOE OSTI2024

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1997

Improvement in Suppression of Pulsed Nd:YAG Laser Light With Iodine Absorption Cells for Filtered Rayleigh Scattering Measurements

Filtered Rayleigh scattering using iodine absorption cells is an effective technique for obtaining density, temperature, and velocity measurements in high speed confined flows. By tuning a single frequency laser to a strong iodine absorption line, stray scattered laser light can be greatly suppressed. For example, the minimum transmission predicted by an iodine absorption model calculation is less than 10(exp -5) at the 18788.44/cm line using a 200 mm absorption cell containing iodine vapor at 0.46 T. Measurements obtained by other researches using a CW Nd:YAG laser agree with the model calculations. However, measurements made by us and by others using Q-switched, injection-seeded, frequency doubled Nd:YAG lasers only show minimum transmission of about 3 x 10(exp -3). This greatly reduces the applicability of the filtered Rayleigh scattering technique using these lasers in experiments having large amounts of stray scattered laser light. The purposes of the present study are to characterize the spectrum of the excess light transmitted by the iodine cell and to make changes to the laser to reduce the transmitted laser light. Transmission data as a function of laser frequency for the iodine absorption line at 18788.44/cm are presented. A planar mirror Fabry-Perot interferometer was used to characterize the frequency spectrum of the light passed through the cell. Measurements taken with the laser tuned to the center of the iodine absorption line show the light transmitted through the iodine cell to have a component with a bandwidth of about 40 GHz. This is probably caused by other modes in the laser that exist in spite of the single frequency injection beam. A second broadband component was also observed, possibly caused by the laser flash lamps or by fluorescence. An intracavity etalon was installed in the laser oscillator cavity to suppress the 40 GHz component. Measurements taken with the etalon tuned to the injection frequency showed a reduction in the transmitted laser light. This improvement allows the iodine cell to block significantly more of the stray laser light in filtered Rayleigh scattering experiments. Examples are given of filtered Rayleigh scattering measurements showing the effect of the etalon on measurements taken in a Mach 3 flow in the NASA Lewis 4 inch by 10 inch supersonic wind tunnel.

Seasholtz, Richard G.↗

NASA NTRS2007

Titan Orbiter Aerorover Mission with Enceladus Science (TOAMES)

Cassini and Huygens have made exciting discoveries at Titan and Enceladus, and at the same time made us aware of how little we understand about these bodies. For example, the source, and/or recycling mechanism, of methane in Titan's atmosphere is still puzzling. Indeed, river beds (mostly dry) and lakes have been spotted, and occasional clouds have been seen, but the physics to explain the observations is still mostly lacking, since our "image" of Titan is still sketchy and quite incomplete. Enceladus, only -500 km in extent, is even more puzzling, with its fiery plumes of vapor, dust and ice emanating from its south polar region, "feeding" Saturn's E ring. Long term variability of magnetospheric plasma, neutral gas, E-ring ice grain density, radio emissions, and corotation of Saturn's planetary magnetic field in response to Enceladus plume activity are of great interest for Saturn system science. Both Titan and Enceladus are bodies of considerable astrobiological interest in view of high organic abundances at Titan and potential subsurface liquid water at Enceladus. We propose to develop a new mission to Titan and Enceladus, the Titan Orbiter Aerorover Mission with Enceladus Science (TOAMES), to address these questions using novel new technologies. TOAMES is a multi-faceted mission that starts with orbit insertion around Saturn using aerobraking with Titan's extended atmosphere. We then have an orbital tour around Saturn (for 1-2 years) and close encounters with Enceladus, before it goes into orbit around Titan (via aerocapture). During the early reconnaissance phase around Titan, perhaps 6 months long, the orbiter will use altimetry, radio science and remote sensing instruments to measure Titan's global topography, subsurface structure and atmospheric winds. This information will be used to determine where and when to release the Aerorover, so that it can navigate safely around Titan and identify prime sites for surface sampling and analysis. In situ instruments will sample the upper atmosphere which may provide the seed population for the complex organic chemistry on the surface. The Aerorover will probably use a "hot air" Montgolfier balloon concept using the waste heat from the MMRTG 1-2 kwatts. New technologies will need to be developed and miniaturization will be required to maintain functionality while controlling mass, power and cost. Duty cycling will be used. The Aerorover will have all the instruments needed to sample Titan's atmosphere and surface with possible methane lakes-rivers. It will e.g., use multi-spectral imagers and for last 6 months of mission, balloon payload will land on surface at predetermined site to take core samples of the surface and use seismometers to help probe the interior. All remote (and active) sensors on the orbiter will share a - 1 meter telescope, called MIDAS (Multiple Instrument Distributed Aperture Sensor). MIDAS observations in stable orbit at Titan can provide full global maps of Titan's surface and could additionally provide long term observations of the Saturn system including Enceladus for extended mission phases over many years, potentially for decades. Experience from the Hubble Space Telescope has shown strong public interest and commitment to exciting generational missions.

Sittler, Edward C.↗

DOE OSTI2025

Medium-range order and compositional correlation in metallic glasses

The compositional atomic ordering in metallic glasses was studied by simulation focusing on the medium-range order (MRO). Many metallic alloy liquids and glasses show MRO characterized by the oscillations in the atomic pair-distribution function (PDF) beyond the first peak, which decay exponentially with distance. To study the effects of the local chemical order on MRO, we examine the compositionally resolved PDF and its MRO for models of various binary metallic alloy glasses. We show that compositional ordering is limited mostly to the nearest-neighbor atoms and the MRO is largely independent of the compositional order. For some elements that strongly repel each other in the alloy, a second MRO periodicity is observed owing to the distinct correlations among them. These results are discussed in light of the idea that the MRO oscillations in the PDF describe the correlations in the atomic density fluctuations, rather than the detailed local atomic structure.

Atomic structure↗

DOE OSTI2020

Steric effects in light-induced solvent proton abstraction

The significance of solvent structural factors in the excited-state proton transfer (ESPT) reactions of Schiff bases with alcohols is reported here. In this work, we use the super photobase FR0-SB and a series of primary, secondary, and tertiary alcohol solvents to illustrate the steric issues associated with solvent to photobase proton transfer. Steady-state and time-resolved fluorescence data show that ESPT occurs readily for primary alcohols, with a probability proportional to the relative –OH concentration. For secondary alcohols, ESPT is greatly diminished, consistent with the barrier heights obtained using quantum chemistry calculations. ESPT is not observed in the tertiary alcohol. We explain ESPT using a model involving an intermediate hydrogen-bonded complex where the proton is “shared” by the Schiff base and the alcohol. The formation of this complex depends on the ability of the alcohol solvent to achieve spatial proximity to and alignment with the FR0-SB* imine lone pair stabilized by the solvent environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1961

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗

DOE OSTI2020

Etching of Si 3 N 4 induced by electron beam plasma from hollow cathode plasma in a downstream reactive environment

An etching system based on the interaction of electrons extracted from a direct current (DC) hollow-cathode (HC) Ar plasma and injected towards a Si3N4 covered silicon substrate located in the downstream reactive environment created by an additional remote CF 4 /O 2 plasma source was developed and evaluated. By controlling the properties of the injected beam electrons, this approach allows to deliver energy to a surface functionalized by exposure to reactive species, and initiate surface etching. The energy of the primary beam electrons is controlled by the acceleration voltage relative to the HC discharge. Ar atoms flow from the high-pressure HC discharge into the low pressure downstream reactive environment in the process chamber. For an acceleration voltage greater than the ionization potential of Ar and/or process gas species, the energetic primary beam electrons produce a secondary plasma in the process chamber, and can also cause additional dissociation. We have characterized the properties of the secondary plasma and also surface etching of Si 3 N 4 as a function of process parameters, including acceleration voltage (0-80 V), discharge current of the HC discharge (1-2 A), pressure (3.5-20 mTorr), source to substrate distance (1.5 to 5 cm), and feed gas composition (20% and 80% O2 in CF 4 /O 2 ). The electron energy probability function (EEPF) measured with a Langmuir probe about 2.5 cm below the extraction ring suggests several major groups of electrons for this situation, including high energy primary beam electrons with an energy that varies as the acceleration voltage is changed and low energy electrons produced by beam electron-induced ionization of the Ar gas in the process chamber. When a remote CF 4 /O 2 plasma is additionally coupled to the process chamber, Si 3 N 4 surfaces can be functionalized, and by varying the energy of the beam electrons, Si 3 N 4 etching can be induced by electron-neutral synergy effect with plasma-surface interaction. For conditions without beam electron injection, the remote plasma etching rate (ER) of Si 3 N 4 depends strongly on the O 2 concentration in the CF 4 /O 2 processing gas mixture and can be suppressed for O 2 -rich process conditions by the formation of a SiONF passivation layer on the Si 3 N 4 surface. The combination of the HC electron beam source with the remote plasma source (HCEB-RP) makes it possible to induce Si 3 N 4 etching for O 2 -rich remote plasma conditions where remote plasma by itself produces negligible Si 3 N 4 etching. The electron enhanced etching of Si 3 N 4 depends strongly on O 2 /CF 4 mixing ratio reflecting changing arrival rates of O and F species at the surface. Optical emission spectroscopy was used to estimate the ratio of gas phase F and O densities and found to be controlled by the gas mixing ratio and independent of HC EB operating conditions. At this time the detailed sequence of events operative in the etching mechanism is unclear. While the increase of the electron energy is ultimately responsible for initiating surface etching, at present we cannot rule out a role of ions from the simultaneously produced secondary plasma in plasma surface interaction mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2024

Quantification of Uncertainty and Risk Sensitivity for Safety of Emerging Operations

The growing need to develop and deploy small unmanned aerial vehicles (sUAVs) for various applications in the airspace necessitates reliable tools to accurately predict the flight trajectories of the sUAVs. The knowledge of the predicted trajectories help decision makers anticipate potential conflict, assess the risk, and take appropriate risk mitigation actions. In addition, uncertainties in vehicle models, weather, and controller action further highlights the need for reliable prediction tools. In this project, the application of mixed sparse grid-based quadrature and generalized polynomial chaos(gPC) expansion method for uncertainty quantification and collision assessment in air traffic consisting of fixed-wing small unmanned aerial vehicles (sUAV) was studied. From the results obtained, it can be concluded that this provides a reliable framework to carry out quantitative conflict assessment in an unmanned air traffic, which when employed, can improve the functionalities of the unmanned traffic management system. It was observed that the results from the gPC expansion framework developed in the project can be utilized to conduct rapid probabilistic collision assessment for near real-time unmanned traffic management in the airspace. From the vehicle models, position updates, and wind-field data, a priori gPC based 3-σcon-fidence ellipses can provide estimates of potential conflict at some future instants. The computational costs scaled linearly when the uncertain inputs were fewer. Further, the largest allowable distribution of para-metric uncertainties that leads to the smallest risk of collision in traffic of small unmanned aerial vehicles could be calculated. The time of closest approach between two sUAVs can be established paving way for development of proactive mitigation strategies. The separation between the sUAVs was found to be most significantly affected by uncertainties in the maximum available thrusts, zero-lift drag coefficients, and wing planform areas of the sUAVs. The study of uncertain wind-fields indicated that a heterogeneous traffic mix resulted in an increased probability of conflict. Increased measurement update rate reduced the uncertain-ties in the trajectories of the vehicles, further reducing the probability of conflict but rapid updates of all vehicles in the airspace poses a stringent communication limitation. The gPC framework also provided the means to analyze vehicle impact (crash region) due to loss of control resulting from actuator failure in sUAS traffic, essentially to predict impact and crash zones for representative vehicles. The predicted regions when compared with non-participant density, provides a means to develop an early mitigation strategy, should the sUAV detect an imminent actuator failure.

Rajnish Bhusal↗

NASA NTRS2008

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Aerosols and especially their effect on clouds are one of the key components of the climate system and the hydrological cycle [Ramanathan et al., 2001]. Yet, the aerosol effect on clouds remains largely unknown and the processes involved not well understood. A recent report published by the National Academy of Science states "The greatest uncertainty about the aerosol climate forcing - indeed, the largest of all the uncertainties about global climate forcing - is probably the indirect effect of aerosols on clouds [NRC, 2001]." The aerosol effect on clouds is often categorized into the traditional "first indirect (i.e., Twomey)" effect on the cloud droplet sizes for a constant liquid water path [Twomey, 1977] and the "semi-direct" effect on cloud coverage [e.g., Ackerman et al ., 2001]." Enhanced aerosol concentrations can also suppress warm rain processes by producing a narrow droplet spectrum that inhibits collision and coalescence processes [e.g., Squires and Twomey, 1961; Warner and Twomey, 1967; Warner, 1968; Rosenfeld, 19991. The aerosol effect on precipitation processes, also known as the second type of aerosol indirect effect [Albrecht, 1989], is even more complex, especially for mixed-phase convective clouds. Table 1 summarizes the key observational studies identifying the microphysical properties, cloud characteristics, thermodynamics and dynamics associated with cloud systems from high-aerosol continental environments. For example, atmospheric aerosol concentrations can influence cloud droplet size distributions, warm-rain process, cold-rain process, cloud-top height, the depth of the mixed phase region, and occurrence of lightning. In addition, high aerosol concentrations in urban environments could affect precipitation variability by providing an enhanced source of cloud condensation nuclei (CCN). Hypotheses have been developed to explain the effect of urban regions on convection and precipitation [van den Heever and Cotton, 2007 and Shepherd, 2005]. Please see Tao et al. (2007) for more detailed description on aerosol impact on precipitation. Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (GCE) model. Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and summertime convection over a mid-latitude continent with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. The impact of atmospheric aerosol concentration on cloud and precipitation will be investigated.

Tao, Wei-Kuo↗

NASA NTRS2013

Enigma Version 12

Enigma Version 12 software combines model building, animation, and engineering visualization into one concise software package. Enigma employs a versatile user interface to allow average users access to even the most complex pieces of the application. Using Enigma eliminates the need to buy and learn several software packages to create an engineering visualization. Models can be created and/or modified within Enigma down to the polygon level. Textures and materials can be applied for additional realism. Within Enigma, these models can be combined to create systems of models that have a hierarchical relationship to one another, such as a robotic arm. Then these systems can be animated within the program or controlled by an external application programming interface (API). In addition, Enigma provides the ability to use plug-ins. Plugins allow the user to create custom code for a specific application and access the Enigma model and system data, but still use the Enigma drawing functionality. CAD files can be imported into Enigma and combined to create systems of computer graphics models that can be manipulated with constraints. An API is available so that an engineer can write a simulation and drive the computer graphics models with no knowledge of computer graphics. An animation editor allows an engineer to set up sequences of animations generated by simulations or by conceptual trajectories in order to record these to highquality media for presentation. Enigma Version 12 Lyndon B. Johnson Space Center, Houston, Texas 28 NASA Tech Briefs, September 2013 Planetary Protection Bioburden Analysis Program NASA's Jet Propulsion Laboratory, Pasadena, California This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware.

Shores, David↗

DOE OSTI2020

Electric and magnetic dipole strength in 112,114,116,118,120,124 Sn

Background: There is renewed interest in electric dipole strength distributions for a variety of reasons including the extraction of the dipole polarizability related to properties of the symmetry energy and a measure for the neutron skin thickness, understanding the structure of low-energy $\textit{E}$1 strength in nuclei with neutron excess, and establishing the systematics of the isovector giant dipole resonance (IVGDR). Inelastic proton scattering at energies of a few hundred MeV and very forward angles including 0° has been established as a tool for the study of electric and magnetic dipole strength distributions in nuclei. Purpose: The present work aims at a systematic investigation of the electric and magnetic dipole strength distributions in the chain of stable even-mass tin isotopes. Methods: Inelastic proton scattering experiments were performed at the Research Center for Nuclear Physics, Osaka, with a 295-MeV beam covering laboratory angles 0°–6° and excitation energies 6–22 MeV. Cross sections due to $\textit{E}$1 and $\textit{M}$1 excitations were extracted with a multipole decomposition analysis (MDA) and then converted to reduced transition probabilities with the “virtual photon method” for $\textit{E}$1 and the “unit cross section method” for $\textit{M}$1 excitations, respectively. Including a theory-aided correction for the high-excitation-energy region not covered experimentally, the electric dipole polarizability was determined from the $\textit{E}$1 strength distributions. Results: In this work, total photoabsorption cross sections derived from the $\textit{E}$1 and $\textit{M}$1 strength distributions show significant differences compared to those from previous ($\textit{γ, xn}$) experiments in the energy region of the IVGDR. The widths of the IVGDR deduced from the present data with a Lorentz parametrization show an approximately constant value of about 4.5 MeV in contrast to the large variations between isotopes observed in previous work. The IVGDR centroid energies are in good correspondence to expectations from empirical systematics of their mass dependence. Furthermore, a study of the dependence of the IVGDR energies on bulk matter properties is presented. The $\textit{E}$1 strengths below neutron threshold show fair agreement in this work with results from ($\textit{γ, γ'}$) experiments on 112,116,120,124 Sn in the energy region between 6 and 7 MeV, where also isoscalar $\textit{E}$1 strength was found for 124 Sn. At higher excitation energies, large differences are observed, pointing to a different nature of the excited states with small ground-state branching ratios. The isovector spin-$\textit{M}$1 strengths exhibit a broad distribution between 6 and 12 MeV in all studied nuclei. Conclusions: The present results contribute to the solution of a variety of nuclear structure problems including the systematics of the energy and width of the IVGDR, the structure of low-energy $\textit{E}$1 strength in nuclei, new constraints to energy density functionals (EDFs) aiming at a systematic description of the dipole polarizability across the nuclear chart, from which properties of the symmetry energy can be derived, and the systematics of the isovector spin-$\textit{M}$1 strength in heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗