SEARCH · Search NASA

Results for “rational”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Numerical mixing calculations of confined reacting jet flows in a cylindrical duct

The results reported in this paper describe some of the main flow characteristics and NOx production results which develop in the mixing process in a constant cross-sectional cylindrical duct. A 3-dimensional numerical model has been used to predict the mixing flow field and NOx characteristics in a mixing section of an RQL combustor. Eighteen configurations have been analyzed in a circular geometry in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying three parameter: (1) jet-to-mainstream momentum-flux ration (J), (2) orifice shape or orifice aspect ratio, and (3) slot slant angle. The results indicate that the mixing flow field and NOx production significantly vary with the value of the jet penetration and subsequently, slanting elongated slots generally improve the NOx production at high J conditions. Round orifices produce low NOx at low J due to the strong jet penetration. The NOx production trends do not correlate with the mixing non-uniformity parameters described herein.

Oechsle, Victor L.↗

NASA NTRS1995

Seamless data-range change using punctured convolutional codes for time-varying signal-to-noise ratios

In a time-varying signal-to-noise ration (SNR) environment, symbol rate is often changed to maximize data return. However, the symbol-rate change has some undesirable effects, such as changing the transmission bandwidth and perhaps causing the receiver symbol loop to lose lock temporarily, thus losing some data. In this article, we are proposing an alternate way of varying the data rate without changing the symbol rate and, therefore, the transmission bandwidth. The data rate change is achieved in a seamless fashion by puncturing the convolutionally encoded symbol stream to adapt to the changing SNR environment. We have also derived an exact expression to enumerate the number of distinct puncturing patterns. To demonstrate this seamless rate change capability, we searched for good puncturing patterns for the Galileo (14,1/4) convolutional code and changed the data rates by using the punctured codes to match the Galileo SNR profile of November 9, 1997. We show that this scheme reduces the symbol-rate changes from nine to two and provides a comparable data return in a day and a higher symbol SNR during most of the day.

Feria, Y.↗

NASA NTRS1994

A dynamic localization model with stochastic backscatter

The modeling of subgrid scales in large-eddy simulation (LES) has been rationalized by the introduction of the dynamic localization procedure. This method allows one to compute rather than prescribe the unknown coefficients in the subgrid-scale model. Formally, the LES equations are supposed to be obtained by applying to the Navier-Stokes equations a 'grid filter' operation. Though the subgrid stress itself is unknown, an identity between subgrid stresses generated by different filters has been derived. Although preliminary tests of the Dynamic Localization Model (DLM) with k-equation have been satisfactory, the use of a negative eddy viscosity to describe backscatter is probably a crude representation of the physics of reverse transfer of energy. Indeed, the model is fully deterministic. Knowing the filtered velocity field and the subgrid-scale energy, the subgrid stress is automatically determined. We know that the LES equations cannot be fully deterministic since the small scales are not resolved. This stems from an important distinction between equilibrium hydrodynamics and turbulence. In equilibrium hydrodynamics, the molecular motions are also not resolved. However, there is a clear separation of scale between these unresolved motions and the relevant hydrodynamic scales. The result of molecular motions can then be separated into an average effect (the molecular viscosity) and some fluctuations. Due to the large number of molecules present in a box with size of the order of the hydrodynamic scale, the ratio between fluctuations and the average effect should be very small (as a result of the 'law of large numbers'). For that reason, the hydrodynamic balance equations are usually purely deterministic. In turbulence, however, there is no clear separation of scale between small and large eddies. In that case, the fluctuations around a deterministic eddy viscosity term could be significant. An eddy noise would then appear through a stochastic term in the subgrid-scale model and could be the source of backscatter.

Carati, Daniele↗

NASA NTRS1995

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗

NASA NTRS1993

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

NASA NTRS1993

TIGER: A user-friendly interactive grid generation system for complicated turbomachinery and axis-symmetric configurations

The issue of time efficiency in grid generation is addressed by developing a user friendly graphical interface for interactive/automatic construction of structured grids around complex turbomachinery/axis-symmetric configurations. The accuracy of geometry modeling and its fidelity is accomplished by adapting the nonuniform rational b-spline (NURBS) representation. A customized interactive grid generation code, TIGER, has been developed to facilitate the grid generation process for complicated internal, external, and internal-external turbomachinery fields simulations. The FORMS Library is utilized to build user-friendly graphical interface. The algorithm allows a user to redistribute grid points interactively on curves/surfaces using NURBS formulation with accurate geometric definition. TIGER's features include multiblock, multiduct/shroud, multiblade row, uneven blade count, and patched/overlapping block interfaces. It has been applied to generate grids for various complicated turbomachinery geometries, as well as rocket and missile configurations.

Shih, Ming H.↗

NASA NTRS1995

On the breakup of viscous liquid threads

A one-dimensional model evolution equation is used to describe the nonlinear dynamics that can lead to the breakup of a cylindrical thread of Newtonian fluid when capillary forces drive the motion. The model is derived from the Stokes equations by use of rational asymptotic expansions and under a slender jet approximation. The equations are solved numerically and the jet radius is found to vanish after a finite time yielding breakup. The slender jet approximation is valid throughout the evolution leading to pinching. The model admits self-similar pinching solutions which yield symmetric shapes at breakup. These solutions are shown to be the ones selected by the initial boundary value problem, for general initial conditions. Further more, the terminal state of the model equation is shown to be identical to that predicted by a theory which looks for singular pinching solutions directly from the Stokes equations without invoking the slender jet approximation throughout the evolution. It is shown quantitatively, therefore, that the one-dimensional model gives a consistent terminal state with the jet shape being locally symmetric at breakup. The asymptotic expansion scheme is also extended to include unsteady and inerticial forces in the momentum equations to derive an evolution system modelling the breakup of Navier-Stokes jets. The model is employed in extensive simulations to compute breakup times for different initial conditions; satellite drop formation is also supported by the model and the dependence of satellite drop volumes on initial conditions is studied.

Papageorgiou, Demetrios T.↗

NASA NTRS1995

The National Grid Project: A system overview

The National Grid Project (NGP) is a comprehensive numerical grid generation software system that is being developed at the National Science Foundation (NSF) Engineering Research Center (ERC) for Computational Field Simulation (CFS) at Mississippi State University (MSU). NGP is supported by a coalition of U.S. industries and federal laboratories. The objective of the NGP is to significantly decrease the amount of time it takes to generate a numerical grid for complex geometries and to increase the quality of these grids to enable computational field simulations for applications in industry. A geometric configuration can be discretized into grids (or meshes) that have two fundamental forms: structured and unstructured. Structured grids are formed by intersecting curvilinear coordinate lines and are composed of quadrilateral (2D) and hexahedral (3D) logically rectangular cells. The connectivity of a structured grid provides for trivial identification of neighboring points by incrementing coordinate indices. Unstructured grids are composed of cells of any shape (commonly triangles, quadrilaterals, tetrahedra and hexahedra), but do not have trivial identification of neighbors by incrementing an index. For unstructured grids, a set of points and an associated connectivity table is generated to define unstructured cell shapes and neighboring points. Hybrid grids are a combination of structured grids and unstructured grids. Chimera (overset) grids are intersecting or overlapping structured grids. The NGP system currently provides a user interface that integrates both 2D and 3D structured and unstructured grid generation, a solid modeling topology data management system, an internal Computer Aided Design (CAD) system based on Non-Uniform Rational B-Splines (NURBS), a journaling language, and a grid/solution visualization system.

Gaither, Adam↗

NASA NTRS1995

Generating grids directly on CAD database surfaces using a parametric evaluator approach

A very important, but often overlooked step in grid generation is acquiring a suitable geometry definition of the vehicle to be analyzed. In the past, geometry was usually obtained by generating a number of cross-sections of each component. A number of recent efforts have focussed on non-uniform rational B-spline surfaces (NURBS) to provide as single type of analytic surface to deal with inside the grid generator. This approach has required the development of tools to read other types of surfaces and convert them, either exactly or by approximation, into a NURBS surface. This paper describes a more generic parametric evaluator approach, which does not rely on a particular surface type internal to the grid generation system and is less restrictive in the number of surface types that can be represented exactly. This approach has been implemented in the McDonnell Douglas grid generation system, MACGS, and offers direct access to all types of surfaces from a Unigraphics part file.

Gatzhe, Timothy D.↗

NASA NTRS1995

Grid generation on trimmed Bezier and NURBS quilted surfaces

This paper presents some recently added capabilities to RAGGS, Rockwell Automated Grid Generation System. Included are the trimmed surface handling and display capability and structures and unstructured grid generation on trimmed Bezier and NURBS (non-uniform rational B-spline surfaces) quilted surfaces. Samples are given to demonstrate the new capabilities.

Woan, Chung-Jin↗

NASA NTRS1995

NASA-IGES Translator and Viewer

NASA-IGES Translator (NIGEStranslator) is a batch program that translates a general IGES (Initial Graphics Exchange Specification) file to a NASA-IGES-Nurbs-Only (NINO) file. IGES is the most popular geometry exchange standard among Computer Aided Geometric Design (CAD) systems. NINO format is a subset of IGES, implementing the simple and yet the most popular NURBS (Non-Uniform Rational B-Splines) representation. NIGEStranslator converts a complex IGES file to the simpler NINO file to simplify the tasks of CFD grid generation for models in CAD format. The NASA-IGES Viewer (NIGESview) is an Open-Inventor-based, highly interactive viewer/ editor for NINO files. Geometry in the IGES files can be viewed, copied, transformed, deleted, and inquired. Users can use NIGEStranslator to translate IGES files from CAD systems to NINO files. The geometry then can be examined with NIGESview. Extraneous geometries can be interactively removed, and the cleaned model can be written to an IGES file, ready to be used in grid generation.

Chou, Jin J.↗

NASA NTRS1995

Proton affinity determinations and proton-bound dimer structure indications in C2 to C15, (alpha),(omega)-alkyldiamines

The 'kinetic method' was used to determine the proton affinity (PA) of a,coalkyldiamines from collision induced dissociation (CID) studies of protonated heterodimers. These PA values were consistently lower than those reported in the proton affinity scale. The apparent discrepancy was rationalized in terms of differences in the conformation of the protonated diamine monomers. The minimum energy species, formed by equilibrium proton transfer processes, have a cyclic conformation and the ion charge is shared by both amino-groups which are bridged by the proton. On the other hand, the species formed through dissociation of protonated dimers have a linear structure and the charge is localized on one of the amino-groups. Thus, the difference in the PA values obtained by both methods is a measure of the additional stability acquired by the protonated diamines through cyclization and charge delocalization. The major collision dissociation pathway of the protonated diamine monomers involved elimination of an ammonia moiety. Other reactions observed included loss of the second amino-group and several other bond cleavages. CID of the protonated dimers involved primarily formation of a protonated monomer through cleavage of the weaker hydrogen bond and subsequently loss of ammonia at higher collision energies. As observed from the CID studies, doubly charged ions were also formed from the diamines under conditions of the electrospray ionization.

Karpas, Z.↗

NASA NTRS1994

Factors affecting the estimate of primary production from space

Remote sensing of primary production in the euphotic zone has been based mostly on visible-band and water-leaving radiance measured with the coastal zone color scanner. There are some robust, simple relationships for calculating integral production based on surface measurements, but they also require knowledge for photoadaptive parameters such as maximum photosynthesis which currently cannot be obtained from spave. A 17,000-station data set is used to show that space-based estimates of maximum photosynthesis could improve predictions of psi, the water column light utiliztion index, which is an important term in many primary productivity models. Temperature is also examined as a factor for predicting hydrographic structure and primary production. A simple model is used to relate temperature and maximum photosynthesis; the model incorporates (1) the positive relationship between maximum photosynthesis and temperature and (2) the strongly negative relationship between temperature and nitrate in the ocean (which directly affects maximum growth rates via nitrogen limitation). Since these two factors relate to carbon and nitrogen, 'balanced carbon/nitrogen assimilation' was calculated using the Redfield ratio, It is expected that the relationship between maximum balanced carbon assimilation versus temperature is concave-down, with the peak dependent on nitrate uptake kinetics, temperature-nitrate relationships,a nd the carbon chlorophyll ration. These predictions were compared with the sea truth data. The minimum turnover time for nitrate was also calculated using this approach. Lastly, sea surface temperature gradients were used to predict the slope of isotherms (a proxy for the slope of isopycnals in many waters). Sea truth data show that at size scales of several hundred kilometers, surface temperature gradients can provide information on the slope of isotherms in the top 200 m of the water column. This is directly relevant to the supply of nutrients into the surface mixed layer, which is useful for predicting integral biomass and primary production.

Balch, W. M.↗

NASA NTRS1994

Heterogeneous conversion of N2O5 to HNO3 in the post-Mount Pinatubo eruption stratosphere

Simultaneous stratospheric volume mixing ration (VMR) profiles of dinitrogen pentoxide (N2O5) and nitric acid (HNO3) at sunrise between 25 deg N and 15 deg S latitude and profiles of HNO3 at sunset between 42 deg S and 53 deg S latitude have been derived from 0.01/cm resolution infrared solar occultation spectra recorded 9.5 months after the massive eruption of the Mount Pinatubo volcano in the Philippine Islands. The measurements were obtained by the atmospheric trace molecule spectroscopy (ATMOS) Fourier transform spectrometer during the ATLAS 1 shuttle mission (March 24 to April 2, 1992). The measured HNO3 VMRs are higher at all altitudes and latitudes than corresponding values measured by the limb infrared monitor of the stratosphere (LIMS) instrument during the same season in 1979, when the aerosol loading was near background levels. The largest relative increase in the HNO3 VMR occurred near the equator at 30-km altitude, where the ATMOS/ATLAS 1 values are about a factor of 2 higher than the LIMS measurements. Two-dimensional model calculations show that the increase in HNO3 and the ATMOS/ATLAS 1 measurement of a steep decrease in the N2O5 VMR below 30 km can be explained by the enhanced conversion of N2O5 to HNO3 on the surfaces of the Mount Pinatubo sulfate aerosols. Our profile results demonstrate the global impact of the N2O5 + H2O yields 2HNO3 heterogeneous reaction in altering the partitioning of stratospheric odd nitrogen after a major volcanic eruption.

Rinsland, C. P.↗

NASA NTRS1994

Aubrite basalt vitrophyres: High sulfur silicate melts and a snapshot of aubrite formation

Two aubrite basalt vitrophyre clasts have been found within AMNH thin sections from the Parsa EH3 chondrite and the Khor Temiki aubrite. Polished sections of the Parsa Aubrite Inclusion (PAI) and the Khor Temiki Inclusion (KTI) were studied by optical, electron probe microanalysis (EPMA), and scanning electron microscopy (SEM) techniques with broad-beam and low absorbed EPMA currents used to minimize glass volatile loss. Some data have previously been reported for PAI and KTI may possibly correlate to a previously reported inclusion in Khor Tiimiki. In polished sections, PAI and KTI are approximately equal 4 mm in diameter and contain a large volume of glass. The clasts have similar textural characteristics and are akin to lunar vitrophyre textures. The glasses have high alkali rhyodacitic compositions Al-though PAI is peraluminous, KTI is significantly peralkaline. Additionally, the glasses have elevated sulfur concentrations that are extremely high by geochemical standards. SEM examination for beam overlap of microscopic CaS, FeS, and (Mg, Mn, Fe) S inclusions showed no such contamination. Furthermore, homogeneity of glass S content and low FeO contents help rule out contamination. Materials research data show that under reducing conditions alumino-silicate melts can dissolve up to several weight percent sulfur in the absence of Fe. The high S and alkali contents, the lack of associated high shock features, and the rationalized phase equilibria suggest that PAI and KTI are igneous melting products of an E-chondrite-like source material. Although large-scale impact melting cannot totally be ruled out, the above observations eliminate the possibility of in-situ shock melting.

Fogel, R. A.↗

NASA NTRS1993

A highly accurate Voigt function algorithm

A complex Voigt lineshape algorithm is presented whose maximum relative error over the complex plane is less than 1 x 10(exp -8). The algorithm consists of series, rational approximations and Gauss-Hermite integrations which makes it suitable as a general purpose software module for a wide variety of uses, including a Voigt function standard.

Shippony, Z.↗

NASA NTRS1994

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

NASA NTRS1995

Development of improved space sampling strategies for ocean chemical properties: Total carbon dioxide and dissolved nitrate

Large-scale ocean observing programs such as the Joint Global Ocean Flux Study (JGOFS) and the World Ocean Circulation Experiment (WOCE) today, must face the problem of designing an adequate sampling strategy. For ocean chemical variables, the goals and observing technologies are quite different from ocean physical variables (temperature, salinity, pressure). We have recently acquired data on the ocean CO2 properties on WOCE cruises P16c and P17c that are sufficiently dense to test for sampling redundancy. We use linear and quadratic interpolation methods on the sampled field to investigate what is the minimum number of samples required to define the deep ocean total inorganic carbon (TCO2) field within the limits of experimental accuracy (+/- 4 micromol/kg). Within the limits of current measurements, these lines were oversampled in the deep ocean. Should the precision of the measurement be improved, then a denser sampling pattern may be desirable in the future. This approach rationalizes the efficient use of resources for field work and for estimating gridded (TCO2) fields needed to constrain geochemical models.

Goyet, Catherine↗