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,369 records · Page 76

Integrated Locomotor Function Tests for Countermeasure Evaluation

Following spaceflight crewmembers experience locomotor dysfunction due to inflight adaptive alterations in sensorimotor function. Countermeasures designed to mitigate these postflight gait alterations need to be assessed with a new generation of tests that evaluate the interaction of various sensorimotor sub-systems central to locomotor control. The goal of the present study was to develop new functional tests of locomotor control that could be used to test the efficacy of countermeasures. These tests were designed to simultaneously examine the function of multiple sensorimotor systems underlying the control of locomotion and be operationally relevant to the astronaut population. Traditionally, gaze stabilization has been studied almost exclusively in seated subjects performing target acquisition tasks requiring only the involvement of coordinated eye-head movements. However, activities like walking involve full-body movement and require coordination between lower limbs and the eye-head-trunk complex to achieve stabilized gaze during locomotion. Therefore the first goal of this study was to determine how the multiple, interdependent, full-body sensorimotor gaze stabilization subsystems are functionally coordinated during locomotion. In an earlier study we investigated how alteration in gaze tasking changes full-body locomotor control strategies. Subjects walked on a treadmill and either focused on a central point target or read numeral characters. We measured: temporal parameters of gait, full body sagittal plane segmental kinematics of the head, trunk, thigh, shank and foot, accelerations along the vertical axis at the head and the shank, and the vertical forces acting on the support surface. In comparison to the point target fixation condition, the results of the number reading task showed that compensatory head pitch movements increased, peak head acceleration was reduced and knee flexion at heel-strike was increased. In a more recent study we investigated the adaptive remodeling of the full-body gaze control systems following exposure to visual-vestibular conflict. Subjects walked on a treadmill before and after a 30- minute exposure to 0.5X minifying during which self-generated sinusoidal vertical head rotations were performed while seated. Following exposure to visual-vestibular conflict subjects showed a restriction in compensatory head movements, increased knee and ankle flexion after heel-strike and a decrease in the rate of body loading during the rapid weight transfer phase after the heel strike event. Taken together, results from both studies provide evidence that the full body contributes to gaze stabilization during locomotion, and that different functional elements are responsive to changes in visual task constraints and are subject to adaptive alterations following exposure to visual-vestibular conflict. This information provides the basis for the design of a new generation of integrative tests that incorporate the evaluation of multiple neural control systems relevant to astronaut operational performance.

Bloomberg, J. J.↗

The effects of Kelvin-Helmholtz instability on resonance absorption layers in coronal loops

One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.

Karpen, Judith T.↗

Free vibration analysis of spinning structural systems.

This article presents an efficient digital computer procedure, along with the complete listing of the associated computer program, which may be conveniently utilized for the accurate solution of a wide range of practical eigenvalue problems. Important applications of the present work are envisaged in the natural frequency analysis of spinning structures discretized by the finite element technique, and in the determination of transfer functions associated with the dynamic blocks of control systems of spacecraft utilizing gas jets or reaction wheels for attitude control, as well as of spin-stabilized and dual-spin-stabilized satellites. The validity of the Sturm sequence property is first established for the related matrix formulation involving Hermitian and real symmetric, positive-definite matrices, both being usually of highly banded configuration. A numerically stable algorithm based on the Sturm sequence method is then developed which fully exploits the banded form of the associated matrices.

Gupta, K. K.↗

An implicit algorithm for the conservative transonic full potential equation using an arbitrary mesh

A new, implicit approximate factorization (AF) algorithm designed to solve the conservative full-potential equation for the transonic flow past arbitrary airfoils has been developed. The new algorithm uses an upwind bias of the density coefficient to provide stability in supersonic regions. This allows the simple two- and three-banded matrix form of the AF scheme to be retained over the entire flow field, even in regions of supersonic flow. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Airfoil pressure distributions have been computed and are in good agreement with independent results.

Holst, T. L.↗

Tradeoff studies in multiobjective insensitive design of airplane control systems

A computer aided design method for multiobjective parameter-insensitive design of airplane control systems is described. Methods are presented for trading off nominal values of design objectives against sensitivities of the design objectives to parameter uncertainties, together with guidelines for designer utilization of the methods. The methods are illustrated by application to the design of a lateral stability augmentation system for two supersonic flight conditions of the Shuttle Orbiter. Objective functions are conventional handling quality measures and peak magnitudes of control deflections and rates. The uncertain parameters are assumed Gaussian, and numerical approximations of the stochastic behavior of the objectives are described. Results of applying the tradeoff methods to this example show that stochastic-insensitive designs are distinctly different from deterministic multiobjective designs. The main penalty for achieving significant decrease in sensitivity is decreased speed of response for the nominal system.

Schy, A. A.↗

A class of optimum digital phase locked loops

This paper presents a class of optimum digital filters for digital phase locked loops, for the important case in which the maximum update rate of the loop filter and numerically controlled oscillator (NCO) is limited. This case is typical when the loop filter is implemented in a microprocessor. In these situations, pure delay is encountered in the loop transfer function and thus the stability and gain margin of the loop are of crucial interest. The optimum filters designed for such situations are evaluated in terms of their gain margin for stability, dynamic error, and steady-state error performance. For situations involving considerably high phase dynamics an adaptive and programmable implementation is also proposed to obtain an overall optimum strategy.

Kumar, R.↗

Navier-Stokes simulation of the crossflow instability in swept-wing flows

The computational modeling of the transition process characteristic of flows over swept wings are described. Specifically, the crossflow instability and crossflow/T-S wave interactions are analyzed through the numerical solution of the full three-dimensional Navier-Stokes equations including unsteadiness, curvature, and sweep. This approach is chosen because of the complexity of the problem and because it appears that linear stability theory is insufficient to explain the discrepancies between different experiments and between theory and experiments. The leading edge region of a swept wing is considered in a three-dimensional spatial simulation with random disturbances as the initial conditions. The work has been closely coordinated with the experimental program of Professor William Saric, examining the same problem. Comparisons with NASA flight test data and the experiments at Arizona State University were a necessary and an important integral part of this work.

Reed, Helen L.↗

Convective instability in differentially rotating disks

A normal mode analysis for nonaxisymmetric perturbations in a thin, differentially rotating disk with a vertical structure that is isothermal and convectively unstable is performed. The vertical gravity is assumed to be external and constant. The perturbation scale is assumed to be much shorter than the radius of the disk but comparable to or less than the thickness. The initial value problem is formulated in shearing coordinates. Dispersion relations are obtained for the three limiting cases of zero shear, axisymmetric perturbations, and small radial wavelengths. The full effects of shear are studied by integrating numerically the initial value problem. Nonaxisymmetric local Fourier modes are found to have a radial wavenumber that increases linearly with time in proportion to the shear times the azimuthal wavenumber. While Coriolis forces exert stabilizing effects on the convective modes, reducing their growth rate and the range of unstable wavelengths, shear has destablizing effects inasmuch as it reduces the epicyclic frequency at a given angular velocity. In a Keplerian disk, perturbations with azimuthal wavelengths about 2 times smaller than vertical wavelengths grow exponentially.

Ryu, Dongsu↗

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗

Influence of Two-Phase Thermocapillary Flow on Liquid Retention in Microscopic Pores

An important feature of screened propellant acquisition devices is the retention capability or maximum maintainable pressure difference across the porous barrier separating the liquid and gas. Previous experiments with liquid hydrogen showed a marked reduction in retention when the tank containing the device was pressurized with hydrogen vapor. These tests, however did not indicate any appreciable degradation in retention with helium pressurization or direct heating through the screen. The objective of this article is to determine if the thermocapillary convection arising from phase change in the microscopic pores of such screens could cause these disparities in performance. A numerical model of flow in a single pore suggests that the thermocapillary-induced gradient in liquid pressure along the surface can strongly affect surface morphology. In an evaporative environment, this gradient exerts a stabilizing influence on surface curvature, and preserves the momentum balance between the liquid and gas. With condensation, it causes a force imbalance and a destabilizing suction in the middle of the pore that reduces retention. Results also indicate that introducing an inert gas, such as helium, suppresses this retention loss mechanism by lowering thermocapillary circulation and its associated interfacial pressure gradient.

Schmidt, George R.↗

Fast, Conservative Algorithm for Solving the Transonic Full-Potential Equation

A fast, fully implicit approximate factorization algorithm designed to solve the conservative, transonic, full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., rotated differencing), thereby greatly enhancing the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted, finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line over-relaxation algorithms.

Holst, Terry L.↗

Hyper-X Research Vehicle (HXRV) Experimental Aerodynamics Test Program Overview

This paper provides an overview of the experimental aerodynamics test program to ensure mission success for the autonomous flight of the Hyper-X Research Vehicle (HXRV). The HXRV is a 12-ft long, 2700 lb lifting body technology demonstrator designed to flight demonstrate for the first time a fully airframe integrated scramjet propulsion system. Three flights are currently planned, two at Mach 7 and one at Mach 10, beginning in the fall of 2000. The research vehicles will be boosted to the prescribed scramjet engine test point where they will separate from the booster, stabilize. and initiate engine test. Following 5+ seconds of powered flight and 15 seconds of cowl-open tares, the cowl will close and the vehicle will fly a controlled deceleration trajectory which includes numerous control doublets for in-flight aerodynamic parameter identification. This paper reviews the preflight testing activities, wind tunnel models, test rationale. risk reduction activities, and sample results from wind tunnel tests supporting the flight trajectory of the HXRV from hypersonic engine test point through subsonic flight termination.

Holland, Scott D.↗

Hyper-X Research Vehicle (HXRV) Experimental Aerodynamics Test Program Overview

This paper provides an overview of the experimental aerodynamics test program to ensure mission success for the autonomous flight of the Hyper-X Research Vehicle (HXRV). The HXRV is a 12-ft long, 2700 lb lifting body technology demonstrator designed to flight demonstrate for the first time a fully airframe integrated scramjet propulsion system. Three flights are currently planned, two at Mach 7 and one at Mach 10, beginning in the fall of 2000. The research vehicles will be boosted to the prescribed scramjet engine test point where they will separate from the booster, stabilize, and initiate engine test. Following 5+ seconds of powered flight and 15 seconds of cow-open tares, the cowl will close and the vehicle will fly a controlled deceleration trajectory which includes numerous control doublets for in-flight aerodynamic parameter identification. This paper reviews the preflight testing activities, wind tunnel models, test rationale, risk reduction activities, and sample results from wind tunnel tests supporting the flight trajectory of the HXRV from hypersonic engine test point through subsonic flight termination.

Holland, Scott D.↗

Crew Restraint Design for the International Space Station

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. Another ISS task that requires special consideration with respect to restraints is robotic teleoperation. The Robot Systems Technology Branch at the NASA Johnson Space Center is developing a humanoid robot astronaut, or Robonaut. It is being designed to perform extravehicular activities (EVAs) in the hazardous environment of space. An astronaut located inside the ISS will remotely operate Robonaut through a telepresence control system. Essentially, the robot mimics every move the operator makes. This requires the operator to be stable enough to prevent inadvertent movements, while allowing the flexibility to accomplish the controlled movements of the robot. Some type of special purpose restraint will be required to operate Robonaut and similar devices.

Norris, Lena↗

Experimental Optimization of a Free-to-Rotate Wing for Small UAS

This paper discusses an experimental investigation conducted to optimize a free-to-rotate wing for use on a small unmanned aircraft system (UAS). Although free-to-rotate wings have been used for decades on various small UAS and small manned aircraft, little is known about how to optimize these unusual wings for a specific application. The paper discusses some of the design rationale of the basic wing. In addition, three main parameters were selected for "optimization", wing camber, wing pivot location, and wing center of gravity (c.g.) location. A small apparatus was constructed to enable some simple experimental analysis of these parameters. A design-of-experiment series of tests were first conducted to discern which of the main optimization parameters were most likely to have the greatest impact on the outputs of interest, namely, some measure of "stability", some measure of the lift being generated at the neutral position, and how quickly the wing "recovers" from an upset. A second set of tests were conducted to develop a response-surface numerical representation of these outputs as functions of the three primary inputs. The response surface numerical representations are then used to develop an "optimum" within the trade space investigated. The results of the optimization are then tested experimentally to validate the predictions.

Logan, Michael J.↗

Potential vorticity structure of Titan's polar vortices from Cassini CIRS observations

The Cassini mission has provided the best opportunity to date to extensively study the seasonal variation in Titan’s atmosphere, with observations spanning almost half a Titan year (𝐿𝑠 = 293 − 93°). An important feature in the Titan middle-atmosphere is the formation of a polar vortex. Observations have shown that an initially well-developed northern vortex enriched with trace gas species gradually breaks down after spring equinox as a new vortex emerges in southern winter. Here we use Cassini CIRS observations to derive the temperature and composition of the middle-atmosphere. We use the gradient wind equation to first estimate the mean zonal winds, and then the Potential Vorticity (PV) throughout Titan’s atmosphere over the timespan of the Cassini mission. PV is a useful diagnostic quantity for studying the dynamics of polar vortices because it is materially conserved for adiabatic and frictionless flows, and can be inverted to find all other dynamical fields. Our results show the formation of a strong zonal jet in the winter hemisphere, with wind velocities reaching 220 m/s, which is consistent with previous studies. An annular PV structure is also observed over the winter poles, whereby a ring of PV encircles a local minima over the pole. Such distributions are often found to be unstable without a restoring force, yet they are seen here in numerous observations in both the northern and southern hemispheres. A comparison with the annular Martian vortices shows that latent heat release from condensation or subsidence-induced adiabatic heating may explain the origin and stability of the annulus. Finally, we investigate the evolution of the size of the vortices and the role of strong PV gradients as a dynamical mixing barrier for trace gas species across the vortex edge. We find that longer lived gases are less confined to the vortex than those with shorter photochemical lifetimes.

Conor Andrew Nixon↗

Ten Years of SNPP VIIRS Reflective Solar Bands On-Orbit Calibration and Performance

The Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (SNPP) has successfully operated on-orbit for nearly ten years since its launch in October 2011, continuously making global observations and improving studies of changes in the Earth’s climate and environment. VIIRS has 22 spectral bands, among which 14 are reflective solar bands (RSBs) covering a spectral range from 0.41 to 2.25 μm. The SNPP VIIRS RSBs are primarily calibrated by the onboard solar diffuser (SD), with its on-orbit degradation tracked by an onboard SD stability monitor (SDSM). The near-monthly scheduled lunar observations, together with the sensor responses over stable ground targets, have contributed to the sensor’s mission-long on-orbit calibration and characterization. Numerous improvements have been made in the RSB calibration methodology since SNPP VIIRS was launched, and the RSB calibration has reached a mature stage after almost ten years of on-orbit operation. SNPP is a joint NASA/NOAA mission and there are two teams, the NASA VIIRS Calibration Support Team (VCST) and the NOAA VIIRS Sensor Data Record Team, which are dedicated to SNPP VIIRS on-orbit calibration. In this paper, we focus on the calibration performed by the NASA VCST. The SNPP VIIRS RSB calibration methodologies used to produce the calibration coefficient look up tables for the latest NASA Level 1B Collection 2 products are reviewed and the calibration improvements incorporated in this collection are described. Recent calibration changes include the removal of image striping caused by non-uniform degradation of the SD, improvements to the method for combining lunar and SD data, mitigation of the effects due a recent anomaly in the SD measurements, estimation of the SD degradation beyond 935 nm, and fitting strategy improvements for look-up table delivery. Overall, the SNPP VIIRS RSBs have performed well since its launch and continue to meet design specifications.

SNPP↗

Improvements of the Load Schedule for the Machine Calibration of a Strain-Gage Balance

The load schedule for the calibration of a six-component force balance in a calibration machine was improved. Now, single-component loads are repeated in regular intervals during the calibration. This approach has several advantages. First, the number of single-component loads increases to about twenty-two percent of all loads and load combinations. Consequently, more accurate numerical estimates of the primary bridge sensitivities can be obtained if global regression is used for the analysis of the calibration data. In addition, single-component repeats make it possible to track the stability of the applied loads during the calibration process. Finally, interactions of single-component repeats can be compared with interactions that are observed during the application of manual loads to the balance. Machine calibration and manual data sets of two force balances are used to illustrate benefits of the new load schedule. It is shown in the examples how differences between the observed interactions of machine calibration and manual data can be quantified. The suggested improvements can also be implemented in the load schedule for the machine calibration of a moment or direct-read balance as long as single-component loads are included that are described in the design load format of the balance.

wind tunnel test↗