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 271 records · Page 15

Research study on stabilization and control: Modern sampled data control theory

A numerical analysis of spacecraft stability parameters was conducted. The analysis is based on a digital approximation by point by point state comparison. The technique used is that of approximating a continuous data system by a sampled data model by comparison of the states of the two systems. Application of the method to the digital redesign of the simplified one axis dynamics of the Skylab is presented.

Kuo, B. C.↗

Stability of polar rings around Neptune

An analytical and numerical study of the stability of polar rings around Neptune is presented. The stability proofs are based on various methods used to study gas disks in galaxies. It is shown that stable polar rings can exist despite energy dissipation by collisions between particles. Also, four equilibrium orientations which pass nearly over the pole of Neptune are found, two of which are stable in the presence of dissipation.

Dobrovolskis, Anthony R.↗

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from a source surface, given the proper boundary conditions. The results of this work; describing the mathematical formalism that was developed, the numerical techniques employed, and the stability criteria developed for these numerical schemes are presented. An analytical solution is used for a test case; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent ( 5%).

Wu, S. T.↗

Numerical Boundary Condition Procedures

Topics include numerical procedures for treating inflow and outflow boundaries, steady and unsteady discontinuous surfaces, far field boundaries, and multiblock grids. In addition, the effects of numerical boundary approximations on stability, accuracy, and convergence rate of the numerical solution are discussed.

Source record↗

Error propagation in the numerical solutions of the differential equations of orbital mechanics

The relationship between the eigenvalues of the linearized differential equations of orbital mechanics and the stability characteristics of numerical methods is presented. It is shown that the Cowell, Encke, and Encke formulation with an independent variable related to the eccentric anomaly all have a real positive eigenvalue when linearized about the initial conditions. The real positive eigenvalue causes an amplification of the error of the solution when used in conjunction with a numerical integration method. In contrast an element formulation has zero eigenvalues and is numerically stable.

Bond, V. R.↗

Influence of exit impedance on finite difference solutions of transient acoustic mode propagation in ducts

The time-dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard-wall duct without flow and with spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the steady solutions of the acoustic field. The transient method did not converge to the steady-state solution for cutoff acoustic duct modes. This has implications as to its use in a variable-area duct, where modes may become cutoff in the smal-area portion of the duct. For single cutoff mode propagation the steady-state impedance boundary condition produced acoustic reflections during the initial transient that caused finite instabilities in the numerical calculations. The stability problem is resolved by reformulating the exit boundary condition. Example calculations show good agreement with exact analytical and numerical results for forcing frequencies above, below, and nearly at the cutoff frequency.

Baumeister, K. J.↗

Energy management during the space shuttle transition.

An approach to calculating optimal, gliding flight paths of the type associated with the space shuttle's transition from entry to cruising flight is presented. Kinetic energy and total energy (per unit weight) replace velocity and time in the dynamic equations, reducing the dimension and complexity of the problem. The capability for treating integral and terminal penalties (as well as Mach number effects) is retained in the numerical optimization; hence, stability and control boundaries can be observed as trajectories to the desired final energy, flight path angle, and range are determined. Numerical results show that the 'jump' to the 'front-side of the L/D curve' need not be made until the end of the transition and that the dynamic model provides a conservative range estimate. Alternatives for real-time trajectory control are discussed.

Stengel, R. F.↗

Energy management during the space shuttle transition

An approach to calculating optimal, gliding flight paths of the type associated with the space shuttle's transition from entry to cruising flight is presented. Kinetic energy and total energy (per unit weight) replace velocity and time in the dynamic equations, reducing the dimension and complexity of the problem. The capability for treating integral and terminal penalties (as well as Mach number effects) is retained in the numerical optimization; hence, stability and control boundaries can be observed as trajectories to the desired final energy, flight path angle, and range are determined. Numerical results show that the jump to the front-side of the L/D curve need not be made until the end of the transition and that the dynamic model provides a conservative range estimate. Alternatives for real time trajectory control are discussed.

Stengel, R. F.↗

Numerical Simulation And Experimental Investigation Of The Lift-Off And Blowout Of Enclosed Laminar Flames

Enclosed flames are found in many industrial applications such as power plants, gas-turbine combustors and jet engine afterburners. A better understanding of the burner stability limits can lead to development of combustion systems that extend the lean and rich limits of combustor operations. This paper reports a fundamental study of the stability limits of co-flow laminar jet diffusion flames. A numerical study was conducted that used an adaptive mesh refinement scheme in the calculation. Experiments were conducted in two test rigs with two different fuels and diluted with three inert species. The numerical stability limits were compared with microgravity experimental data. Additional normal-gravity experimental results were also presented.

Venuturmilli, Rajasekhar↗

Two-dimensional stability of laminar flames

The stability of laminar flames was studied numerically and the dependence of stability on finite rate chemistry with low activation energy and variable thermodynamic and transport properties was addressed. The calculations show that activation energy and details of chemistry play a minor role in altering the linear neutral stability results from asymptotic analysis. Variable specific heat makes a marginal change to the stability; variable transport properties, on the other hand, tend to substantially enhance the stability from a critical wave number of about 0.50 to 0.20. Also, the effects of variable properties tend to nullify the effects of nonunity Lewis number. When the Lewis number of a single species is different from unity, as is true in a hydrogen-air premixed flame, the stability results remain close to that of unity Lewis number.

Mukunda, H. S.↗

Partially Grooved Domain Stabilization Structures for Vertical Bloch Line Memory

Bias field stability ranges were measured and numerically simulated for magnetic domains ingarnets stabilized by partially grooved rectangular and ring grooves. Simulation results agreefavorably with experimental results when finite slope effects of the groove walls are included. Asbias fields increase, rectangular and ring domains both destabilize through stripe head recession. Asbias fields decrease, destabilization in rectangular domains occurs by runout, while destabilization inring domains occurs by midstripe domain buckling. While ring domains are stable at lower biasfields than rectangular domains, bias field stability ranges are approximately equal. Hence, for thesame partial grooving depth, rectangular domains are preferred because they offer higher storagedensity potential in Vertical Block Line (VBL) storage arrays as long as bit propagation margins atstripe ends are sufficient.

Katti, R. R.↗

Construction of Low Dissipative High Order Well-Balanced Filter Schemes for Non-Equilibrium Flows

The goal of this paper is to generalize the well-balanced approach for non-equilibrium flow studied by Wang et al. [26] to a class of low dissipative high order shock-capturing filter schemes and to explore more advantages of well-balanced schemes in reacting flows. The class of filter schemes developed by Yee et al. [30], Sjoegreen & Yee [24] and Yee & Sjoegreen [35] consist of two steps, a full time step of spatially high order non-dissipative base scheme and an adaptive nonlinear filter containing shock-capturing dissipation. A good property of the filter scheme is that the base scheme and the filter are stand alone modules in designing. Therefore, the idea of designing a well-balanced filter scheme is straightforward, i.e., choosing a well-balanced base scheme with a well-balanced filter (both with high order). A typical class of these schemes shown in this paper is the high order central difference schemes/predictor-corrector (PC) schemes with a high order well-balanced WENO filter. The new filter scheme with the well-balanced property will gather the features of both filter methods and well-balanced properties: it can preserve certain steady state solutions exactly; it is able to capture small perturbations, e.g., turbulence fluctuations; it adaptively controls numerical dissipation. Thus it shows high accuracy, efficiency and stability in shock/turbulence interactions. Numerical examples containing 1D and 2D smooth problems, 1D stationary contact discontinuity problem and 1D turbulence/shock interactions are included to verify the improved accuracy, in addition to the well-balanced behavior.

Wang, Wei↗

Langley Stability and Transition Analysis Code (LASTRAC) Version 1.2 User Manual

LASTRAC is a general-purposed, physics-based transition prediction code released by NASA for Laminar Flow Control studies and transition research. The design and development of the LASTRAC code is aimed at providing an engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. It was written from scratch based on the state-of-the-art numerical methods for stability analysis and modern software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory or linear parabolized stability equations method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. This document describes the governing equations, numerical methods, code development, detailed description of input/output parameters, and case studies for the current release of LASTRAC.

Chang, Chau-Lyan↗

Stability Limits and Dynamics of Nonaxisymmetric Liquid Bridges

Theoretical and experimental investigation of the stability of nonaxisymmetric and axisymmetric bridges contained between equal and unequal radii disks as a function of Bond and Weber number with emphasis on the transition from unstable axisymmetric to stable nonaxisymmetric shapes. Numerical analysis of the stability of nonaxisymmetric bridges between unequal disks for various orientations of the gravity vector Experimental and numerical investigation of bridge stability (nonaxisymmetric and axisymmetric), large amplitude (nonaxisymmetric) oscillations and breaking.

Alexander, J. Iwan D.↗

The Langley Stability and Transition Analysis Code (LASTRAC) : LST, Linear and Nonlinear PSE for 2-D, Axisymmetric, and Infinite Swept Wing Boundary Layers

During the past two decades, our understanding of laminar-turbulent transition flow physics has advanced significantly owing to, in a large part, the NASA program support such as the National Aerospace Plane (NASP), High-speed Civil Transport (HSCT), and Advanced Subsonic Technology (AST). Experimental, theoretical, as well as computational efforts on various issues such as receptivity and linear and nonlinear evolution of instability waves take part in broadening our knowledge base for this intricate flow phenomenon. Despite all these advances, transition prediction remains a nontrivial task for engineers due to the lack of a widely available, robust, and efficient prediction tool. The design and development of the LASTRAC code is aimed at providing one such engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. LASTRAC was written from scratch based on the state-of-the-art numerical methods for stability analysis and modem software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory (LST) or linear parabolized stability equations (LPSE) method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. Coupled with the built-in receptivity model that is currently under development, the nonlinear PSE method offers a synergistic approach to predict transition onset for a given disturbance environment based on first principles. This paper describes the governing equations, numerical methods, code development, and case studies for the current release of LASTRAC. Practical applications of LASTRAC are demonstrated for linear stability calculations, N-factor transition correlation, non-linear breakdown simulations, and controls of stationary crossflow instability in supersonic swept wing boundary layers.

Chang, Chau-Lyan↗

Remarks on disputed numerical results in compressible boundary-layer stability theory

Past work on the influence of Mach number on the viscous and inviscid instability of flat-plate boundary layers is reviewed, and new spatial calculations are presented. These calculations support the previous view that viscosity is only stabilizing for both two- and three-dimensional first-mode waves above M1 = 3.0, and for second-mode waves at all Mach numbers. It is concluded that the calculations of Wazzan, Taghavi, and Keltner that show viscous instability at M1 = 6.0 for first-mode 50 deg waves, and at M1 = 3.0 for two-dimensional second-mode waves, are not correct.

Mack, L. M.↗