SEARCH · Search NASA

Results for “first order”

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 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗

NASA NTRS1994

Errors induced by the neglect of polarization in radiance calculations for Rayleigh-scattering atmospheres

Although neglecting polarization and replacing the rigorous vector radiative transfer equation by its approximate scalar counterpart has no physical background, it is a widely used simplification when the incident light is unpolarized and only the intensity of the reflected light is to be computed. We employ accurate vector and scalar multiple-scattering calculations to perform a systematic study of the errors induced by the neglect of polarization in radiance calculations for a homogeneous, plane-parallel Rayleigh-scattering atmosphere (with and without depolarization) above a Lambertian surface. Specifically, we calculate percent errors in the reflected intensity for various directions of light incidence and reflection, optical thicknesses of the atmosphere, single-scattering albedos, depolarization factors, and surface albedos. The numerical data displayed can be used to decide whether or not the scalar approximation may be employed depending on the parameters of the problem. We show that the errors decrease with increasing depolarization factor and/or increasing surface albedo. For conservative or nearly conservative scattering and small surface albedos, the errors are maximum at optical thicknesses of about 1. The calculated errors may be too large for some practical applications, and, therefore, rigorous vector calculations should be employed whenever possible. However, if approximate scalar calculations are used, we recommend to avoid geometries involving phase angles equal or close to 0 deg and 90 deg, where the errors are especially significant. We propose a theoretical explanation of the large vector/scalar differences in the case of Rayleigh scattering. According to this explanation, the differences are caused by the particular structure of the Rayleigh scattering matrix and come from lower-order (except first-order) light scattering paths involving right scattering angles and right-angle rotations of the scattering plane.

Mishchenko, M. I.↗

NASA NTRS1995

Blending Velocities In Task Space In Computing Robot Motions

Blending of linear and angular velocities between sequential specified points in task space constitutes theoretical basis of improved method of computing trajectories followed by robotic manipulators. In method, generalized velocity-vector-blending technique provides relatively simple, common conceptual framework for blending linear, angular, and other parametric velocities. Velocity vectors originate from straight-line segments connecting specified task-space points, called "via frames" and represent specified robot poses. Linear-velocity-blending functions chosen from among first-order, third-order-polynomial, and cycloidal options. Angular velocities blended by use of first-order approximation of previous orientation-matrix-blending formulation. Angular-velocity approximation yields small residual error, quantified and corrected. Method offers both relative simplicity and speed needed for generation of robot-manipulator trajectories in real time.

Volpe, Richard A.↗

NASA NTRS2003

Some Advanced Concepts in Discrete Aerodynamic Sensitivity Analysis

An efficient incremental iterative approach for differentiating advanced flow codes is successfully demonstrated on a two-dimensional inviscid model problem. The method employs the reverse-mode capability of the automatic differentiation software tool ADIFOR 3.0 and is proven to yield accurate first-order aerodynamic sensitivity derivatives. A substantial reduction in CPU time and computer memory is demonstrated in comparison with results from a straightforward, black-box reverse-mode applicaiton of ADIFOR 3.0 to the same flow code. An ADIFOR-assisted procedure for accurate second-rder aerodynamic sensitivity derivatives is successfully verified on an inviscid transonic lifting airfoil example problem. The method requires that first-order derivatives are calculated first using both the forward (direct) and reverse (adjoinct) procedures; then, a very efficient noniterative calculation of all second-order derivatives can be accomplished. Accurate second derivatives (i.e., the complete Hesian matrices) of lift, wave drag, and pitching-moment coefficients are calculated with respect to geometric shape, angle of attack, and freestream Mach number.

Taylor, Arthur C., III↗

NASA NTRS1980

Computer program documentation for the dynamic analysis of a noncontacting mechanical face seal

A computer program is presented which achieves a numerical solution for the equations of motion of a noncontacting mechanical face seal. The flexibly-mounted primary seal ring motion is expressed by a set of second order differential equations for three degrees of freedom. These equations are reduced to a set of first order equations and the GEAR software package is used to solve the set of first order equations. Program input includes seal design parameters and seal operating conditions. Output from the program includes velocities and displacements of the seal ring about the axis of an inertial reference system. One example problem is described.

Auer, B. M.↗

NASA NTRS1987

An analytical study of the effects of transverse shear deformation and anisotropy on buckling loads of laminated cylinders

Buckling loads of thick-walled orthotropic and anisotropic simply supported circular cylinders are predicted using a higher-order transverse-shear deformation theory. A comparison of buckling loads predicted by the conventional first-order transverse-shear deformation theory and the higher-order theory show that the additional allowance for transverse shear deformation has a negligible effect on the predicted buckling loads of medium-thick metallic isotropic cylinders. However, the higher-order theory predicts buckling loads which are significantly lower than those predicted by the first-order transverse-shear deformation theory for certain short, thick-walled cylinders which have low through-the-thickness shear moduli. A parametric study of the effects of ply orientation on the buckling load of axially compressed cylinders indicates that laminates containing 45 degree plies are most sensitive to transverse-shear deformation effects. Interaction curves for buckling loads of cylinders subjected to axial compressive and external pressure loadings indicate that buckling loads due to external pressure loadings are as sensitive to transverse-shear deformation effects as buckling loads due to axial compressive loadings. The effects of anisotropy are important over a much wider range of cylinder geometries than the effects of transverse shear deformation.

Jegley, Dawn C.↗

NASA NTRS1996

Thermohydrodynamic Analysis of Cryogenic Liquid Turbulent Flow Fluid Film Bearings

Computational programs developed for the thermal analysis of tilting and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings are described. The motion of a cryogenic liquid on the thin film annular region of a fluid film bearing is described by a set of mass and momentum conservation, and energy transport equations for the turbulent bulk-flow velocities and pressure, and accompanied by thermophysical state equations for evaluation of the fluid material properties. Zeroth-order equations describe the fluid flow field for a journal static equilibrium position, while first-order (linear) equations govern the fluid flow for small amplitude-journal center translational motions. Solution to the zeroth-order flow field equations provides the bearing flow rate, load capacity, drag torque and temperature rise. Solution to the first-order equations determines the rotordynamic force coefficients due to journal radial motions.

Luis San Andres↗

NASA NTRS1977

High-energy heavy-ion scattering and the optical phase shift function

Within the framework of the Glauber approximation, an exact expression for the optical phase shift function is derived which takes into account the center-of-mass correlation function. We calculate total and differential cross sections using the first-order term for the shift function and compare them with the usual (first-order) optical limit results and, for the lighter nuclei, with the exact Glauber results. The nucleus-nucleus elastic scattering amplitude does not exhibit the large-q divergence which characterizes the usual optical limit. The obtained results improve the calculated total and differential cross sections dramatically for light nuclei, and significantly for medium and heavy nuclei.

Franco, V.↗

NASA NTRS1982

Manufacturing in space: Fluid dynamics numerical analysis

Numerical computations were performed for natural convection in circular enclosures under various conditions of acceleration. It was found that subcritical acceleration vectors applied in the direction of the temperature gradient will lead to an eventual state of rest regardless of the initial state of motion. Supercritical acceleration vectors will lead to the same steady state condition of motion regardless of the initial state of motion. Convection velocities were computed for acceleration vectors at various angles of the initial temperature gradient. The results for Rayleigh numbers of 1000 or less were found to closely follow Weinbaum's first order theory. Higher Rayleigh number results were shown to depart significantly from the first order theory. Supercritical behavior was confirmed for Rayleigh numbers greater than the known supercritical value of 9216. Response times were determined to provide an indication of the time required to change states of motion for the various cases considered.

Robertson, S. J.↗

NASA NTRS1996

Three-Dimensional Simulations of Marangoni-Benard Convection in Small Containers by the Least-Squares Finite Element Method

This paper reports a numerical study of the Marangoni-Benard (MB) convection in a planar fluid layer. The least-squares finite element method (LSFEM) is employed to solve the three-dimensional Stokes equations and the energy equation. First, the governing equations are reduced to be first-order by introducing variables such as vorticity and heat fluxes. The resultant first-order system is then cast into a div-curl-grad formulation, and its ellipticity and permissible boundary conditions are readily proved. This numerical approach provides an equal-order discretization for velocity, pressure, vorticity, temperature, and heat conduction fluxes, and therefore can provide high fidelity solutions for the complex flow physics of the MB convection. Numerical results reported include the critical Marangoni numbers (M(sub ac)) for the onset of the convection in containers with various aspect ratios, and the planforms of supercritical MB flows. The numerical solutions compared favorably with the experimental results reported by Koschmieder et al..

Yu, Sheng-Tao↗

NASA NTRS1954

Experimental Determination of Linear Dynamics of Two-Spool Turbojet Engines

Transfer functions descriptive of the response of most engine variables were determined from transient data that were obtained from approximate step inputs in fuel flow and in exhaust-nozzle area. The speed responses of both spools to fuel flow and to turbine-inlet temperature appeared as identical first-order lags. Response to exhaust-nozzle area was characterized by a first-order lag response of the outer-spool speed, accompanied by virtually no change in inner-spool speed.

Novik, David↗

NASA NTRSDate not supplied

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗

NASA NTRS1982

On the dynamic collapse of a column impacting a rigid surface

Results are presented of an analytical investigation on the dynamic collapse of an elastic periodically supported column. The column has an attached mass at one end and impacts a rigid surface with prescribed velocity and angle of incidence at the other end. A first-order approximate nonlinear analysis is developed in which it is assumed that only first-order nonlinear terms need be retained. Differential equations are derived and solved numerically by explicit time integration in order to determine specific ranges of the basic nondimensional parameters for which the response characteristics are markedly different. To assess modeling sophistication relative to prediction accuracy, a comparison is made of the results from this approximate analysis with those of a finite element computer code based on a convected coordinate formulation.

Housner, J. M.↗

NASA NTRS1987

Analysis of implicit local linearization techniques for upwind and TVD algorithms

An attempt is made to investigate local time linearization techniques for implicit flux-difference splitting and flux-vector splitting schemes in the simplest settings (i.e., first-order spatial schemes and one-dimensional Euler flows). It is noted that first-order spatial schemes provide the simplest examples of schemes which are collective extensions of scalar TVD schemes. Simple analytical results concerning the local linearizations are highlighted and subsequently verified using a numerical fixed-point analysis on selected problems. It is noted that while primary emphasis is on asymptotic behavior, many of the results have implications for time-accurate calculations as well.

Barth, Timothy J.↗

NASA NTRS1990

Fringe formation in dual-hologram interferometry

Reference-fringe formation in nondiffuse dual-hologram interferometry is described by combining a first-order geometrical hologram treatment with interference fringes generated by two point sources. The first-order imaging relationships can be used to describe reference-fringe patterns for the geometry of the dual-hologram interferometry. The process can be completed without adjusting the two holograms when the reconstructing wavelength is less than the exposing wavelength, and the process is found to facilitate basic intereferometer adjustments.

Burner, A. W.↗

NASA NTRS1990

Thickness and orientational design for a maximum stiff membrane

Recent results from sensitivity analysis for strain energy with anisotropic elasticity are applied to thickness and orientational design of laminated membranes. Primarily, the first order gradients of the total elastic energy are used in an optimality criteria based method. This traditional method is shown to give slow convergence with respect to design parameters, although the convergence of strain energy is very good. To get a deeper insight into this rather general characteristic, second order derivatives are included and it is shown how they can be obtained by first order sensitivity analysis. Examples of only thickness design, only orientational design, and combined thickness--orientational design are presented.

Pedersen, Pauli↗

NASA NTRS2000

Model of ASTM Flammability Test in Microgravity: Iron Rods

There is extensive qualitative results from burning metallic materials in a NASA/ASTM flammability test system in normal gravity. However, this data was shown to be inconclusive for applications involving oxygen-enriched atmospheres under microgravity conditions by conducting tests using the 2.2-second Lewis Research Center (LeRC) Drop Tower. Data from neither type of test has been reduced to fundamental kinetic and dynamic systems parameters. This paper reports the initial model analysis for burning iron rods under microgravity conditions using data obtained at the LERC tower and modeling the burning system after ignition. Under the conditions of the test the burning mass regresses up the rod to be detached upon deceleration at the end of the drop. The model describes the burning system as a semi-batch, well-mixed reactor with product accumulation only. This model is consistent with the 2.0-second duration of the test. Transient temperature and pressure measurements are made on the chamber volume. The rod solid-liquid interface melting rate is obtained from film records. The model consists of a set of 17 non-linear, first-order differential equations which are solved using MATLAB. This analysis confirms that a first-order rate, in oxygen concentration, is consistent for the iron-oxygen kinetic reaction. An apparent activation energy of 246.8 kJ/mol is consistent for this model.

Steinberg, Theodore A↗

NASA NTRS1991

Integrated analysis and design for controlled nonlinear multibody systems

A hybrid approach for integrated analysis and design of nonlinear multibody systems is proposed. Based on the approach, a general-purpose design tool is developed for nonlinear dynamic systems subjected to nonlinear design constraints. For analysis purposes, second-order nonlinear equations of motion of the system are automatically generated using general-purpose multibody formulations. Once the equations are solved, they are written in first-order form to take advantage of the first-order formulations of design sensitivity analysis for dynamic systems. A nonlinear programming technique is then used to optimize the nonlinear systems, such that a nonlinear cost function is minimized and performance constraints are satisfied. The approach proved to be very general and useful for design and analysis of controlled multibody systems. The approach is applied to the design of passive dynamic controllers for slewing control of multibody systems.

Wu, Shih-Chin↗