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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1982

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.↗

NASA NTRS1983

The cosmic ray abundances of the platinum-lead elements as measured on HEAO-3

The relative abundances of elements in the charge ranges of Z = 75-79 (platinum) and Z = 80-83 (lead) should be a sensitive indication of the contributions of the r- and s-processes in nucleosynthesis. Data from the HEAO 3 Heavy Nuclei Experiment are used to establish abundances, relative to iron, of these elements in the cosmic radiation, as well as the ratio of 'secondary' elements, in the Z = 62-74 range, to the primary lead-platinum elements. These results appear to suggest that either the source abundances are deficient in s-process elements or that they are not organized solely by first ionization potential. In addition, present propagation models can adequately represent the relative abundances of primary and secondary elements.

Fixsen, D. J.↗

NASA NTRS1993

Are cosmic rays effective for ionization of the solar nebula?

In this paper, we argue that the effectiveness of cosmic rays to ionize the bulk of the nebular gas may be further impaired by the influence of the magnetic field on the propagation of cosmic rays. When cosmic rays enter the nebular disk they ionize the gas and make the dynamo generation of magnetic fields possible. However, once magnetic fields are embedded in the nebular gas, the upcoming cosmic rays can no longer penetrate directly into the nebular disk because they start to interact with the magnetic field and lose their energy before propagating significantly toward the midplane. That, in turn, undercuts the ionization source within the bulk of the gas stopping the dynamo action. Nebular dynamo models ignored this back reaction of magnetic fields on cosmic rays. We calculate this back reaction effect, but for the sake of mathematical simplicity, we ignore the effect of magnetic field weakening due to diminishing ionization by cosmic rays.

Dolginov, A. Z.↗

DOE OSTI2021

Programmatic Advantages of Linear Equivalent Seismic Models

Underground explosions nonlinearly deform the surrounding earth material and can interact with the free surface to produce spall. However, at typical seismological observation distances the seismic wavefield can be accurately modeled using linear approximations. Although nonlinear algorithms can accurately simulate very near field ground motions, they are computationally expensive and potentially unnecessary for far field wave simulations. Conversely, linearized seismic wave propagation codes are orders of magnitude faster computationally and can accurately simulate the wavefield out to typical observational distances. Thus, devising a means of approximating a nonlinear source in terms of a linear equivalent source would be advantageous both for scenario modeling and for interpretation of seismic source models that are based on linear, far-field approximations. This allows fast linear seismic modeling that still incorporates many features of the nonlinear source mechanics built into the simulation results so that one can have many of the advantages of both types of simulations without the computational cost of the nonlinear computation. In this report we first show the computational advantage of using linear equivalent models, and then discuss how the near-source (within the nonlinear wavefield regime) environment affects linear source equivalents and how well we can fit seismic wavefields derived from nonlinear sources.

58 GEOSCIENCES↗

NASA NTRS2024

Exploring the ENSO modulation of the QBO periods with GISS E2.2 models

Observational studies have shown that the El Niño–Southern Oscillation (ENSO) exerts an influence on the Quasi-Biennial Oscillation (QBO). The downward propagation of the QBO tends to speed up and slow down during El Niño and La Niña, respectively. Recent results from general circulation models have indicated that the ENSO modulation of the QBO requires a relatively high horizontal resolution, and that it does not show up in the climate models with parameterized but temporally constant gravity wave sources. Here, we demonstrate that the NASA Goddard Institute for Space Studies (GISS) E2.2 models can capture the observed ENSO modulation of the QBO period with a horizontal resolution of 2 ∘ latitude by 2.5 ∘ longitude but with its gravity wave sources being parameterized interactively. This is because El Niño events lead to more vigorous gravity wave sources generating more absolute momentum fluxes over the equatorial belt, as well as less filtering of these waves into the tropical lower stratosphere through a weakening of the Walker circulation. Various components of the ENSO system, such as the sea surface temperatures, the convective activities, and the Walker circulation, are intimately involved in the generation and propagation of parameterized gravity waves, through which ENSO modulates the QBO period in GISS E2.2 models.

El Niño–Southern Oscillation↗

DOE OSTI2025

Logical error rates for the surface code under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions

With fault-tolerant quantum computing (FTQC) on the horizon, it is critical to understand sources of logical errors in plausible hardware implementations of quantum error-correcting codes. Detailed error modeling of computational instructions on particular FTQC architectures will enable the better prediction of error propagation in FT-encoded quantum circuits while revealing where greater attention is needed in hardware design. In this work, we consider logical error rates for the surface code implemented on a hypothetical grid-based trapped-ion quantum charge-coupled device architecture. Specifically, we construct logical channels for the idling surface code and examine its diamond error under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions. We include the coherent dephasing noise that is known to accumulate during physical qubit idling and transport in these systems, determining idling and transport durations using the time-resolved output of an open-source trapped-ion surface code compiler. To estimate expectation values of logical Pauli observables following hardware circuits containing non-Clifford sources of noise, we utilize a Monte Carlo technique to sample from an underlying quasiprobability distribution of Clifford circuits that we independently simulate in a phase-sensitive fashion. We verify error suppression up to code distance 𝑑 = 11 at coherent dephasing rates near and below those of current-generation trapped-ion quantum computers and find that logical error rates align with those of analogous fully stochastic simulations in this regime. Exploring higher dephasing rates at 𝑑 = 3−5, we find evidence for growing coherent rotations about all three logical Pauli axes, increased diagonal logical error process matrix elements relative to those of stochastic simulations, and a reduced dephasing rate threshold. Overall, our work paves a way toward realistic hardware emulation of small fault-tolerant quantum processes, e.g., members of an FTQC instruction set.

Quantum benchmarking↗

NASA NTRS2012

The Scaling of Broadband Shock-Associated Noise with Increasing Temperature

A physical explanation for the saturation of broadband shock-associated noise (BBSAN) intensity with increasing jet stagnation temperature has eluded investigators. An explanation is proposed for this phenomenon with the use of an acoustic analogy. For this purpose the acoustic analogy of Morris and Miller is examined. To isolate the relevant physics, the scaling of BBSAN at the peak intensity level at the sideline ( = 90 degrees) observer location is examined. Scaling terms are isolated from the acoustic analogy and the result is compared using a convergent nozzle with the experiments of Bridges and Brown and using a convergent-divergent nozzle with the experiments of Kuo, McLaughlin, and Morris at four nozzle pressure ratios in increments of total temperature ratios from one to four. The equivalent source within the framework of the acoustic analogy for BBSAN is based on local field quantities at shock wave shear layer interactions. The equivalent source combined with accurate calculations of the propagation of sound through the jet shear layer, using an adjoint vector Green s function solver of the linearized Euler equations, allows for predictions that retain the scaling with respect to stagnation pressure and allows for the accurate saturation of BBSAN with increasing stagnation temperature. This is a minor change to the source model relative to the previously developed models. The full development of the scaling term is shown. The sources and vector Green s function solver are informed by steady Reynolds-Averaged Navier-Stokes solutions. These solutions are examined as a function of stagnation temperature at the first shock wave shear layer interaction. It is discovered that saturation of BBSAN with increasing jet stagnation temperature occurs due to a balance between the amplification of the sound propagation through the shear layer and the source term scaling.A physical explanation for the saturation of broadband shock-associated noise (BBSAN) intensity with increasing jet stagnation temperature has eluded investigators. An explanation is proposed for this phenomenon with the use of an acoustic analogy. For this purpose the acoustic analogy of Morris and Miller is examined. To isolate the relevant physics, the scaling of BBSAN at the peak intensity level at the sideline psi = 90 degrees) observer location is examined. Scaling terms are isolated from the acoustic analogy and the result is compared using a convergent nozzle with the experiments of Bridges and Brown and using a convergent-divergent nozzle with the experiments of Kuo, McLaughlin, and Morris at four nozzle pressure ratios in increments of total temperature ratios from one to four. The equivalent source within the framework of the acoustic analogy for BBSAN is based on local field quantities at shock wave shear layer interactions. The equivalent source combined with accurate calculations of the propagation of sound through the jet shear layer, using an adjoint vector Green s function solver of the linearized Euler equations, allows for predictions that retain the scaling with respect to stagnation pressure and allows for the accurate saturation of BBSAN with increasing stagnation temperature. This is a minor change to the source model relative to the previously developed models. The full development of the scaling term is shown. The sources and vector Green s function solver are informed by steady Reynolds-Averaged Navier-Stokes solutions. These solutions are examined as a function of stagnation temperature at the first shock wave shear layer interaction. It is discovered that saturation of BBSAN with increasing jet stagnation temperature occurs due to a balance between the amplification of the sound propagation through the shear layer and the source term scaling.

Miller, Steven A.↗

NASA NTRS2004

Diffuse Galactic Continuum Gamma Rays. A Model Compatible with EGRET Data and Cosmic-ray Measurements

We present a study of the compatibility of some current models of the diffuse Galactic continuum gamma-rays with EGRET data. A set of regions sampling the whole sky is chosen to provide a comprehensive range of tests. The range of EGRET data used is extended to 100 GeV. The models are computed with our GALPROP cosmic-ray propagation and gamma-ray production code. We confirm that the "conventional model" based on the locally observed electron and nucleon spectra is inadequate, for all sky regions. A conventional model plus hard sources in the inner Galaxy is also inadequate, since this cannot explain the GeV excess away from the Galactic plane. Models with a hard electron injection spectrum are inconsistent with the local spectrum even considering the expected fluctuations; they are also inconsistent with the EGRET data above 10 GeV. We present a new model which fits the spectrum in all sky regions adequately. Secondary antiproton data were used to fix the Galactic average proton spectrum, while the electron spectrum is adjusted using the spectrum of diffuse emission it- self. The derived electron and proton spectra are compatible with those measured locally considering fluctuations due to energy losses, propagation, or possibly de- tails of Galactic structure. This model requires a much less dramatic variation in the electron spectrum than models with a hard electron injection spectrum, and moreover it fits the y-ray spectrum better and to the highest EGRET energies. It gives a good representation of the latitude distribution of the y-ray emission from the plane to the poles, and of the longitude distribution. We show that secondary positrons and electrons make an essential contribution to Galactic diffuse y-ray emission.

Strong, Andrew W.↗

NASA NTRS1985

Parametrized energy spectrum of cosmic-ray protons with kinetic energies down to 1 GeV

A new estimation of the interstellar proton spectrum is made in which the source term of primary protons is taken from shock acceleration theory and the cosmic ray propagation calculation is based on a proposed nonuniform galactic disk model.

Tan, L. C.↗

NASA NTRS1984

The plasma wave environment of an auroral arc. II - ULF waves on an auroral arc boundary

On March 9, 1978, a sounding rocket launched from Poker Flat, Alaska, at 2200 LT, made a four-component measurement of a 0.5 Hz hydromagnetic wave as the payload crossed the poleward boundary of a quiet homogeneous auroral arc. An energy flux of about 10 to the -6th W/sq m was observed propagating upward with a left-handed polarization within the arc, and a flux 6 times greater was observed propagating downward with a right-handed polarization on the arc boundary. The waves were identified as shear mode Alfven waves. Various models for the source of the free energy are discussed with the conclusion that the most likely production mechanism was either the electromagnetic or electrostatic Kelvin-Helmholtz instability.

Gelpi, C. G.↗

NASA NTRS2007

Experimental Measurements of Two-dimensional Planar Propagating Edge Flames

The study of edge flames has received increased attention in recent years. This work reports the results of a recent study into two-dimensional, planar, propagating edge flames that are remote from solid surfaces (called here, free-layer flames, as opposed to layered flames along floors or ceilings). They represent an ideal case of a flame propagating down a flammable plume, or through a flammable layer in microgravity. The results were generated using a new apparatus in which a thin stream of gaseous fuel is injected into a low-speed laminar wind tunnel thereby forming a flammable layer along the centerline. An airfoil-shaped fuel dispenser downstream of the duct inlet issues ethane from a slot in the trailing edge. The air and ethane mix due to mass diffusion while flowing up towards the duct exit, forming a flammable layer with a steep lateral fuel concentration gradient and smaller axial fuel concentration gradient. We characterized the flow and fuel concentration fields in the duct using hot wire anemometer scans, flow visualization using smoke traces, and non-reacting, numerical modeling using COSMOSFloWorks. In the experiment, a hot wire near the exit ignites the ethane air layer, with the flame propagating downwards towards the fuel source. Reported here are tests with the air inlet velocity of 25 cm/s and ethane flows of 967-1299 sccm, which gave conditions ranging from lean to rich along the centerline. In these conditions the flame spreads at a constant rate faster than the laminar burning rate for a premixed ethane air mixture. The flame spread rate increases with increasing transverse fuel gradient (obtained by increasing the fuel flow rate), but appears to reach a maximum. The flow field shows little effect due to the flame approach near the igniter, but shows significant effect, including flow reversal, well ahead of the flame as it approaches the airfoil fuel source.

Villa-Gonzalez, Marcos↗

NASA NTRS2017

Modeling Solar Energetic Particle Events Using ENLIL Heliosphere Simulations

Solar energetic particle (SEP) event modeling (SEPMOD) has gained renewed attention in part because of the availability of a decade of multipoint measurements from STEREO (Solar TErrestrial RElations Observatory) and L1 (Lagrangian point 1) spacecraft at 1 AU (Astronomical Unit). These observations are coupled with improving simulations of the geometry and strength of heliospheric shocks obtained by using coronagraph images to send erupted material into realistic solar wind backgrounds. The STEREO and ACE (Aerosol, Cloud systems, ocean Ecosystems) measurements in particular have highlighted the sometimes surprisingly widespread nature of SEP events. It is thus an opportune time for testing SEP models, which typically focus on protons approximately 1-100 megaelectronvolts, toward both physical insight to these observations and potentially useful space radiation environment forecasting tools. Some approaches emphasize the concept of particle acceleration and propagation from close to the Sun, while others emphasize the local field line connection to a traveling, evolving shock source. Among the latter is the previously introduced SEPMOD treatment, based on the widely accessible and well-exercised WSA-ENLIL (Wang-Sheeley-Arge-ENLIL)-cone model. SEPMOD produces SEP proton time profiles at any location within the ENLIL domain. Here we demonstrate a SEPMOD version that accommodates multiple, concurrent shock sources occurring over periods of several weeks. The results illustrate the importance of considering longer-duration time periods and multiple CME (Coronal Mass Ejection) contributions in analyzing, modeling, and forecasting SEP events.

Luhmann, J. G.↗

DOE OSTI2021

Bayesian learning of orthogonal embeddings for multi-fidelity Gaussian Processes

Uncertainty propagation in complex engineering systems often poses significant computational challenges related to modeling and quantifying probability distributions of model outputs, as those emerge as the result of various sources of uncertainty that are inherent in the system under investigation. Gaussian Processes regression (GPs) is a robust meta-modeling technique that allows for fast model prediction and exploration of response surfaces. Multi-fidelity variations of GPs further leverage information from cheap and low fidelity model simulations in order to improve their predictive performance on the high fidelity model. In order to cope with the high volume of data required to train GPs in high dimensional design spaces, a common practice is to introduce latent design variables that are typically projections of the original input space to a lower dimensional subspace, and therefore substitute the problem of learning the initial high dimensional mapping, with that of training a GP on a low dimensional space. Here in this paper, we present a Bayesian approach to identify optimal transformations that map the input points to low dimensional latent variables. The \projection" mapping consists of an orthonormal matrix that is considered a priori unknown and needs to be inferred jointly with the GP parameters, conditioned on the available training data. The proposed Bayesian inference scheme relies on a two-step iterative algorithm that samples from the marginal posteriors of the GP parameters and the projection matrix respectively, both using Markov Chain Monte Carlo (MCMC) sampling. In order to take into account the orthogonality constraints imposed on the orthonormal projection matrix, a Geodesic Monte Carlo sampling algorithm is employed, that is suitable for exploiting probability measures on manifolds. We extend the proposed framework to multi-fidelity models using GPs including the scenarios of training multiple outputs together. We validate our framework on three synthetic problems with a known lower-dimensional subspace. The benefits of our proposed framework, are illustrated on the computationally challenging aerodynamic optimization of a last-stage blade for an industrial gas turbine, where we study the effect of an 85-dimensional shape parameterization of a three-dimensional airfoil on two output quantities of interest, specifically on the aerodynamic efficiency and the degree of reaction

42 ENGINEERING↗

NASA NTRS1988

Evaluation of the communications impact of a low power arcjet thruster

The interaction of a 1 kW arcjet thruster plume with a communications signal is evaluated. A two-parameter, source flow equation has been used to represent the far flow field distribution of the arcjet plume in a realistic spacecraft configuration. Modelling the plume as a plasma slab, the interaction of the plume with a 4 GHz communications signal is then evaluated in terms of signal attenuation and phase shift between transmitting and receiving antennas. Except for propagation paths which pass very near the arcjet source, the impacts to transmission appear to be negligible. The dominant signal loss mechanism is refraction of the beam rather than absorption losses due to collisions. However, significant reflection of the signal at the sharp vacuum-plasma boundary may also occur for propagation paths which pass near the source.

Carney, Lynnette M.↗

NASA NTRS1988

Evaluation of the communications impact of a low power arcjet thruster

Carney, Lynnette M.↗

NASA NTRS2006

Using CO2:CO Correlations to Improve Inverse Analyses of Carbon Fluxes

Observed correlations between atmospheric concentrations of CO2 and CO represent potentially powerful information for improving CO2 surface flux estimates through coupled CO2-CO inverse analyses. We explore the value of these correlations in improving estimates of regional CO2 fluxes in east Asia by using aircraft observations of CO2 and CO from the TRACE-P campaign over the NW Pacific in March 2001. Our inverse model uses regional CO2 and CO surface fluxes as the state vector, separating biospheric and combustion contributions to CO2. CO2-CO error correlation coefficients are included in the inversion as off-diagonal entries in the a priori and observation error covariance matrices. We derive error correlations in a priori combustion source estimates of CO2 and CO by propagating error estimates of fuel consumption rates and emission factors. However, we find that these correlations are weak because CO source uncertainties are mostly determined by emission factors. Observed correlations between atmospheric CO2 and CO concentrations imply corresponding error correlations in the chemical transport model used as the forward model for the inversion. These error correlations in excess of 0.7, as derived from the TRACE-P data, enable a coupled CO2-CO inversion to achieve significant improvement over a CO2-only inversion for quantifying regional fluxes of CO2.

Palmer, Paul I.↗

DOE OSTI2021

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗

DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗