Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source term calculation”

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 415 records · Page 23

Spacecraft camera image registration

A system for achieving spacecraft camera (1, 2) image registration comprises a portion external to the spacecraft and an image motion compensation system (IMCS) portion onboard the spacecraft. Within the IMCS, a computer (38) calculates an image registration compensation signal (60) which is sent to the scan control loops (84, 88, 94, 98) of the onboard cameras (1, 2). At the location external to the spacecraft, the long-term orbital and attitude perturbations on the spacecraft are modeled. Coefficients (K, A) from this model are periodically sent to the onboard computer (38) by means of a command unit (39). The coefficients (K, A) take into account observations of stars and landmarks made by the spacecraft cameras (1, 2) themselves. The computer (38) takes as inputs the updated coefficients (K, A) plus synchronization information indicating the mirror position (AZ, EL) of each of the spacecraft cameras (1, 2), operating mode, and starting and stopping status of the scan lines generated by these cameras (1, 2), and generates in response thereto the image registration compensation signal (60). The sources of periodic thermal errors on the spacecraft are discussed. The system is checked by calculating measurement residuals, the difference between the landmark and star locations predicted at the external location and the landmark and star locations as measured by the spacecraft cameras (1, 2).

Kamel, Ahmed A.↗

Low-frequency modes and nonbarotropic effects in pseudo-Newtonian accretion disks

A recently developed formalism is used to reexamine the question of the existence of hydrodynamical modes that pulsate with very low frequencies in the inner regions of accretion disks. The formalism is valid in an exact sense for the adiabatic pulsations of rotating Newtonian fluids that are generally nonbarotropic (such as those with 'nonadiabatic temperature gradients,' for example), and hence its application in the present context represents an improvement over previous analyses that are more approximate. The formalism is applied to thin non-self-gravitating disks, with the gravitational potential of the central source modified in the usual way in order to simulate relativistic effects. In the barotropic limit, the analyses indicate that in many cases nearly Keplerian disks exhibit nonaxisymmetric modes of pulsation that are trapped in the inner disk regions, with pulsation periods much longer than the dynamical timescale. These results are similar to those of earlier calculations that assume disks pulsate without changing the temperature distribution. A method is developed for including lowest order nonbarotropic effects. Previous analyses have been incapable of accurately treating the nonbarotropic regime. The application of the present method to the low-frequency modes reveals that, due to unexpected cancellations among terms, the nonbarotropic correction to the pusation frequency omega is only of order tilde-omega(sub BV exp 2) omega, where tilde-omega(sub BV) is the appropriate dimensionless Brunt-Vaisala frequency. This correction is much smaller than the expected correction of order tilde-omega(sub BV) Omega, where Omega is the rotation angular velocity. The important conclusion drawn from this is that nonbarotropic corrections are generally small and hence that low-frequency modes persist into the nonbarotropic regime. For disk temperatures appropriate to X-ray emission, the adiabatic frequencies of trapped modes are of the same order as the frequencies associated with quasi-periodic oscillations (QPOs) of certain low-mass X-ray binaries and other X-ray sources. It is suggested that further consideration be given to the possiblity that modulation of emission processes by these modes can explain at least some of the QPO phenomena.

Ipser, James R.↗

A comparison of Lorentz, planetary gravitational, and satellite gravitational resonances

We consider a charged dust grain whose orbital motion is dominated by a planet's point-source gravity, but perturbed by higher-order terms in the planet's gravity field as well as by the Lorentz force arising from an asymmetric planetary magnetic field. Perturbations to Keplerian orbits due to a nonspherical gravity field are expressed in the traditional way: in terms of a disturbing function which can be expanded in a series of spherical harmonics (W. M. Kaula, 1966). In order to calculate the electromagnetic perturbation, we first write the Lorentz force in terms of the orbital elements and then substitute it into Gauss' perturbation equations. We use our result to derive strengths of Lorentz resonances and elucidate their properties. In particular, we compare Lorentz resonances to two types of gravitational resonances: those arising from periodic tugs of a satellite and those due to the attraction of an arbitrarily shaped planet. We find that Lorentz resonances share numerous properties with their gravitational counterparts and show, using simple physical arguments, that several of these patterns are fundamental, applying not only to our expansions, but to all quantities expressed in terms of orbital elements. Some of these patterns have been previously called 'd'Alembert rules' for satellite resonances. Other similarities arise because, to first-order in the perturbing force, the three problems share an integral of the motion. Yet there are also differences; for example, first-order inclination resonances exist for perturbations arising from planetary gravity and from the Lorentz force, but not for those due to an orbiting satellite. Finally, we provide a heuristic treatment of a particle's orbital evolution under the influence of drag and resonant forces. Particles brought into mean-motion resonances experience either trapping or resonant 'jumps,' depending on the direction from which the resonance is approached. We show that this behavior does not depend on the details of the perturbing force but rather is fundamental to all mean-motion resonances.

Hamilton, Douglas P.↗

BISON Simulated and Experimental Fission Product Release Comparisons from Reradiated AGR-3/4 Compacts During High Temperature Heating Tests

The fuel performance modeling code BISON was used to predict the release of fission products iodine-131 (131I), xenon-133 (133Xe), and krypton-85 (85Kr) from four re-irradiated AGR-3/4 fuel compacts containing tristructural isotropic (TRISO) coated particles during high-temperature isothermal heating tests. The AGR-3/4 fuel compacts were irradiated in the Advanced Test Reactor (ATR) as part of the third and fourth series of planned experiments to support the Advanced Gas Reactor (AGR) Program. They were subsequently stored and re-irradiated in the Neutron Radiography (NRAD) reactor for approximately five days and then stored for another five to eight days before being subjected to isothermal heating tests in the Fuel Accident Condition Simulation (FACS furnace) for 200 to 300 hours at temperatures between 1000°C and 1600°C to evaluate fission product release at elevated temperatures. New nuclide-specific fission product source term models for the three nuclides of interest were developed using the reactor multiphysics code Griffin and implemented into BISON to support this work. The new source term models were incorporated into coupled compact- and particle-scale BISON simulations, which predict spatially- and temporally-resolved radionuclide generation, radioactive decay, transport, and release throughout the entire irradiation history, including the initial ATR irradiation, NRAD re-irradiations, FACS heating tests, and intermediate periods spent in storage. The experimentally measured fission product release from the heating tests were compared to modeling release predictions calculated by BISON to evaluate how well the code compares to experimental results. Overall, the experimental measured and BISON predicted comparative release results varied but generally agreed to within 5 particle equivalents. Comparative release results identified general observations to take into consideration to help refine future models and reduce uncertainties associated with both the measurement results and predictive results. This includes developing new uranium oxycarbide (UCO) specific kernel diffusivities for the three isotopes examined to more accurately reflect the material properties of the fuel form. Deriving new diffusivities will aid in producing a more informed BISON model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An Appraisal of the Performance and Characteristics of Summary Rays Calculated for the SALSA3D Traveltime Dataset

The SALSA3D tomographic model provides a crucial community resource for improving the quality (in terms of both accuracy and precision) of predictions of the traveltimes of seismic waves and therefore improving our ability to locate anthropogenic or natural seismic events. Constructing the requisite tomographic model requires addressing the challenges implied by a massive and growing dataset of traveltime measurements. This study explores one approach to tackle this challenge: the use of summary rays, which average traveltime measurements from sources within evenly spaced cells, thereby eliminating redundant data.

58 GEOSCIENCES↗

The specific ( α ,n) production rate for 234 U in UF 6

Natural, depleted, and low enriched uranium hexafluoride is manufactured and transported on an industrial scale to support commercial nuclear power plant operation. One-way material production and flow is verified by nuclear safeguards inspectorates is by neutron counting. The dominant neutron source term is ( α ,n) production due to 234 U which gets enriched/depleted along with 235 U. We describe recent efforts to accurately determine the ( α ,n)/s/g of 234 U in UF 6 . This piece of nuclear data is of direct value to the nuclear safeguards community in the science-based verification and assay of bulk uranium hexafluoride. It also provides a benchmark on codes that calculate yield curves for various compounds from cross section data, and hence, provides a scale or normalization constraint on evaluated basic cross section data. Finally, in turn, evaluated basic data allows the yield from other materials to be calculated; for instance, moist uranyl fluoride, which is of interest in holdup measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Direct Determination of Absolute Absorption Cross Sections at the L-Edge of Dilute Mn Complexes in Solution Using a Transmission Flatjet

The 3d transition metals play a pivotal role in many charge transfer processes in catalysis and biology. X-ray absorption spectroscopy at the L-edge of metal sites probes metal 2p–3d excitations, providing key access to their valence electronic structure, which is crucial for understanding these processes. In this paper, we report L-edge absorption spectra of Mn II (acac) 2 and Mn III (acac) 3 complexes in solution, utilizing a liquid flatjet for X-ray absorption spectroscopy in transmission mode. With this, we derive absolute absorption cross-sections for the L-edge transitions with peak magnitudes as large as 12 and 9 Mb for Mn II (acac) 2 and Mn III (acac) 3 , respectively. We provide insight into the electronic structure with ab initio restricted active space calculations of these L-edge transitions, reproducing the experimental spectra with excellent agreement in terms of shapes, relative energies, and relative intensities for the two complexes. Crystal field multiplet theory is used to assign spectral features in terms of the electronic structure. Comparison to charge transfer multiplet calculations reveals the importance of charge transfer in the core-excited final states. Finally, on the basis of our experimental observations, we extrapolate the feasibility of 3d transition metal L-edge absorption spectroscopy using the liquid flatjet approach in probing highly dilute biological solution samples and possible extensions to table-top soft X-ray sources.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predictive onboard flow control for packet switching satellites

We outline two alternate approaches to predicting the onset of congestion in a packet switching satellite, and argue that predictive, rather than reactive, flow control is necessary for the efficient operation of such a system. The first method discussed is based on standard, statistical techniques which are used to periodically calculate a probability of near-term congestion based on arrival rate statistics. If this probability exceeds a present threshold, the satellite would transmit a rate-reduction signal to all active ground stations. The second method discussed would utilize a neural network to periodically predict the occurrence of buffer overflow based on input data which would include, in addition to arrival rates, the distributions of packet lengths, source addresses, and destination addresses.

Bobinsky, Eric A.↗

Predictive onboard flow control in packet switching satellites

We outline two alternate approaches to predicting the onset of congestion in a packet switching satellite, and argue that predictive, rather than reactive, flow control is necessary for the efficient operation of such a system. The first method discussed is based on standard, statistical techniques which are used to periodically calculate a probability of near-term congestion based on arrival rate statistics. If this probability exceeds a present threshold, the satellite would transmit a rate-reduction signal to all active ground stations. The second method discussed would utilize a neural network to periodically predict the occurrence of buffer overflow based on input data which would include, in addition to arrival rates, the distributions of packet lengths, source addresses, and destination addresses.

Bobinsky, E. A.↗

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Melting behavior and phase relations of lunar samples

Cooling rate studies of 12002 were conducted and the results interpreted in terms of the crystallization history of this rock and certain other picritic Apollo 12 samples. Calculations of liquid densities and viscosities during crystallization, crystal settling velocities, and heat loss by the parent rock body are discussed, as are petrographic studies of other Apollo 12 samples. The process of magmatic differentiation that must have accompanied the early melting and chemical fractionation of the moon's outer layers was investigated. The source of regions of both high- and low-titanium mare basalts were also studied.

Hays, J. F.↗

Multistage Simulations of the GE90 Turbine

The average passage approach has been used to analyze three multistage configurations of the GE90 turbine. These are a high pressure turbine rig, a low pressure turbine rig and a full turbine configuration comprising 18 blade rows of the GE90 engine at takeoff conditions. Cooling flows in the high pressure turbine have been simulated using source terms. This is the first time a dual-spool cooled turbine has been analyzed in 3D using a multistage approach. There is good agreement between the simulations and experimental results. Multistage and component interaction effects are also presented. The parallel efficiency of the code is excellent at 87.3% using 121 processors on an SGI Origin for the 18 blade row configuration. The accuracy and efficiency of the calculation now allow it to be effectively used in a design environment so that multistage effects can be accounted for in turbine design.

Turner, Mark G.↗

Refined interpretation of electron temperature response to neutral beam injection at DIII-D

Accurate particle and power deposition profiles of neutral beam injection (NBI) are essential to transport studies, and that information is usually acquired through Monte Carlo simulations with a given collisional model. The deposition process of the energetic beam particles leads to the informative electron temperature (T e ) evolution trajectory, which can be captured by electron cyclotron emission (ECE) system due to its good spatial and temporal resolution. Previously, some work has been done to interpret the T e responses to the pulsed NBI as a linear heating source with Fourier-based techniques, although that approach fell short when the fast ion slowing-down time becomes significant (~100 ms). It has been observed in DIII-D that the modulated NBI pulses (10–50 Hz) reduce local core T e values ~0.1 keV through cold electron dilution in high-T e (>2 keV) plasmas alongside accumulative heating. Here, a novel approach to interpret the T e response to NBI was developed by linearizing and modeling the detailed T e evolution trajectory using coherently averaged ECE data based on the different time scales of the terms in the local power and particle balance equations. The technique does not require absolute calibrations of ECE and is independent of collisional models. The resulting beam deposition profiles show good consistency and reasonable agreement with Monte Carlo calculations based on the atomic data from the Atomic Data and Analysis Structure (ADAS). Local electron density response measured by Thomson scattering (TS) also suggests the same features when the beam pulse is large enough for that diagnostic to resolve. The remaining discrepancies are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

L Band Brightness Temperature from Forest: Comparison of Approximate Techniques

In this paper, three approximate physical microwave radiometry models have been used to calculate brightness temperatures from a forest canopy at L-band. These models are (1) tau-omega model (zero order scattering approximation to radiative transfer equations), (2) successive order of scattering model up to first order (first order scattering approximation to the radiative transfer equations), and (3) Peak technique utilizing the active solution obtained from the Distorted Born Approximation (DBA). These models are physically-based and treat vegetation as a layer of discrete scatterers over a rough surface. Vegetation components within the canopy are represented by canonical shapes such as dielectric discs and cylinders. The tau-omega model is based on a zero-order solution to the radiative transfer (RT) equations. The model ignores scattering except for the effect of the scatterers in the attenuation of the emission through the vegetation. Application of the tau-omega model to data acquired during airborne and ground-based campaigns over the years has solidified scientific understanding of microwave interactions with different landscapes. In particular, shrubland, grasslands, agricultural crops, and light to moderate vegetation have been investigated. Its applicability to areas with a significant tree fraction is unknown. The first order scattering model is based on an iterative solution of the RT equation up to the first order. The first order solution is obtained by substituting the zeroth-order solution into the scattering source term and then solving the resulting radiative transfer equations. This formulation adds a new scattering term to the tau-omega model. It represents emission by particles in the layer and emission by the ground that is scattered once by particles in the layer. The resulting model represents an improvement over the standard zero-order solution (the tau-omega model) since it accounts for the scattered vegetation and ground radiation that can have a pronounced effect on the observed brightness temperature. The third model is based on the Peake formulation in conjunction with the DBA. The procedure for calculation of forest emission is accomplished by first calculating the bistatic scattering cross section for each type of scatterer, then by using the DBA to calculate specular albedo of the ground and the diffused albedo of the layer. Once the albedos are determined, Peake s principle relating active and passive problems can be used to determine the effective emissivity of the forest layer.

Kurum, Mehmet↗

The De-Icing Comparison Experiment (D-ICE): a study of broadband radiometric measurements under icing conditions in the Arctic

Abstract. Surface-based measurements of broadband shortwave (solar) and longwave (infrared) radiative fluxes using thermopile radiometers are made regularly around the globe for scientific and operational environmental monitoring. The occurrence of ice on sensor windows in cold environments – whether snow, rime, or frost – is a common problem that is difficult to prevent as well as difficult to correct in post-processing. The Baseline Surface Radiation Network (BSRN) community recognizes radiometer icing as a major outstanding measurement uncertainty. Towards constraining this uncertainty, the De-Icing Comparison Experiment (D-ICE) was carried out at the NOAA Atmospheric Baseline Observatory in Utqiaġvik (formerly Barrow), Alaska, from August 2017 to July 2018. The purpose of D-ICE was to evaluate existing ventilation and heating technologies developed to mitigate radiometer icing. D-ICE consisted of 20 pyranometers and 5 pyrgeometers operating in various ventilator housings alongside operational systems that are part of NOAA's Barrow BSRN station and the US Department of Energy Atmospheric Radiation Measurement (ARM) program North Slope of Alaska and Oliktok Point observatories. To detect icing, radiometers were monitored continuously using cameras, with a total of more than 1 million images of radiometer domes archived. Ventilator and ventilator–heater performance overall was skillful with the average of the systems mitigating ice formation 77 % (many >90 %) of the time during which icing conditions were present. Ventilators without heating elements were also effective and capable of providing heat through roughly equal contributions of waste energy from the ventilator fan and adiabatic heating downstream of the fan. This provided ∼0.6 ∘C of warming, enough to subsaturate the air up to a relative humidity (with respect to ice) of ∼105 %. Because the mitigation technologies performed well, a near complete record of verified ice-free radiometric fluxes was assembled for the duration of the campaign. This well-characterized data set is suitable for model evaluation, in particular for the Year of Polar Prediction (YOPP) first Special Observing Period (SOP1). We used the data set to calculate short- and long-term biases in iced sensors, finding that biases can be up to +60 W m−2 (longwave) and −211 to +188 W m−2 (shortwave). However, because of the frequency of icing, mitigation of ice by ventilators, cloud conditions, and the timing of icing relative to available sunlight, the biases in the monthly means were generally less than the aggregate uncertainty attributed to other conventional sources in both the shortwave and longwave.

54 ENVIRONMENTAL SCIENCES↗

Compton and synchrotron processes in spherically-symmetric non-thermal sources

A framework is developed for the accurate calculation of Compton and synchrotron processes in a special class of spherically-symmetric sources. The models considered take a distribution in magnetic field values and high-energy electron density that varies smoothly from central values to values at the outer boundary of the source. Further, the assumption is made that the magnetic field has local disorder (isotropy) and the high-energy electron distribution has local isotropy. Then, in terms of parameters of the source, the synchrotron brightness distribution, total synchrotron flux, Compton-synchrotron flux, and synchrotron self-absorption turnover are all computed accurately. Results are given with a view toward application to the analysis of compact non-thermal sources.

Gould, R. J.↗

Area estimation of environmental phenomena from NOAA-n satellite data

A technique for documenting changes in size of NOAA-n pixels in order to calibrate the data for use in performing area calculations is described. Based on Earth-satellite geometry, a function for calculating the effective pixel size, measured in terms of ground area, on any given pixel was derived. The equation is an application of the law of sines plus an arclength formula. Effective pixel dimensions for NOAA 6 and 7 satellites for all pixels between nadir and the extreme view angles are presented. The NOAA 6 data were used to estimate the areas of several lakes, with an accuracy within 5%. Sources of error are discussed.

Tappan, G.↗

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗