Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantification of LEU Holdup using gamma ray imaging and inverse transport solver

Holdup is the residual amount of special nuclear material (SNM) remaining in a processing facility after the bulk materials have been cleaned out. In commercial uranium processing facilities, quantification of holdup is a major challenge because of the highly variable shapes and sizes of the deposits. Any method that attempts to generalize and calibrate deposit shapes in order to quantify holdup will be prone to high uncertainties. Uncertainties on the order of ±50% are typical in holdup results. In international safeguards applications, a ±50% uncertainty can result in a large amount of material unaccounted for (MUF) thereby increasing the difficulty of detecting material diversion and facility misuse. An imaging-based methodology has been developed with the objective of significantly reducing this uncertainty by using the true deposit shape, instead of relying on oversimplified geometric assumptions. The project is a collaboration between ORNL, Y-12, and the University of Tennessee, Knoxville, TN. Uranium sources of known masses were measured using the Germanium Gamma-ray Imager (GeGI), a high-resolution imaging spectrometer, creating a pixelated map for each spectral bin. Two different gamma imaging methods are employed in this work: coded aperture imaging and Compton imaging. A validated MonteCarlo model of the detector has been developed using the GEANT4 code for determining the intrinsic response of the detector, its enclosure, and the coded aperture mask. An inverse transport solver based on the Markov Chain Monte-Carlo approach known as Differential Evolution Adaptive Metropolis (DREAM) is employed to use the measurement data from the image pixels (coded aperture or Compton) to solve for the mass of 235 U in the deposit. A reliable method based on the DREAM solver has been developed to flag the infinite thickness condition of a uranium deposit. The project team is working towards improving the image reconstruction for Compton imaging so that a better localization of the source can be achieved. Besides treating the coded aperture and Compton imaging methods independently, the project is also evaluating a combined method that uses the Compton scatter data from a coded aperture measurement. GEANT4 simulations are being performed to evaluate the combined approach. The impact on the DREAM optimization as the source thickness progressively approaches infinite thickness is being evaluated. A number of uranium sources available at ORNL have been measured, and the DREAM results have been tested and validated for the coded aperture imaging. A similar effort will be carried out to validate the Compton based method once the development of algorithms for better localization are complete. The imaging based quantification is very amenable to unattended monitoring of holdup accumulation at key measurement points. A proof of concept measurement has been completed to demonstrate this capability The current work used the high energy resolution imager GeGI. However, the approach and methodologies are applicable to other imagers such as the cadmium zin telluride (CZT) based imager manufactured by H3D, Inc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Spectroscopic Analysis of Milky Way Outer Halo Satellites: Aquarius II and Boötes II

In this paper, we present a chemical and kinematic analysis of two ultrafaint dwarf galaxies (UFDs), Aquarius II (Aqu II) and Boötes II (Boo II), using Magellan/IMACS spectroscopy. We present the largest sample of member stars for Boo II (12), and the largest sample of red giant branch members with metallicity measurements for Aqu II (eight). In both UFDs, over 80% of targets selected based on Gaia proper motions turned out to be spectroscopic members. In order to maximize the accuracy of stellar kinematic measurements, we remove the identified binary stars and RR Lyrae variables. For Aqu II, we measure a systemic velocity of -65.3 ± 1.8 km s-1 and a metallicity of [Fe/H] = $-2.57^{+0.17}_{-0.17}$. When compared with previous measurements, these values display a ~6 km s -1 difference in radial velocity and a decrease of 0.27 dex in metallicity. Similarly for Boo II, we measure a systemic velocity of $-130.4^{+1.4}_{-1.1}$ km s -1 , more than 10 km s -1 different from the literature, a metallicity almost 1 dex smaller at [Fe/H] = $-2.71^{+0.11}_{-0.10}$, and a velocity dispersion 3 times smaller at σ v hel = $2.9^{+1.6}_{-1.2}$km s -1 . Additionally, we derive systemic proper-motion parameters and model the orbits of both UFDs. Finally, we highlight the extremely dark-matter-dominated nature of Aqu II and compute the J-factor for both galaxies to aid searches of dark matter annihilation. Despite the small size and close proximity of Boo II, it is an intermediate target for the indirect detection of dark matter annihilation due to its low-velocity dispersion and corresponding low dark matter density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Constraining the Twomey effect from satellite observations: issues and perspectives

Abstract. The Twomey effect describes the radiative forcing associated with a change in cloud albedo due to an increase in anthropogenic aerosol emissions. It is driven by the perturbation in cloud droplet number concentration (ΔNd, ant) in liquid-water clouds and is currently understood to exert a cooling effect on climate. The Twomey effect is the key driver in the effective radiative forcing due to aerosol–cloud interactions, but rapid adjustments also contribute. These adjustments are essentially the responses of cloud fraction and liquid water path to ΔNd, ant and thus scale approximately with it. While the fundamental physics of the influence of added aerosol particles on the droplet concentration (Nd) is well described by established theory at the particle scale (micrometres), how this relationship is expressed at the large-scale (hundreds of kilometres) perturbation, ΔNd, ant, remains uncertain. The discrepancy between process understanding at particle scale and insufficient quantification at the climate-relevant large scale is caused by co-variability of aerosol particles and updraught velocity and by droplet sink processes. These operate at scales on the order of tens of metres at which only localised observations are available and at which no approach yet exists to quantify the anthropogenic perturbation. Different atmospheric models suggest diverse magnitudes of the Twomey effect even when applying the same anthropogenic aerosol emission perturbation. Thus, observational data are needed to quantify and constrain the Twomey effect. At the global scale, this means satellite data. There are four key uncertainties in determining ΔNd, ant, namely the quantification of (i) the cloud-active aerosol – the cloud condensation nuclei (CCN) concentrations at or above cloud base, (ii) Nd, (iii) the statistical approach for inferring the sensitivity of Nd to aerosol particles from the satellite data and (iv) uncertainty in the anthropogenic perturbation to CCN concentrations, which is not easily accessible from observational data. This review discusses deficiencies of current approaches for the different aspects of the problem and proposes several ways forward: in terms of CCN, retrievals of optical quantities such as aerosol optical depth suffer from a lack of vertical resolution, size and hygroscopicity information, non-direct relation to the concentration of aerosols, difficulty to quantify it within or below clouds, and the problem of insufficient sensitivity at low concentrations, in addition to retrieval errors. A future path forward can include utilising co-located polarimeter and lidar instruments, ideally including high-spectral-resolution lidar capability at two wavelengths to maximise vertically resolved size distribution information content. In terms of Nd, a key problem is the lack of operational retrievals of this quantity and the inaccuracy of the retrieval especially in broken-cloud regimes. As for the Nd-to-CCN sensitivity, key issues are the updraught distributions and the role of Nd sink processes, for which empirical assessments for specific cloud regimes are currently the best solutions. These considerations point to the conclusion that past studies using existing approaches have likely underestimated the true sensitivity and, thus, the radiative forcing due to the Twomey effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Influence of simple terrain on the spatial variability of a low-level jet and wind farm performance in the AWAKEN field campaign

In wind energy research, scientific challenges are often associated with complex terrain sites, where orography, vegetation, and buildings disrupt flow uniformity. However, even sites characterized as simple terrain can exhibit significant spatial variability in wind speed, particularly during stable boundary layers (SBLs) and low-level jets (LLJs). This study investigates these terrain interactions using both simulations and observations from the American WAKe ExperimeNt (AWAKEN). We employ a multiscale Weather Research and Forecasting (WRF) model simulation, integrating mesoscale forcing in the coarse domains and representing three rows of turbines from the King Plains wind farm as generalized actuator disks (GAD) in the large-eddy simulation (LES) domains. During a nocturnal LLJ event on 3 April 2023, the downstream, wake-affected turbine rows outperformed the upstream, unwaked row by 25 %–51 %. This counterintuitive result arises from terrain-induced streamwise variations in hub-height wind speed of approximately 4 m s −1 over 5 km – equivalent to ∼50 % of the upstream reference speed. This enhancement outweighs the wake-induced reduction in mean wind speed (∼12 %) and global blockage effects reported in the literature (∼1 %–3.4 %). The multiscale simulations capture the intra-farm spatial variability in power performance observed in SCADA data. Terrain-induced vertical displacement of the LLJ, coupled with large wind shear below the jet maximum, drives the substantial streamwise acceleration within the wind farm. These findings underscore the importance of accounting for spatial variability related to terrain, even in simple landscapes, particularly during LLJ conditions. Incorporating such effects into reduced-order modeling frameworks for wind farm design and control could significantly enhance their effectiveness.

17 WIND ENERGY↗
NASA NTRS1971

The three-dimensional steady radial expansion of a viscous gas from a sonic source into a vacuum.

The three-dimensional steady radial expansion of a viscous, heat-conducting, compressible fluid from a spherical sonic source into a vacuum is analyzed using the Navier-Stokes equations as a basis. It is assumed that the model fluid is a perfect gas having constant specific heats, a constant Prandtl number of order unity, and viscosity coefficients varying as a power of the absolute temperature. Limiting forms for the flow variable solutions are studied for the Reynolds number based on the sonic source conditions, going to infinity and the Newtonian parameter both fixed and going to zero. For the case of the viscosity-temperature exponent between .5 and 1, it is shown that the velocity as well as the pressure approach zero as the radial distance goes to infinity. The formulations of the distinct regions that span the domain extending from the sonic source to the vacuum are presented.

Bush, W. B.↗
NASA NTRS1975

The geometric approach to sets of ordinary differential equations and Hamiltonian dynamics

The calculus of differential forms is used to discuss the local integration theory of a general set of autonomous first order ordinary differential equations. Geometrically, such a set is a vector field V in the space of dependent variables. Integration consists of seeking associated geometric structures invariant along V: scalar fields, forms, vectors, and integrals over subspaces. It is shown that to any field V can be associated a Hamiltonian structure of forms if, when dealing with an odd number of dependent variables, an arbitrary equation of constraint is also added. Families of integral invariants are an immediate consequence. Poisson brackets are isomorphic to Lie products of associated CT-generating vector fields. Hamilton's variational principle follows from the fact that the maximal regular integral manifolds of a closed set of forms must include the characteristics of the set.

Estabrook, F. B.↗
NASA NTRS1976

A finite difference method of solving anisotropic scattering problems

A new method of solving radiative transfer problems is described including a comparison of its speed with that of the doubling method, and a discussion of its accuracy and suitability for computations involving variable optical properties. The method uses a discretization in angle to produce a coupled set of first-order differential equations which are integrated between discrete depth points to produce a set of recursion relations for symmetric and anti-symmetric angular sums of the radiation field at alternate depth points. The formulation given here includes depth-dependent anisotropic scattering, absorption, and internal sources, and allows arbitrary combinations of specular and non-Lambertian diffuse reflection at either or both boundaries. Numerical tests of the method show that it can return accurate emergent intensities even for large optical depths. The method is also shown to conserve flux to machine accuracy in conservative atmospheres

Barkstrom, B. R.↗
NASA NTRS1977

Large-scale V/STOL testing

Several facets of large-scale testing of V/STOL aircraft configurations are discussed with particular emphasis on test experience in the Ames 40- by 80-foot wind tunnel. Examples of powered-lift test programs are presented in order to illustrate tradeoffs confronting the planner of V/STOL test programs. It is indicated that large-scale V/STOL wind-tunnel testing can sometimes compete with small-scale testing in the effort required (overall test time) and program costs because of the possibility of conducting a number of different tests with a single large-scale model where several small-scale models would be required. The benefits of both high- and full-scale Reynolds numbers, more detailed configuration simulation, and number and type of onboard measurements increase rapidly with scale. Planning must be more detailed at large scale in order to balance the trade-offs between the increased costs, as number of measurements and model configuration variables increase and the benefits of larger amounts of information coming out of one test.

Koenig, D. G.↗
NASA NTRS1978

Computerized power supply analysis: State equation generation and terminal models

To aid engineers that design power supply systems two analysis tools that can be used with the state equation analysis package were developed. These tools include integration routines that start with the description of a power supply in state equation form and yield analytical results. The first tool uses a computer program that works with the SUPER SCEPTRE circuit analysis program and prints the state equation for an electrical network. The state equations developed automatically by the computer program are used to develop an algorithm for reducing the number of state variables required to describe an electrical network. In this way a second tool is obtained in which the order of the network is reduced and a simpler terminal model is obtained.

Garrett, S. J.↗
NASA NTRS1978

Methods used for Space Shuttle SRB thrust shape design

Optimization of the Space Shuttle trajectory is discussed with reference to the low acceleration profile required for Shuttle missions. Static tests of the nominal flight curve are described in terms of impulse requirements, vacuum thrust, and burn time. Attention is given to BARF (Burning Anomaly Rate Factor), and it is noted that mandrel fabrication is intended to include the flexibility to counter BARF, should it occur. Test results are presented in which both BARF and specific impulse are considered as independent variables. It was found that no erosive burning occurred, BARF did not occur, specific impulse was on the order of 265 sec, and flow anomalies in the star region produced head-to-aft stagnation pressure drops in excess of theoretical predictions. In other areas, good agreement is noted between theoretical prediction and empirical data.

Baker, J.↗
NASA NTRS1979

The rotational motion of an earth orbiting gyroscope according to the Einstein theory of general relativity

The classical Poisson equations of rotational motion are used to study the attitude motions of an earth orbiting, rapidly spinning gyroscope perturbed by the effects of general relativity (Einstein theory). The center of mass of the gyroscope is assumed to move about a rotating oblate earth in an evolving elliptic orbit which includes all first-order oblateness effects produced by the earth. A method of averaging is used to obtain a transformation of variables, for the nonresonance case, which significantly simplifies the Poisson differential equations of motion of the gyroscope. Long-term solutions are obtained by an exact analytical integration of the simplified transformed equations. These solutions may be used to predict both the orientation of the gyroscope and the motion of its rotational angular momentum vector as viewed from its center of mass. The results are valid for all eccentricities and all inclinations not near the critical inclination.

Hoots, F. R.↗
NASA NTRS1981

The ultraviolet variability of the T Tauri star RW Aurigae

Between 1978 and 1979 the visible brightness of RW Aurigae increased by 0.9 mag. During this time (1) CIV and SiIV increased by factors of 2 to 4 while the lower ionization lines remained unchanged; (2) the fluorescent OI line increased by a factor of 8; (3) the shell spectrum changed from emission to absorption; and (4) the ultraviolet continuum brightened by 2.3 mag. On a time scale of a week the continuum varied by as much as 0.8 mag., but the MgII emission lines showed no variability over 10 percent. An active chromosphere, transition region, and envelope cooled by mass loss are hypothesized in order to explain the ultraviolet observations of RW Aur.

Imhoff, C. L.↗
NASA NTRS1984

Integration becomes the name of the OTV game

Since 1981, NASA has been supporting an Orbit Transfer Vehicle (OTV) propulsion technology development program which encompasses the efforts of three major engine manufacturers. A thrust variability ratio of 30:1 has been stipulated for the baseline 520 lbf-sec/lbm engines, in order to yield the versatility needed for low acceleration missions, orbit-transfer missions, and aeromaneuvering tasks. These goals may not be reachable either individually or collectively; the program supports generic technology benefitting all three engine concepts, as well as concept-specific technology. Engine/vehicle integration will determine the final configuration. The three engines under study differ as to turbomachinery types and operating speeds.

Garland, S.↗
NASA NTRS1988

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitary time, G(bar a, bar t), and that of the time when a crack reaches an arbitary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. TAking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗
NASA NTRS1989

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitrary time, G(bar a, bar t), and that of the time when a crack reaches an arbitrary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. Taking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗
NASA NTRS1989

A two dimensional power spectral estimate for some nonstationary processes

A two dimensional estimate for the power spectral density of a nonstationary process is being developed. The estimate will be applied to helicopter noise data which is clearly nonstationary. The acoustic pressure from the isolated main rotor and isolated tail rotor is known to be periodically correlated (PC) and the combined noise from the main and tail rotors is assumed to be correlation autoregressive (CAR). The results of this nonstationary analysis will be compared with the current method of assuming that the data is stationary and analyzing it as such. Another method of analysis is to introduce a random phase shift into the data as shown by Papoulis to produce a time history which can then be accurately modeled as stationary. This method will also be investigated for the helicopter data. A method used to determine the period of a PC process when the period is not know is discussed. The period of a PC process must be known in order to produce an accurate spectral representation for the process. The spectral estimate is developed. The bias and variability of the estimate are also discussed. Finally, the current method for analyzing nonstationary data is compared to that of using a two dimensional spectral representation. In addition, the method of phase shifting the data is examined.

Smith, Gregory L.↗
NASA NTRS1990

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗
NASA NTRS1991

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗