SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2008

Accuracy of Gradient Reconstruction on Grids with High Aspect Ratio

Gradient approximation methods commonly used in unstructured-grid finite-volume schemes intended for solutions of high Reynolds number flow equations are studied comprehensively. The accuracy of gradients within cells and within faces is evaluated systematically for both node-centered and cell-centered formulations. Computational and analytical evaluations are made on a series of high-aspect-ratio grids with different primal elements, including quadrilateral, triangular, and mixed element grids, with and without random perturbations to the mesh. Both rectangular and cylindrical geometries are considered; the latter serves to study the effects of geometric curvature. The study shows that the accuracy of gradient reconstruction on high-aspect-ratio grids is determined by a combination of the grid and the solution. The contributors to the error are identified and approaches to reduce errors are given, including the addition of higher-order terms in the direction of larger mesh spacing. A parameter GAMMA characterizing accuracy on curved high-aspect-ratio grids is discussed and an approximate-mapped-least-square method using a commonly-available distance function is presented; the method provides accurate gradient reconstruction on general grids. The study is intended to be a reference guide accompanying the construction of accurate and efficient methods for high Reynolds number applications

Thomas, James↗

NASA NTRS1988

Analytical approximation of a distorted reflector surface defined by a discrete set of points

Reflector antennas on Earth orbiting spacecrafts generally cannot be described analytically. The reflector surface is subjected to a large temperature fluctuation and gradients, and is thus warped from its true geometrical shape. Aside from distortion by thermal stresses, reflector surfaces are often purposely shaped to minimize phase aberrations and scanning losses. To analyze distorted reflector antennas defined by discrete surface points, a numerical technique must be applied to compute an interpolatory surface passing through a grid of discrete points. In this paper, the distorted reflector surface points are approximated by two analytical components: an undistorted surface component and a surface error component. The undistorted surface component is a best fit paraboloid polynomial for the given set of points and the surface error component is a Fourier series expansion of the deviation of the actual surface points, from the best fit paraboloid. By applying the numerical technique to approximate the surface normals of the distorted reflector surface, the induced surface current can be obtained using physical optics technique. These surface currents are integrated to find the far field radiation pattern.

Acosta, Roberto J.↗

NASA NTRS1979

Convective instability when the temperature gradient and rotation vector are oblique to gravity. I - Fluids without diffusion

A linear stability analysis of fluid layers under uniform rotation (generally oblique to gravity) which possess both vertical and horizontal temperature gradients is made by considering ideal fluids without diffusion within a Boussinesq approximation. This simplified configuration is used to assess the preferred convective modes as a function of latitude on a planet like Jupiter. The tilted rotation vector introduces a preference for roll-like disturbances with north-south orientations, while the horizontal temperature gradient produces a thermal wind shear which favors convective rolls oriented parallel to the flow in an east-west direction. It is found that the horizontal temperature gradient needed to produce a preference for the axisymmetric or east-west rolls increases with an increasing rotation rate and a decreasing latitude. The parameter values for Jupiter are estimated with the use of a simple radiative convective model, indicating a preference for axisymmetric rolls at nearly all latitudes if the convection zone depth is greater than about 200 km below the one atmosphere pressure level and convective roll characteristics which contribute to an equatorial acceleration.

Hathaway, D. H.↗

NASA NTRS1994

An Approximate Solution for the Circumsolar Flow Field of a Sun Moving Through the Local Interstellar Medium (LISM)

The solar system is in relative motion with respect to the ambient interstellar medium. The supersonic solar wind is expected to pass through the termination shock, thus the solar wind plasma eventually has to enter into an asymptotic outflow geometry appropriately adopted to this counterflow situation. Many attempts have been done to simulate the interaction between the solar wind and the LISM numerically. In this paper we generalize a Parker type analytical solution of the counterflow. The idea is to introduce a special kind of compressibility of the solar wind flow. With the assumption that only a transversal component of the density gradient normal to the flow lines exists we are able to calculate a full set of hydrodynamical quantities describing the circumsolar flow field of a Sun moving through the LISM. The equations governing the velocity and density fields lead to analytical solutions which can be taken as good approximations to the more general case of compressible plasma flows.

Ratkiewicz, Romana E.↗

NASA NTRS1994

Low-frequency modes and nonbarotropic effects in pseudo-Newtonian accretion disks

A recently developed formalism is used to reexamine the question of the existence of hydrodynamical modes that pulsate with very low frequencies in the inner regions of accretion disks. The formalism is valid in an exact sense for the adiabatic pulsations of rotating Newtonian fluids that are generally nonbarotropic (such as those with 'nonadiabatic temperature gradients,' for example), and hence its application in the present context represents an improvement over previous analyses that are more approximate. The formalism is applied to thin non-self-gravitating disks, with the gravitational potential of the central source modified in the usual way in order to simulate relativistic effects. In the barotropic limit, the analyses indicate that in many cases nearly Keplerian disks exhibit nonaxisymmetric modes of pulsation that are trapped in the inner disk regions, with pulsation periods much longer than the dynamical timescale. These results are similar to those of earlier calculations that assume disks pulsate without changing the temperature distribution. A method is developed for including lowest order nonbarotropic effects. Previous analyses have been incapable of accurately treating the nonbarotropic regime. The application of the present method to the low-frequency modes reveals that, due to unexpected cancellations among terms, the nonbarotropic correction to the pusation frequency omega is only of order tilde-omega(sub BV exp 2) omega, where tilde-omega(sub BV) is the appropriate dimensionless Brunt-Vaisala frequency. This correction is much smaller than the expected correction of order tilde-omega(sub BV) Omega, where Omega is the rotation angular velocity. The important conclusion drawn from this is that nonbarotropic corrections are generally small and hence that low-frequency modes persist into the nonbarotropic regime. For disk temperatures appropriate to X-ray emission, the adiabatic frequencies of trapped modes are of the same order as the frequencies associated with quasi-periodic oscillations (QPOs) of certain low-mass X-ray binaries and other X-ray sources. It is suggested that further consideration be given to the possiblity that modulation of emission processes by these modes can explain at least some of the QPO phenomena.

Ipser, James R.↗

NASA NTRS1972

Dynamics of stratified liquids in the presence of space charge.

The dispersion characteristics of small signal electrohydrodynamic space charge and polarization waves in continuously and discretely stratified liquids are developed. The liquids, modeled as incompressible, inviscid, and perfectly insulating have distributions of mass density, charge density, permittivity, and electric field. A general set of relations for perturbation field and flow variables defined for a prototype layer are used to describe the dynamics of many superposed planar layers. It is shown that the general relations can be used to represent systems of coupled layers which approximate continuous distributions by a series of step functions. Specific examples of weak gradient and exponential distributions are presented showing that the solution found directly from the distributed theory is approached by the system of coupled layers, if the limit is taken in which the number of layers become large while each layer thickness approaches zero.

Zahn, M.↗

NASA NTRS1995

Effects of pressure gradients on turbulent premixed flames

The influence of a constant acceleration on a turbulent premixed flame is studied by direct numerical simulation. This acceleration induces a mean pressure gradient across the flame brush, leading to a modification of the turbulent flame structure due to differential buoyancy mechanisms between heavy cold fresh and light hot burnt gases. Such a pressure gradient may be encountered in practical applications in ducted flames. A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. A favorable pressure gradient also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. The balance equation for the turbulent flux of the Favre averaged progress variable is also analyzed. The first results show that the fluctuating pressure term, cannot be neglected as generally assumed in models. Simple models assuming that a high mean pressure gradient may only be balanced by the cross-dissipation term seem too approximate. This analysis has to be continued to compare simulation data and closure schemes proposed for the transport equation. The analysis developed by Veynante et al.(1995) has been extended to imposed acceleration and mean pressure gradients. A simple model for the turbulent flux is proposed and validated from simulation data. Then, a modified criterion is derived to delineate between counter-gradient and gradient turbulent diffusion. In fact, counter-gradient diffusion may occur in most practical applications, especially for ducted flames.

Veynante, D.↗

NASA NTRS2000

General Purpose Fortran Program for Discrete-Ordinate-Method Radiative Transfer in Scattering and Emitting Layered Media: An Update of DISORT

This update reports a state-of-the-art discrete ordinate algorithm for monochromatic unpolarized radiative transfer in non-isothermal, vertically inhomogeneous, but horizontally homogeneous media. The physical processes included are Planckian thermal emission, scattering with arbitrary phase function, absorption, and surface bidirectional reflection. The system may be driven by parallel or isotropic diffuse radiation incident at the top boundary, as well as by internal thermal sources and thermal emission from the boundaries. Radiances, fluxes, and mean intensities are returned at user-specified angles and levels. DISORT has enjoyed considerable popularity in the atmospheric science and other communities since its introduction in 1988. Several new DISORT features are described in this update: intensity correction algorithms designed to compensate for the 8-M forward-peak scaling and obtain accurate intensities even in low orders of approximation; a more general surface bidirectional reflection option; and an exponential-linear approximation of the Planck function allowing more accurate solutions in the presence of large temperature gradients. DISORT has been designed to be an exemplar of good scientific software as well as a program of intrinsic utility. An extraordinary effort has been made to make it numerically well-conditioned, error-resistant, and user-friendly, and to take advantage of robust existing software tools. A thorough test suite is provided to verify the program both against published results, and for consistency where there are no published results. This careful attention to software design has been just as important in DISORT's popularity as its powerful algorithmic content.

Tsay, Si-Chee↗

NASA NTRS2008

Validation of Aura Microwave Limb Sounder HCl Measurements

The Earth Observing System (EOS) Microwave Limb Sounder (MLS) aboard the Aura satellite has provided daily global HCl profiles since August 2004. We provide a characterization of the resolution, random and systematic uncertainties, and known issues for the version 2.2 MLS HCl data. The MLS sampling allows for comparisons with many (~1500 to more than 3000) closely matched profiles from the Halogen Occultation Experiment (HALOE) and Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS). These data sets provide HCl latitudinal distributions that are, overall, very similar to those from (coincident) MLS profiles, although there are some discrepancies in the upper stratosphere between the MLS and HALOE gradients. As found in previous work, MLS and ACE HCl profiles agree very well (within approximately 5%, on average), but the MLS HCl abundances are generally larger (by 10-20%) than HALOE HCl. The bias versus HALOE is unlikely to arise mostly from MLS, as a similar systematic bias (of order 15%) is not observed between average MLS and balloon-borne measurements of HCl, obtained over Fort Sumner, New Mexico, in 2004 and 2005. At the largest pressure (147 hPa) for MLS HCl, a high bias (approximately 0.2 ppbv) is apparent in analyses of low to midlatitude data versus in situ aircraft chemical ionization mass spectrometry (CIMS) HCl measurements from the Aura Validation Experiment (AVE) campaigns in 2004, 2005, and 2006; this bias is also observed in comparisons of MLS and aircraftHCl/O3 correlations. Good agreement between MLS and CIMS HCl is obtained at 100 to 68 hPa. The recommended pressure range for MLS HCl is from 100 to 0.15 hPa.

Earth Observing System (EOS)↗

NASA NTRS1977

Introductory Remarks

Advanced transonic configurations, such as the supercritical wing, are inherently more sensitive to Reynolds number than earlier configurations because the pressure recovery gradients imposed on the boundary layer are generally steeper. The results of two-dimensional supercritical airfoil investigations and theoretical calculations show this effect. A technique for approximately simulating full-scale Reynolds number characteristics at present wind-tunnel Reynolds numbers for near-cruise conditions is described. The transition strip, which in the past has been located near the leading edge of the wing, is rearward so that the relative displacement thickness of the boundary layer at the trailing edge of the wing is the same as might be expected on a full-scale configuration with the transition near the leading edge. Two-dimensional wind-tunnel results indicate that the technique provides a very good simulation of airfoil characteristics at full-scale Reynolds number. The need for sorting the effects of Reynolds number and aeroelasticity, which can be done in the NTF, is also addressed.

Whitcomb, R.↗

NASA NTRS1959

Charts and Tables for Estimating the Stability of the Compressible Laminar Boundary Layer with Heat Transfer and Arbitrary Pressure Gradient

The minimum critical Reynolds numbers for the similar solutions of the compressible laminar boundary layer computed by Cohen and Reshotko and also for the Falkner and Skan solutions as recomputed by Smith have been calculated by Lin's rapid approximate method for two-dimensional disturbances. These results enable the stability of the compressible laminar boundary layer with heat transfer and pressure gradient to be easily estimated after the behavior of the boundary layer has been computed by the approximate method of Cohen and Reshotko. The previously reported unusual result (NACA Technical Note 4037) that a highly cooled stagnation point flow is more unstable than a highly cooled flat-plate flow is again encountered. Moreover, this result is found to be part of the more general result that a favorable pressure gradient is destabilizing for very cool walls when the Mach number is less than that for complete stability. The minimum critical Reynolds numbers for these wall temperature ratios are, however, all larger than any value of the laminar-boundary-layer Reynolds number likely to be encountered. For Mach numbers greater than those for which complete stability occurs a favorable pressure gradient is stabilizing, even for very cool walls.

Tetervin, Neal↗

NASA NTRS1982

A nonlocal theory of an electrostatic sinusoidal density drift instability

The stability of space plasmas in which the macroscopic variables vary in the direction x perpendicular to the magnetic field is usually studied by means of the local approximation, in which the fluctuating potential is assumed to be independent of x. To remove this approximation, a nonlocal problem is studied in which the background density, and hence the fluctuating potentials, are periodic in x. From the linear Vlasov/Poisson equations, a set of coupled, linear, homogeneous, algebraic equations relating the Fourier amplitudes of the eigensolution is derived. The equations are solved numerically and a hierarchy of exact eigenmodes characterized by different growth rates and spatial structures is found. At points where the density gradient is locally zero, the mode amplitude is generally several orders of magnitude lower than the peak amplitude. The dependence of these solutions on the parameters which define the background density perturbation is studied and a correspondence between these nonlocal results and those obtained with the local approximation is demonstrated.

Thomsen, M. F.↗

NASA NTRS1976

An iterative procedure for obtaining maximum-likelihood estimates of the parameters for a mixture of normal distributions, 2

The problem of obtaining numerically maximum likelihood estimates of the parameters for a mixture of normal distributions is addressed. In recent literature, a certain successive approximations procedure, based on the likelihood equations, is shown empirically to be effective in numerically approximating such maximum-likelihood estimates; however, the reliability of this procedure was not established theoretically. Here, a general iterative procedure is introduced, of the generalized steepest-ascent (deflected-gradient) type, which is just the procedure known in the literature when the step-size is taken to be 1. With probability 1 as the sample size grows large, it is shown that this procedure converges locally to the strongly consistent maximum-likelihood estimate whenever the step-size lies between 0 and 2. The step-size which yields optimal local convergence rates for large samples is determined in a sense by the separation of the component normal densities and is bounded below by a number between 1 and 2.

Peters, B. C., Jr.↗

NASA NTRS1978

An iterative procedure for obtaining maximum-likelihood estimates of the parameters for a mixture of normal distributions

This paper addresses the problem of obtaining numerically maximum-likelihood estimates of the parameters for a mixture of normal distributions. In recent literature, a certain successive-approximations procedure, based on the likelihood equations, was shown empirically to be effective in numerically approximating such maximum-likelihood estimates; however, the reliability of this procedure was not established theoretically. Here, we introduce a general iterative procedure, of the generalized steepest-ascent (deflected-gradient) type, which is just the procedure known in the literature when the step-size is taken to be 1. We show that, with probability 1 as the sample size grows large, this procedure converges locally to the strongly consistent maximum-likelihood estimate whenever the step-size lies between 0 and 2. We also show that the step-size which yields optimal local convergence rates for large samples is determined in a sense by the 'separation' of the component normal densities and is bounded below by a number between 1 and 2.

Peters, B. C., Jr.↗

NASA NTRS1991

An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations

An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required for the upwind PNS code are approximately equal to an explicit PNS MacCormack's code and existing implicit PNS solvers.

Korte, John J.↗

NASA NTRS1979

Covariance analyses of satellite-derived mesoscale wind fields

Statistical structure functions have been computed independently for nine satellite-derived mesoscale wind fields that were obtained on two different days. Small cumulus clouds were tracked at 5 min intervals, but since these clouds occurred primarily in the warm sectors of midlatitude cyclones the results cannot be considered representative of the circulations within cyclones in general. The field structure varied considerably with time and was especially affected if mesoscale features were observed. The wind fields on the 2 days studied were highly anisotropic with large gradients in structure occurring approximately normal to the mean flow. Structure function calculations for the combined set of satellite winds were used to estimate random error present in the fields. It is concluded for these data that the random error in vector winds derived from cumulus cloud tracking using high-frequency satellite data is less than 1.75 m/s. Spatial correlation functions were also computed for the nine data sets. Normalized correlation functions were considerably different for u and v components and decreased rapidly as data point separation increased for both components. The correlation functions for transverse and longitudinal components decreased less rapidly as data point separation increased.

Maddox, R. A.↗

NASA NTRS1989

New contributions to transit-time damping in multidimensional systems

The existence of two previously unrecognized contributions to transit-time damping in systems of more than one dimension is demonstrated and discussed. It is shown that these contributions cannot be treated by one-dimensional analyses unless it is assumed that the gradient of the field perpendicular to itself always vanishes. Such an assumption is unjustified in general and the new contributions can dominate damping by fast particles in more general situations. Analytic expressions obtained using a Born approximation are found to be in excellent agreement with numerical test-particle calculations of transit-time damping for a variety of field configurations. These configurations include those of a resonance layer and of a spherical wave packet, which approximates a collapsing wave packet in a strongly turbulent plasma. It is found that the fractional power absorption can be strongly enhanced in non-slablike field configurations.

Robinson, P. A.↗

NASA NTRS2006

Small-Scale Gravity Waves in ER-2 MMS/MTP Wind and Temperature Measurements during CRYSTAL-FACE

Lower stratospheric wind and temperature measurements made from NASA's high-altitude ER-2 research aircraft during the CRYSTAL-FACE campaign in July 2002 were analyzed to retrieve information on small scale gravity waves (GWs) at the aircraft's flight level (typically approximately 20 km altitude). For a given flight segment, the S-transform (a Gaussian wavelet transform) was used to search for and identify small horizontal scale GW events, and to estimate their apparent horizontal wavelengths. The horizontal propagation directions of the events were determined using the Stokes parameter method combined with the cross S-transform analysis. The vertical temperature gradient was used to determine the vertical wavelengths of the events. GW momentum fluxes were calculated from the cross S-transform. Other wave parameters such as intrinsic frequencies were calculated using the GW dispersion relation. More than 100GW events were identified. They were generally high frequency waves with vertical wavelength of approximately 5 km and horizontal wavelength generally shorter than 20 km. Their intrinsic propagation directions were predominantly toward the east, whereas their ground-based propagation directions were primarily toward the west. Among the events, approximately 20% of them had very short horizontal wavelength, very high intrinsic frequency, and relatively small momentum fluxes, and thus they were likely trapped in the lower stratosphere. Using the estimated GW parameters and the background winds and stabilities from the NCAR/NCEP reanalysis data, we were able to trace the sources of the events using a simple reverse ray-tracing. More than 70% of the events were traced back to convective sources in the troposphere, and the sources were generally located upstream of the locations of the events observed at the aircraft level. Finally, a probability density function of the reversible cooling rate due to GWs was obtained in this study, which may be useful for cirrus cloud models.

Wang, L.↗