Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source term calculation”

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 505 records · Page 28

'Flux conservation' by a Schwarzschild gravitational lens

It is demonstrated, by an explicit calculation, that a single, isolated Schwarzschild gravitational lens leads to a vanishing net total amplification on its own, without any 'help' from other lenses. It is shown that, in the limit of small M/R (where M is the mass of the lens in geometrical units and R is the distance from the source to the observers), the positive net total amplification, obtained from the classical approximation of the amplification for small impact angle lensing, is balanced by a negative net total amplification from the first-order term in M/R in an approximation of the amplification for large impact angle lensing derived here. The relation of these results to the common interpretation of average flux conservation by gravitational lenses is discussed as due to the collective action of many lenses in a cosmological framework.

Shulami, I.↗

A Numerical Experiment on the Role of Surface Shear Stress in the Generation of Sound

The sound generated due to a localized flow over an infinite flat surface is considered. It is known that the unsteady surface pressure, while appearing in a formal solution to the Lighthill equation, does not constitute a source of sound but rather represents the effect of image quadrupoles. The question of whether a similar surface shear stress term constitutes a true source of dipole sound is less settled. Some have boldly assumed it is a true source while others have argued that, like the surface pressure, it depends on the sound field (via an acoustic boundary layer) and is therefore not a true source. A numerical experiment based on the viscous, compressible Navier-Stokes equations was undertaken to investigate the issue. A small region of a wall was oscillated tangentially. The directly computed sound field was found to to agree with an acoustic analogy based calculation which regards the surface shear as an acoustically compact dipole source of sound.

Shariff, Karim↗

The velocity dispersion of the giant molecular clouds - A viscous origin

The energy source is proposed and the details of the acceleration mechanism for the random motion of the giant molecular clouds (GMCs) in the Galaxy are discussed. Specifically, it is proposed that gravitational scattering of the massive clouds off each other in the differentially rotating Galactic disk constitutes an effective gravitational viscosity, which causes an increase in the random kinetic energy of the GMCs at the expense of their ordered, rotational kinetic energy in the Galactic disk. The energy input due to this gravitational viscosity is calculated, and the other terms contributing to the energy balance of the GMCs are computed. These results are used to determine the steady state cloud velocity dispersion and its dependence on the cloud mass and radius and on the galactocentric radial distance.

Jog, Chanda J.↗

Comparison of Aerosol Single Scattering Albedos Derived By Diverse Techniques in Two North Atlantic Experiments

Aerosol single scattering albedo w (the ratio of scattering to extinction) is important in determining aerosol climatic effects, in explaining relationships between calculated and measured radiative fluxes, and in retrieving aerosol optical depths from satellite radiances. Recently, two experiments in the North Atlantic region, TARFOX and ACE-2, determined aerosol w by a variety of techniques. The techniques included fitting of calculated to measured fluxes; retrievals of w from skylight radiances; best fits of complex refractive index to profiles of backscatter, extinction, and size distribution; and in situ measurements of scattering and absorption at the surface and aloft. Both TARFOX and ACE-2 found a fairly wide range of values for w at midvisible wavelengths, with 0.85 less than wmidvis less than 0.99 for the marine aerosol impacted by continental pollution. Frequency distributions of w could usually be approximated by lognormals in wmax-w, with some occurrence of bimodality, suggesting the influence of different aerosol sources or processing. In both TARFOX and ACE-2, closure tests between measured and calculated radiative fluxes yielded best-fit values of wmidvis of 0.90+/-0.04 for the polluted boundary layer. Although these results have the virtue of describing the column aerosol unperturbed by sampling, they are subject to questions about representativeness and possible artifacts (e.g., unknown gas absorption). The other techniques gave larger values for wmidvis for the polluted boundary layer, with a typical result of wmidvis = 0.95+/-0.04, Current uncertainties in vv are large in terms of climate effects. More tests are needed of the consistency among different methods and of humidification effects on w.

Russell, P. B.↗

Some atmospheric scattering considerations relevant to BATSE: A model calculation

The orbiting Burst and Transient Source Experiement (BATSE) will locate gamma ray burst sources by analysis of the relative numbers of photons coming directly from a source and entering its prescribed array of detectors. In order to accurately locate burst sources it is thus necessary to identify and correct for any counts contributed by events other than direct entry by a mainstream photon. An effort is described which estimates the photon numbers which might be scattered into the BATSE detectors from interactions with the Earth atmosphere. A model was developed which yielded analytical expressions for single-scatter photon contributions in terms of source and satellite locations.

Young, John H.↗

Extension of SCALE/Sampler’s sensitivity analysis

Nuclear data are a major source of uncertainties in reactor physics calculations. The propagation of nuclear data uncertainties to important system responses is instrumental when determining appropriate safety margins in reactor safety analyses. It is also important to understand the major contributors to the observed uncertainties to make recommendations for further measurements and evaluations and aid in the understanding of the studied system. The SCALE code system allows for nuclear data uncertainty analysis based on the random sampling approach as implemented in SCALE’s Sampler sequence. Sampler was recently extended by a sensitivity analysis in terms of the calculation of two correlation-based sensitivity indices. This analysis allows for the identification of the top contributing nuclear reactions to any analyzed output uncertainty. This paper presents the sensitivity indices, along with their interpretation and limitations. It demonstrates the application in an eigenvalue and decay heat analysis for a boiling water reactor fuel assembly.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Increasing phosphorus loss despite widespread concentration decline in US rivers

The loss of phosphorous (P) from the land to aquatic systems has polluted waters and threatened food production worldwide. Systematic trend analysis of P, a nonrenewable resource, has been challenging, primarily due to sparse and inconsistent historical data. Here, we leveraged intensive hydrometeorological data and the recent renaissance of deep learning approaches to fill data gaps and reconstruct temporal trends. We trained a multitask long short-term memory model for total P (TP) using data from 430 rivers across the contiguous United States (CONUS). Trend analysis of reconstructed daily records (1980–2019) shows widespread decline in concentrations, with declining, increasing, and insignificantly changing trends in 60%, 28%, and 12% of the rivers, respectively. Concentrations in urban rivers have declined the most despite rising urban population in the past decades; concentrations in agricultural rivers however have mostly increased, suggesting not-as-effective controls of nonpoint sources in agriculture lands compared to point sources in cities. TP loss, calculated as fluxes by multiplying concentration and discharge, however exhibited an overall increasing rate of 6.5% per decade at the CONUS scale over the past 40 y, largely due to increasing river discharge. Results highlight the challenge of reducing TP loss that is complicated by changing river discharge in a warming climate.

Science & Technology - Other Topics↗

Interface extinction and subsurface peaking of the radiation pattern of a line source

The radiation pattern of a line source lying along the plane interface of two dielectric half-spaces is calculated. It is found that the pattern at the interface has a null (interface extinction); that the pattern in the upper half-space, whose index of refraction is taken to be less than that of the lower half-space, has a single lobe with a maximum normal to the interface; and that the pattern in the lower half-space (subsurface region) has two maxima (peaks) straddling symmetrically a minimum. Interpretation of these results in terms of ray optics, Oseen's extinction theorem, and the Cerenkov effect are given.

Engheta, N.↗

Radionuclide Screening Analysis and Transport Parameters for Pahute Mesa Detonations, Nevada National Security Site

A “screening analysis” is implemented at Pahute Mesa (PM), Nevada National Security Site (NNSS), to determine subsets of the 43 radionuclides listed in the radionuclide inventory of Finnegan et al. (2016) that are of relevance or potential relevance to the hydrologic source term (HST) for assessment of radionuclide transport in groundwater. Consideration is added to how levels of contamination are defined and whether any other radionuclides not in the inventory are relevant to the HST. A model is developed to estimate a range of possible source concentrations in groundwater that account for uncertainty in partitioning into melt glass and sorption into surrounding rock of the exchange zone for the 82 PM underground nuclear tests detonated in vertical shafts. The model also accounts for the varied hydrogeochemical settings. Transport parameters needed for the screening model calculations are developed from databases for hydrogeologic units, chemistry, mineralogy, fracture spacing, fracture aperture, fracture openness, matrix porosity, bulk density, saturation, and alteration. Of key importance is consideration of diagenetic zonation of the mineralogy related to hydrothermal alteration. The screening analysis compares model results with available groundwater radiochemistry data. Consistency between screening model results and available data helps validate the model for application to all 43 radionuclides, most of which have no measurements or only non-detect measurements of concentration in groundwater. The screening analysis determines that ten radionuclides – tritium, Sr-90, I-129, Cs-137, U-232, U-233, U-234, U-238, and the total of Pu-239 and Pu-240 are relevant to the HST. Determination of relevance is primarily based on data and/or model results indicating source concentration exceeding a maximum contaminant level (MCL) out to 100 years since the time of detonation. U-238 is relevant to assessment of contamination relative to MCLs for uranium and gross alpha particle activity. Five other radionuclides – C-14, Cl-36, Tc-99, Np-237, and Pu-238 – are determined potentially relevant to the HST based on a 0.1 MCL threshold. The screening analysis recommends additional attention to the natural daughter radionuclides of the uranium series (U-238) and thorium series (Th-232) decay chains, particularly Ra-226 and Ra-228, that were not included the inventory.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Numerical Simulation of Turbulent MHD Flows Using an Iterative PNS Algorithm

A new parabolized Navier-Stokes (PNS) algorithm has been developed to efficiently compute magnetohydrodynamic (MHD) flows in the low magnetic Reynolds number regime. In this regime, the electrical conductivity is low and the induced magnetic field is negligible compared to the applied magnetic field. The MHD effects are modeled by introducing source terms into the PNS equation which can then be solved in a very efficient manner. To account for upstream (elliptic) effects, the flowfields are computed using multiple streamwise sweeps with an iterated PNS algorithm. Turbulence has been included by modifying the Baldwin-Lomax turbulence model to account for MHD effects. The new algorithm has been used to compute both laminar and turbulent, supersonic, MHD flows over flat plates and supersonic viscous flows in a rectangular MHD accelerator. The present results are in excellent agreement with previous complete Navier-Stokes calculations.

Kato, Hiromasa↗

Observational and theoretical evidence in support of a significant in-situ photochemical source of tropospheric ozone

The latitudinal and seasonal variation of ozone in the troposphere is discussed. Of particular interest is the asymmetrical behavior of this gas with respect to the two hemispheres. These asymmetries, when coupled with a diagnostic photochemical model of the troposphere, lends support to the view that ozone cannot be viewed as an inert tracer of stratospheric origin. In the calculations it is noted that it is likely that the budgets of carbon monoxide and tropospheric ozone may be quite dependent on each other and the calculations are discussed in light of the uncertainty which currently exists about representative global tropospheric background concentrations of the nitrogen oxides. In addition, the seasonal variation of excess (C-14)O2 (a stratospheric tracer) is examined and compared with the seasonal ozone variation during the same period of observations at the same location and altitudes. The distinct maxima for ozone found during the summer in the lower troposphere are not present for the (C-14)O2 data. This finding likewise suggests that photochemical processes taking place in the troposphere are an important source term for tropospheric ozone.

Fishman, J.↗

Imaging the sun in hard x rays using Fourier telescopes

For several years, solar flares have been observed with a variety of instruments confirming that tremendous amounts of energy are locally stored in the solar magnetic field and then rapidly released during the life of the flare. In concert with observations, theorists have attempted to describe the means by which these energetic events occur and evolve. Two competing theories have emerged and have stood the test of time. One theory describes the flare in terms of nonthermal, electron beam injection into a thick target while the other uses a thermal approach. Both theories provide results which are reasonably consistent with current observations; but to date, none have been able to provide conclusive evidence as to the validity of either model. Imaging on short time scales (1 s) and/or small size scales (1 arc s) should give definitive answers to these questions. In order to test whether a realistic telescope can indeed discriminate between models, we construct model sources based upon the thermal and the nonthermal models and calculate the emission as a function of time and energy in the range from 10 to 100 keV. In addition, we construct model telescopes representing both the spatial modulation collimator (SMC) and the rotating modulation collimator (RMC) techniques of observation using random photon counting statistics. With these two types of telescopes we numerically simulate the instrument response to the above two model flares to see if there are distinct x-ray signatures which may be discernable. We find that theoretical descriptions of the primary models of solar flares do indeed predict different hard x-ray signatures for 1 sec time scales and at 1-5 arc sec spatial resolution. However, these distinguishing signatures can best be observed early in the impulsive phase and from a position perpendicular to the plane of the loop. Furthermore, we find that Fourier telescopes with reasonable and currently attainable design characteristics can image these signatures and that the same sensitive areas and short temporal integration times relative to source evolution (i.e., 1 s), the RMC and the SMC will both provide about the same performance.

Campbell, J. W.↗

High charge state carbon and oxygen ions in Earth's equatorial quasi-trapping region

Observations of energetic (1.5 - 300 keV/e) medium-to-high charge state (+3 less than or equal to Q less than or equal to +7) solar wind origin C and O ions made in the quasi-trapping region (QTR) of Earth's magnetosphere are compared to ion trajectories calculated in model equatorial magnetospheric magnetic and electric fields. These comparisons indicate that solar wind ions entering the QTR on the nightside as an energetic component of the plasma sheet exit the region on the dayside, experiencing little or no charge exchange on the way. Measurements made by the CHarge Energy Mass (CHEM) ion spectrometer on board the Active Magnetospheric Particle Tracer Explorer/Charge Composition Explorer (AMPTE/CCE) spacecraft at 7 less than L less than 9 from September 1984 to January 1989 are the source of the new results contained herein: quantitative long-term determination of number densities, average energies, energy spectra, local time distributions, and their variation with geomagnetic disturbance level as indexed by Kp. Solar wind primaries (ions with charge states unchanged) and their secondaries (ions with generally lower charge states produced from primaries in the magnetosphere via charge exchange)are observed throughout the QTR and have distinctly different local time variations that persist over the entire 4-year analysis interval. During Kp larger than or equal to 3 deg intervals, primary ion (e.g., O(+6)) densities exhibit a pronounced predawn maximum with average energy minimum and a broad near-local-noon density minimum with average energy maximum. Secondary ion (e.g., O(+5)) densities do not have an identifiable predawn peak, rather they have a broad dayside maximum peaked in local morning and a nightside minimum. During Kp less than or equal to 2(-) intervals, primary ion density peaks are less intense, broader in local time extent, and centered near midnight, while secondary ion density local time variations diminish. The long-time-interval baseline helps to refine and extend previous observations; for example, we show that ionospheric contribution to O(+3)) is negligible. Through comparison with model ion trajectories, we interpret the lack of pronounced secondary ion density peaks colocated with the primary density peaks to indicate that: (1) negligible charge exchange occurs at L greater than 7, that is, solar wind secondaries are produced at L less than 7, and (2) solar wind secondaries do not form a significant portion of the plasma sheet population injected into the QTR. We conclude that little of the energetic solar wind secondary ion population is recirculated through the magnetosphere.

Christon, S. P.↗

Measuring the 511 keV emission in the direction of 1E1740.7-2942 with BATSE

Observations of the 511 keV emission in the direction of 1E 1740.7-2942 (1E) using the earth burst and transient source experiment (BATSE) onboard the Compton Gamma Ray Observatory (CGRO), are presented. The CGRO phase 1 average spectrum of 1E is calculated using a method which assumes that a given source spectrum is the sum of the flux coming directly from the object and the contribution from the surrounding diffuse emission. The 1E light curve is calculated in the 40 to 150 keV range. It presents a constant flux excess of 70 mCrab in comparison with observations from the SIGMA gamma ray telescope onboard the GRANAT observatory. By removing this contribution, the 1E spectral transition from the low state to the high standard state observed by SIGMA is confirmed, and it is shown that the 511 keV flux is independent of the 1E long term evolution from low state to high standard state. It is concluded that the 511 keV emission of (4.2 +/- 1.3) x 140(exp -4) photons/sq cm s observed in the direction of 1E is mainly diffuse and spatially extended.

Wallyn, P.↗

Fatigue analysis of the Spallation Neutron Source 2 MW target design

Upgrades to the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory are under way with two major projects. The first project is the Proton Power Upgrade (PPU), which will double the power capacity of the SNS accelerator system to enable the operation of a future second neutron source target station—the second project. PPU also will increase power to the SNS first target station – which currently operates at 1.4 MW with 1.0 GeV protons – to 2 MW with 1.3 GeV protons. Final design of the PPU 2 MW target module is finished, including all necessary design and analysis calculations. The fatigue evaluation described herein is an assessment of the fitness of the PPU 2 MW target mercury vessel to resist failure from cyclic loading from the pulsed 60 Hz proton beam and thermal changes from disrupted operation. In addition to the reliability predictions, the method, detail, and depth of this fatigue evaluation are superior to those of past target design assessments. Furthermore, the thermal fatigue life is predicted to exceed 3125 thermal cycles per 1250 h of operation for all load cases. The design also meets the (relative) fatigue design goal for combined thermal and pulse loading of greater than 0.5 times the minimum fatigue life calculated for a jet-flow target design operating at 1.4 MW—the most robust design to date in terms of resistance to fatigue failure. From the onset of the PPU, it has been known that pulse loading requires effective application of helium gas injection into the target mercury to reduce beam pressure loading. For the PPU 2 MW target, none of the pulse load cases requires a maximum strain reduction from gas injection of more than 50% (for base material) to meet the fatigue design goal. This level of reduction is achievable from gas injection, according to historical strain measurements from 1.4 MW target operation. Overall, the PPU 2 MW target design is predicted to have superior resistance to fatigue failure compared with past target designs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modeling impurity concentrations in liquid argon detectors

Impurities in noble liquid detectors used for neutrino and dark matter experiments can significantly impact the quality of data. We present an experimentally verified model for describing the dynamics of impurity distributions in liquid argon (LAr) detectors. The model considers sources, sinks, and transport of impurities within and between the gas and liquid argon phases. Measurements of oxygen concentrations in a 20-L LAr multi-purpose test stand are compared to calculations made with this model to show that an accurate description of the concentrations under various operational conditions can be obtained. A result of this analysis is a determination of Henry’s coefficient for oxygen in LAr. Additionally, these calculations also show that some processes have small effects on the impurity dynamics and excluding them yields a solution as a sum of two exponential terms. This solution provides a simple way to extract Henry’s coefficient with negligible approximation error. It is applied to the data and the Henry’s coefficient for oxygen in LAr is obtained as 0.84 +0.09 -0.05 , consistent with literature results. Finally, based on the analysis of the data with the model, we further suggest that, for a large liquid argon detector, barriers to flow (“baffles”) installed in the gas phase to restrict flow can help reduce the ultimate impurity concentration in the LAr.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗

Calculated and measured performance of a 'near ideal' locally reacting duct liner in grazing incidence

The calculated acoustic response of a nearly ideal (locally reacting and linear) duct liner has been compared with its measured performance. The ceramic honeycomb liner structure was such that the local reaction assumption was well satisfied over the frequency range of interest (0.5-3.0 kHz). In addition, the high repeatability of the measured impedance values over a large SPL range verified the assumption of linearity. A finite element algorithm was used to compute the duct liner transmission properties in terms of the liner's normal incidence impedance and the duct termination impedance using a plane wave source. Very good agreement between theory and experiment was observed in general, but required accurate normal impedance values in the high attenuation frequency bands of the liner's response.

Parrott, T. L.↗