Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source Term Analysis”

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 307 records · Page 17

Search for scattered X-ray halos around variable sources - The X-ray halo of Cygnus X-1

The results of a program to search for the presence of halos due to X-ray scattering by interstellar grains in the line-of-sight to variable X-ray sources are reported. As part of this program, four Einstein HRI images of Cyg X-1 were examined. The analysis technique exploits the intrinsic aperiodic variability of this source to map the point response function of the optics. Any scattered halo present will not reflect short-term central source time variability, since such variability is smoothed by differential time delays of order days. Thus, a residual, nonvariable component to the surface brightness distribution (comprising about 12 percent or more of the source flux) is interpreted as a scattered halo. The Cyg X-1 halo is consistent with those of other sources found in previous studies using different techniques. Comparison is also made with a scattering model and, despite uncertainties in source spectrum and distance, reasonable agreement with the observatons is found using a standard interstellar grain model. The potential of X-ray scattering as a probe of the properties of interstellar grains is demonstrated.

Bode, M. F.↗

On the numerical treatment of nonlinear source terms in reaction-convection equations

The objectives of this paper are to investigate how various numerical treatments of the nonlinear source term in a model reaction-convection equation can affect the stability of steady-state numerical solutions and to show under what conditions the conventional linearized analysis breaks down. The underlying goal is to provide part of the basic building blocks toward the ultimate goal of constructing suitable numerical schemes for hypersonic reacting flows, combustions and certain turbulence models in compressible Navier-Stokes computations. It can be shown that nonlinear analysis uncovers much of the nonlinear phenomena which linearized analysis is not capable of predicting in a model reaction-convection equation.

Lafon, A.↗

Graphite Baseline Status

Subtask Relevancy Inherent variation exists in nuclear-grade graphite due to processing techniques and raw material sources Comprehensive graphite properties database does not exist that fully establishes physical and mechanical property relationships Predictability of long-term behavior in a nuclear environment requires full analysis of variability Technical Approach, Accomplishments/Results Perform comprehensive mechanical and physical properties testing to populate a fully-validated database Grade selection is based upon candidate graphites from four international suppliers 15,000+ property values measured to date Physical properties testing that is based upon AGC geometries is a key component to establishing baseline/AGC data link

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Advanced Test Reactor Experiment Safety Analysis

This presentation contains general information about the Advanced Test Reactor experiment safety analysis program. Accidents and dose consequence are discussed in general terms. Pictures of ATR facilities and experiments are included. In cases where experiments are shown, identification numbers have been removed. Department of Energy limits on dose consequence are provided. These limits are found in many non-CUI documents and a derived from publicly available sources. A summary of the presentation content is as follows: -Define Experiment Safety Analysis -Discuss what we are protecting and the applicable limits -Show how safety analysis fits into the experiment process -Explore the scope of experiment activities and associated accidents -Present how we categorize accidents and establish accident frequencies -Identify how controls are implemented and process improvements

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Near-Real-Time Statistical Analysis and Visualization of Streamflow from a Deep-Learning Rainfall-Runoff Model

Near-real-time (NRT) streamflow data are critical importance for timely water resources management. Here, we developed an open-source tool, FlowStats, for NRT streamflow analysis and visualization in Germany, based on NRT meteorological data from the German Weather Service and simulated streamflow from a long short-term memory neural network (LSTM). The LSTM model achieved very good overall performance, median NSE of 0.80 for the test period across 1,479 catchments. FlowStats provides options for deriving various streamflow statistics, from normal and abnormal streamflow detection to drought and flood analyses. An example analysis from FlowStats revealed widespread below-normal to extreme low-flow conditions across Germany from March to May 2025, which weakened from June to September 2025. Drought analysis for September 2025 highlighted severe to extreme drought conditions in northwestern Germany, while flood classifications indicated that high-flow events occurred in southwestern Germany. FlowStats can be used for various hydrological assessments to support water resources management.

Hydrological modeling↗

Unveiling the chemical fingerprint of phosphorus-rich stars: I. In the infrared region of APOGEE-2

The origin of phosphorus, one of the essential elements for life on Earth, is currently unknown. Prevalent models of Galactic chemical evolution (GCE) systematically underestimate the amount of P compared to observations, especially at low metallicities. The recently discovered P-rich ([P/Fe] ≳ 1.2 dex) and metal-poor ([Fe/H] ≃ –1.0 dex) giants further challenge the GCE models, calling current theories on stellar nucleosynthesis into question. Since the observed low-mass giants are not expected to produce their high P contents themselves, our primary goal is to find clues on their progenitor or polluter. By increasing the number of known P-rich stars, we aim to narrow down a statistically reliable chemical abundance pattern that defines these peculiar stars. In this way, we place more robust constraints on the nucleosynthetic mechanism that causes the unusually high P abundances. In the long term, identifying the progenitor of the P-rich stars may contribute to the search for the source of P in our Galaxy. We performed a detailed chemical abundance analysis based on the high-resolution near-infrared (H band) spectra from the latest data release (DR17) of the APOGEE-2 survey. Employing the BACCHUS code, we measured the abundances of 13 elements in the inspected sample, which is mainly composed of a recent collection of Si-enhanced giants. We also analyzed the orbital motions and compared the abundance results to possible nucleosynthetic formation scenarios, and also to detailed GCE models. These models were produced with the OMEGA+ chemical evolution code, using four different massive star yield sets to investigate different scenarios for massive star evolution. We enlarged the sample of confirmed P-rich stars from 16 to a group of 78 giants, which represents the largest sample of P-rich stars to date. The sample includes the first detection of a P-rich star in a Galactic globular cluster. Significant enhancements in O, Al, Si, and Ce, as well as systematic correlations among the studied elements, unveil the unique chemical fingerprint of the P-rich stars. In contrast, the high [Mg/Fe] and [(C+N)/Fe] found in some of the P-rich stars with respect to P-normal stars is not confirmed over the full sample because of the current uncertainties. Strikingly, the strong overabundance in the α-element Si is accompanied by normal Ca and S abundances. This is at odds with current stellar nucleosynthesis models of massive stars. Our analysis of the orbital motion showed that the P-rich stars do not belong to a locally specific population in the Galaxy. In addition, we confirm that the majority of the sample stars are not part of binary systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Uranium Immobility in Tims Branch Wetland

The objective of this study was to address an item identified by the Low-Level Waste Disposal Facility Federal Review Group (LFRG) panel during their review of the Savannah River Site’s Composite Analysis. Their concern was that the Tims Branch wetland, which has elevated levels of uranium (U) originating from site operations between 1952 to 1985, may be a secondary source term, posing a human health risk. This report briefly describes three studies dealing with: 1) quantifying the mass of U remaining in the wetland, 2) quantifying the tendency of U to release from Tims Branch sediment, and 3) describing the distribution of U in aqueous and suspended solids in Tims Branch water during episodic storm events. The Tims Branch wetland has been extremely effective for the last 50 years at immobilizing U released to the system. A detailed map including 270,000 gamma spectra of the U distribution in the wetland indicated that 84% of the initial 43,500 kg U introduced in the system remain there after 50 years (Kaplan et al., 2020). The effectiveness of immobilizing U is supported by stream water measurements collected over the last several decades indicating that U concentrations have been consistently below drinking water regulatory limits (i.e., the EPA Maximum Contaminant Level, 30 pg/L (SRNS, 2019)). Furthermore, using contaminated Tims Branch sediments, in situ desorption Kd values (ratio of U sediment/U Water) were extremely high, ranging from 104to 106 L/kg, suggesting the U is strongly bound to these wetland sediments. Finally, a detailed study of Tims Branch stream water properties during episodic rain events indicated that essentially all the U transported in this stream was associated with suspended solids (Batson et al., 1996). Together these findings indicate that the Tims Branch wetland is extremely effective in its present state at immobilizing U and that the small masses of U released downstream will be predominantly in association with suspended solids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

User’s Manual for RESRAD-RDD&IND Code Version 2: Vol. 2—User’s Guide for RESRAD-RDD&IND Code

Version 2.0 of the RESRAD-RDD&IND computer code is designed to support the implementation of protective action guides (PAGs) after a nuclear emergency incident including a radiological dispersal device (RDD) and/or an improvised nuclear device (IND) incident (EPA 2017). Eight different group types, addressing various decisions, are available for selection. The RESRAD-RDD&IND code calculates radiological doses, stay times, etc., for the selected group that the user wishes to focus on. (That is, the results for all the groups are not calculated simultaneously, and the input for those other groups do not matter, although some parameter values are shared between groups.) Version 2.0 has a user-friendly interface so that the RESRAD-RDD&IND code can be used with minimal training. For example, the user can select the major characteristics of the problem-event type, source term, and decision type from the left side of the interface and then calculate the results with the default assumptions for the exposure scenarios. More in-depth analysis would include specifying site-specific exposure scenario characteristics in the right side of the interface. The procedures for data entry and results viewing are self-explanatory. This is because common window maneuvering features and text instructions were incorporated in the interface design. General and context-specific help are available to aid users entering parameter values, as well. The RESRAD-RDD&IND computer code gives the user the option to select either an RDD or IND incident for analysis. For an RDD event analysis, 11 radionuclides (Am-241, Cf-252, Cm-244, Co-60, Cs-137, Ir-192, Po-210, Pu-238, Pu-239, Ra-226, and Sr-90) are included. These 11 radionuclides are the radionuclides most likely used for an RDD. More than 90 radionuclides can be selected for an IND event analysis. Initial default concentrations are provided for 44 radionuclides for a uranium-fueled IND event. These 44 radionuclides are those that would contribute significantly to the radiation dose associated with a uranium-fueled bomb detonation. The radionuclides generated from ingrowth of these 44 initial radionuclides are also automatically included in the analysis. Pu-239, Cs-134m, Ru-105, and Rb-89 and their progeny can be selected for analysis if they are detected and their concentrations are determined. This user’s guide, which is Volume 2 of the User’s Manual for RESRAD-RDD&IND Code Version 2, provides instructions to users on how to install the RESRAD-RDD&IND code, navigate the interface, and use the various features, including those discussed above, to set up an analysis and view/print the results in text outputs. Volume 1 of the User’s Manual for RESRAD-RDD&IND Code Version 2 (Yu et al. 2026), which contains descriptions of the methodology and theoretical basis for dose modeling and the mathematical equations implemented in the code, can be accessed and viewed through the Help menu in the code or can be downloaded from the RESRAD website (https://resrad.evs.anl.gov).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Design of optimum acoustic treatment for rectangular ducts with flow

A design optimization technique for acoustic treatment in rectangular ducts with uniform mean flow is presented. The technique is based on the acoustic wave solution in terms of series of characteristic duct modes. The analysis allows multiple axial treatment sections along the length of the duct and requires a known modal characterization of the sound source. Conditions of acoustic pressure and acoustic velocity continuity are used to match modal solutions at planes of impedance discontinuity in the duct. Experimental techniques for obtaining this modal characterization are presented. Using duct modes measured at the source plane, the optimization technique is exercised to design an optimized single element liner in a case without mean flow, and optimized single and dual element liners in cases with mean flow. The validity of the program for predicting noise suppression is demonstrated by comparing analytical predictions with measured data for several (non-optimum) cases. Application to treatment design in turbomachinery exhaust ducts is considered.

Motsinger, R. E.↗

The analysis of delays in simulator digital computing systems. Volume 1: Formulation of an analysis approach using a central example simulator model

The effects of spurious delays in real time digital computing systems are examined. Various sources of spurious delays are defined and analyzed using an extant simulator system as an example. A specific analysis procedure is set forth and four cases are viewed in terms of their time and frequency domain characteristics. Numerical solutions are obtained for three single rate one- and two-computer examples, and the analysis problem is formulated for a two-rate, two-computer example.

Heffley, R. K.↗

Helicopter main-rotor noise: Determination of source contributions using scaled model data

Acoustic data from a test of a 40 percent model MBB BO-105 helicopter main rotor are scaled to equivalent full-scale flyover cases. The test was conducted in the anechoic open test section of the German-Dutch Windtunnel (DNW). The measured data are in the form of acoustic pressure time histories and spectra from two out-of-flow microphones underneath and foward of the model. These are scaled to correspond to measurements made at locations 150 m below the flight path of a full-scale rotor. For the scaled data, a detailed analysis is given for the identification in the data of the noise contributions from different rotor noise sources. Key results include a component breakdown of the noise contributions, in terms of noise criteria calculations of a weighted sound pressure level (dBA) and perceived noise level (PNL), as functions of rotor advance ratio and descent angle. It is shown for the scaled rotor that, during descent, impulsive blade-vortex interaction (BVI) noise is the dominant contributor to the noise. In level flight and mild climb, broadband blade-turbulent wake interaction (BWI) noise is dominant due to the absence of BVI activity. At high climb angles, BWI is reduced and self-noise from blade boundary-layer turbulence becomes the most prominent.

Brooks, Thomas F.↗

Amazon Basin ozone and aerosol - Wet season observations

The ozone and aerosol number density data obtained from aircraft flights over the Amazon Basin (AB) during April and May of 1987 are summarized with particular attention given to the general characteristics of ozone and aerosol for wet season (WS), compared to dry season (DS), meteorology. It is concluded that the AB provides favorable conditions for the production of hydroxyl radical and subsequent chemical and photochemical reactions which generate or consume key atmospheric trace species. It is determined that a negative ozone gradient exists within the mixed layer (ML). WSML ozone is about 7 ppbv lower than that of DS. The net sink characteristics of the WSML are seen throughout the troposphere of the AB in that the ozone at 3-4 km altitude is typically 15-25 ppbv as compared to DS values of 30-35 ppbv. In terms of aerosol source characteristics of the AB, ML aerosols are 5-10 times higher than observed in the troposphere with ML values of 100-200 aerosols/cu cm. An analysis of aerosol size distribution data is presented.

Gregory, Gerald L.↗

Analysis and interpretation of Wide Field Planetary Camera outgas data collected in the temperature range from -20 C to -100 C

The data from the Wide Field Planetary Camera (WFPC) are analyzed to determine the absolute source-emission parameters and contaminant reemission parameters of outgassing mass accumulation. The data correspond to the mass accumulated on a -100 C quartz crystal microbalance and then desorbed at -50 C as a function of time. The parameters investigated include contaminant mass, first-order rate constants, and characteristic energies. The data are modeled in terms of a minimum number of species characterized by source-rate and reemission-rate parameters. The modeling technique is found to produce adequate fits, and the contaminants that accumulate at -100 C and desorb at -50 C do not accumulate at -50 C. The data analysis is found to be an adequate technique for assessing outgassing parameters for the prediction of WFPC internal contamination.

Barengoltz, Jack B.↗

The Severe Weather Outbreak of 10 November 2002: Lightning and Radar Analysis of Storms in the Deep South

On the afternoon and evening of 10 November 2002, the Midwest and Deep South were struck by a major outbreak of severe storms that produced some 80 tornadoes. In terms of number of tornadoes, this was the largest outbreak in the United States since November 1992. Some 32 of the tornadoes occurred in Tennessee, Mississippi, Alabama and Georgia, including several long-track killers. We use the North Alabama Lightning Mapping Array (LMA) and other data sources to perform a comprehensive analysis of the structure and evolution of the outbreak. Most of the Southern tornadoes occurred in isolated, fast-moving supercell storms that formed in warm, moist air ahead of a major cold front. Storms tended to form in lines parallel to storm cell motion, resulting in many communities being hit multiple times by severe storms that evening. Supercells in Tennessee produced numerous strong tornadoes with short to medium-length track paths, while the supercells further south produced several very long-track tornadoes. Radar data indicate that the Tennessee storms tended to split frequently, apparently limiting their ability to sustain long-lived tornadoes, while storms further south split at most one time. The differences between these storms appear to be related to the presence of stronger jetstream winds in Tennessee relative to those present in Mississippi, Alabama and Georgia. LMA-derived flash rates associated with most of the supercell storm cores were about 1-2 flashes per second. Rapid increases in lightning rates (or "jumps") occurred prior to tornado touchdown in many instances. Lightning "holes" (lightning-free regions associated with the echo-free vault) occurred in two of the Tennessee supercells. The complexity of the relationship between lightning and storm severity is revealed by the behavior of one Alabama supercell, which produced a peak flash rate of nearly 14 flashes per second, well after the end of its long-track tornado, while interacting and ultimately merging with a daughter supercell on its southwest flank. Close examination of this powerful storm indicates that its prodigious flash rate was the result of strong flash activity over an unusually large area, rather than a concentrated core of extremely high flash rate activity.

Buechler, D. E.↗

Dose Rate Analysis of the WCS Consolidated Interim Storage Facility

This report describes a confirmatory calculation for a proposed consolidated interim storage facility (CISF) to support the US Nuclear Regulatory Commission’s review of the license application submitted by the Interim Storage Partners. The scope of this report is limited to an assessment of the annual dose from direct radiation associated with Phase 1 of the Waste Control Specialists (WCS) CISF and the determination of the minimum controlled area boundary based on the 10 CFR 72.104 limit of 25 mrem annual dose to the whole body. The confirmatory dose rate calculations used the source term and shielding calculation capabilities of the SCALE 6.2.3 computer code system. The calculation method employed in the confirmatory calculations uses a detailed Monte Carlo radiation transport simulation from source to dose rate. This method differs from the two-step method used in the safety analysis report (SAR), which requires determination of the photon and neutron energy and angular distributions on the cask external surface and then the use of this surface source in a new radiation transport calculation to determine dose rate as a function of distance from the storage facility. A simulation was made using the complete site geometry (all casks present) but with only one cask containing source, and a total of 467 independent calculations were performed to obtain the dose rate maps produced by each storage cask. Basic input data and assumptions used in the SAR (e.g., cask design parameters and design basis assembly characteristics) were used in this confirmatory calculation for consistency with the SAR. Local atmospheric conditions for skyshine calculations were simulated using US standard atmosphere data, which describe average air temperature, pressure, and density as a function of altitude. The determined minimum required distances from the WCS Phase 1 CISF to a contour defined by the 25 mrem annual dose limit were approximately 623 m (0.387 miles) in the N-NE direction, 594 m (0.369 miles) in the SSW direction, and 533 m (0.331 miles) in the E-SE direction and the W-NW direction. The annual dose results of the confirmatory calculation are in relatively good agreement with the annual dose values presented in the WCS SAR. The annual dose produced by the Phase 1 CISF at the nearest site boundary, 0.75 miles away from the CISF, was estimated by this study to be approximately 6.02E-02 ± 2.0E-04 mrem, which is significantly smaller than the 25 mrem annual dose limit provided in 10 CFR 72.104. The SAR annual dose rate value for the nearest site boundary was 7.52E-2 mrem, which included 7.73E-03 mrem due to postulated leakage of the FO-, FC-, and FF-canisters. A comparison with the dose rate values reported in the SAR for 16 detector locations, D1 through D16, on the protected area boundary is also included in this report. This comparison shows that the SAR determined higher annual dose values for detectors D1 through D14 by 12% to 45% and lower annual dose values for detectors D15 and D16 by 18% and 8%, respectively, as compared with this ORNL confirmatory calculation. Sensitivity studies have been performed to support a comprehensive review of the application. In these studies, the effects on dose rate of the modeled air volume, density, and humidity were analyzed. The sensitivity calculations showed that the air volume used in the Monte Carlo simulations was adequate and that a larger air volume would not further increase the dose rate due to skyshine contributions. The sensitivity analysis showed that simulation of local atmospheric conditions is important for the accuracy of dose rate estimates at large distances from the facility. Air humidity was demonstrated to have negligible effects on dose rate values for the simulated local atmospheric conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗

WSF PISA example for SAFER pilot

On October 15, 2021, the Radioactive and Hazardous Waste Management Facility Manager declared a Potential Inadequacy in the Safety Analysis (PISA) due to a new information relating to an error discovered in calculation AB-WSF-20-002. Calculation AB-WSF-20-002 is used to determine source terms and heat release rates for vehicle and aircraft impacts with fuel fires involving TRU waste containers. The source term for an aircraft’s impact has different source terms factors (i.e., DR, ARF, RF, LPF) that are applied to multiple categories of drums impacted in the accident scenario based on the physical stresses of the aircraft impact and involvement in the subsequent fuel pool fire. The source terms of the multiple categories of affected drums are summed to determine the final source term for the accident scenario. Certain categories of affected drums are modeled to lose their lid, causing material to eject and burn unconfined on the ground. For these specific categories, a lower ARF value was incorrectly applied in AB-WSF-20-002 (i.e., an ARF value of 1E-3 was used instead of the DOE-STD-5506-2007 directed ARF value of 1E-2). For these categories, the ARF is applied to a relatively small number of low-activity drums (i.e., 17 drums at 2.9 PE-Ci each) compared to the total number of drums affected (i.e., 166 drums totaling 621.8 PE-Ci). Therefore, this error is not expected to significantly increase final dose for the accident scenario. No immediate actions or compensatory measures are required to maintain the facility in a safe condition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparison and analysis of aircraft measurements and mesoscale atmospheric chemistry model simulations of tropospheric ozone

The Regional Acid Deposition Model (RADM) has been applied to several of the field experiments which were part of the Acid Models Operational and Diagnostic Evaluation Study (Acid MODES). The experiment which was of particular interest with regards to ozone photochemistry involved horizontal zig-zag flight patterns (ZIPPER) over an area from the eastern Ohio River valley to the Adirondacks of New York. Model simulations by both the standard resolution RADM (delta x = 80 km) and the nested grid RADM (delta x = 26.7 km) compare well to measurements in the low emission regions in central Pennsylvania and upstate New York, but underestimate in the high emission upper Ohio River valley. The nested simulation does considerably better, however, than the coarse grid simulation in terms of horizontal pattern and concentration magnitudes. Analysis of NO(x) and HO(x) concentrations and photochemical products rates of ozone show that the model's response to large point source emissions is very unsystematic both spatially and temporally. This is due to the models instability to realistically simulate the small scale (subgrid) gradients in precursor concentrations in and around large point source plumes.

Pleim, Jonathan E.↗

A Multipole Expansion Method for Analyzing Lightning Field Changes

Changes in the surface electric field are frequently used to infer the locations and magnitudes of lightning-caused changes in thundercloud charge distributions. The traditional procedure is to assume that the charges that are effectively deposited by the flash can be modeled either as a single point charge (the Q model) or a point dipole (the P model). The Q model has four unknown parameters and provides a good description of many cloud-to-ground (CG) flashes. The P model has six unknown parameters and describes many intracloud (IC) discharges. In this paper we introduce a new analysis method that assumes that the change in the cloud charge can be described by a truncated multipole expansion, i.e., there are both monopole and dipole terms in the unknown source distribution, and both terms are applied simultaneously. This method can be used to analyze CG flashes that are accompanied by large changes in the cloud dipole moment and complex IC discharges. If there is enough information content in the measurements, the model can also be generalized to include quadrupole and higher order terms. The parameters of the charge moments are determined using a dme-dimensional grid search in combination with a linear inversion, and because of this, local minima in the error function and the associated solution ambiguities are avoided. The multipole method has been tested on computer-simulated sources and on natural lightning at the NASA Kennedy Space Center and U.S. Air Force Eastern Range.

Koshak, William J.↗