Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed sources”

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 145 records · Page 8

Applying Machine Learning and Bayesian Inference to Identify and Locate Moving Anthropogenic Sources Using Distributed Acoustic Sensing Data

Distributed acoustic sensing (DAS) systems, which use existing telecommunication fibers, offer high‐resolution capabilities ideal for recording anthropogenic sources. However, the complexity of urban environments and the large amount of data recorded by DAS require automated methods to efficiently detect and categorize anthropogenic sources. Here, we evaluate how well three machine learning models (k‐nearest neighbor [k‐NN], convolutional neural networks, and recurrent‐convolutional neural networks) can identify various anthropogenic sources recorded by DAS. Our findings reveal that both k‐NN and neural network methods perform well in high signal‐to‐noise ratio (SNR) settings. However, their accuracy decreases at SNRs <4. We also use Kalman filtering, a form of Bayesian inference, on backprojected locations of these sources to recover locations that generally fall within standard smartphone Global Positioning System errors. By combining machine learning and Kalman filter results, we calculate a multidimensional model of moving anthropogenic sources. These results demonstrate the potential of DAS data in urban seismology for accurately identifying and locating such sources. Depending on the research objectives, these sources can be further studied or filtered out to improve the quality of seismic data for earthquake studies. Such methods provide a valuable tool for urban seismology and seismic hazard analysis.

Luckie, Thomas William [Sandia National Laboratori↗

Theory on acoustic sources

A theory is described for the radiation emission emission from acoustic multipole sources. The sources can be stationary or moving at speeds including supersonic and experience stationary or moving disturbances. The effect of finite source distributions and disturbances is investigated as well as the manner in which they interact. Distinction is made between source distributions that responsed as a function of time and those that respond as a function of space.

Wright, S. E.↗

Numerical convergence and validation of the DIMP inverse particle transport model

The data integration with modeled predictions (DIMP) model is a promising inverse radiation transport method for solving the special nuclear material (SNM) holdup problem. Unlike previous methods, DIMP is a completely passive nondestructive assay technique that requires no initial assumptions regarding the source distribution or active measurement time. DIMP predicts the most probable source location and distribution through Bayesian inference and quasi-Newtonian optimization of predicted detector re-sponses (using the adjoint transport solution) with measured responses. DIMP performs well with for-ward hemispherical collimation and unshielded measurements, but several considerations are required when using narrow-view collimated detectors. DIMP converged well to the correct source distribution as the number of synthetic responses increased. DIMP also performed well for the first experimental validation exercise after applying a collimation factor, and sufficiently reducing the source search vol-ume's extent to prevent the optimizer from getting stuck in local minima. DIMP's simple point detector response function (DRF) is being improved to address coplanar false positive/negative responses, and an angular DRF is being considered for integration with the next version of DIMP to account for highly collimated responses. Overall, DIMP shows promise for solving the SNM holdup inverse problem, especially once an improved optimization algorithm is implemented.

42 ENGINEERING↗

GEANT4 investigation of 40 K and 222 Rn interference effects on actinide detection sensitivities of CEMRC BEGe lung detectors

Internal dosimetry measurements of trans-uranic (TRU) radionuclides within the human body can be accomplished by γ-ray spectroscopy using high-purity germanium detectors. While measurements are typically conducted inside shielded chambers in order to minimize the background detector signal, the presence of radioactive 40 K within the body cannot be removed, which produces a high-energy γ ray at 1461 keV. The 40 K, in addition to airborne 222 Rn, contributes to the radiation field within the chamber that increases to the minimum detectable activity that can be achieved in practice. The objective of the present work is to estimate the contributions from 40 K internal to the human body and airborne 222 Rn to the background signal in a lung counter germanium detector array using GEANT4 simulations, with a focus on TRU isotope energy region (17–59 keV). The investigations consisted of three main simulations: 1. Simulating the γ-ray flux at the geometric center of the CEMRC counting chamber arising from a 25 nCi source of 40 K from (a) a point source, and (b) a source distributed throughout a BOMAB phantom. 2. Simulating the response of the CEMRC BEGe detectors arising from a 25 nCi source of 40 K from (a) a point source, and (b) a source distributed throughout a BOMAB phantom. 3. Simulating the response of the CEMRC BEGe detectors to the environmental background in the counting chamber, consisting of airborne 222 Rn (0–4 pCi/L) and 40 K (0–160 nCi) internal to the human body. The results are used to investigate the 40 K background contribution to a routine spectrum recorded by the lung detectors at the CEMRC Internal Dosimetry laboratory over a typical half hour counting interval used in the dosimetry of radiological workers. The minimum detectable activity was calculated for various 222 Rn and 40 K activities and compared to experimentally determined values. The contribution of cosmic rays was not included in this investigation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

The Diversity of Refractory Organic Material in Comets

Organic matter exists in comets (most notably in 81P/Wild 2 [Stardust], 67P/Churyomov- Gerasimenko (67P/C-G) [Rosetta], chondritic porous IDPS, and UCAMMs) and in primitive carbonaceous chondrites that likely retain some chemistry that reflects an origin in the prenatal cold molecular cloud (Alexander+2017). Heavy isotopic enrichments, 15N/14N and possibly D/H, signify preserved molecular cloud organics. In the cold outer disk, if grains are lofted above the disk mid-plane then organics likely experience significant UV processing (Ciesla+2012). In remote sensing of comet comae, organics in the dust are considered refractory or semirefractory. Semi-refractory organics have limited comae lifetimes and produce distributed sources of molecules (H2CO and CO). Rosetta's close passes of 67P/C-G's nucleus (10-15 km) reveals a distributed source of glycine, methyl amine and ethylamine (Altwegg+2016). Cometary samples and primitive meteorites have two types of organic matter: (1) acid-insoluble organic matter (IOM), which is a macromolecular polymer with a mixture of aromatic and aliphatic moieties, and (2) labile, soluble organics, which includes the amino acids, such as glycine (Stardust, Elsila+2009). Meteoritic IOM is robust, withstanding experimental temperatures of 1200 K (Dobrica+2011, Cody+2008). Nanoglobules are a type of IOM; they have a distinct physical structure, but often share the same chemistry as the other IOM from the same meteorite. Moderate-sized PAHs (20 C-atoms) are detected in Stardust samples (Clemett+2010). Refractory organic IOM is ubiquitous yet has a great diversity of abundances between cometary samples. IOM is in primitive chondrites, 67P/C-G (Rosetta), 81P/Wild 2 (Stardust), 1P/Halley, 26P/Grigg-Skjellerup, UCAMMs, anhydrous IDPs, and in chondritic porous IDPs (CP IDPs) and larger cluster IDPs (e.g., Fray+2016, Fomenkova+94, Busemann+ 09, Dobrica+2011, Dobrica+2012). 81P's refractory organic matter is of two types (De Gregorio+2011): nanoglobules of highly aromatic refractory organic matter and polyaromatic carbonyl-containing organic matter, which is similar to IOM in primitive meteorites and IDPs. Fray+2017 estimate that 50% of carbon in 67P/C-G is in IOM. 67P/C-G's organics appear to lack the soluble organic matter, aliphatic carbon, amino acids, and PAHs (Fray+2016). Other notable aspects of the diversity in IOM in cometary samples are the ranges of atomic ratios of N/C, O/C, and H/C, and the range of isotopic enrichments of 15N/14N and D/H. Aqueous and thermal processing on asteroids changes the balance of soluble to insoluble organics, and may be important for diversifying the range of OM delivered to Earth.

Wooden, Diane↗

Recent advances in antenna synthesis.

A technique for the synthesis of planar aperture antennas has been developed. It applies equally well to continuous or array type apertures. Patterns with shaped main beams and complex sidelobe structures can be synthesized. The method involves applying corrections iteratively and converging toward the desired pattern. The iteration process terminates when the synthesized pattern approximates the desired pattern to within a specified tolerance. A computer program has been written to perform all calculations. The source distribution which is required to produce the synthesized pattern is provided as output. The method is such that patterns requiring complex source distributions will not be synthesized (such as superdirective patterns).

Stutzman, W. L.↗

Waste Characterization Methods Based on Gamma-ray Imaging Spectrometers - 19454

Characterization of LLW, ILW, U, NORM and TENORM is critical to the nuclear industry's bottom line. Accurate measurements enable the disposal of waste at the lowest possible cost. In the past the most accurate technology available involved using HPGe detectors and extensive system modeling to estimate the total source term given what the user knows about the geometry of the source and any shielding between the source and the detectors. Now there are a variety of imaging detectors that are available commercially, including CZT and HPGe detectors that achieve energy resolution ≤1% FWHM at 662 keV. These devices couple the traditional high resolution spectroscopic analysis with next generation gamma-ray imaging information to more accurately characterize the radiation source term. Uncertainties in the source location will result in significant error in the characterization of an object's radioactivity. It is helpful if the object can be spun, such as a barrel, so that any asymmetry in the source distribution is averaged out from the detector's point of view. Even in the case of a spinning barrel it is important to determine if the radiation is concentrated near the center of the barrel or distributed near the edge. The closer one measures to a barrel, which usually dictates the time required to acquire data, the more these spatial variations in the source distribution will affect the calculated activity of each isotope. This paper will focus on showing the use of CZT imaging spectrometers to characterize complicated radiation source terms. The radiation images are used to determine where the source term should be defined in a radiation transport model that calculates the expected number of photopeak counts for each isotope detected. This is used to calculate the activity of the source. The transport model is used to determine the attenuation of the radiation source term through distance and shielding material, but the detection efficiency is based solely on radiation measurement data. These measurements are taken with a robotic source positioning device that allows for detailed directional efficiency calibration. It is a fully automated system where a detector is loaded into the device and all of the efficiency curves are captured as a function of direction by creating new measurements for each location the robot positions the source. These measurements are tagged with the true source location and the efficiency curve is calculated for each measurement. This is used to create a response function matrix as a function of energy and incident direction which is used to convert the direction dependent flux on the detector head determined from the simulation model to the expected number of counts per unit activity of each isotope of interest. The results of these source term reconstructions in both the spatial and spectral domain will be shown. The impact of improper specification of the source term will be shown for a few standard waste characterization scenarios to demonstrate the importance of proper source location determination. Finally, the underlying calibration method and some sample efficiency curves as a function of direction will be shown to demonstrate the importance of calibrating out directional response functions in imaging devices (which often have a strong forward sensitivity bias). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Self-consistent current sheet structures in the quiet-time magnetotail

The structure of the quiet-time magnetotail is studied using a test particle simulation. Vlasov equilibria are obtained in the regime where v(D) = E(y) c/B(z) is much less than the ion thermal velocity and are self-consistent in that the current and magnetic field satisfy Ampere's law. Force balance between the plasma and magnetic field is satisfied everywhere. The global structure of the current sheet is found to be critically dependent on the source distribution function. The pressure tensor is nondiagonal in the current sheet with anisotropic temperature. A kinetic mechanism is proposed whereby changes in the source distribution results in a thinning of the current sheet.

Holland, Daniel L.↗

A matrix-inversion method for gamma-source mapping from gamma-count data

In a previous paper it was proposed that a simple matrix inversion method could be used to extract source distributions from gamma-count maps, using simple models to calculate the response matrix. The method was tested using numerically generated count maps. In the present work a 100 kBq Co{sup 60} source has been placed on a gridded surface and the count rate measured using a NaI scintillation detector. The resulting map of gamma counts was used as input to the matrix inversion procedure and the source position recovered. A multi-source array was simulated by superposition of several single-source count maps and the source distribution was again recovered using matrix inversion. The measurements were performed for several detector heights. The effects of uncertainties in source-detector distances on the matrix inversion method are also examined. The results from this work give confidence in the application of the method to practical applications, such as the segregation of highly active objects amongst fuel-element debris. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Saturn Neutron Exosphere as Source for Inner and Innermost Radiation Belts

Energetic proton and electron measurements by the ongoing Cassini orbiter mission are expanding our knowledge of the highest energy components of the Saturn magnetosphere in the inner radiation belt region after the initial discoveries of these belts by the Pioneer 11 and Voyager 2 missions. Saturn has a neutron exosphere that extends throughout the magnetosphere from the cosmic ray albedo neutron source at the planetary main rings and atmosphere. The neutrons emitted from these sources at energies respectively above 4 and 8 eV escape the Saturn system, while those at lower energies are gravitationally bound. The neutrons undergo beta decay in average times of about 1000 seconds to provide distributed sources of protons and electrons throughout Saturn's magnetosphere with highest injection rates close to the Saturn and ring sources. The competing radiation belt source for energetic electrons is rapid inward diffusion and acceleration of electrons from the middle magnetosphere and beyond. Minimal losses during diffusive transport across the moon orbits, e.g. of Mimas and Enceladus, and local time asymmetries in electron intensity, suggest that drift resonance effects preferentially boost the diffusion rates of electrons from both sources. Energy dependences of longitudinal gradient-curvature drift speeds relative to the icy moons are likely responsible for hemispheric differences (e.g., Mimas, Tethys) in composition and thermal properties as at least partly produced by radiolytic processes. A continuing mystery is the similar radial profiles of lower energy (<10 MeV) protons in the inner belt region. Either the source of these lower energy protons is also neutron decay, but perhaps alternatively from atmospheric albedo, or else all protons from diverse distributed sources are similarly affected by losses at the moon' orbits, e.g. because the proton diffusion rates are extremely low. Enceladus cryovolcanism, and radiolytic processing elsewhere on the icy moon and ring surfaces, are additional sources of protons via ionization and charge exchange from breakup of water molecules. But one must then account somehow for local acceleration to the observed keV-MeV energies, since moon sweeping and E-ring absorption would remove protons diffusing inward from the middle magnetosphere. Although the main rings block further inward diffusion from the inner radiation belts, the exospheric neutron-decay source, combined with much slower diffusion of protons relative to electrons, may produce an innermost radiation belt in the gap between the upper atmosphere and the D-ring. This innermost belt will first be explored in-situ during the final proximal orbits of the Cassini mission.

Cooper, John↗

Jet Noise Diagnostics Supporting Statistical Noise Prediction Methods

The primary focus of my presentation is the development of the jet noise prediction code JeNo with most examples coming from the experimental work that drove the theoretical development and validation. JeNo is a statistical jet noise prediction code, based upon the Lilley acoustic analogy. Our approach uses time-average 2-D or 3-D mean and turbulent statistics of the flow as input. The output is source distributions and spectral directivity. NASA has been investing in development of statistical jet noise prediction tools because these seem to fit the middle ground that allows enough flexibility and fidelity for jet noise source diagnostics while having reasonable computational requirements. These tools rely on Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) solutions as input for computing far-field spectral directivity using an acoustic analogy. There are many ways acoustic analogies can be created, each with a series of assumptions and models, many often taken unknowingly. And the resulting prediction can be easily reverse-engineered by altering the models contained within. However, only an approach which is mathematically sound, with assumptions validated and modeled quantities checked against direct measurement will give consistently correct answers. Many quantities are modeled in acoustic analogies precisely because they have been impossible to measure or calculate, making this requirement a difficult task. The NASA team has spent considerable effort identifying all the assumptions and models used to take the Navier-Stokes equations to the point of a statistical calculation via an acoustic analogy very similar to that proposed by Lilley. Assumptions have been identified and experiments have been developed to test these assumptions. In some cases this has resulted in assumptions being changed. Beginning with the CFD used as input to the acoustic analogy, models for turbulence closure used in RANS CFD codes have been explored and compared against measurements of mean and rms velocity statistics over a range of jet speeds and temperatures. Models for flow parameters used in the acoustic analogy, most notably the space-time correlations of velocity, have been compared against direct measurements, and modified to better fit the observed data. These measurements have been extremely challenging for hot, high speed jets, and represent a sizeable investment in instrumentation development. As an intermediate check that the analysis is predicting the physics intended, phased arrays have been employed to measure source distributions for a wide range of jet cases. And finally, careful far-field spectral directivity measurements have been taken for final validation of the prediction code. Examples of each of these experimental efforts will be presented. The main result of these efforts is a noise prediction code, named JeNo, which is in middevelopment. JeNo is able to consistently predict spectral directivity, including aft angle directivity, for subsonic cold jets of most geometries. Current development on JeNo is focused on extending its capability to hot jets, requiring inclusion of a previously neglected second source associated with thermal fluctuations. A secondary result of the intensive experimentation is the archiving of various flow statistics applicable to other acoustic analogies and to development of time-resolved prediction methods. These will be of lasting value as we look ahead at future challenges to the aeroacoustic experimentalist.

Bridges, James E.↗

Phased-Array Measurements of Single Flow Hot Jets

A 16 microphone phased-array system has been successfully applied to measure jet noise source distributions. In this study, a round convergent nozzle was tested at various hot and cold flow conditions: acoustic Mach numbers are between 0.35 and 1.6 and static temperature ratios are varied from cold to 2.7. The classical beamforming method was applied on narrowband frequencies. From the measured source distributions locations of peak strength were tracked and found to be very consistent between adjacent narrowband frequencies. In low speed heated and unheated jets, the peak source locations vary smoothly from the nozzle exit to downstream as the frequency is decreased. When the static temperature ratio was kept constant, the peak source position moved downstream with increasing acoustic Mach number for the Strouhal numbers smaller than about 1.5. It was also noted that the peak source locations of low frequencies occur farther downstream than the end of potential core.

Bridges, James↗

Methods for Detector Placement and Analysis of Criticality Accident Alarm Systems

Determining the optimum placement to minimize the number of detectors for a criticality accident alarm system (CAAS) in a large manufacturing facility is a complex problem. There is typically a target for the number of detectors that can be used over a given zone of the facility. A study to optimize detector placement typically begins with some initial guess at the placement of the detectors and is followed by either predictive calculations of accidents at specific locations or adjoint calculations based on preferred detector locations. Within an area of a facility, there may be a large number of potential criticality accident sites. For any given placement of the detectors, the list of accident sites can be reduced to a smaller number of locations at which accidents may be difficult for detectors to detect. Developing the initial detector placement and determining the list of difficult accident locations are both based on the practitioner's experience. Simulations following fission particles released from an accident location are called 'forward calculations.' These calculations can be used to answer the question 'where would an alarm be triggered?' by an accident at a specified location. Conversely, 'adjoint calculations' start at a detector site using the detector response function as a source and essentially run in reverse. These calculations can be used to answer the question 'where would an accident be detected?' by a specified detector location. If the number of accidents, P, is much less than the number of detectors, Q, then forward simulations may be more convenient and less time-consuming. If Q is large or the detectors are not placed yet, then a mesh tally of dose observed by a detector at any location must be computed over the entire zone. If Q is much less than P, then adjoint calculations may be more efficient. Adjoint calculations employing a mesh tally can be even more advantageous because they do not rely on a list of specific difficult-to-detect accident sites, which may not have included every possible accident location. Analog calculations (no biasing) simply follow particles naturally. For sparse buildings and line-of-sight calculations, analog Monte Carlo (MC) may be adequate. For buildings with internal walls or large amounts of heavy equipment (dense geometry), variance reduction may be required. Calculations employing the CADIS method use a deterministic calculation to create an importance map and a matching biased source distribution that optimize the final MC to quickly calculate one specific tally. Calculations employing the FW-CADIS method use two deterministic calculations (one forward and one adjoint) to create an importance map and a matching biased source distribution that are designed to make the MC calculate a mesh tally with more uniform uncertainties in both high-dose and low-dose areas. Depending on the geometry of the problem, the number of detectors, and the number of accident sites, different approaches to CAAS placement studies can be taken. These are summarized in Table I. SCALE 6.1 contains the MAVRIC sequence, which can be used to perform any of the forward-based approaches outlined in Table I. For analog calculations, MAVRIC simply calls the Monaco MC code. For CADIS and FW-CADIS, MAVRIC uses the Denovo discrete ordinates (SN) deterministic code to generate the importance map and biased source used by Monaco. An adjoint capability is currently being added to Monaco and should be available in the next release of SCALE. An adjoint-based approach could be performed with Denovo alone - although fine meshes, large amounts of memory, and long computation times may be required to obtain accurate solutions. Coarse-mesh SN simulations could be employed for adjoint-based scoping studies until the adjoint capability in Monaco is complete. CAAS placement studies, especially those dealing with mesh tallies, require some extra utilities to aid in the analysis. Detectors must receive a minimum dose rate in order to alarm; therefore, a simple yes/no plot could be more useful to the analyst than a standard dose rate contour plot. Alarm systems that require several detectors to alarm simultaneously for a given accident would need to combine several yes/no plots in order to show accident sites where multiple detectors would be triggered. This could require a plot of 'areas of coverage' that is mapped over the building geometry. Several new utilities (which will be part of SCALE 6.2) were created to help with CAAS analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Hexavalent Chromium Sources and Distribution in California Groundwater

Groundwater resources in California represent a confluence of high-risk factors for hexavalent chromium contamination as a result of industrial activities, natural geology, and, potentially, land use. Here, we examine state-wide links in California between groundwater Cr(VI) concentrations and chemicals that provide signatures for source attribution. In environmental monitoring wells, Cr(VI) had the highest co-occurrence and also clustered with 1,4-dioxane and several chlorinated hydrocarbons indicative of the metal plating industry. Additionally, hotspots of Cr(VI) co-occurring with bromoform result from volatile organic compound remediation using in situ chemical oxidation that inadvertently oxidizes naturally occurring Cr(III). In groundwater supply wells, which are typically free of industrial inputs, Cr(VI) correlates with dichlorodiphenyldichloroethylene (DDE), vanadium, and ammonia and clusters with nitrate and dissolved oxygen, suggesting potential links between agricultural activities and Cr(VI). Specific controls on Cr(VI) vary substantially by region: from the metal plating industry around Los Angeles and the San Francisco Bay areas to natural redox conditions along flow paths in the Mojave Desert and to correlations with agricultural practices in the Central Valley of California. Here, while industrial uses of Cr lead to the most acute cases of groundwater Cr(VI) contamination, oxidation of naturally occurring Cr affects a larger area, more wells, and a greater number of people throughout California.

54 ENVIRONMENTAL SCIENCES↗

An adaptive fast multipole accelerated Poisson solver for complex geometries

Here, we present a fast, direct and adaptive Poisson solver for complex two-dimensional geometries based on potential theory and fast multipole acceleration. More precisely, the solver relies on the standard decomposition of the solution as the sum of a volume integral to account for the source distribution and a layer potential to enforce the desired boundary condition. The volume integral is computed by applying the FMM on a square box that encloses the domain of interest. For the sake of efficiency and convergence acceleration, we first extend the source distribution (the right-hand side in the Poisson equation) to the enclosing box as a $C$ 0 function using a fast, boundary integral-based method. We demonstrate on multiply connected domains with irregular boundaries that this continuous extension leads to high accuracy without excessive adaptive refinement near the boundary and, as a result, to an extremely efficient “black box” fast solver.

97 MATHEMATICS AND COMPUTING↗

Transfer function analysis of thermospheric perturbations

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Mayr, H. G.↗

Jet-Surface Interaction Noise Models Comparison to TMP17 and Next Steps

A set of empirical jet-surface interaction noise models, developed for single-stream round nozzles exhausting over a simple surface in a static ambient, are evaluated for use in more realistic applications that include multi-stream nozzle systems with chevrons and dihedral surface geometries. The simple-single-stream models have several advantages when used in system-level noise studies: they are robust, they are quickly computed, and they are generally applicable to a wide range of configurations. However, these models rely on a fixed noise source distribution that does not account for the presence of chevrons at the nozzle exit. The effect of this omission on the prediction accuracy is shown by this comparison to experimental data. This motivates discussion of future work to develop a more robust modeling framework that can account for generalized source distributions and more complex airframe geometries. Re: Top-Mounted Propulsion (TMP17) Test Plans.

Propulsion Noise↗

Comparison of ADVANTG Simulations to Simple Shielding Measurements - Paper 36

Monte Carlo (MC) radiation transport simulations, like those performed using Monte Carlo N-Particle (MCNP), can be highly accurate but require time to reduce the statistical uncertainties in the calculated detector response(s). Many different 'variance reduction' techniques devised over the years increase the efficiency of MC calculations, but these methods can be more art than science. Many of these methods are problem specific or require careful tuning of parameters for good performance. One of the best variance reduction methods, in terms of both reliable performance increases and ease of use, is the use of weight windows. Weight windows are derived from an importance map in space and energy that allows the transport code to focus more simulation time on particles that have a higher likelihood of being important to the final detector response. Less time is spent simulating particles that have little chance of making a contribution to the detector response. MCNP already includes the abilities to generate and use weight windows. Unfortunately, generating weight windows can take a very long time, perhaps longer than the original problem. Oak Ridge National Laboratory has led the development of hybrid methods using a discrete ordinates adjoint calculation to quickly develop an importance map (weight windows) for use in a Monte Carlo code. The Consistent Adjoint Driven Importance Sampling (CADIS) method uses an adjoint calculation to develop both an importance map and a biased source distribution that work together. This method is reasonably automatic, is easy to use, and has been implemented in the Automated Variance reduction Generator (ADVANTG) package that works with MCNP5. ADVANTG reads an MCNP input file and a small amount of extra user information and then runs an adjoint deterministic calculation using Denovo. The importance map and biased source distribution derived from the CADIS method can then be incorporated back into the original MCNP input file to enhance the figure-of-merit (FOM) of the calculation. The FOM is a measure of the efficiency of the calculation using the amount of run time and corresponding level of statistical uncertainty in the detector result. ADVANTG-enhanced MCNP input files have been shown to increase the FOM of many types of problems, obtaining the same level of statistical uncertainty in a shorter time compared to standard MCNP. ADVANTG can also generate variance reduction parameters (importance map and biased source) using the Forward-Weighted CADIS (FW-CADIS) method for optimizing mesh tallies. The CADIS method creates variance reduction parameters that optimize a tally's total quantity (flux, dose rate, interaction rate, etc.) and will not generally increase performance for all parts of a segmented tally. To optimize a segmented tally, the FW-CADIS method uses two discrete ordinates calculations to make variance reduction parameters that help compute a segmented tally with uniform relative uncertainties in each bin. This report will demonstrate that using ADVANTG for neutron and photon radiation shielding problems will (1) preserve the standard MCNP results and (2) obtain those results either with less statistical uncertainty or in shorter times than the standard MCNP calculation. The four demonstration problems in this report are based on simple neutron or photon transmission measurements through different materials. All of the problems here used the CADIS method and still saw improvement across the flux tallies segmented by particle energy. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗