Search NASA⌕ Search

SEARCH · Search NASA

Results for “NUMERICAL FLUX FUNCTION”

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 NTRS2009

Barium Depletion in Hollow Cathode Emitters

The effect of tungsten erosion, transport and redeposition on the operation of dispenser hollow cathodes was investigated in detailed examinations of the discharge cathode inserts from an 8200 hour and a 30,352 hour ion engine wear test. Erosion and subsequent re-deposition of tungsten in the electron emission zone at the downstream end of the insert reduces the porosity of the tungsten matrix, preventing the ow of barium from the interior. This inhibits the interfacial reactions of the barium-calcium-aluminate impregnant with the tungsten in the pores. A numerical model of barium transport in the internal xenon discharge plasma shows that the barium required to reduce the work function in the emission zone can be supplied from upstream through the gas phase. Barium that flows out of the pores of the tungsten insert is rapidly ionized in the xenon discharge and pushed back to the emitter surface by the electric field and drag from the xenon ion flow. This barium ion flux is sufficient to maintain a barium surface coverage at the downstream end greater than 0.6, even if local barium production at that point is inhibited by tungsten deposits. The model also shows that the neutral barium pressure exceeds the equilibrium vapor pressure of the impregnant decomposition reaction over much of the insert length, so the reactions are suppressed. Only a small region upstream of the zone blocked by tungsten deposits is active and supplies the required barium. These results indicate that hollow cathode failure models based on barium depletion rates in vacuum dispenser cathodes are very conservative.

plasma discharges↗
NASA NTRS1988

Hydrogen line and continuum emission in young stellar objects. I - Excitation model

Two mechanisms that populate the n = 2 level of hydrogen after the Lyman continuum is depleted are identified. They are ionization of N I from its excited states, followed by charge-exchange between N II and H I, and Ly-beta line wing absorption. Both processes involve absorption of the sub-Lyman continuum between 11 and 13.6 eV. With population in level n = 2 thus maintained, the strong Brackett line fluxes observed are then produced as a result of Balmer photoionization. The magnitude of the sub-Lyman continuum dictates the fraction of Balmer continuum that will be absorbed. Numerical calculations for four young stellar objects with luminosity ranging from 10 to 10,000 solar are performed, and it is concluded that this two-step process of sub-Lyman continuum absorption followed by Balmer photoionization can account for the great majority of observed Brackett line fluxes. The location and mass of the emitting as as determined from the Brackett line fluxes are reported. The Br-alpha luminosity is calculated as a function of the mass loss rate.

Kwan, John↗
NASA NTRS1995

Bleed Roughness Mechanism

Suction applied through wall bleed is expected to reduce boundary-layer thickness and lead to a profile that can resist separation. However, experiments with rows of suction holes or slots have shown that the benefits of suction are not realized under certain circumstances. Lee, Sloan and Paynter attributed this to a 'bleed roughness', and modeled the phenomena using a modified inner-length scale with convective lag for the eddy viscosity. Test information was required to determine the roughness values at each Mach number and bleed rate. The present work used numerical calculation, subsonic and supersonic experimentation, and physical modeling were used to explore bleed roughness. It was found that the phenomena can exist in the absence of bleed, that it is a function of geometry and Mach number, and that it does not appear in low subsonic flow. The inner-length scale/convective lag approach was found to be inconsistent with recent physical observation. Alternate steady models for the roughness effect disagreed with mass flux trends and the choked flow limit, and suggested that non-steady mechanisms are present. While new supersonic experiments did exhibit low-frequency plenum oscillations, high-frequency unsteadiness associated with the orifice flow scale appears to be the cause of bleed roughness. In fact, boundary-layer calculations with an oscillatory orifice flow produced the same effect in the time-averaged downstream profile for incompressibile flow. A new empirical turbulence-model was proposed which increased the eddy viscosity throughout the entire boundary-layer by a constant which was calibrated using the zero-bleed profile. Preliminary calculations agreed with the experimental data. The model significantly reduced the test information required to determine roughness values for turbulence modeling, and provided prediction capability for supersonic bleed flows. An abstract of a thesis entitled 'Boundary-Layer Bleed Roughness' is included.

Russell, David A.↗
DOE OSTI2024

A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure Based on Data-Driven Shape Function

Here, a new variable Eddington factor (VEF) model is presented for nonlinear problems of thermal radiative transfer (TRT). The VEF model is data-driven and acts on known (a-priori) radiation-diffusion solutions for material temperatures in the TRT problem. A linear auxiliary problem is constructed for the radiative transfer equation (RTE) whose emission source and opacities are evaluated at these known material temperatures. The solution to this RTE approximates the specific intensity distribution in phase-space and time. It is applied as a shape function to define the Eddington tensor for the presented VEF model. The shape function computed via the auxiliary RTE problem will capture some degree of transport effects within the TRT problem. The VEF moment equations closed with this approximate Eddington tensor will thus carry with them these captured transport effects. In this study, the temperature data comes from multigroup P 1 , P 1/3 , and flux-limited diffusion radiative transfer models. The proposed VEF model can be interpreted as a transport-corrected diffusion reduced-order model. Numerical results are presented on the Fleck-Cummings test problem which models a supersonic wavefront of radiation. The VEF model is shown to improve accuracy by 1–2 orders of magnitude compared to the considered radiation-diffusion model solutions to the TRT problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

A parametric heat transfer study for cryogenic ball bearings in SSME HPOTP

A numerical modeling is to examine the effects of coolant convective heat transfer coefficient and frictional heating on the local temperature characteristics of a ball element in Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump (HPOTP) bearing. The present modeling uses a control-volume based, finite-difference method to solve the non-dimensionalized heat conduction equation in spherical coordinate system. The dimensionless temperature is found as a function of Biot number, heat flux ratio between the two race contacts, and location in the ball. The current results show that, for a given cooling capability, the ball temperature generally increases almost linearly with the heat input from the race-contacts. This increase is always very high at one of the two contacts. An increase in heat transfer coefficient generally reduces the ball temperature and alleviates the temperature gradient, except for the regions very close to the race contacts. For a 10-fold increase of heat transfer coefficient, temperature decrease is 35 percent for the average over entire ball, and 10 percent at the inner-race contact. The corresponding change of temperature gradient displays opposing trends between the regions immediately adjacent to the contacts and the remaining portion of the ball. The average temperature gradient in the vicinity of both contacts increases approximately 70 to 100 percent. A higher temperature gradient produces excessive thermal stress locally which may be detrimental to the material integrity. This, however, is the only unfavorable issue for an increase of heat transfer coefficient.

Chyu, Mingking K.↗
NASA NTRS1997

Some Experiences with High Order Discretization Schemes for Compressible Fluid Flow

It is widely recognized that the structure and complexity of compressible fluid flow places high demands on numerical discretization techniques for the fluid flow equations. Fluid flows arising in external aerodynamics often contain both flow field discontinuities and fluid boundary-layers. Both must be accurately resolved to provide useful information to aerodynamic design and analysis engineers. These accuracy requirements motivated the present author to examine a class of finite-volume techniques on arbitrary triangulated domains based on linear or quadratic reconstruction of integral-averaged data followed by upwind flux function evaluation and small time evolution. More recently, we have considered some new upwind techniques which yield compact discretizations while maintaining higher order accuracy. In the mini-symposium talk we will discuss both of these techniques as well as demonstrate the relative merits of each method by computing a number of aerodynamic flows containing shock waves and boundary-layers.

Barth, Timothy J.↗
NASA NTRS1960

Heat Diffusion in Gases, Including Effects of Chemical Reaction

The diffusion of heat through gases is treated where the coefficients of thermal conductivity and diffusivity are functions of temperature. The diffusivity is taken proportional to the integral of thermal conductivity, where the gas is ideal, and is considered constant over the temperature interval in which a chemical reaction occurs. The heat diffusion equation is then solved numerically for a semi-infinite gas medium with constant initial and boundary conditions. These solutions are in a dimensionless form applicable to gases in general, and they are used, along with measured shock velocity and heat flux through a shock reflecting surface, to evaluate the integral of thermal conductivity for air up to 5000 degrees Kelvin. This integral has the properties of a heat flux potential and replaces temperature as the dependent variable for problems of heat diffusion in media with variable coefficients. Examples are given in which the heat flux at the stagnation region of blunt hypersonic bodies is expressed in terms of this potential.

Hansen, C. Frederick↗
NASA NTRSDate not supplied

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗
NASA NTRS1998

A Two-Fluid, MHD Coronal Model

We describe first results from a numerical two-fluid MHD model of the global structure of the solar corona. The model is two-fluid in the sense that it accounts for the collisional energy exchange between protons and electrons. As in our single-fluid model, volumetric heat and momentum sources are required to produce high speed wind from coronal holes, low speed wind above streamers, and mass fluxes similar to the empirical solar wind. By specifying different proton and electron heating functions we obtain a high proton temperature in the coronal hole and a relatively low proton temperature in the streamer (in comparison with the electron temperature). This is consistent with inferences from SOHO/UVCS, and with the Ulysses/SWOOPS proton and electron temperature measurements which we show from the fast latitude scan. The density in the coronal hole between 2 solar radii and 5 solar radii (2RS and 5RS) is similar to the density reported from SPARTAN 201-01 measurements by Fisher and Guhathakurta. The proton mass flux scaled to 1 AU is 2.4 x 10(exp 8)/sq cm s, which is consistent with Ulysses observations. Inside the closed field region, the density is sufficiently high so that the simulation gives equal proton and electron temperatures due to the high collision rate. In open field regions (in the coronal hole and above the streamer) the proton and electron temperatures differ by varying amounts. In the streamer, the temperature and density are similar to those reported empirically by Li et al and the plasma beta is larger than unity everywhere above approx. 1.5 R(sub s), as it is in all other MHD coronal streamer models.

Suess, Steven T.↗
NASA NTRS1999

A Two-Fluid, MHD Coronal Model

We describe first results from a numerical two-fluid MHD model of the global structure of the solar Corona. The model is two-fluid in the sense that it accounts for the collisional energy exchange between protons and electrons. As in our single-fluid model, volumetric heat and Momentum sources are required to produce high speed wind from Corona] holes, low speed wind above streamers, and mass fluxes similar to the empirical solar wind. By specifying different proton and electron heating functions we obtain a high proton temperature in the coronal hole and a relatively low proton temperature above the streamer (in comparison with the electron temperature). This is consistent with inferences from SOHO/UltraViolet Coronagraph Spectrometer instrument (UVCS), and with the Ulysses/Solar Wind Observations Over the Poles of the Sun instrument (SWOOPS) proton and electron temperature measurements which we show from the fast latitude scan. The density in the coronal hole between 2 and 5 solar radii (2 and 5 R(sub S)) is similar to the density reported from SPARTAN 201.-01 measurements by Fisher and Guhathakurta [19941. The proton mass flux scaled to 1 AU is 2.4 x 10(exp 8)/sq cm s, which is consistent with Ulysses observations. Inside the closed field region, the density is sufficiently high so that the simulation gives equal proton and electron temperatures due to the high collision rate. In open field regions (in the coronal hole and above the streamer) the proton and electron temperatures differ by varying amounts. In the streamer the temperature and density are similar to those reported empirically by Li et al. [1998], and the plasma beta is larger than unity everywhere above approx. 1.5 R(sub S), as it is in all other MHD coronal streamer models [e.g., Steinolfson et al., 1982; also G. A. Gary and D. Alexander, Constructing the coronal magnetic field, submitted to Solar Physics, 1998].

Suess, S. T.↗
NASA NTRS1973

Preliminary study of a wall stabilized constructed arc

An iterative, implicit, finite-difference numerical technique is described which is suitable for obtaining solutions to the governing equations for a gas flowing in an axially symmetric constricted-arc heater. The method is shown to provide adequate solutions for three cases of simple pipe flows found in the literature, and for flow in a constricted-arc heater by direct comparison with experimental data. The comparison with arc-heater data includes static pressure, arc voltage, and wall heat flux, all as a function of axial location, and a radial temperature profile at one axial station. The arc-heater data were taken with air as the test gas at a heater inlet pressure approximately 0.40 atm and two flow rates of 2.2 and 4.8 g/sec. The arc currents investigated were between 377 and 584 amperes.

Graves, R. A.↗
NASA NTRS1975

Computational procedure for three-dimensional boundary layers on aircraft and aerospace vehicles

An implicit numerical procedure is developed to solve the full three-dimensional laminar, boundary layer equations over realistic aircraft and spacecraft configurations. This procedure utilizes splined functions in the direction normal to the surface and finite differences in the cross flow direction. Calculation of off-pitch plane attachment and detachment lines can be calculated. The solution procedure permits the incorporation of a wide variety of turbulent models and boundary conditions including; entropy layer procedures, specified surface temperatures or heat flux distributions, surface catalysis, mass addition, and the coupling to both three-dimensional inviscid codes and surface ablation codes. In addition, reasonable computing times are obtained.

Kendall, R. M.↗
NASA NTRS1983

On the stellar rotation-activity connection

The relationship between rotation rates and surface activity in late-type dwarf stars is explored in a survey of recent theoretical and observational studies. Current theoretical models of stellar-magnetic-field production and coronal activity are examined, including linear kinematic dynamo theory, nonlinear dynamos using approximations, and full numerical simulations of the MHD equations; and some typical results are presented graphically. The limitations of the modeling procedures and the constraints imposed by the physics are indicated. The statistical techniques used in establishing correlations between various observational parameters are analyzed critically, and the methods developed for quasar luminosity functions by Avni et al. (1980) are used to evaluate the effects of upper detection bounds, incomplete samples, and missing data for the case of rotation and X-ray flux data.

Rosner, R.↗
NASA NTRS1993

Effects of a parallel electric field and the geomagnetic field in the topside ionosphere on auroral and photoelectron energy distributions

The consequences of electric field acceleration and an inhomogeneous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one-dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogeneous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of O(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function is investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogeneous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Min, Q.-L.↗
NASA NTRS1993

Effects of a Parallel Electric Field and the Geomagnetic Field in the Topside Ionosphere on Auroral and Photoelectron Energy Distributions

The consequences of electric field acceleration and an inhomogencous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one- dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogencous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of 0(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function in investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogencous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Min, Q.-L.↗
NASA NTRS1999

Impact of Assimilating GOES-Derived Land Surface Variables into the PSU/NCAR MM5

Land surface processes are known to have a profound impact on the overlying atmosphere over a wide range of spatial and temporal scales. Many atmospheric numerical models include special parameterizations to improve the specification and partitioning of surface fluxes which are critical to the accurate prediction of warm season boundary layer behavior, organized mesoscale circulations, and convective precipitation. However, the added degrees of freedom resulting from the inclusion of vegetation and soil schemes require the specification of additional surface parameters such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, and the vertical distribution of soil moisture. As satellite data have become more readily available in recent years, many investigations have attempted to use these new measurements to infer missing components of the surface energy budget. Sensitivity studies have shown land-skin temperature (LST) tendencies during the mid-morning hours are strongly sensitive to the surface moisture availability (a function of soil wetness and vegetation) and less sensitive to other parameters such as surface roughness. Based upon results from these and other studies, developed a simple technique that dynamically assimilates Geostationary Operational Environmental Satellite (GOES) derived land-surface products into the surface energy budget of a mesoscale model. The purpose of this paper is to demonstrate that assimilating the GOES satellite data has the potential to improve the representation of land surface characteristics within the model without prior knowledge of the land surface characteristics. The assimilation technique is presented in Section 2 and the numerical experiments are detailed in Section 3. Preliminary results and conclusions are presented in Sections 4 and 5, respectively.

Lapenta, William M.↗
NASA NTRS1993

A simulated Antarctic fast ice ecosystem

A 2D numerical ecosystem model of Antarctic land fast ice is developed to elucidate the primary production with the Antarctic sea ice zone. The physical component employs atmospheric data to simulate congelation ice growth, initial brine entrapment, desalination, and nutrient flux. The biological component is based on the concept of a maximum temperature-dependent algal growth rate which is reduced by limitations imposed from insufficient light or nutrients, as well as suboptimal salinity. Preliminary simulations indicate that, during a bloom, microalgae are able to maintain their vertical position relative to the lower congelation ice margin and are not incorporated into the crystal matrix as the ice sheet thickens. It is inferred that land fast sea ice contains numerous microhabitats that are functionally distinct based upon the unique set of processes that control microalgal growth and accumulation within each.

Arrigo, Kevin R.↗
NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗