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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

An unsteady lifting surface method for single rotation propellers

The mathematical formulation of a lifting surface method for evaluating the steady and unsteady loads induced on single rotation propellers by blade vibration and inflow distortion is described. The scheme is based on 3-D linearized compressible aerodynamics and presumes that all disturbances are simple harmonic in time. This approximation leads to a direct linear integral relation between the normal velocity on the blade (which is determined from the blade geometry and motion) and the distribution of pressure difference across the blade. This linear relation is discretized by breaking the blade up into subareas (panels) on which the pressure difference is treated as approximately constant, and constraining the normal velocity at one (control) point on each panel. The piece-wise constant loads can then be determined by Gaussian elimination. The resulting blade loads can be used in performance, stability and forced response predictions for the rotor. Mathematical and numerical aspects of the method are examined. A selection of results obtained from the method is presented. The appendices include various details of the derivation that were felt to be secondary to the main development in Section 1.

Williams, Marc H.↗
NASA NTRS1990

Analytical model for orbital debris environmental management

A differential equation, also referred to as the PIB (particle-in-a-box) model, expressing the time rate of change of the number of objects in orbit, is developed, and its applicability is illustrated. The model can be used as a tool for the assessment of LEO environment stability, and as a starting point for the development of numerical evolutionary models. Within the context of the model, evolutionary scenarios are examined, and found to be sensitive to the growth rate. It is determined that the present environment is slightly unstable to catastrophic growth, and that the number of particles on orbit will continue to increase until approximately 2250-2350 AD, with a maximum of 2,000,000. The model is expandable to the more realistic (complex) case of multiple species in a multiple-tier system.

Talent, David L.↗
NASA NTRS1991

Centralized, decentralized, and independent control of a flexible manipulator on a flexible base

The dynamics and control of a flexible manipulator arm with payload mass on a flexible base in space are considered. The controllers are provided by one torquer at the center of the base and one torquer at the connection joint of the robot and the base. The nonlinear dynamics of the system is modeled by applying the finite element method and Lagrangian formula. Three control strategies are considered and compared, i.e., centralized control, decentralized control, and independent control. All these control designs are based on the linear quadratic regulator theory. A mathematical decomposition is used in the decentralization process so that the coupling between the subsystems is weak, while a physical decomposition is used in the independent control design process. For both the decentralized and the independent controls, the stability of the overall linear system is checked before a numerical simulations is initiated. Two numerical examples show that the response of the independent control system are close to those of the centralized control system, while the responses of the decentralized control system are not.

Li, Feiyue↗
NASA NTRS1992

The influence of bulges on boundary-layer instability

Local disturbances caused by a spanwise surface corrugation affect the position of the boundary-layer transition, and so the drag, of an object. This premature transition from laminar to turbulent flow is often associated with a separation of the laminar boundary-layer from its surface. Also the roughness-induced separation bubble provides an important link between the pressure and velocity fluctuations in the environment and the development of the disturbance in the laminar boundary-layer, i.e., the receptivity problem. To investigate the influence of a laminar separation bubble on boundary-layer instability, a separated flow generated by a velocity gradient over a flat plate was analyzed by direct numerical simulation using finite-difference solutions of the Navier-Stokes equations. The bubble acts as a strong amplifier of the instability waves and a highly nonlinear flow field is shown to develop downstream of the bubble. Consequently, the results of the direct numerical simulation differ noticeably from those of the classical linear stability theory proving the fact that the nonparallel effects together with the nonlinear interactions are crucial to this flow development. In the present paper, the effect of physical perturbations such as humps and hollows on boundary-layer instability is analyzed. This problem has been considered theoretically by several researchers (e.g., Nayfeh et al., 1987 and 1990; Cebeci et al., 1988). They used linear stability theory in their approach which does not include the nonparallel nor the nonlinear effects. Therefore, to account for these important effects in studying flow over humps and hollows the direct simulation technique is being implemented in generalized coordinates.

Elli, S.↗
NASA NTRS1992

On the structure of cellular solutions in Rayleigh-Benard-Marangoni flows in small-aspect-ratio containers

Multiple steady flow patterns occur in surface-tension/buoyancy-driven convection in a liquid layer heated from below (Rayleigh-Benard-Marangoni flows). Techniques of numerical bifurcation theory are used to study the multiplicity and stability of two-dimensional steady flow patterns (rolls) in rectangular small-aspect-ratio containers as the aspect ratio is varied. For pure Marangoni flows at moderate Biot and Prandtl number, the transitions occurring when paths of codimension 1 singularities intersect determine to a large extent the multiplicity of stable patterns. These transitions also lead, for example, to Hopf bifurcations and stable periodic flows for a small range in aspect ratio. The influence of the type of lateral walls on the multiplicity of steady states is considered. 'No-slip' lateral walls lead to hysteresis effects and typically restrict the number of stable flow patterns (with respect to 'slippery' sidewalls) through the occurrence of saddle node bifurcations. In this way 'no-slip' sidewalls induce a selection of certain patterns, which typically have the largest Nusselt number, through secondary bifurcation.

Dijkstra, Henk A.↗
NASA NTRS1992

Viscous and inviscid instabilities of a trailing vortex

The linear stability of the trailing line vortex model of Batchelor (1964) is studied using a spectral collocation and matrix eigenvalue method. The entire unstable region in the swirl/axial wavenumber parameter space is mapped out for various azimuthal wavenumbers for both the inviscid and viscous stability problem. The results of the study provide a direct numerical validation of the large-azimuthal-wavenumber asymptotic analysis of Leibovich and Stewartson (1983). It is shown that accurate results are obtained up to azimuthal wavenumbers of 10,000 and greater, and the agreement with the asymptotic theory is excellent.

Mayer, Ernst W.↗
NASA NTRS1998

Comparison of Graphite Fabric Reinforced PMR-15 and Avimid N Composites After Long Term Isothermal Aging at Various Temperatures

Extensive effort is currently being expended to demonstrate the feasibility of using high-performance, polymer-matrix composites as engine structural materials over the expected operating lifetime of the aircraft, which can extend from 18,000 to 30,000 hr. The goal is to develop light-weight, high-strength, and high-modulus materials for use in higher temperature sections of advanced 21 st century aircraft propulsion systems. To accomplish this goal, it is necessary to pursue the development of thermal and mechanical durability models for graphite-fiber-reinforced, polymer-matrix composites. Numerous investigations have been reported regarding the thermo-oxidative stability (TOS) of the polyimide PMR-15 (1-5). A significant amount of this work has been directed at edge and geometry effects, reinforcement fiber influences, and empirical modeling of high-temperature weight loss behavior. It is yet to be determined if the information obtained from the PMR-15 composite tests is applicable to other polyimide-matrix composites. The condensation-curing polymer Avimid N is another advanced composite material often considered for structural applications at high temperatures. Avimid N has better thermo-oxidative stability than PMR-15 (6), but the latter is more easily processed. The most comprehensive study of the thermo-oxidative stability of Avimid N neat resin and composites at 371 (infinity)C is found in Salin and Seferis (7). The purposes of the work described herein were to compare the thermal aging behavior of these two matrix polymers and to determine the reasons for and the consequences of the difference in thermal durability. These results might be of some use in improving future polymer development through the incorporation of the desirable characteristics of both polyimides.

Bowles, Kenneth J.↗
NASA NTRS1999

Studies of Premixed Laminar and Turbulent Flames at Microgravity

Several topics relating to combustion limits in premixed flames at reduced gravity have been studied. These topics include: (1) flame balls; (2) numerical simulation of flame ball and planar flame structure and stability; (3) experimental simulation of buoyancy effects in premixed flames using aqueous autocatalytic reactions; and (4) premixed flame propagation in Hele-Shaw cells.

Abid, M.↗
NASA NTRS1999

Finite Differences and Collocation Methods for the Solution of the Two Dimensional Heat Equation

In this paper we combine finite difference approximations (for spatial derivatives) and collocation techniques (for the time component) to numerically solve the two dimensional heat equation. We employ respectively a second-order and a fourth-order schemes for the spatial derivatives and the discretization method gives rise to a linear system of equations. We show that the matrix of the system is non-singular. Numerical experiments carried out on serial computers, show the unconditional stability of the proposed method and the high accuracy achieved by the fourth-order scheme.

Kouatchou, Jules↗
NASA NTRS1999

Studies of Premixed Laminar and Turbulent Flames at Microgravity

Several topics relating to combustion limits in premixed flames at reduced gravity have been studied. These topics include: (1) flame balls; (2) numerical simulation of flame ball and planar flame structure and stability; (3) experimental simulation of buoyancy effects in premixed flames using aqueous autocatalytic reactions; and (4) premixed flame propagation in Hele-Shaw cells.

Abid, M.↗
NASA NTRS2002

Lattice Boltzmann Equation On a 2D Rectangular Grid

We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitras~y inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.

Bouzidi, MHamed↗
NASA NTRS2006

Large-Scale Hybrid Motor Testing

Hybrid rocket motors can be successfully demonstrated at a small scale virtually anywhere. There have been many suitcase sized portable test stands assembled for demonstration of hybrids. They show the safety of hybrid rockets to the audiences. These small show motors and small laboratory scale motors can give comparative burn rate data for development of different fuel/oxidizer combinations, however questions that are always asked when hybrids are mentioned for large scale applications are - how do they scale and has it been shown in a large motor? To answer those questions, large scale motor testing is required to verify the hybrid motor at its true size. The necessity to conduct large-scale hybrid rocket motor tests to validate the burn rate from the small motors to application size has been documented in several place^'^^.^. Comparison of small scale hybrid data to that of larger scale data indicates that the fuel burn rate goes down with increasing port size, even with the same oxidizer flux. This trend holds for conventional hybrid motors with forward oxidizer injection and HTPB based fuels. While the reason this is occurring would make a great paper or study or thesis, it is not thoroughly understood at this time. Potential causes include the fact that since hybrid combustion is boundary layer driven, the larger port sizes reduce the interaction (radiation, mixing and heat transfer) from the core region of the port. This chapter focuses on some of the large, prototype sized testing of hybrid motors. The largest motors tested have been AMROC s 250K-lbf thrust motor at Edwards Air Force Base and the Hybrid Propulsion Demonstration Program s 250K-lbf thrust motor at Stennis Space Center. Numerous smaller tests were performed to support the burn rate, stability and scaling concepts that went into the development of those large motors.

Story, George↗
NASA NTRS2009

Parametric Optimization of Ares I Propellant Slosh Characteristics Using Frequency Response Criteria

A novel technique for developing propellant slosh damping requirements with respect to the stability characteristics of large flexible launch vehicles is presented. A numerical algorithm is devised which allows an automated software program to rapidly converge to pseudo-optimal solutions that minimize required propellant slosh damping for multiple tanks while maintaining constraints on the frequency response characteristics of a particular open-loop plant transfer function. An implementation of the algorithm using a high-order linear model of the Ares I plant dynamics considers all relevant dynamic interactions of flexible body modes, propellant slosh, and nozzle inertia effects. A high-resolution propellant damping requirements table is produced that can be used for baffle design. The method is demonstrated to provide exceptional speed and accuracy when compared with the alternative human-in-the-loop approach.

Orr, Jeb S.↗
DOE OSTI2021

Borehole breakout modeling in arkose and granite rocks

The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.

58 GEOSCIENCES↗
NASA NTRS1990

Instability of a supersonic shock free elliptic jet

This paper presents a comparison of the measured and the computed spatial stability properties of an aspect ratio 2 supersonic shock free elliptic jet. The shock free nature of the elliptic jet provides an ideal test of validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with linear hydrodynamic stability theory. Both aerodynamic and acoustic data were measured. The data are used to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. A hybrid numerical scheme is applied to solve the Rayleigh problem governing the inviscid linear spatial stability of the jet. The measured mean velocity profiles are used to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles used in the stability analysis. Computational results are presented for several modes of instability at two jet cross sections. The acoustic measurements show that a varicose instability is the jet's perferred mode of motion. The stability analysis predicts that the Strouhal number varies linearly as a function of axial distance in the jet's initial mixing region, which is in good qualitative agreement with previous measurements.

Baty, Roy S.↗
DOE OSTI2025

Innovative Symmetrical Electrolyte Architecture Enables Ultra-Flat and Thermal Resilient Protonic Ceramic Electrochemical Cells

Protonic ceramic electrochemical cells (PCECs) represent a promising class of solid‐state energy conversion devices capable of high‐efficiency hydrogen production and power generation. However, the practical deployment of planar PCECs is fundamentally constrained by severe structural deformation and mechanical failure during fabrication, stemming from asymmetric shrinkage between the thin electrolyte and the thick NiO‐based support layer. Here, in this work, a functionally integrated, symmetry‐engineered double‐sided electrolyte (DE) design is unveiled, which not only suppresses thermally induced curvature but also unlocks significant gains in electrochemical performance and stability. This architecture intrinsically balances shrinkage dynamics across the cell bilaterally, enabling the fabrication of ultra‐flat 5 × 5 cm 2 cells with sub‐100 µm thickness variation. A numerical solid mechanics simulation is introduced to investigate and interpret this achievement. Beyond structural advantages, the DE configuration enhances the cell operational stability, delivering a low open‐circuit voltage degradation of 9.5 mV/100 h across 80 thermal cycles. This work establishes a compelling paradigm wherein architectural symmetry directly translates to both mechanical fidelity and functional enhancement, offering a promising route toward PCECs scale‐up.

08 - HYDROGEN↗
DOE OSTI2024

Wind Turbine Aeroelastic Stability in OpenFAST

Wind turbines are growing in size and increasingly suffer from aeroelastic instabilities. Unfortunately, numerical models often show inconsistent results during verification studies. We address this gap by first introducing novel linearization capabilities within the open-source aero-hydro-servo-elastic framework OpenFAST. Next, a code-to-code benchmark study is presented that compares modal parameters between OpenFAST and HAWCStab2 for a land-based version of the International Energy Agency 15-MW reference wind turbine modeled with quasi-steady aerodynamics. The two solvers are in strong agreement except for discrepancies in the second rotor flapwise modes. The differences are attributed to the torsional flexibility of the tower, which is assumed torsionally stiff in the OpenFAST model. Work is ongoing to close this modeling gap. The aeroelastic stability of a low-specific-power land-based wind turbine is also investigated. The impact of design choices is discussed, high-lighting how narrow the margins are between a stable design and an unstable design.

17 WIND ENERGY↗