Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete linearization method”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Total metals, carbon, nitrogen & anion concentration data; Slate River & East River floodplains, Crested Butte, CO; May 2022-October 2022

This data package includes processed and undiluted measurements for metal, total carbon, total nitrogen, and anion concentrations from pore water (groundwater) and surface water samples from the Slate River and East River floodplains of Crested Butte, CO, focus field sites for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and October of 2022. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS (metals):Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA. The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solutions (SPEX Certiprep, Metuchen, NJ). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace elements in water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. Lastly, a suitable internal standard (usually Rh, In, Ga or Ge) is added using on-line addition into the sample line and mixing tee.Analysis by Shimadzu TOC-L (TOC/TN):The TOC-L system is a combustion technique where liquid samples are injected and combusted into CO2 for carbon detection by non-dispersive infrared (NDIR) and NO for detection by chemiluminescence. A calibration curve using five standard solutions between 0.1 and 7 ppm for carbon and 0.05 and 3.5 ppm for nitrogen is made for each type of measurement with a linearity >0.99. All samples, standards, and QC’s are prepared in 24mL scintillation vials that have been baked for 4hrs at 475 Cº and made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.Analysis by Ion Chromatography (Anions):The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with sodium carbonate eluent. A calibration curve using five standard solutions between 5 and 250 umol/L is made with a linearity >0.99. Standards and QC’s are prepared in 15mL polypropylene conical tubes, pipetted along with the samples into 1.5mL polypropylene vials. Dilutions are made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Changes in ribbon synapses and rough endoplasmic reticulum of rat utricular macular hair cells in weightlessness

This study combined ultrastructural and statistical methods to learn the effects of weightlessness on rat utricular maculae. A principle aim was to determine whether weightlessness chiefly affects ribbon synapses of type II cells, since the cells communicate predominantly with branches of primary vestibular afferent endings. Maculae were microdissected from flight and ground control rat inner ears collected on day 13 of a 14-day spaceflight (F13), landing day (R0) and day 14 postflight (R14) and were prepared for ultrastructural study. Ribbon synapses were counted in hair cells examined in a Zeiss 902 transmission electron microscope. Significance of synaptic mean differences was determined for all hair cells contained within 100 section series, and for a subset of complete hair cells, using SuperANOVA software. The synaptic mean for all type II hair cells of F13 flight rats increased by 100%, and that for complete cells by 200%. Type I cells were less affected, with synaptic mean differences statistically insignificant in complete cells. Synapse deletion began within 8 h upon return to Earth. Additionally, hair cell laminated rough endoplasmic reticulum of flight rats was reversibly disorganized on R0. Results support the thesis that synapses in type II hair cells are uniquely affected by altered gravity. Type II hair cells may be chiefly sensors of gravitational and type I cells of translational linear accelerations.

short duration↗
NASA NTRS2007

Analysis of Crossflow Transition Flight Experiment aboard the Pegasus Launch Vehicle

The Pegasus wing-glove flight experiment was designed to provide crossflow transition data at high Mach numbers, specifically to help validate stability based predictions for transition onset in a flight environment. This paper provides an analysis of the flight experiment, with emphasis on computational results for crossflow disturbances and the correlation of disturbance growth factors with in-flight transition locations via the e(sup N) method. Implications of the flight data for attachment line stability are also examined. Analysis of the thermocouple data reveals that transition (from turbulent to laminar flow) was first detected during the ascending flight of the rocket when the free stream Mach number exceeded about 4. Therefore, computations have been performed for flight Mach numbers of 4.13, 4.35, 4.56 and 4.99. Due to continually decreasing unit Reynolds number at higher altitudes, the entire wing-glove boundary layer became laminar at the highest flight Mach number computed above. In contrast, the boundary layer flow over the inboard tile region remained transitional up to and somewhat beyond the time of laminarization over the instrumented glove region. Linear stability predictions confirmed that the tile boundary layer is indeed more unstable to crossflow disturbances than the much colder stainless steel glove boundary layer. The transition locations based on thermocouple data from both the glove and the tile regions are found to correlate with stationary-crossflow N-factors within the range of 7 to 12.4 and with traveling mode N-factors between 7.6 and 14.1. Data from the thermocouples and hot film sensors indicates that transition from turbulent to laminar flow (i.e., laminarization) at a fixed point over the glove is generally completed within a flight time interval of 3 seconds. However, the times at which transition begins and ends as inferred from the hot film sensors are found to differ by about 2 seconds from the corresponding estimates based on the thermocouple data.

Malik, Mujeeb R.↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Grossman, Bernard↗
DOE OSTI2019

Non-linear evolution in QCD at high-energy beyond leading order

The standard formulation of the high-energy evolution in perturbative QCD, based on the Balitsky-Kovchegov equation, is known to suffer from severe instabilities associated with radiative corrections enhanced by double transverse logarithms, which occur in all orders starting with the next-to-leading one. Over the last years, several methods have been devised to resum such corrections by enforcing the time-ordering of the successive gluon emissions. We observe that the instability problem is not fully cured by these methods: various prescriptions for performing the resummation lead to very different physical results and thus lack of predictive power. We argue that this problem can be avoided by using the rapidity of the dense target (which corresponds to Bjorken x) instead of that of the dilute projectile as an evolution time. This automatically ensures the proper time-ordering and also allows for a direct physical interpretation of the results. We explicitly perform this change of variables at NLO. We observe the emergence of a new class of double logarithmic corrections, potentially leading to instabilities, which are however less severe, since disfavoured by the typical BK evolution for “dilute-dense” scattering. We propose several prescriptions for resumming these new double-logarithms to all orders and find only little scheme dependence: different prescriptions lead to results which are consistent to each other to the accuracy of interest. We restore full NLO accuracy by completing one of the resummed equations (non-local in rapidity) with the remaining NLO corrections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

How Useful are Conventional I–Vs for Performance Calibration of Single- and Two-Junction Perovskite Solar Cells? A Statistical Analysis of Performance Data on ≈200 Cells from 30 Global Sources

As perovskite photovoltaics (PV) advance from the laboratory to commercial prototypes, their accurate and reliable performance testing is becoming increasingly important. The well-documented dynamic response of perovskite solar cells to an external applied voltage has led to the development of steady-state performance measurement methods; however, these methods have not been widely adopted by the perovskite PV community. A key reason for this is that steady-state measurement methods take tens of minutes to complete, as opposed to conventional "fast" current-voltage (I-V) measurements usually lasting a few seconds. Fast I-Vs arise from a snapshot, almost always not a steady-state condition of the device; however, given their widespread use, the question arises: how do performance parameters of perovskite PV compare when measured with fast I-V and with a steady-state method? Results compiled from approximately 200 perovskite PV cells, including single junction, and two-terminal perovskite-perovskite and perovskite-Si tandems, show that fast I-Vs can provide a useful measure of the open-circuit voltage of the devices, while the short-circuit current and the overall efficiency can be widely misestimated. Here, the implications of these findings on performance testing protocols are discussed and possible options for fast and accurate testing of perovskite PV are proposed.

14 SOLAR ENERGY↗
DOE OSTI2019

The International Linear Collider. A Global Project

A large, world-wide community of physicists is working to realise an exceptional physics program of energy-frontier, electron-positron collisions with the International Linear Collider (ILC). This program will begin with a central focus on high-precision and model-independent measurements of the Higgs boson couplings. This method of searching for new physics beyond the Standard Model is orthogonal to and complements the LHC physics program. The ILC at 250 GeV will also search for direct new physics in exotic Higgs decays and in pair-production of weakly interacting particles. Polarised electron and positron beams add unique opportunities to the physics reach. The ILC can be upgraded to higher energy, enabling precision studies of the top quark and measurement of the top Yukawa coupling and the Higgs self-coupling. The key accelerator technology, superconducting radio-frequency cavities, has matured. Optimised collider and detector designs, and associated physics analyses, were presented in the ILC Technical Design Report, signed by 2400 scientists. There is a strong interest in Japan to host this international effort. A detailed review of the many aspects of the project is nearing a conclusion in Japan. Now the Japanese government is preparing for a decision on the next phase of international negotiations, that could lead to a project start within a few years. The potential timeline of the ILC project includes an initial phase of about 4 years to obtain international agreements, complete engineering design and prepare construction, and form the requisite international collaboration, followed by a construction phase of 9 years.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

Vivid Motor Imagery as an Adaptation Method for Head Turns on a Short-Arm Centrifuge

Artificial gravity (AG) has been proposed as a potential countermeasure to the debilitating physiological effects of long duration space flight. The most economical means of implementing AG may be through the use of a short-radius (2m or less) centrifuge. For such a device to produce gravitational forces comparable to those on earth requires rotation rates in excess of 20 revolutions per minute (rpm). Head turns made out of the plane of rotation at these rates, as may be necessary if exercise is combined with AG, result in cross-coupled stimuli (CCS) that cause adverse side effects including motion sickness, illusory sensations of motion, and inappropriate eye movements. Recent studies indicate that people can adapt to CCS and reduce these side effects by making multiple head turns during centrifuge sessions conducted over consecutive days. However, about 25% of the volunteers for these studies have difficulty tolerating the CCS adaptation paradigm and often drop out due to motion sickness symptoms. The goal of this investigation was to determine whether vivid motor imagery could be used as a pseudostimulus for adapting subjects to this unique environment. Twenty four healthy human subjects (14 males, 10 females), ranging in age from 21 to 48 years (mean 33, sd 7 years) took part in this study. The experimental stimuli were produced using the NASA JSC short-arm centrifuge (SAC). Subjects were oriented supinely on this device with the nose pointed toward the ceiling and head centered on the axis of rotation. Thus, centrifuge rotation was in the body roll plane. After ramp-up the SAC rotated clockwise at a constant rate of 23 rpm, producing a centrifugal force of approximately 1 g at the feet. Semicircular canal CCS were produced by having subjects make yaw head turns from the nose up (NU) position to the right ear down (RED) position and from RED to NU. Each head turn was completed in about one second, and a 30 second recovery period separated consecutive head movements. Participants were randomly assigned to one of three groups (n=8 per group): physical adapters (PA), mental adapters (MA), or a control group (CG). Each subject participated in a one hour test session on each of three consecutive days. Each test session consisted of an initial (preadaptation) period during which the subject performed six CCS maneuvers in the dark, followed by an adaptation period with internal lighting on the centrifuge, and a final (postadaptation) period during which six more CCS maneuvers were performed in the dark. For the PA group, the adaptation period consisted of performing 30 additional CCS maneuvers in the light. For the MA and CG group the centrifuge was ramped down to 0 rpm after the pre-adaptation period and ramped back up to 23 rpm before the post-adaptation period. For the both of these groups, the adaptation period consisted of making 30 CCS maneuvers in the light with the centrifuge stationary (so no cross-coupling occurred). MA group subjects were instructed to vividly imagine the provocative sensations produced by the preadaptation CCS maneuvers in terms of magnitude, duration, and direction of illusory body tilt, as well as any accompanying levels of motion sickness. CG group subjects were asked to answer low imagery content questions (trivial pursuit) during each adaptation period head turn. During the 30 second recovery following each head turn, psychophysical data were collected including self reports of motion sickness, magnitude and direction estimates of illusory body tilt, and the overall duration of these sensations. A multilevel mixed effects linear regression analysis performed on all response variables indicated that all three groups experienced some psychophysical adaptation across the three test sessions. For illusory tilt magnitude, the PA group exhibited the most overall adaptation, followed by the MA group, and the CG group. The slopes of these adaptation trajectories by group over day were significantly diffent from one another. For the perceived duration of sensations, the CG group again exhibited the least amount of adaptation. However, the rates of adaptation of the PA and the MA groups were indistinguishable, suggesting that the imagined pseudostimulus appeared to be just as effective a means of adaptation as the actual stimulus. The MA group's rate of adaptation to motion sickness symptoms was also comparable to the PA group. The use of vivid motor imagery may be an effective method for adapting to the illusory sensations and motion sickness symptoms produced by cross-coupled stimuli. For space-based AG applications, this technique may prove quite useful in retaining astronauts considered highly susceptible to motion sickness as it reduces the number of actual CCS required to attain adaptation.

Newby, N. J.↗
NASA NTRS2004

Comparisons of the Solutions for the Category 3--Problem 2: Cascade-Gust Interaction

Six different solutions were submitted for this benchmark problem. These were obtained using a variety of methods that can be conveniently categorized in two main groups, a nonlinear time-domain group and a linearized frequency-domain group. The first includes solutions submitted by (1) Hixon, (2) Nallasamy et. al, (3) Shieh et. al, and (4) Wang et. al, and the second includes solutions submitted by (5) Coupland and (6) Serrano et. al. Methods (1) and (2) use sixth order compact differencing schemes and the rest are essentially second order in space. With the exception of the solution submitted by Shieh et. al, all are individually discussed in great detail in the workshop proceedings. Comparisons of the submitted solutions with the benchmark solution are presented below. Due to differences in the level of solution detail provided to the author by the participants, the comparisons do not always include results from all submissions. It should be noted at the outset that, since the benchmark solution itself was numerically computed, the comparisons are somewhat subjective. In order to provide maximum latitude for the participants of the workshop, no restrictions were placed on the type of method that could be used to solve the problem. Neither was there were any stipulations to use a particular grid topology or grid density. Therefore, without a detailed study of the critical features of the computed solutions, it is not possible to make concrete statements about the relative merits of one method over another. Such a study is beyond the scope of the current exercise, especially since complete flowfield details were not provided to the author by all participants. Instead a package, containing the information about the benchmark solution (both the steady and unsteady parts of it), is included on the proceedings CD should the authors who submitted solutions for this problem wish to examine in detail the benchmark solution and compare their results to it.

Envia, Edmane↗
NASA NTRS2021

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗
DOE OSTI2021

The PAU Survey: an improved photo- z sample in the COSMOS field

Here we present – and make publicly available – accurate and precise photometric redshifts in the ACS footprint from the COSMOS field for objects with i AB ≤ 23. The redshifts are computed using a combination of narrow-band photometry from PAUS, a survey with 40 narrow bands spaced at $100\,\mathring{\rm A}$ intervals covering the range from 4500 to $8500\,\mathring{\rm A}$, and 26 broad, intermediate, and narrow bands covering the UV, visible and near-infrared spectrum from the COSMOS2015 catalogue. We introduce a new method that models the spectral energy distributions as a linear combination of continuum and emission-line templates and computes its Bayes evidence, integrating over the linear combinations. The correlation between the UV luminosity and the $\mathrm{O\,{\small II}}$ line is measured using the 66 available bands with the zCOSMOS spectroscopic sample, and used as a prior which constrains the relative flux between continuum and emission-line templates. The flux ratios between the $\mathrm{O\,{\small II}}$ line and H α , H β and $\mathrm{O\,{\small III}}$ are similarly measured and used to generate the emission-line templates. Comparing to public spectroscopic surveys via the quantity Δ z ≡ (z photo – z spec )/(1 + z spec ), we find the photometric redshifts to be more precise than previous estimates, with σ 68 (Δ z ) ≈ (0.003, 0.009) for galaxies at magnitude i AB ~ 18 and i AB ~ 23, respectively, which is three times and 1.66 times tighter than COSMOS2015. Additionally, we find the redshifts to be very accurate on average, yielding a median of the Δ z distribution compatible with |median(Δ z )| ≤ 0.001 at all redshifts and magnitudes considered. Both the added PAUS data and new methodology contribute significantly to the improved results. The catalogue produced with the technique presented here is expected to provide a robust redshift calibration for current and future lensing surveys, and allows one to probe galaxy formation physics in an unexplored luminosity-redshift regime, thanks to its combination of depth, completeness, and excellent redshift precision and accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Quantifying Rates of Erosion Using Historical Aerial Photos - 19478

Erosion can play an important role in the fate and transport of waste at sites where disposal of long-lived waste is anticipated. The statistical characterization of key processes related to erosion is essential to the understanding of site stability through time. In this work multi-temporal historical aerial images are analyzed in conjunction with orthophotography and Lidar data to develop probability distributions for variables representing important erosion processes. Scans of photographic negatives were obtained for multiple photographs and several negatives were scanned with a calibrated photogrammetric scanner by Geomni (formerly Blue Skies Consulting) in Belen, New Mexico. Other negatives were provided by USDA or scanned on a high-resolution, large format media scanner in the United States National Archives Cartographic Research Room, College Park, Maryland (1939). An additional set of aerial photographs was obtained from the USGS EROS Data Center web site. For other years, photographic prints were available for aerial photography. Several prints were scanned with a large-format media scanner at 1200 dpi. Ground control points for selected digital images used in the stereoscopic analyses were identified in ArcMap. Ground control points were collected by matching features which could be identified in the aerial photography as well as the orthophotography or shaded relief imagery from the Lidar dataset. Such features included fence line intersections, road intersections, railroad bridges, and structures, depending on the age of the photography. Many old fence lines and abandoned roads remained visible in the shaded relief imagery while they were completely obscured by tree canopy in modern photography. Ground control elevations were sampled from the Lidar-based DEM or directly from point clouds. The image coordinates and matching three-dimensional projected coordinates of the ground control points were exported to IDL where parameters for a three-dimensional direct linear transformation (DLT) were obtained. Built-in singular value decomposition routines were used to solve an overdetermined set of equations for the eleven parameters of the DLT for a particular historical photograph. A solution to the DLT parameters was obtained for each historical aerial photograph in which features of interest were identified using the stereoscope. Ultimately, this process resulted in the identification of locations for features like gully heads using both the Lidar dataset and projection of the estimated location from the historical aerial photo under consideration back to the Lidar dataset. The distance between the estimate of the location from the historical image and that of the Lidar dataset was the estimated distance the feature has moved. This distance was then divided by the number of years between the Lidar dataset and the year of the historical aerial image of interest to get a rate of movement through time. This was done for several dozen points on each historical aerial image. There was inevitable uncertainty associated with the estimation of points from the historical aerial imagery. This includes estimation of changes in locations of features like gully heads through time. In this work a method to estimate the uncertainty associated with measurements of locations of features using the stereoscopic analysis is developed and applied. For features with locations in space that do not change (e.g., road intersections, etc.) the difference in measurements of locations from the Lidar dataset and that of the estimate from the historical aerial image provides an observation of the measurement error. These points used to estimate the measurement error are referred to as control points. For each of the historical aerial images, a different measurement error distribution is estimated using control points. These measurement errors are combined with the estimated distances from the stereoscopic analyses to produce distributions of gully retreat rates and hillslope erosion rates used along with output from Erosion Working Group (EWG) Landscape Evolution Models (LEMs) to provide statistical probability distributions for input into the PPA Model. Specifically, gully head retreat rate and gully widening rate were characterized using statistical probability distributions. The central limit theorem was then applied to develop rates that represent long-term averages that are used to inform a PPA Model built in GoldSim. The importance of these rates with respect to characterization of future risk can then assessed using sensitivity analyses. This provides an approach to determine if subsequent reduction of the uncertainty associated with the probability distributions of erosion rates will correspond to meaningful reduction in the estimation of future risk from the site. Ultimately this is a general method that utilizes historical aerial imagery to estimate rates of erosion with uncertainty. This information can be used to inform PPA models to project future risks from a given site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics

Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 4: Maintenance document (version 3.0)

The Maintenance Document Version 3.0 is a guide to the PAN AIR software system, a system which computes the subsonic or supersonic linear potential flow about a body of nearly arbitrary shape, using a higher order panel method. The document describes the overall system and each program module of the system. Sufficient detail is given for program maintenance, updating, and modification. It is assumed that the reader is familiar with programming and CRAY computer systems. The PAN AIR system was written in FORTRAN 4 language except for a few CAL language subroutines which exist in the PAN AIR library. Structured programming techniques were used to provide code documentation and maintainability. The operating systems accommodated are COS 1.11, COS 1.12, COS 1.13, and COS 1.14 on the CRAY 1S, 1M, and X-MP computing systems. The system is comprised of a data base management system, a program library, an execution control module, and nine separate FORTRAN technical modules. Each module calculates part of the posed PAN AIR problem. The data base manager is used to communicate between modules and within modules. The technical modules must be run in a prescribed fashion for each PAN AIR problem. In order to ease the problem of supplying the many JCL cards required to execute the modules, a set of CRAY procedures (PAPROCS) was created to automatically supply most of the JCL cards. Most of this document has not changed for Version 3.0. It now, however, strictly applies only to PAN AIR version 3.0. The major changes are: (1) additional sections covering the new FDP module (which calculates streamlines and offbody points); (2) a complete rewrite of the section on the MAG module; and (3) strict applicability to CRAY computing systems.

Purdon, David J.↗
DOE OSTI2022

Emergence of local scaling relations in adsorption energies on high-entropy alloys

Alloying has been proposed to circumvent scaling relations between the adsorption energies thus allowing for the complete optimization of multistep reactions. Herein the fidelity of scaling rules on high-entropy alloy (HEA) surfaces is assessed focusing on hydrogen-containing molecules, *AH x for A = C and N (x = 0, 1, 2, 3), A = S (x = 0, 1, 2) and A = O (x = 0, 1). Using an adsorbate- and site-specific deep learning model to rapidly compute the adsorption energies on CoMoFeNiCu HEA surfaces, the energies of *AH x and *A are shown to be linearly correlated if *A and *AH x have identical adsorption site symmetry. However, a local linear dependence emerges between the configuration-averaged adsorption energies irrespective of the site symmetry. Although these correlations represent a weaker form of the scaling relationships, they are sufficient to prohibit the optimization of multistep reactions. The underpinning of this behavior is twofold (1) the nearsightedness principle and (2) the narrow distribution of the adsorption energies around the mean-field value. While the nearsightedness is general for all electronic systems, the second criterion applies in HEAs with relatively strong reactive elements. The present findings strongly suggest that alloys may not generally enable the breaking of scaling relationships.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗
DOE OSTI2022

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗