Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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,315 records · Page 73

Research in Parallel Algorithms and Software for Computational Aerosciences

Phase I is complete for the development of a Computational Fluid Dynamics parallel code with automatic grid generation and adaptation for the Euler analysis of flow over complex geometries. SPLITFLOW, an unstructured Cartesian grid code developed at Lockheed Martin Tactical Aircraft Systems, has been modified for a distributed memory/massively parallel computing environment. The parallel code is operational on an SGI network, Cray J90 and C90 vector machines, SGI Power Challenge, and Cray T3D and IBM SP2 massively parallel machines. Parallel Virtual Machine (PVM) is the message passing protocol for portability to various architectures. A domain decomposition technique was developed which enforces dynamic load balancing to improve solution speed and memory requirements. A host/node algorithm distributes the tasks. The solver parallelizes very well, and scales with the number of processors. Partially parallelized and non-parallelized tasks consume most of the wall clock time in a very fine grain environment. Timing comparisons on a Cray C90 demonstrate that Parallel SPLITFLOW runs 2.4 times faster on 8 processors than its non-parallel counterpart autotasked over 8 processors.

Domel, Neal D.↗

Hover Testing of the NASA/Army/MIT Active Twist Rotor Prototype Blade

Helicopter rotor individual blade control promises to provide a mechanism for increased rotor performance and reduced rotorcraft vibrations and noise. Active material methods, such as piezoelectrically actuated trailing-edge flaps and strain-induced rotor blade twisting, provide a means of accomplishing individual blade control without the need for hydraulic power in the rotating system. Recent studies have indicated that controlled strain induced blade twisting can be attained using piezoelectric active fiber composite technology. In order to validate these findings experimentally, a cooperative effort between NASA Langley Research Center, the Army Research Laboratory, and the MIT Active Materials and Structures Laboratory has been developed. As a result of this collaboration an aeroelastically-scaled active-twist model rotor blade has been designed and fabricated for testing in the heavy gas environment of the Langley Transonic Dynamics Tunnel (TDT). The results of hover tests of the active-twist prototype blade are presented in this paper. Comparisons with applicable analytical predictions of active-twist frequency response in hovering flight are also presented.

Matthew L Wilbur↗

Validation of Airborne FMCW Radar Measurements of Snow Thickness Over Sea Ice in Antarctica

Antarctic sea ice and its snow cover are integral components of the global climate system, yet many aspects of their vertical dimensions are poorly understood, making their representation in global climate models poor. Remote sensing is the key to monitoring the dynamic nature of sea ice and its snow cover. Reliable and accurate snow thickness data are currently a highly sought after data product. Remotely sensed snow thickness measurements can provide an indication of precipitation levels, predicted to increase with effects of climate change in the polar regions. Airborne techniques provide a means for regional-scale estimation of snow depth and distribution. Accurate regional-scale snow thickness data will also facilitate an increase in the accuracy of sea ice thickness retrieval from satellite altimeter freeboard estimates. The airborne data sets are easier to validate with in situ measurements and are better suited to validating satellite algorithms when compared with in situ techniques. This is primarily due to two factors: better chance of getting coincident in situ and airborne data sets and the tractability of comparison between an in situ data set and the airborne data set averaged over the footprint of the antennas. A 28-GHz frequency modulated continuous wave (FMCW) radar loaned by the Center for Remote Sensing of Ice Sheets to the Australian Antarctic Division is used to measure snow thickness over sea ice in East Antarctica. Provided with the radar design parameters, the expected performance parameters of the radar are summarized. The necessary conditions for unambiguous identification of the airsnow and snowice layers for the radar are presented. Roughnesses of the snow and ice surfaces are found to be dominant determinants in the effectiveness of layer identification for this radar. Finally, this paper presents the first in situ validated snow thickness estimates over sea ice in Antarctica derived from an FMCW radar on a helicopterborne platform.

sea ice↗

Use of an OSSE to Evaluate Background Error Covariances Estimated by the 'NMC Method'

The NMC method has proven utility for prescribing approximate background-error covariances required by variational data assimilation systems. Here, untunedNMCmethod estimates are compared with explicitly determined error covariances produced within an OSSE context by exploiting availability of the true simulated states. Such a comparison provides insights into what kind of rescaling is required to render the NMC method estimates usable. It is shown that rescaling of variances and directional correlation lengths depends greatly on both pressure and latitude. In particular, some scaling coefficients appropriate in the Tropics are the reciprocal of those in the Extratropics. Also, the degree of dynamic balance is grossly overestimated by the NMC method. These results agree with previous examinations of the NMC method which used ensembles as an alternative for estimating background-error statistics.

Data Assimilation↗

Earth-Moon Libration Point Orbit Stationkeeping: Theory, Modeling and Operations

Collinear Earth-Moon libration points have emerged as locations with immediate applications. These libration point orbits are inherently unstable and must be maintained regularly which constrains operations and maneuver locations. Stationkeeping is challenging due to relatively short time scales for divergence effects of large orbital eccentricity of the secondary body, and third-body perturbations. Using the Acceleration Reconnection and Turbulence and Electrodynamics of the Moon's Interaction with the Sun (ARTEMIS) mission orbit as a platform, the fundamental behavior of the trajectories is explored using Poincare maps in the circular restricted three-body problem. Operational stationkeeping results obtained using the Optimal Continuation Strategy are presented and compared to orbit stability information generated from mode analysis based in dynamical systems theory.

Folta, David C.↗

Recent Improvements and Verification of A Full Body Model in Opensim

BACKGROUND: The dynamic feasibility [1,2] criterion, that the subject’s Center of Pressure (COP) be located within the Base of Support (BOS) which outlines the feet, has aided in assessing the stability of human motion recorded on earth while performing the recorded tasks in lunar gravity or as countermeasures exercises on a vibration isolation and stabilization system in microgravity. The convex hull of the BOS on the platform under the subject’s feet was estimated using virtual markers on the feet of the scaled subject model. The COP was calculated using the ground reaction forces and moments determined from motion capture data with biomechanical modeling tools [3]. Occasionally, large-amplitude oscillatory spikes or “artifacts” were observed in the subject’s linear and angular momentum derivatives, affecting some COP data derived from motion capture. The purpose of this investigation was to assess and improve the accuracy of model scaling and BOS estimation as well as to determine the efficacy of model adjustments in mitigating artifacts influencing motion capture-derived ground reaction force and COP results. METHODS: To aid evaluation of proposed process and model updates, motion capture data were collected for two subjects during unit test and range of motion trials, lunar tasks, and countermeasures exercise motions. Markers were added to the full body Plug-in Gait marker set [4] during data collection. New markers were placed on the front, back, sides, and top of the head to improve scaling using distances between marker pairs. Medial elbow markers were added to stabilize the upper arm during OpenSim Inverse Kinematics (IK) [5]. Finally, markers were added on the outer edge of the heels and on the outside edges of the first and last toes on each foot. These additional foot markers were made available to test new automated foot scaling techniques and to calculate the error between the subject’s estimated and recorded BOS. The modified unscaled OpenSim Full Body Rajagopal Model [6,7] was adjusted using some previously investigated techniques [8] to mitigate rapid shifts in joint angles occurring during IK, as these were found to cause the spike artifacts observed in subsequent stages of analysis. Since OpenSim models use Euler angles and rotation sequences, the arm axes of rotation were adjusted, and the pelvis order of rotation was changed to minimize the likelihood of encountering “gimbal lock” during common human motion. The model clavicle, arm, elbow, wrist, pelvis, and ankle angle limits were adjusted to better accommodate the full human range of motion seen in exercise and lunar data. The shoulder joint center was calculated using a “pivoting” algorithm [9], and both shoulder joint center and upper arm markers were included during IK to provide additional shoulder stability on a case-by-case basis. The quality of IK results was assessed by three criteria: minimizing error between recorded motion capture markers and model markers, checking for reasonable rates of change in joint angles between fames (i.e., no IK artifacts), and ensuring the absence of spikes in the inertial forces and angular momentum derivatives calculated using a custom OpenSim plugin [10]. RESULTS: The additional markers placed on the subject during data collection allowed the head and feet to be scaled more accurately using distances between new marker pairs. Scaling with BOS markers placed on the subject and removing the limit on subtalar angle resulted in more accurate BOS determination. Unrealistically large changes in joint angles between frames could be reduced by including clavicle, sternum, and medial elbow markers during IK. In cases with large arm ranges of motion, results could be further improved by running IK using medial elbow and virtual shoulder joint center markers. Model adjustments significantly improved the IK results affecting COP calculation and increased the accuracy of BOS estimation.

C A Bell↗

On the Connection Between Microbursts and Nonlinear Electronic Structures in Planetary Radiation Belts

Using a dynamical-system approach, we have investigated the efficiency of large-amplitude whistler waves for causing microburst precipitation in planetary radiation belts by modeling the microburst energy and particle fluxes produced as a result of nonlinear wave-particle interactions. We show that wave parameters, consistent with large amplitude oblique whistlers, can commonly generate microbursts of electrons with hundreds of keV-energies as a result of Landau trapping. Relativistic microbursts (greater than 1 MeV) can also be generated by a similar mechanism, but require waves with large propagation angles Theta (sub k)B greater than 50 degrees and phase-speeds v(sub phi) greater than or equal to c/9. Using our result for precipitating density and energy fluxes, we argue that holes in the distribution function of electrons near the magnetic mirror point can result in the generation of double layers and electron solitary holes consistent in scales (of the order of Debye lengths) to nonlinear structures observed in the radiation belts by the Van Allen Probes. Our results indicate a relationship between nonlinear electrostatic and electromagnetic structures in the dynamics of planetary radiation belts and their role in the cyclical production of energetic electrons (E greater than or equal to 100 keV) on kinetic timescales, which is much faster than previously inferred.

plasmas↗

Solid-state Precipitation in Molecular Dynamics: KMC-MD Hybrid Simulations

A critical drawback of molecular dynamics (MD) simulations when modeling physical properties is the limitation to short timescales. Solid state precipitation is one behavior that is not significantly observable at timescales accessible to MD, and thus has been difficult to simulate at an MD scale. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation during MD simulation implemented in LAMMPS. This method incorporates the atomic neighbor hopping behavior utilized by KMC simulations into a MD simulation with a variable local lattice configuration. By accelerating nearest-neighbor atomic hopping with this technique, solid-state precipitation can be observed at an MD timescale. Utilizing theory based on persistent-embryo methods for liquid-solid precipitation in MD we demonstrate the calculation of critical nucleus size for the gamma’ precipitation out of gamma phase in the Ni-Al system.

Molecular Dynamics↗

Computational Fluid Dynamic (CFD) analysis of axisymmetric plume and base flow of film/dump cooled rocket nozzle

Film/dump cooling a rocket nozzle with fuel rich gas, as in the National Launch System (NLS) Space Transportation Main Engine (STME), adds potential complexities for integrating the engine with the vehicle. The chief concern is that once the film coolant is exhausted from the nozzle, conditions may exist during flight for the fuel-rich film gases to be recirculated to the vehicle base region. The result could be significantly higher base temperatures than would be expected from a regeneratively cooled nozzle. CFD analyses were conduced to augment classical scaling techniques for vehicle base environments. The FDNS code with finite rate chemistry was used to simulate a single, axisymmetric STME plume and the NLS base area. Parallel calculations were made of the Saturn V S-1 C/F1 plume base area flows. The objective was to characterize the plume/freestream shear layer for both vehicles as inputs for scaling the S-C/F1 flight data to NLS/STME conditions. The code was validated on high speed flows with relevant physics. This paper contains the calculations for the NLS/STME plume for the baseline nozzle and a modified nozzle. The modified nozzle was intended to reduce the fuel available for recirculation to the vehicle base region. Plumes for both nozzles were calculated at 10kFT and 50kFT.

Tucker, P. K.↗

Investigation of the Physical Processes Governing Large-scale Tracer Transport in the Stratosphere and Troposphere

This report reviews the second year of a three-year research program to investigate the physical mechanisms which underlie the transport of trace constituents in the stratosphere- troposphere system. The primary scientific goal of the research is to identify the processes which transport air masses within the lower stratosphere, particularly between the tropics and middle latitudes. The SASS program seeks to understand the impact of the present and future fleets of conventional jet traffic on the upper troposphere and lower stratosphere, while complementary airborne observations under UARP seek to understand the complex interactions of dynamical and chemical processes that affect the ozone layer. The present investigation contributes to the goals of each of these by diagnosing the history of air parcels intercepted by NASA research aircraft in UARP and AEAP campaigns.

Selkirk, Henry B.↗

Transport Aircraft System Identification from Wind Tunnel Data

Recent studies have been undertaken to investigate and develop aerodynamic models that predict aircraft response in nonlinear unsteady flight regimes for transport configurations. The models retain conventional static and rotary dynamic terms but replace conventional acceleration terms with more general indicial functions. In the Integrated Resilient Aircraft Controls project of the NASA Aviation Safety Program one aspect of the research is to apply these current developments to transport configurations to facilitate development of advanced controls technology. This paper describes initial application of a more general modeling methodology to the NASA Langley Generic Transport Model, a sub-scale flight test vehicle.

Murphy, Patrick C.↗

Solid-state Precipitation in Molecular Dynamics: kMC-MD Hybrid Simulations

A critical drawback of molecular dynamics (MD) simulations when modeling physical properties is the limitation to short timescales. Solid state precipitation is one behavior that is not significantly observable at timescales accessible to MD, and thus has been difficult to simulate at an MD scale. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation during MD simulation implemented in LAMMPS. This method incorporates the atomic neighbor hopping behavior utilized by KMC simulations into a MD simulation with a variable local lattice configuration. By accelerating nearest-neighbor atomic hopping with this technique, solid-state precipitation can be observed at an MD timescale. Utilizing theory based on persistent-embryo methods for liquid-solid precipitation in MD we demonstrate the calculation of critical nucleus size for the gamma’ precipitation out of gamma phase in the Ni-Al system.

Molecular Dynamics↗

Dissipational galaxy formation - Confrontation with observations

An exploration is presented of the hypothesis that a protogalaxy consists of an aggregate of interacting gas clouds which undergo mergers with neighboring systems, as envisaged by both the hierarchical clustering and fragmentation schemes of galaxy formation. Both gaseous dissipation and violet relaxation play fundamental roles in this galaxy formation model, in order to account for such diverse structural and dynamical properties of spheroidal galaxies as velocity anisotropy and metallicity gradients. Protogalaxy mergers during the initial stages of galaxy clustering can account for the observed spatial distribution of spiral, S0, and elliptical galaxies, and galaxy formation can occur slowly and at late epochs, since the time scale for disk formation is less than about 10 to the 10th years.

Silk, J.↗

Automatic control of cryogenic wind tunnels

Inadequate Reynolds number similarity in testing of scaled models affects the quality of aerodynamic data from wind tunnels. This is due to scale effects of boundary-layer shock wave interaction which is likely to be severe at transonic speeds. The idea of operation of wind tunnels using test gas cooled to cryogenic temperatures has yielded a quantrum jump in the ability to realize full scale Reynolds number flow similarity in small transonic tunnels. In such tunnels, the basic flow control problem consists of obtaining and maintaining the desired test section flow parameters. Mach number, Reynolds number, and dynamic pressure are the three flow parameters that are usually required to be kept constant during the period of model aerodynamic data acquisition. The series of activity involved in modeling, control law development, mechanization of the control laws on a microcomputer, and the performance of a globally stable automatic control system for the 0.3-m Transonic Cryogenic Tunnel (TCT) are discussed. A lumped multi-variable nonlinear dynamic model of the cryogenic tunnel, generation of a set of linear control laws for small perturbation, and nonlinear control strategy for large set point changes including tunnel trajectory control are described. The details of mechanization of the control laws on a 16 bit microcomputer system, the software features, operator interface, the display and safety are discussed. The controller is shown to provide globally stable and reliable temperature control to + or - 0.2 K, pressure to + or - 0.07 psi and Mach number to + or - 0.002 of the set point value. This performance is obtained both during large set point commands as for a tunnel cooldown, and during aerodynamic data acquisition with intrusive activity like geometrical changes in the test section such as angle of attack changes, drag rake movements, wall adaptation and sidewall boundary-layer removal. Feasibility of the use of an automatic Reynolds number control mode with fixed Mach number control is demonstrated.

Balakrishna, S.↗

Introduction to Penelope

A formal program verification is a (mathematical) proof that a program executed according to its intended model meets some specification. This proves that the algorithm defined by the program is correct in the precise technical sense of being consistent with a particular specification. A program correct in this sense is free from a large and important class of errors, even though its behavior may still produce unintended results--either because the implementation of the programming language itself does not match the model of execution, or because the specification does not correctly express the user's intentions. Penelope is a prototype system for interactively developing and verifying programs that are written in a rich subset of sequential Ada. Penelope can be used to develop a program and its correctness proof incrementally, and in concert with one another. Incrementality is used in a number of ways to help make verification more tractable and more productive. For example, if an already-verified program is modified, one can attempt to prove the modified version by replaying and modifying the original verification. Penelope's specification language, Larch/Ada, belongs to the family of Larch interface languages. Larch/Ada scales up properly, in the sense that it is demonstrably sound to decompose a system hierarchically and reason locally about the implementation of each piece. Penelope has been applied in various demonstration projects--for specification (guidance control, distributed operating systems), verification (of off-the-shelf code), and formal development (by non-expert as well as expert users). Some features of Penelope have been embodied in Ada Wise, a lint-like non-interactive tool that warns of the potential for certain dynamic semantic errors in Ada programs.

Guaspari, David↗

Weather from the Stratosphere?

Is the stratosphere, the atmospheric layer between about 10 and 50 km, important for predicting changes in weather and climate? The traditional view is that the stratosphere is a passive recipient of energy and waves from weather systems in the underlying troposphere, but recent evidence suggests otherwise. At a workshop in Whistler, British Columbia (1), scientists met to discuss how the stratosphere responds to forcing from below, initiating feedback processes that in turn alter weather patterns in the troposphere. The lowest layer of the atmosphere, the troposphere, is highly dynamic and rich in water vapor, clouds, and weather. The stratosphere above it is less dense and less turbulent (see the figure). Variability in the stratosphere is dominated by hemispheric-scale changes in airflow on time scales of a week to several months. Occasionally, however, stratospheric air flow changes dramatically within just a day or two, with large-scale jumps in temperature of 20 K or more. The troposphere influences the stratosphere mainly through atmospheric waves that propagate upward. Recent evidence shows that the stratosphere organizes this chaotic wave forcing from below to create long-lived changes in the stratospheric circulation. These stratospheric changes can feed back to affect weather and climate in the troposphere.

Baldwin, Mark P.↗

Similarity Rules for Scaling Solar Sail Systems

Future science missions will require solar sails on the order of 200 square meters (or larger). However, ground demonstrations and flight demonstrations must be conducted at significantly smaller sizes, due to limitations of ground-based facilities and cost and availability of flight opportunities. For this reason, the ability to understand the process of scalability, as it applies to solar sail system models and test data, is crucial to the advancement of this technology. This paper will approach the problem of scaling in solar sail models by developing a set of scaling laws or similarity criteria that will provide constraints in the sail design process. These scaling laws establish functional relationships between design parameters of a prototype and model sail that are created at different geometric sizes. This work is applied to a specific solar sail configuration and results in three (four) similarity criteria for static (dynamic) sail models. Further, it is demonstrated that even in the context of unique sail material requirements and gravitational load of earth-bound experiments, it is possible to develop appropriate scaled sail experiments. In the longer term, these scaling laws can be used in the design of scaled experimental tests for solar sails and in analyzing the results from such tests.

Canfield, Stephen L.↗

High Resolution Global Climate Modeling with GEOS-5: Intense Precipitation, Convection and Tropical Cyclones on Seasonal Time-Scales.

In 2008 the World Modeling Summit for Climate Prediction concluded that "climate modeling will need-and is ready-to move to fundamentally new high-resolution approaches to capitalize on the seamlessness of the weather-climate continuum." Following from this, experimentation with very high-resolution global climate modeling has gained enhanced priority within many modeling groups and agencies. The NASA Goddard Earth Observing System model (GEOS-5) has been enhanced to provide a capability for the execution at the finest horizontal resolutions POS,SIOle with a global climate model today. Using this high-resolution, non-hydrostatic version of GEOS-5, we have developed a unique capability to explore the intersection of weather and climate within a seamless prediction system. Week-long weather experiments, to mUltiyear climate simulations at global resolutions ranging from 3.5- to 14-km have demonstrated the predictability of extreme events including severe storms along frontal systems, extra-tropical storms, and tropical cyclones. The primary benefits of high resolution global models will likely be in the tropics, with better predictions of the genesis stages of tropical cyclones and of the internal structure of their mature stages. Using satellite data we assess the accuracy of GEOS-5 in representing extreme weather phenomena, and their interaction within the global climate on seasonal time-scales. The impacts of convective parameterization and the frequency of coupling between the moist physics and dynamics are explored in terms of precipitation intensity and the representation of deep convection. We will also describe the seasonal variability of global tropical cyclone activity within a global climate model capable of representing the most intense category 5 hurricanes.

Putnam, WilliamM.↗