SEARCH · Engineering Papers
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.
As an aid in the numerical construction of Hill functions and their derivatives, an algorithm using local coordinates and an expansion in Legendre polynomials is proposed. The algorithm is shown to possess sufficient stability, and the orthogonality of the Legendre polynomials simplifies the computation when the Ritz-Galerkin technique is used.
The dynamics of a helicopter blade in forward flight are described by a system of linear differential equations with periodic coefficients. The stability of this periodic aeroelastic system is determined, using multivariable Floquet-Liapunov theory. The transition matrix at the end of the period is evaluated by: (1) direct numerical integration, and (2) a new, approximate method, which consists in approximating a periodic function by a series of step functions. The numerical accuracy and efficiency of the methods is compared, and the second method is shown to be superior by far. Results illustrating the effect of the periodic coefficients and various blade parameters are presented.
Within the context of the restricted problem of three bodies the effects caused by varying the mass ratio of the primaries and the eccentricity of their orbits, upon periodic orbits of the infinitesimal mass which are numerical continuations of circular orbits in the ordinary problem of two bodies are shown. A recursive power series technique is used to numerically integrate the equations of motion as well as the first variational equations in order to generate a two parameter family of perodic orbits and identify the linear stability characteristics. Seven such families are investigated with equally spaced mass ratios from 0.0 to 1.0 and eccentricities of the orbits of the primaries in a range 0.0 to 0.6. Stable orbits are associated with large distances of the infinitesimal mass from the perturbing primary, nearly circular motion of the primaries, and small mass ratios of the primaries. Unstable orbits for the infinitesimal mass are associated with small distances from the perturbing primary, highly elliptic orbits of the primaries and large mass ratios.
A bump-on-tail unstable reduced velocity distribution, constructed from data obtained at the upstream boundary of the electron foreshock by the GSFC electron spectrometer experiment on the ISEE-1 satellite, is used as the initial plasma state for a numerical integration of the 1D-Vlasov-Maxwell system of equations. The integration is carried through the growth of the instability, beyond its saturation, and well into the stabilized plasma regime. A power spectrum computed for the electric field of the stabilized plasma is dominated by a narrow peak at the Bohm-Gross frequency of the unstable field mode but also contains significant power at the harmonics of the Bohm-Gross frequency. The harmonic power is in sharp peaks which are split into closely spaced doublets. The fundamental peak at the Bohm-Gross frequency is split into a closely spaced triplet. The mechanism for excitation of the second harmonic is shown to be second order wave-wave coupling.
A bump-on-tail unstable reduced velocity distribution, constructed from data obtained at the upstream boundary of the electron foreshock by the GSFC electron spectrometer experiment on the ISEE-1 satellite, is used as the initial plasma state for a numerical integration of the 1D-Vlasov-Maxwell system of equations. The integration is carried through the growth of the instability, beyond its saturation, and well into the stabilized plasma regime. A power spectrum computed for the electric field of the stabilized plasma is dominated by a narrow peak at the Bohm-Gross frequency of the unstable field mode but also contains significant power at the harmonics of the Bohm-Gross frequency. The harmonic power is in sharp peaks which are split into closely spaced doublets. The fundamental peak at the Bohm-Gross frequency is split into a closely spaced triplet. The mechanism for excitation of the second harmonic is shown to be second order wave-wave coupling. Previously announced in STAR as N83-17315
Two families of symmetric periodic orbits of the planar, general, three-body problem are presented. The masses of the three bodies include ratios equal to the sun-Jupiter-Saturn system and the periods of the orbits of Jupiter and Saturn are in a 2:5 resonance. The (linear) stability of the orbits are studied in relation to eccentricity and mass variations. The generation of the two families of periodic orbits follows a systematic approach and employs (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to the elliptic restricted problem and from the circular and elliptic problems to the general problem through bifurcation phenomena relating the three dynamical systems. The approach also provides insight into the evolutionary process of periodic orbits continued from the restricted problems to the general problem.
The formation of a coronal cavity and its relation to a quiescent prominence is studied theoretically. The stability of condensation modes of a plasma in the coronal streamer model (Steinolfson et al., 1982) is considered using a two-dimensional time-dependent ideal MHD numerical simulation. It is found that a plasma with beta = 0.5 is unstable but one with beta = 4 is stable because the density enhancement of the plasma trapped by the closed fields increases with the strength of the magnetic field. The means by which condensation modes can produce a coronal cavity and/or initiate the formation of a prominence (depending on the field configuration) are discussed. It is argued that prominence and cavity material is all supplied from the chromospheric level in the form of spicules.
Primitive variables with central differencing on a staggered grid are used in the present, factored ADI finite-difference scheme for artificial compressibility method solution of the incompressible Navier-Stokes equations, leading to a close coupling between velocity and pressure that both enhances stability and eliminates the need for artificial damping. Computational efficiency is enhanced through the use of a spatially variable, fixed Courant number-based time-step. The numerical results obtained for a driven cavity at Re of 10,000, with local cell Re as high as 100, exhibits no flow variable spatial oscillations on a 40 x 40 stretched grid solution.
In order to assess the state of the art in transonic flow disciplines and to glimpse at future directions, NASA-Langley held a Transonic Symposium. Emphasis was placed on steady, three dimensional external, transonic flow and its simulation, both numerically and experimentally. The symposium included technical sessions on wind tunnel and flight experiments; computational fluid dynamic applications; inviscid methods and grid generation; viscous methods and boundary layer stability; and wind tunnel techniques and wall interference. This, being volume 1, is unclassified.
The stream function equation is solved for steady two-dimensional (and axisymmetric) rotational flows. Both finite differences and finite volumes discretization techniques are studied, using generalized body fitted coordinates and unstructured staggered grids, respectively. For inviscid transonic flows, a new artificial viscosity scheme which does not produce any artificial vorticity is introduced, for the stability of the mixed flow calculations and for capturing shocks. The solution of Euler equations, in primitive variables, are also considered. The effects of the artificial viscosity and numerical boundary conditions on the total enthalpy and the vorticity distributions are demonstrated.
Consideration is given to the effect of compressibility on a plane mixing layer that is a prototype free shear layer, amenable to study by numerical simulation and experiment. The full time-dependent compressible Navier-Stokes equations are solved numerically for a temporally evolving mixing layer employing a mixed spectral and high-order finite difference method. Simulations with random initial conditions confirm the prediction of linear stability theory that at high Mach numbers oblique waves grow faster than two-dimensional waves. Simulations are presented of the nonlinear temporal evolution of the most rapidly amplified linear instability waves.
A frequency-response method for determining the critical control-gearing and hunting oscillations of airplanes with automatic pilots is presented. The method is graphical and has several advantages over the standard numerical procedure based on Routh's discriminant. The chief advantage of the method is that direct use can be made of the measured response characteristics of the automatic pilot. This feature is especially useful in determining the existence, amplitude, and frequency of the hunting oscillations that may be present when the automatic pilot has nonlinear dynamic characteristics. Several examples are worked out to illustrate the application of the frequency-response method in determining the effect of automatic-pilot lag or lead on critical control gearing and in determining the amplitude and frequency hunting. It is shown that the method may be applied to the case of a control geared to airplane motions about two axes.
We examine the evolution of baroclinic vortices in a time-dependent, nonlinear numerical model of a Jovian atmosphere. The model uses a normal-mode expansion in the vertical, using the barotropic and first two baroclinic modes. Results for the stability of baroclinic vortices on an f plane in the absence of a mean zonal flow are similar to results of Earth vortex models, although the presence of a fluid interior on the Jovian planets shifts the stability boundaries to smaller length scales. The presence of a barotropic mean zonal flow in the interior stabilizes vortices against instability and significantly modifies the finite amplitude form of baroclinic instabilities. The effect of a zonal flow on a form of barotropic instability produces periodic oscillations in the latitude and longitude of the vortex as observed at the level of the cloud tops. This instability may explain some, but not all, observations of longitudinal oscillations of vortices on the outer planets. Oscillations in aspect ratio and orientation of stable vortices in a zonal shear flow are observed in this baroclinic model, as in simpler two-dimensional models. Such oscillations are also observed in the atmospheres of Jupiter and Neptune. The meridional propagation and decay of vortices on a beta plane is inhibited by the presence of a mean zonal flow. The direction of propagation of a vortex relative to the mean zonal flow depends upon the sign of the meridional potential vorticity gradient; combined with observations of vortex drift rates, this may provide a constraint on model assumption for the flow in the deep interior of the Jovian planets.
This paper has a dual purpose: it presents a multidomain Chebyshev method for the solution of the two-dimensional reactive compressible Navier-Stokes equations, and it reports the results of the application of this code to the numerical simulations of high Mach number reactive flows in recessed cavity. The computational method utilizes newly derived interface boundary conditions as well as an adaptive filtering technique to stabilize the computations. The results of the simulations are relevant to recessed cavity flameholders.
A forced oscillation test of the Orion Crew Module (CM) was conducted in the Langley 20-Foot Vertical Spin Tunnel. The objective of the test was to quantify the rate damping characteristics of the CM-drogue chute system. Numerous configurations were tested to measure the influence of the chutes on the CM dynamic aerodynamics and, conversely, the influence of the CM on drogue performance. Results show that the CM-drogue system is well-damped at all combinations of frequency, amplitude, and Strouhal number. The wake of the CM significantly reduces the drogue chute riser line force, and the drogues have little upstream influence on the CM aerodynamics. These results are being used to improve simulation model fidelity of CM flight with drogues deployed, which has been identified by the project as key to a successful Orion Critical Design Review.
Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.
A computer program is developed to solve the compressible, laminar boundary-layer equations for two-dimensional flow, axisymmetric flow, and quasi-three-dimensional flows including the flow along the plane of symmetry, flow along the leading-edge attachment line, and swept-wing flows with a conical flow approximation. The finite-difference numerical procedure used to solve the governing equations is second-order accurate. The flow over a wide range of speed, from subsonic to hypersonic speed with perfect gas assumption, can be calculated. Various wall boundary conditions, such as wall suction or blowing and hot or cold walls, can be applied. The results indicate that this boundary-layer code gives velocity and temperature profiles which are accurate, smooth, and continuous through the first and second normal derivatives. The code presented herein can be coupled with a stability analysis code and used to predict the onset of the boundary-layer transition which enables the assessment of the laminar flow control techniques. A user's manual is also included.
A hyperbolic initial-boundary-value problem can be approximated by a system of ordinary differential equations (ODEs) by replacing the spatial derivatives by finite-difference approximations. The resulting system of ODEs is called a semidiscrete approximation. A complication is the fact that more boundary conditions are required for the spatially discrete approximation than are specified for the partial differential equation. Consequently, additional numerical boundary conditions are required and improper treatment of these additional conditions can lead to instability. For a linear initial-boundary-value problem (IBVP) with homogeneous analytical boundary conditions, the semidiscrete approximation results in a system of ODEs of the form du/dt = Au whose solution can be written as u(t) = exp(At)u(O). Lax-Richtmyer stability requires that the matrix norm of exp(At) be uniformly bounded for O less than or = t less than or = T independent of the spatial mesh size. Although the classical Lax-Richtmyer stability definition involves a conventional vector norm, there is no known algebraic test for the uniform boundedness of the matrix norm of exp(At) for hyperbolic IBVPs. An alternative but more complicated stability definition is used in the theory developed by Gustafsson, Kreiss, and Sundstrom (GKS). The two methods are compared.