Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

ENSO Diversity Simulated in a Revised Cane-Zebiak Model

The El Niño-Southern Oscillation (ENSO) phenomenon features rich sea surface temperature (SST) spatial pattern variations dominated by the Central Pacific (CP) and Eastern Pacific (EP) patterns during its warm phase. Understanding such ENSO pattern diversity has been a subject under extensive research activity. To provide a framework for unveiling the fundamental dynamics of ENSO diversity, an intermediate coupled model based on the Cane-Zebiak-type framework, named RCZ, is established in this study. Compared with the original Cane-Zebiak model, RCZ consists of revised model formulation and well-tuned parameterization schemes. All model components are carefully validated against the observations via the standalone mode, in which the observed anomalous SST (wind stress) forcing is prescribed to drive the atmospheric (oceanic) component. The superiority of RCZ’s model components over those in the original Cane-Zebiak model is evidenced by their better performance in simulating the observations. Coupled simulation with RCZ satisfactorily reproduces aspects of the observed ENSO characteristics, including the spatial pattern, phase-locking, amplitude asymmetry, and, particularly, ENSO diversity/bi-modality. RCZ serves as a promising tool for studying dynamics of ENSO diversity as it resolves most of the relevant processes proposed in the literature, including atmospheric nonlinear convective heating, oceanic nonlinear dynamical heating, and the ENSO/westerly wind burst interaction.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1972

Temperature compensation of light emitting diodes and photodiodes.

Improved temperature compensation of the output from a GaAs LED can be obtained by connecting the diode in series with the emitter of an n-p-n transistor whose base circuit incorporates a high-gain operational amplifier and a parallel thermistor-resistor network. The temperature dependence of several photodiodes is also reported along with a practical temperature compensation scheme. Finally, the temperature compensation of a photon coupled pair consisting of a GaAs LED and a silicon photodiode is described. The combination provides a stability, linearity, and repeatability of about 1% over the upper 60% of the LED forward current range in the temperature interval from 15 to 45 C.

Burrous, C. N.↗
NASA NTRS1979

Space Shuttle status

Hardware considerations and the orbital flight test program of the Space Shuttle are discussed, including reference to the Space Shuttle main engine, payload bay doors, TV cameras for flight inspection, and launch crew escape schemes. Particular attention is given to the thermal protection system consisting of 32,000 tiles composed of approximately 3% pure silica fibers and 97% air, which has a design requirement of keeping the basic aluminum skin below 350 F for a surface temperature of 2700 F.

Slayton, D. K.↗
NASA NTRS1981

E-beam generated holographic masks for optical vector-matrix multiplication

An optical vector matrix multiplication scheme that encodes the matrix elements as a holographic mask consisting of linear diffraction gratings is proposed. The binary, chrome on glass masks are fabricated by e-beam lithography. This approach results in a fairly simple optical system that promises both large numerical range and high accuracy. A partitioned computer generated hologram mask was fabricated and tested. This hologram was diagonally separated outputs, compact facets and symmetry about the axis. The resultant diffraction pattern at the output plane is shown. Since the grating fringes are written at 45 deg relative to the facet boundaries, the many on-axis sidelobes from each output are seen to be diagonally separated from the adjacent output signals.

Arnold, S. M.↗
NASA NTRS1990

Compiler-directed cache management in multiprocessors

The necessity of finding alternatives to hardware-based cache coherence strategies for large-scale multiprocessor systems is discussed. Three different software-based strategies sharing the same goals and general approach are presented. They consist of a simple invalidation approach, a fast selective invalidation scheme, and a version control scheme. The strategies are suitable for shared-memory multiprocessor systems with interconnection networks and a large number of processors. Results of trace-driven simulations conducted on numerical benchmark routines to compare the performance of the three schemes are presented.

Cheong, Hoichi↗
NASA NTRS1991

Three dimensional unstructured multigrid for the Euler equations

The three-dimensional Euler equations are solved on unstructured tetrahedral meshes using a multigrid strategy. The driving algorithm consists of an explicit vertex-based finite-element scheme, which employs an edge-based data-structure to assemble the residuals. The multigrid approach employs a sequence of independently generated coarse and fine meshes to accelerate the convergence to steady-state of the fine grid solution. Variables, residuals and corrections are passed back and forth between the various grids of the sequence using linear interpolation. The addresses and weights for interpolation are determined in a preprocessing stage using an efficient graph traversal algorithm. The preprocessing operation is shown to require a negligible fraction of the CPU time required by the overall solution procedure, while gains in overall solution efficiencies greater than an order of magnitude are demonstrated on meshes containing up to 350,000 vertices. Solutions using globally regenerated fine meshes as well as adaptively refined meshes are given.

Mavriplis, D. J.↗
NASA NTRS2006

An Integrated Snow Radiance and Snow Physics Modeling Framework for Cold Land Surface Modeling

Recent developments in forward radiative transfer modeling and physical land surface modeling are converging to allow the assembly of an integrated snow/cold lands modeling framework for land surface modeling and data assimilation applications. The key elements of this framework include: a forward radiative transfer model (FRTM) for snow, a snowpack physical model, a land surface water/energy cycle model, and a data assimilation scheme. Together these form a flexible framework for self-consistent remote sensing and water/energy cycle studies. In this paper we will describe the elements and the integration plan. Each element of this framework is modular so the choice of element can be tailored to match the emphasis of a particular study. For example, within our framework, four choices of a FRTM are available to simulate the brightness temperature of snow: Two models are available to model the physical evolution of the snowpack and underlying soil, and two models are available to handle the water/energy balance at the land surface. Since the framework is modular, other models-physical or statistical--can be accommodated, too. All modules will operate within the framework of the Land Information System (LIS), a land surface modeling framework with data assimilation capabilities running on a parallel-node computing cluster at the NASA Goddard Space Flight Center. The advantages of such an integrated modular framework built on the LIS will be described through examples-e.g., studies to analyze snow field experiment observations, and simulations of future satellite missions for snow and cold land processes.

Kim, Edward J.↗
NASA NTRS2007

Smooth Phase Interpolated Keying

Smooth phase interpolated keying (SPIK) is an improved method of computing smooth phase-modulation waveforms for radio communication systems that convey digital information. SPIK is applicable to a variety of phase-shift-keying (PSK) modulation schemes, including quaternary PSK (QPSK), octonary PSK (8PSK), and 16PSK. In comparison with a related prior method, SPIK offers advantages of better performance and less complexity of implementation. In a PSK scheme, the underlying information waveform that one seeks to convey consists of discrete rectangular steps, but the spectral width of such a waveform is excessive for practical radio communication. Therefore, the problem is to smooth the step phase waveform in such a manner as to maintain power and bandwidth efficiency without incurring an unacceptably large error rate and without introducing undesired variations in the amplitude of the affected radio signal. Although the ideal constellation of PSK phasor points does not cause amplitude variations, filtering of the modulation waveform (in which, typically, a rectangular pulse is converted to a square-root raised cosine pulse) causes amplitude fluctuations. If a power-efficient nonlinear amplifier is used in the radio communication system, the fluctuating-amplitude signal can undergo significant spectral regrowth, thus compromising the bandwidth efficiency of the system. In the related prior method, one seeks to solve the problem in a procedure that comprises two major steps: phase-value generation and phase interpolation. SPIK follows the two-step approach of the related prior method, but the details of the steps are different. In the phase-value-generation step, the phase values of symbols in the PSK constellation are determined by a phase function that is said to be maximally smooth and that is chosen to minimize the spectral spread of the modulated signal. In this step, the constellation is divided into two groups by assigning, to information symbols, phase values that result in equal numbers of clockwise and counter-clockwise phase rotations for equally likely symbols. The purpose served by assigning phase values in this way is to prevent unnecessary generation of spectral lines and prevent net shifts of the carrier signal. In the phase-interpolation step, the smooth phase values are interpolated over a number, n, of consecutive symbols (including the present symbol) by means of an unconventional spline curve fit.

Borah, Deva K.↗
NASA NTRS2010

Development and Application of Agglomerated Multigrid Methods for Complex Geometries

We report progress in the development of agglomerated multigrid techniques for fully un- structured grids in three dimensions, building upon two previous studies focused on efficiently solving a model diffusion equation. We demonstrate a robust fully-coarsened agglomerated multigrid technique for 3D complex geometries, incorporating the following key developments: consistent and stable coarse-grid discretizations, a hierarchical agglomeration scheme, and line-agglomeration/relaxation using prismatic-cell discretizations in the highly-stretched grid regions. A signi cant speed-up in computer time is demonstrated for a model diffusion problem, the Euler equations, and the Reynolds-averaged Navier-Stokes equations for 3D realistic complex geometries.

Nishikawa, Hiroaki↗
NASA NTRS2004

Increasing Accuracy in Computed Inviscid Boundary Conditions

A technique has been devised to increase the accuracy of computational simulations of flows of inviscid fluids by increasing the accuracy with which surface boundary conditions are represented. This technique is expected to be especially beneficial for computational aeroacoustics, wherein it enables proper accounting, not only for acoustic waves, but also for vorticity and entropy waves, at surfaces. Heretofore, inviscid nonlinear surface boundary conditions have been limited to third-order accuracy in time for stationary surfaces and to first-order accuracy in time for moving surfaces. For steady-state calculations, it may be possible to achieve higher accuracy in space, but high accuracy in time is needed for efficient simulation of multiscale unsteady flow phenomena. The present technique is the first surface treatment that provides the needed high accuracy through proper accounting of higher-order time derivatives. The present technique is founded on a method known in art as the Hermitian modified solution approximation (MESA) scheme. This is because high time accuracy at a surface depends upon, among other things, correction of the spatial cross-derivatives of flow variables, and many of these cross-derivatives are included explicitly on the computational grid in the MESA scheme. (Alternatively, a related method other than the MESA scheme could be used, as long as the method involves consistent application of the effects of the cross-derivatives.) While the mathematical derivation of the present technique is too lengthy and complex to fit within the space available for this article, the technique itself can be characterized in relatively simple terms: The technique involves correction of surface-normal spatial pressure derivatives at a boundary surface to satisfy the governing equations and the boundary conditions and thereby achieve arbitrarily high orders of time accuracy in special cases. The boundary conditions can now include a potentially infinite number of time derivatives of surface-normal velocity (consistent with no flow through the boundary) up to arbitrarily high order. The corrections for the first-order spatial derivatives of pressure are calculated by use of the first-order time derivative velocity. The corrected first-order spatial derivatives are used to calculate the second- order time derivatives of velocity, which, in turn, are used to calculate the corrections for the second-order pressure derivatives. The process as described is repeated, progressing through increasing orders of derivatives, until the desired accuracy is attained.

Dyson, Roger↗
DOE OSTI2022

A simple improved low temperature correction for the hierarchical equations of motion

The study of open system quantum dynamics has been transformed by the hierarchical equations of motion (HEOM) method, which gives the exact dynamics for a system coupled to a harmonic bath at arbitrary temperature and system–bath coupling strength. However, in its standard form, this method is only consistent with the weak-coupling quantum master equation at all temperatures when many auxiliary density operators are included in the hierarchy, even when low temperature corrections are included. Here, we propose a new low temperature correction scheme for the termination of the hierarchy based on Zwanzig projection, which alleviates this problem and restores consistency with the weak-coupling master equation with a minimal hierarchy. The utility of the new correction scheme is demonstrated on a range of model systems, including the Fenna–Matthews–Olson complex. The new closure is found to improve convergence of the HEOM even beyond the weak-coupling limit and is very straightforward to implement in existing HEOM codes.

Fay, Thomas P. (ORCID:000000030625731X)↗
NASA NTRS2009

WRF Simulations of the 20-22 January 2007 Snow Events over Eastern Canada: Comparison with In-Situ and Satellite Observations

One of the grand challenges of the Global Precipitation Measurement (GPM) mission is to improve cold season precipitation measurements in middle and high latitudes through the use of high-frequency passive microwave radiometry. For this, the Weather Research and Forecasting (WRF) model with the Goddard microphysics scheme is coupled with a satellite data simulation unit (WRF-SDSU) that has been developed to facilitate over-land snowfall retrieval algorithms by providing a virtual cloud library and microwave brightness temperature (Tb) measurements consistent with the GPM Microwave Imager (GMI). This study tested the Goddard cloud microphysics scheme in WRF for two snowstorm events, a lake effect and a synoptic event, that occurred between 20 and 22 January 2007 over the Canadian CloudSAT/CALIPSO Validation Project (C3VP) site in Ontario, Canada. The 24h-accumulated snowfall predicted by the WRF model with the Goddard microphysics was comparable to the observed accumulated snowfall by the ground-based radar for both events. The model correctly predicted the onset and ending of both snow events at the CARE site. WRF simulations capture the basic cloud properties as seen by the ground-based radar and satellite (i.e., CloudSAT, AMSU-B) observations as well as the observed cloud streak organization in the lake event. This latter result reveals that WRF was able to capture the cloud macro-structure reasonably well.

Shi, J. J.↗
NASA NTRS1977

Interpolation schemes for orbit determination with the global positioning system

This paper demonstrates that the Global Positioning System (GPS) and simple interpolation schemes can be utilized to satisfy typical orbit determination demands of applications satellites. The complete GPS consists of 24 satellites and permits a position fix of a user satellite at any arbitrary instant. The interpolation formulae used in this report fit generalized Keplerian orbits through a number of position fixes. For a given accuracy level the telemetry requirements as measured by the time intervals between position fixes and the computational load of an interpolation as measured by the number of points through which a generalized Keplerian orbit is fitted vary inversely. A set of possible compromises between these two factors is presented. The Phase I GPS consists of 6 satellites and permits a position fix of a user satellite just 36% of the time. It is shown that with this system more sophisticated interpolation schemes which model atmospheric drag and higher degree terms of the earth's gravity field must be employed to obtain accurate orbits.

Argentiero, P.↗
NASA NTRS2009

LDPC-PPM Coding Scheme for Optical Communication

In a proposed coding-and-modulation/demodulation-and-decoding scheme for a free-space optical communication system, an error-correcting code of the low-density parity-check (LDPC) type would be concatenated with a modulation code that consists of a mapping of bits to pulse-position-modulation (PPM) symbols. Hence, the scheme is denoted LDPC-PPM. This scheme could be considered a competitor of a related prior scheme in which an outer convolutional error-correcting code is concatenated with an interleaving operation, a bit-accumulation operation, and a PPM inner code. Both the prior and present schemes can be characterized as serially concatenated pulse-position modulation (SCPPM) coding schemes. Figure 1 represents a free-space optical communication system based on either the present LDPC-PPM scheme or the prior SCPPM scheme. At the transmitting terminal, the original data (u) are processed by an encoder into blocks of bits (a), and the encoded data are mapped to PPM of an optical signal (c). For the purpose of design and analysis, the optical channel in which the PPM signal propagates is modeled as a Poisson point process. At the receiving terminal, the arriving optical signal (y) is demodulated to obtain an estimate (a^) of the coded data, which is then processed by a decoder to obtain an estimate (u^) of the original data.

Barsoum, Maged↗
NASA NTRS1987

A multiprocessing architecture for real-time monitoring

A multiprocessing architecture for performing real time monitoring and analysis using knowledge-based problem solving techniques is discussed. To handle asynchronous inputs and perform in real time, the system consists of three or more separate processes which run concurrently on one or more processors and communicate via a message passing scheme. The Data Management Process gathers, compresses, scales and sends the incoming telemetry data to other tasks. The Inference Process consists of a proprietary high performance inference engine that runs at 1000 rules per second using telemetry data to perform real time analysis on the state and health of the Space Telescope. The multiprocessing architecture has been interfaced to a simulator and is able to process the incoming telemetry in real time.

Laffey, Thomas J.↗
DOE OSTI2023

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Direct adaptive control of a PUMA 560 industrial robot

The implementation and experimental validation of a new direct adaptive control scheme on a PUMA 560 industrial robot is described. The testbed facility consists of a Unimation PUMA 560 six-jointed robot and controller, and a DEC MicroVAX II computer which hosts the Robot Control C Library software. The control algorithm is implemented on the MicroVAX which acts as a digital controller for the PUMA robot, and the Unimation controller is effectively bypassed and used merely as an I/O device to interface the MicroVAX to the joint motors. The control algorithm for each robot joint consists of an auxiliary signal generated by a constant-gain Proportional plus Integral plus Derivative (PID) controller, and an adaptive position-velocity (PD) feedback controller with adjustable gains. The adaptive independent joint controllers compensate for the inter-joint couplings and achieve accurate trajectory tracking without the need for the complex dynamic model and parameter values of the robot. Extensive experimental results on PUMA joint control are presented to confirm the feasibility of the proposed scheme, in spite of strong interactions between joint motions. Experimental results validate the capabilities of the proposed control scheme. The control scheme is extremely simple and computationally very fast for concurrent processing with high sampling rates.

Seraji, Homayoun↗
NASA NTRS1977

Composite methods for hyperbolic equations

A composite approximation procedure combining the properties of the Lax-Wendroff and leapfrog algorithms is proposed for solving hyperbolic equations. For a one-dimensional equation, a three-step approximation consisting of a two-step Richtmeyer method followed by a leapfrog step is considered. This is a two-level scheme, so all difficulties, including storage requirements, associated with the three-level leapfrog are eliminated. For two-dimensional problems a generalization of the preceding method is used consisting of a rotated Richtmeyer method followed by a modified leapfrog step. It is found that the composite schemes are effective in reducing oscillations and nonlinear instabilities that affect the leapfrog method. The dissipation in the composite schemes is much less than in the Richtmeyer algorithm, and hence can be used for long term integrations.

Turkel, E.↗