Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local time stepping”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case

Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.

Goodrich, John W.↗
NASA NTRS1994

Large-eddy simulation of a plane wake

Previously the theoretical development leading to the dynamic localization model (DLM) for large-eddy simulation (LES) was presented. The method has been successfully applied to isotropic turbulence, channel flow, and the flow over a backward-facing step. Here we apply the model to the computation of the temporally developing place wake. The two main objectives of this project are: (1) Use the model to perform an LES of a time developing plane wake and compare the results with direction numerical simulation (DNS) data to see if important statistical measures can be readily predicted, and to provide a relative evaluation of the several versions of the model in terms of predictive capability and cost; and (2) If the tests in (1) show that the model generates reliable predictions, then use the LES to study various aspects of the physics of turbulent wakes and mixing layers.

Ghosal, Sandip↗
NASA NTRS2000

Plasma Sheet Source and Loss Processes

Data from the TIMAS ion mass spectrometer on the Polar satellite, covering 15 ev/e to 33 keV/e in energy and essentially 4(pi) in view angles, are used to investigate the properties of earthward (sunward) field-aligned flows of ions, especially protons, in the plasma sheet-lobe transition region near local midnight. A total of 142 crossings of this region are analyzed at 12-sec time resolution, all in the northern hemisphere, at R(SM) approx. 4 - 7 R(sub E), and most (106) in the poleward (sunward) direction. Earthward proton flows are prominent in this transition region (greater than 50% of the time), typically appearing as sudden "blasts" with the most energetic protons (approx. 33 keV) arriving first with weak flux, followed by protons of decreasing energy and increasing flux until either: (1) a new "blast" appears, (2) the flux ends at a sharp boundary, or (3) the flux fades away within a few minutes as the mean energy drops to a few keV. Frequent step-like changes (less than 12 sec) of the flux suggest that perpendicular gradients on the scale of proton gyroradii are common. Peak flux is similar to central plasma sheet proton flux (10(exp 5) - 10(exp 6)/[cq cm sr sec keV/e] and usually occurs at E approx. 4 - 12 keV. Only the initial phase of each "blast" (approx. 1 min) displays pronounced field-alignment of the proton velocity distribution, consistent with the time-of-flight separation of a more or less isotropic source distribution with df/d(nu) less than 0. The dispersive signatures are often consistent with a source at R(SM) less than or equal to 30 R(sub E). No systematic latitudinal velocity dispersion is found, implying that the equatorial plasma source is itself convecting. In short, the proton "blasts" appear as sudden local expansions of central plasma sheet particles along reconfigured ("dipolarized") magnetic field lines.

Lennartsson, O. W.↗
NASA NTRS1980

A theory for the origin of a self-replicating chemical system. I - Natural selection of the autogen from short, random oligomers

A general theory is presented for the origin of a self-replicating chemical system, termed an autogen, which is capable of both crude replication and translation (protein synthesis). The theory requires the availability of free energy and monomers to the system, a significant background low-yield synthesis of kinetically stable oligopeptides and oligonucleotides, the localization of the oligomers, crude oligonucleotide selectivity of amino acids during oligopeptide synthesis, crude oligonucleotide replication, and two short peptide families which catalyze replication and translation, to produce a localized group of at least one copy each of two protogenes and two protoenzymes. The model posits a process of random oligomerization, followed by the random nucleation of functional components and the rapid autocatalytic growth of the functioning autogen to macroscopic amounts, to account for the origin of the first self-replicating system. Such a process contains steps of such high probability and short time periods that it is suggested that the emergence of an autogen in a laboratory experiment of reasonable time scale may be possible.

White, D. H.↗
NASA NTRS2002

Procedure for Tooth Contact Analysis of a Face Gear Meshing With a Spur Gear Using Finite Element Analysis

A procedure was developed to perform tooth contact analysis between a face gear meshing with a spur pinion using finite element analysis. The face gear surface points from a previous analysis were used to create a connected tooth solid model without gaps or overlaps. The face gear surface points were used to create a five tooth face gear Patran model (with rim) using Patran PCL commands. These commands were saved in a series of session files suitable for Patran input. A four tooth spur gear that meshes with the face gear was designed and constructed with Patran PCL commands. These commands were also saved in a session files suitable for Patran input. The orientation of the spur gear required for meshing with the face gear was determined. The required rotations and translations are described and built into the session file for the spur gear. The Abaqus commands for three-dimensional meshing were determined and verified for a simplified model containing one spur tooth and one face gear tooth. The boundary conditions, loads, and weak spring constraints were determined to make the simplified model work. The load steps and load increments to establish contact and obtain a realistic load was determined for the simplified two tooth model. Contact patterns give some insight into required mesh density. Building the two gears in two different local coordinate systems and rotating the local coordinate systems was verified as an easy way to roll the gearset through mesh. Due to limitation of swap space, disk space and time constraints of the summer period, the larger model was not completed.

Bibel, George↗
NASA NTRS2017

Low-Thrust Trajectory Optimization with Simplified SQP Algorithm

The problem of low-thrust trajectory optimization in highly perturbed dynamics is a stressing case for many optimization tools. Highly nonlinear dynamics and continuous thrust are each, separately, non-trivial problems in the field of optimal control, and when combined, the problem is even more difficult. This paper de-scribes a fast, robust method to design a trajectory in the CRTBP (circular restricted three body problem), beginning with no or very little knowledge of the system. The approach is inspired by the SQP (sequential quadratic programming) algorithm, in which a general nonlinear programming problem is solved via a sequence of quadratic problems. A few key simplifications make the algorithm presented fast and robust to initial guess: a quadratic cost function, neglecting the line search step when the solution is known to be far away, judicious use of end-point constraints, and mesh refinement on multiple shooting with fixed-step integration.In comparison to the traditional approach of plugging the problem into a “black-box” NLP solver, the methods shown converge even when given no knowledge of the solution at all. It was found that the only piece of information that the user needs to provide is a rough guess for the time of flight, as the transfer time guess will dictate which set of local solutions the algorithm could converge on. This robustness to initial guess is a compelling feature, as three-body orbit transfers are challenging to design with intuition alone. Of course, if a high-quality initial guess is available, the methods shown are still valid.We have shown that endpoints can be efficiently constrained to lie on 3-body repeating orbits, and that time of flight can be optimized as well. When optimizing the endpoints, we must make a trade between converging quickly on sub-optimal endpoints or converging more slowly on end-points that are arbitrarily close to optimal. It is easy for the mission design engineer to adjust this trade based on the problem at hand.The biggest limitation to the algorithm at this point is that multi-revolution transfers (greater than 2 revolutions) do not work nearly as well. This restriction comes in because the relationship between node 1 and node N becomes increasingly nonlinear as the angular distance grows. Trans-fers with more than about 1.5 complete revolutions generally require the line search to improve convergence. Future work includes: Comparison of this algorithm with other established tools; improvements to how multiple-revolution transfers are handled; parallelization of the Jacobian computation; in-creased efficiency for the line search; and optimization of many more trajectories between a variety of 3-body orbits.

Parrish, Nathan L.↗
NASA NTRS1974

The effect of yttrium and thorium on the oxidation behavior of Ni-Cr-Al alloys

The effect of quaternary additions of 0.5% Y, 0.5 and 1.0% Th to a base alloy of Ni-10CR-5Al on the oxidation behavior and mechanism was studied during oxidation in air over the range of 1000 to 1200 C. The presence of yttrium decreased the oxidation kinetics slightly, whereas, the addition of thorium caused a slight increase. Oxide scale adherence was markedly improved by the addition of the quaternary elements. Although a number of oxides formed on yttrium containing alloys, quantitative X-ray diffraction clearly showed that the rate-controlling step was the diffusion of aluminum through short circuit paths in a thin layer of alumina that formed parabolically with time. Although the scale adherence of the yttrium containing alloy was considerably better than the base alloys, spalling did occur that was attributed to the formation of the voluminous YAG particles which grew in a mushroom-like manner, lifting the protective scale off the subrate locally. The YAG particles formed primarily at grain boundaries in the substrate in which the yttrium originally existed as YNi9.

Kumar, A.↗
NASA NTRS1993

Monoclonal Antibody Testing for Cancer Metastasis

Malignant cells are characterized by the ability to invade surrounding normal tissues. Tumor invasion is abetted by proteolytic enzymes that have been correlated with recurrent disease and metastasis. These enzymes are involved in a cascade of proteolytic interactions with other enzymes and inhibitors which allow cancer cells to dissolve surrounding extracellular matrix, thereby enabling the cells to rapidly invade adjacent tissues and migrate to metastatic sites distant from the primary tumor. Among these proteases are the plasminogen activators (PA), collagenase IV, faminase, and in some cases cathepsin D, which together mediate key steps in the invasion process of metastasis. Cells which have the selective advantage for invasion and metastasis are those capable of regulating their proteolytic activity and proliferation. Cells in the process of invasion would be probably down-regulated for proliferation, but subsequent to attachment and adhesion at a distant site, would then be in a proliferative mode, up-regulating DNA replication. Urokinase (uPA) can be present in the tissues in several molecular forms. The inactive proenzyme is a single chain protein (scuPA) that is cleaved at Lys. 158 to form the double chain, high molecular weight active form (HMW-uPA) of 54 kD. A low molecular weight form (LMW-uPA) can also be produced by cleavage of the HMW-U PA at Lys. 135 - Lys. 136 giving a 35 kD active enzyme. Recently, it has been shown that the HMW active form of urokinase, bound to the tumor cell membrane, is responsible for the local lysis of the extracellular matrix, hence the tissue invasion mechanism for metastasis (Andreasen et al, 19861. Receptor- (membrane) bound uPA is twice as efficient (catalytically) as free fluid-phase uPA. Tho unbound uPA and the LMW form is not responsible for most of the local dissolution of extracellular matrix in the immediate vicinity of the metastatic tumor cell. High levels of urokinase (greater than 3.49 ng/mg of total protein) extracted from breast tumor tissues have recently been shown, together with plasminogen activator inhibitor 1 (PAI-1), to be a good prognostic indicator for high risk of recurrence and shorter patient survival times. In this project, we have attempted to develop immunocytochemical methodologies for the clinical assessment of the expression of urokinase plasminogen activator, which has been implicated to be important for initial steps in tumor invasion, and to relate it to cell proliferation and DNA replication at the single-cell level.

Source record↗
NASA NTRS1992

Divergence preserving discrete surface integral methods for Maxwell's curl equations using non-orthogonal unstructured grids

Several new discrete surface integral (DSI) methods for solving Maxwell's equations in the time-domain are presented. These methods, which allow the use of general nonorthogonal mixed-polyhedral unstructured grids, are direct generalizations of the canonical staggered-grid finite difference method. These methods are conservative in that they locally preserve divergence or charge. Employing mixed polyhedral cells, (hexahedral, tetrahedral, etc.) these methods allow more accurate modeling of non-rectangular structures and objects because the traditional stair-stepped boundary approximations associated with the orthogonal grid based finite difference methods can be avoided. Numerical results demonstrating the accuracy of these new methods are presented.

Madsen, Niel K.↗
NASA NTRS1989

Thickening oscillation of a delta wing using Navier-Stokes and Navier-displacement equations

The problem of unsteady, supersonic, locally conical, vortical flow around a delta wing undergoing thickening oscillation is solved using the unsteady, thin-layer Navier-Stokes equations and the unsteady linearized Navier-displacement equations. The unsteady, thin-layer Navier-Stokes equations are solved using the implicit approximate-factorization finite-volume scheme to compute the conservative components of the flow vector field. With the conservative components known at any time step, the linearized, Navier-displacement equations are solved using the alternating, direction-implicit scheme to obtain the grid points displacements due to known displacement boundary conditions. A grid-displacements limiter, in the form of a low mesh Reynolds number, is used to limit grid-folding in regions of highly reversed flow.

Chuang, Hsin-Kung A.↗
NASA NTRS1986

Mass and energy balance of the cold Io torus

A new model of the cold Io torus is described. Ions and energy are injected into the system by independent processes so that the mass balance is isolated from the energy balance. The primary source of energy is local ionization and acceleration of hot pickup ions resulting from charge exchange between thermal ions and an extended cloud of Iogenic sulfur and oxygen atoms. The primary energy loss mechanism of the plasma is collisionally excited line emission at optical wavelengths. The primary ion source is radial diffusion inward from the hot torus on a time scale of 140-710 days. The primary ion loss mechanism is a novel two-step enhanced recombination mechanism involving charge exchange between thermal ions and an extended cloud of neutral SO2 molecules, followed by rapid dissociative recombination of the resultant molecular ion. The model provides a self-consistent solution which reconciles a number of diverse observations with known physical processes.

Moreno, M. A.↗
NASA NTRSDate not supplied

Stability Analysis of the Flow over a Swept Forward-Facing Step using PIV Base Flowsin a Non-Orthogonal Coordinate System

Understanding the flow of a crossflow-vortex-dominated boundary layer over a forward-facing step excrescence is necessary in order to mitigate an early transition scenario on the wings and tail of commercial aircraft. By performing BiGlobal stability analysis on the flowfield measured with high-resolution, stereographic Particle Image Velocimetry (PIV), previous work has shown that a family of unstable disturbances exists in the direct downstream vicinity of a supercritical, forward-facing step. The goal of the present paper is to improve the previously used stability approach on two fronts: 1) the stability problem is formulated in a non-orthogonal coordinate system and 2) the contribution of the out-of-plane base-flow derivatives is accounted for locally. The resulting eigen solutions feature significantly stronger growth rates and eigenfunctions that are localized above regions of reverse flow. The changes in the problem formulation furthermore establish a significant improvement in the comparison of the stability solutions with the results from the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content. A large destabilizing effect by the out-of-plane base-flow derivatives is determined to be justified, despite the assumption that these derivatives are small, because the effect can be reconstructed upon using an eigenvalue-correction formula that assumes the responsible terms are infinitesimally small. The unstable perturbation mechanisms are demonstrated to have a convective nature by assessing the relation between the frequency and the out-of-plane wavenumber, which indicates that their group speed does not approach zero. The last 2 facts, that 1) the perturbations are convective and 2) that the out-of-plane base-flow-derivative terms are small, remove any qualitative suspicions that the perturbation problem is not conducive to parabolization. This opens the path to analyzing this problem with a plane-marching approach.

Boundary- Layer Transition↗
NASA NTRS1989

Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a tubine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of independence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

Choo, Yung K.↗
NASA NTRS1989

Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a turbine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of dependence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

Choo, Yung K.↗
NASA NTRS1998

Interfacing with Legacy using Remote Method Invocation

The assignment described was enough to make a neophyte Java developer bolt for the door: provide a remote method for use by an applet which invokes a native method that wraps a function in an existing legacy library. The purpose of the remote method is to return an instance of a class object whose contents reflect the data structure returned by the legacy function. While embroiled in implementation, I would have spent the time wading through their JNI use group archive as well, but I couldn't seem to locate one. Subsequently, I made the decision to try to document my findings in order to assist others. Before we start on the class design, let's look at what the existing legacy code does. The C function to be called, Get-Legacy-Data, consists of two steps: an ASII file is read from the local disk and its contents are parsed into a Legacy_Type structure whose address is passed as an argument by the caller. The legacy code was compiled into a shared object library, legacy. so, using the IRIX 6.2 compiler and then loaded onto the Web server, a Silicon Graphics Indy station loaded with the IRIX 6.4 operating system. As far as the class design is concerned, the first thing required is a class to act as a template for the data structure returned by the legacy function. This class, JLegacy, declares a series of public instance variables which correspond to the members of Legacy_Type and provides a parameterless constructor. This constructor is never called, not even by the native method which allocates the object for return to the remote method. Next, the remote interface declaration for the remote object must be defined. In order for JLegacyRO to implement getJLegacy, JLegacyRO must interface with the existing legacy code through a native method, getn. getn is declared in the JLegacyRO class but implemented in C, just like the legacy code. getn returns a JLegacy instance and is declared static since its implementation is the same for all instances of the JLegacyRO class.

Howard, Scott M.↗
NASA NTRS1998

Impact of the 1997-1998 El-Nino of Regional Hydrology

The 1997-1998 El-Nino brought with it a range of severe local-regional hydrological phenomena. Record high temperatures and extremely dry soil conditions in Texas is an example of this regional effect. The El-Nino and La-Nina change the continental weather patterns considerably. However, connections between continental weather anomalies and regional or local anomalies have not been established to a high degree of confidence. There are several unique features of the recent El-Nino and La-Nina. Due to the recognition of the present El-Nino well in advance, there have been several coupled model studies on global and regional scales. Secondly, there is a near real-time monitoring of the situation using data from satellite sensors, namely, SeaWIFS, TOVS, AVHRR and GOES. Both observations and modeling characterize the large scale features of this El-Nino fairly well. However the connection to the local and regional hydrological phenomenon still needs to be made. This paper will use satellite observations and analysis data to establish a relation between local hydrology and large scale weather patterns. This will be the first step in using satellite data to perform regional hydrological simulations of surface temperature and soil moisture.

Lakshmi, Venkataraman↗
NASA NTRS2011

Calibrated Noise Measurements with Induced Receiver Gain Fluctuations

The lack of well-developed techniques for modeling changing statistical moments in our observations has stymied the application of stochastic process theory in science and engineering. These limitations were encountered when modeling the performance of radiometer calibration architectures and algorithms in the presence of non stationary receiver fluctuations. Analyses of measured signals have traditionally been limited to a single measurement series. Whereas in a radiometer that samples a set of noise references, the data collection can be treated as an ensemble set of measurements of the receiver state. Noise Assisted Data Analysis is a growing field of study with significant potential for aiding the understanding and modeling of non stationary processes. Typically, NADA entails adding noise to a signal to produce an ensemble set on which statistical analysis is performed. Alternatively as in radiometric measurements, mixing a signal with calibrated noise provides, through the calibration process, the means to detect deviations from the stationary assumption and thereby a measurement tool to characterize the signal's non stationary properties. Data sets comprised of calibrated noise measurements have been limited to those collected with naturally occurring fluctuations in the radiometer receiver. To examine the application of NADA using calibrated noise, a Receiver Gain Modulation Circuit (RGMC) was designed and built to modulate the gain of a radiometer receiver using an external signal. In 2010, an RGMC was installed and operated at the National Institute of Standards and Techniques (NIST) using their Noise Figure Radiometer (NFRad) and national standard noise references. The data collected is the first known set of calibrated noise measurements from a receiver with an externally modulated gain. As an initial step, sinusoidal and step-function signals were used to modulate the receiver gain, to evaluate the circuit characteristics and to study the performance of a variety of calibration algorithms. The receiver noise temperature and time-bandwidth product of the NFRad are calculated from the data. Statistical analysis using temporal-dependent calibration algorithms reveals that the natural occurring fluctuations in the receiver are stationary over long intervals (100s of seconds); however the receiver exhibits local non stationarity over the interval over which one set of reference measurements are collected. A variety of calibration algorithms have been applied to the data to assess algorithms' performance with the gain fluctuation signals. This presentation will describe the RGMC, experiment design and a comparative analysis of calibration algorithms.

Racette, Paul↗
NASA NTRS2004

Time-Extended Policies in Mult-Agent Reinforcement Learning

Reinforcement learning methods perform well in many domains where a single agent needs to take a sequence of actions to perform a task. These methods use sequences of single-time-step rewards to create a policy that tries to maximize a time-extended utility, which is a (possibly discounted) sum of these rewards. In this paper we build on our previous work showing how these methods can be extended to a multi-agent environment where each agent creates its own policy that works towards maximizing a time-extended global utility over all agents actions. We show improved methods for creating time-extended utilities for the agents that are both "aligned" with the global utility and "learnable." We then show how to crate single-time-step rewards while avoiding the pi fall of having rewards aligned with the global reward leading to utilities not aligned with the global utility. Finally, we apply these reward functions to the multi-agent Gridworld problem. We explicitly quantify a utility's learnability and alignment, and show that reinforcement learning agents using the prescribed reward functions successfully tradeoff learnability and alignment. As a result they outperform both global (e.g., team games ) and local (e.g., "perfectly learnable" ) reinforcement learning solutions by as much as an order of magnitude.

Tumer, Kagan↗