Search NASA⌕ Search

SEARCH · Search NASA

Results for “source term”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Magnetic dipole moment of a spherical shell with TRM acquired in a field of internal origin

The acquisition of thermoremanent magnetization (TRM) by a cooling spherical shell is studied for internal magnetizing dipole fields, using Runcorn's (1975) theorems on magnetostatics. If the shell cools progressively inward, inner regions acquire TRM in a net field composed of the dipole source term plus a uniform field due to the outer magnetized layers. In this case, the global dipole moment and external remanent field are nonzero when the whole shell has cooled below the Curie point and the source dipole has disappeared. The remanent field outside the shell is found to depend on the thickness, radii, and cooling rate of the shell, as well as the coefficient of TRM and the intensity of the magnetizing field. Some implications for the moon's remanent dipole moment are discussed.

Srnka, L. J.↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
NASA NTRS2006

Land Water Storage within the Congo Basin Inferred from GRACE Satellite Gravity Data

GRACE satellite gravity data is used to estimate terrestrial (surface plus ground) water storage within the Congo Basin in Africa for the period of April, 2002 - May, 2006. These estimates exhibit significant seasonal (30 +/- 6 mm of equivalent water thickness) and long-term trends, the latter yielding a total loss of approximately 280 km(exp 3) of water over the 50-month span of data. We also combine GRACE and precipitation data set (CMAP, TRMM) to explore the relative contributions of the source term to the seasonal hydrological balance within the Congo Basin. We find that the seasonal water storage tends to saturate for anomalies greater than 30-44 mm of equivalent water thickness. Furthermore, precipitation contributed roughly three times the peak water storage after anomalously rainy seasons, in early 2003 and 2005, implying an approximately 60-70% loss from runoff and evapotranspiration. Finally, a comparison of residual land water storage (monthly estimates minus best-fitting trends) in the Congo and Amazon Basins shows an anticorrelation, in agreement with the 'see-saw' variability inferred by others from runoff data.

Crowley, John W.↗
NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗
NASA NTRS1985

Interpretation of ozone temperature correlations. II - Analysis of SBUV ozone data

In an effort to understand the behavior of ozone in the transition region, the terms in the perturbation continuity equation have been evaluated from solar backscattered ultraviolet ozone data, National Meteorological Center temperature data, and geostrophic winds derived from the temperature data at 2 mbar, winter 1978-1979. A linear parameterization is used to estimate the chemical source term. The perturbation ozone field calculated by integrating the net forcing evaluated from the data compares well with the measured perturbation ozone field. There are occasions during the winter when the sum of the dynamic terms is small due to cancellation, and the ozone and temperature perturbations are anticorrelated. The -1 correlation between the ozone and temperature perturbations does not necessarily imply that chemical forcing is dominant.

Douglass, A. R.↗
NASA NTRSDate not supplied

Integral Velocity Sampling for Unsteady Rotor Models on Cartesian Meshes

An unsteady body-force rotor model utilizing integral velocity sampling is developed for Cartesian meshes with a focus on accurate performance predictions. In this approach, the rotors are modeled as body-force source terms, and integral sampling is used to define the local freestream velocity vector. The method is verified in two dimensions for airfoils, providing accurate angle of attack and freestream Mach number estimates even on relatively coarse meshes, and demonstrating excellent mesh convergence across a variety of source region thicknesses. Extensions of the model to three dimensions for both fixed and rotating wings are considered, with numerical experiments demonstrating the need for a spanwise-continuous representation of the induced velocity or angle of attack. The model provides reasonable spanwise lift distributions without any explicit tip-loss model. Comparisons to high-fidelity simulations and experimental data show excellent agreement for predictions of an isolated DJI Phantom 3 propeller in hover. Finally, simulations of the full aircraft in edgewise forward flight also compare well with wind tunnel data across a variety of pitch angles and propeller speeds.

ARMD↗
NASA NTRS1979

Scattering and distortion of the unsteady motion on transversely sheared mean flows

It is shown that the pressure and velocity fluctuations of the unsteady motion on a transversely sheared mean flow can be expressed entirely in terms of the derivatives of two potential functions. One of these is a convected quantity that can be specified as a boundary condition and is related to a transverse component of the upstream velocity field. The other can be determined by solving an inhomogeneous wave equation whose source term is also a convected quantity that can be specified as a boundary condition in any given problem. The general theory is used to study the interaction of an unsteady flow with a semi-infinite plate embedded in a shear layer. The acoustic field produced by this interaction is calculated in the limits of low and high frequency. The results are compared with experimental one-third octave sound pressure level radiation patterns. The agreement is found to be excellent, especially in the low frequency range, where the mean-flow and convective effects are shown to have a strong influence on the directivity of the sound.

Goldstein, M. E.↗
NASA NTRS1993

Examination of nonequilibrium effects in an ionized nitrogen flow

A thermo-chemical nonequilibrium ionized nitrogen flow in a shock tube is examined. A one-dimensional computational fluid dynamics code has been developed to study the flowfield incorporating a novel treatment of the electron energy equation. In the present approach, the electron pressure is included in the electron energy flux. In contrast, previous work has removed this term from the flux to facilitate the flux splitting. The code uses a five-species gas model for ionized nitrogen characterized by translational-rotational, vibrational, and electron-electronic temperatures. The results give good agreement with experimental data except for an inadequacy in the electronic source terms. A small difference is observed between the new and old splittings which increases with higher levels of ionization.

Hatfield, John A.↗
NASA NTRS2003

Various Numerical Applications on Tropical Convective Systems Using a Cloud Resolving Model

In recent years, increasing attention has been given to cloud resolving models (CRMs or cloud ensemble models-CEMs) for their ability to simulate the radiative-convective system, which plays a significant role in determining the regional heat and moisture budgets in the Tropics. The growing popularity of CRM usage can be credited to its inclusion of crucial and physically relatively realistic features such as explicit cloud-scale dynamics, sophisticated microphysical processes, and explicit cloud-radiation interaction. On the other hand, impacts of the environmental conditions (for example, the large-scale wind fields, heat and moisture advections as well as sea surface temperature) on the convective system can also be plausibly investigated using the CRMs with imposed explicit forcing. In this paper, by basically using a Goddard Cumulus Ensemble (GCE) model, three different studies on tropical convective systems are briefly presented. Each of these studies serves a different goal as well as uses a different approach. In the first study, which uses more of an idealized approach, the respective impacts of the large-scale horizontal wind shear and surface fluxes on the modeled tropical quasi-equilibrium states of temperature and water vapor are examined. In this 2-D study, the imposed large-scale horizontal wind shear is ideally either nudged (wind shear maintained strong) or mixed (wind shear weakened), while the minimum surface wind speed used for computing surface fluxes varies among various numerical experiments. For the second study, a handful of real tropical episodes (TRMM Kwajalein Experiment - KWAJEX, 1999; TRMM South China Sea Monsoon Experiment - SCSMEX, 1998) have been simulated such that several major atmospheric characteristics such as the rainfall amount and its associated stratiform contribution, the Qlheat and Q2/moisture budgets are investigated. In this study, the observed large-scale heat and moisture advections are continuously applied to the 2-D model. The modeled cloud generated from such an approach is termed continuously forced convection or continuous large-scale forced convection. A third study, which focuses on the respective impact of atmospheric components on upper Ocean heat and salt budgets, will be presented in the end. Unlike the two previous 2-D studies, this study employs the 3-D GCE-simulated diabatic source terms (using TOGA COARE observations) - radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the Ocean mixed-layer (OML) model.

Shie, C.-L.↗
NASA NTRS2017

Identification of Spurious Signals from Permeable Ffowcs Williams and Hawkings Surfaces

Integral forms of the permeable surface formulation of the Ffowcs Williams and Hawkings (FW-H) equation often require an input in the form of a near field Computational Fluid Dynamics (CFD) solution to predict noise in the near or far field from various types of geometries. The FW-H equation involves three source terms; two surface terms (monopole and dipole) and a volume term (quadrupole). Many solutions to the FW-H equation, such as several of Farassat's formulations, neglect the quadrupole term. Neglecting the quadrupole term in permeable surface formulations leads to inaccuracies called spurious signals. This paper explores the concept of spurious signals, explains how they are generated by specifying the acoustic and hydrodynamic surface properties individually, and provides methods to determine their presence, regardless of whether a correction algorithm is employed. A potential approach based on the equivalent sources method (ESM) and the sensitivity of Formulation 1A (Formulation S1A) is also discussed for the removal of spurious signals.

Lopes, Leonard V.↗
NASA NTRS2004

High-Resolution Simulation of Hurricane Bonnie (1998): Water Budget - Part 2

The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weak updrafts, indicative of a dominance of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where. cooling from melting enhances condensation. Unlike squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance stratiform precipitation development through additional growth by vapor deposition and aggregation. In addition, artificial source terms for cloud and precipitation mass associated with setting negative mixing ratios, produced as a result of numerical advection errors, to zero are described. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Braun, Scott A.↗
NASA NTRS2005

The Water Budget of a Simulated Hurricane

The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weaker updrafts, indicative of a larger role of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where cooling from melting enhances condensation. Unlike in squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance in situ stratiform precipitation development through additional growth by vapor deposition and aggregation. We examine artificial source terms for cloud and precipitation mass associated with setting to zero negative mixing ratios that arise from numerical advection errors. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation and 15-20% of the precipitation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Braun, S.↗
NASA NTRS2001

Water Vapor Tracers as Diagnostics of the Regional Hydrologic Cycle

Numerous studies suggest that local feedback of surface evaporation on precipitation, or recycling, is a significant source of water for precipitation. Quantitative results on the exact amount of recycling have been difficult to obtain in view of the inherent limitations of diagnostic recycling calculations. The current study describes a calculation of the amount of local and remote geographic sources of surface evaporation for precipitation, based on the implementation of three-dimensional constituent tracers of regional water vapor sources (termed water vapor tracers, WVT) in a general circulation model. The major limitation on the accuracy of the recycling estimates is the veracity of the numerically simulated hydrological cycle, though we note that this approach can also be implemented within the context of a data assimilation system. In the WVT approach, each tracer is associated with an evaporative source region for a prognostic three-dimensional variable that represents a partial amount of the total atmospheric water vapor. The physical processes that act on a WVT are determined in proportion to those that act on the model's prognostic water vapor. In this way, the local and remote sources of water for precipitation can be predicted within the model simulation, and can be validated against the model's prognostic water vapor. As a demonstration of the method, the regional hydrologic cycles for North America and India are evaluated for six summers (June, July and August) of model simulation. More than 50% of the precipitation in the Midwestern United States came from continental regional sources, and the local source was the largest of the regional tracers (14%). The Gulf of Mexico and Atlantic regions contributed 18% of the water for Midwestern precipitation, but further analysis suggests that the greater region of the Tropical Atlantic Ocean may also contribute significantly. In most North American continental regions, the local source of precipitation is correlated with total precipitation. There is a general positive correlation between local evaporation and local precipitation, but it can be weaker because large evaporation can occur when precipitation is inhibited. In India, the local source of precipitation is a small percentage of the precipitation owing to the dominance of the atmospheric transport of oceanic water. The southern Indian Ocean provides a key source of water for both the Indian continent and the Sahelian region.

Bosilovich, Michael G.↗
NASA NTRS2021

Preventing Anomalous Torques in Circumbinary Accretion Simulations

Numerical experiments are the primary method of studying the evolution of circumbinary disks due to the strong nonlinearities involved. Many circumbinary simulations also require the use of numerical mass sinks: source terms that prevent gas from unphysically accumulating around the simulated point masses by removing gas at a given rate. However, special care must be taken when drawing physical conclusions from such simulations to ensure that results are not biased by numerical artifacts. We demonstrate how improved sink methods reduce some of these biases, using simulations of aspect ratio 0.1 accretion disks around binaries with mass ratios between 0.1 and 1. We show that sink terms that do not reduce the angular momentum of gas relative to the accreting object (1) reduce the dependence on the sink rate of quantities such as the torque on the binary, distribution of accretion between binary components, and evolution of the binary semimajor axis; (2) reduce the degree to which the sink rate affects the structure of the accretion disks around each binary component; (3) alter the variability of accretion onto the binary, regularizing it in time. We also investigate other potential sources of systematic error, such as gravitational softening and simplifications to the viscous stress tensor. Because of the strong dependence of binary orbital evolution on both the torque and the distribution of mass between binary components, the sink method influences the orbital evolution of the binary at all mass ratios, with mass ratios below ∼0.3 most affected.

Alexander Dittmann↗
NASA NTRS1996

Theoretical Basis for Finite Difference Extrapolation of Sonic Boom Signatures

Calculation of sonic boom signatures for aircraft has traditionally followed the methods of Whitham' and Walkden. The wave disturbance generated by the vehicle is obtained by area rule linearized supersonic flow methods, which yield a locally axisymmetric asymptotic solution. This solution is acoustic in nature, i.e., first order in disturbance quantities, and corresponds to ray acoustics. Cumulative nonlinear distortion of the signature is incorporated by using this solution to adjust propagation speed to first order, thus yielding a solution second order in disturbance quantities. The effects of atmospheric gradients are treated by Blokhintzov's method of geometrical acoustics. Both nonlinear signature evolution and ray tracing are applied as if the pressure field very close to the vehicle were actually that given by the source term (the 'F-function') of the asymptotic linearized flow solution. The viewpoint is thus that the flow solution exists at a small radius near the vehicle, and may be treated as an input to an extrapolation procedure consisting of ray tracing and nonlinear aging. The F-function is often regarded as a representation of a near-field pressure signature, and it is common for computational implementations to treat it interchangeably with the pressure signature. There is a 'matching radius' between the source function and the subsequent propagation extrapolation. This viewpoint has been supported by wind tunnel tests of simple models, and very typically yields correct results for actual flight vehicles. The assumption that the F-function and near-field signature are interchangeable is generally not correct. The flowfield of a vehicle which is not axisymmetric contains crossflow components which are very significant at small radii and less so at larger distances. From an acoustical viewpoint, the crossflow is equivalent to source diffraction portions of the wave field. Use of the F-function as a near field signature effectively assumes that the diminution of the crossflow/diffraction component may be applied all at once at the matching radius noted above. This approximation, though not rigorously validated, is responsible for the usual correct far-field results. On the other hand, if an actual near-field signature (either from wind tunnel or CFD data) is used at a starting point rather than one based on th effective source distribution, the predicted far-field signature is generally wrong.

Plotkin, Kenneth J.↗
NASA NTRS2002

Modeling Aerodynamically Generated Sound of Helicopter Rotors

A great deal of progress has been made in the modeling of aerodynamically generated sound of rotors over the past decade. Although the modeling effort has focused on helicopter main rotors, the theory is generally valid for a wide range of rotor configurations. The Ffowcs Williams Hawkings (FW-H) equation has been the foundation for much of the development. The monopole and dipole source terms of the FW-H equation account for the thickness and loading noise, respectively. Bladevortex-interaction noise and broadband noise are important types of loading noise, hence much research has been directed toward the accurate modeling of these noise mechanisms. Both subsonic and supersonic quadrupole noise formulations have been developed for the prediction of high-speed impulsive noise. In an effort to eliminate the need to compute the quadrupole contribution, the FW-H equation has also been utilized on permeable surfaces surrounding all physical noise sources. Comparisons of the Kirchhoff formulation for moving surfaces with the FW-H equation have shown that the Kirchhoff formulation for moving surfaces can give erroneous results for aeroacoustic problems. Finally, significant progress has been made incorporating the rotor noise models into full vehicle noise prediction tools.

Brentner, Kenneth S.↗
NASA NTRS2018

High-resolution Wave Propagation Method for Stratified Flows

The implementation of the multidimensional f-waves Riemann solver for the time-dependent, three-dimensional, nonhydrostatic, meso- and microscale atmospheric flows is described in detail. The Riemann solver employs flux-based wave decomposition (f-waves) for the calculation of Godunov fluxes in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The scheme incorporates the source term due to gravity without introducing discretization errors which is an important property in the context of atmospheric flows. The resulting flow solver is conservative, accurate, stable, and well-balanced. The implementation of the solver is evaluated using benchmark test cases for atmospheric dynamics.

Riemann problem↗