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 181 records · Page 10

Comparison between pressure gradient method and MAC method on high Re calculation

A cavity flow driven by shear and buoyancy forces is used as a test problem in the application of a nonstaggered pressure gradient (PG) method in solving the two-dimensional incompressible Navier-Stokes equations. Twelve finite differencing schemes are used to solve the cavity flow problem. The schemes consist of various combinations of grid arrangements, upwinding treatments, and conservativeness of convection terms. An artificial source term is introduced, and the solutions are compared with those obtained by the conventional marker-and-cell (MAC) method. The comparisons favor the PG method. Numerical results obtained by the twelve schemes are compared with exact solutions in order to assess the stability and accuracy of each scheme.

Tan, C.-H.↗

Variations in eddy diffusion and associated transport processes

Variation in the eddy diffusion coefficient constitutes a momentum source. This variation near the turbopause affects both the composition and the wind and temperature fields of the thermosphere. Two types of calculations are carried out, for horizontal winds considered self-consistently in a three-dimensional dynamic model and horizontal winds artificially suppressed to simulate a one-dimensional model. Momentum transfer associated with variations in eddy diffusion coefficients represents the source function in the analysis presented. It is shown that for long-term variations in the thermosphere, the compositional dependence on variations in eddy diffusion coefficients can be entirely different in one and three dimensional models with horizontal winds included. Horizontal winds can greatly reduce the effects of eddy diffusion in general and the ratio between He and O amplitudes in particular.

Mayr, H. G.↗

Radiation properties and emissivity parameterization of high level thin clouds

To parameterize emissivity of clouds at 11 microns, a study has been made in an effort to understand the radiation field of thin clouds. The contributions to the intensity and flux from different sources and through different physical processes are calculated by using the method of successive orders of scattering. The effective emissivity of thin clouds is decomposed into the effective absorption emissivity, effective scattering emissivity, and effective reflection emissivity. The effective absorption emissivity depends on the absorption and emission of the cloud; it is parameterized in terms of optical thickness. The effective scattering emissivity depends on the scattering properties of the cloud; it is parameterized in terms of optical thickness and single scattering albedo. The effective reflection emissivity follows the similarity relation as in the near infrared cases. This is parameterized in terms of the similarity parameter and optical thickness, as well as the temperature difference between the cloud and ground.

Wu, M.-L. C.↗

Co-Seismic Mass Dislocation and its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the shaking that is the earthquake, leaves behind permanent (step-function-like) dislocations in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field (in terms of spherical harmonic Stokes coefficients). The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results based on Chao & Gross (1987). The calculation uses the normal mode summation scheme, applied to nearly twenty thousand major earthquakes that occurred during 1976-2002, according to source mechanism solutions given by the Harvard Central Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies. For example, earthquakes conspire to decrease J2 and J22 while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to nudge the Earth rotation pole towards approximately 140 degrees E, roughly opposite to the observed polar drift direction. The geophysical significance and implications will be further studied.

Chao, B. F.↗

Tonlé Sap Food Security and Agriculture: Evaluating the Effects of Land Use and Hydrological Change on Ecosystem Vitality using Remotely-Sensed Data in the Tonlé Sap Lake Basin

Tonlé Sap Lake, the largest lake in Southeast Asia, is a critical source of fish and freshwater resources for the region. The health of this freshwater system is under pressure from accelerating dam construction, intensifying agriculture, deforestation, and changing climate patterns, forcing tradeoffs between immediate food security and the long-term vitality and productivity of the ecosystem. Efficient freshwater system monitoring is crucial to navigating these challenges. In collaboration with Conservation International, the Cambodian Ministry of Water Resources and Meteorology, and the Tonlé Sap Authority, we developed and tested remotely-sensed proxies for sub-indicators of the Freshwater Health Index (FHI), which is typically calculated using in situ datasets. We used landcover datasets derived from Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI), PROBA-V Vegetation sensor (VGT), Sentinel-2 Multispectral Imager (MSI), Advanced Very High Resolution Radiometer (AVHRR), and Envisat Medium Resolution Imaging Spectrometer (MERIS) as inputs to calculate land cover naturalness and bank modification. Additionally, we created a lake-level time series using a collection of altimetry data sources to estimate deviation from natural flow. We observed a decrease in landcover naturalness and a breakdown in the volume and regularity of annual lake levels from 2000-2020, reflecting increased pressure on water supply and agricultural productivity. At least 8% of forested areas in the basin were lost and rice harvest intensity increased over the course of the study period. These results will help our partners make informed decisions regarding freshwater management. Furthermore, our remotely-sensed FHI analysis can be replicated in other regions, providing decision makers with a snapshot of freshwater health in data-scarce environments.

Marco Vallejos↗

Investigations of medium wavelength magnetic anomalies in the eastern Pacific using MAGSAT data

The computer program package for dealing with MAGSAT data was refined and a program listing and deck of cards was forwarded for use by interested investigators. Instead of doing equivalent source calculations, it is proposed to use spherical harmonics to model the observed data since spherical harmonic functions are orthogonal and addition of higher degree and order harmonics can be done without recalculation of the harmonics of lower degree and order already calculated. These functions are much more easy to manipulate in that the height of the term of the observation can be taken care of with great facility. Downward calculation (after suitable filtering) can also be performed very simply.

Harrison, C. G. A.↗

Noniterative three-dimensional grid generation using parabolic partial differential equations

A new algorithm for generating three-dimensional grids has been developed and implemented which numerically solves a parabolic partial differential equation (PDE). The solution procedure marches outward in two coordinate directions, and requires inversion of a scalar tridiagonal system in the third. Source terms have been introduced to control the spacing and angle of grid lines near the grid boundaries, and to control the outer boundary point distribution. The method has been found to generate grids about 100 times faster than comparable grids generated via solution of elliptic PDEs, and produces smooth grids for finite-difference flow calculations.

Edwards, T. A.↗

Infrared spectra of rotating protostars

Earlier calculations of the infrared emission expected from stars in the process of being made are corrected to include the most important observable effects of rotation and generalized. An improved version of the spherical model of a previous paper is developed, and the corresponding emergent spectral energy distributions are calculated for the theoretically expected mass infall rate in the cores of cool and quiescent molecular clouds. The dust grain opacity model and the temperature profile parameterization are improved. It is shown that the infrared spectrum of the IRAS source 04264+2426, which is associated with a Herbig-Haro object, can be adequately represented in terms of a rotating and accreting protostar. This strengthens the suggestion that collimated outflows in young stellar objects originate when a stellar wind tries to emerge and reverse the swirling pattern of infall which gave birth to the central star.

Adams, F. C.↗

Atomic data from the iron project. 3: Rate coefficients for electron impact excitation of boron-like ions: Ne VI, Mg VIII, Al IX, Si X, S XII, Ar XIV, Ca XVI and Fe XXII

Collison strengths and maxwellian averaged rate coefficients have been calculated for the 105 transitions among all 15 fine structure levels of the 8 LS terms 2s(sup 2) 2 P(P-2(sup 0 sub 1/, 3/2)), 2s2p(sup 2)(P-4(sub 1/2,3/2,5/2), D-2(sub 3/2, 5/2), S-2(sub 1/2), P-2(sub 1/2, 3/2)), 2p(sup 3)(S-4(sup 0)(sub 3/2), D-2(sup 0 sub 3/2, 5/2), P-2(sup 0 sub 1/2, 3/2)) in highly- charged B-like Ne, Mg, Al, Si, S, Ar, Ca and Fe. Rate coefficients have been tabulated at a wide range of temperatures, depending on the ion charge and abundance in plasma sources. Earlier work for O IV has also been extended to include the high temperature range. A brief discussion of the calculations, sample results, and comparison with earlier works is also given. While much of the new data should be applicable to UV spectral diagnostics, the new rates for the important ground state fine structure transition P-2(sup 0 sub 1/2)-P-2(sup 0 sub 3/2) should result in significant revision of the IR cooling rates in plasmas where B-like ions are prominent constituents, since the new rate coefficients are generally higher by several factors compared with the older data.

Zhang, Hong Lin↗

Variational approach to probabilistic finite elements

Probabilistic finite element method (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties, and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

Variational approach to probabilistic finite elements

Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

NASA Glenn Coefficients for Calculating Thermodynamic Properties of Individual Species

This report documents the library of thermodynamic data used with the NASA Glenn computer program CEA (Chemical Equilibrium with Applications). This library, containing data for over 2000 solid, liquid, and gaseous chemical species for temperatures ranging from 200 to 20,000 K, is available for use with other computer codes as well. The data are expressed as least-squares coefficients to a seven-term functional form for C((sup o)(sub p)) (T) / R with integration constants for H (sup o) (T) / RT and S(sup o) (T) / R. The NASA Glenn computer program PAC (Properties and Coefficients) was used to calculate thermodynamic functions and to generate the least-squares coefficients. PAC input was taken from a variety of sources. A complete listing of the database is given along with a summary of thermodynamic properties at 0 and 298.15 K.

Bonnie J McBride↗

Model 'zero-age' lunar thermal profiles resulting from electrical induction

Thermal profiles for the moon are calculated under the assumption that a pre-main-sequence T-Tauri-like solar wind excites both transverse magnetic and transverse electric induction while the moon is accreting. A substantial initial temperature rise occurs, possibly of sufficient magnitude to cause subsequent early extensive melting throughout the moon in conjunction with nominal long-lived radioactives. In these models, accretion is an unimportant direct source of thermal energy but is important because even small temperature rises from accretion cause significant changes in bulk electrical conductivity. Induction depends upon the radius of the moon, which we take to be accumulating while it is being heated electrically. The 'zero-age' profiles calculated in this paper are proposed as initial conditions for long-term thermal evolution of the moon.

Herbert, F.↗

The extension of the concept of the cosmic-ray path-length distribution to nonrelativistic energies

The concept of the cosmic-ray path-length distribution is examined. The corresponding cosmic-ray propagation calculational procedure has been justified theoretically at relativistic energies (Ginzburg and Syrovatskii, 1964) where the effects of ionization energy loss are negligible. The present paper extends the use of the path-length distribution concept in cosmic-ray propagation calculations to nonrelativistic energies. Sufficient constraints to effect this extension are presented. The solution of the cosmic-ray propagation equations in terms of a Green's function approach is also investigated and is used to provide a formulation of the path-length distribution at nonrelativistic as well as relativistic energies in terms of the cosmic-ray source distribution and the propagation characteristics of the interstellar medium. The leaky-box model of cosmic-ray propagation is also examined.

Lezniak, J. A.↗

Biological effectiveness of high-energy protons - Target fragmentation

High-energy protons traversing tissue produce local sources of high-linear-energy-transfer ions through nuclear fragmentation. The contribution of these target fragments to the biological effectiveness of high-energy protons using the cellular track model is examined. The effects of secondary ions are treated in terms of the production collision density using energy-dependent parameters from a high-energy fragmentation model. Calculations for mammalian cell cultures show that at high dose, at which intertrack effects become important, protons deliver damage similar to that produced by gamma rays, and with fragmentation the relative biological effectiveness (RBE) of protons increases moderately from unity. At low dose, where sublethal damage is unimportant, the contribution from target fragments dominates, causing the proton effectiveness to be very different from that of gamma rays with a strongly fluence-dependent RBE. At high energies, the nuclear fragmentation cross sections become independent of energy. This leads to a plateau in the proton single-particle-action cross section, below 1 keV/micron, since the target fragments dominate.

Cucinotta, Francis A.↗

Calculation of the three-dimensional flow field in supersonic inlets at angle of attack using a bicharacteristic method with discrete shock wave fitting

An analysis is presented for calculating the flow field in supersonic mixed-compression aircraft inlets operating at angle of attack. The flow field is computed by a steady three-dimensional bicharacteristic method. The bow shock wave and the reflected internal shock wave system are computed by a three-dimensional discrete shock wave fitting procedure. Viscous and thermal diffusion may be included as source terms in the bicharacteristic method. A production type computer program capable of determining the flow field in a variety of axisymmetric mixed-compression supersonic inlets is available. The results of the present analysis agree well with those produced by the two-dimensional method of characteristics when axisymmetric flow fields are computed. For three-dimensional flow fields, the results of the present analysis agree well with experimental data except in regions of high viscous interaction and boundary layer removal. The present analysis does not compute the boundary layer, nor does it account for boundary layer bleed.

Vadyak, J.↗

On the transport and propagation of relativistic electrons in galaxies

Analytic solutions to the time-dependent equation describing the transport of relativistic electrons that are being convected and diffused, and may be either gaining or losing energy, are described. The calculations yield a body of solutions from which specific ones can be chosen for the modelling of synchrotron radiation in the Galaxy and in other galaxies. The solutions allow for (1) spatially varying convection velocity; (2) an energy and spatially dependent diffusion coefficient; (3) adiabatic acceleration effects; (4) a term describing spatially dependent synchrotron loss and/or Fermi acceleration; (5) an injection source for the electrons that may be impulsive, steady-state, or variable-injection; (6) a spatial injection or spatially more diffuse source; (7) arbitrary energy dependence of the injection source; and (8) spatial and energy boundaries where the spatial or energy form of the convection, diffusion, and energy supply terms in the basic equation can change character.

Lerche, I.↗

Magnetic confinement of cosmic clouds

The role of the magnetic field in the confinement or compression of interstellar gas clouds is reconsidered. The virial theorem for an isolated magnetized cloud in the presence of distant magnetic sources is reformulated in terms of moments of the internal and external currents, and an equilibrium condition is derived. This condition is applied to the interaction between isolated clouds for the simple- and artificial-case in which the field of each cloud is a dipole. With the simplest of statistical assumptions, the probability of any given cloud being compressed is calculated as about 10 percent, the magnetic field acting as a medium which transmits the kinetic pressure between clouds. Even when compression occurs the magnetic pressure 1/2 B-squared may decrease on leaving the cloud surface.

Azar, Michel↗