Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time discretization”

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 721 records · Page 40

Accelerated viscoelastic characterization of T300-5208 graphite-epoxy laminates

A viscoelastic response scheme for the accelerated characterization of polymer-based composite laminates in applied to T300/5208 graphite/epoxy. The response of uni-directional specimens is modeled. The transient component of the viscoelastic creep compliance is assumed to follow a power law approximation. A recursive relationship is developed, based upon the Schapery single-integral equation, which allows approximation of a continuous time-varying uniaxial load using discrete steps in stress. The viscoelastic response of T300/5208 to transverse normal and shear stresses is determined unsing 90 deg and 10 deg off-axis tensile specimens. In each case the seven viscoelastic material parameters required in the analysis are determined experimentally using short-term creep and creep recovery tests. It is shown that an accurate measure of the power law exponent is crucial for accurate long-term prediction. A short term test cycle selection procedure is proposed, which should provide useful guidelines for the evaluation of other viscoelastic materials.

Tuttle, M. E.↗

Control-System Design Program

Control-theory design package, Optimal Regulator Algorithms for Control of Linear Systems (ORACLS), developed to aid in design of controllers and optimal filters for systems modeled by linear, time-invariant differential and difference equations. Optimal linear quadratic regulator theory, Linear-Quadratic-Gaussian (LQG) problem, most widely accepted method of determining optimal control policy. Provides for solution to time-in-variant continuous or discrete LQG problems. Attractive to control-system designer providing rigorous tool for dealing with multi-input and multi-output dynamic systems in continuous and discrete form. CDO version written in FORTRAN IV. VAX version written in FORTRAN 77.

Frisch, Harold P.↗

Analysis of the convergence history of flow through nozzles with shocks

Acceleration techniques such as Wynn's (1986) epsilon algorithm and analysis techiques such as eigensystem analysis are used here to study numerically the convergence properties of an iterative scheme applied to the quasi-one-dimensional Euler and Navier-Stokes equations for flow through nozzles with shocks. The convergence and stability properties are studied by analyzing the dependence of convergence of the code on the discretization technique, boundary conditions, time-step, number of grid points, and the physics of the problem.

Cheer, A. Y.↗

Pulse shaping and extraction of information from ultrasonic reflections in composite materials

The application of ultrasonic inspection to the nondestructive testing of composite materials is examined with reference to specific examples including monitoring of epoxy curing, evaluation of impact damage in a graphite/epoxy laminate, and monitoring of a repaired area in a composite wing plate. Signal processing methods are described which improve time resolution, facilitating identification of discrete echoes corresponding to events in the material.

Kishoni, Doron↗

Delay compensation in integrated communication and control systems. II - Implementation and verification

The implementation and verification of the delay-compensation algorithm are addressed. The delay compensator has been experimentally verified at an IEEE 802.4 network testbed for velocity control of a DC servomotor. The performance of the delay-compensation algorithm was also examined by combined discrete-event and continuous-time simulation of the flight control system of an advanced aircraft that uses the SAE (Society of Automotive Engineers) linear token passing bus for data communications.

Luck, Rogelio↗

Nonclassical aileron buzz in transonic flow

A computational study of inviscid, transonic aileron and trailing-edge buzz instabilities is presented. A mixed Eulerian-Lagrangian formulation is used to model the fluid-structure system and to obtain a system of space-discretized equations that is time-marched to simulate the aeroelastic behavior of the wing-aileron system. Results obtained suggest that shock-induced separation may not be an essential driving force behind all buzz phenomena. Several examples are shown where the shock motion interacts with the aileron motion to extract energy from the flow. If the trailing-edge region is sufficiently flexible and the shocks are at the trailing edge, a trailing-edge buzz instability appears possible.

Bendiksen, Oddvar O.↗

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗

A generic efficient adaptive grid scheme for rocket propulsion modeling

The objective of this research is to develop an efficient, time-accurate numerical algorithm to discretize the Navier-Stokes equations for the predictions of internal one-, two-dimensional and axisymmetric flows. A generic, efficient, elliptic adaptive grid generator is implicitly coupled with the Lower-Upper factorization scheme in the development of ALUNS computer code. The calculations of one-dimensional shock tube wave propagation and two-dimensional shock wave capture, wave-wave interactions, shock wave-boundary interactions show that the developed scheme is stable, accurate and extremely robust. The adaptive grid generator produced a very favorable grid network by a grid speed technique. This generic adaptive grid generator is also applied in the PARC and FDNS codes and the computational results for solid rocket nozzle flowfield and crystal growth modeling by those codes will be presented in the conference, too. This research work is being supported by NASA/MSFC.

Mo, J. D.↗

The onset of Wolf-Rayet wind outflow and the nature of the hot component in the symbiotic nova PU Vulpecula

We have analyzed temporal variations in the far ultraviolet He II (1640), Si IV (1393, 1402), and C IV (1548, 1550) line profiles in eight high dispersion, International Ultraviolet Explorer Short Wavelength Prime spectra of the symbiotic nova PU Vul by comparatively examining these profiles on a common velocity scale. We see clear evidence of the onset of a Wolf-Rayet-like wind outflow from the bloated, contracting white dwarf hot component with terminal velocity of approximately equals -550 to -600 km/s. We have quantitatively analyzed the complicated He II (1640) emission region for the first time and show that the discrete absorption features seen in the He II region occur at precisely the same velocites in each spectrum, thus demonstrating that the absorbing source is steady and not affected by any orbital motion. We demonstrate that there is an underlying He II wind emission feature whose true shape is hidden by superposed absorption due to the foreground red giant wind flowing in front of the white dwarf and abscuring the white dwarf's wind outflow. We present synthetic spectra of He II emission behind an absorbing slab with u = 20 km/s, T = 5000 K, and column densities in the range N = 1 x 10(exp 22) and 1 x 10(exp 23)/sq cm which explain these absorptions. Our analysis of the Si IV and C IV resonance doublets, in velocity space, reveal temporal variations in the profile between 1987 and 1991 with the emergence of clear P Cygni profiles in Si IV by 1990. A nebular emission feature in C III 1909 also appears in the most recent spectra (e.g., SW42538H) while it was absent or extremely weak in the earliest spectra (e.g., SW36332H), thus strengthening evidence that the nebular emission, as seen in permitted and semiforbidden lines, intensities in step with the onset of the hot, fast, wind outflow. We also report the first detection of narrow interstellar (circumbinary shell?) absorption lines near -1 km/s, most strongly in Al III (1854, 1862) and Si IV (1392, 1402). We have carried out a rough quantitative analysis of the He II wind emission by using the theoretical He II Wolf-Rayet profiles of Hamann & Schmutz (1987). We obtain a lower limit to the He II net emission equivalent width of approximately 1 A, a hot component temperature in the range 25,000 to 35,000 K, a hot component radius in the range 5 solar radius to 30 solar radius, a maximum wind velocity of approximately equals -600 km/s and a rough upper limit to the mass-loss rate of less than 1 x 10(exp -5) solar mass/yr. To our knowledge, this is the first quantitative wind analysis, albeit crude, to be carried out for the hot component of a symbiotic nova or symbiotic variable.

Sion, Edward M.↗

Large-eddy simulation of a backward facing step flow using a least-squares spectral element method

We report preliminary results obtained from the large eddy simulation of a backward facing step at a Reynolds number of 5100. The numerical platform is based on a high order Legendre spectral element spatial discretization and a least squares time integration scheme. A non-reflective outflow boundary condition is in place to minimize the effect of downstream influence. Smagorinsky model with Van Driest near wall damping is used for sub-grid scale modeling. Comparisons of mean velocity profiles and wall pressure show good agreement with benchmark data. More studies are needed to evaluate the sensitivity of this method on numerical parameters before it is applied to complex engineering problems.

Chan, Daniel C.↗

Comparison of Several Dissipation Algorithms for Central Difference Schemes

Several algorithms for introducing artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical results are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme.

Swanson, R. C.↗

Improved Results for Route Planning in Stochastic Transportation Networks

In the bus network problem, the goal is to generate a plan for getting from point X to point Y within a city using buses in the smallest expected time. Because bus arrival times are not determined by a fixed schedule but instead may be random. the problem requires more than standard shortest path techniques. In recent work, Datar and Ranade provide algorithms in the case where bus arrivals are assumed to be independent and exponentially distributed. We offer solutions to two important generalizations of the problem, answering open questions posed by Datar and Ranade. First, we provide a polynomial time algorithm for a much wider class of arrival distributions, namely those with increasing failure rate. This class includes not only exponential distributions but also uniform, normal, and gamma distributions. Second, in the case where bus arrival times are independent and geometric discrete random variable,. we provide an algorithm for transportation networks of buses and trains, where trains run according to a fixed schedule.

Boyan, Justin↗

On Some Numerical Dissipation Schemes

Several schemes for introducing an artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar dissipation and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical solutions are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme. The coarse-grid accuracy for the original CUSP scheme is improved by modifying the limiter function used with the scheme, giving comparable accuracy to that obtained with the MATD scheme. The modifications reduce the background dissipation and provide control over the regions where the scheme can become first order.

Swanson, R. C.↗

Oh! I Slipped the Surly Bounds of Earth....and Ran into Space Weather!

Over the past decade the concept of space weather has been introduced and matured in both the scientific community and popular press. Likewise the concept of space climatology recently also is being advanced. Closely linked to these concepts are their impacts on ground- and space-based technological systems; one such system commonly mentioned is manned space flight exemplified by the Space Shuttle and International Space Station (ISS). From a manned space flight perspective, space weather and space climatology have significant effects on the amount of radiation exposure received by humans in space from the ambient high-energy charged particles present in interplanetary space and trapped in the geomagnetosphere. Whereas the impact of space weather for most technological systems is usually discrete and well correlated in time, the principle impact of space weather and space climatology is to increase the probability of latent cancer formation in thetraveler cohort. In this regard, while space weather may be the dominating factor for a given mission, over the life of a long-term program such as the Space Shuttle or ISS space climatology is the controlling factor of latent cancer risk. Human radiation exposure enhancements associated with space weather disturbances has been a concern among scientist and mission controllers since the inception of manned spaceflight nearly forty years ago. This led NASA to develop, in conjunction with the Environmental Science Services Administration s Space Disturbance Forecast Center and the USAF/AWS, the Solar Particle Alert Network (SPAN)-the foundation of an initial U.S. space weather monitoring and forecasting service. Since Apollo, routine space flight operations have evolved to include the use of space weather and climatology data provided through a world-wide network of operational space weather data services to predict and recommend actions to minimize astronaut radiation exposures. NASA Space Radiation Analysis Group (SRAG) flight controllers use real-time space weather data to detect and assess the impact of solar particle events, outer electron belt enhancements, the formation of pseudo-stable additional trapped radiation belts, and the solar cycle modulation of trapped radiation belts and galactic cosmic rays. Energetic particle data from GOES spacecraft are automatically ingested from NOAA Space Environment Center data servers and used to drive a model for the estimating the exposure to astronauts from solar particle events. While adequate for current manned space flight support, the existing operational space weather support system requires improvements to address the anticipated evolution in both the character of manned missions as well as space flight operations management. Necessary space weather data improvements include: reliably available (near) real-time space weather data on a fixed schedule via redundant access methods that support autonomous data acquisition by computer systems behind enterprise firewalls; and rapid transition of promising research sensors into operational systems.

Golightly, Michael J.↗

Rover Graphical Simulator

Rover Graphical Simulator (RGS) is a package of software that generates images of the motion of a wheeled robotic exploratory vehicle (rover) across terrain that includes obstacles and regions of varying traversability. The simulated rover moves autonomously, utilizing reasoning and decision-making capabilities of a fuzzy-logic navigation strategy to choose its path from an initial to a final state. RGS provides a graphical user interface for control and monitoring of simulations. The numerically simulated motion is represented as discrete steps with a constant time interval between updates. At each simulation step, a dot is placed at the old rover position and a graphical symbol representing the rover is redrawn at the new, updated position. The effect is to leave a trail of dots depicting the path traversed by the rover, the distances between dots being proportional to the local speed. Obstacles and regions of low traversability are depicted as filled circles, with buffer zones around them indicated by enclosing circles. The simulated robot is equipped with onboard sensors that can detect regional terrain traversability and local obstacles out to specified ranges. RGS won the NASA Group Achievement Award in 2002.

Bon, Bruce↗

Evaluating the Application of Multi-Satellite Observations in Hydrologic Modeling

When monitoring local or regional hydrosphere dynamics for applications such as agricultural productivity or drought and flooding events, it is necessary to have accurate, high-resolution estimates of terrestrial water and energy storages. Though in-situ observations provide reliable estimates of hydrologic states and fluxes, they are only capable of accurately capturing the dynamics at relatively discrete points in space and time, which makes them inadequate for characterizing the variability of the water budget across scales. In contrast, satellite-based remote sensing is ideal for providing observations of hydrological states and fluxes because it provides spatially-distributed observations at spatial and temporal scales required for regional land surface process modeling. Due to the continued progress in algorithm development and emerging satellite technology, we now have near-real time monitoring of several components of the water cycle including precipitation, soil moisture, lake and river height, terrestrial water storage, snow cover, and evapotranspiration. As these data become more readily available, their application to hydrologic modeling is becoming more common, however there remains little consensus on the most appropriate method for optimal integration and evaluation in regard to hydrological applications. Here we present two case studies operationally applying several remotely sensed products from AMSR-E, GRACE, and MODIS and discuss assimilation strategies, ease of integration and interpretation, and methods for quantifying the success of the application methodology.

Bolten, John↗

SpF: Enabling Petascale Performance for Pseudospectral Dynamo Models

Pseudospectral (PS) methods possess a number of characteristics (e.g., efficiency, accuracy, natural boundary conditions) that are extremely desirable for dynamo models. Unfortunately, dynamo models based upon PS methods face a number of daunting challenges, which include exposing additional parallelism, leveraging hardware accelerators, exploiting hybrid parallelism, and improving the scalability of global memory transposes. Although these issues are a concern for most models, solutions for PS methods tend to require far more pervasive changes to underlying data and control structures. Further, improvements in performance in one model are difficult to transfer to other models, resulting in significant duplication of effort across the research community.We have developed an extensible software framework for pseudospectral methods called SpF that is intended to enable extreme scalability and optimal performance. High-level abstractions provided by SpF unburden applications of the responsibility of managing domain decomposition and load balance while reducing the changes in code required to adapt to new computing architectures. The key design concept in SpF is that each phase of the numerical calculation is partitioned into disjoint numerical kernels that can be performed entirely in-processor. The granularity of domain-decomposition provided by SpF is only constrained by the data-locality requirements of these kernels. SpF builds on top of optimized vendor libraries for common numerical operations such as transforms, matrix solvers, etc., but can also be configured to use open source alternatives for portability. SpF includes several alternative schemes for global data redistribution and is expected to serve as an ideal testbed for further research into optimal approaches for different network architectures.In this presentation, we will describe the basic architecture of SpF as well as preliminary performance data and experience with adapting legacy dynamo codes. We will conclude with a discussion of planned extensions to SpF that will provide pseudospectral applications with additional flexibility with regard to time integration, linear solvers, and discretization in the radial direction.

Pseudospectral (PS)↗

Applicability of Spectral Decomposition by Varimax-Rotated, Principal Component Analysis to the Surface Biology and Geology (SBG) VNIR Mission Concept

Cyanobacterial and Harmful Algal Blooms (CyanoHABs) are a growing concern in coastal and inland waters. But, spectral interference from multiple constituents in optically complex waters can hamper application of remote sensing using traditional image processing methods. The Kent State University (KSU) spectral decomposition method can be applied to multispectral and hyperspectral remote sensing images (e.g. HICO and the NASA Glenn HSI2) to partition and identify signals related to cyanobacteria, algae, pigment degradation products and suspended sediment in each pixel. Fundamental to the use of remote sensing data is the ability to extract independent signals from correlated hyperspectral VNIR data cubes. The Kent State University varimax-rotated, principal component analysis method (VPCA) is important to integrate into the SBG VNIR mission concept because it provides greater specificity, a software-based SNR boost relative to hardware performance, and can assist with Cal/Val, Modeling and Applications. We present examples of the hyperspectral application of the KSU VPCA method with relevance to SBG. The information extracted by VPCA can be validated spectrally or spatially with laboratory and/or in situ sensors, which capture spatial or time series of information at discrete points within remote sensing images. Comparisons show hyperspectral sensors extract more components than multispectral ones, but more independent information can be extracted from multispectral sensors by VPCA than traditional band ratio approaches. The spectral decomposition method is capable of enhancing the signal to noise ratio (SNR) of the NASA Glenn, second-generation hyperspectral imager by a factor of 7x to 20x, with a spectral reproducibility of ±3%. The spectral decomposition method, when compared against existing remote sensing monitoring methods exhibits both greater specificity and a lower detection limit. The method has been validated with multispectral images in Lake Erie to quantify the Microcystis CyanoHAB and from the Indian River Lagoon, Florida to quantify the Brown Tide resulting from A. lagunesnsis. Field operations in the Western Basin of Lake Erie were conducted using a bbe Fluoroprobe to collect vertical profiles and horizontal tows along a transect from the Toledo to the Detroit Lighthouse during coincident satellite overpasses. Extraction of pixel values from the MODIS Aqua sensor yields agreement between in situ field and lab-based measures of cyanobacterial, cryptophyte, diatoms and green algae, suspended sediment and pigment degradation products with R2>0.8.

hyperspectral↗