Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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 937 records · Page 52

Error analysis for pulsed maneuvers of a dual-spin spacecraft

The Galileo is a dual-spin interplanetary spacecraft which is scheduled for launch toward Jupiter in October of 1989. The spacecraft has twelve thrusters, located on the rotor, which can be operated in continuous or pulsed modes. In this paper probabilistic error models are presented for the pulsed maneuvers which are used to change the spacecraft's velocity. These models take into account all error sources which might affect the maneuver accuracies, such as wobble, nutation, plume impingement, thrust level variation and misalignment, gyro drift, burn timing errors and mass property uncertainties. The analytic models are sufficiently general so that they may be applied to all spin and dual-spin stabilized spacecraft. When these models are applied to the case of the Galileo spacecraft, the numerical results demonstrate that the spacecraft can achieve the maneuver accuracies required in order to successfully reach Jupiter and perform its scientific mission.

Longuski, James M.↗

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. The present results indicate that viscous effects are of primary importance in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure and tubulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. Two numerical techniques were tested: vector extrapolation, to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices which arise in chemically reacting flows.

Yungster, Shaye↗

The three-dimensional evolution of a plane mixing layer. Part 1: The Kelvin-Helmholtz roll-up

The Kelvin Helmholtz roll up of three dimensional, temporally evolving, plane mixing layers were simulated numerically. All simulations were begun from a few low wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity profile. The spanwise disturbance wavelength was taken to be less than or equal to the streamwise wavelength associated with the Kelvin Helmholtz roll up. A standard set of clean structures develop in most of the simulations. The spanwise vorticity rolls up into a corrugated spanwise roller, with vortex stretching creating strong spanwise vorticity in a cup shaped region at the vends of the roller. Predominantly streamwise rib vortices develop in the braid region between the rollers. For sufficiently strong initial three dimensional disturbances, these ribs collapse into compact axisymmetric vortices. The rib vortex lines connect to neighboring ribs and are kinked in the opposite direction of the roller vortex lines. Because of this, these two sets of vortex lines remain distinct. For certain initial conditions, persistent ribs do not develop. In such cases the development of significant three dimensionality is delayed. When the initial three dimensional disturbance energy is about equal to, or less than, the two dimensional fundamental disturbance energy, the evolution of the three dimensional disturbance is nearly linear (with respect to the mean and the two dimensional disturbances), at least until the first Kelvin Helmholtz roll up is completed.

Rogers, Michael M.↗

Real parameter margin computation for uncertain structural dynamic systems

This paper is concerned with the problem of computing real parameter margins for stabilized, structural dynamic systems with the masses, spring constants, and damping constants as uncertain parameters. Numerical algorithms for uncertain systems with multilinear parameters are investigated. It is shown that for a certain class of structural dynamic systems with small passive damping (e.g., flexible structures in space), the computational complexity of the problem can be avoided by modeling the system as a conservative plant, without loss of any practical significance.

Wie, Bong↗

New Directions in Phthalocyanine Pigments

The objectives were the following: (1) investigation of the synthesis of new phthalocyanines; (2) characterization of the new phthalocyanines synthesized; (3) investigate the properties of the newly synthesized phthalocyanines with emphasis on UV protection of plastics and coatings; and (4) utilize quantum mechanics to evaluate the structural relationships with possible properties and synthetic approaches. The proposed research targeted the synthesis of phthalocyanines containing an aromatic bridge between two phthalocyanine rings. The goal was to synthesize pigments which would protect plastics when exposed to the photodegradation effects of the sun in space. The stability and extended conjugation of the phthalocyanines offer a unique opportunity for energy absorption and numerous radiative and non-radiative energy loss mechanisms. Although the original targeted phthalocyanines were changed early in the project, several new and unique phthalocyanine compounds were prepared. The basic goals of this work were met and some unique and unexpected outcomes of the work were the result of the integral use of quantum mechanics and molecular modeling with the synthetic effort.

Vandemark, Michael R.↗

Dynamics and Structure of Weakly-Strained Flames at Normal and Microgravity

Experimental and numerical studies were conducted for weakly-strained, laminar premixed flames. The dynamic response and stability of such flames was assessed for a large number of mixtures. A new technique is proposed for the direct experimental determination of laminar flame speeds at the limit of near-zero strain rate.

Vagelopoulos, Christine M.↗

Noise from Supersonic Coaxial Jets: Normal Velocity Profile - Part 2

Instability waves have been established as noise generators in supersonic jets. Recent analysis of these slowly diverging jets has shown that these instability waves radiate noise to the far field when the waves have components with phase velocities that are supersonic relative to the ambient speed of sound. This instability wave noise generation model has been applied to supersonic jets with a single shear layer and is now applied to supersonic coaxial jets with two initial shear layers. In this paper the case of coaxial jets with normal velocity profiles is considered, where the inner jet stream velocity is higher than the outer jet stream velocity. To provide mean flow profiles at all axial locations, a numerical scheme is used to calculate the mean flow properties. Calculations are made for the stability characteristics in the coaxial jet shear layers and the noise radiated from the instability waves for different operating conditions with the same total thrust, mass flow and exit area as a single reference jet. The effects of changes in the velocity ratio, the density ratio and the area ratio are each considered independently.

Dahl, M. D.↗

MHD Turbulence Sheared in Fixed and Rotating Frames

We consider homogeneous turbulence in a conducting fluid that is exposed to a uniform external magnetic field while being sheared in fixed and rotating frames. We take both the frame-rotation axis and the applied magnetic field to be aligned in the direction normal to the plane of the mean shear. Here a systematic parametric study is carried out in a series of Direct Numerical Simulations (DNS) in order to clarify the main effects determining the structural anisotropy and stability of the flow.

Kassinos, S. C.↗

Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains

Borrowing from techniques developed for conservation law equations, numerical schemes which discretize the Hamilton-Jacobi (H-J), level set, and Eikonal equations on triangulated domains are presented. The first scheme is a provably monotone discretization for certain forms of the H-J equations. Unfortunately, the basic scheme lacks proper Lipschitz continuity of the numerical Hamiltonian. By employing a virtual edge flipping technique, Lipschitz continuity of the numerical flux is restored on acute triangulations. Next, schemes are introduced and developed based on the weaker concept of positive coefficient approximations for homogeneous Hamiltonians. These schemes possess a discrete maximum principle on arbitrary triangulations and naturally exhibit proper Lipschitz continuity of the numerical Hamiltonian. Finally, a class of Petrov-Galerkin approximations are considered. These schemes are stabilized via a least-squares bilinear form. The Petrov-Galerkin schemes do not possess a discrete maximum principle but generalize to high order accuracy.

Barth, Timothy J.↗

Receptivity of Hypersonic Boundary Layers to Acoustic and Vortical Disturbances (Invited)

Boundary-layer receptivity to two-dimensional acoustic and vortical disturbances for hypersonic flows over two-dimensional and axi-symmetric geometries were numerically investigated. The role of bluntness, wall cooling, and pressure gradients on the receptivity and stability were analyzed and compared with the sharp nose cases. It was found that for flows over sharp nose geometries in adiabatic wall conditions the instability waves are generated in the leading-edge region and that the boundary layer is much more receptive to slow acoustic waves as compared to the fast waves. The computations confirmed the stabilizing effect of nose bluntness and the role of the entropy layer in the delay of boundary layer transition. The receptivity coefficients in flows over blunt bodies are orders of magnitude smaller than that for the sharp cone cases. Wall cooling stabilizes the first mode strongly and destabilizes the second mode. However, the receptivity coefficients are also much smaller compared to the adiabatic case. The adverse pressure gradients increased the unstable second mode regions.

Balakumar, P.↗

Recent Checks on the Radiometric and Spatial Calibration of AIRS In-Orbit

The Atmospheric Infrared Sounder (AIRS) on the EOS Aqua Spacecraft was launched on May 4, 2002 and is currently fully operational. AIRS acquires hyperspectral infrared radiances in 2378 channels ranging in wavelength from 3.7-15.4 um with spectral resolution of better than 1200, and spatial resolution of 13.5 km with global daily coverage. The AIRS was designed to measure temperature and water vapor profiles for improvement in weather forecast and improved parameterization of climate processes. Currently the AIRS Level 1B Radiance Products are assimilated by NWP centers worldwide and have shown considerable forecast improvement. Although the calibration of AIRS (< 200 mK 3 sigma) is sufficient for data assimilation into Numerical Weather Prediction (NWP) models, long term trends of Earth’s climate require radiances with stability approaching 10 mK/year, and absolute accuracies better than 100 mK. This investigation uses views of space during roll maneuvers of the Aqua spacecraft to calibrate the mirror emission (one of the largest error sources for AIRS) and reduce the residual errors in cold scenes. We also present results of a secondary study that uses MODIS data to determine the alignment of the AIRS boresight. In this study we match AIRS and MODIS data and iterate on the assumed boresight to find the minimum difference in signal. In this way we are able to confirm the boresight projections determined shortly after launch.

Pagano, Thomas S.↗

On the linear stability of compressible plane Couette flow

The linear stability of compressible plane Couette flow is investigated. The correct and proper basic velocity and temperature distributions are perturbed by a small amplitude normal mode disturbance. The full small amplitude disturbance equations are solved numerically at finite Reynolds numbers, and the inviscid limit of these equations is then investigated in some detail. It is found that instability can occur, although the stability characteristics of the flow are quite different from unbounded flows. The effects of viscosity are also calculated, asymptotically, and shown to have a stabilizing role in all the cases investigated. Exceptional regimes to the problem occur when the wavespeed of the disturbances approaches the velocity of either of the walls, and these regimes are also analyzed in some detail. Finally, the effect of imposing radiation-type boundary conditions on the upper (moving) wall (in place of impermeability) is investigated, and shown to yield results common to both bounded and unbounded flows.

Duck, Peter W.↗

The stability of a compressible stratified shear layer

The stability of a shear layer under the effect of gravity is investigated using the compressible magnetohydrodynamic (MHD) equations, including an effective gravity term to represent the curvature effects of the flow and magnetic field line geometry. A general eigenmode equation is derived for a two-dimensional MHD fluid, and an energy-principle analysis to explain the effect of compressibility on the critical Richardson number is presented. For the case of a hyperbolic tangent shear flow and exponential density profile, it was found that, in the Boussinesq approximation, the compressibility raises the critical Richardson number from 1/4 to as much as 1/2, with the exact value depending on the value of the magnetic field at infinity. Under approximation of a strong asymptotic magnetic field, without invoking the Boussinesq approximation, it is shown both analytically and numerically that the density gradient terms cause the shear instability to be dispersive. The long-wavelength stability boundary for the Richardson number J = 0 is characterized by a normalized phase velocity c =

Wang, Z.↗

On the wall-normal velocity of the compressible boundary-layer equations

Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x,y) plane to a computational (xi,eta) plane in which the evolution of the flow is 'slow' in the time-like xi direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem for a wide class of non-similar boundary-layer flows, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently nonlinear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. In light of recent research which shows mean-flow non-parallelism to significantly influence the stability of high-speed compressible flows, the contribution of the wall-normal velocity in the analysis of stability should not be routinely neglected. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally-accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly-accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which is extraordinarily smooth and accurate, and which satisfies the continuity equation nearly to machine precision. These qualities make the method well suited to the computation of the non-parallel mean flows needed by spatial direct numerical simulations (DNS) and parabolized stability equation (PSE) approaches to the analysis of stability.

Pruett, C. David↗

Higher-order accurate Osher schemes with application to compressible boundary layer stability

Two fourth order accurate Osher schemes are presented which maintain higher order accuracy on nonuniform grids. They use either a conservative finite difference or finite volume discretization. Both methods are successfully used for direct numerical simulations of flat plate boundary layer instability at different Mach numbers. Results of growth rates of Tollmien-Schlichting waves compare well with direct simulations of incompressible flow and for compressible flow with results obtained by solving the parabolic stability equations.

Van Der Vegt, J. J. W.↗

Higher-order accurate Osher schemes with application to compressible boundary layer stability

Two fourth order accurate Osher schemes are presented which maintain higher order accuracy on nonuniform grids. They use either a conservative finite difference or finite volume discretization. Both methods are successfully used for direct numerical simulations of flat plate boundary layer instability at different Mach numbers. Results of growth rates of Tollmien-Schlichting waves compare well with direct simulations of incompressible flow and for compressible flow with results obtained by solving the parabolic stability equations.

Vandervegt, J. J. W.↗

Effect of periodic accelerations on interface stability in a multilayered fluid configuration

The increasing number of research opportunities in a microgravity environment will benefit not only fundamental studies in fluid dynamics, but also technological applications such as those involving materials processing. In particular, fluid configurations that involve fluid-fluid interfaces would occur in a variety of experimental investigations. This work investigates the stability of a configuration involving fluid-fluid interfaces in the presence of a time-dependent (periodic) forcing. The fluid configuration is multilayered and infinite in extent. The analysis is linear and inviscid, and the acceleration vector is oriented perpendicular to each interface. A Floquent analysis is employed, and the resulting algebraic eigensystem is truncated. Nondimensional parameters appear in the algebraic system. A numerical study is performed to elucidate the regions of instability and the effect of parameter variation on the fluid configuration stability.

Lyell, M. J.↗

Research study on stabilization and control: Modern sampled-data control theory. Design of the large space telescope system

Conditions of self-sustained oscillations in a two-axis model of the nonlinear LST system are studied. The describing function of the CMG frictional nonlinearity of the LST system is used for the analysis, as well as continuous-data and discrete-data models of the simplified LST control system. A numerical-iterative method is described for the analysis of the two-axis system. Approximation methods and the direct plotting of the stability equation are implemented in the study. It is shown that although the dynamics of the two axes are identical, the amplitudes of self-sustained oscillations in the two axes may in principle be different. Analysis shows that the LST systems are of equal amplitudes but with 180-degree phase shift.

Kuo, B. C.↗