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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Numerical simulation of the non-Newtonian mixing layer

This work is a continuing effort to advance our understanding of the effects of polymer additives on the structures of the mixing layer. In anticipation of full nonlinear simulations of the non-Newtonian mixing layer, we examined in a first stage the linear stability of the non-Newtonian mixing layer. The results of this study show that, for a fluid described by the Oldroyd-B model, viscoelasticity reduces the instability of the inviscid mixing layer in a special limit where the ratio (We/Re) is of order 1 where We is the Weissenberg number, a measure of the elasticity of the flow, and Re is the Reynolds number. In the present study, we pursue this project with numerical simulations of the non-Newtonian mixing layer. Our primary objective is to determine the effects of viscoelasticity on the roll-up structure. We also examine the origin of the numerical instabilities usually encountered in the simulations of non-Newtonian fluids.

Azaiez, Jalel↗
NASA NTRS1977

Unsteady aerodynamic modeling for arbitrary motions

A study is presented on the unsteady aerodynamic loads due to arbitrary motions of a thin wing and their adaptation for the calculation of response and true stability of aeroelastic modes. In an Appendix, the use of Laplace transform techniques and the generalized Theodorsen function for two-dimensional incompressible flow is reviewed. New applications of the same approach are shown also to yield airloads valid for quite general small motions. Numerical results are given for the two-dimensional supersonic case. Previously proposed approximate methods, starting from simple harmonic unsteady theory, are evaluated by comparison with exact results obtained by the present approach. The Laplace inversion integral is employed to separate the loads into 'rational' and 'nonrational' parts, of which only the former are involved in aeroelastic stability of the wing. Among other suggestions for further work, it is explained how existing aerodynamic computer programs may be adapted in a fairly straightforward fashion to deal with arbitrary transients.

Edwards, J. W.↗
NASA NTRS1971

Analysis of cavitation bubble dynamics in a liquid

General differential equations governing the dynamics of the cavitation bubbles in a liquid were derived. With the assumption of spherical symmetry the governing equations were simplified. Closed form solutions were obtained for simple cases, and numerical solutions were calculated for complicated ones. The growth and the collapse of the bubble were analyzed, oscillations of the bubbles were studied, and the stability of the cavitation bubbles were investigated. The results show that the cavitation bubbles are unstable, and the oscillation is not sinusoidal.

Fontenot, L. L.↗
NASA NTRS1983

Computation of viscous flows over airfoils, including separation, with a coupling approach

Viscous incompressible flows over single or multiple airfoils, with or without separation, were computed using an inviscid flow calculation, with modified boundary conditions, and by a method providing calculation and coupling for boundary layers and wakes, within conditions of strong viscous interaction. The inviscid flow is calculated with a method of singularities, the numerics of which were improved by using both source and vortex distributions over profiles, associated with regularity conditions for the fictitious flows inside of the airfoils. The viscous calculation estimates the difference between viscous flow and inviscid interacting flow, with a direct or inverse integral method, laminar or turbulent, with or without reverse flow. The numerical method for coupling determines iteratively the boundary conditions for the inviscid flow. For attached viscous layers regions, an underrelaxation is locally calculated to insure stability. For separated or separating regions, a special semi-inverse algorithm is used. Comparisons with experiments are presented.

Leballeur, J. C.↗
NASA NTRS1988

Stability analysis of large space structure control systems with delayed input

Large space structural systems, due to their inherent flexibility and low mass to area ratio, are represented by large dimensional mathematical models. For implementation of the control laws for such systems a finite amount of time is required to evaluate the control signals; and this time delay may cause instability in the closed loop control system that was previously designed without taking the input delay into consideration. The stability analysis of a simple harmonic oscillator representing the equation of a single mode as a function of delay time is treated analytically and verified numerically. The effect of inherent damping on the delay is also analyzed. The control problem with delayed input is also formulated in the discrete time domain.

Reddy, A. S. S. R.↗
NASA NTRS1988

Stability analysis of large space structure control systems with delayed input

Large space structural systems, due to their inherent flexibility and low mass to area ratio, are represented by large dimensional mathematical models. For implementation of the control laws for such systems a finite amount of time is required to evaluate the control signals; and this time delay may cause instability in the closed loop control system that was previously designed without taking the input delay into consideration. The stability analysis of a simple harmonic oscillator representing the equation of a single mode as a function of delay time is analyzed analytically and verified numerically. The effect of inherent damping on the delay is also analyzed. The control problem with delayed input is also formulated in the discrete time domain.

Reddy, A. S. S. R.↗
NASA NTRS1989

A simple element for aeroelastic analysis of undamaged and damaged wings

The effects of material damage on the aeroelastic response of an anisotropic composite wing are investigated analytically. The wing is modeled as a Timoshenko beam of arbitrary cross section, applying a 24-DOF one-dimensional beam element and Cartesian coordinates. The formulations for the element stiffness matrix, the strain-displacement relations, the box-beam constitutive relations, the contributions of stringers and web, and the warping correction factor are outlined; the application of strip theory in deriving the aerodynamic model is briefly explained; and numerical results are presented in extensive graphs. It is shown that asymmetry and damage can induce bending-stretching coupling phenomena and thereby affect the stability of the aircraft either positively or negatively. The need to take such effects into account in the design of advanced composite wings is stressed.

Kapania, Rakesh K.↗
NASA NTRS1989

On an efficient and accurate method to integrate restricted three-body orbits

This work is a quantitative analysis of the advantages of the Bulirsch-Stoer (1966) method, demonstrating that this method is certainly worth considering when working with small N dynamical systems. The results, qualitatively suspected by many users, are quantitatively confirmed as follows: (1) the Bulirsch-Stoer extrapolation method is very fast and moderately accurate; (2) regularization of the equations of motion stabilizes the error behavior of the method and is, of course, essential during close approaches; and (3) when applicable, a manifold-correction algorithm reduces numerical errors to the limits of machine accuracy. In addition, for the specific case of the restricted three-body problem, even a small eccentricity for the orbit of the primaries drastically affects the accuracy of integrations, whether regularized or not; the circular restricted problem integrates much more accurately.

Murison, Marc A.↗
NASA NTRS1990

A multigrid solver for semi-implicit global shallow-water models

A multigrid solver is developed for the discretized two-dimensional elliptic equation on the sphere that arises from a semiimplicit time discretization of the global shallow-water equations. Different formulations of the semiimplicit scheme result in variable-coefficient Helmholtz-type equations for which no fast direct solvers are available. The efficiency of the multigrid solver is optimal, in the sense that the total operation count is proportional to the number of unknowns. Numerical experiments using initial data derived from actual 300-mb height and wind velocity fields indicate that the present model has very good accuracy and stability properties.

Barros, Saulo R. M.↗
NASA NTRS1990

Instabilities in rotating channel flow

A combination of linear, weakly nonlinear, and secondary stability theories has been used to study the transition process in rotating channel flow. Results for parts of the theory have been confirmed by direct numerical simulations. Preliminary results indicate that the theory can help explain some experimental observations regarding the evolution of streamwise vortices.

Ng, Lian↗
NASA NTRS1991

An unconditionally stable staggered algorithm for transient finite element analysis of coupled thermoelastic problems

An unconditionally stable second order accurate implicit-implicit staggered procedure for the finite element solution of fully coupled thermoelasticity transient problems is proposed. The procedure is stabilized with a semi-algebraic augmentation technique. A comparative cost analysis reveals the superiority of the proposed computational strategy to other conventional staggered procedures. Numerical examples of one and two-dimensional thermomechanical coupled problems demonstrate the accuracy of the proposed numerical solution algorithm.

Farhat, C.↗
NASA NTRS1991

An unconditionally stable staggered algorithm for transient finite element analysis of coupled thermoelastic problems

An unconditionally stable second order accurate implicit-implicit staggered procedure for the finite element solution of fully coupled thermoelasticity transient problems is proposed. The procedure is stabilized with a semi-algebraic augmentation technique. A comparative cost analysis reveals the superiority of the proposed computational strategy to other conventional staggered procedures. Numerical examples of one- and two-dimensional thermomechanical coupled problems demonstrate the accuracy of the proposed numerical solution algorithm.

Farhat, Charbel↗
NASA NTRS1997

Noise from Supersonic Coaxial Jets: Inverted Velocity Profile - Part 3

The instability wave noise generation model is used to study the instability waves in the two shear layers of an inverted velocity profile, supersonic, coaxial jet and the noise radiated from the dominant wave. The inverted velocity profile jet has a high speed outer stream surrounding a low speed inner stream and the outer shear layer is always larger than the inner shear layer. The jet mean flows are calculated numerically. The operating conditions are chosen to exemplify the effect of the coaxial jet outer shear layer initial spreading rates. 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. Results for inverted velocity profile jets indicate that relative maximum instability wave amplitudes and far field peak noise levels can be reduced from that of the reference jet by having higher spreading rates for the outer shear layer, low velocity ratios, and outer streams hotter than the inner stream.

Dahl, Milo D.↗
NASA NTRS1980

Implicit Finite-Difference Simulations of Three-Dimensional Compressible Flow

An implicit finite-difference procedure for unsteady three-dimensional flow capable of handling arbitrary geometry through the use of general coordinate transformations is described. Viscous effects are optionally incorporated with a "thin-layer" approximation of the Navier-Stokes equations. An implicit approximate factorization technique is employed so that the small grid sizes required for spatial accuracy and viscous resolution do not impose stringent stability limitations. Results obtained from the program include transonic inviscid or viscous solutions about simple body configurations. Comparisons with existing theories and experiments are made. Numerical accuracy and the effect of three-dimensional coordinate singularities are also discussed.

Pulliam, Thomas H.↗
NASA NTRS2003

The Effect of Stabilization Heat Treatments on the Tensile and Creep Behavior of an Advanced Nickel-Based Disk Alloy

As part of NASA s Advanced Subsonic Technology Program, a study of stabilization heat treatment options for an advanced nickel-base disk alloy, ME 209, was performed. Using a simple, physically based approach, the effect of stabilization heat treatments on tensile and creep properties was analyzed in this paper. Solutions temperature, solution cooling rate, and stabilization temperature/time were found to have a significant impact on tensile and creep properties. These effects were readily quantified using the following methodology. First, the effect of solution cooling rate was assessed to determine its impact on a given property. The as-cooled property was then modified by using two multiplicative factors which assess the impact of solution temperature and stabilization parameters. Comparison of experimental data with predicted values showed this physically based analysis produced good results that rivaled the statistical analysis employed, which required numerous changes in the form of the regression equation depending on the property and temperature in question. As this physically based analysis uses the data for input, it should be noted that predictions which attempt to extrapolate beyond the bounds of the data must be viewed with skepticism. Future work aimed at expanding the range of the stabilization/aging parameters explored in this study would be highly desirable, especially at the higher solution cooling rates.

Gayda, John↗
NASA NTRS2014

All-Digital Baseband 65nm PLL/FPLL Clock Multiplier using 10-cell Library

PLLs for clock generation are essential for modern circuits, to generate specialized frequencies for many interfaces and high frequencies for chip internal operation. These circuits depend on analog circuits and careful tailoring for each new process, and making them fault tolerant is an incompletely solved problem. Until now, all digital PLLs have been restricted to sampled data DSP techniques and not available for the highest frequency baseband applications. This paper presents the design and preliminary evaluation of an all-digital baseband technique built entirely with an easily portable 10-cell digital library. The library is also described, as it aids in research and low volume design porting to new processes. The advantages of the digital approach are the wide variety of techniques available to give varying degrees of fault tolerance, and the simplicity of porting the design to new processes, even to exotic processes that may not have analog capability. The only tuning parameter is digital gate delay. An all-digital approach presents unique problems and standard analog loop stability design criteria cannot be directly used. Because of the quantization of frequency, there is effectively infinite gain for very small loop error feedback. The numerically controlled oscillator (NCO) based on a tapped delay line cannot be reliably updated while a pulse is active in the delay line, and ordinarily does not have enough frequency resolution for a low-jitter output.

Shuler, Robert L., Jr.↗
NASA NTRS2014

ALL-Digital Baseband 65nm PLL/FPLL Clock Multiplier Using 10-Cell Library

PLLs for clock generation are essential for modern circuits, to generate specialized frequencies for many interfaces and high frequencies for chip internal operation. These circuits depend on analog circuits and careful tailoring for each new process, and making them fault tolerant is an incompletely solved problem. Until now, all digital PLLs have been restricted to sampled data DSP techniques and not available for the highest frequency baseband applications. This paper presents the design and preliminary evaluation of an all-digital baseband technique built entirely with an easily portable 10-cell digital library. The library is also described, as it aids in research and low volume design porting to new processes. The advantages of the digital approach are the wide variety of techniques available to give varying degrees of fault tolerance, and the simplicity of porting the design to new processes, even to exotic processes that may not have analog capability. The only tuning parameter is digital gate delay. An all-digital approach presents unique problems and standard analog loop stability design criteria cannot be directly used. Because of the quantization of frequency, there is effectively infinite gain for very small loop error feedback. The numerically controlled oscillator (NCO) based on a tapped delay line cannot be reliably updated while a pulse is active in the delay line, and ordinarily does not have enough frequency resolution for a low-jitter output.

Schuler, Robert L., Jr.↗
NASA NTRS2019

An Integrated Approach to Weather Radar Calibration and Monitoring Using Ground Clutter and Satellite Comparisons

The stability and accuracy of weather radar reflectivity calibration are imperative for quantitative applications, such as rainfall estimation, severe weather monitoring and nowcasting, and assimilation in numerical weather prediction models. Various radar calibration and monitoring techniques have been developed, but only recently have integrated approaches been proposed, that is, using different calibration techniques in combination. In this paper the following three techniques are used: 1) ground clutter monitoring, 2) comparisons with spaceborne radars, and 3) the self-consistency of polarimetric variables. These techniques are applied to a C-band polarimetric radar (CPOL) located in the Australian tropics since 1998. The ground clutter monitoring technique is applied to each radar volumetric scan and provides a means to reliably detect changes in calibration, relative to a baseline. It is remarkably stable to within a standard deviation of 0.1 dB (decibels). To obtain an absolute calibration value, CPOL observations are compared to spaceborne radars on board TRMM (Tropical Rainfall Measuring Mission) and GPM (Global Precipitation Measurement) using a volume-matching technique. Using an iterative procedure and stable calibration periods identified by the ground echoes technique, we improve the accuracy of this technique to about 1 dB. Finally, we review the self-consistency technique and constrain its assumptions using results from the hybrid TRMM-GPM and ground echo technique. Small changes in the self-consistency parameterization can lead to 5 dB of variation in the reflectivity calibration. We find that the drop-shape model of Brandes et al. with a standard deviation of the canting angle of 12 degrees best matches our dataset.

Louf, Valentin↗