Search NASA⌕ Search

SEARCH · Search NASA

Results for “source terms”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

A Fast Method of Deriving the Kirchhoff Formula for Moving Surfaces

The Kirchhoff formula for a moving surface is very useful in many wave propagation problems, particularly in the prediction of noise from rotating machinery. Several publications in the last two decades have presented derivations of the Kirchhoff formula for moving surfaces in both time and frequency domains. Here we present a method originally developed by Farassat and Myers in time domain that is both simple and direct. It is based on generalized function theory and the useful concept of imbedding the problem in the unbounded three-dimensional space. We derive an inhomogeneous wave equation with the source terms that involve Dirac delta functions with their supports on the moving data surface. This wave equation is then solved using the simple free space Green's function of the wave equation resulting in the Kirchhoff formula. The algebraic manipulations are minimal and simple. We do not need the Green's theorem in four dimensions and there is no ambiguity in the interpretation of any terms in the final formulas. Furthermore, this method also gives the simplest derivation of the classical Kirchhoff formula which has a fairly lengthy derivation in physics and applied mathematics books. The Farassat-Myers method can be used easily in frequency domain.

Farassat, F.↗
NASA NTRS2002

A New Time Domain Formulation for Broadband Noise Predictions

A new analytic result in acoustics called "Formulation 1B," proposed by Farassat, is used to compute the loading 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. The formulation contains a far field surface integral that depends on the time derivative and the surface gradient of the pressure on the airfoil, as well as a contour integral on the boundary of the airfoil surface. As a first test case, the new formulation is used to compute the noise radiated from a flat plate, moving through a sinusoidal gust of constant frequency. The unsteady surface pressure for this test case is analytically specied from a result based on linear airfoil theory. This test case is used to examine the velocity scaling properties of Formulation 1B and to demonstrate its equivalence to Formulation 1A of Farassat. The new acoustic formulation, again with an analytic surface pressure, is then used to predict broadband noise radiated from an airfoil immersed in homogeneous, isotropic turbulence. The results are compared with experimental data previously reported by Paterson and Amiet. Good agreement between predictions and measurements is obtained. Finally, an alternative form of Formulation 1B is described for statistical analysis of broadband noise.

Casper, Jay H.↗
NASA NTRS2009

Performance of MODIS Thermal Emissive Bands On-orbit Calibration Algorithms

Two nearly identical copies of the Moderate Resolution Imaging Spectroradiometer (MODIS) are currently operated on-board the Terra and Aqua spacecrafts, launched in December 1999 and May 2002, respectively. Together, they have produced an unprecedented amount of science data products, which are widely used for the studies of changes in the Earth's system of land, oceans, and atmosphere. MODIS is a cross-track scanning radiometer, which uses a two-sided scan mirror and collects data continuously over a wide scan angle range (+/-55 degree relative to the instrument nadir) each scan of 1.47 seconds. It has 36 spectral bands with wavelengths ranging from visible (VIS) to long-wave infrared (LWIR). MODIS bands 1-19 and 26 are the reflective solar bands (RSB) and bands 20-25 and 27-36 are the thermal emissive bands (TEB). MODIS was developed and designed with improvements made over its heritage sensors (such as AVHRR and Landsat) and, in particular, with more stringent calibration requirements. Because of this, MODIS was built with a set of state-of-art on-board calibrators (OBC), which include a solar diffuser (SD), a solar diffuser stability monitor (SDSM), a blackbody (BB), a spectroradiometric calibration assembly (SRCA), and a space view (SV) port. With the exception of view angle differences, MODIS OBC measurements and the Earth View (EV) observations are made via the same optical path. MODIS TEB have a total of 160 individual TEB detectors (10 per band), which are located on two cold focal plane assemblies (CFPA). For nominal on-orbit operation, the CFPA temperature is controlled at 83K via a passive radiative cooler. For the TEB, the calibration requirements at specified typical scene radiances are less than or equal to 1% with an exception for the fire detection (low gain) band. MODIS TEB on-orbit calibration is performed on a scan-by-scan basis using a quadratic calibration algorithm, and data collected from sensor responses to the onboard BB and SV. The BB serves as the thermal calibration source and the SV provides measurements for the sensor's background and offsets. MODIS on-board BB is a v-grooved plate with its temperature measured using 12 platinum resistive thermistors (PRT) uniformly embedded in the BB substrate. All the BB thermistors were characterized pre-launch with reference to the NIST temperature standards. Unlike typical BB operations in many heritage sensors, which have no temperature control capability, the MODIS on-board BB can be operated at any temperatures between instrument ambient (about 270K) and 315K and can also be varied continuously within this range. This feature has significantly enhanced the MODIS' capability of tracking and updating the TEB nonlinear calibration coefficients over its entire mission. Following a brief description of MODIS TEB on-orbit calibration methodologies and its onboard BB operational activities, this paper provides a comprehensive performance assessment of MODIS TEB quadratic calibration algorithm. It examines the scan-by-scan, orbit-by-orbit, daily, and seasonal variations of detector responses and associated impact due changes in the CFPA and instrument temperatures. Specifically, this paper will analyze the contribution by each individual thermal emissive source term (BB, scan cavity, and scan mirror), the impact on the Level 1 B data product quality due to pre-launch and on-orbit calibration uncertainties. A comparison of Terra and Aqua TEB on-orbit performance, lessons learned, and suggestions for future improvements will also be made.

Xiong, Xiaoxiong↗
NASA NTRS2010

Meeting Radiation Protection Requirements and Reducing Spacecraft Mass - A Multifunctional Materials Approach

Both crew and radio-sensitive systems, especially electronics must be protected from the effects of the space radiation environment. One method of mitigating this radiation exposure is to use passive-shielding materials. In previous vehicle designs such as the International Space Station (ISS), materials such as aluminum and polyethylene have been used as parasitic shielding to protect crew and electronics from exposure, but these designs add mass and decrease the amount of usable volume inside the vehicle. Thus, it is of interest to understand whether structural materials can also be designed to provide the radiation shielding capability needed for crew and electronics, while still providing weight savings and increased useable volume when compared against previous vehicle shielding designs. In this paper, we present calculations and analysis using the HZETRN (deterministic) and FLUKA (Monte Carlo) codes to investigate the radiation mitigation properties of these structural shielding materials, which includes graded-Z and composite materials. This work is also a follow-on to an earlier paper, that compared computational results for three radiation transport codes, HZETRN, HETC, and FLUKA, using the Feb. 1956 solar particle event (SPE) spectrum. In the following analysis, we consider the October 1989 Ground Level Enhanced (GLE) SPE as the input source term based on the Band function fitting method. Using HZETRN and FLUKA, parametric absorbed doses at the center of a hemispherical structure on the lunar surface are calculated for various thicknesses of graded-Z layups and an all-aluminum structure. HZETRN and FLUKA calculations are compared and are in reasonable (18% to 27%) agreement. Both codes are in agreement with respect to the predicted shielding material performance trends. The results from both HZETRN and FLUKA are analyzed and the radiation protection properties and potential weight savings of various materials and materials lay-ups are compared.

Atwell, William↗
NASA NTRS2009

Coupled Particle Transport and Pattern Formation in a Nonlinear Leaky-Box Model

Effects of particle-particle coupling on particle characteristics in nonlinear leaky-box type descriptions of the acceleration and transport of energetic particles in space plasmas are examined in the framework of a simple two-particle model based on the Fokker-Planck equation in momentum space. In this model, the two particles are assumed coupled via a common nonlinear source term. In analogy with a prototypical mathematical system of diffusion-driven instability, this work demonstrates that steady-state patterns with strong dependence on the magnetic turbulence but a rather weak one on the coupled particles attributes can emerge in solutions of a nonlinearly coupled leaky-box model. The insight gained from this simple model may be of wider use and significance to nonlinearly coupled leaky-box type descriptions in general.

Barghouty, A. F.↗
NASA NTRS2010

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↗
NASA NTRS2006

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↗
NASA NTRS2004

Further Improvement in 3DGRAPE

3DGRAPE/AL:V2 denotes version 2 of the Three-Dimensional Grids About Anything by Poisson's Equation with Upgrades from Ames and Langley computer program. The preceding version, 3DGRAPE/AL, was described in Improved 3DGRAPE (ARC-14069) NASA Tech Briefs, Vol. 21, No. 5 (May 1997), page 66. These programs are so named because they generate volume grids by iteratively solving Poisson's Equation in three dimensions. The grids generated by the various versions of 3DGRAPE have been used in computational fluid dynamics (CFD). The main novel feature of 3DGRAPE/AL:V2 is the incorporation of an optional scheme in which anisotropic Lagrange-based trans-finite interpolation (ALBTFI) is coupled with exponential decay functions to compute and blend interior source terms. In the input to 3DGRAPE/AL:V2 the user can specify whether or not to invoke ALBTFI in combination with exponential-decay controls, angles, and cell size for controlling the character of grid lines. Of the known programs that solve elliptic partial differential equations for generating grids, 3DGRAPE/AL:V2 is the only code that offers a combination of speed and versatility with most options for controlling the densities and other characteristics of grids for CFD.

Alter, Stephen↗
NASA NTRS2011

Models of the Solar Atmospheric Response to Flare Heating

I will present models of the solar atmospheric response to flare heating. The models solve the equations of non-LTE radiation hydrodynamics with an electron beam added as a flare energy source term. Radiative transfer is solved in detail for many important optically thick hydrogen and helium transitions and numerous optically thin EUV lines making the models ideally suited to study the emission that is produced during flares. I will pay special attention to understanding key EUV lines as well the mechanism for white light production. I will also present preliminary results of how the model solar atmosphere responds to Fletcher & Hudson type flare heating. I will compare this with the results from flare simulations using the standard thick target model.

Allred, Joel↗
NASA NTRS2012

The High-Resolution Wave-Propagation Method Applied to Meso- and Micro-Scale Flows

The high-resolution wave-propagation method for computing the nonhydrostatic atmospheric flows on meso- and micro-scales is described. The design and implementation of the Riemann solver used for computing the Godunov fluxes is discussed in detail. The method uses a flux-based wave decomposition in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The two advantages of the technique are: 1) the need for an explicit definition of the Roe matrix is eliminated and, 2) the inclusion of source term due to gravity does not result in discretization errors. The resulting flow solver is conservative and able to resolve regions of large gradients without introducing dispersion errors. The methodology is validated against exact analytical solutions and benchmark cases for non-hydrostatic atmospheric flows.

Ahmad, Nashat N.↗
NASA NTRS2012

Magnetic Field Effects on Plasma Plumes

Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results

Ebersohn, F.↗
NASA NTRS2012

Radiation Coupling with the FUN3D Unstructured-Grid CFD Code

The HARA radiation code is fully-coupled to the FUN3D unstructured-grid CFD code for the purpose of simulating high-energy hypersonic flows. The radiation energy source terms and surface heat transfer, under the tangent slab approximation, are included within the fluid dynamic ow solver. The Fire II flight test, at the Mach-31 1643-second trajectory point, is used as a demonstration case. Comparisons are made with an existing structured-grid capability, the LAURA/HARA coupling. The radiative surface heat transfer rates from the present approach match the benchmark values within 6%. Although radiation coupling is the focus of the present work, convective surface heat transfer rates are also reported, and are seen to vary depending upon the choice of mesh connectivity and FUN3D ux reconstruction algorithm. On a tetrahedral-element mesh the convective heating matches the benchmark at the stagnation point, but under-predicts by 15% on the Fire II shoulder. Conversely, on a mixed-element mesh the convective heating over-predicts at the stagnation point by 20%, but matches the benchmark away from the stagnation region.

Wood, William A.↗
NASA NTRS2012

Performance Impact of Deflagration to Detonation Transition Enhancing Obstacles

A sub-model is developed to account for the drag and heat transfer enhancement resulting from deflagration-to-detonation (DDT) inducing obstacles commonly used in pulse detonation engines (PDE). The sub-model is incorporated as a source term in a time-accurate, quasi-onedimensional, CFD-based PDE simulation. The simulation and sub-model are then validated through comparison with a particular experiment in which limited DDT obstacle parameters were varied. The simulation is then used to examine the relative contributions from drag and heat transfer to the reduced thrust which is observed. It is found that heat transfer is far more significant than aerodynamic drag in this particular experiment.

Paxson, Daniel E.↗
NASA NTRS2012

Arctic Ocean Freshwater: How Robust are Model Simulations

The Arctic freshwater (FW) has been the focus of many modeling studies, due to the potential impact of Arctic FW on the deep water formation in the North Atlantic. A comparison of the hindcasts from ten ocean-sea ice models shows that the simulation of the Arctic FW budget is quite different in the investigated models. While they agree on the general sink and source terms of the Arctic FW budget, the long-term means as well as the variability of the FW export vary among models. The best model-to-model agreement is found for the interannual and seasonal variability of the solid FW export and the solid FW storage, which also agree well with observations. For the interannual and seasonal variability of the liquid FW export, the agreement among models is better for the Canadian Arctic Archipelago (CAA) than for Fram Strait. The reason for this is that models are more consistent in simulating volume flux anomalies than salinity anomalies and volume-flux anomalies dominate the liquid FW export variability in the CAA but not in Fram Strait. The seasonal cycle of the liquid FW export generally shows a better agreement among models than the interannual variability, and compared to observations the models capture the seasonality of the liquid FW export rather well. In order to improve future simulations of the Arctic FW budget, the simulation of the salinity field needs to be improved, so that model results on the variability of the liquid FW export and storage become more robust.

Jahn, A.↗
NASA NTRS2014

Metal Hydrides, MOFs, and Carbon Composites as Space Radiation Shielding Mitigators

Recently, metal hydrides and MOFs (Metal-Organic Framework/microporous organic polymer composites - for their hydrogen and methane storage capabilities) have been studied with applications in fuel cell technology. We have investigated a dual-use of these materials and carbon composites (CNT-HDPE) to include space radiation shielding mitigation. In this paper we present the results of a detailed study where we have analyzed 64 materials. We used the Band fit spectra for the combined 19-24 October 1989 solar proton events as the input source term radiation environment. These computational analyses were performed with the NASA high energy particle transport/dose code HZETRN. Through this analysis we have identified several of the materials that have excellent radiation shielding properties and the details of this analysis will be discussed further in the paper.

Atwell, William↗
NASA NTRS2013

On Spurious Numerics in Solving Reactive Equations

The objective of this study is to gain a deeper understanding of the behavior of high order shock-capturing schemes for problems with stiff source terms and discontinuities and on corresponding numerical prediction strategies. The studies by Yee et al. (2012) and Wang et al. (2012) focus only on solving the reactive system by the fractional step method using the Strang splitting (Strang 1968). It is a common practice by developers in computational physics and engineering simulations to include a cut off safeguard if densities are outside the permissible range. Here we compare the spurious behavior of the same schemes by solving the fully coupled reactive system without the Strang splitting vs. using the Strang splitting. Comparison between the two procedures and the effects of a cut off safeguard is the focus the present study. The comparison of the performance of these schemes is largely based on the degree to which each method captures the correct location of the reaction front for coarse grids. Here "coarse grids" means standard mesh density requirement for accurate simulation of typical non-reacting flows of similar problem setup. It is remarked that, in order to resolve the sharp reaction front, local refinement beyond standard mesh density is still needed.

Numerics↗
NASA NTRS2013

The North-South Asymmetry of the Heliospheric Current Sheet: Results of an MHD Simulation

A displacement of the heliospheric current sheet (HCS) south of the helioequator by approx.10deg was proposed by Simpson et al. (1996) as a possible explanation of the north-south asymmetry in the galactic cosmic rays observed by Ulysses during its first fast transit in 1994-1995. The idea was not supported by magnetic field measurements on Ulysses and, on this ground, was dismissed by Simpson et al. (1996). In addition, Erdos & Balogh (1998) argued that any north-south symmetry was unlikely as there should be flux balance between the magnetic sectors of opposite polarity. Nonetheless, many in the scientific community have accepted the original suggestion of Simpson et al. (1996) that a displacement of the HCS was responsible for the cosmic ray asymmetry. In this paper, using a magnetohydrodynamic model of the solar corona and solar wind that includes both dipole and quadrupole magnetic source terms, we show that a north-south asymmetry of the magnetic field on the Sun does not give rise to a displacement of the HCS. The lack of displacement of the HCS results from a latitudinal redistribution of magnetic flux near the Sun where the plasma beta much < 1. The latitudinal redistribution is a direct consequence of the magnetic field gradient between pole and equator. Near the Sun, the latitudinal gradient in magnetic field generates meridional flows directed equatorward that tend to relax the gradient in the magnetic field (to make it more latitude-independent) as heliocentric distance increases. If there is an asymmetry between north and south magnetic field strength then the meridional flows are also asymmetric (i.e., stronger in the hemisphere of stronger magnetic field). Because the magnetic fluxes (positive and negative) in the hemispheres must be equal, the redistribution shifts the HCS into balance by approx. 16 R(solar mass). At larger distances, where the magnetic field is relatively weak (beta much > 1), the HCS can be displaced if there is a difference in total pressure between the hemispheres.

solar wind↗
NASA NTRS2013

The Scaling of Broadband Shock-Associated Noise with Increasing Temperature

A physical explanation for the saturation of broadband shock-associated noise (BBSAN) intensity with increasing jet stagnation temperature has eluded investigators. An explanation is proposed for this phenomenon with the use of an acoustic analogy. To isolate the relevant physics, the scaling of BBSAN peak intensity level at the sideline observer location is examined. The equivalent source within the framework of an acoustic analogy for BBSAN is based on local field quantities at shock wave shear layer interactions. The equivalent source combined with accurate calculations of the propagation of sound through the jet shear layer, using an adjoint vector Green's function solver of the linearized Euler equations, allows for predictions that retain the scaling with respect to stagnation pressure and allows for saturation of BBSAN with increasing stagnation temperature. The sources and vector Green's function have arguments involving the steady Reynolds- Averaged Navier-Stokes solution of the jet. It is proposed that saturation of BBSAN with increasing jet temperature occurs due to a balance between the amplication of the sound propagation through the shear layer and the source term scaling.

Miller, Steven A. E.↗