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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Coarse Grid Modeling of Turbine Film Cooling Flows Using Volumetric Source Terms

The recent trend in numerical modeling of turbine film cooling flows has been toward higher fidelity grids and more complex geometries. This trend has been enabled by the rapid increase in computing power available to researchers. However, the turbine design community requires fast turnaround time in its design computations, rendering these comprehensive simulations ineffective in the design cycle. The present study describes a methodology for implementing a volumetric source term distribution in a coarse grid calculation that can model the small-scale and three-dimensional effects present in turbine film cooling flows. This model could be implemented in turbine design codes or in multistage turbomachinery codes such as APNASA, where the computational grid size may be larger than the film hole size. Detailed computations of a single row of 35 deg round holes on a flat plate have been obtained for blowing ratios of 0.5, 0.8, and 1.0, and density ratios of 1.0 and 2.0 using a multiblock grid system to resolve the flows on both sides of the plate as well as inside the hole itself. These detailed flow fields were spatially averaged to generate a field of volumetric source terms for each conservative flow variable. Solutions were also obtained using three coarse grids having streamwise and spanwise grid spacings of 3d, 1d, and d/3. These coarse grid solutions used the integrated hole exit mass, momentum, energy, and turbulence quantities from the detailed solutions as volumetric source terms. It is shown that a uniform source term addition over a distance from the wall on the order of the hole diameter is able to predict adiabatic film effectiveness better than a near-wall source term model, while strictly enforcing correct values of integrated boundary layer quantities.

Heidmann, James D.↗
NASA NTRS1975

Influence of grazing flow on duct wall normal impedances

A theory is developed for evaluating the normal impedance of an orifice connecting a cylindrical cavity to a flow duct. The impedance of a piston velocity source at the orifice is calculated in terms of the cavity modes which are modified by the presence of a vortex flow. The strength of the vortex is related to the mean flow in the duct. Results are presented showing the variation of the normal impedance spectrum with the cavity vortex. The trends of these variations are confirmed by comparison with measurements.

Mungur, P.↗
NASA NTRS1995

Numerical investigations of low-density nozzle flow by solving the Boltzmann equation

A two-dimensional finite-difference code to solve the BGK-Boltzmann equation has been developed. The solution procedure consists of three steps: (1) transforming the BGK-Boltzmann equation into two simultaneous partial differential equations by taking moments of the distribution function with respect to the molecular velocity u(sub z), with weighting factors 1 and u(sub z)(sup 2); (2) solving the transformed equations in the physical space based on the time-marching technique and the four-stage Runge-Kutta time integration, for a given discrete-ordinate. The Roe's second-order upwind difference scheme is used to discretize the convective terms and the collision terms are treated as source terms; and (3) using the newly calculated distribution functions at each point in the physical space to calculate the macroscopic flow parameters by the modified Gaussian quadrature formula. Repeating steps 2 and 3, the time-marching procedure stops when the convergent criteria is reached. A low-density nozzle flow field has been calculated by this newly developed code. The BGK Boltzmann solution and experimental data show excellent agreement. It demonstrated that numerical solutions of the BGK-Boltzmann equation are ready to be experimentally validated.

Deng, Zheng-Tao↗
NASA NTRS2013

Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Yee, Helen M. C.↗
NASA NTRS1978

Photoelectrons in the upper atmosphere - A formulation incorporating effects of transport

A technique for solving the problem of photoelectron spatial transport and spectral degradation is presented. The developed model, which allows for a relatively complete description of the relevant transport and collisional physics, can be applied to study processes at either low or high altitude, and may be used to investigate interhemispheric transport. The technique involves dividing the spatial dimension into a set of discrete ordinates and calculating the photoelectron intensity for a specific energy at each point in space and at a set of angle points. The selected range of the energy spectrum was covered completely by a set of energy bins, and the photoelectron intensity in the highest energy bin was calculated first by assuming that the only source term was derived from photoionization. The method is similar to a method used in radiative transfer theory. Spatial and angular redistribution of electrons play a dominant role at heights above 400 km.

Prather, M. J.↗
NASA NTRS1990

The sound of moving bodies

The importance of the quadrupole source term in the Ffowcs, Williams, and Hawkings (FWH) equation was addressed. The quadrupole source contains fundamental components of the complete fluid mechanics problem, which are ignored only at the risk of error. The results made it clear that any application of the acoustic analogy should begin with all of the source terms in the FWH theory. The direct calculation of the acoustic field as part of the complete unsteady fluid mechanics problem using CFD is considered. It was shown that aeroelastic calculation can indeed be made with CFD codes. The results indicate that the acoustic field is the most susceptible component of the computation to numerical error. Therefore, the ability to measure the damping of acoustic waves is absolutely essential both to develop acoustic computations. Essential groundwork for a new approach to the problem of sound generation by moving bodies is presented. This new computational acoustic approach holds the promise of solving many problems hitherto pushed aside.

Brentner, Kenneth Steven↗
NASA NTRS1992

A k-omega-multivariate beta PDF for supersonic combustion

In an attempt to study the interaction between combustion and turbulence in supersonic flows, an assumed PDF has been employed. This makes it possible to calculate the time average of the chemical source terms that appear in the species conservation equations. In order to determine the averages indicated in an equation, two transport equations, one for the temperature (enthalpy) variance and one for Q, are required. Model equations are formulated for such quantities. The turbulent time scale controls the evolution. An algebraic model similar to that used by Eklund et al was used in an attempt to predict the recent measurements of Cheng et al. Predictions were satisfactory before ignition but were less satisfactory after ignition. One of the reasons for this behavior is the inadequacy of the algebraic turbulence model employed. Because of this, the objective of this work is to develop a k-omega model to remedy the situation.

Alexopoulos, G. A.↗
NASA NTRS1989

Determination of near and far field acoustics for advanced propeller configurations

A method has been studied for predicting the acoustic field of the SR-3 transonic propfan using flow data generated by two versions of the NASPROP-E computer code. Since the flow fields calculated by the solvers include the shock-wave system of the propeller, the nonlinear quadrupole noise source term is included along with the monopole and dipole noise sources in the calculation of the acoustic near field. Acoustic time histories in the near field are determined by transforming the azimuthal coordinate in the rotating, blade-fixed coordinate system to the time coordinate in a nonrotating coordinate system. Fourier analysis of the pressure time histories is used to obtain the frequency spectra of the near-field noise.

Korkan, K. D.↗
NASA NTRS1988

The exact calculation of quadrupole sources for some incompressible flows

This paper is concerned with the application of the acoustic analogy of Lighthill to the acoustic and aerodynamic problems associated with moving bodies. The Ffowcs Williams-Hawkings equation, which is an interpretation of the acoustic analogy for sound generation by moving bodies, manipulates the source terms into surface and volume sources. Quite often in practice the volume sources, or quadrupoles, are neglected for various reasons. Recently, Farassat, Long and others have attempted to use the FW-H equation with the quadrupole source and neglected to solve for the surface pressure on the body. The purpose of this paper is to examine the contribution of the quadrupole source to the acoustic pressure and body surface pressure for some problems for which the exact solution is known. The inviscid, incompressible, 2-D flow, calculated using the velocity potential, is used to calculate the individual contributions of the various surface and volume source terms in the FW-H equation. The relative importance of each of the sources is then assessed.

Brentner, Kenneth S.↗
NASA NTRS1989

Acoustic design considerations: Review of rotor acoustic sources

It is not sufficient to optimize a rotor design in terms of a single noise level calculated for a single flight condition and a single measurement location. The various noise sources, their frequency content, amplitude, and directivity as a function of operating condition must be considered. A summary of the frequency ranges, directivity patterns and the most important operational and design parameters for major rotor noise sources is presented. It is difficult to generalize design requirements for rotor noise because the acoustic output varies so widely depending on the noise source, flight condition, measurement location, and frequency range. However, assuming the rotor must lift a fixed nominal payload and operate over a wide range of flight conditions, three general design guidelines can be stated: (1) minimize tip Mach number; (2) minimize blade thickness in the tip region; and (3) minimize gradients in the spanwise lift distribution in the tip region. Constraints on blade thickness, maximum values for hover tip Mach number, advancing tip Mach number and spanwise lift coefficient gradient will be specified during the aerodynamic, dynamic and structural optimization process. The rotor noise sources to be considered include the low frequency loading and thickness noise, and the higher frequency noise due to blade-vortex interactions (BVI). The analyses to be employed will include the comprehensive rotor analysis and design program CAMRAD and the rotor noise prediction program WOPWOP.

Martin, Ruth M.↗
NASA NTRS1979

Precipitation fluxes of energetic electrons at Jupiter - An estimated upper limit

Divine's (1976) model for the observed energetic electron fluxes in the inner Jovian magnetosphere is used to calculate space densities, and the combined source and loss term is evaluated for steady-state radial diffusion conserving the first and second adiabatic invariants. Upper-limit estimates for the ionospheric precipitation fluxes at various values of L and over various electron energy ranges are derived by assuming that contributions to the source term are negligible and that particle precipitation due to pitch-angle diffusion in the loss cone is the main contributor to the loss term. The energy deposition rates corresponding to these precipitated particles are estimated, along with the bremsstrahlung X-ray emissions that should result. It is found that the X-ray emissions should be 5 to 6 orders of magnitude below the reported upper observational limit.

Thomsen, M. F.↗
NASA NTRS1997

Computation of Reacting Flows in Combustion Processes

The main objective of this research was to develop an efficient three-dimensional computer code for chemically reacting flows. The main computer code developed is ALLSPD-3D. The ALLSPD-3D computer program is developed for the calculation of three-dimensional, chemically reacting flows with sprays. The ALL-SPD code employs a coupled, strongly implicit solution procedure for turbulent spray combustion flows. A stochastic droplet model and an efficient method for treatment of the spray source terms in the gas-phase equations are used to calculate the evaporating liquid sprays. The chemistry treatment in the code is general enough that an arbitrary number of reaction and species can be defined by the users. Also, it is written in generalized curvilinear coordinates with both multi-block and flexible internal blockage capabilities to handle complex geometries. In addition, for general industrial combustion applications, the code provides both dilution and transpiration cooling capabilities. The ALLSPD algorithm, which employs the preconditioning and eigenvalue rescaling techniques, is capable of providing efficient solution for flows with a wide range of Mach numbers. Although written for three-dimensional flows in general, the code can be used for two-dimensional and axisymmetric flow computations as well. The code is written in such a way that it can be run in various computer platforms (supercomputers, workstations and parallel processors) and the GUI (Graphical User Interface) should provide a user-friendly tool in setting up and running the code.

Keith, Theo G., Jr.↗
NASA NTRS1983

A series expansion of the acoustic power radiated from planar sources

A series expansion in ascending powers of the wavenumber k is derived for the acoustic power delivered by baffled or unbaffled planar sources. This series provides a relatively simple means of derving expressions for the power radiated by a baffled source with a known velocity distribution and can be used for unbaffled plates when the velocity field outside the plate is also known. The terms in the series are calculated from the moments of this velocity distribution in the plane containing the source. If these moments are written as derivaties in wavenumber space, it is shown that a MacLaurin expansion of the Fourier transformed velocity provides an easy technique for computing the first few terms of the acoustic power. Examples are provided for baffled, rectangular plates with various boundary conditions. The arbirarily shaped plate with free boundaries is particularly interesting. It is proven that the volume flow across it surface must be zero and as a result corner and edge mode radiation cannot exist for this kind of source.

Willams, E. G.↗
NASA NTRS1994

'Spherical' disks: Moving toward a unified source model for L1551

To predict the effects of a disk on the spectral energy distribution of a deeply embedded protostar, we construct disk models with power-law temperature distributions T is proportional to r(exp -q). We then use the spherically averaged disk emission as the central source for a spherical envelope, hence the term, 'spherical' disk. We then calculate the predicted spectral energy distribution of the disk and envelope, using a spherically symmetric radiative transport code. Applying this procedure to L1551 IRS 5, we find that the predicted far-infrared flux is not very sensitive to the nature of the central source. The best source model is consistent with the far-infrared emission arising from the infalling region in an 'inside-out' collapse model, independent of the nature of the central source. Disk models are superior to the star-only model when we try to match millimeter interferometer data. While disks with various q can reproduce the observed 2.7 mm interferometer flux, only an active disk (q = 0.5) can produce enough emission in a region small enough to match the observed 2.7 mm visibilities. However, if the disk is backwarmed by the envelope, even purely reprocessing disks can meet this constraint. All types of backwarmed disks are virtually in distinguishable in their millimeter properties. We find that all reasonable envelope models are sufficiently opaque in the mid-infrared to attenuate any disk model to a level well below the observations, unless the ratio of the mid-infrared to far-infrared dust opacities is similar to that of the dust opacities advocated by Mathis, Mezger, & Panagia (1983).

Butner, Harold M.↗
NASA NTRS1984

Implications of Observed Frequencies of Solar P Modes

A preliminary comparison is made of the observed frequencies of 5 min P modes with theoretical frequencies for a traditional solar model. The differences between observations and theory can be understood qualitatively in terms of two separate sources of error in the frequency calculation, one near the solar surface and the other at the base of the convection zone. There is no indication of errors in the deep interior of the model.

Christensen-Dalsgaard, J.↗
NASA NTRS1993

A nonequilibrium model for a moderate pressure hydrogen microwave discharge plasma

This document describes a simple nonequilibrium energy exchange and chemical reaction model to be used in a computational fluid dynamics calculation for a hydrogen plasma excited by microwaves. The model takes into account the exchange between the electrons and excited states of molecular and atomic hydrogen. Specifically, electron-translation, electron-vibration, translation-vibration, ionization, and dissociation are included. The model assumes three temperatures, translational/rotational, vibrational, and electron, each describing a Boltzmann distribution for its respective energy mode. The energy from the microwave source is coupled to the energy equation via a source term that depends on an effective electric field which must be calculated outside the present model. This electric field must be found by coupling the results of the fluid dynamics and kinetics solution with a solution to Maxwell's equations that includes the effects of the plasma permittivity. The solution to Maxwell's equations is not within the scope of this present paper.

Scott, Carl D.↗
NASA NTRS1991

Systems engineering studies of lunar base construction

Many ingenious concepts have been proposed for lunar base construction, but few systematic studies exist which relate time-consistent lunar base construction technologies and the choice of lunar base approach with the long-term SEI objectives - i.e., lunar indigenous base construction and Mars Exploration equipment development. To fill this gap, CSC has taken a two-pronged approach. First, the Center undertook basic geotechnical investigations of lunar soil, fabrication of a scale prototype of a lunar construction crane, a multi-robot construction team laboratory experiment, and a preliminary design of lunar base structures. Second, during Jun. and Jul. 1991 two lunar base construction systems engineering studies were accomplished - a 'near term lunar base' study, and a 'far-term lunar base' study. The goals of these studies were to define the major lunar base construction research problems in consistent technology/construction frameworks, and to define design requirements for construction equipment such as a lunar crane and a regolith mover. The 'near-term lunar base' study examined three different construction concepts for a lunar base comprised of pre-fabricated, pre-tested, Space Station Freedom-type modules, which would be covered with regolith shielding. Concept A used a lunar crane for unloading and transportation; concept B, a winch and cart; and concept C, a walker to move the modules from the landing site to the base site and assemble them. To evaluate the merits of each approach, calculations were made of mass efficiency measure, source mass, reliability, far-term base mass, Mars base mass, and base assembly time. The model thus established was also used to define the requirements for crane speed and regolith mover m(sup 3)/sec rates. A major problem addressed is how to 'mine' the regolith and stack it over the habitats as shielding. To identify when the cost of using indigenous lunar materials to construct the base exceeds the cost of development and delivery of the equipment for processing lunar materials, a study of construction of a candidate sintered regolith 'far term lunar base' was undertaken. A technique was devised for casting slabs of sintered (basaltic) regolith and assembling these into a hemispherical (or geodesic) dome. The major problem occurs with the inner liner. At 14.7 psi and 20 percent oxygen internal atmosphere, the entire structure is in tension, even with the regolith load. Also, another study has indicated that at 14.7 psi major resupply of air will be needed because of leakage, and astronauts may have to engage in extensive pre-breathing and post-breathing for extravehicular activity (EVA) tasks, thus detracting from useful mission work time. An alternative is to operate part of the base at, say, 5 psi and 70 percent oxygen, or to equip the astronauts with hard suits at 8.3 psi or greater. All of these choices directly influence base design and construction techniques.

Morgenthaler, George W.↗
NASA NTRS1987

Ring vortex/cylinder sound production revisited

The acoustic radiation produced by the passage of a vortex ring near a circular cylinder is calculated utilizing the Coriolis acceleration of the vortex as the source term and the low frequency Green's function approach. Necessary integrals are evaluated exactly via complex integration techniques. The results exhibit good agreement (on a Pascal basis) with experimental data. Of particular interest is the rotation of the directivity pattern during this transient encounter and the extreme sensitivity of the results to the distance at which the vortex passes by the cylinder.

Hardin, Jay C.↗