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,279 records · Page 71

Drogue Parachute Effects on the Orion Crew Module Stability

A forced oscillation test of the Orion Crew Module (CM) was conducted in the Langley 20-Foot Vertical Spin Tunnel. The objective of the test was to quantify the rate damping characteristics of the CM-drogue chute system. Numerous configurations were tested to measure the influence of the chutes on the CM dynamic aerodynamics and, conversely, the influence of the CM on drogue performance. Results show that the CM-drogue system is well-damped at all combinations of frequency, amplitude, and Strouhal number. The wake of the CM significantly reduces the drogue chute riser line force, and the drogues have little upstream influence on the CM aerodynamics. These results are being used to improve simulation model fidelity of CM flight with drogues deployed, which has been identified by the project as key to a successful Orion Critical Design Review.

Aubuchon, Vanessa V.↗

BLSTA: A boundary layer code for stability analysis

A computer program is developed to solve the compressible, laminar boundary-layer equations for two-dimensional flow, axisymmetric flow, and quasi-three-dimensional flows including the flow along the plane of symmetry, flow along the leading-edge attachment line, and swept-wing flows with a conical flow approximation. The finite-difference numerical procedure used to solve the governing equations is second-order accurate. The flow over a wide range of speed, from subsonic to hypersonic speed with perfect gas assumption, can be calculated. Various wall boundary conditions, such as wall suction or blowing and hot or cold walls, can be applied. The results indicate that this boundary-layer code gives velocity and temperature profiles which are accurate, smooth, and continuous through the first and second normal derivatives. The code presented herein can be coupled with a stability analysis code and used to predict the onset of the boundary-layer transition which enables the assessment of the laminar flow control techniques. A user's manual is also included.

Wie, Yong-Sun↗

Stability of semidiscrete approximations for hyperbolic initial-boundary-value problems: An eigenvalue analysis

A hyperbolic initial-boundary-value problem can be approximated by a system of ordinary differential equations (ODEs) by replacing the spatial derivatives by finite-difference approximations. The resulting system of ODEs is called a semidiscrete approximation. A complication is the fact that more boundary conditions are required for the spatially discrete approximation than are specified for the partial differential equation. Consequently, additional numerical boundary conditions are required and improper treatment of these additional conditions can lead to instability. For a linear initial-boundary-value problem (IBVP) with homogeneous analytical boundary conditions, the semidiscrete approximation results in a system of ODEs of the form du/dt = Au whose solution can be written as u(t) = exp(At)u(O). Lax-Richtmyer stability requires that the matrix norm of exp(At) be uniformly bounded for O less than or = t less than or = T independent of the spatial mesh size. Although the classical Lax-Richtmyer stability definition involves a conventional vector norm, there is no known algebraic test for the uniform boundedness of the matrix norm of exp(At) for hyperbolic IBVPs. An alternative but more complicated stability definition is used in the theory developed by Gustafsson, Kreiss, and Sundstrom (GKS). The two methods are compared.

Warming, Robert F.↗

International Space Station Crew Restraint Design

With permanent human presence onboard the International Space Station (ISS), crews will be living and working in microgravity, dealing with the challenges of a weightless environment. In addition, the confined nature of the spacecraft environment results in ergonomic challenges such as limited visibility and access to the activity areas, as well as prolonged periods of unnatural postures. Without optimum restraints, crewmembers may be handicapped for performing some of the on-orbit tasks. Currently, many of the tasks on ISS are performed with the crew restrained merely by hooking their arms or toes around handrails to steady themselves. This is adequate for some tasks, but not all. There have been some reports of discomfort/calluses on the top of the toes. In addition, this type of restraint is simply insufficient for tasks that require a large degree of stability. Glovebox design is a good example of a confined workstation concept requiring stability for successful use. They are widely used in industry, university, and government laboratories, as well as in the space environment, and are known to cause postural limitations and visual restrictions. Although there are numerous guidelines pertaining to ventilation, seals, and glove attachment, most of the data have been gathered in a 1-g environment, or are from studies that were conducted prior to the early 1980 s. Little is known about how best to restrain a crewmember using a glovebox in microgravity. In 2004, The Usability Testing and Analysis Facility (UTAF) at the NASA Johnson Space Center completed development/evaluation of several design concepts for crew restraints to meet the various needs outlined above. Restraints were designed for general purpose use, for teleoperation (Robonaut) and for use with the Life Sciences Glovebox. All design efforts followed a human factors engineering design lifecycle, beginning with identification of requirements followed by an iterative prototype/test cycle. Anthropometric modeling was completed for the 5th percentile Asian female and the 95th percentile American male for all restraints. A series of three evaluations was performed onboard NASA's reduced gravity aircraft (KC-135). For all evaluations, participants performed representative tasks while being videotaped, and then completed a questionnaire following each flight day. The questionnaire included ratings scales and free format questions to assess topics such as comfort, stability provided, flexibility provided, etc. Results from the three flight evaluations are being used to develop the human factors design requirements for crew restraint concepts for 1) general purpose restraints, 2) teleoperation restraints and 3) glovebox restraints. The poster presentation will describe the detailed methodology used, results from each of the three evaluations, and the resulting human factors recommendations for the design of these restraints.

Whitmore, M.↗

Numerical analysis of confined turbulent flow

The considered investigation is concerned with the development of an efficient computational method for obtaining a physical understanding of an internal turbulent field. The employed approach makes use of a 'two equation' type model for the turbulence to obtain the numerical solution of a two-dimensional confined turbulent flow. The mean flow governing equations are considered along with the governing equation of the mean temperature and concentrations, and the boundary conditions. The numerical procedure for solving the turbulent flow is discussed, taking into account an approximation to the nonlinear terms, and the inner and outer coupling. Attention is given to a stability convergence analysis, the stability characteristics, and computational examples.

Lin, A.↗

Theoretical study of multiple equilibria in simple axisymmetric tropical circulations

The possibility that the asymmetric part of the atmospheric circulation can possess multiple equilibrium states is examined using a two-layer axisymmetric model involving balance equations on an equatorial beta plane. Mountains are excluded from consideration and a Newtonian cooling formulation represents thermal forcing. A temperature maximum at 25 deg N is selected to simulate summer conditions in the Northern Hemisphere. Steady-state solutions obtained are investigated for stability with regard to first and second y-mode perturbations. A single stable mode is found, together with two other quasi-stable states. Attention is given to numerically modeling multiple equilibria in symmetric circulations, and one steady-state is determined for the two-layer model. A model employing primitive equations with the Boussinesq approximation is also examined, and it also furnishes only one steady state. The reasons for the lack of multiple steady-states as derived by the models are discussed.

Goswami, B. N.↗

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗

Flexible system model reduction and control system design based upon actuator and sensor influence functions

A model reduction technique based on aggregation with respect to sensor and actuator influence functions rather than modes is presented for large systems of coupled second-order differential equations. Perturbation expressions which can predict the effects of spillover on both the reduced-order plant model and the neglected plant model are derived. For the special case of collocated actuators and sensors, these expressions lead to the derivation of constraints on the controller gains that are, given the validity of the perturbation technique, sufficient to guarantee the stability of the closed-loop system. A case study demonstrates the derivation of stabilizing controllers based on the present technique. The use of control and observation synthesis in modifying the dimension of the reduced-order plant model is also discussed. A numerical example is provided for illustration.

Yam, Yeung↗

Development of high resolution simulations of the atmospheric environment using the MASS model

Numerical simulations were performed with a very high resolution (7.25 km) version of the MASS model (Version 4.0) in an effort to diagnose the vertical wind shear and static stability structure during the Shuttle Challenger disaster which occurred on 28 January 1986. These meso-beta scale simulations reveal that the strongest vertical wind shears were concentrated in the 200 to 150 mb layer at 1630 GMT, i.e., at about the time of the disaster. These simulated vertical shears were the result of two primary dynamical processes. The juxtaposition of both of these processes produced a shallow (30 mb deep) region of strong vertical wind shear, and hence, low Richardson number values during the launch time period. Comparisons with the Cape Canaveral (XMR) rawinsonde indicates that the high resolution MASS 4.0 simulation more closely emulated nature than did previous simulations of the same event with the GMASS model.

Kaplan, Michael L.↗

Further development towards a new virtual-pulse time integral methodology for general nonlinear transient thermal analysis

New developments describing the theoretical basis towards an effective virtual-pulse time integral procedure are presented for general nonlinear transient heat transfer problems. To validate the proposed methodology of computation, numerical test examples are given and comparisons are drawn with the implicit second-order accurate Crank-Nicolson technique. The proposed method shows superior or improved accuracy and stability characteristics for the models tested.

Tamma, K. K.↗

On the stability of nongyrotropic ion populations - A first (analytic and simulation) assessment

The wave and dispersion equations for perturbations propagating parallel to an ambient magnetic field in magnetoplasmas with nongyrotropic ion populations show, in general, the occurrence of coupling between the parallel (left- and right-hand circularly polarized electromagnetic and longitudinal electrostatic) eigenmodes of the associated gyrotropic medium. These interactions provide a means to driving linearly one mode with free-energy sources of other modes in homogeneous media. Different types of nongyrotropy bring about distinct classes of coupling. The stability of a hydrogen magnetoplasma with anisotropic, nongyrotropic protons that only couple the electromagnetic modes to each other is investigated analytically (via solution of the derived dispersion equation) and numerically (via simulation with a hybrid code). Nongyrotropy enhances growth and enlarges the unstable spectral range relative to the corresponding gyrotropic situation. The relevance of the properties of nongyrotropic populations to space plasma environments is also discussed.

Brinca, A. L.↗

Control of Combustion-Instabilities Through Various Passive Devices

It is well known that under some operating conditions, rocket engines (using solid or liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and shutdown. The sources of these instabilities are diverse and are dependent on fuel, chamber geometry and various upstream sources such as pumps, valves and injection mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system [1, 2]. Giammar and Putnam [3] performed a comprehensive study of noise generated by gasfired industrial burners and made several key observations; flow noise was sometimes more intense than combustion roar, which tended to have a characteristic frequency spectrum. Turbulence was amplified by the flame. The noise power varied directly with combustion intensity and also with the product of pressure drop and heat release rate. Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the combustion chamber. This is important, since it quantified how much of the noise from an engine originates in the combustor. A physical interpretation of the interchange of energy between sound waves and unsteady heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et al [6] extended Rayleigh s criterion to describe the interaction of unsteady combustion with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and energy conservation equations in the pre- and post-flame zones were matched by making several assumptions about the combustion process. They concluded that changes in boundary conditions affect the energy balance of acoustic waves in the combustor. Abouseif et al [7] also solved the one-dimensional flow equations, but they used a onestep reaction to evaluate the unsteady heat release rate by relating it to temperature and velocity perturbations. Their analysis showed that oscillations arise from coupling between entropy waves produced at the flame and pressure waves originating from the nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity and pressure oscillations. Conservation equations were expressed in integral form and solutions for the acoustic wave equations and complex frequencies were obtained. The imaginary part of the frequency indicated stability regions of the flame. Activation energy asymptotics together with a one-step reaction were used by McIntosh [9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional flames. He found that the flame stability was altered by the upstream acoustic feedback. Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional longitudinal, combustion instabilities. However, numerical diffusion was not completely eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure perturbations with premixed flames. He used complex chemistry to study responses of pressure perturbations in one-dimensional combustors. His results indicated that reflected and transmitted waves differed significantly from incident waves.

Frendi, Kader↗

Managing the Risk for Early Onset Osteoporosis in Long-Duration Astronauts Due to Spaceflight

Early Onset Osteoporosis is probably the most recognized but poorly understood long-term health risk due to spaceflight. Osteoporosis management is primarily prophylactic and clinical interventions rely upon the ability to predict fractures which is currently determined by surrogate measures of bone strength. The RMAT for Early Onset Osteoporosis identified some open issues related to the fact that long-duration astronauts compose a unique group of subjects for which clinical approaches for osteoporosis management do not apply. Long-duration astronauts are healthy, young (25 to 55 years of age), predominantly male, and physical fit relative to the typical osteoporosis patient. Moreover, during prolonged space missions (typically 6-month missions) the skeleton not only adapts to weightlessness, but is influenced by numerous risk factors induced by operational constraints, e.g., inability to maintain preflight weight-bearing and aerobic activities, sub-optimal dietary intake (e.g., high sodium content for food stability, lack of fresh fruit and vegetables), suppression of vitamin D metabolism by uv shielding, and remote medicine care. Moreover, adaptation results in novel changes to astronauts bones that cannot be detected by current medically-useful measures. Consequently, a panel of clinicians (recognized leaders and policy-makers in osteoporosis) was convened to review the dataset of bone measures and bone loss risk factors in long-duration astronauts. Driven by the queries in the RMAT, the panel was charged to determine 1) if an intervention is required to prevent this risk, 2) what type and at what time would intervention be optimal, 3) what is the clinical trigger that would require a medical response from flight surgeons and 4) how should research data be used in the clinical care of astronauts. Hence, the RMAT determined that a bone health policy need to be formulated specific for this unique cohort subjected to a novel skeletal condition

Sibonga, Jean D.↗

Pd/CeO2/SiC Chemical Sensors

The incorporation of nanostructured interfacial layers of CeO2 has been proposed to enhance the performances of Pd/SiC Schottky diodes used to sense hydrogen and hydrocarbons at high temperatures. If successful, this development could prove beneficial in numerous applications in which there are requirements to sense hydrogen and hydrocarbons at high temperatures: examples include monitoring of exhaust gases from engines and detecting fires. Sensitivity and thermal stability are major considerations affecting the development of high-temperature chemical sensors. In the case of a metal/SiC Schottky diode for a number of metals, the SiC becomes more chemically active in the presence of the thin metal film on the SiC surface at high temperature. This increase in chemical reactivity causes changes in chemical composition and structure of the metal/SiC interface. The practical effect of the changes is to alter the electronic and other properties of the device in such a manner as to degrade its performance as a chemical sensor. To delay or prevent these changes, it is necessary to limit operation to a temperature <450 C for these sensor structures. The present proposal to incorporate interfacial CeO2 films is based partly on the observation that nanostructured materials in general have potentially useful electrical properties, including an ability to enhance the transfer of electrons. In particular, nanostructured CeO2, that is CeO2 with nanosized grains, has shown promise for incorporation into hightemperature electronic devices. Nanostructured CeO2 films can be formed on SiC and have been shown to exhibit high thermal stability on SiC, characterized by the ability to withstand temperatures somewhat greater than 700 C for limited times. The exchanges of oxygen between CeO2 and SiC prevent the formation of carbon and other chemical species that are unfavorable for operation of a SiC-based Schottky diode as a chemical sensor. Consequently, it is anticipated that in a Pd/CeO2/SiC Schottky diode, the nanostructured interfacial CeO2 layer would contribute to thermal stability and, by contributing to transfer of electrons, would also contribute to sensitivity.

Lu, Weijie↗

An Experimental and Numerical Study of Roughness-Induced Instabilities in a Mach 3.5 Boundary Layer

Progress on a joint experimental and numerical study of laminar-to-turbulent transition induced by an isolated roughness element in a high-speed laminar boundary layer is reported in this paper. The numerical analysis suggests that transition is driven by the instability of high- and low-speed streaks embedded in the wake of the isolated roughness element. In addition, spatial stability analysis revealed that the wake flow supports multiple modes (even and odd) of convective instabilities that experience strong enough growth to cause transition. The experimental measurements, which included hot-wire and pitot-probe surveys, confirmed the existence of embedded high- and low-speed streaks in the roughness wake. Furthermore, the measurements indicate the presence of both even and odd modes of instability, although their relative magnitude depends on the specifics of the roughness geometry and flow conditions (e.g., the value of Re(sub kk) or k/delta. For the two test cases considered in the measurements (Re(sub kk) values of 462 and 319), the even mode and the odd mode were respectively dominant and appear to play a primary role in the transition process.

Kegerise, Michael A.↗

A Flight-Calibrated Methodology for Determination of Cassini Thruster On-Times for Reaction Wheel Biases

The Cassini spacecraft, the largest and most complex interplanetary spacecraft ever built, continues to undertake unique scientific observations of planet Saturn, Titan, Enceladus, and other moons of the ring world. In order to maintain a stable attitude during the course of its mission, this three-axis stabilized spacecraft uses two different control systems: the Reaction Control System (or RCS) and the Reaction Wheel Assembly (RWA) control system. In the course of its mission, Cassini performs numerous reaction wheel momentum biases (or unloads) using its reaction control thrusters. The use of the RCS thrusters often imparts undesired velocity changes (delta Vs) on the spacecraft and it is crucial for Cassini navigation and attitude control teams to be able to, quickly but accurately, predict the hydrazine usage and delta V vector in Earth Mean Equatorial (J2000) inertial coordinates for reaction wheel bias events, without actually having to spend time and resources simulating the event in a dynamic or hardware-in-the-loop simulation environments. The flight-calibrated methodology described in this paper, and the ground software developed thereof, are designed to provide the RCS thruster on-times, with acceptable accuracy and without any form of dynamic simulation, for reaction wheel biases, along with the hydrazine usage and the delta V in EME-2000 inertial frame.

Reaction Control System↗

On the dependence of travelling wave stability on basic state wave phase speed

The dependence of stability on basic state wave phase speed is examined in a nondivergent barotropic model on a sphere for planetary scale waves with phase speeds typical of waves observed in the atmosphere. Triad interactions are examined analytically and compared to numerical results from a full spectral model. A number of triads may result in growth of the perturbation components for each basic state wave; for each triad there is a basic state wave phase speed where it is resonant, identified as the point where the critical amplitude for instability to that perturbation becomes zero. Critical amplitudes for instability obtained using the full spectral model generally agree well with triad results. Since the basic state wave phase speed determines which triad will grow, spatial structure and critical amplitudes for instability for the growing disturbances depend strongly on basic state wave phase speed. The results of this idealized study suggest that phase speed may be an important factor in determining the stability of planetary scale waves in the atmosphere.

Manney, Gloria L.↗