SEARCH · Engineering Papers

Results for “balance function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,549 records · Page 86

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Seal Technology for Liquid Oxygen (LOX) Turbopumps

Two types of advanced seals for liquid oxygen (LOX) turbopumps were investigated. One was a spiral-groove face seal whose function is to seal high-pressure LOX at the impeller end of the turbopump. The other was a floating-ring, Rayleigh-step, helium buffered seal used to prevent LOX ingress to the turbine side of the unit. For each seal type, two sizes were investigated (50 and 20 mm). A turbine-driven test rig was designed and manufactured, and a test program was completed on the 50 mm floating-ring, Rayleigh-step, helium buffered seal. Significant results were: vaporization in the flow path could cause failure by overheating; therefore, the spiral-groove pumping portion of the seal that provides the fluid film must circulate fluid without disruption if vaporization occurs in the sealing dam. This is successfully accomplished by a pressure-balanced spiral-groove concept that is described. The spiral-groove configuration is affected by turbulence in the fluid film and pressure drops due to fluid inertia at sudden contractions. The net results of these effects are deep grooves, large operating films, and high power loss when compared against seals operating with laminar films. Turbulence and inertia are induced by the high-density and low-viscosity characteristics of LOX. The program clearly pointed out the need to consider system environmental factors such as thermal and centrifugal distortions and rotor vibrations in the seal design.

Shapiro, Wilbur↗

NASA NTRS2014

Understanding Differences in Chemistry Climate Model Projections of Stratospheric Ozone

Chemistry climate models (CCMs) are used to project future evolution of stratospheric ozone as concentrations of ozone-depleting substances (ODSs) decrease and greenhouse gases increase, cooling the stratosphere. CCM projections exhibit not only many common features but also a broad range of values for quantities such as year of ozone return to 1980 and global ozone level at the end of the 21st century. Multiple linear regression is applied to each of 14 CCMs to separate ozone response to ODS concentration change from that due to climate change. We show that the sensitivity of lower stratospheric ozone to chlorine change Delta Ozone/Delta inorganic chlorine is a near-linear function of partitioning of total inorganic chlorine into its reservoirs; both inorganic chlorine and its partitioning are largely controlled by lower stratospheric transport. CCMs with best performance on transport diagnostics agree with observations for chlorine reservoirs and produce similar ozone responses to chlorine change. After 2035, differences in Delta Ozone/Delta inorganic chlorine contribute little to the spread in CCM projections as the anthropogenic contribution to inorganic chlorine becomes unimportant. Differences among upper stratospheric ozone increases due to temperature decreases are explained by differences in ozone sensitivity to temperature change Delta Ozone/Delta T due to different contributions from various ozone loss processes, each with its own temperature dependence. Ozone decrease in the tropical lower stratosphere caused by a projected speedup in the Brewer-Dobson circulation may or may not be balanced by ozone increases in the middle- and high-latitude lower stratosphere and upper troposphere. This balance, or lack thereof, contributes most to the spread in late 21st century projections.

Stratospheric Transport↗

DOE OSTI2024

Spread spectrum time domain reflectometry (SSTDR) and frequency domain reflectometry (FDR) cable inspection using machine learning

Cables are initially qualified for nuclear power plant use for 40 years. As plants extend their operating license to 60 and 80 years, justification for continued cable use must shift to a condition-based approach since it is cost prohibitive to completely replace cables that are likely still capable of performing their design function. The Pacific Northwest National Laboratory (PNNL) Accelerated and Real Time Experimental Nodal Analysis (ARENA) cable motor test bed was used to test the response of a commercial spread spectrum time domain reflectometry (SSTDR) system, a laboratory instrument software-controlled SSTDR, and a vector network analyzer-based frequency domain reflectometry (FDR) system to various cable anomalies. The three instrument systems were able to interrogate cables over a range of frequency bandwidths that can be helpful for human data analysis. Data were subjected to supervised and unsupervised machine learning (ML) analyses to distinguish normal undamaged cable responses from anomalous cable responses. Both supervised and unsupervised ML approaches produced encouraging results with an undamaged/anomalous prediction accuracy from 0.69% to 0.87%. Recommendations for further development and field implementation include increased and more balanced sample sets particularly including more training data.

SSTDR, FDR, Reflectometry, Machine Learning, ARENA↗

NASA NTRS1980

Pioneer Venus Sounder Probe Solar Flux Radiometer

The Solar Flux Radiometer aboard the Pioneer Venus Sounder Probe operated successfully during its descent through the atmosphere of Venus. The instrument measured atmospheric radiance over the spectral range from 400 to 1800 nm as a function of altitude. Elevation and azimuthal measurements on the radiation field were made with five optical channels. Twelve filtered Si and Ge photovoltaic detectors were maintained near 30 C with a phase-change material. The detector output currents were processed with logarithmic transimpedance converters and digitized with an 11-bit A/D converter. Atmospheric sampling in both elevation and azimuth was done according to a Gaussian integration scheme. The serial output data averaged 20 bits/sec, including housekeeping (sync, spin period, sample timing and mode). The data were used to determine the deposition of solar energy in the atmosphere of Venus between 67 km and the surface along with upward and downward fluxes and radiances with an altitude resolution of several hundred meters. The results allow for more accurate modeling of the radiation balance of the atmosphere than previously possible.

Tomasko, M. G.↗

NASA NTRS1983

Creep testing of foil-gage metals at elevated temperature using an automated data acquisition system

A method is being developed to obtain creep data on foil gage metals at elevated temperatures using an automated data acquisition system in conjunction with a mechanically counter balanced extensometer. The automated system components include the Hewlett-Packard (HP) 9845A desktop computer, the HP 3455A digital voltmeter and the HP 3495A scanner. Software for test monitoring and data collection was developed; data manipulation, including curve plotting was done with a HP regression analysis software package. Initial creep tests were conducted on .003 in. thick foil specimens of Ti-6A1-4V at temperatures of 800 F and 1000 F and at stress levels of 25 ksi and 45 ksi. For comparison, duplicate tests were run on .049 in. thick specimens sheet of the same alloy. During testing, the furnace and specimen temperature, bridge voltage, strain and load output were automatically monitored and recorded at predetermined intervals. Using the HP regression analysis program, recorded strain output was plotted as a function of time. These resultant creep curves indicate that, under similar conditions of temperature and stress, foil gage specimens exhibit a higher creep rate than sheet specimens.

Blackburn, L. B.↗

NASA NTRS1996

MSC/NASTRAN DMAP Alter Used for Closed-Form Static Analysis With Inertia Relief and Displacement-Dependent Loads

Solving for the displacements of free-free coupled systems acted upon by static loads is a common task in the aerospace industry. Often, these problems are solved by static analysis with inertia relief. This technique allows for a free-free static analysis by balancing the applied loads with the inertia loads generated by the applied loads. For some engineering applications, the displacements of the free-free coupled system induce additional static loads. Hence, the applied loads are equal to the original loads plus the displacement-dependent loads. A launch vehicle being acted upon by an aerodynamic loading can have such applied loads. The final displacements of such systems are commonly determined with iterative solution techniques. Unfortunately, these techniques can be time consuming and labor intensive. Because the coupled system equations for free-free systems with displacement-dependent loads can be written in closed form, it is advantageous to solve for the displacements in this manner. Implementing closed-form equations in static analysis with inertia relief is analogous to implementing transfer functions in dynamic analysis. An MSC/NASTRAN (MacNeal-Schwendler Corporation/NASA Structural Analysis) DMAP (Direct Matrix Abstraction Program) Alter was used to include displacement-dependent loads in static analysis with inertia relief. It efficiently solved a common aerospace problem that typically has been solved with an iterative technique.

Source record↗

NASA NTRS2011

Current Sheets in the Heliosheath: Voyager 1, 2009

We identified all of the current sheets for which we have relatively complete and accurate magnetic field (B) data from Voyager 1 (V1) from days of year (DOYs) 1 to 331, 2009, which were obtained deep in the heliosheath between 108.5 and 111.8 AU. Three types of current sheets were found: (1) 15 proton boundary layers (PBLs), (2) 10 and 3 magnetic holes and magnetic humps, respectively, and (3) 3 sector boundaries. The magnetic field strength changes across PBL, and the profile B(t) is linearly related to the hyperbolic tangent function, but the direction of B does not change. For each of the three sector boundaries, B rotated in a plane normal to the minimum variance direction, and the component of B along the minimum variance direction was zero within the uncertainties, indicating that the sector boundaries were tangential discontinuities. The structure of the sector boundaries was not as simple as that for PBLs. The average thickness of magnetic holes and humps (approx.30 RL) was twice that of the PBLs (approx.15 RL). The average thickness of the current sheets associated with sector boundaries was close to the thickness of the PBLs. Our observations are consistent with the hypothesis that magnetic holes and humps are solitons, which are initiated by the mirror mode instability, and evolve by nonlinear kinetic plasma processes to pressure balanced structures maintained by magnetization currents and proton drift currents in the gradients of B.

Burlaga, L. F.↗

NASA NTRS2012

Intercomparison of MODIS Albedo Retrievals and In Situ Measurements Across the Global FLUXNET Network

Surface albedo is a key parameter in the Earth's energy balance since it affects the amount of solar radiation directly absorbed at the planet surface. Its variability in time and space can be globally retrieved through the use of remote sensing products. To evaluate and improve the quality of satellite retrievals, careful intercomparisons with in situ measurements of surface albedo are crucial. For this purpose we compared MODIS albedo retrievals with surface measurements taken at 53 FLUXNET sites that met strict conditions of land cover homogeneity. A good agreement between mean yearly values of satellite retrievals and in situ measurements was found (R(exp 2)= 0.82). The mismatch is correlated to the spatial heterogeneity of surface albedo, stressing the relevance of land cover homogeneity when comparing point to pixel data. When the seasonal patterns of MODIS albedo is considered for different plant functional types, the match with surface observation is extremely good at all forest sites. On the contrary, in non-forest sites satellite retrievals underestimate in situ measurements across the seasonal cycle. The mismatch observed at grasslands and croplands sites is likely due to the extreme fragmentation of these landscapes, as confirmed by geostatistical attributes derived from high resolution scenes.

Cescatti, Alessandro↗

NASA NTRS2024

NASA Engineering and Safety Center Lunar Rover Design Concepts Assessments

Mass is a significant risk to programs and projects as they transition from formulation to implementation, especially in larger human space systems where delivery mass or volume can be a constrained commodity. • Concepts developed without an adequate systems-engineering basis, including design and operations experience, may miss significant functionality and subsequent mass required for integration or operations. • Further uncertainty can be associated with not fully understanding design best practices and standards that drive mass, such as design for minimum risk or fault tolerance.1 • Finally, the appropriate systems engineering mass properties management rigor and technical discipline insight are required to set the mass baseline, including appropriate mass growth allowance (MGA) to ensure program success. Mass is a key quantity that should be constantly monitored by the systems engineer and the program/project management stakeholder to ensure mission compatibility throughout the project life cycle. In addition, mass is a key programmatic performance metric monitored by the NASA Chief Financial Officer for current and future program cost estimating. NASA Exploration Systems sought an assessment of reasonableness with respect to different potential rover concepts that balances mass needs and human-class cargo lander capabilities. Appropriate systems engineering mass properties management rigor and technical discipline insight were used, and are required to set the mass baseline, including appropriate MGA to ensure program success. Two independent Lunar rover concepts were evaluated, with a goal to understand concept credibility and the Lunar rover designs potential extensibility for Mars surface operations. A notional generic rover concept is shown in Figure 1.

Systems Engineering↗

NASA NTRS2007

Optimisation of the Future Routine Orbit for Mars Express

Mars Express (MEX), the first planetary mission of the European Space Agency (ESA), reached Mars on December 25th 2003. Since then it is performing routine operations. Its operational phase had to cover one Martian year, with the possibility of an extension for a second Martian year (i.e. until November 2007). The end of the mission extension is approaching but, given the good health of the payload instruments and the high science return of the mission, there is a strong will to achieve further extensions. Mars Express is also seen as an important asset, capable to provide relay functions for future Martian missions. The ESA Science Program Committee has recently approved a second extension of the MEX mission until May 2009 and even further extensions are possible. Mars Express has an eccentric quasi-polar orbit with a period of approximately 6.72 hours and a pericentre height of about 300 km. Science observations are mainly performed at pericentre (but not only). In addition the orbit has a resonance of 11 revolutions per 3 Martian days. This means that ground tracks corresponding to orbits separated by 11 revolutions are adjacent, such that a given area can be covered by the on-board camera without leaving gaps. The J2 effect of Mars causes a drift of both ascending node and argument of pericentre. The drift of argument of pericentre makes it possible to observe periodically all Mars latitudes from close distance. Illumination conditions at pericentre are influenced by both the drift of the argument of pericentre and the drift of ascending node, as well as by the rotation of Mars around the Sun. The original MEX routine orbit was optimized for the duration of the nominal mission and extension, such that it produced a balanced share of day-side observations (for the optical instruments) and night-side observations (for the radar). The orbit was thus not optimized for the time beyond the assumed extension. Indeed, the evolution of the ascending node and argument of pericentre would cause in the following years a drift of the pericentre towards night-side observation conditions, hence uninteresting for the optical instruments. In order to prevent this an optimisation process for the future routine orbit has taken place. The share between day-side and night-side observations can be controlled by adjusting the drift of argument of pericentre and ascending node. This can in particular be achieved by changing the semimajor axis, eccentricity and/or inclination. A change of inclination is inefficient compared to a change in semimajor axis and eccentricity, and has therefore been discarded. An in-plane maneuvre can be performed to change both semi-major axis and eccentricity, and thus the period of the orbit. Although an apocentre manoeuvre is cheaper in terms of deltaV, it would result in raising the pericentre height, which is unfavourable for close observations. Hence a pericentre manoeuvre is proposed, which will increase the apocentre height. A repeat cycle is still required to allow mapping areas with adjacent ground tracks, so the change of semimajor axis must result in a new resonance. Resonances 18:5, 25:7 and 7:2 have been considered as potential candidates. The resulting long term evolution of the observation conditions has been analysed. Finally it has been decided to perform a change of orbit to reach the 18:5. Another aspect of the optimisation process is the control of the ground track. The previous MEX reference trajectory included regular maneuvres at every apocentre in order to adjust the orbital period, such that the separation of the ground tracks would be optimal, regardless of the latitude of pericentre. The implementation of the actual delatVs on-board was done partly by optimizing the attitude of reaction-wheel desaturation activities. Despite of it, this strategy has a significant propellant cost, because it prevents to optimize reaction wheel de-saturation activities to minimize propellant consumption. Therore, with the aim at preserving propellant resources for a long time extension it has been agreed to stop the ground-track control. This requires now a more accurate science operation planning, with improved attitude pointing control. Finally, the approach to phase Mars Express to provide back-up relay functions for NASA Phoenix landing is explained. In the context of the routine trajectory optimisation a new requirement for close fly-bys at Phobos, with different observation geometries, has been specified. The approach to fulfill this requirement is explained.

Carranza, Manuel↗

NASA NTRS1981

The ionization state in a gas with a non-Maxwellian electron distribution

The inferred degree of ionization of a gas is often used in astrophysics as a diagnostic of the gas temperature. In the solar transition region and corona, in the outer atmospheres of cool stars, and in some portions of the interstellar medium), photoionization can be neglected, and the ionization state is fixed by the balance between ion-electron collisional ionization and dielectronic and/or radiative recombination. Under these conditions, higher degrees of ionization result from higher energy ion-electron collisions which are common in a high temperature gas. Actually, ionization occurs through collisions with electrons that have kinetic energies greater than the ionization potential of the given ion, and so the ionization rate depends on to the number of such high-energy electrons in the tail of the electron velocity distribution. High-velocity electrons move across large distances between effective coulomb collisions, and, in a strong temperature or density gradient, the tail can be overpopulated relative to Maxwell-Boltzmann distribution of equivalent energy density. Thus, the ionization rate can also be greatly increased. These effects for a parameterized form of the electron distribution function with an enhanced high-velocity tail, namely the kappa distribution are illustrated.

Owocki, S. P.↗

NASA NTRS2009

Regression Model Optimization for the Analysis of Experimental Data

A candidate math model search algorithm was developed at Ames Research Center that determines a recommended math model for the multivariate regression analysis of experimental data. The search algorithm is applicable to classical regression analysis problems as well as wind tunnel strain gage balance calibration analysis applications. The algorithm compares the predictive capability of different regression models using the standard deviation of the PRESS residuals of the responses as a search metric. This search metric is minimized during the search. Singular value decomposition is used during the search to reject math models that lead to a singular solution of the regression analysis problem. Two threshold dependent constraints are also applied. The first constraint rejects math models with insignificant terms. The second constraint rejects math models with near-linear dependencies between terms. The math term hierarchy rule may also be applied as an optional constraint during or after the candidate math model search. The final term selection of the recommended math model depends on the regressor and response values of the data set, the user s function class combination choice, the user s constraint selections, and the result of the search metric minimization. A frequently used regression analysis example from the literature is used to illustrate the application of the search algorithm to experimental data.

Ulbrich, N.↗

NASA NTRS2015

Lumped Parameter Models of the Central Nervous System for VIIP Research

INTRODUCTION: Current long-duration missions to the International Space Station and future exploration-class missions beyond low-Earth orbit, such as to Mars and asteroids, expose astronauts to increased risk of Visual Impairment and Intracranial Pressure (VIIP) syndrome [1]. It has been hypothesized that the headward shift of cerebral spinal fluid (CSF) and blood in microgravity may cause significant elevation of intracranial pressure (ICP), which in turn induces VIIP syndrome through biomechanical pathways [1, 2]. However, there is insufficient evidence to confirm this hypothesis. In this light, we are developing lumped-parameter models of fluid transport in the central nervous system (CNS) as a means to simulate the influence of microgravity on ICP. The CNS models will also be used in concert with the lumped parameter and finite element models of the eye described in the realted IWS abstracts submitted by Nelson et al., Feola et al. and Ethier et al. METHODS: We have developed a nine compartment CNS model (Figure 1) capable of both time-dependent and steady state fluid transport simulations, based on the works of Stevens et al. [3]. The breakdown of compartments within the model includes: vascular (3), CSF (2), brain (1) and extracranial (3). The boundary pressure in the Central Arteries [A] node is prescribed using an oscillating pressure function PA(t) simulating the carotid pulsatile pressure wave as developed by Linninger et al. [4]. For each time step, pressures are integrated through time using an adaptive-timestep 4th and 5th order Runga-Kutta solver. Once pressures are found, constitutive equations are used to solve for flowrates (Q) between each compartment. In addition to fluid flow between the different compartments, compliance (C) interactions between neighboring compartments are represented. We are also developing a second CNS model based on the works of Linninger et al. [4] which takes a more granular approach to represent the interactions of the intracranial and spinal compartments with the inclusion of arteries, arterioles, capillaries, venules, veins, venous sinus, and ventricles. The flow through the arteries, veins and CSF compartments are governed by continuity, momentum and distensibility balance equations. Furthermore, unlike the Stevens et al. approach, the Monro-Kellie doctrine of constant cranial volume and the bi-phasic nature of the brain parenchyma are implemented. These features appear to be more consistent with the physiologic and anatomical behavior of the CNS, and follow a modeling philosophy similar to the lumped parameter eye model that is intended to be integrated with the CNS model. However, Linninger’s approach has never been implemented to include hydrostatic gradient and microgravity simulation capabilities. Therefore, we aim at implement this modeling approach for spaceflight simulations and assess its overall applicability to VIIP research. OBJECTIVES: We will present verification and validation test results for both models, as well as head-to-head comparison to explore their strengths and limitations with respect to mathematical implementation and physiological significance for VIIP research. In doing so, we hope to provide some guidance to the HRP research community on how to appropriately leverage lumped parameter models for space biomedical research.

Vera, J.↗

DOE OSTI2022

A Mountaintop View Requires Minimal Sorting: A Faster Contour Tree Algorithm

Consider a scalar field f : M → R, where M is a triangulated simplicial mesh in R d . A level set, or contour, at value v is a connected component of f –1 (v). As v is changed, these contours change topology, merge into each other, or split. Contour trees are concise representations of f that track this contour behavior. The vertices of these trees are the critical points of f, where the gradient is zero. The edges represent changes in the topology of contours. It is a fundamental data structure in data analysis and visualization, and there is significant previous work (both theoretical and practical) on algorithms for constructing contour trees. Suppose M has n vertices, N facets, and t critical points. A classic result of Carr, Snoeyink, and Axen (2000) gives an algorithm that takes O(n log n+Nα(N)) time (where α(·) is the inverse Ackermann function). A further improvement to O(t log t + N) time was given by Chiang et al. All these algorithms involve a global sort of the critical points, a significant computational bottleneck. Unfortunately, lower bounds of Ω(t log t) also exist. We present the first algorithm that can avoid the global sort and has a refined time complexity that depends on the contour tree structure. Intuitively, if the tree is short and fat, we get significant improvements in running time. For a partition of the contour tree into a set of descending paths, P, our algorithm runs in O($\Sigma$ pϵP |p| log |p| + tα(t) + N). This is at most O(t log D + N), where D is the diameter of the contour tree. Moreover, it is O(tα(t) + N) for balanced trees, a significant improvement over the previous complexity. Our algorithm requires numerous ideas: partitioning the contour tree into join and split trees, a local growing procedure to iteratively build contour trees, and the use of heavy path decompositions for the time complexity analysis. There is a crucial use of a family of binomial heaps to maintain priorities, ensuring that any comparison made is between comparable nodes of the contour tree. We also prove lower bounds showing that the $\Sigma$ pϵP |p| log |p| complexity is inherent to computing contour trees.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Climate adaptation and sustainability in switchgrass: exploring plant-microbe-soil interactions across continental scale environmental gradients

Less carbon-intensive energy sources are needed to reduce greenhouse gas emissions and their predicted role in climate change. There is growing interest in the potential of biofuels for meeting this need. A critical question is whether large-scale biofuel production can be sustainable over the time scales needed to mitigate our carbon debt from fossil fuel consumption. The carbon balance and ultimately the sustainability of biofuel feedstock production is the result of complex climate-coupled interactions between carbon fixation, sequestration, and release through combustion. Similarly, the long-term productivity of biofuels depends on the environmental factors limiting plant growth. These factors are often related to soil resources which involve complex interactions at the plant-microbe-soil interface impacting their availability and cycling. Our collaborative project addressed sustainable switchgrass (Panicum virgatum) production by exploring Plant Systems, Plant-Microbiome Interactions, and Ecosystem Processes through the integrating lens of Multi-Scale Modeling. Our research was based on detailed characterization of genetically diverse switchgrass genotypes planted in common gardens across a continental latitudinal gradient. The underlying theme of our Plant Systems research was the use of locally adapted plant material to explore plant function, to understand the mechanistic basis of environmental interactions, and to discover the plant genes important for adaptation and sustainability in the face of climate change. Our Plant-Microbiome Interaction project characterized the microbial communities associated with switchgrass using genomic tools. Our Ecosystem Processes research focused on carbon cycle responses at the ecosystem level using stand level plantings. Finally, our Multi-Scale Modeling helped to define conditions of a sustainable biofuel system and identify key tradeoffs between genetic diversity, productivity, and ecosystem services. Genome-wide association analyses were used to identify alleles that contribute to successful establishment and biomass production across North America. Together, our work provided a baseline analyses of the potential of switchgrass as a biofuel feedstock. Our project resulted in a number of successful outcomes. First, we were successful in collecting switchgrass germplasm across the species range, propagating the material, and establishing common garden experiments across the species range. In collaboration with DOE JGI, we successfully assembled the first tetraploid switchgrass genome and published this resource with an analyses of the genetic basis local adaptation from our gardens (Lowry et al. 2019, Lovell et al. 2021). The gardens were used to characterize the genetic architecture for a number of important plant phenotypes. Our project also conducted extensive sampling and sequencing to characterize the bacterial and fungal associates of switchgrass roots and leaves. We showed that host genotype, location, and harvesting practices can play a role in microbiome assembly (Singer et al. 2019 & 2022, Van Wallendael et al. 2020 & 2022, Edwards et al. 2023). Our ecosystem processes work created baseline dataset of carbon and nutrient cycling in realistic stand plantings of switchgrass. Data from this experiment provided new insight into the role of plant traits, phenology, and local environments in ecosystem processes like soil respiration, net-ecosystem exchange, and dynamics of soil and plant nutrients (Ricketts et al. 2023). Finally, our crop modelling experiments help to characterize the sensitivity of common modeling frameworks to parameters, identify key limiters of productivity across large geographic scales, and leverage patterns of local adaptation in prediction. Ultimately, these studies help to identify critical plant-microbe-soil traits that may be manipulated, through breeding or agronomic management, to improve the sustainability of biofuel feedstocks.

09 BIOMASS FUELS↗

NASA NTRS1995

Upper Thermosphere Winds and Temperatures in the Geomagnetic Polar Cap: Solar Cycle, Geomagnetic Activity, and Interplanetary Magnetic Field Dependencies

Ground-based Fabry-Perot interferometers located at Thule, Greenland (76.5 deg. N, 69.0 deg. W, lambda = 86 deg.) and at Sondre Stromfjord, Greenland (67.0 deg. N, 50.9 deg. W, lambda = 74 deg.) have monitored the upper thermospheric (approx. 240-km altitude) neutral wind and temperature over the northern hemisphere geomagnetic polar cap since 1983 and 1985, respectively. The thermospheric observations are obtained by determining the Doppler characteristics of the (OI) 15,867-K (630.0-nm) emission of atomic oxygen. The instruments operate on a routine, automatic, (mostly) untended basis during the winter observing seasons, with data coverage limited only by cloud cover and (occasional) instrument failures. This unique database of geomagnetic polar cap measurements now extends over the complete range of solar activity. We present an analysis of the measurements made between 1985 (near solar minimum) and 1991 (near solar maximum), as part of a long-term study of geomagnetic polar cap thermospheric climatology. The measurements from a total of 902 nights of observations are compared with the predictions of two semiempirical models: the Vector Spherical Harmonic (VSH) model of Killeen et al. (1987) and the Horizontal Wind Model (HWM) of Hedin et al. (1991). The results are also analyzed using calculations of thermospheric momentum forcing terms from the Thermosphere-ionosphere General Circulation Model TGCM) of the National Center for Atmospheric Research (NCAR). The experimental results show that upper thermospheric winds in the geomagnetic polar cap have a fundamental diurnal character, with typical wind speeds of about 200 m/s at solar minimum, rising to up to about 800 m/s at solar maximum, depending on geomagnetic activity level. These winds generally blow in the antisunward direction, but are interrupted by episodes of modified wind velocity and altered direction often associated with changes in the orientation of the Interplanetary Magnetic Field (IMF). The central polar cap (greater than approx. 80 magnetic latitude) antisunward wind speed is found to be a strong function of both solar and geomagnetic activity. The polar cap temperatures show variations in both solar and geomagnetic activity, with temperatures near 800 K for low K(sub p) and F(sub 10.7) and greater than about 2000 K for high K(sub p) and F(sub 10.7). The observed temperatures are significantly greater than those predicted by the mass spectrometer/incoherent scatter model for high activity conditions. Theoretical analysis based on the NCAR TIGCM indicates that the antisunward upper thermospheric winds, driven by upstream ion drag, basically 'coast' across the polar cap. The relatively small changes in wind velocity and direction within the polar cap are induced by a combination of forcing terms of commensurate magnitude, including the nonlinear advection term, the Coriolis term, and the pressure gradient force term. The polar cap thennospheric thermal balance is dominated by horizontal advection, and adiabatic and thermal conduction terms.

Killeen, T. L.↗

DOE OSTI2022

Evaluated 238 U(n,f) Average Prompt Fission Neutron Multiplicities Including the CGMF Model

This report documents an evaluation of the average prompt fission neutron multiplicity, $\overline{v}_p$, of 238 U from 800 keV to MeV. This evaluation had to be re-done from “scratch” as the input to previous $\overline{v}_p$ evaluations, specifically ENDF/B-VIII.0, was not found. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. The new evaluated 238 U $\overline{v}_p$ based on only experimental data differs distinctly from ENDF/B-VIII.0 $\overline{v}_p$ from 2 to 4.5 MeV, and from 6 to 7 MeV, and is otherwise similar. The difference from 2 to 4.5 MeV stems from the fact that ENDF/B-VIII.0 was tweaked in this energy range to data of Frehaut, while two other, equally trustworthy, data sets would indicate an evaluated 238 U $\overline{v}_p$ that is up to 2% higher. Also, second chance fission in ENDF/B-VIII.0 was smoothed over from 6–7 MeV. Another major difference to ENDF/B-VIII.0 is that one of the evaluations presented here includes model information from the Hauser-Feshbach fission fragment decay code CGMF, while ENDF/B-VIII.0 is based purely on experimental data. CGMF links several fission quantities with each other; $\overline{v}_p$ is predicted by assumptions made on, e.g., pre-neutron emission yields as a function of mass, the total kinetic energy, or spin and parity of fission fragments. This allows to validate the new 238 U $\overline{v}_p$ by using CGMF parameters obtained from fitting to experimental 238 U $\overline{v}_p$ to predict yields as a function of mass, the average total kinetic energy, or the mean energy of the prompt fission neutron spectrum. These model-predicted values can then be compared to experimental and evaluated data. The model-predicted fission-observable values using evaluated parameters obtained here are reasonably close to experimental data for some observables, but are farther away from experimental data related to TKE observables. In addition to that, the evaluated 238 U(n,f) $\overline{v}_p$ shows similar deviations from ENDF/B-VIII.0 as for the evaluation with only experimental data. This difference is expected to lead to changes in simulated effective neutron multiplication factor, $k_{eff}$ of ICSBEP critical assemblies that are sensitive to 238 U in the fast range (BigTen, Flattop, Flattop-Pu). These changes in $k_{eff}$ need to be counter-balanced. Chi-Nu PFNS experimental data are expected to be released in the next few months that might lead to the needed changes in the PFNS. Until then, we hold off in benchmarking the new 238 U(n,f) $\overline{v}_p$ as well as submitting it to ENDF/B-VIII.1. Also, new high-precision 238 U $\overline{v}_p$ are expected to be measured by the CEA in the next two years that will shed further light on question on 238 U $\overline{v}_p$ from 2–4.5 and 6–7 MeV.

238U↗

NASA NTRS2012

Extension of the Time-Spectral Approach to Overset Solvers for Arbitrary Motion

Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.

Leffell, Joshua Isaac↗