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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Delayed Gamma Signature Accounting for Field Line Transport of Debris Ions after a Nuclear Explosion

Delayed gammas (DGs) are produced by the decay of radioactive fission fragments after a nuclear explosion. Currently, the time and energy profile of the delayed gamma signal is propagated in DIORAMA to the GEO satellite sensors by assuming a point source at the burst location for the production of the delayed gammas. However, in reality the production of delayed gammas is spread over some spatial extent. First, the debris ions are confined within an expanding region by the magnetic field and/or surrounding plasma. This expanding confinement bubble eventually collapses, at which point the debris ions can be trapped on magnetic field lines and transported away from the original burst location. In this work the effects of debris ion transport on the associated time profile of delayed gammas at a GEO satellite is presented. The first step was the development of the Delayed Gamma Source Model (DGSM), which is based on the analytical Electron Source Model (ESM) code developed at LANL that models the spatial distribution of energetic electrons produced after a high-altitude nuclear explosion (HANE) for use as a source term in LANL artificial belt diffusion code. This model creates a spatial grid in latitude, longitude, and altitude of the distribution of radioactive debris ions. The output of this code is combined with a separate analysis code also described, which produces the DG time profile signature at a GEO sensor faceplate when accounting for the spatial spread of delayed gamma production during field line transport. Some simple case studies are presented to show how the DG signature at GEO calculated with this higher-fidelity model and a simple point-source model, as is currently assumed in DIORAMA, compare.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

Modeling the dust size distribution in comets with dust fragmentation

A hydrodynamic model was developed of a spherically symmetric dusty gas flow in a cometary atmosphere assuming a single fluid, inviscid, perfect gas. The hydrodynamics for gas and dust, which involves the gas drag force (momentum transfer), heat exchange between gas and dust, photodissociation energy for H2O gas, and radiative heating and cooling terms for dust particles, are solved using the Gear method for stiff, coupled differential equations. Calculations were done with a dust size distribution for radii alpha = 0.01 micron to 10 cm with densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g/cu cm and a total dust-to-gas mass ratio chi = 1 was adopted. There are indications from in situ observations that dust particle fragment into smaller ones. Fragmentation of dust particles was incorporated into the model. This is done by adding source and sink terms in the continuity equations for the dust. Lifetimes for the decay of dust particles were assumed as a function of particle size. It is also assumed that dust particles always fragment only into the next smaller size.

Konno, Ichishiro↗
NASA NTRS1975

Short-term time variability of Cygnus X-1

Autocorrelation functions, power-density spectra, and cross-correlation functions were calculated for 71 observations of Cyg X-1 by the Uhuru satellite. The average autocorrelation function and power-density spectrum were found to be independent of energy and consistent with a shot-noise source with a characteristic time constant of 0.45 s. The data in different energy bands were found to be 100 percent correlated.

Weisskopf, M. C.↗
NASA NTRS1974

The generation of gravitational waves. 1. Weak-field sources: A plug-in-and-grind formalism

A plug-in-and-grind formalism is derived for calculating the gravitational waves emitted by any system with weak internal gravitational fields. If the internal fields have negligible influence on the system's motions, then the formalism reduces to standard linearized theory. Whether or not gravity affects the motions, if the motions are slow and internal stresses are weak, then the new formalism reduces to the standard quadrupole-moment formalism. In the general case the new formalism expresses the radiation in terms of a retarded Green's function for slightly curved spacetime, and then breaks the Green's-function integral into five easily understood pieces: direct radiation, produced directly by the motions of the sources; whump radiation, produced by the the gravitational stresses of the source; transition radiation, produced by a time-changing time delay (Shapiro effect) in the propagation of the nonradiative, 1/r field of the source; focussing radiation produced when one portion of the source focusses, in a time-dependent way, the nonradiative field of another portion of the source, and tail radiation, produced by backscatter of the nonradiative field in regions of focussing.

Thorne, K. S.↗
DOE OSTI2023

A general predictive methodology for fuel-mixture properties up to supercritical conditions

A predictive thermodynamic model is utilized for the calculation of fuel properties of oxymethylene dimethyl ethers (OME 3–4 ), surrogates for gasoline, diesel and aviation fuel, as well as alcohol blends with gasoline and diesel. The alcohols used for these blends are methanol, ethanol, propanol, butanol and pentanol; their mixing ratio ranges from 10 to 50% by volume. The model is based on the Perturbed-Chain Statistical Association Fluid Theory (PC-SAFT) equation of state (EoS) and Vapor Liquid Equilibrium (VLE) calculations at constant temperature, density and composition. The model includes the association term, with the assumption of two association sites (2B scheme), to enable the modeling of alcohols. The pure-component parameters are estimated based on the Group Contribution (GC) method of various sources, as well as a parametrization model specifically designed for the case of OME 3–4 . The results of the computational model for the density, vapor pressure and distillation curves at various conditions, including high-pressure, high-temperature (HPHT), are compared to experimental and computational data available in the literature. In the cases where no measurements are available for the surrogates, experimental data for the corresponding target fuel are used, taking into consideration the inherent deviation in properties between real and surrogate fuel. Overall, the results are in good agreement with the data from the literature, with the average deviation not exceeding 12% for temperature (Kelvin) on the distillation curves, 10% for density and 46% for vapor pressure and the general trend being captured successfully. The use of different pure component parameter estimation techniques can further improve the prediction quality in the cases of OME3–4 and the aviation fuel surrogate, especially for the vapor pressure, leading to an average deviation lower than 18%. These results demonstrate the predictive capabilities of the model, which extend to a wide range of fuel types and pressure/temperature conditions. Through this investigation, the present work aims to establish the limits of applicability of this thermodynamic property prediction methodology.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

(U) A Theoretical Study of Asay Foil Trajectories (V.02)

We consider the trajectory of a generic Asay foil ejecta momentum diagnostic sensor, for scenarios where ejecta are produced at a planar surface and fly ballistically through a perfect vacuum to the sensor. To do so, we build upon a previously established mathematical formalism derived for the analytic study of stationary sensors (i.e., piezopins). First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source areal mass function (“source model”). This defines an integro-differential equation (IDE) for the foil trajectory. When ejecta production is instantaneous - as is generally assumed in momentum diagnostic data analyses - the IDE leads to an implicit and easily calculable closed-form solution for the foil trajectory in a perfect system, as long as the ejecta particle velocity distribution is twice-integrable. General properties of the instant-production trajectory solution indicate the existence of a boundary condition the particle velocity distribution must satisfy in order for the analytically predicted foil trajectories to be compatible with certain features commonly observed in foil data. This boundary condition is identical to one derived previously from a consideration of piezopin data. Armed with the analytic solution for instant production, we also consider various techniques used to extract time-dependent cumulative ejecta masses from foil trajectories, and derive an expression for the error imposed by using an approximated equation of motion. This analytic trajectory solution furthermore makes it possible to examine the common practice of presenting inferred cumulative ejecta masses as a function of a normalized implied velocity; we derive conditions under which this methodology is and is not meaningful. We also propose a strategy for extending the instant-production trajectory solution to time-dependent source functions.

42 ENGINEERING↗
NASA NTRS1973

Ionospheric scintillation studies

The diffracted field of a monochromatic plane wave was characterized by two complex correlation functions. For a Gaussian complex field, these quantities suffice to completely define the statistics of the field. Thus, one can in principle calculate the statistics of any measurable quantity in terms of the model parameters. The best data fits were achieved for intensity statistics derived under the Gaussian statistics hypothesis. The signal structure that achieved the best fit was nearly invariant with scintillation level and irregularity source (ionosphere or solar wind). It was characterized by the fact that more than 80% of the scattered signal power is in phase quadrature with the undeviated or coherent signal component. Thus, the Gaussian-statistics hypothesis is both convenient and accurate for channel modeling work.

Rino, C. L.↗
NASA NTRS1980

Gravitational waves from rotating and precessing rigid bodies. II - General solutions and computationally useful formulas

A rigid, freely precessing Newtonian body emits gravitational radiation. In this paper the classical-mechanics results for free precession which are needed in order to calculate the weak-field slow-motion quadrupole-moment gravitational waves are reviewed. Within that formalism, algorithms for computing the exact gravitational power radiated and waveforms produced by arbitrary rigid-body freely precessing source are given. Also presented are the dominant terms in series expansions of the waveforms for the case of an almost-spherical object precessing with a small wobble angle. These series expansions, which retain the precise frequency dependence of the waves, may be useful for gravitational astronomers when freely precessing sources begin to be observed.

Zimmermann, M.↗
NASA NTRS1987

Sensitivity studies and a simple ozone perturbation experiment with a truncated two-dimensional model of the stratosphere

The 1-1/2-D model of Holton (1986), which is actually a highly truncated two-dimensional model, describes latitudinal variations of tracer mixing ratios in terms of their projections onto second-order Legendre polynomials. The present study extends the work of Holton by including tracers with photochemical production in the stratosphere (O3 and NOy). It also includes latitudinal variations in the photochemical sources and sinks, improving slightly the calculated global mean profiles for the long-lived tracers studied by Holton and improving substantially the latitudinal behavior of ozone. Sensitivity tests of the dynamical parameters in the model are performed, showing that the response of the model to changes in vertical residual meridional winds and horizontal diffusion coefficients is similar to that of a full two-dimensional model. A simple ozone perturbation experiment shows the model's ability to reproduce large-scale latitudinal variations in total ozone column depletions as well as ozone changes in the chemically controlled upper stratosphere.

Stordal, Frode↗
NASA NTRS1990

Near-field noise of a single-rotation propfan at an angle of attack

The near field noise characteristics of a propfan operating at an angle of attack are examined utilizing the unsteady pressure field obtained from a 3-D Euler simulation of the propfan flowfield. The near field noise is calculated employing three different procedures: a direct computation method in which the noise field is extracted directly from the Euler solution, and two acoustic-analogy-based frequency domain methods which utilize the computed unsteady pressure distribution on the propfan blades as the source term. The inflow angles considered are -0.4, 1.6, and 4.6 degrees. The results of the direct computation method and one of the frequency domain methods show qualitative agreement with measurements. They show that an increase in the inflow angle is accompanied by an increase in the sound pressure level at the outboard wing boom locations and a decrease in the sound pressure level at the (inboard) fuselage locations. The trends in the computed azimuthal directivities of the noise field also conform to the measured and expected results.

Nallasamy, M.↗
NASA NTRS1990

Near-field noise of a single rotation propfan at an angle of attack

The near-field noise characteristics of a propfan operating at an angle of attack are examined utilizing the unsteady pressure field obtained from a three-dimensional Euler simulation of the propfan flowfield. The near-field noise is calculated employing three different procedures: a direct computation method in which the noise field is extracted directly from the Euler solution, and two acoustic-analogy-based frequency domain methods which utilize the computed unsteady pressure distribution on the propfan blades as the source term. The inflow angles considered are -0.4, 1.6, and 4.6 degrees. The results of the direct computation method and one of the frequency domain methods show qualitative agreement with measurements. They show that an increase in the inflow angle is accompanied by an increase in the sound pressure level at the outboard wing boom locations and a decrease in the sound pressure level at the (inboard) fuselage locations. The trends in the computed azimuthal directivities of the noise field also conform to the measured and expected results.

Nallasamy, M.↗
NASA NTRS1977

Solar seismology. I - The stability of the solar p-modes

The stability of the radial p-modes of the sun is investigated by computing nonadiabatic eigenvalues and eigenfunctions for a solar envelope model which extends from an inner radius of about 0.3 solar radius out to an optical depth of about 0.0003. The calculations take into account in a crude fashion the response of the convective flux to the oscillation. The dynamical effect of turbulence in the convection zone is parametrized in terms of a turbulent shear viscosity. The results show that if damping by turbulent viscosity is neglected, all modes with periods longer than 6 minutes are unstable. The familiar kappa-mechanism, which operates in the H ionization-H(-) opacity region, is the dominant source of driving of the oscillations. Modes with periods shorter than 6 minutes are stabilized by radiative damping in the solar atmosphere. When turbulent dissipation of pulsational energy is included, all modes are predicted to be stable. However, the margin of stability is very small. In view of the large uncertainty that must be assigned to the estimate of turbulent damping, it is concluded that theoretical calculations cannot unequivocally resolve the question of the stability of the solar p-modes.

Goldreich, P.↗
DOE OSTI2023

A geospatial risk analysis graphical user interface for identifying hazardous chemical emission sources

Background: Performing back trajectory and forward trajectory using the Hybrid Single-Particle Lagrangian Integrated Trajectory Model (HYSPLIT) is a reliable approach for assessing particle transport after release among mid-field atmospheric models. HYSPLIT has an externally facing online interface that allows non-expert users to run the model trajectories without requiring extensive training or programming. However, the existing HYSPLIT interface is limited if simulations have a large amount of meteorological data and timesteps that are not coincident. The objective of this study is to design and develop a more robust tool to rapidly evaluate hazard transport conditions and to perform risk analysis, while still maintaining an intuitive and user-friendly interface. Methods: HYSPLIT calculates forward and backward trajectories of particles based on wind speed, wind direction, and the corresponding location, timestamp, and Pasquill stability classes of the regions of the atmosphere in terms of the wind speed, the amount of solar radiation, and the fractional cloud cover. The computed particle transport trajectories, combined with the online Proton Transfer Reaction-Mass Spectrometry (PTR-MS) data (https://figshare.com/articles/dataset/ARL_Data_from_PROS_station_at_Hanford_site/19993964), can be used to identify and quantify the sources and affected area of the hazardous chemicals’ emission using the potential source distribution function (PSDF). PSDF is an improved statistical function based on the well-known potential source contribution function (PSCF) in establishing the air pollutant source and receptor relationship. Performing this analysis requires a range of meteorological and pollutant concentration measurements to be statistically meaningful. The existing HYSPLIT graphical user interface (GUI) does not easily permit computations of trajectories of a dataset of meteorological data in high temporal frequency. To improve the performance of HYSPLIT computations from a large dataset and enhance risk analysis of the accidental release of material at risk, a geospatial risk analysis tool (GRAT-GUI) is created to allow large data sets to be processed instantaneously and to provide ease of visualization. Results: The GRAT-GUI is a native desktop-based application and can be run in any Windows 10 system without any internet access requirements, thus providing a secure way to process large meteorological datasets even on a standalone computer. GRAT-GUI has features to import, integrate, and convert meteorological data with various formats for hazardous chemical emission source identification and risk analysis as a self-explanatory user interface. The tool is available at https://figshare.com/articles/software/GRAT/19426742.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Numerical study of hydrogen-air supersonic combustion by using elliptic and parabolized equations

The two-dimensional Navier-Stokes and species continuity equations are used to investigate supersonic chemically reacting flow problems which are related to scramjet-engine configurations. A global two-step finite-rate chemistry model is employed to represent the hydrogen-air combustion in the flow. An algebraic turbulent model is adopted for turbulent flow calculations. The explicit unsplit MacCormack finite-difference algorithm is used to develop a computer program suitable for a vector processing computer. The computer program developed is then used to integrate the system of the governing equations in time until convergence is attained. The chemistry source terms in the species continuity equations are evaluated implicitly to alleviate stiffness associated with fast chemical reactions. The problems solved by the elliptic code are re-investigated by using a set of two-dimensional parabolized Navier-Stokes and species equations. A linearized fully-coupled fully-implicit finite difference algorithm is used to develop a second computer code which solves the governing equations by marching in spce rather than time, resulting in a considerable saving in computer resources. Results obtained by using the parabolized formulation are compared with the results obtained by using the fully-elliptic equations. The comparisons indicate fairly good agreement of the results of the two formulations.

Chitsomboon, T.↗
DOE OSTI2023

Convective and Turbulent Motions in Nonprecipitating Cu. Part III: Characteristics of Turbulence Motions

Velocity field in a nonprecipitating Cu under BOMEX conditions, simulated by SAM with 10-m resolution and spectral bin microphysics is separated into the convective part and the turbulent part, using a wavelet filtering. In Part II of the study properties of convective motions of this Cu were investigated. Here in Part III of the study, the parameters of cloud turbulence are calculated in the cloud updraft zone at different stages of cloud development. The main points of this study are (i) application of a fine-scale LES model of a single convective cloud allowed a direct estimation of turbulence parameters using the resolved flow in the cloud and (ii) the separation of the resolved flow into the turbulence flow and the nonturbulence flow allowed us to estimate different turbulent parameters with sufficient statistical accuracy. We calculated height and time dependences of the main turbulent parameters such as turbulence kinetic energy (TKE), spectra of TKE, dissipation rate, and the turbulent coefficient. It was found that the main source of turbulence in the cloud is buoyancy whose contribution is described by the buoyancy production term (BPT). The shear production term (SPT) increases with height and reaches its maximum near cloud top, and so does BPT. In agreement with the behavior of BPT and SPT, turbulence in the lower cloud part (below the inversion level) is weak and hardly affects the processes of mixing and entrainment. The fact that BPT is larger than SPT determines many properties of cloud turbulence. For instance, the turbulence is nonisotropic, so the vertical component of TKE is substantially larger than the horizontal components. Another consequence of the fact that BPT is larger than STP manifests itself in the finding that the turbulence spectrum largely obeys the -11/5 Bolgiano–Obukhov scaling. The classical Kolmogorov -5/3 scaling dominates for the low part of a cloud largely at the dissolving stage of cloud evolution. Using the spectra obtained we evaluated an “effective” dissipation rate which increases with height from nearly zero at cloud base up to 20 cm 2 s -3 near cloud top. The coefficient of turbulent diffusion was found to increase with height and ranged from 5 m 2 s -1 near cloud base to 25 m 2 s -1 near cloud top. In conclusion, the possible role of turbulence in the process of lateral entrainment and mixing is discussed.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Trailing Edge Noise Prediction Based on a New Acoustic Formulation

A new analytic result in acoustics called 'Formulation 1B,' proposed by Farassat, is used to compute broadband trailing edge noise from an unsteady surface pressure distribution on a thin airfoil in the time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with the loading source term, and has been shown in previous research to provide time domain predictions of broadband noise that are in excellent agreement with experiment. Furthermore, this formulation lends itself readily to rotating reference frames and statistical analysis of broadband trailing edge noise. Formulation 1B is used to calculate the far field noise radiated from the trailing edge of a NACA 0012 airfoil in low Mach number flows, using both analytical and experimental data on the airfoil surface. The results are compared to analytical results and experimental measurements that are available in the literature. Good agreement between predictions and measurements is obtained.

Casper, J.↗
NASA NTRS2017

Direct Estimation of the Global Distribution of Vertical Velocity Within Cirrus Clouds

Cirrus clouds determine the radiative balance of the upper troposphere and the transport of water vapor across the tropopause.The representation of vertical wind velocity, W, in atmospheric models constitutes the largest source of uncertainty in the calculation of the cirrus formation rate. Using global atmospheric simulations with a spatial resolution of 7 kilometers we obtain for the first time a direct estimate of the distribution of W at the scale relevant for cirrus formation, validated against long-term observations at two different ground sites. The standard deviation in W, sigma (sub W), varies widely over the globe with the highest values resulting from orographic uplift and convection, and the lowest occurring in the Arctic. Globally about 90 of the simulated sigma (sub W) values are below 0.1 meters per second and about one in 10 (sup 4) cloud formation events occur in environments with sigma (sub W) greater than 0.8 meters per second. Combining our estimate with reanalysis products and an advanced cloud formation scheme results in lower homogeneous ice nucleation frequency than previously reported, and a decreasing average ice crystal concentration with decreasing temperature. These features are in agreement with observations and suggest that the correct parameterization of sigma (sub W) is critical to simulate realistic cirrus properties.

Radiative Balance↗
DOE OSTI2021

Uranium Attenuated by a Wetland 50 Years after Release into a Stream

: Wetlands have several important roles in the hydrological cycle, including maintaining water quality by removing surface and groundwater contaminants. Over time, the wetlands themselves can become contaminated, posing a secondary environmental threat. The objective of this study was to calculate the inventory of uranium (U) remaining in a wetland 50 years after the Fuel Fabrication Facility on the Savannah River Site located in South Carolina released 43.5 Mg of U into the nearby environment. Over 232 700 gamma spectra and their associated global positioning system (GPS) coordinates were collected and collated into a map of the contaminated land area. Five core samples were also collected that contained U concentrations as high as 14 099 mg/kg (background levels are about 2.7 mg/kg U). The contaminated area was 278 000 m2 , and it contained 36.2 Mg U, about 83% of the U released. About 80% of the U in the wetland was concentrated in a former beaver pond, a 73 000 m2 area (26% of the contaminated area). This contaminated wetland area was almost 2 km from the source, indicating that it comprised unique hydro-biogeochemical properties for immobilizing the released U. To the best of our knowledge, this is the first data-rich study to quantify the long-term effectiveness of a wetland to immobilize inorganic contaminants. Significant environmental changes to the system, such as those associated with hydrology, forest fires, or anthropogenic land use, may alter the complex hydro-biogeochemical interactions necessary for the long-term immobilization of the U.

Contaminant inventory↗