Search NASA⌕ Search

SEARCH · Search NASA

Results for “Taylor approximation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Approach for Input Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

An implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for quasi 3-D Euler CFD code is presented. Given uncertainties in statistically independent, random, normally distributed input variables, first- and second-order statistical moment procedures are performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, these moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗
NASA NTRS1976

Polynominal Interpolation Methods for Viscous Flow Calculations

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and by Hermitian (Taylor series) finite-difference discretization. The similarities and special features of these different developments are discussed. The governing systems apply for both uniform and variable meshes. Hybrid schemes resulting from two different polynomial approximations for the first and second derivatives lead to a nonuniform mesh extension of the so-called compact or Pad? difference technique (Hermite 4). A variety of fourth-order methods are described and the Hermitian approach is extended to sixth-order (Hermite 6). The appropriate spline boundary conditions are derived for all procedures. For central finite differences, this leads to a two-point, second-order accurate generalization of the commonly used three-point end-difference formula. Solutions with several spline and Hermite procedures are presented for the boundary layer equations, with and without mass transfer, and for the incompressible viscous flow in a driven cavity. Divergence and nondivergence equations are considered for the cavity. Among the fourth-order techniques, it is shown that spline 4 has the smallest truncation error. The spline 4 procedure generally requires one-quarter the number of mesh points in a given coordinate direction as a central finite-difference calculation of equal accuracy. The Hermite 6 procedure leads to remarkably accurate boundary layer solutions.

Rubin, S. G.↗
NASA NTRS1993

Solution of the linear shallow water equations by the fourth-order leapfrog scheme

Numerical schemes of the first and second order of approximation introduce numerical distortion when the wave propagation over a long distance is investigated. To alleviate this problem, the fourth-order leapfrog scheme is constructed. The standard leapfrog method is based on the truncated Taylor series expansion which depicts an error proportional to the second-order terms. In the proposed method the numerical solution is corrected for these terms. The space and time corrections work well in diminishing numerical dispersion and dissipation.

Kowalik, Z.↗
NASA NTRS1978

Formation of a rotational accretion column

We consider the problem of the flow produced by the accretion of gas from a rotating cloud onto a point gravitational source. It is shown that, for a cloud initially in uniform rotation, the effect of accretion is to produce a Taylor column in the cloud parallel to the rotation axis. Only gas within the column is accreted; the gas outside the column undergoes oscillatory flow, which is probably subject to shear instability. The column forms in approximately one-half rotation period from the time at which accretion begins. Its initial radius is determined by the accretion rate and the angular velocity of the cloud. Application of the results to the formation of Jupiter by accretion is discussed.

Cassen, P.↗
NASA NTRS1987

Full CI benchmark calculations on CH3

Full CI calculations have been performed on the CH3 radical. The full CI results are compared to those obtained using CASSCF/multireference CI and coupled-pair functional methods, both at the equilibrium CH distance and at geometries with the three CH bonds extended. In general, the performance of the approximate methods is similar to that observed in calculations on other molecules in which one or two bonds were stretched.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1977

Coarse-grained basalt 71597 - A product of partial olivine accumulation

Rock 71597, a coarse-grained plagioclase-poikilitic basalt that contains 15.8% MgO and 19.3% modal olivine, is described. Olivine and ilmenite crystals occur in two size generations, and large subequant olivine crystals, small 'matrix' olivines, and large ilmenite crystals are characterized. Since large olivine and ilmenite crystals with skeletal habits suggestive of crystallization under relatively rapid cooling are found in a rock whose texture is consonant with slower cooling, it is suggested that the large skeletal crystals formed toward the more rapidly cooled margin of a flow and accumulated by gravitational crystal settling toward the flow's interior. In support of this hypothesis, calculations indicate that the observed olivine compositions are approximately 10 mol% too Fe-rich to have crystallized from a melt of composition corresponding to that of 71597.

Warner, R. D.↗
NASA NTRS2001

Approach for Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

This paper presents an implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for a quasi 1-D Euler CFD (computational fluid dynamics) code. Given uncertainties in statistically independent, random, normally distributed input variables, a first- and second-order statistical moment matching procedure is performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, the moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗
NASA NTRS1994

On the Lighthill relationship and sound generation from isotropic turbulence

In 1952, Lighthill developed a theory for determining the sound generated by a turbulent motion of a fluid. With some statistical assumptions, Proudman applied this theory to estimate the acoustic power of isotropic turbulence. Recently, Lighthill established a simple relationship that relates the fourth-order retarded time and space covariance of his stress tensor to the corresponding second-order covariance and the turbulent flatness factor, without making statistical assumptions for a homogeneous turbulence. Lilley revisited Proudman's work and applied the Lighthill relationship to evaluate directly the radiated acoustic power from isotropic turbulence. After choosing the time separation dependence in the two-point velocity time and space covariance based on the insights gained from direct numerical simulations, Lilley concluded that the Proudman constant is determined by the turbulent flatness factor and the second-order spatial velocity covariance. In order to estimate the Proudman constant at high Reynolds numbers, we analyzed a unique data set of measurements in a large wind tunnel and atmospheric surface layer that covers a range of the Taylor microscale based on Reynolds numbers 2.0 x 10(exp 3) less than or equal to R(sub lambda) less than or equal to 12.7 x 10(exp 3). Our measurements demonstrate that the Lighthill relationship is a good approximation, providing additional support to Lilley's approach. The flatness factor is found between 2.7 - 3.3 and the second order spatial velocity covariance is obtained. Based on these experimental data, the Proudman constant is estimated to be 0.68 - 3.68.

Zhou, YE↗
NASA NTRS1991

Helium abundances on the moon: Assumptions and estimates

Nuclear energy is a highly desirable source of energy, and He-3 is the most prized of the fusion reactants. As the Wisconsin Group has emphasized, He-3 may be the only true economic ore on the Moon. The lack of a shielding atmosphere on the Moon permits solar-wind alpha particles to impinge upon the lunar regolith and become implanted into the various solid components. In particular, large quantities of helium (5 to 50 ppm) are presented. The measured parameter of I(sub s)/FeO, a direct indicator of maturity and exposure age, can be used as a first approximation to predict the abundances of many solar-wind components in the soils. However, because ilmenite has a much higher retentivity for helium than the other phases, the TiO2 contents of the soils are better indicators of helium contents (Taylor, Space 90). High-Ti mare bassalt regions, such as at the Apollo 17 locale, appear to be the best areas for He mining (15 to 50 ppm He(sub T)), versus 3 to 9 ppm in the Highlands. However, the relationships between I(sub s)/FeO, TiO2 and He-3 contents are complicated - e.g., many of the most He-rich soils are immature to submature. The amount of He-3 in the regolith of the moon is estimated at 220,000 tons in the outer 2 m of the Maria.

Taylor, Lawrence A.↗
NASA NTRS1995

An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case

Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.

Goodrich, John W.↗
NASA NTRSDate not supplied

Higher-Order Approximations for Stabilizing Zero-Energy Modes in Peridynamics Crystal Plasticity Models with Large Horizon Interactions

The non-ordinary state-based peridynamics theory combines non-local dynamic techniques with a desirable correspondence material principle, allowing for the use of continuum mechanics constitutive models. Such an approach presents a unique capability for solving problems involving discontinuities (e.g., strain localization, fracture, and fragmentation). However, the correspondence-based peridynamics models often suffer from zero-energy mode instabilities in numerical implementation, primarily due to the approximations of the non-local deformation gradient tensor. This paper focuses on a computational scheme for eliminating the zero-energy mode oscillations using a choice of influence functions that improve the truncation error in a higher-order Taylor series expansion of the deformation gradient. The novelty here is a tensor-based derivation of the linear constraint equations, which can be used to systematically identify the particle interaction weight functions for various user-specified horizon radii. In this paper, the proposed higher-order stabilization scheme is demonstrated for multi-dimensional examples involving polycrystalline and composite microstructures, along with comparisons against conventional finite element methods. The proposed stabilization scheme is shown to be highly effective in suppressing the spurious zero-energy mode oscillations in all numerical examples while enabling efficient simulations of strain localizations across material interfaces.

Non-Ordinary State-Based Peridynamics↗
NASA NTRS1991

An efficient method for estimating neighboring steady-state numerical solutions to the Euler equations

The paper concentrates on a linear approximation method for predicting the changes occurring in steady-state numerical solutions of the Euler equations as a consequence of small changes in the independent variables which control the problem. The importance of proper boundary-condition treatment and other issues concerning the problem are covered along with the importance of proper algorithm selection for a fully supersonic inviscid flow. The method is applied to a subsonic nozzle involving variation of the pressure on the outflow boundary and to a supersonic inlet involving variation of the inflow Mach number. In the subsonic test case, the comparisons between the predicted and conventional numerical solutions are shown to be good, while in the supersonic test case, the agreement between the approximation method and conventional numerical solution starts out well but rapidly degenerates at some point in the flowfield as the perturbation of the boundary conditions is increased.

Taylor, Arthur C., III↗
NASA NTRS2014

Big Impacts and Transient Oceans on Titan

We have studied the thermal consequences of very big impacts on Titan [1]. Titan's thick atmosphere and volatile-rich surface cause it to respond to big impacts in a somewhat Earth-like manner. Here we construct a simple globally-averaged model that tracks the flow of energy through the environment in the weeks, years, and millenia after a big comet strikes Titan. The model Titan is endowed with 1.4 bars of N2 and 0.07 bars of CH4, methane lakes, a water ice crust, and enough methane underground to saturate the regolith to the surface. We assume that half of the impact energy is immediately available to the atmosphere and surface while the other half is buried at the site of the crater and is unavailable on time scales of interest. The atmosphere and surface are treated as isothermal. We make the simplifying assumptions that the crust is everywhere as methane saturated as it was at the Huygens landing site, that the concentration of methane in the regolith is the same as it is at the surface, and that the crust is made of water ice. Heat flow into and out of the crust is approximated by step-functions. If the impact is great enough, ice melts. The meltwater oceans cool to the atmosphere conductively through an ice lid while at the base melting their way into the interior, driven down in part through Rayleigh-Taylor instabilities between the dense water and the warm ice. Topography, CO2, and hydrocarbons other than methane are ignored. Methane and ethane clathrate hydrates are discussed quantitatively but not fully incorporated into the model.

methane↗
NASA NTRS1983

Some spatial and temporal relationships between lightning and storm structure and evolution

A number of systems were used to study relationships between lightning and other storm parameters. Instrumentation includes a system for locating cloud to ground lightning strike points, an acoustic lightning mapper and very high frequencies (VHF) (30 to 80 MHz) lightning mapper to study the structure and location of lightning, and S band radar to provide other storm data. Several characteristics of the reconstructed lightning have been relatively constant: (1) lightning structure is usually predominately horizontal; (2) lightning tends to occur preferrentially in certain altitude ranges that are only a few kilometers thick and above the O C isotherm; (3) it tends to occur near, but not necessarily inside, cores of high reflectivity (approximately or 45 dBz), and it often connects with regions of relatively weak reflectivity downstream from stronger reflectivity and updraft cores. Other relationships are also discussed.

Macgorman, D. R.↗
NASA NTRS1984

Helicopter rotor blade design for minimum vibration

The importance of blade design parameters in rotor vibratory response and the design of a minimum vibration blade based upon this understanding are examined. Various design approaches are examined for a 4 bladed articulated rotor operating at a high speed flight condition. Blade modal shaping, frequency placement, structural and aerodynamic coupling, and intermodal cancellation are investigated to systematically identify and evaluate blade design parameters that influence blade airloads, blade modal response, hub loads, and fuselage vibration. The relative contributions of the various components of blade force excitation and response to the vibratory hub loads transmitted to the fuselage are determined in order to isolate primary candidates for vibration alleviation. A blade design is achieved which reduces the predicted fuselage vibration from the baseline blade by approximately one half. Blade designs are developed that offer significant reductions in vibration (and fatigue stresses) without resorting to special vibration alleviation devices, radical blade geometries, or weight penalties.

Taylor, R. B.↗
NASA NTRS1994

Aeroacoustic computation of gust-blade interaction

To better understand and address the challenges faced in computing the acoustics of flow fields, test problems must be considered. In the present study, the sound radiated by the interaction of a flat plate with an oncoming gust containing a two component, mean velocity is computed. The gust has a uniform mean flow in x with Mach number M(infinity) equal to 0.5. The gust's mean velocity in y is of smaller amplitude and is given by: v = 0.1 sin(pi/8(x/M(sub infinity) - t)). This problem has been posed for an upcoming ICASE/LaRC workshop on benchmark problems in computational aeroacoustics. A plate with a length of 30 units in x is used. The plate is assumed to be infinitesimally thin and is centered at the origin. All variables are made dimensionless using the scales specified. Acoustic quantities are obtained by numerically integrating the linearized Euler equations. Integration is performed on the computational domain -100.0 less than or equal to x less than or equal to 100.0, -100.0 less than or equal to y less than or equal to 100.0, using unit length grid spacing in x and in y. An integration scheme is sought which will provide accurate solution to the small quantities of interest at a minimal computational expense. Results indicate that with the given discretization a scheme of minimal fourth order accuracy might be adequate to approximate the waves within the given flow. Thus, a variation of the MacCormack scheme with fourth order accuracy in space and second order accuracy in time was chosen. A scheme with sixth order accuracy in space has also been implemented and results compared with those of the fourth order accurate scheme. To ensure no mass flux, zero normal velocity is assigned at the plate. This condition will induce a discontinuity in the pressure across the plate location. Values for the perturbation pressure p' along the surface of the plate are obtained using a one-sided, third order Taylor expansion, such that p'(sub y) = O. In accordance with the Kutta condition, perturbation pressure at the trailing edge is assigned to zero. In the far field, radiation boundary conditions have been implemented. The effectiveness of the far field conditions are validated by computing in a larger computational domain and comparing the results. Early in time, sound waves begin to radiate from the plate. A Doppler effect is observed. After the initial transients disappear, the strongest waves leave the trailing edge at an approximate 45 degree angle. The intensity pattern of pressure fluctuations shows five lobes (of increasing magnitude with increasing downstream direction) emerging from the plate. Undesirable short wave contaminants are observed in the computed pressure distribution along the plate surface. For a more accurate solution at small scales, a more refined discretization will be required.

Martin, James E.↗
NASA NTRS1991

Planetary surface structure and evolution of Mars

The surface of the planet Mars is characterized by two different hemispheres: old densely cratered structures in the Southern Uplands, and sparsely cratered younger topographically lower regions covering approximately the northern third of the planet's surface. The model for explaining these global surface structures is characterized by the following features: (1) homogeneous accretion of Mars; (2) formation of a metal melt layer; (3) northward translation of the undifferentiated core due to a Rayleigh-Taylor instability; (4) vigorous convection in the southern parts and formation of the Southern Uplands' crust; (5) fragmentation of the primordial core, slowly dissolving rockbergs, beginning of Tharsis uplift and volcanism; (6) formation of the Northern Lowlands' crust from only weakly differentiated silicatic material; and (7) reaching of the present state with symmetrically placed core and further thermal evolution.

Franck, Siegfried↗
NASA NTRS1971

Standard deviation of vertical two-point longitudinal velocity differences in the atmospheric boundary layer.

Statistical estimates of wind shear in the planetary boundary layer are important in the design of V/STOL aircraft, and for the design of the Space Shuttle. The data analyzed in this study consist of eleven sets of longitudinal turbulent velocity fluctuation time histories digitized at 0.2 sec intervals with approximately 18,000 data points per time history. The longitudinal velocity fluctuations were calculated with horizontal wind and direction data collected at the 18-, 30-, 60-, 90-, 120-, and 150-m levels. The data obtained confirm the result that Eulerian time spectra transformed to wave-number spectra with Taylor's frozen eddy hypothesis possess inertial-like behavior at wave-numbers well out of the inertial subrange.

Fichtl, G. H.↗