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 433 records · Page 24

Limits on coronal reflection using high-frequency solar oscillations

Acoustic waves in the Sun with frequencies above about 5.3 mHz can propagate in the chromosphere. We examine imaged solar intensity data for evidence of reflection of these waves in the upper chromosphere, where the temperature increases by a large factor over a short distance. Our method is to compare the observed and theoretically derived frequency spacings between peaks in the power spectrum. We find that our theoretical frequencies provide the best fit to the data when the reflection in the upper atmosphere is eliminated. In particular, the model of Kumar (1993b), which includes the source depth, and radiative damping, in the calculation of power spectra but ignores chromospheric reflection, gives peak frequencies that are in good agreement with the observations. For acoustic waves of frequency greater than 6 mHz we put an upper limit to the reflectivity of chromosphere and corona, using our method, of about 10%. At a given spherical harmonic degree, the frequency spacing between peaks in the data generally decreases with increasing frequency, because the lower turning point of the waves is moving inward. However, between 5 and 5.5 mHz the frequency spacing increases slightly. This feature is probably associated with the acoustic cutoff frequency in the solar atmosphere, i.e., it indicates a transition from trapped waves to propagating waves. We are able to reproduce the observed behavior by a crude modeling of the solar atmosphere. Further study of these peaks should provide an independent way of exploring the mean structure of the solar atmosphere, particularly around the temperature minimum region.

Kumar, P.↗

The ionization of the diffuse ionized gas

The diffuse ionized gas (DIG), responsible for producing the faint but pervasive H(sub alpha), requires 15% of the power of all Galactic O stars for its ionization. No other source of ionization seems practical. The spectrum is strong in low stages of ionization (N II, SII), weak in (O III), and very weak in (O I) lambda 6300, at least in one well-observed direction, all significantly different from denser H II regions. We give low-density, low-excitation photoionization models that explain the observed spectrum. The lambda 6300 observation eliminates the simplest models for the DIG (Mathis 1986), in which neutral H extends beyond the edge of the very dilute stellar radiation field. Our present models include two components. One, representing the edges of interstellar H I clouds, extends to the point where H becomes neutral. In the second, the fraction of H(sup o) is not allowed to exceed 0.05 to 0.1. Both have very low values of the ionization parameter, or ratio of the number densities of ionizing photons to electrons. The ionization parameter required by our models is shown to be compatible with observed values of photon and electron densities in the diffuse interstellar medium. Interstellar dust is not important in either model. Predictions of the model are that the Galaxy is leaking about 4% of the ionizing radiation of at least the local O stars into the intergalactic medium, that (O II) lambda 3727 approximately equals 1.1 H(sub alpha), and that He(sup +)/H (sup +)) approximately equals 0.6 He/H. A major difficul ty of this picture is having ionizing radiation propagate from the ionizing stars in the plane of the Galaxy to the DIG. We suggest, as do Miller & Cox (1993), that the ionizing radiation from O stars produces holes or tunnels in the distribution of neutral H as seen from their positions, while we see a sky covered with H I sheets or filaments because we are not near a source of ionizing radiation. With our models, several interfaces between ionized and neutral gas are allowed before the observational limit on lambda 6300 is exceeded. Another difficulty is the gas-phase N(S(sup +))/N(H(sup o)) ratio toward the high-latitude star HD 93521 (Spitzer & Fitzpatrick 1993). We require (S/H) to be somewhat lower than solar, either because of depletion onto grains or a Galactic abundance gradient, so that the S(sup +) can be produced in each of separate H(sup o) and H(sup +) regions.

Domgorgen, H.↗

Testing Special Relativity at High Energies with Astrophysical Sources

Since the group of Lorentz boosts is unbounded, there is a question as to whether Lorentz invariance (LI) holds to infinitely short distances. However, special and general relativity may break down at the Planck scale. Various quantum gravity scenarios such as loop quantum gravity, as well as some forms of string theory and extra dimension models may imply Lorentz violation (LV) at ultrahigh energies. The Gamma-Ray Large Area Space Telescope (GLAST), to be launched in mid-December, will measure the spectra of distant extragalactic sources of high energy gamma-rays, particularly active galactic nuclei and gamma-ray bursts. GLAST can look for energy-dependent gamma-ray propagation effects from such sources as a signal of Lorentz invariance violation. These sources may also exhibit the high energy cutoffs predicted to be the result of intergalactic annihilation interactions with low energy photons having a flux level as determined by various astronomical observations. With LV the threshold for such interactions can be significantly raised, changing the predicted absorption turnover in the observed spectrum of the sources. Stecker and Glashow have shown that the existence such absorption features in the spectra of extragalactic sources puts constraints on LV. Such constraints have important implications for some quantum gravity and large extra dimension models. Future spaceborne detectors dedicated to measuring gamma-ray polarization can look for birefringence effects as a possible signal of loop quantum gravity. A very small LV may also result in the modification or elimination of the GZK effect, thus modifying the spectrum of ultrahigh energy cosmic rays. This possibility can be explored with ground-based arrays such as Auger or with a space based detector system such as the proposed OWL satellite mission.

Stecker, F. W.↗

Production and propagation of particles with A greater than 81 in the Galaxy.

A detailed mathematical model is developed to describe the transformation of the charge composition of superheavy (A greater than 81) cosmic rays by spallation on interstellar hydrogen in the Galaxy. If one takes a single source which suddenly injects particles with an r-process charge spectrum t years in the past, the calculated relative abundance ratios in four atomic mass groups are consistent with measured values for t equal to 250,000 years. We further conclude that for any initial source distribution the upper limit on t is 1 m.y.

Kaiser, T. B.↗

The unusual radio outburst of Nova Vulpeculae 1984 No. 2

Multifrequency radio monitoring of Nova Vulpeculae 1984 No. 2 has revealed a unique radio light curve, exhibiting a strong outburst which precedes the appearance of normal radio emission from the principal ejecta of the nova by at least 100 days. The early emission is extremely optically-thick and has brightness temperature in excess of 100,000 K. A model is discussed in which the radiation is produced by a strong shock propagating outward through the principal ejecta of the nova, as a result of an interaction with a later, high velocity wind from the central source. It is shown that the general features of the radio light curve can be explained by the presence of a central wind with a mass loss rate of about 0.00001 solar masses/yr, lasting for a period of 200 to 300 days after the optical outburst. The first radio map of nova ejecta shortly after outburst is presented. Comparison of the angular expansion rate to the ejection velocity implies a distance of 3.6 kpc, and preliminary analysis of the main outburst emission indicates a mass and kinetic energy of the ejecta respectively of about 0.0008 solar masses and 8 x 10 to the 45th erg.

Taylor, A. R.↗

The solar wind - Advances in our knowledge through two solar cycles

As the Pioneer and Voyager spacecraft have moved outward they have gradually unfolded a view of distant regions of the heliosphere. Information on the solar wind velocity, density and temperature as a function of distance out to more than 40 AU has been gathered. Meanwhile the description of the solar wind has evolved. Long-standing questions on the sources of the wind causing geomagnetic activity were clarified by the discovery of coronal holes and coronal mass ejections. The propagation of the resultant solar wind disturbances through the heliosphere has been studied using both observations and models. Plasma physical processes have been studied. This review focuses on the development of the concepts that have been used to describe the solar wind in the three dimensional heliosphere over the last two solar cycles. Collisionless shocks, transient disturbances in space, disturbance propagation and the distant solar wind are discussed.

Feynman, Joan↗

Comparison of distributed acceleration and standard models of cosmic-ray transport

Recent cosmic-ray abundance measurements for elements in the range 3 less than or equal to Z less than or equal to 28 and energies 10 MeV/n less than or equal to E less than or equal to 1 TeV/n have been analyzed with computer transport modeling. About 500 elemental and isotopic measurements have been explored in this analysis. The transport code includes the effects of ionization losses, nuclear spallation reactions (including those of secondaries), all nuclear decay modes, stripping and attachment of electrons, escape from the Galaxy, weak reacceleration and solar modulation. Four models of reacceleration (with several submodels of various reacceleration strengths) were explored. A chi (exp 2) analysis show that the reacceleration models yield at least equally good fits to the data as the standard propagation model. However, with reacceleration, the ad hoc assumptions of the standard model regarding discontinuities in the energy dependence of the mean path length traversed by cosmic rays, and in the momentum spectrum of the cosmic-ray source spectrum are eliminated. Futhermore, the difficulty between rigidity dependent leakage and energy independent anisotropy below energies of 10(exp 14) eV is alleviated.

Letaw, J. R.↗

Physics of Acoustic Radiation from Jet Engine Inlets

Numerical simulations of acoustic radiation from a jet engine inlet are performed using advanced computational aeroacoustics (CAA) algorithms and high-quality numerical boundary treatments. As a model of modern commercial jet engine inlets, the inlet geometry of the NASA Source Diagnostic Test (SDT) is used. Fan noise consists of tones and broadband sound. This investigation considers the radiation of tones associated with upstream propagating duct modes. The primary objective is to identify the dominant physical processes that determine the directivity of the radiated sound. Two such processes have been identified. They are acoustic diffraction and refraction. Diffraction is the natural tendency for an acoustic wave to follow a curved solid surface as it propagates. Refraction is the turning of the direction of propagation of sound waves by mean flow gradients. Parametric studies on the changes in the directivity of radiated sound due to variations in forward flight Mach number and duct mode frequency, azimuthal mode number, and radial mode number are carried out. It is found there is a significant difference in directivity for the radiation of the same duct mode from an engine inlet when operating in static condition and in forward flight. It will be shown that the large change in directivity is the result of the combined effects of diffraction and refraction.

Tam, Christopher K. W.↗

Rapid Prediction of Installed Jet Noise from RANS

A new method of computing jet noise, called mSrc, was developed on the general principles of acoustic analogies. In the method, the problem of translating turbulent flow energy into acoustic energy at a far-field observer is broken into two parts, a calculation of acoustic source strengths and then their propagation. The acoustic sources are related to turbulent quantities in the jet plume in a robust manner. The propagation, which is more properly computed using a Green's function accounting for nonuniform speed of sound and solid surfaces, is instead modeled using commonly observed features of jet acoustic far field directivity, and by diffraction barrier theory for surfaces. The mSrc method does not require gradients of the predicted flow field, allowing it to make use of robust unstructured RANS CFD methods, including embedded boundary codes. Such codes do not require specification of surface meshes, and auto-refine their grid to resolve flow gradients, putting resolution where it is required without a priori user input. The ability of mSrc to use such radically unstructured flow input results in an efficient method of estimating noise from jet flows from complex nozzles installed on aircraft. Many validation cases are presented to demonstrate the accuracy and range of applicability of the mSrc method for representative jet noise applications.

Bridges, James↗

Rapid Prediction of Installed Jet Noise From RANS

A new method of computing jet noise, called mSrc, was developed on the general principles of acoustic analogies. In the method, the problem of translating turbulent flow energy into acoustic energy at a far-field observer is broken into two parts, a calculation of acoustic source strengths and then their propagation. The acoustic sources are related to turbulent quantities in the jet plume in a robust manner. The propagation, which is more properly computed using a Green's function accounting for nonuniform speed of sound and solid surfaces, is instead modeled using commonly observed features of jet acoustic far field directivity, and by diffraction barrier theory for surfaces. The mSrc method does not require gradients of the predicted flow field, allowing it to make use of robust unstructured RANS CFD methods, including embedded boundary codes. Such codes do not require specification of surface meshes, and auto-refine their grid to resolve flow gradients, putting resolution where it is required without a priori user input. The ability of mSrc to use such radically unstructured flow input results in an efficient method of estimating noise from jet flows from complex nozzles installed on aircraft. Many validation cases are presented to demonstrate the accuracy and range of applicability of the mSrc method for representative jet noise applications.

Bridges, James↗

Sources of Uncertainty in the Prediction of LAI / fPAR from MODIS

To explicate the sources of uncertainty in the prediction of biophysical variables over space, consider the general equation: where z is a variable with values on some nominal, ordinal, interval or ratio scale; y is a vector of input variables; u is the spatial support of y and z ; x and u are the spatial locations of y and z , respectively; f is a model and B is the vector of the parameters of this model. Any y or z has a value and a spatial extent which is called its support. Viewed in this way, categories of uncertainty are from variable (e.g. measurement), parameter, positional. support and model (e.g. structural) sources. The prediction of Leaf Area Index (LAI) and the fraction of absorbed photosynthetically active radiation (fPAR) are examples of z variables predicted using model(s) as a function of y variables and spatially constant parameters. The MOD15 algorithm is an example of f, called f(sub 1), with parameters including those defined by one of six biome types and solar and view angles. The Leaf Canopy Model (LCM)2, a nested model that combines leaf radiative transfer with a full canopy reflectance model through the phase function, is a simpler though similar radiative transfer approach to f(sub 1). In a previous study, MOD15 and LCM2 gave similar results for the broadleaf forest biome. Differences between these two models can be used to consider the structural uncertainty in prediction results. In an effort to quantify each of the five sources of uncertainty and rank their relative importance for the LAI/fPAR prediction problem, we used recent data for an EOS Core Validation Site in the broadleaf biome with coincident surface reflectance, vegetation index, fPAR and LAI products from the Moderate Resolution Imaging Spectrometer (MODIS). Uncertainty due to support on the input reflectance variable was characterized using Landsat ETM+ data. Input uncertainties were propagated through the LCM2 model and compared with published uncertainties from the MOD15 algorithm.

Dungan, Jennifer L.↗

An Extension of the Truncated Ingard-Myers Impedance Boundary Condition for High Mach Number Grazing Flows

The Ingard-Myers impedance boundary condition is widely recognized to be prone to hydrodynamic Kelvin-Helmholtz-type instability, primarily due to its use of a vortex sheet to model the flow at the boundary. Recently, a stabilized approximation of the Ingard-Myers condition, called the Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC), has been introduced. This approximation offers a good representation of the Ingard-Myers condition for grazing flows characterized by low to mid subsonic Mach numbers. This paper explores an extension of the TIMIBC, referred to as TIMIBC-ext, designed for liners in grazing flows with high subsonic Mach numbers. The TIMIBC-ext introduces a tunable parameter and demonstrates that by judiciously selecting the parameter value, the accuracy of the TIMIBC can be improved for either upstream or downstream propagating waves in grazing flows with high subsonic Mach numbers. Consequently, the TIMIBC-ext represents an enhancement to the TIMIBC when the location of the source relative to the liner (upstream or downstream) is known. Time domain implementations of the TIMIBC-ext are also provided, employing a multipole expansion model for the impedance function. A numerical example is provided to assess the performance of the TIMIBC-ext, comparing computational results with measurements from a recent set of measurements in the NASA Langley Grazing Flow Impedance Tube.

time domain acoustic scattering↗

Prediction Skill of the MJO Teleconnection Signals in the NASA GEOS Subseasonal Reforecasts

Tropical-extratropical teleconnections are considered key to advancing subseasonal prediction. The Madden Julian oscillation (MJO), characterized by large scale convective envelopes propagating along the tropical Indo-Pacific sector, is known to modulate midlatitude circulation and associated weather patterns. Although there is a general consensus on the MJO's influence on the midlatitude circulation, which is thought to be due to modulations of the North Atlantic Oscillation (NAO) and the Pacific North America (PNA) pattern, relatively less is known about the predictability of these teleconnection signals in dynamical forecast models. The composite evolution of the midlatitude circulation anomalies and associated wave train structure as the delayed response to tropical heating are reported in many studies that have examined reanalyses and long model simulations. However, it is yet to be determined whether they lend any beneficial subseasonal forecast skill, especially to weekly mean surface temperature and precipitation over North America. Investigating useful predictable signals from the MJO teleconnections is also complicated by the fact that the MJO is a moving heat source with an approximate periodicity of 30-60 days, and that the structure and amplitude of the midlatitude response can be sensitive to the longitudinal positioning of the heating anomaly as well as the propagation speed of the MJO. The objective of this study is to investigate the impact of MJO teleconnections on forecast accuracy at 2-3 week lead over North America, with an emphasis on the above-mentioned lesser known aspects of these teleconnections. To this end, we utilize a suite of subseasonal reforecasts performed with the latest NASA GEOS-5 seasonal-to-subseasonal (S2S) system. These reforecasts were performed as part of the NOAA SubX project, wherein the NASA GEOS-5 atmosphere-ocean coupled model was run at degree horizontal resolution, initialized every 5 days for the period 1999-2016. The GEOS-5 model shows skillful predictions of the MJO, with the correlation coefficient based on the real-time multivariate MJO (RMM) index staying at or above 0.5 up to forecast lead 26-36 days. The system is thus a useful tool for investigating MJO teleconnection processes.

Achuthavarier, Deepthi↗

Geologic Framework Model for the Dry Alluvium Geology (DAG) Experiment Testbed Yucca Flat, Nevada National Security Site

Geologic framework models (GFMs) provide a methodology for integrating geology into, and thus geologically informing, other modeling and simulation activities. GFMs provide a threedimensional (3-D), geology-based digital framework for developing and parametrizing meshes and evaluating model and simulation results. This report describes a 3-D GFM constructed for the Source Physics Experiment Phase II Dry Alluvium Geology test series located in Yucca Flat at the Nevada National Security Site. The geology in the Yucca Flat region is complex and diverse, which creates challenges to modeling seismic wave propagation from SPE tests. The Yucca Flat GFM helps address these challenges by providing the 3-D distribution of relevant geologic features and physical properties necessary to more effectively model seismic wave propagation. The GFM includes 7 model layers and 48 faults that cut and offset the layers. An appendix is included that provides quantitative data on physical properties for each model layer.

58 GEOSCIENCES↗

Geologic framework model for the Dry Alluvium Geology (DAG) experiment testbed, Yucca Flat, Nevada National Security Site

Geologic framework models (GFMs) provide a methodology for integrating geology into, and thus geologically informing, other modeling and simulation activities. GFMs provide a three-dimensional (3-D), geology-based digital framework for developing and parametrizing meshes and evaluating model and simulation results. This report describes a 3-D GFM constructed for the Source Physics Experiment Phase II Dry Alluvium Geology test series located in Yucca Flat at the Nevada National Security Site. The geology in the Yucca Flat region is complex and diverse, which creates challenges to modeling seismic wave propagation from SPE tests. The Yucca Flat GFM helps address these challenges by providing the 3-D distribution of relevant geologic features and physical properties necessary to more effectively model seismic wave propagation. The GFM includes 7 model layers and 48 faults that cut and offset the layers. An appendix is included that provides quantitative data on physical properties for each model layer.

58 GEOSCIENCES↗

Ion Channel Laser Based on Direct Laser Acceleration of a Shaped Beam Driver (Final Technical Report)

The main focus of our research was to understand the synergies between an electron bunch and a laser pulse when the two co-propagate through the plasma. We have clearly demonstrated using theoretical and computational modeling that the propagation distance of both the bunch and the laser pulse could be extended. The ability of the combined bunch/pulse system to extend the propagation distance and the size of the plasma bubble enables highly-efficient sources of relativistic electrons. As pointed out by the Plasma Decadal Study, such electron sources can be used for generating extremely bright x-rays for a variety of applications can serve as probes and diagnostics for other plasma experiments: high energy density (HED) sciences, inertial confinement fusion (ICF), and potentially Fusion Materials and Technology (FM&T). Development of advanced diagnostics of plasma-based accelerators is yet another key area identified by recent reports. Even broader security and medical applications of compact accelerator-based radiation sources, such as very high energy electron (VHEE) sources for FLASH radiobiology, have been identified by a recent multi-agency panel. A number of key technical issues were considered and successfully resolved during the course of the grant. Those include: beam loading effect of both the driver and witness bunches on the wake, laser channeling by the bunch, and the direct laser acceleration (DLA) of the driver bunch by the laser pulse. To demonstrate clear synergy, we were able to ascertain that the total energy gain of a witness bunch in the wake of the combined bunch/laser complex is large than the sum of the energy gains in the wake of the laser pulse alone, and the driver bunch alone. We have also demonstrated that not only the energy gain is improved, but the energy spread is not sacrificed. To carry out these simulations, we have to develop a range of in-house computational tools, including a fully-3D code WAND-PIC.

43 PARTICLE ACCELERATORS↗

Temporal Explosion Source Processes of Declared Nuclear Tests in the Democratic People’s Republic of Korea

In this work we highlight a preliminary temporal source analysis of the six declared Democratic People's Republic of Korea (DPRK) nuclear tests. We use regional seismic data to estimate relative source time functions (RSTFs) via iterative time-domain deconvolution (Ammon, 2006; Pippin, 2022) of vertical-component ground motions recorded within 2000 km of the source region. Since RSTFs are ideally independent of site and propagation effects, their amplitude spectrum is equivalent to the source spectral ratio, but they also retain phase information. We compare observed RSTFs (in the time and frequency domains) with synthetic RSTFs derived from the Mueller & Murphy (1971) explosion source model. The resolution of these time functions varies, however, we generally obtain high-quality results within the limitations of the recording broadband instrumentation. The results indicate that this method effectively preserves source time-history information that can be used for temporal analysis of remote nuclear explosions. This preliminary analysis is intended to assess the viability of using time-domain deconvolution methods for extracting temporal source information.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Quantum Machine Learning with SQUID

In this work we present the Scaled QUantum IDentifier (SQUID), an open-source framework for exploring hybrid Quantum-Classical algorithms for classification problems. The classical infrastructure is based on PyTorch and we provide a standardized design to implement a variety of quantum models with the capability of back-propagation for efficient training. We present the structure of our framework and provide examples of using SQUID in a standard binary classification problem from the popular MNIST dataset. In particular, we highlight the implications for scalability for gradient-based optimization of quantum models on the choice of output for variational quantum models.

97 MATHEMATICS AND COMPUTING↗