Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source and propagation modeling”

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 289 records · Page 16

Actinide cosmic ray chronometers - Relative abundances and the cosmic ray lifetime

Measurements of ultraheavy cosmic rays will provide radioactive chronometers, including the actinides, which are expected to be present in the source. These primary chronometers differ from the Be-10 secondary chronometer in total inelastic cross-section and decay half-life. A propagation code, which includes nuclear fragmentation, radioactive decay, and ionization energy-loss of cosmic rays is being used to investigate the variations in the fluxes and mean ages which would be expected for various models of cosmic ray propagation. Some preliminary results in the actinide region are presented here which indicate that measurements of the relative abundances of the actinides in the cosmic rays will be very useful for understanding source abundances but less useful for studying propagation effects.

Margolis, S. H.↗

Uncertainty propagation in a stratospheric model. I - Development of a concise stratospheric model. II - Monte Carlo analysis of imprecisions due to reaction rates

The paper discusses the development of a concise stratospheric model which uses iteration to obtain coupling between interacting species. The one-dimensional, steady-state, diurnally-averaged model generates diffusion equations with appropriate sources and sinks for species odd oxygen, H2O, H2, CO, N2O, odd nitrogen, CH4, CH3Cl, CCl4, CF2Cl2, CFCl3, and odd chlorine. The model evaluates steady-state perturbations caused by injections of chlorine and NO(x) and may be used to predict ozone depletion. The model is used in a Monte Carlo study of the propagation of reaction-rate imprecisions by calculating an ozone perturbation caused by the addition of chlorine. Since the model is sensitive to only 10 of the more than 50 reaction rates considered, only about 1000 Monte Carlo cases are required to span the space of possible results.

Rundel, R. D.↗

Comprehensive framework for data-driven model form discovery of the closure laws in thermal-hydraulics codes

The two-phase two-fluid model is a basis of many thermal-hydraulics codes used in design, licensing, and safety considerations of nuclear power plants. Thermal-hydraulics codes rely on the closure laws to close the system of conservation equations and describe the interactions between phases. These laws, derived from years of experimental investigations, are semi-empirical correlations that lack generality and have a limited range of applicability. Increase of computational power, availability of new experiments, and development of high-fidelity simulations has increased the number of validation data. The discrepancies between the code predictions and the validation data are a great source of knowledge. Missing physics that are not included in the model but are important for the considered phenomena can be discovered by propagating the information from the experimental results through the model. Furthermore, physics-discovered data-driven model form (P3DM) methodology integrates available integral effect tests and separate effects tests to determine the necessary corrections to the model form of the closure laws. In contrast to existing calibration techniques, the methodology modifies the functional form of the closure laws. Based on the functional form of the correction, the missing physics that were not included in the original model can be discovered. The methodology provides the alternative to the machine learning approach, in which the model is discovered in the form of the intractable black-box relation. In this work, the methodology was applied to the CTF subchannel code to improve the prediction of the two-phase flow phenomena.

42 ENGINEERING↗

Dynamic Rupture Simulations of Caldera Collapse Earthquakes: Effects of Wave Radiation, Magma Viscosity, and Evidence of Complex Nucleation at Kilauea 2018

All instrumented basaltic caldera collapses have generated M w > 5 very long period earthquakes. However, previous studies of source dynamics have been limited to lumped models treating the caldera block as rigid, leaving open questions related to how ruptures initiate and propagate around the ring fault, and the seismic expressions of those dynamics. Here, we present the first 3D numerical model capturing the nucleation and propagation of ring fault rupture, the mechanical coupling to the underlying viscoelastic magma, and the associated seismic wavefield. We demonstrate that seismic radiation, neglected in previous models, acts as a damping mechanism reducing coseismic slip by up to half, with effects most pronounced for large magma chamber volume/ring fault radius or highly compliant crust/compressible magma. Viscosity of basaltic magma has negligible effect on collapse dynamics. In contrast, viscosity of silicic magma significantly reduces ring fault slip. We use the model to simulate the 2018 Kilauea caldera collapse. Three stages of collapse, characterized by ring fault rupture initiation and propagation, deceleration of the downward-moving caldera block and magma column, and post-collapse resonant oscillations, in addition to chamber pressurization, are identified in simulated and observed (unfiltered) near-field seismograms. A detailed comparison of simulated and observed displacement waveforms corresponding to collapse earthquakes with hypocenters at various azimuths of the ring fault reveals a complex nucleation phase for earthquakes initiated on the northwest. Our numerical simulation framework will enhance future efforts to reconcile seismic and geodetic observations of caldera collapse with conceptual models of ring fault and magma chamber dynamics.

58 GEOSCIENCES↗

Seismoacoustic analysis of underground explosions using the Rayleigh integral

SUMMARY A coupled seismoacoustic model is developed for the analysis of acoustic signals produced by underground explosive events with an aim to develop a means of improving estimated depth and yield for explosion monitoring. A ground spall model is used to predict surface motion characteristics produced by an underground explosion and the Rayleigh integral is applied to relate the surface motion to the acoustic signal some distance from surface ground zero. The low-frequency component of the ground motion associated with the prolonged free fall of lofted material during spall is found to dominate the acoustic signal propagating away from surface ground zero at shallow angles. The model is applied to observed ground motion and acoustic signals recorded during the Source Physics Experiment (SPE) with promising results. In addition to accurately predicting characteristics of the observed acoustic signals from several of the SPE events, the model provides a means of explaining the lack of signals observed during several events in the SPE due to the directionality of the higher frequency acoustic signals associated with the uplift and closure components as well as the lack of a spatially localized, longer duration dwell.

58 GEOSCIENCES↗

SIRU utilization. Volume 1: Theory, development and test evaluation

The theory, development, and test evaluations of the Strapdown Inertial Reference Unit (SIRU) are discussed. The statistical failure detection and isolation, single position calibration, and self alignment techniques are emphasized. Circuit diagrams of the system components are provided. Mathematical models are developed to show the performance characteristics of the subsystems. Specific areas of the utilization program are identified as: (1) error source propagation characteristics and (2) local level navigation performance demonstrations.

Musoff, H.↗

Measurement of the cosmic-ray iron spectrum between 60 and 200 GeV per nucleon

A measurement of the spectral index of Galactic cosmic-ray (GCR) iron has been made using a high-energy gas Cerenkov spectrometer. The spectral index of GCR iron is found to be 2.56 + or - 0.11 in the energy range 57-200 GeV/ nucleon. This result indicates that the source spectrum of GCR iron is similar to that of other primary GCR nuclei and is consistent with the simplest models of GCR propagation.

Esposito, Joseph A.↗

Microwave Driven Actuators Power Allocation and Distribution

Design, fabrication and test of a power allocation and distribution (PAD) network for microwave driven actuators is presented in this paper. Development of a circuit that would collect power from a rectenna array amplify and distribute the power to actuators was designed and fabricated for space application in an actuator array driven by a microwave. A P-SPICE model was constructed initially for data reduction purposes, and was followed by a working real-world model. A voltage up - converter (VUC) is used to amplify the voltage from the individual rectenna. The testing yielded a 26:1 voltage amplification ratio with input voltage at 9 volts and a measured output voltage 230VDC. Future work includes the miniaturization of the circuitry, the use of microwave remote control, and voltage amplification technology for each voltage source. The objective of this work is to develop a model system that will collect DC voltage from an array of rectenna and propagate the voltage to an array of actuators.

Forbes, Timothy↗

The elemental and isotopic composition of Galactic cosmic-ray nuclei from scandium through nickel

Measurements of the relative elemental and isotopic abundances of iron-group Galactic cosmic rays at energies of about 325 MeV per nucleon have been made. The source abundance ratio of Ni-60/Ni-58 is 1.07 +/- 0.39, which is a factor of 2.8 +/- 1.0 larger than the solar system value. Our measurements imply the presence of Co-59 at the source, which can be reconciled with the predictions of conventional nucleosynthesis models if there exists a time delay of more than about 100,000 yr between nucleosynthesis and acceleration. Most of the Mn-54 produced by spallation during cosmic-ray propagation in the Galaxy is found to have decayed to Fe-54, indicating a confinement time of greater than 2 Myr. The source ratio of Fe-54/Fe-56 corrected for the Mn-54 decay is 0.046+/- 0.020, which is consistent with the solar system value of 0.063.

Leske, Richard A.↗

The theory of magnetohydrodynamic wave generation by localized sources. I - General asymptotic theory

The magnetohydrodynamic wave emission from several localized, periodic, kinematically specified fluid velocity fields are calculated using Lighthill's method for finding the far-field wave forms. The waves propagate through an isothermal and uniform plasma with a constant B field. General properties of the energy flux are illustrated with models of pulsating flux tubes and convective rolls. Interference theory from geometrical optics is used to find the direction of minimum fast-wave emission from multipole sources and slow-wave emission from discontinuous sources. The distribution of total flux in fast and slow waves varies with the ratios of the source dimensions l to the acoustic and Alfven wavelengths.

Collins, William↗

VpROM: a novel variational autoencoder-boosted reduced order model for the treatment of parametric dependencies in nonlinear systems

Reduced Order Models (ROMs) are of considerable importance in many areas of engineering in which computational time presents difficulties. Established approaches employ projection-based reduction, such as Proper Orthogonal Decomposition. The limitation of the linear nature of such operators is typically tackled via a library of local reduction subspaces, which requires the assembly of numerous local ROMs to address parametric dependencies. Our work attempts to define a more generalisable mapping between parametric inputs and reduced bases for the purpose of generative modeling. We propose the use of Variational Autoencoders (VAEs) in place of the typically utilised clustering or interpolation operations, for inferring the fundamental vectors, termed as modes, which approximate the manifold of the model response for any and each parametric input state. The derived ROM still relies on projection bases, built on the basis of full-order model simulations, thus retaining the imprinted physical connotation. However, it additionally exploits a matrix of coefficients that relates each local sample response and dynamics to the global phenomena across the parametric input domain. The VAE scheme is utilised for approximating these coefficients for any input state. This coupling leads to a high-precision low-order representation, which is particularly suited for problems where model dependencies or excitation traits cause the dynamic behavior to span multiple response regimes. Moreover, the probabilistic treatment of the VAE representation allows for uncertainty quantification on the reduction bases, which may then be propagated to the ROM response. The performance of the proposed approach is validated on an open-source simulation benchmark featuring hysteresis and multi-parametric dependencies, and on a large-scale wind turbine tower characterised by nonlinear material behavior and model uncertainty.

Conditional VAEs↗

Measurements of the spatial structure and directivity of 100 KeV photon sources in solar flares using PVO and ISEE-3 spacecraft

The objective of this grant was to measure the spatial structure and directivity of the hard X-ray and low energy gamma-ray (100 keV-2 MeV) continuum sources in solar flares using stereoscopic observations made with spectrometers aboard the Pioneer Venus Orbiter (PVO) and Third International Sun Earth Explorer (ISEE-3) spacecraft. Since the hard X-ray emission is produced by energetic electrons through the bremsstrahlung process, the observed directivity can be directly related to the 'beaming' of electrons accelerated during the flare as they propagate from the acceleration region in the corona to the chromosphere/transition region. Some models (e.g., the thick-target model) predict that most of the impulsive hard X-ray/low energy gamma-ray source is located in the chromosphere, the effective height of the X-ray source above the photosphere increasing with the decrease in the photon energy. This can be verified by determining the height-dependence of the photon source through stereoscopic observations of those flares which are partially occulted from the view of one of the two spacecraft. Thus predictions about beaming of electrons as well as their spatial distributions could be tested through the analysis proposed under this grant.

Anderson, Kinsey A.↗

Full Scale Rotor Aeroacoustic Predictions and the Link to Model Scale Rotor Data

The NASA Aeroacoustic Prediction System (NAPS) is used to establish a link between model-scale and full-scale rotor predictions and is partially validated against measured wind tunnel and flight aeroacoustic data. The prediction approach of NAPS couples a comprehensive rotorcraft analysis with acoustic source noise and propagation codes. The comprehensive analysis selected for this study is CAMRAD-II, which provides the performance/trim/wake solution for a given rotor or flight condition. The post-trim capabilities of CAMRAD-II are used to compute high-resolution sectional airloads for the acoustic tone noise analysis, WOPMOD. The tone noise is propagated to observers on the ground with the propagation code, RNM (Rotor Noise Model). Aeroacoustic predictions are made with NAPS for an isolated rotor and compared to results of the second Harmonic Aeroacoustic Rotor Test (HART-II) program, which tested a 40% dynamically and Mach-scaled BO-105 main rotor at the DNW. The NAPS is validated with comparisons for three rotor conditions: a baseline condition and two Higher Harmonic Control (HHC) conditions. To establish a link between model and full-scale rotor predictions, a full-scale BO-105 main rotor input deck for NAPS is created from the 40% scale rotor input deck. The full-scale isolated rotor predictions are then compared to the model predictions. The comparisons include aerodynamic loading, acoustic levels, and acoustic pressure time histories for each of the three conditions. With this link established, full-scale predictions are made for a range of descent flight conditions and compared with measured trends from the recent Rotorcraft Operational Noise Abatement Procedures (RONAP) flight test conducted by DLR and ONERA. Additionally, the effectiveness of two HHC conditions from the HART-II program is demonstrated for the full-scale rotor in flight.

Boyd, D. Douglas, Jr.↗

Near-global summer circulation response to the spring surface temperature anomaly in Tibetan Plateau –– the GEWEX/LS4P first phase experiment

Subseasonal to seasonal (S2S) prediction of droughts and floods is one of the major challenges of weather and climate prediction. Recent studies suggest that the springtime land surface temperature/subsurface temperature (LST/SUBT) over the Tibetan Plateau (TP) can be a new source of S2S predictability. The project “Impact of Initialized Land Surface Temperature and Snowpack on Subseasonal to Seasonal Prediction (LS4P)” was initiated to study the impact of springtime LST/SUBT anomalies over high mountain areas on summertime precipitation predictions. The present work explores the simulated global scale response of the atmospheric circulation to the springtime TP land surface cooling by 16 current state-of-the-art Earth System Models (ESMs) participating in the LS4P Phase I (LS4P-I) experiment. The LS4P-I results show, for the first time, that springtime TP surface anomalies can modulate a persistent quasi-barotropic Tibetan Plateau-Rocky Mountain Circumglobal (TRC) wave train from the TP via the northeast Asia and Bering Strait to the western part of the North America, along with the springtime westerly jet from TP across the whole North Pacific basin. The TRC wave train modulated by the TP thermal anomaly play a critical role on the early summer surface air temperature and precipitation anomalies in the regions along the wave train, especially over the northwest North America and the southern Great Plains. The participant models that fail in capturing the TRC wave train greatly under-predict climate anomalies in reference to observations and the successful models. These results suggest that the TP LST/SUBT anomaly via the TRC wave train is the first order source of the S2S variability in the regions mentioned. Furthermore, the TP surface temperature anomaly can influence the Southern Hemispheric circulation by generating cross-equator wave trains. However, the simulated propagation pathways from the TP into the Southern Hemisphere show large inter-model differences. More dynamical understanding of the TRC wave train as well as its cross-equator propagation into the Southern Hemisphere will be explored in the newly launched LS4P phase II experiment.

54 ENVIRONMENTAL SCIENCES↗

Gravity Waves and Mesospheric Clouds in the Summer Middle Atmosphere: A Comparison of Lidar Measurements and Ray Modeling of Gravity Waves Over Sondrestrom, Greenland

We conducted gravity wave ray-tracing experiments within an atmospheric region centered near the ARCLITE lidar system at Sondrestrom, Greenland (67N, 310 deg E), in efforts to understand lidar observations of both upper stratospheric gravity wave activity and mesospheric clouds during August 1996 and the summer of 2001. The ray model was used to trace gravity waves through realistic three-dimensional daily-varying background atmospheres in the region, based on forecasts and analyses in the troposphere and stratosphere and climatologies higher up. Reverse ray tracing based on upper stratospheric lidar observations at Sondrestrom was also used to try to objectively identify wave source regions in the troposphere. A source spectrum specified by reverse ray tracing experiments in early August 1996 (when atmospheric flow patterns produced enhanced transmission of waves into the upper stratosphere) yielded model results throughout the remainder of August 1996 that agreed best with the lidar observations. The model also simulated increased vertical group propagation of waves between 40 km and 80 km due to intensifying mean easterlies, which allowed many of the gravity waves observed at 40 km over Sondrestrom to propagate quasi-vertically from 40-80 km and then interact with any mesospheric clouds at 80 km near Sondrestrom, supporting earlier experimentally-inferred correlations between upper stratospheric gravity wave activity and mesospheric cloud backscatter from Sondrestrom lidar observations. A pilot experiment of real-time runs with the model in 2001 using weather forecast data as a low-level background produced less agreement with lidar observations. We believe this is due to limitations in our specified tropospheric source spectrum, the use of climatological winds and temperatures in the upper stratosphere and mesosphere, and missing lidar data from important time periods.

LIDAR (LIGHT DETECTION AND RANGING)↗

Field-independent source localization of Neptune's radio bursts

During the Voyager 2 encounter with Neptune, a narrowbanded bursty radio component was observed between 500 and 1326 kHz by the Planetary Radio Astronomy instrument. Based on the emission occurrence pattern, the radio source has been localized without the explicit use of the Neptunian offset-tilted dipole magnetic field model, which is accurate only at distances greater than 4 R(N) (Neptune radii) from the planet. Only assumptions based upon the general nature of radio wave propagation in planetary magnetospheres were used. A number of different candidate radial positions were sampled. For example, at 1.5 R(N), the derived source location was positioned only about 10 deg from the south magnetic pole. The radiation from this source was beamed into a cone of 77.5 + or - 6.3 deg half-angle that was tilted about 10 deg from the radial direction to the north-northeast. At other sampled radial positions, similar source locations were obtained. Due to its proximity to the south magnetic pole, the kilometric emission radio source is believed to be associated with an active auroral region, similar in nature to those found at earth and Saturn.

Farrell, W. M.↗

Observation of high-energy neutrinos from the Galactic plane

The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth’s atmosphere, is unknown. Because of deflection by interstellar magnetic fields, cosmic rays produced within the Milky Way arrive at Earth from random directions. However, cosmic rays interact with matter near their sources and during propagation, which produces high-energy neutrinos. We searched for neutrino emission using machine learning techniques applied to 10 years of data from the IceCube Neutrino Observatory. By comparing diffuse emission models to a background-only hypothesis, we identified neutrino emission from the Galactic plane at the 4.5σ level of significance. The signal is consistent with diffuse emission of neutrinos from the Milky Way but could also arise from a population of unresolved point sources.

Science & Technology - Other Topics↗

A thundercloud electric field sounding - Charge distribution and lightning

An instrumented free balloon measured electric fields and field changes as it rose through a thundercloud above Langmuir Laboratory, New Mexico. The variation of the electric field with altitude implied that the cloud contained negative space charge of density -0.6 to -4 nC/cu m between 5.5 and 8.0 km MSL. The environmental temperature at these levels ranged from -5 to -20 C. The measurements imply that the areal extent of this negative charge center was significantly greater than that of the cloud's intense precipitation shafts. At altitudes greater than 8 km, the instrument ascended past net positive charge. In addition, positive space charge adjacent to the earth's surface (concentration 0.6 nC/cu m and in the lowest portion of the cloud (1.0 nC/cu m) is inferred from the measurements. Electric field changes from intracloud lightning were interpreted by using a simple model for the developing streamer of the initial phase. Thunder source reconstructions provided estimates for the orientation of lightning channels. Seven 'streamers' so analyzed propagated on the average, at 50,000 m/s and carried a current of 390 A. The mean charge dissipated during a flash was 30 C.

Weber, M. E.↗