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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

On the Computation of Sound by Large-Eddy Simulations

The effect of the small scales on the source term in Lighthill's acoustic analogy is investigated, with the objective of determining the accuracy of large-eddy simulations when applied to studies of flow-generated sound. The distribution of the turbulent quadrupole is predicted accurately, if models that take into account the trace of the SGS stresses are used. Its spatial distribution is also correct, indicating that the low-wave-number (or frequency) part of the sound spectrum can be predicted well by LES. Filtering, however, removes the small-scale fluctuations that contribute significantly to the higher derivatives in space and time of Lighthill's stress tensor T(sub ij). The rms fluctuations of the filtered derivatives are substantially lower than those of the unfiltered quantities. The small scales, however, are not strongly correlated, and are not expected to contribute significantly to the far-field sound; separate modeling of the subgrid-scale density fluctuations might, however, be required in some configurations.

Piomelli, Ugo↗
NASA NTRS1996

Computational Modeling of Multi-Phase/Multi-Species Flows with Applications to Liquid Rocket Engines

Accurate prediction of all physical phenomena in a combustion chamber is essential for better understanding of the system performance. Atomization, evaporation, combustion, chemical kinetics, and turbulence are those processes of great importance that need to be well understood. Processes involving the liquid phase in a combustion chamber will be further complicated under supercritical conditions. More advanced and accurate numerical techniques are required to extend our understanding of the above phenomena. A computer program for multi-species/multi-phase flow was developed for NASA/MSFC in 1992. This code, called Liquid Thrust Chamber Performance (LTCP) program takes an Eulerian- Eulerian approach and is based on the Total Variation Diminishing (TVD) technique with Lax-Friedrichs upwind method. Under the NASA/ASEE SFFP the LTCP code was used to predict the performance characteristics of several engines that were of particular interest to NASA. This code was also successful in a combustion detonation study. Converting the program to the PC platform was accomplished which extends usability and makes it available to a wider range of users. The Eulerian formulation of the liquid phase provides a suitable model that can be extended to include combustion modeling under supercritical conditions. The results have been compared against the ones of other codes and available measured data. The algorithm proved to be robust and efficient for problems with stiff source terms.

Navaz, Homayun K.↗
NASA NTRS1998

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗
NASA NTRS1999

Effects of G-Jitter on Directional Solidification of a Binary Alloy

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermo-solutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter. A background gravity of 1 microgram has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Benjapiyaporn, C.↗
NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗
NASA NTRS1998

Inner Magnetospheric Superthermal Electron Transport: Photoelectron and Plasma Sheet Electron Sources

Two time-dependent kinetic models of superthermal electron transport are combined to conduct global calculations of the nonthermal electron distribution function throughout the inner magnetosphere. It is shown that the energy range of validity for this combined model extends down to the superthermal-thermal intersection at a few eV, allowing for the calculation of the en- tire distribution function and thus an accurate heating rate to the thermal plasma. Because of the linearity of the formulas, the source terms are separated to calculate the distributions from the various populations, namely photoelectrons (PEs) and plasma sheet electrons (PSEs). These distributions are discussed in detail, examining the processes responsible for their formation in the various regions of the inner magnetosphere. It is shown that convection, corotation, and Coulomb collisions are the dominant processes in the formation of the PE distribution function and that PSEs are dominated by the interplay between the drift terms. Of note is that the PEs propagate around the nightside in a narrow channel at the edge of the plasmasphere as Coulomb collisions reduce the fluxes inside of this and convection compresses the flux tubes inward. These distributions are then recombined to show the development of the total superthermal electron distribution function in the inner magnetosphere and their influence on the thermal plasma. PEs usually dominate the dayside heating, with integral energy fluxes to the ionosphere reaching 10(exp 10) eV/sq cm/s in the plasmasphere, while heating from the PSEs typically does not exceed 10(exp 8) eV/sq cm/s. On the nightside, the inner plasmasphere is usually unheated by superthermal electrons. A feature of these combined spectra is that the distribution often has upward slopes with energy, particularly at the crossover from PE to PSE dominance, indicating that instabilities are possible.

Khazanov, G. V.↗
NASA NTRS1999

Streamwise Vorticity Generation in Laminar and Turbulent Jets

Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Generation mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical (DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching occurs, RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds number turbulent jets produce axis-switching even for cases with several instability modes in the mixing layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the instantaneous results shows that, complex interactions of gradients of the normal and shear Reynolds stresses are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations confirm these results. k - epsilon turbulence model computations fail to reproduce the phenomenon, whereas algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are computed explicitly, succeeded in reproducing the phenomenon accurately.

Demuren, Aodeji O.↗
NASA NTRS1999

Multistage Simulations of the GE90 Turbine

The average passage approach has been used to analyze three multistage configurations of the GE90 turbine. These are a high pressure turbine rig, a low pressure turbine rig and a full turbine configuration comprising 18 blade rows of the GE90 engine at takeoff conditions. Cooling flows in the high pressure turbine have been simulated using source terms. This is the first time a dual-spool cooled turbine has been analyzed in 3D using a multistage approach. There is good agreement between the simulations and experimental results. Multistage and component interaction effects are also presented. The parallel efficiency of the code is excellent at 87.3% using 121 processors on an SGI Origin for the 18 blade row configuration. The accuracy and efficiency of the calculation now allow it to be effectively used in a design environment so that multistage effects can be accounted for in turbine design.

Turner, Mark G.↗
NASA NTRS1998

Inner Magnetospheric Superthermal Electron Transport: Photoelectron and Plasma Sheet Electron Sources

Two time-dependent kinetic models of superthermal electron transport are combined to conduct global calculations of the nonthermal electron distribution function throughout the inner magnetosphere. It is shown that the energy range of validity for this combined model extends down to the superthermal-thermal intersection at a few eV, allowing for the calculation of the entire distribution function and thus an accurate heating rate to the thermal plasma. Because of the linearity of the formulas, the source terms are separated to calculate the distributions from the various populations, namely photoelectrons (PEs) and plasma sheet electrons (PSEs). These distributions are discussed in detail, examining the processes responsible for their formation in the various regions of the inner magnetosphere. It is shown that convection, corotation, and Coulomb collisions are the dominant processes in the formation of the PE distribution function, and that PSEs are dominated by the interplay between the drift terms. Of note is that the PEs propagate around the nightside in a narrow channel at the edge of the plasmasphere as Coulomb collisions reduce the fluxes inside of this and convection compresses the flux tubes inward. These distributions are then recombined to show the development of the total superthermal electron distribution function in the inner magnetosphere and their influence on the thermal plasma. PEs usually dominate the dayside heating, with integral energy fluxes to the ionosphere reaching 10(exp 10) eV/sq cm/s in the plasmasphere, while heating from the PSEs typically does not exceed 10(exp 8)eV/sq cm/s. On the nightside, the inner plasmasphere is usually unheated by superthermal electrons. A feature of these combined spectra is that the distribution often has upward slopes with energy, particularly at the crossover from PE to PSE dominance, indicating that instabilities are possible.

Khazanov, G. V.↗
NASA NTRS1999

A One Dimensional, Time Dependent Inlet/Engine Numerical Simulation for Aircraft Propulsion Systems

The NASA Lewis Research Center (LeRC) and the Arnold Engineering Development Center (AEDC) have developed a closely coupled computer simulation system that provides a one dimensional, high frequency inlet/engine numerical simulation for aircraft propulsion systems. The simulation system, operating under the LeRC-developed Application Portable Parallel Library (APPL), closely coupled a supersonic inlet with a gas turbine engine. The supersonic inlet was modeled using the Large Perturbation Inlet (LAPIN) computer code, and the gas turbine engine was modeled using the Aerodynamic Turbine Engine Code (ATEC). Both LAPIN and ATEC provide a one dimensional, compressible, time dependent flow solution by solving the one dimensional Euler equations for the conservation of mass, momentum, and energy. Source terms are used to model features such as bleed flows, turbomachinery component characteristics, and inlet subsonic spillage while unstarted. High frequency events, such as compressor surge and inlet unstart, can be simulated with a high degree of fidelity. The simulation system was exercised using a supersonic inlet with sixty percent of the supersonic area contraction occurring internally, and a GE J85-13 turbojet engine.

Doug Garrard↗
NASA NTRS1999

Numerical Modeling of Three-Dimensional Fluid Flow with Phase Change

We present a numerical method to compute phase change dynamics of three-dimensional deformable bubbles. The full Navier-Stokes and energy equations are solved for both phases by a front tracking/finite difference technique. The fluid boundary is explicitly tracked by discrete points that are connected by triangular elements to form a front that is used to keep the stratification of material properties sharp and to calculate the interfacial source terms. Two simulations are presented to show robustness of the method in handling complex phase boundaries. In the first case, growth of a vapor bubble in zero gravity is studied where large volume increase of the bubble is managed by adaptively increasing the front resolution. In the second case, growth of a bubble under high gravity is studied where indentation at the rear of the bubble results in a region of large curvature which challenges the front tracking in three dimensions.

Esmaeeli, Asghar↗
NASA NTRS2001

An Overview of In-Stu Treatability Studies at Marshall Space Flight Center, Huntsville, Alabama

Marshall Space Flight Center (MSFC) is located in Huntsville, Alabama (north-central Alabama), on approximately 1,840 acres near the center of the U.S. Army's Redstone Arsenal (RSA). MSFC is the National Aeronautics and Space Administration's (NASA's) principal propulsion development center. Its scientists, engineers, and support personnel play a major role in the National Space Transportation System by managing space shuttle mission activities, including the microgravity laboratory. In addition, MSFC will be a significant contributor to several of NASA's future programs, including the Reusable Launch Vehicle (X-33), International Space Station, and Advanced X-ray Astrophysics Facility, as well as research on a variety of space science applications. MSFC has been used to develop, test and manufacture space vehicles and components since 1960, when civilian rocketry and missile activities were transferred from RSA to MSFC. In 1994, MSFC was placed on the National Priority List for the management of hazardous waste sites, under the requirements of the Comprehensive Environmental Response, Compensation, and Liability Act (CERCLA). One requirement of the CERCLA program is to evaluate the nature and extent of environmental contamination resulting from identified CERCLA sites, assess the public health and environmental risks associated with the identified contamination, and identify potential remedial actions. A CERCLA remedial investigation (RI) for the groundwater system has identified at least five major plumes of chlorinated volatile organic compounds (CVOCs) in the groundwater beneath the facility. These plumes are believed to be the result of former management practices at 14 main facility locations (termed "source areas") where CVOCs were released to the subsurface. Trichloroethene (TCE) is the predominant CVOC and is common to all the plumes. Perchloroethene (PCE) also exists in two of the plumes. In addition to TCE and PCE, carbon tetrachloride and 1,1,2,2-tetrachloroethane are contained in one of the plumes. The CVOCs are believed to exist as dense non-aqueous phase liquids (DNAPLs) beneath many of the source areas.

McElroy, Bill↗
NASA NTRS2001

Methanol Uptake by Low Temperature Aqueous Sulfuric Acid Solutions

The global methanol budget is currently unbalanced, with source terms significantly larger than the sinks terms. To evaluate possible losses of gaseous methanol to sulfate aerosols, the solubility and reactivity of methanol in aqueous sulfuric acid solutions representative of upper tropospheric and lower stratospheric aerosols is under investigation. Methanol will partition into sulfate aerosols according to its Henry's law solubility. Using standard uptake techniques in a Knudsen cell reactor, we have measured the effective Henry's law coefficient, H*, for cold (196 - 220 K) solutions ranging between 45 and 70 wt % H2SO4. We have found that methanol solubility ranges from approx. 10(exp 5) - 10(exp 7) M/atm for UT/LS conditions. Solubility increases with decreasing temperature and with increasing sulfuric acid content. Although methanol is slightly more soluble than are acetone and formaldehyde, current data indicate that uptake by clean aqueous sulfuric acid particles will not be a significant sink for methanol in the UT/LS. These solubility measurements include uptake due to physical solvation and any rapid equilibria which are established in solution. Reaction between primary alcohols and sulfuric acid does occur, leading to the production of alkyl sulfates. Literature values for the rate of this reaction suggest that formation of CH3OSO3H is not significant over our experimental time scale for solutions below 80 wt % H2SO4. To confirm this directly, results obtained using a complementary equilibrium measurement technique will also be presented.

Iraci, L. T.↗
NASA NTRS2000

Effects of Space Environment on Flow and Concentration During Directional Solidification

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we numerically investigate the effects of g-jitter on directional solidification. A background gravity of 1 micro-g has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Benjapiyaporn, C.↗
NASA NTRS2000

Effect of g-jitter on Directional Solidification of a Binary Alloy

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter: the small, rapid fluctuations in gravitational acceleration which may be experienced in an orbiting space vehicle. A background gravity of 1 micro-g has been assumed, and new results for the effects of orientation angle of the periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Santiviriyapanich, P.↗
NASA NTRS2001

Applications of the Method of Space-Time Conservation Element and the Solution Element to Unsteady Chemically Reactive Flows

This document reports the conclusion and findings of our research activities for this grant. The goal of the project is the development and application of the method of Space-Time Conservation Element and Solution Element, or the CE/SE method, to simulate chemically reacting flows. The product of this project will be a high-fidelity, time-accurate flow solver analyzing unsteady flow fields advanced propulsion concepts, including the low-emission turbojet engine combustion and flow fields of the Pulse Detonation Engines (PDE). Based on the documents and computer software of the CE/SE method that we have received from the CE/SE working group at NASA Lewis, we have focused our research effort on addressing outstanding technical issues related to the extension of the CE/SE method for unsteady, chemically reactive flows. In particular, we have made progresses in the following three aspects: (1) Derivation of the governing equations for reacting flows; (2) Numerical treatments of stiff source terms; and (3) Detailed simulations of ZND detonation waves.

Yu, Sheng-Tao↗
NASA NTRS2001

A Numerical Study of Tropical Sea-Air Interactions Using a Cloud Resolving Model Coupled with an Ocean Mixed-Layer Model

Coupling a cloud resolving model (CRM) with an ocean mixed layer (OML) model can provide a powerful tool for better understanding impacts of atmospheric precipitation on sea surface temperature (SST) and salinity. The objective of this study is twofold. First, by using the three dimensional (3-D) CRM-simulated (the Goddard Cumulus Ensemble model, GCE) diabatic source terms, radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the OML model, the respective impact of individual component on upper ocean heat and salt budgets are investigated. Secondly, a two-way air-sea interaction between tropical atmospheric climates (involving atmospheric radiative-convective processes) and upper ocean boundary layer is also examined using a coupled two dimensional (2-D) GCE and OML model. Results presented here, however, only involve the first aspect. Complete results will be presented at the conference.

Shie, Chung-Lin↗
NASA NTRS2001

Io's Magnetospheric Interaction: An MHD Model with Day-Night Asymmetry

In this paper we present the results of all improved three-dimensional MHD model for Io's interaction with Jupiter's magnetosphere. We have included the day-night asymmetry into the spatial distribution of our mass-loading, which allowed us to reproduce several smaller features or the Galileo December 1995 data set. The calculation is performed using our newly modified description of the pick-up processes that accounts for the effects of the corotational electric field existing in the Jovian magnetosphere. This change in the formulation of the source terms for the MHD equations resulted in significant improvements in the comparison with the Galileo measurements. We briefly discuss the limitations of our model and possible future improvements.

Kabin, K.↗