Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized gradient approximations”

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 109 records · Page 6

Modeling of a Non-Dilute Alloy Solidification Under Terrestrial and Microgravity Conditions

A general finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe from the melt. 2-D and 3-D models with and without Oberbeck-Boussinesq (O-B) approximation were used. Verification of the 3-D model was based on experimental data. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. Double-diffusive convection in the melt is caused by radial temperature gradients and by rejection at the solid-melt interface of the heavier HgTe-rich solute. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. Models based on complete set of Navier-Stokes equations for compressible media are superior to the O-B approximation in microgravity, because they provide more accurate estimation of convection.

Bune, Andris V.↗

Morphological Evolution of Directional Solidification Interfaces in Microgravity: An Analysis of Model Experiments Performed on the International Space Station

A series of experiments performed using the Pore Formation and Mobility Investigation (PFMI) apparatus within the glovebox facility (GBX) on board the International Space Station (ISS) has provided video images of the morphological evolution of a three-dimensional interface in a diffusion controlled regime. The experimental samples were prepared on ground by filling glass tubes, 1 cm ID and approximately 30 cm in length, with "alloys" of succinonitrile (SCN) and water in an atmosphere of nitrogen at 450 millibar pressure. The compositions of the samples processed and analyzed are 0.25,0.5 and 1.0 wt% water. Experimental processing parameters of temperature gradient and translation speed, as well as camera settings, were remotely monitored and manipulated from the ground Telescience Center (TSC) at the Marshall !3pace Flight Center. During the experiments, the sample was first subjected to a unidirectional melt back, generally at 10 microns per second, with a constant temperature gradient ahead of the melting interface. Following the melt back, the interface was allowed to stabilize before translation is initiated. The temperatures in the sample were monitored by six in situ thermocouples and the position is monitored by an optical linear encoder. For the experiments performed and analyzed, the gradients ranged from 2.5 - 3.3 K/mm and the initial pulling velocities ranged from 0.7 micron per second to 1 micron per second with subsequent transition velocities of up to 100 microns per second. The data provided by the PFMI for analysis includes near-real-time (NRT) video captured on the ground during the experiment runs, ISS Video Tape Recorder (VTR) data dumped from the VTR at the end of the experiment run and recorded on the ground, telemetry data including temperature and position measurements, and limited flight HI-8 tapes in 2 camera views of experiment runs for which tapes have been returned to the investigators from ISS. Because of limited down mass from the ISS, the majority of the initial analysis has been performed using the NRT and VTR video data but will be supplemented with the HI-8 video as it becomes available. hardware and procedures required to prepare samples for correlation to flight samples is described. Using this ground-based system, a series of experiments have been performed for direct comparison with the flight data. The results of these comparisons as well as implications to future microgravity experiments are presented and discussed. A ground-based thin-sample directional solidification system as well as all associated

Strutzenberg, Louise L.↗

Forecasting Coronal Mass Ejections from Magnetograms

We report further results from our ongoing assessment of magnetogram-based measures of active-region nonpotentiality (magnetic shear and twist) and size as predictors of coronal mass ejections (CMEs). From a set of 36 vector magnetograms of bipolar active regions (Falconer, Moore, & Gary 2004, ApJ, submitted), we have found that (1) each of our 5 measures of active-region nonpotentiality has a 75-80% success rate in predicting whether an active region will produce a CME within 2 days after the magnetogram, (2) the correlation of each of these measures with CME production in this window is statistically significant (confidence level greater than 95%), (3) our measure of active-region size has a 65% success rate in predicting CMEs in this window, but (4) the correlation between active-region size and CME productivity is not statistically significant (confidence level approximately 80%). As part of the work under our pending 2003 LWS TR&T proposal, we will double our sample to approximately 70 active regions in order to demonstrate a statistically significant correlation between active-region size and CME productivity, and to determine whether this correlation is a byproduct of any Sun-produced correlation between magnetic twist and size of active regions. Since the 2002 LWS Science Workshop, we have found a measure of active-region nonpotentiality (the length of strong-gradient neutral line) that can be measured from a line-of-sight magnetogram, and we have generalized this measure and the corresponding measure from a vector magnetogram (the length of strong-shear neutral line) so that they can be applied to multiple-bipole complex active regions as well as to bipolar active regions. Preliminary results will be presented from application of these two generalized measures to our sample of bipolar active regions and to a new sample of multiple-bipole active regions.

Falconer, D. A.↗

An Analytic Function of Lunar Surface Temperature for Exospheric Modeling

We present an analytic expression to represent the lunar surface temperature as a function of Sun-state latitude and local time. The approximation represents neither topographical features nor compositional effects and therefore does not change as a function of selenographic latitude and longitude. The function reproduces the surface temperature measured by Diviner to within +/-10 K at 72% of grid points for dayside solar zenith angles of less than 80, and at 98% of grid points for nightside solar zenith angles greater than 100. The analytic function is least accurate at the terminator, where there is a strong gradient in the temperature, and the polar regions. Topographic features have a larger effect on the actual temperature near the terminator than at other solar zenith angles. For exospheric modeling the effects of topography on the thermal model can be approximated by using an effective longitude for determining the temperature. This effective longitude is randomly redistributed with 1 sigma of 4.5deg. The resulting ''roughened'' analytical model well represents the statistical dispersion in the Diviner data and is expected to be generally useful for future models of lunar surface temperature, especially those implemented within exospheric simulations that address questions of volatile transport.

moon↗

Turbulent transport observed just above the Amazon forest

Observations of turbulent heat, moisture, and momentum transport made at two levels, approximately 5 and 10 m above the Amazon rain forest canopy, are presented. Data acquired at 10 Hz included variances and some mixed third moments of vertical velocity, temperature, and humidity. Two features of the data appear to question the displacement height hypothesis: (1) the characteristic dissipation length scale in the near-canopy layer varied between 20 m in stable conditions to approximately 150 m during afternoon convective conditions, generally larger scales than would be expected, and (2) no appreciable differences in dissipation scales was seen at two observed levels. Heat budgets on selected days show that frequent periods with negative heat flux concurrent with continuing positive moisture flux occur in early afternoon, and this is believed to indicate the patchy nature of canopy-atmosphere coupling. Vertical velocity skewness was observed to be negative on three successive days and exhibited a sharp positive gradient. Time series of some of the terms in the turbulence budgets of vertical velocity variance and kinematic heat flux are presented.

Fitzjarrald, David R.↗

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.↗

Near-symmetric instability for general Prandtl number

An analysis of the stability of a simple baroclinic flow to perturbations whose horizontal wave vector lies in and near the vertical plane containing the density gradient has been performed. The Ekman number (E) and 'azimuthal' wavenumber (alpha) were both assumed much less than 1, and expansions about these quantities were performed for the eigenmodes and Richardson number (R). Agreement with the previous analysis of Busse and Chen (1981) that the correction to the critical R was O(alpha) for Prandtl number (P) away from unity was obtained. An expression for the correction to the critical R for arbitrary P, and an approximate expression for P = 1 were obtained as new results. The correction for P = 1 has the same sign as that for P less than 1, and is O(alpha E exp 2/3). This result compares well with the numerical results of Miller and Antar (1986) for small E.

Reynolds, Nathaniel D.↗

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗

Three-dimensional inviscid analysis of radial turbine flow and a limited comparison with experimental data

The three-dimensional inviscid DENTON code is used to analyze flow through a radial-inflow turbine rotor. Experimental data from the rotor are compared with analytical results obtained by using the code. The experimental data available for comparison are the radial distributions of circumferentially averaged values of absolute flow angle and total pressure downstream of the rotor exit. The computed rotor-exit flow angles are generally underturned relative to the experimental values, which reflect the boundary-layer separation at the trailing edge and the development of wakes downstream of the rotor. The experimental rotor is designed for a higher-than-optimum work factor of 1.126 resulting in a nonoptimum positive incidence and causing a region of rapid flow adjustment and large velocity gradients. For this experimental rotor, the computed radial distribution of rotor-exit to turbine-inlet total pressure ratios are underpredicted due to the errors in the finite-difference approximations in the regions of rapid flow adjustment, and due to using the relatively coarser grids in the middle of the blade region where the flow passage is highly three-dimensional. Additional results obtained from the three-dimensional inviscid computation are also presented, but without comparison due to the lack of experimental data. These include quasi-secondary velocity vectors on cross-channel surfaces, velocity components on the meridional and blade-to-blade surfaces, and blade surface loading diagrams. Computed results show the evolution of a passage vortex and large streamline deviations from the computational streamwise grid lines. Experience gained from applying the code to a radial turbine geometry is also discussed.

Choo, Y. K.↗

Three-dimensional inviscid analysis of radial-turbine flow and a limited comparison with experimental data

The three-dimensional inviscid DENTON code is used to analyze flow through a radial-inflow turbine rotor. Experimental data from the rotor are compared with analytical results obtained by using the code. The experimental data available for comparison are the radial distributions of circumferentially averaged values of absolute flow angle and total pressure downstream of the rotor exit. The computed rotor-exit flow angles are generally underturned relative to the experimental values, which reflect the boundary-layer separation at the trailing edge and the development of wakes downstream of the rotor. The experimental rotor is designed for a higher-than-optimum work factor of 1.126 resulting in a nonoptimum positive incidence and causing a region of rapid flow adjustment and large velocity gradients. For this experimental rotor, the computed radial distribution of rotor-exit to turbine-inlet total pressure ratios are underpredicted due to the errors in the finite-difference approximations in the regions of rapid flow adjustment, and due to using the relatively coarser grids in the middle of the blade region where the flow passage is highly three-dimensional. Additional results obtained from the three-dimensional inviscid computation are also presented, but without comparison due to the lack of experimental data. These include quasi-secondary velocity vectors on cross-channel surfaces, velocity components on the meridional and blade-to-blade surfaces, and blade surface loading diagrams. Computed results show the evolution of a passage vortex and large streamline deviations from the computational streamwise grid lines. Experience gained from applying the code to a radial turbine geometry is also discussed.

Choo, Y. K.↗

Fronts propagating with curvature dependent speed: Algorithms based on Hamilton-Jacobi formulations

New numerical algorithms are devised (PSC algorithms) for following fronts propagating with curvature-dependent speed. The speed may be an arbitrary function of curvature, and the front can also be passively advected by an underlying flow. These algorithms approximate the equations of motion, which resemble Hamilton-Jacobi equations with parabolic right-hand-sides, by using techniques from the hyperbolic conservation laws. Non-oscillatory schemes of various orders of accuracy are used to solve the equations, providing methods that accurately capture the formation of sharp gradients and cusps in the moving fronts. The algorithms handle topological merging and breaking naturally, work in any number of space dimensions, and do not require that the moving surface be written as a function. The methods can be used also for more general Hamilton-Jacobi-type problems. The algorithms are demonstrated by computing the solution to a variety of surface motion problems.

Osher, Stanley↗

Experimental Investigation of Effects of Moderate Sideslip on the Flow Fields near a 45 Degree Swept-wing-fuselage Combination at Low Speed

The flow fields near a 45 degree swept-wing-fuselage combination at moderate angles of sideslip (plus-or-minus 8 degrees), as determined experimentally at low speed, are presented as variations with chordwise distance for various spanwise and vertical locations and angles of attack. The results indicated that for positions close to the fuselage (on and near the plane of symmetry) changes in the angle of sideslip caused large changes in the flow-field characteristics and particularly in the local angles of sideslip, which in some cases were nearly double the static angle of sideslip. In general, the effects of changing the angle of sideslip on the flow-field characteristics for all of the outboard underwing locations were qualitatively similar, although conditions at the more inboard and outboard locations were somewhat more severe for lifting conditions than at the one-half semispan location. The chordwise gradients in the flow parameters for the underwing locations were more severe than for the fuselage locations although the effect of changing the angle of sideslip was less severe, in that the incremental changes in the local angles of sideslip were approximately equal to the static angle of sideslip. Flow conditions near the wing tip were found to be critically dependent on vertical location, with the largest sideslip-induced variations occurring at the nearest vertical locations. The results also indicated that for the outboard underwing locations the wing was the predominant factor in disturbing the field of flow for the conditions investigated.

Alford, William J , Jr↗

The role of large-scale eddies in the climate equilibrium. Part 2: Variable static stability

Lorenz's two-level model on a sphere is used to investigate how the results of Part 1 are modified when the interaction of the vertical eddy heat flux and static stability is included. In general, the climate state does not depend very much on whether or not this interaction is included, because the poleward eddy heat transport dominates the eddy forcing of mean temperature and wind fields. However, the climatic sensitivity is significantly affected. Compared to two-level model results with fixed static stability, the poleward eddy heat flux is less sensitive to the meridional temperature gradient and the gradient is more sensitive to the forcing. For example, the logarithmic derivative of the eddy flux with respect to the gradient has a slope that is reduced from approximately 15 on a beta-plane with fixed static stability and approximately 6 on a sphere with fixed static stability, to approximately 3 to 4 in the present model. This last result is more in line with analyses from observations. The present model also has a stronger baroclinic adjustment than that in Part 1, more like that in two-level beta-plane models with fixed static stability, that is, the midlatitude isentropic slope is very insensitive to the forcing, the diabatic heating, and the friction, unless the forcing is very weak.

Zhou, Shuntai↗

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗

Learning and tuning fuzzy logic controllers through reinforcements

A new method for learning and tuning a fuzzy logic controller based on reinforcements from a dynamic system is presented. In particular, our Generalized Approximate Reasoning-based Intelligent Control (GARIC) architecture: (1) learns and tunes a fuzzy logic controller even when only weak reinforcements, such as a binary failure signal, is available; (2) introduces a new conjunction operator in computing the rule strengths of fuzzy control rules; (3) introduces a new localized mean of maximum (LMOM) method in combining the conclusions of several firing control rules; and (4) learns to produce real-valued control actions. Learning is achieved by integrating fuzzy inference into a feedforward network, which can then adaptively improve performance by using gradient descent methods. We extend the AHC algorithm of Barto, Sutton, and Anderson to include the prior control knowledge of human operators. The GARIC architecture is applied to a cart-pole balancing system and has demonstrated significant improvements in terms of the speed of learning and robustness to changes in the dynamic system's parameters over previous schemes for cart-pole balancing.

Berenji, Hamid R.↗

Chemical and physical gradients along the OMC-1 ridge

We present a survey of the distribution of 20 chemical and isotopic molecular species along the central ridge of the Orion molecular cloud from 6' north to 6' south of BN-KL observed with the QUARRY focal plane array on the FCRAO 14 m telescope, which provides an angular resolution of approximately 50" in the 3 mm wavelength region. We use standard tools of multivariate analysis for a systematic investigation of the similarities and differences among the maps of integrated intensities of the 32 lines observed. The maps fall in three broad classes: first, those strongly peaked toward BN-KL; second, those having rather flat distributions along the ridge; and third, those with a clear north-south gradient or contrast. We identify six positions or regions where we calculate relative abundances. Line velocities and line widths indicate that the optically thin lines generally trace the same volume of dense gas, except in the molecular bar, where C18O, C34S, H13CO+, CN, C2H, SO, and C3H2 have velocities characteristic of the bar itself, whereas the emission from other detected species is dominated by the background cloud. The strongest abundance variations in our data are the well-known enhancements seen in HCN, CH3OH, HC3N, and SO toward BN-KL and, less strongly, toward the Orion-South outflow 1'.3S. The principal result of this study is that along the extended quiescent ridge the chemical abundances, within factors of 3-4, exhibit an impressive degree of uniformity. The northern part of the ridge has a chemistry closest to that found in quiescent dense clouds. While temperature and density are similar around the northern radical-ion peak near 3'.5N and in the southern core near 4'.2S, some abundances, in particular, those of the ions HCO+ and N2H+, are significantly lower toward 4'.2S. The areas near 4.'2S and the molecular bar itself around (1'.7E, 2'.4S) stand out with peculiar and similar properties probably caused by stronger UV fields penetrating deeper into the clumpy molecular gas. This leads to higher electron abundances and thereby reduced abundances of the ions, as well as a lack of complex molecules.

NASA Discipline Exobiology↗

Learning and tuning fuzzy logic controllers through reinforcements

This paper presents a new method for learning and tuning a fuzzy logic controller based on reinforcements from a dynamic system. In particular, our generalized approximate reasoning-based intelligent control (GARIC) architecture (1) learns and tunes a fuzzy logic controller even when only weak reinforcement, such as a binary failure signal, is available; (2) introduces a new conjunction operator in computing the rule strengths of fuzzy control rules; (3) introduces a new localized mean of maximum (LMOM) method in combining the conclusions of several firing control rules; and (4) learns to produce real-valued control actions. Learning is achieved by integrating fuzzy inference into a feedforward neural network, which can then adaptively improve performance by using gradient descent methods. We extend the AHC algorithm of Barto et al. (1983) to include the prior control knowledge of human operators. The GARIC architecture is applied to a cart-pole balancing system and demonstrates significant improvements in terms of the speed of learning and robustness to changes in the dynamic system's parameters over previous schemes for cart-pole balancing.

Berenji, Hamid R.↗

Analysis of the Radar Reflectivity of Aircraft Vortex Wakes

Radar has been proposed as a way to track wake vortices to reduce aircraft spacing and tests have revealed radar echoes from aircraft wakes in clear air. The results are always interpreted qualitatively using Tatarski's theory of weak scattering by isotropic atmospheric turbulence. The goal of the present work was to predict the value of the radar cross-section (RCS) using simpler models. This is accomplished in two steps. First, the refractive index is obtained. Since the structure of the aircraft wakes is different from atmospheric turbulence, three simple mechanisms specific to vortex wakes are considered: (1) Radial density gradient in a two-dimensional vortex, (2) three-dimensional fluctuations in the vortex cores, and (3) Adiabatic transport of the atmospheric fluid in a two-dimensional oval surrounding the pair of vortices. The index of refraction is obtained more precisely for the two-dimensional mechanisms than for the three-dimensional ones. In the second step, knowing the index of refraction, a scattering analysis is performed. Tatarski's weak scattering approximation is kept but the usual assumptions of a far-field and a uniform incident wave are dropped. Neither assumption is generally valid for a wake that is coherent across the radar beam. For analytical insight, a simpler approximation that invokes, in addition to weak scattering, the far-field and wide cylindrical beam assumptions, is also developed and compared with the more general analysis. The predicted RCS values for the oval surround the vortices (mechanism C) agree with the experiments of Bilson conducted over a wide range of frequencies. However, the predictions have a cut-off away from normal incidence which is not present in the measurements. Estimates suggest that this is due to turbulence in the baroclinic vorticity generated at the boundary of the oval. The reflectivity of a vortex itself (mechanism A) is comparable to that of the oval (mechanism C) but cuts-off at frequencies lower than those considered in all the experiments to date. The RCS of a vortex happens to peak at the frequency (about 49 MHz) where atmospheric radars (known as ST radars) operate and so the present prediction could be verified in the future. Finally , we suggest that hot engine exhaust could increase RCE by 40 db and reveal vortex circulation, provided its mixing with the surroundings is prevented in the laminarising flow of the vortices.

Shariff, Karim↗