Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source term estimation”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

High speed turboprop aeroacoustic study (single rotation). Volume 1: Model development

A frequency-domain noncompact-source theory for the steady loading and volume-displacement (thickness) noise of high speed propellers has been developed and programmed. Both near field and far field effects have been considered. The code utilizes blade surface pressure distributions obtained from three-dimensional nonlinear aerodynamic flow field analysis programs as input for evaluating the steady loading noise. Simplified mathematical models of the velocity fields induced at the propeller disk by nearby wing and fuselage surfaces and by angle-of-attack operation have been developed to provide estimates of the unsteady loading imposed on the propeller by these potential field type interactions. These unsteady blade loadings have been coupled to a chordwise compact propeller unsteady loading noise model to provide predictions of unsteady loading noise caused by these installation effects. Finally, an analysis to estimate the corrections to be applied to the free-field noise predictions in order to arrive at the measurable fuselage sound pressure levels has been formulated and programmed. This analysis considers the effects of fuselage surface reflection and diffraction together with surface boundary layer refraction. The steady loading and thickness model and the unsteady loading model have been verified using NASA-supplied data for the SR-2 and SR-3 model propfans. In addition, the steady loading and thickness model has been compared with data from the SR-6 model propfan. These theoretical models have been employed in the evaluation of the SR-7 powered Gulfstream aircraft in terms of noise characteristics at representative takeoff, cruise, and approach operating conditions. In all cases, agreement between theory and experiment is encouraging.

Whitfield, C. E.↗
NASA NTRS1998

Envelope structure on 700 AU scales and the molecular outflows of low-mass young stellar objects

Aperture synthesis observations of HCO+ J = 1-0, 13CO 1-0, and C18O 1-0 obtained with the Owens Valley Millimeter Array are used to probe the small-scale (5" approximately 700 AU) structure of the molecular envelopes of a well-defined sample of nine embedded low-mass young stellar objects in Taurus. The interferometer results can be understood in terms of: (1) a core of radius approximately or less than 1000 AU surrounding the central star, possibly flattened and rotating; (2) condensations scattered throughout the envelope that may be left over from the inhomogeneous structure of the original cloud core or that may have grown during collapse; and (3) material within the outflow or along the walls of the outflow cavity. Masses of the central cores are 0.001-0.1 M (solar), and agree well with dust continuum measurements. Averaged over the central 20" (3000 AU) region, an HCO+ abundance of 4 x 10(-8) is inferred, with a spread of a factor of 3 between the different sources. Reanalysis of previously presented single-dish data yields an HCO+ abundance of (5.0 +/- 1.7) x 10(-9), which may indicate an average increase by a factor of a few on the smaller scales sampled by the interferometer. Part of this apparent abundance variation could be explained by contributions from extended cloud emission to the single-dish C18O lines, and uncertainties in the assumed excitation temperatures and opacities. The properties of the molecular envelopes and outflows are further investigated through single-dish observations of 12CO J = 6-5, 4-3, and 3-2, 13CO 6-5 and 3-2, and C18O 3-2 and 2-1, obtained with the James Clerk Maxwell and IRAM 30 m telescopes, along with the Caltech Submillimeter Observatory. Ratios of the mid-J CO lines are used to estimate the excitation temperature, with values of 25-80 K derived for the gas near line centre. The outflow wings show a similar range, although Tex is enhanced by a factor of 2-3 in at least two sources. In contrast to the well-studied L1551 IRS 5 outflow, which extends over 10' (0.4 pc), seven of the remaining eight sources are found to drive 12CO 3-2 outflows over < or = 1' (0.04 pc); only L1527 IRS has a well-developed outflow of some 3'(0.12 pc). Estimates are obtained for the outflow kinetic luminosity, Lkin, and the flow momentum rate, FCO, applying corrections for line opacity and source inclination. The flow force FCO correlates with the envelope mass and with the 2.7 mm flux of the circumstellar disk. Only a weak correlation is seen with Lbol, while none is found with the relative age of the object as measured by integral Tmb(HCO+ 3-2)dV/Lbol. These trends support the hypothesis that outflows are driven by accretion through a disk, with a global mass infall rate determined by the mass and density of the envelope. The association of compact HCO+ emission with the walls of the outflow cavities indicates that outflows in turn influence the appearance of the envelopes. It is not yet clear, however, whether they are actively involved in sweeping up envelope material, or merely provide a low-opacity pathway for heating radiation to reach into the envelope.

NASA Discipline Exobiology↗
NASA NTRS2010

Premonsoon Aerosol Characterization and Radiative Effects Over the Indo-Gangetic Plains: Implications for Regional Climate Warming

The Himalayas have a profound effect on the South Asian climate and the regional hydrological cycle, as it forms a barrier for the strong monsoon winds and serves as an elevated heat source, thus controlling the onset and distribution of precipitation during the Indian summer monsoon. Recent studies have suggested that radiative heating by absorbing aerosols, such as dust and black carbon over the Indo-Gangetic Plains (IGP) and slopes of the Himalayas, may significantly accelerate the seasonal warming of the Hindu Kush-Himalayas-Tibetan Plateau (HKHT) and influence the subsequent evolution of the summer monsoon. This paper presents a detailed characterization of aerosols over the IGP and their radiative effects during the premonsoon season (April-May-June) when dust transport constitutes the bulk of the regional aerosol loading, using ground radiometric and spaceborne observations. During the dust-laden period, there is a strong response of surface shortwave flux to aerosol absorption indicated by the diurnally averaged forcing efficiency of -70 W/sq m per unit optical depth. The simulated aerosol single-scattering albedo, constrained by surface flux and aerosol measurements, is estimated to be 0.89+/- 0.01 (at approx.550 nm) with diurnal mean surface and top-of-atmosphere forcing values ranging from -11 to -79.8 W/sq m and +1.4 to +12 W/sq m, respectively, for the premonsoon period. The model-simulated solar heating rate profile peaks in the lower troposphere with enhanced heating penetrating into the middle troposphere (5-6 km), caused by vertically extended aerosols over the IGP with peak altitude of approx.5 km as indicated by spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization observations. On a long-term climate scale, our analysis, on the basis of microwave satellite measurements of tropospheric temperatures from 1979 to 2007, indicates accelerated annual mean warming rates found over the Himalayan-Hindu Kush region (0.21 C/decade+/-0.08 C/decade) and underscores the potential role of enhanced aerosol solar absorption in the maximum warming localized over the western Himalayas (0.26 C/decade f 0.09 C/decade) that significantly exceed the entire HKHT and global warming rates. We believe the accelerated warming rates reported here are critical to both the South Asian summer monsoon and hydro-glaciological resource variability in the Himalayan-Hindu Kush snowpack and therefore to the densely populated downstream regions.

Gautam, Ritesh↗
NASA NTRS2007

An Experimental System for a Global Flood Prediction: From Satellite Precipitation Data to a Flood Inundation Map

Floods impact more people globally than any other type of natural disaster. It has been established by experience that the most effective means to reduce the property damage and life loss caused by floods is the development of flood early warning systems. However, advances for such a system have been constrained by the difficulty in estimating rainfall continuously over space (catchment-. national-, continental-. or even global-scale areas) and time (hourly to daily). Particularly, insufficient in situ data, long delay in data transmission and absence of real-time data sharing agreements in many trans-boundary basins hamper the development of a real-time system at the regional to global scale. In many countries around the world, particularly in the tropics where rainfall and flooding co-exist in abundance, satellite-based precipitation estimation may be the best source of rainfall data for those data scarce (ungauged) areas and trans-boundary basins. Satellite remote sensing data acquired and processed in real time can now provide the space-time information on rainfall fluxes needed to monitor severe flood events around the world. This can be achieved by integrating the satellite-derived forcing data with hydrological models, which can be parameterized by a tailored geospatial database. An example that is a key to this progress is NASA's contribution to the Tropical Rainfall Measuring Mission (TRMM), launched in November 1997. Hence, in an effort to evolve toward a more hydrologically-relevant flood alert system, this talk articulates a module-structured framework for quasi-global flood potential naming, that is 'up to date' with the state of the art on satellite rainfall estimation and the improved geospatial datasets. The system is modular in design with the flexibility that permits changes in the model structure and in the choice of components. Four major components included in the system are: 1) multi-satellite precipitation estimation; 2) characterization of land surface including digital elevation from NASA SRTM, topography-derived hydrologic parameters such as flood direction. flow accumulation, basin, and river network etc.; 3) spatially distributed hydrological models to infiltrate rainfall and route overland runoff; and 4) an implementation interface to relay thc input data to the models and display the flood inundation results to the users and decision-makers. Early results appear reasonable in terms of location and frequency of events. Case studies of this experimental system are evaluated with surface runoff data and other river monitoring systems. such as Dartmouth Flood Observatory's "Surface Water Watch" array of river reaches that are measured daily via other satellite remote sensing data. A major outcome of this progress will be the availability of a global overview of flood alerts that should consequently improve the performance of Decision Support System. We expect these developments in utilization of satellite remote sensing technology to offer a practical solution to the challenge of building a cost-effective early warning system for data scarce and under-developed areas.

Adler, Robert↗
DOE OSTI2023

Surface Complexation/Ion Exchange Hybrid Model for Radionuclide Sorption to Clay Minerals (M4SF-23LL010301062)

This progress report (Level 4 Milestone Number M4SF-23LL010301062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite International Collaborations Activity Number SF-23LL01030106. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. The focus of this milestone is the establishment of international collaborations for surface complexation modeling and the associated impacts of unlocking larger, community-based datasets. More specifically, we are developing a database framework for Spent Fuel and Waste and Science Technology (SFWST) that is aligned with the Helmholtz Zentrum Dresden Rossendorf (HZDR) sorption database development group in support of the database needs of the SFWST program. In our FY22 effort, we described a detailed analysis of U(VI) sorption to quartz through both traditional surface complexation modeling and through a hybrid ML framework. In FY23, effort was placed on publication of these results and expansion of the LLNL surface complexation and ion exchange database (L-SCIE) in order to assess mineral-based radionuclide retardation under a wider variety of geochemical conditions (e.g., ionic strength, varying electrolyte compositions). Efforts were initiated to expand L-SCIE to include radionuclide surface complexation and ion exchange to clays that are relevant to subsurface geochemical processes occurring at nuclear waste repositories. In particular, a large source of sorption data for clays resides at the Paul Scherrer Institute (PSI) (work primarily by Bradbury and Baeyens) and we initiated discussions on how to retrieve those data and apply FAIR principles to those datasets. In addition to L-SCIE development, two hybrid models that incorporate AI/ML were investigated and compared to discern the most promising approaches for accurate and precise estimations of radionuclide retardation. Key considerations for future model development include (1) the ability to reduce computational burden on determining retardation coefficients for PA and (2) the ability to quantify and predict radionuclide-mineral partitioning at a more efficient, rapid pace due to automated workflows. Upon the careful consideration of the most effective modeling approaches, we are identifying ways to implement these approaches into PA. Ultimately, the data science-based workflows will provide a major incentive for other institutions to adopt a FAIR-formatted, interoperable database. LLNL will play a key role in disseminating sorption data and acting as good data stewards by updating the database in a consistent format and assessing the quality of the newly assimilated data in an organized fashion. To this end, all data and workflows are open access and made available on the LLNL Seaborg research website (https://seaborg.llnl.gov/resources/geochemical-databases-modeling-codes).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS2020

What Can We Learn from One Billion Ground System Log Messages?

Shortage of log-based data in a ground system they have traditionally been the under achievers in a satellite ground system. This is due to several factors: Once log messages scroll out of view on the TTC event console window they are soon forgotten. Application and system log files are scattered across directories within a system, across a multitude of servers, and across one or more databases making access cumbersome. Typical tools to perform log file content searching are generally crude and typically only employed as part of trouble-shooting exercises.As we move towards satellite constellations and fleets and add even more status information, the number of messages keeps growing. One mission now estimates that they could generate 3,000,000 messages per day 1 billion per year - for the life of their mission. What to do with those 1 billion messages? That is the challenge. With the recent technological advances in the management of large data sets, text-based processing, and data analytics, there are now capabilities that we can provide to the ground system engineers and satellite operators to address what we postulate are missed opportunities. Advanced real-time log analysis can allow us to be less reactionary in favor of being more proactive. Analytics goals include the ability to: Identify root cause of unexpected events, failures or error conditions enabled by correlating disparate data. Detect security breaches attempts before they are successful. Help admins ensure IT resources continue running optimally. Identify trends and patterns that may indicate impending failures or error conditions for valuable assets before they happen. Compare satellites in a fleet or constellation in terms of number of alarms, number of command sent to them, etc.. Answer questions like "Are the operations support needs increasing over the past year?" or "Have we seen this combination of alarm conditions before?" But really, once the tools are readily available the users will start realizing what can be done with their new powers. In this presentation we will show the results of analyzing millions of actual mission operations log messages, how the results can be displayed to the user, and how new products now available as open source can be applied to the challenges of large scale time-tagged text-based mission operations messages. Flight operations team members believe that this is a powerful new option for how they assess overall system and space asset health. Technical descriptions of the design, tools, and storage will be provided. One billion messages? Bring'em on!

Orsborne, Sharon↗
NASA NTRS2013

Evidence For Quasi-Periodic X-ray Dips From An Ultraluminous X-ray Source: Implications for the Binary Motion

We report results from long-term (approx.1240 days) X-ray (0.3-8.0 keV) monitoring of the ultraluminous X-ray source NGC 5408 X-1 with the Swift/X-Ray Telescope. Here we expand on earlier work by Strohmayer (2009) who used only a part of the present data set. Our primary results are: (1) the discovery of sharp, quasi-periodic, energy-independent dips in the X-ray intensity that recur on average every 243 days, (2) the detection of an energy dependent (variability amplitude decreases with increasing energy), quasi-sinusoidal X-ray modulation with a period of 112.6 +/- 4 days, the amplitude of which weakens during the second half of the light curve, and (3) spectral evidence for an increase in photoelectric absorption during the last continuous segment of the data. We interpret the X-ray modulations within the context of binary motion in analogy to that seen in high-inclination accreting X-ray binaries. If correct, this implies that NGC 5408 X-1 is in a binary with an orbital period of 243 +/- 23 days, in contrast to the 115.5 day quasi-sinusoidal period previously reported by Strohmayer (2009). We discuss the overall X-ray modulation within the framework of accretion via Roche-lobe overflow of the donor star. In addition, if the X-ray modulation is caused by vertically structured obscuring material in the accretion disk, this would imply a high value for the inclination of the orbit. A comparison with estimates from accreting X-ray binaries suggests an inclination > or approx.70deg. We note that, in principle, a precessing accretion disk could also produce the observed X-ray modulations.

Pasham, Dheeraj R.↗
DOE OSTI2021

Active virus-host interactions at sub-freezing temperatures in Arctic peat soil

Abstract Background Winter carbon loss in northern ecosystems is estimated to be greater than the average growing season carbon uptake and is primarily driven by microbial decomposers. Viruses modulate microbial carbon cycling via induced mortality and metabolic controls, but it is unknown whether viruses are active under winter conditions (anoxic and sub-freezing temperatures). Results We used stable isotope probing (SIP) targeted metagenomics to reveal the genomic potential of active soil microbial populations under simulated winter conditions, with an emphasis on viruses and virus-host dynamics. Arctic peat soils from the Bonanza Creek Long-Term Ecological Research site in Alaska were incubated under sub-freezing anoxic conditions with H 2 18 O or natural abundance water for 184 and 370 days. We sequenced 23 SIP-metagenomes and measured carbon dioxide (CO 2 ) efflux throughout the experiment. We identified 46 bacterial populations (spanning 9 phyla) and 243 viral populations that actively took up 18 O in soil and respired CO 2 throughout the incubation. Active bacterial populations represented only a small portion of the detected microbial community and were capable of fermentation and organic matter degradation. In contrast, active viral populations represented a large portion of the detected viral community and one third were linked to active bacterial populations. We identified 86 auxiliary metabolic genes and other environmentally relevant genes. The majority of these genes were carried by active viral populations and had diverse functions such as carbon utilization and scavenging that could provide their host with a fitness advantage for utilizing much-needed carbon sources or acquiring essential nutrients. Conclusions Overall, there was a stark difference in the identity and function of the active bacterial and viral community compared to the unlabeled community that would have been overlooked with a non-targeted standard metagenomic analysis. Our results illustrate that substantial active virus-host interactions occur in sub-freezing anoxic conditions and highlight viruses as a major community-structuring agent that likely modulates carbon loss in peat soils during winter, which may be pivotal for understanding the future fate of arctic soils' vast carbon stocks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Warming of the Willamette River, 1850–present: the effects of climate change and river system alterations

Abstract. Using archival research methods, we recovered and combined data from multiple sources to produce a unique, 140-year record of daily water temperature (Tw) in the lower Willamette River, Oregon (1881–1890, 1941–present). Additional daily weather and river flow records from the 1850s onwards are used to develop and validate a statistical regression model of Tw for 1850–2020. The model simulates the time-lagged response of Tw to air temperature and river flow and is calibrated for three distinct time periods: the late 19th, mid-20th, and early 21st centuries. Results show that Tw has trended upwards at 1.1 ∘C per century since the mid-19th century, with the largest shift in January and February (1.3 ∘C per century) and the smallest in May and June (∼ 0.8 ∘C per century). The duration that the river exceeds the ecologically important threshold of 20 ∘C has increased by about 20 d since the 1800s, to about 60 d yr−1. Moreover, cold-water days below 2 ∘C have virtually disappeared, and the river no longer freezes. Since 1900, changes are primarily correlated with increases in air temperature (Tw increase of 0.81 ± 0.25 ∘C) but also occur due to alterations in the river system such as depth increases from reservoirs (0.34 ± 0.12 ∘C). Managed release of water affects Tw seasonally, with an average reduction of up to 0.56 ∘C estimated for September. River system changes have decreased variability (σ) in daily minimum Tw by 0.44 ∘C, increased thermal memory, reduced interannual variability, and reduced the response to short-term meteorological forcing (e.g., heat waves). These changes fundamentally alter the response of Tw to climate change, posing additional stressors on fauna.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Dynamical Model for the Zodiacal Cloud and Sporadic Meteors

The solar system is dusty, and would become dustier over time as asteroids collide and comets disintegrate, except that small debris particles in interplanetary space do not last long. They can be ejected from the solar system by Jupiter, thermally destroyed near the Sun, or physically disrupted by collisions. Also, some are swept by the Earth (and other planets), producing meteors. Here we develop a dynamical model for the solar system meteoroids and use it to explain meteor radar observations. We find that the Jupiter Family Comets (JFCs) are the main source of the prominent concentrations of meteors arriving to the Earth from the helion and antihelion directions. To match the radiant and orbit distributions, as measured by the Canadian Meteor Orbit Radar (CMOR) and Advanced Meteor Orbit Radar (AMOR), our model implies that comets, and JFCs in particular, must frequently disintegrate when reaching orbits with low perihelion distance. Also, the collisional lifetimes of millimeter particles may be longer (approx. > 10(exp 5) yr at 1 AU) than postulated in the standard collisional models (approx 10(exp 4) yr at 1 AU), perhaps because these chondrule-sized meteoroids are stronger than thought before. Using observations of the Infrared Astronomical Satellite (IRAS) to calibrate the model, we find that the total cross section and mass of small meteoroids in the inner solar system are (1.7-3.5) 10(exp 11) sq km and approx. 4 10(exp 19) g, respectively, in a good agreement with previous studies. The mass input required to keep the Zodiacal Cloud (ZC) in a steady state is estimated to be approx. 10(exp 4)-10(exp 5) kg/s. The input is up to approx 10 times larger than found previously, mainly because particles released closer to the Sun have shorter collisional lifetimes, and need to be supplied at a faster rate. The total mass accreted by the Earth in particles between diameters D = 5 micron and 1 cm is found to be approx 15,000 tons/yr (factor of 2 uncertainty), which is a large share of the accretion flux measured by the Long Term Duration Facility (LDEF). Majority of JFC particles plunge into the upper atmosphere at <15 km/s speeds, should survive the atmospheric entry, and can produce micrometeorite falls. This could explain the compositional similarity of samples collected in the Antarctic ice and stratosphere, and those brought from comet Wild 2 by the Stardust spacecraft. Meteor radars such as CMOR and AMOR see only a fraction of the accretion flux (approx 1- 10% and approx 10-50%, respectively), because small particles impacting at low speeds produce ionization levels that are below these radars detection capabilities.

Nesvorny, David↗
NASA NTRS2017

How Long Do Satellites Need to Overlap? Evaluation of Climate Data Stability from Overlapping Satellite Records

Sensors on satellites provide unprecedented understanding of the Earth's climate system by measuring incoming solar radiation, as well as both passive and active observations of the entire Earth with outstanding spatial and temporal coverage. A common challenge with satellite observations is to quantify their ability to provide well-calibrated, long-term, stable records of the parameters they measure. Ground-based intercomparisons offer some insight, while reference observations and internal calibrations give further assistance for understanding long-term stability. A valuable tool for evaluating and developing long-term records from satellites is the examination of data from overlapping satellite missions. This paper addresses how the length of overlap affects the ability to identify an offset or a drift in the overlap of data between two sensors. Ozone and temperature data sets are used as examples showing that overlap data can differ by latitude and can change over time. New results are presented for the general case of sensor overlap by using Solar Radiation and Climate Experiment (SORCE) Spectral Irradiance Monitor (SIM) and Solar Stellar Irradiance Comparison Experiment (SOLSTICE) solar irradiance data as an example. To achieve a 1 % uncertainty in estimating the offset for these two instruments' measurement of the Mg II core (280 nm) requires approximately 5 months of overlap. For relative drift to be identified within 0.1 %/yr uncertainty (0.00008 W/sq m/nm/yr), the overlap for these two satellites would need to be 2.5 years. Additional overlap of satellite measurements is needed if, as is the case for solar monitoring, unexpected jumps occur adding uncertainty to both offsets and drifts; the additional length of time needed to account for a single jump in the overlap data may be as large as 50 % of the original overlap period in order to achieve the same desired confidence in the stability of the merged data set. Results presented here are directly applicable to satellite Earth observations. Approaches for Earth observations offer additional challenges due to the complexity of the observations, but Earth observations may also benefit from ancillary observations taken from ground-based and in situ sources. Difficult choices need to be made when monitoring approaches are considered; we outline some attempts at optimizing networks based on economic principles. The careful evaluation of monitoring overlap is important to the appropriate application of observational resources and to the usefulness of current and future observations.

climate data records↗
NASA NTRS1986

Space infrared telescope pointing control system. Infrared telescope tracking in the presence of target motion

The use of charge-coupled-devices, or CCD's, has been documented by a number of sources as an effective means of providing a measurement of spacecraft attitude with respect to the stars. A method exists of defocussing and interpolation of the resulting shape of a star image over a small subsection of a large CCD array. This yields an increase in the accuracy of the device by better than an order of magnitude over the case when the star image is focussed upon a single CCD pixel. This research examines the effect that image motion has upon the overall precision of this star sensor when applied to an orbiting infrared observatory. While CCD's collect energy within the visible spectrum of light, the targets of scientific interest may well have no appreciable visible emissions. Image motion has the effect of smearing the image of the star in the direction of motion during a particular sampling interval. The presence of image motion is incorporated into a Kalman filter for the system, and it is shown that the addition of a gyro command term is adequate to compensate for the effect of image motion in the measurement. The updated gyro model is included in this analysis, but has natural frequencies faster than the projected star tracker sample rate for dim stars. The system state equations are reduced by modelling gyro drift as a white noise process. There exists a tradeoff in selected star tracker sample time between the CCD, which has improved noise characteristics as sample time increases, and the gyro, which will potentially drift further between long attitude updates. A sample time which minimizes pointing estimation error exists for the random drift gyro model as well as for a random walk gyro model.

Powell, J. D.↗
DOE OSTI2023

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Analysis of the NASA White Sands Test Facility (WSTF) Test System for Friction-Ignition of Metallic Materials

Friction is a known ignition source for metals in oxygen-enriched atmospheres. The test system developed by the NASA White Sands Test Facility in response to ASTM G-94 has been used successfully to determine the relative ignition from friction of numerous metallic materials and metallic materials pairs. These results have been ranked in terms of a pressure-velocity product (PV) as measured under the prescribed test conditions. A high value of 4.1(exp 8) watts per square meter for Inconel MA 754 is used to imply resistance to friction ignition, whereas a low value of 1.04(exp 8) watts per square meter for stainless steel 304 is taken as indicating material susceptible to friction ignition. No attempt has been made to relate PV values to other material properties. This work reports the analysis of the WSTF friction-ignition test system for producing fundamental properties of metallic materials relating to ignition through friction. Three materials, aluminum, titanium, and nickel were tested in the WSTF frictional ignition instrument system under atmospheres of oxygen or nitrogen. Test conditions were modified to reach a steady state of operation, that is applied, the force was reduced and the rotational speed was reduced. Additional temperature measurements were made on the stator sample. The aluminum immediately galled on contact (reproducible) and the test was stopped. Titanium immediately ignited as a result of non-uniform contact of the stator and rotor. This was reproducible. A portion of the stator sampled burned, but the test continued. Temperature measurements on the stator were used to validate the mathematical model used for estimating the interface (stator/rotor) temperature. These interface temperature measurements and the associate thermal flux into the stator were used to distinguish material-phase transitions, chemical reaction, and mechanical work. The mechanical work was used to analyze surface asperities in the materials and to estimate a coefficient of fiction. The coefficient of fiction was analyzed in terms of material properties that is, hardness, Young's modulus and elasticity/plasticity of the material.

Shoffstall, Michael S.↗
DOE OSTI2024

Denoising Seismic Waveforms Using a Wavelet-Transform-Based Machine-Learning Method

Seismic waveform data recorded at stations can be thought of as a superposition of the signal from a source of interest and noise from other sources. Frequency‐based filtering methods for waveform denoising do not result in desired outcomes when the targeted signal and noise occupy similar frequency bands. Recently, denoising techniques based on deep‐learning convolutional neural networks (CNNs), in which a recorded waveform is decomposed into signal and noise components, have led to improved results. These CNN methods, which use short‐time Fourier transform representations of the time series, provide signal and noise masks for the input waveform. These masks are used to create denoised signal and designaled noise waveforms, respectively. However, advancements in the field of image denoising have shown the benefits of incorporating discrete wavelet transforms (DWTs) into CNN architectures to create multilevel wavelet CNN (MWCNN) models. The MWCNN model preserves the details of the input due to the good time–frequency localization of the DWT. In this report we use a data set of over 382,000 constructed seismograms recorded by the University of Utah Seismograph Stations network to compare the performance of CNN and MWCNN‐based denoising models. Evaluation of both models on constructed test data shows that the MWCNN model outperforms the CNN model in the ability to recover the ground‐truth signal component in terms of both waveform similarity and preservation of amplitude information. Model evaluation of real‐world data shows that both the CNN and MWCNN models outperform standard band‐pass filtering (BPF; average improvement in signal‐to‐noise ratio of 9.6 and 19.7 dB, respectively, with respect to BPF). Evaluation of continuous data suggests the MWCNN denoiser can improve both signal detection capabilities and phase arrival time estimates.

58 GEOSCIENCES↗
NASA NTRS2012

Applications of MODIS Fluorescent Line Height Measurements to Monitor Water Quality Trends and Algal Bloom Activity

Recent advances in satellite and airborne remote sensing, such as improvements in sensor and algorithm calibrations, processing techniques and atmospheric correction procedures have provided for increased coverage of remote-sensing, ocean-color products for coastal regions. In particular, for the Moderate Resolution Imaging Spectrometer (MODIS) sensor calibration updates, improved aerosol retrievals and new aerosol models has led to improved atmospheric correction algorithms for turbid waters and have improved the retrieval of ocean color in coastal waters. This has opened the way for studying ocean phenomena and processes at finer spatial scales, such as the interactions at the land-sea interface, trends in coastal water quality and algal blooms. Human population growth and changes in coastal management practices have brought about significant changes in the concentrations of organic and inorganic, particulate and dissolved substances entering the coastal ocean. There is increasing concern that these inputs have led to declines in water quality and have increase local concentrations of phytoplankton, which cause harmful algal blooms. In two case studies we present MODIS observations of fluorescence line height (FLH) to 1) assess trends in water quality for Tampa Bay, Florida and 2) illustrate seasonal and annual variability of algal bloom activity in Monterey Bay, California as well as document estuarine/riverine plume induced red tide events. In a comprehensive analysis of long term (2003-2011) in situ monitoring data and satellite imagery from Tampa Bay we assess the validity of the MODIS FLH product against chlorophyll-a and a suite of water quality parameters taken in a variety of conditions throughout a large optically complex estuarine system. A systematic analysis of sampling sites throughout the bay is undertaken to understand how the relationship between FLH and in situ chlorophyll-a responds to varying conditions and to develop a near decadal trend in water quality changes. In situ monitoring locations that correlated well with satellite imagery were in depths greater than seven meters and were located over five kilometers from shore. Water quality parameter of total nitrogen, phosphorous, turbidity and biological oxygen demand had high correlations with these sites, as well. Satellite FLH estimates show improving water quality from 2003-2007 with a slight decline up through 2011. Dinoflagellate blooms in Monterey Bay, California (USA) have recently increased in frequency and intensity. Nine years of MODIS FLH observations are used to describe the annual and seasonal variability of bloom activity within the Bay. Three classes of MODIS algorithms were correlated against in situ chlorophyll measurements. The FLH algorithm provided the most robust estimate of bloom activity. Elevated concentrations of phytoplankton were evident during the months of August-November, a period during which increased occurrences of dinoflagellate blooms have been observed in situ. Seasonal patterns of FLH show the on- and offshore movement of areas of high phytoplankton biomass between oceanographic seasons. Higher concentrations of phytoplankton are also evident in the vicinity of the land-based nutrient sources and outflows, and the cyclonic bay-wide circulation can transport these nutrients to the northern Bay bloom incubation region. Both of these case studies illustrate the utility MODIS FLH observations in supporting management decisions in coastal and estuarine waters.

Fischer, Andrew↗
DOE OSTI2022

Experimental study of the Si 30 ( He 3 , d ) P 31 reaction and thermonuclear reaction rate of Si 30 ( p , γ ) P 31

Background: Abundance anomalies in some globular clusters, such as the enhancement of potassium and the depletion of magnesium, can be explained in terms of an earlier generation of stars polluting the presently observed ones. It was shown that the potential range of temperatures and densities of the polluting sites depends on the strength of a few critical reaction rates. The 30 Si(p,γ) 31 P reaction has been identified as one of these important reactions. Purpose: The key ingredient for evaluating the thermonuclear 30 Si(p,γ) 31 P reaction rate is the strength of the resonances which, at low energy, are proportional to their proton width. Therefore, the goal of this work is to determine the proton widths of unbound 31 P states. Method: States in 31 P were studied at the Maier-Leibnitz-Laboratorium using the one-proton 30 Si( 3 He,d) 31 P transfer reaction. Deuterons were detected with the high resolution one quadrupole three dipoles (Q3D) magnetic spectrometer. Angular distribution and spectroscopic factors were extracted for 27 states, and proton widths and resonance strengths were calculated for the unbound states. Results: Several 31 P unbound states have been observed for the first time in a one-proton transfer reaction. Above 20 MK, the 30 Si(p,γ) 31 P reaction rate is now entirely estimated from the observed properties of 31 P states. The reaction rate uncertainty from all resonances other than the E$^{c.m.}_{r}$ = 149 keV resonance has been reduced down to less than a factor of 2 above that temperature. The unknown spin and parity of the E$^{c.m.}_{r}$ = 149 keV resonance dominates the uncertainty in the rate in the relevant temperature range. Conclusion: Here, the remaining source of uncertainty on the 30 Si(p,γ) 31 P reaction rate comes from the unknown spin and parity of the E$^{c.m.}_{r}$ = 149 keV resonance which can change the reaction rate by a factor of 10 in the temperature range of interest.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2007

Advanced Systems for Monitoring Underwater Sounds

The term "Passive Acoustic Monitoring System" (PAMS) describes a developmental sensing-and-data-acquisition system for recording underwater sounds. The sounds (more precisely, digitized and preprocessed versions from acoustic transducers) are subsequently analyzed by a combination of data processing and interpretation to identify and/or, in some cases, to locate the sources of those sounds. PAMS was originally designed to locate the sources such as fish of species that one knows or seeks to identify. The PAMS unit could also be used to locate other sources, for example, marine life, human divers, and/or vessels. The underlying principles of passive acoustic sensing and analyzing acoustic-signal data in conjunction with temperature and salinity data are not new and not unique to PAMS. Part of the uniqueness of the PAMS design is that it is the first deep-sea instrumentation design to provide a capability for studying soniferous marine animals (especially fish) over the wide depth range described below. The uniqueness of PAMS also lies partly in a synergistic combination of advanced sensing, packaging, and data-processing design features with features adapted from proven marine instrumentation systems. This combination affords a versatility that enables adaptation to a variety of undersea missions using a variety of sensors. The interpretation of acoustic data can include visual inspection of power-spectrum plots for identification of spectral signatures of known biological species or artificial sources. Alternatively or in addition, data analysis could include determination of relative times of arrival of signals at different acoustic sensors arrayed at known locations. From these times of arrival, locations of acoustic sources (and errors in those locations) can be estimated. Estimates of relative locations of sources and sensors can be refined through analysis of the attenuation of sound in the intervening water in combination with water-temperature and salinity data acquired by instrumentation systems other than PAMS. A PAMS is packaged as a battery-powered unit, mated with external sensors, that can operate in the ocean at any depth from 2 m to 1 km. A PAMS includes a pressure housing, a deep-sea battery, a hydrophone (which is one of the mating external sensors), and an external monitor and keyboard box. In addition to acoustic transducers, external sensors can include temperature probes and, potentially, underwater cameras. The pressure housing contains a computer that includes a hard drive, DC-to- DC power converters, a post-amplifier board, a sound card, and a universal serial bus (USB) 4-port hub.

Lane, Michael↗