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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Potential Flow Interactions With Directional Solidification

The effect of convective melt motion on the growth of morphological instabilities in crystal growth has been the focus of many studies in the past decade. While most of the efforts have been directed towards investigating the linear stability aspects, relatively little attention has been devoted to experimental and numerical studies. In a pure morphological case, when there is no flow, morphological changes in the solid-liquid interface are governed by heat conduction and solute distribution. Under the influence of a convective motion, both heat and solute are redistributed, thereby affecting the intrinsic morphological phenomenon. The overall effect of the convective motion could be either stabilizing or destabilizing. Recent investigations have predicted stabilization by a flow parallel to the interface. In the case of non-parallel flows, e.g., stagnation point flow, Brattkus and Davis have found a new flow-induced morphological instability that occurs at long wavelengths and also consists of waves propagating against the flow. Other studies have addressed the nonlinear aspects (Konstantinos and Brown, Wollkind and Segel)). In contrast to the earlier studies, our present investigation focuses on the effects of the potential flow fields typically encountered in Hele-Shaw cells. Such a Hele-Shaw cell can simulate a gravity-free environment in the sense that buoyancy-driven convection is largely suppressed, and hence negligible. Our interest lies both in analyzing the linear stability of the solidification process in the presence of potential flow fields, as well as in performing high-accuracy nonlinear simulations. Linear stability analysis can be performed for the flow configuration mentioned above. It is observed that a parallel potential flow is stabilizing and gives rise to waves traveling downstream. We have built a highly accurate numerical scheme which is validated at small amplitudes by comparing with the analytically predicted results for the pure morphological case. We have been able to observe nonlinear effects at larger times. Preliminary results for the case when flow is imposed also provide good validation at small amplitudes.

Buddhavarapu, Sudhir S.↗
NASA NTRS1988

Influence of heating rate on the condensational instability

Analysis and numerical simulation are used to determine the effect that various heating rates have on the linear and nonlinear evolution of a typical plasma within a solar magnetic flux tube subject to the condensational instability. It is found that linear stability depends strongly on the heating rate. The results of numerical simulations of the nonlinear evolution of the condensational instability in a solar magnetic flux tube are presented. Different heating rates lead to quite different nonlinear evolutions, as evidenced by the behavior of the global internal energy.

Dahlburg, R. B.↗
NASA NTRS1984

Boundary-layer linear stability theory

Most fluid flows are turbulent rather than laminar and the reason for this was studied. One of the earliest explanations was that laminar flow is unstable, and the linear instability theory was first developed to explore this possibility. A series of early papers by Rayleigh produced many notable results concerning the instability of inviscid flows, such as the discovery of inflectional instability. Viscosity was commonly thought to act only to stabilize the flow, and flows with convex velocity profiles appeared to be stable. The investigations that led to a viscous theory of boundary layer instability was reported. The earliest application of linear stability theory to transition prediction calculated the amplitude ratio of the most amplified frequency as a function of Reynolds number for a Blasius boundary layer, and found that this quantity had values between five and nine at the observed Ret. The experiment of Schubauer and Skramstad (1947) completely reversed the prevailing option and fully vindicated the Gottingen proponents of the theory. This experiment demonstrated the existence of instability waves in a boundary layer, their connection with transition, and the quantitative description of their behavior by the theory of Tollmien and Schlichting. It is generally accepted that flow parameters such as pressure gradient, suction and heat transfer qualitatively affect transition in the manner predicted by the linear theory, and in particular that a flow predicted to be stable by the theory should remain laminar. The linear theory, in the form of the e9, or N-factor is today in routine use in engineering studies of laminar flow. The stability theory to boundary layers with pressure gradients and suction was applied. The only large body of numerical results for exact boundary layer solutions before the advent of the computer age by calculating the stability characteristics of the Falkner-Skan family of velocity profiles are given. When the digital computer reached a stage of development which permit the direct solution of the primary differential equations, numerical results were obtained from the linear theory during the next 10 years for many different boundary layer flows: three dimensional boundary layers; free convention boundary layers; compressible boundary layers; boundary layers on compliant walls; a recomputation of Falkner-Skan flows; unsteady boundary layers; and heated wall boundary layers.

Leslie M Mack↗
NASA NTRS1972

Equations of motion for the X-14 Aircraft, phase 2 study

A study of the control and power requirements of the X-14 VTOL aircraft is presented. The complete equations of motion for X-14 are derived. The fundamental assumption is that the aircraft is a single rigid body. The equations of motion are derived with respect to a set of axes fixed to the aircraft. Additional assumptions used are that any wind disturbances are irrotational that the twin engines used on the aircraft rotate in the same direction at the same speed and that the engine exhaust is diverted by means of vanes to provide a direction varying thrust vector. The equations obtained are subsequently linearized about various reference conditions and numerical values for the trim parameters and the stability derivatives at these conditions are tabulated.

Hoffman, M. A.↗
NASA NTRS1974

Research study on stabilization and control modern sampled-data control theory

The methods of continuous and discrete describing function analysis were applied to predicting the existence of self-sustained oscillations in the single-axis model of the large space telescope system with nonlinear control moment gyroscope friction characteristics. It is shown that the stability equations may be solved by a numerical-iterative technique using the describing function analysis, instead of the usual graphical methods. The numerical method is found to be effective in leading to a convergent solution rapidly, with an appropriate guess of the initial condition.

Kuo, B. C.↗
NASA NTRS1975

The Cepheid mode problem

The modal behavior of RR Lyrae stars and low mass Cepheids has been investigated using a nonlinear relaxation technique and stability analysis. The advantages of this type of numerical approach in investigations of preferred mode of pulsation are discussed. The results obtained for both classes of variable stars are quite similar: first harmonic pulsation toward the blue, fundamental pulsation toward the red, and mixed-mode behavior at the extreme red edge of the instability strip. In addition, stars near the center of the strip can pulsate in either the fundamental or the first harmonic mode. Possible implications for observational results, including the Oosterhoff dichotomy of globular clusters and the beat Cepheids, are discussed.

Stellingwerf, R. F.↗
NASA NTRS1975

Measurement updating using the U-D factorization

A new mechanization of the Kalman updating algorithm based on a U-D factorization of the estimate error covariance is introduced. Efficient and stable updating recursions are developed for the unit upper triangular factor U and the diagonal factor D, treating only the parameter estimation problem. Properties of the factorization update performed here include efficient one point at a time processing that requires little more computation than the optimal but numerically unstable conventional Kalman measurement update algorithm, and stability that compares with the square root filter.

Bierman, G. J.↗
NASA NTRS1977

Remote sounding of cloudy atmospheres. II - Multiple cloud formations

The dual frequency-range principle developed in Part I (Chahine, 1974) for infrared remote sounding of atmospheric temperature profiles in the presence of a single cloud layer is extended here to the case of multiple cloud formations. The approach requires no a priori knowledge of the spectral properties of the clouds or the number of cloud layers in the fields of view. The method of solution requires measurements over adjacent fields of view and leads to the determination of the clear-column atmospheric temperature profiles with the same degree of accuracy and vertical resolution permitted under cloudless conditions. Numerical verifications are carried out to illustrate the stability and accuracy of the method using simulated radiance data, from the 4.3 and 15 micron CO2 bands in the terrestrial atmosphere, in the presence of up to three cloud layers in the fields of view.

Chahine, M. T.↗
NASA NTRS1983

On the solution of elastic-plastic static and dynamic postbuckling collapse of general structure

Many investigations have considered structural collapse from strictly the transient point of view. While such an approach is ideally correct, certain difficulties have to be overcome in its implementation. The present investigation is concerned with the development of self-adaptive algorithms which make it possible to conduct the analysis of both static elastic and elastic-plastic postbuckling, as well as static loading to the onset of buckling followed by subsequent dynamic postbuckling. The approach employed to solve the static portion of loading is to extend the constrained Incremental Newton-Raphson (INR) algorithm by incorporating elastic-plastic constitutive characterizations. Large deformation moderate strain theory is adopted to establish the overall strategy. Attention is given to governing field equations, aspects of algorithmic development, and numerical experiments conducted to illustrate the efficiency and stability of the developed schemes.

Padovan, J.↗
NASA NTRS1983

Linear density drift instabilities in very low beta plasmas A different approach

This paper reports a study of the linear Vlasov electromagnetic dispersion equation for density drift instabilities in very low beta plasmas. A uniform magnetic field and a weak uniform density gradient are assumed. This paper differs from most other studies of this topic in three important ways: (1) no low frequency or long wavelength approximations are made, (2) no gauge condition is imposed and (3) a modification of the local approximation is used which is argued to be less arbitrary than the usual local approximation. Numerical solution of the resulting dispersion equation yields stabilization of the universal density drift instability at values of beta somewhat higher than those obtained from the local approximation, and near maximum growth shows no evidence of coupling between the universal and Alfven modes.

Gary, S. P.↗
NASA NTRS1986

Two time scale output feedback regulation for ill-conditioned systems

Issues pertaining to the well-posedness of a two time scale approach to the output feedback regulator design problem are examined. An approximate quadratic performance index which reflects a two time scale decomposition of the system dynamics is developed. It is shown that, under mild assumptions, minimization of this cost leads to feedback gains providing a second-order approximation of optimal full system performance. A simplified approach to two time scale feedback design is also developed, in which gains are separately calculated to stabilize the slow and fast subsystem models. By exploiting the notion of combined control and observation spillover suppression, conditions are derived assuring that these gains will stabilize the full-order system. A sequential numerical algorithm is described which obtains output feedback gains minimizing a broad class of performance indices, including the standard LQ case. It is shown that the algorithm converges to a local minimum under nonrestrictive assumptions. This procedure is adapted to and demonstrated for the two time scale design formulations.

Calise, A. J.↗
NASA NTRS1987

Decoupled direct method for sensitivity analysis in combustion kinetics

An efficient, decoupled direct method for calculating the first order sensitivity coefficients of homogeneous, batch combustion kinetic rate equations is presented. In this method the ordinary differential equations for the sensitivity coefficients are solved separately from , but sequentially with, those describing the combustion chemistry. The ordinary differential equations for the thermochemical variables are solved using an efficient, implicit method (LSODE) that automatically selects the steplength and order for each solution step. The solution procedure for the sensitivity coefficients maintains accuracy and stability by using exactly the same steplengths and numerical approximations. The method computes sensitivity coefficients with respect to any combination of the initial values of the thermochemical variables and the three rate constant parameters for the chemical reactions. The method is illustrated by application to several simple problems and, where possible, comparisons are made with exact solutions and those obtained by other techniques.

Radhakrishnan, Krishnan↗
NASA NTRS1986

Stability with large step sizes for multistep discretizations of stiff ordinary differential equations

One-leg and multistep discretizations of variable-coefficient linear systems of ODEs having both slow and fast time scales are investigated analytically. The stability properties of these discretizations are obtained independent of ODE stiffness and compared. The results of numerical computations are presented in tables, and it is shown that for large step sizes the stability of one-leg methods is better than that of the corresponding linear multistep methods.

Majda, George↗
NASA NTRS1989

Time domain numerical calculations of unsteady vortical flows about a flat plate airfoil

A time domain numerical scheme is developed to solve for the unsteady flow about a flat plate airfoil due to imposed upstream, small amplitude, transverse velocity perturbations. The governing equation for the resulting unsteady potential is a homogeneous, constant coefficient, convective wave equation. Accurate solution of the problem requires the development of approximate boundary conditions which correctly model the physics of the unsteady flow in the far field. A uniformly valid far field boundary condition is developed, and numerical results are presented using this condition. The stability of the scheme is discussed, and the stability restriction for the scheme is established as a function of the Mach number. Finally, comparisons are made with the frequency domain calculation by Scott and Atassi, and the relative strengths and weaknesses of each approach are assessed.

Hariharan, S. I.↗
NASA NTRS1991

Application of the boundary element method to transient heat conduction

An advanced boundary element method (BEM) is presented for the transient heat conduction analysis of engineering components. The numerical implementation necessarily includes higher-order conforming elements, self-adaptive integration and a multiregion capability. Planar, three-dimensional and axisymmetric analyses are all addressed with a consistent time-domain convolution approach, which completely eliminates the need for volume discretization for most practical analyses. The resulting general purpose algorithm establishes BEM as an attractive alternative to the more familiar finite difference and finite element methods for this class of problems. Several detailed numerical examples are included to emphasize the accuracy, stability and generality of the present BEM. Furthermore, a new efficient treatment is introduced for bodies with embedded holes. This development provides a powerful analytical tool for transient solutions of components, such as casting moulds and turbine blades, which are cumbersome to model when employing the conventional domain-based methods.

Dargush, G. F.↗
NASA NTRS1991

Numerical simulation of shock-induced combustion/detonation in a premixed H2-air mixture using Navier-Stokes equations

A numerical study was conducted to address the structural stability of the oblique detonation wave. An attempt was made to isolate the structural instability predicted by Buckmaster's (1990) linear stability analysis. For this study a detailed viscous flowfield with finite-rate chemistry was computed past a 20-deg wedge. The finite-rate chemistry was modeled by using a seven-species and seven-reaction model. Instability modes were isolated by analyzing the Fourier power spectrum of the H2O mass fraction at selected sample stations. No unstable temporal modes with dominant frequency were isolated. It was concluded that the instability which is predicted by the Buckmaster work was a high-frequency low-amplitude phenomenon. For engineering purposes, the oblique detonation wave is a stable phenomenon as long as sufficient levels of overdrive are present.

Singh, D. J.↗