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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Asymptotic behavior of modulated Taylor-Couette flows with a crystalline inner cylinder

The linear stability of a modulated Taylor-Couette system when the inner cylindrical boundary consists of a crystalline solid-liquid interface is considered. Both experimentally and in numerical calculations it is found that the two-phase system is significantly less stable than the analogous rigid-walled system for materials with moderately large Prandtl numbers. A numerical treatment based on Floquet theory is described, which gives results that are in good agreement with preliminary experimental findings. In addition, this instability is further examined by carrying out a formal asymptotic expansion of the solution in the limit of large Prandtl number. In this limit the Floquet analysis is considerably simplified, and the linear stability of the modulated system can be determined to leading order through a conventional stability analysis, without recourse to Floquet theory. The resulting simplified problem is then studied for both the narrow gap geometry and for the case of a finite gap. It is surprising that the determination of the linear stability of the two-phase system is considerably simpler than that of the rigid-walled system, despite the complications introduced by the presence of the crystal-melt interface.

Braun, R. J.↗

NASA NTRS1997

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗

NASA NTRS1995

Finite element buckling and postbuckling solutions for multilayered composite panels

A study is made of the buckling and postbuckling responses of flat, unstiffened composite panels subjected to various combinations of mechanical and thermal loads. The analysis is based on a first-order shear deformation von Karman-type plate theory. A mixed formulation is used with the fundamental unknowns consisting of the strain components, stress resultants and the generalized displacements of the plate. The stability boundary, postbuckling response and the sensitivity coefficients are evaluated. The sensitivity coefficients measure the sensitivity of the buckling and postbuckling responses to variations in the different lamination and material parameters of the panel. Numerical results are presented for both solid panels and panels with central circular cutouts. The results show the effects of the variations in the fiber orientation angels, aspect ratio of the panel, and the hole diameter (for panels with cutouts) on the stability boundary, postbuckling response and sensitivity coefficients.

Noor, Ahmed K.↗

NASA NTRS1993

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗

NASA NTRS1983

Viscous, resistive magnetohydrodynamic stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one-dimensional magnetohydrodynamic quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds-like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds-like numbers. Marginal stability curves, growth rates versus Reynolds-like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result that appears general is that instability has been found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three-dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

NASA NTRS1991

The L sub 1 finite element method for pure convection problems

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Jiang, Bo-Nan↗

NASA NTRS1990

Numerical simulation of transition in a rotating disk flow

An investigation of the laminar-turbulence transition process in a three-dimensional boundary layer on a rotating disk is initiated. Preliminary results from a direct numerical simulation using spectral collocation method are presented. The simulation results in the linear regime is tested against linear stability theory. Saturation of the primary disturbance energy suggests the importance of modal interaction in initiating transition.

Balachandar, S.↗

NASA NTRS1991

Optimal discrete-time dynamic output-feedback design - A w-domain approach

An alternative method for optimal digital control design is described in this paper. The method is based on the usage of the w-transform and has many attractive design features. One of these is its immediate connection with frequency loop-shaping techniques that are now popular and effective for multivariable control synthesis in continuous-time domain. Furthermore, any design algorithms originally developed for continuous-time systems can now be immediately extended to the discrete-time domain. The main results presented in this paper are the exact problem formulation and solution of an optimal discrete-time dynamic output-feedback design in the w-domain involving a quadratic performance index to random disturbances. In addition, necessary conditions for optimality are obtained for the numerical solution of the optimal output-feedback compensator design. A numerical example is presented illustrating its application to the design of a low-order dynamic compensator in a stability augmentation system of a commercial transport.

Ha, Cheolkeun↗

NASA NTRS1991

The L sub 1 finite element method for pure convection problems

Jiang, Bo-Nan↗

NASA NTRS1995

The mesoscale forcing of a midlatitude upper-tropospheric jet streak by a simulated convective system. 1: Mass circulation and ageostrophic processes

The mutual forcing of a midlatitude upper-tropospheric jet streak by organized mesoscale adiabatic and diabatic processes within a simulated convective system (SCS) is investigated. Using isentropic diagnostics, results from a three-dimensional numerical simulation of an SCS are examined to study the isallobaric flow field, modes of dominant ageostrophic motion, and stability changes in relation to the mutual interdependence of adiabatic processes and latent heat release. Isentropic analysis affords an explicit isolation of a component of isallobaric flow associated with diabatic processes within the SCS. Prior to convective development within the simulations, atmospheric destabilization occurs through adiabatic ageostrophic mass adjustment and low-level convergence in association with the preexisting synoptic-scale upper-tropospheric jet streak. The SCS develops in a baroclinic zone and quickly initiates a vigorous mass circulation. By the mature stage, a pronounced vertical couplet of low-level convergence and upper-level mass divergence is established, linked by intense midtropospoheric diabatic heating. Significant divergence persists aloft for several hours subsequent to SCS decay. The dominant role of ageostrophic motion within which the low-level mass convergence develops is the adiabatic isallobaric component, while the mass divergence aloft develops principally through the diabatic isallobaric component. Both compnents are intrinsically linked to the convectively forced vertical mass transport. The inertial diabatic ageostrophic component is largest near the level of maximum heating and is responsible for the development of inertial instability to the north of SCS, resulting in this quadrant being preferred for outflow. The inertial advective component, the dominant term that produces the new downstream wind maximum, rapidly develops north of the SCS and through mutual adjustment creates the baroclinic support for the new jet streak.

Wolf, Bart J.↗

NASA NTRS2001

Experimental and Computational Study of Trapped Vortex Combustor Sector Rig with Tri-pass Diffuser

The Trapped Vortex Combustor (TVC) potentially offers numerous operational advantages over current production gas turbine engine combustors. These include lower weight, lower pollutant emissions, effective flame stabilization, high combustion efficiency, excellent high altitude relight capability, and operation in the lean burn or RQL (Rich burn/Quick mix/Lean burn) modes of combustion. The present work describes the operational principles of the TVC, and provides detailed performance data on a configuration featuring a tri-pass diffusion system. Performance data include EINOx (NO(sub x) emission index) results for various fuel-air ratios and combustor residence times, combustion efficiency as a function of combustor residence time, and combustor lean blow-out (LBO) performance. Computational fluid dynamics (CFD) simulations using liquid spray droplet evaporation and combustion modeling are performed and related to flow structures observed in photographs of the combustor. The CFD results are used to understand the aerodynamics and combustion features under different fueling conditions. Performance data acquired to date are favorable in comparison to conventional gas turbine combustors. Further testing over a wider range of fuel-air ratios, fuel flow splits, and pressure ratios is in progress to explore the TVC performance. In addition, alternate configurations for the upstream pressure feed, including bi-pass diffusion schemes, as well as variations on the fuel injection patterns, are currently in test and evaluation phases.

Hendricks, Robert C.↗

NASA NTRS2004

Experimental and Computational Study of Trapped Vortex Combustor Sector Rig With Tri-Pass Diffuser

The Trapped Vortex Combustor (TVC) potentially offers numerous operational advantages over current production gas turbine engine combustors. These include lower weight, lower pollutant emissions, effective flame stabilization, high combustion efficiency, excellent high altitude relight capability, and operation in the lean burn or RQL modes of combustion. The present work describes the operational principles of the TVC, and extends diffuser velocities toward choked flow and provides system performance data. Performance data include EINOx results for various fuel-air ratios and combustor residence times, combustion efficiency as a function of combustor residence time, and combustor lean blow-out (LBO) performance. Computational fluid dynamics (CFD) simulations using liquid spray droplet evaporation and combustion modeling are performed and related to flow structures observed in photographs of the combustor. The CFD results are used to understand the aerodynamics and combustion features under different fueling conditions. Performance data acquired to date are favorable compared to conventional gas turbine combustors. Further testing over a wider range of fuel-air ratios, fuel flow splits, and pressure ratios is in progress to explore the TVC performance. In addition, alternate configurations for the upstream pressure feed, including bi-pass diffusion schemes, as well as variations on the fuel injection patterns, are currently in test and evaluation phases.

Hendricks, R. C.↗

NASA NTRS2008

Ares-I Bending Filter Design using a Constrained Optimization Approach

The Ares-I launch vehicle represents a challenging flex-body structural environment for control system design. Software filtering of the inertial sensor output is required to ensure adequate stable response to guidance commands while minimizing trajectory deviations. This paper presents a design methodology employing numerical optimization to develop the Ares-I bending filters. The design objectives include attitude tracking accuracy and robust stability with respect to rigid body dynamics, propellant slosh, and flex. Under the assumption that the Ares-I time-varying dynamics and control system can be frozen over a short period of time, the bending filters are designed to stabilize all the selected frozen-time launch control systems in the presence of parameter uncertainty. To ensure adequate response to guidance command, step response specifications are introduced as constraints in the optimization problem. Imposing these constrains minimizes performance degradation caused by the addition of the bending filters. The first stage bending filter design achieves stability by adding lag to the first structural frequency to phase stabilize the first flex mode while gain stabilizing the higher modes. The upper stage bending filter design gain stabilizes all the flex bending modes. The bending filter designs provided here have been demonstrated to provide stable first and second stage control systems in both Draper Ares Stability Analysis Tool (ASAT) and the MSFC MAVERIC 6DOF nonlinear time domain simulation.

Hall, Charles↗

NASA NTRS2010

Ares I Flight Control System Design

The Ares I launch vehicle represents a challenging flex-body structural environment for flight control system design. This paper presents a design methodology for employing numerical optimization to develop the Ares I flight control system. The design objectives include attitude tracking accuracy and robust stability with respect to rigid body dynamics, propellant slosh, and flex. Under the assumption that the Ares I time-varying dynamics and control system can be frozen over a short period of time, the flight controllers are designed to stabilize all selected frozen-time launch control systems in the presence of parametric uncertainty. Flex filters in the flight control system are designed to minimize the flex components in the error signals before they are sent to the attitude controller. To ensure adequate response to guidance command, step response specifications are introduced as constraints in the optimization problem. Imposing these constraints minimizes performance degradation caused by the addition of the flex filters. The first stage bending filter design achieves stability by adding lag to the first structural frequency to phase stabilize the first flex mode while gain stabilizing the higher modes. The upper stage bending filter design gain stabilizes all the flex bending modes. The flight control system designs provided here have been demonstrated to provide stable first and second stage control systems in both Draper Ares Stability Analysis Tool (ASAT) and the MSFC 6DOF nonlinear time domain simulation.

Jang, Jiann-Woei↗

NASA NTRS2019

Application of System Identification to Parachute Modeling

Parachute models are used in numerous flight simulation tools to predict a wide range of parachute flight performance characteristics (e.g., parachute inflation loads, parachute stability and dynamics, vehicle touchdown conditions, and, ultimately, the safety and survivability of the system using the parachute). The current state of the art in developing parachute models is to initially estimate the parachute characteristics based on the parachute geometry and historical data and then add increased model fidelity based on data from wind tunnel and/or flight tests. This approach, however, can be deficient in identifying which parachute states (e.g., angle of attack, sideslip, angular rates, flyout angles, descent rate, dynamic pressure, proximity to other parachutes) are responsible for the parachute motion, and the relationship between those states and the forces on the parachute.

Murri, Daniel G.↗

NASA NTRS2001

The Behavior of Methane-Air Partially Premixed Flames Under Normal- and Zero-G Conditions

Partially premixed flames (PPFs) represent a class of hybrid flames containing multiple reaction zones. These flames are established when less than stoichiometric quantity of oxidizer is molecularly mixed with the fuel stream before entering the reaction zone where additional oxidizer is available for complete combustion. This mode of combustion can be used to exploit the advantages of both nonpremixed and premixed flames regarding operational safety, lower pollutant emissions and flame stabilization. A double flame containing a fuel-rich premixed reaction zone, which is anchored by a nonpremixed reaction zone, is one example of a partially premixed flame. A triple flame is also a PPF that contains three reaction zones, namely, a fuel-rich premixed zone, a fuel-lean premixed zone, and a nonpremixed reaction zone. Herein we focus on two aspects of our investigation, one involving the development of optical diagnostics that can be used on a microgravity rig, which has been recently fabricated, and the other on the numerically predicted differences between normal- and zero-gravity PPFs. Both the measurements and simulations examine the detailed structure of methane-air PPFs stabilized on a Wolfhard-Parker slot burner.

Puri, Ishwar K.↗

NASA NTRS1994

Analytical investigation of the dynamics of tethered constellations in Earth orbit, phase 2

This final report covers nine years of research on future tether applications and on the actual flights of the Small Expendable Deployment System (SEDS). Topics covered include: (1) a description of numerical codes used to simulate the orbital and attitude dynamics of tethered systems during station keeping and deployment maneuvers; (2) a comparison of various tethered system simulators; (3) dynamics analysis, conceptual design, potential applications and propagation of disturbances and isolation from noise of a variable gravity/microgravity laboratory tethered to the Space Station; (4) stability of a tethered space centrifuge; (5) various proposed two-dimensional tethered structures for low Earth orbit for use as planar array antennas; (6) tethered high gain antennas; (7) numerical calculation of the electromagnetic wave field on the Earth's surface on an electrodynamically tethered satellite; (8) reentry of tethered capsules; (9) deployment dynamics of SEDS-1; (10) analysis of SEDS-1 flight data; and (11) dynamics and control of SEDS-2.

Lorenzini, Enrico C.↗

NASA NTRS1986

Numerical and physical aspects of aerodynamic flows III; Proceedings of the Third Symposium, California State University, Long Beach, CA, January 21-24, 1985

The present symposium considers hydrodynamic stability and turbulent transition, transition calculations in three-dimensional flows, a quasi-simultaneous finite difference approach for strongly interacting flows, the significance of the thin layer Navier-Stokes approximation, unsteady airfoil boundary layers, predictions and experiments on airfoils at low Reynolds numbers, and a comparison of interactive boundary layer and thin layer Navier-Stokes procedures. Also discussed are a viscous-inviscid interaction method for computing unsteady transonic separation, massive separation and dynamic stall on a cusped trailing edge airfoil, the computation of turbulent separated flows over wings, an iterative scheme for three-dimensional transonic flows, and the computation of three-dimensional flows with shock wave/boundary layer interaction. Many of the papers in this symposium were abstracted previously (cf., A85-42951).

Cebeci, T.↗