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

Analysis and Improvement of Upwind and Centered Schemes on Quadrilateral and Triangular Meshes

Second-order accurate upwind and centered schemes are presented in a framework that facilitates their analysis and comparison. The upwind scheme employed consists of a reconstruction step (MUSCL approach) followed by an upwind step (Roe's flux-difference splitting). The two centered schemes are of Lax-Friedrichs (L-F) type. They are the nonstaggered versions of the Nessyahu-Tadmor (N-T) scheme and the CE/SE method with epilson = 1/2. The upwind scheme is extended to the case of two spatial dimensions (2D) in a straightforward manner. The N-T and CE/SE schemes are extended in a manner similar to the 2D extensions of the CE/SE scheme by Wang and Chang for a triangular mesh and by Zhang, Yu, and Chang for a quadrilateral mesh. The slope estimates, however, are simplified. Fourier stability and accuracy analyses are carried out for these schemes for the standard 1D and the 2D quadrilateral mesh cases. In the nonstandard case of a triangular mesh, the triangles must be paired up when analyzing the upwind and N-T schemes. An observation resulting in an extended N-T scheme which is faster and uses only one-third of the storage for flow data compared with the CE/SE method is presented. Numerical results are shown. Other improvements to the schemes are discussed.

Huynh, H. T.↗

Investigation of Advanced Counterrotation Blade Configuration Concepts for High Speed Turboprop Systems. Task 3: Advanced Fan Section Grid Generator Final Report and Computer Program User's Manual

A procedure is studied for generating three-dimensional grids for advanced turbofan engine fan section geometries. The procedure constructs a discrete mesh about engine sections containing the fan stage, an arbitrary number of axisymmetric radial flow splitters, a booster stage, and a bifurcated core/bypass flow duct with guide vanes. The mesh is an h-type grid system, the points being distributed with a transfinite interpolation scheme with axial and radial spacing being user specified. Elliptic smoothing of the grid in the meridional plane is a post-process option. The grid generation scheme is consistent with aerodynamic analyses utilizing the average-passage equation system developed by Dr. John Adamczyk of NASA Lewis. This flow solution scheme requires a series of blade specific grids each having a common axisymmetric mesh, but varying in the circumferential direction according to the geometry of the specific blade row.

Crook, Andrew J.↗

Compiler-assisted multiple instruction rollback recovery using a read buffer

Multiple instruction rollback (MIR) is a technique that has been implemented in mainframe computers to provide rapid recovery from transient processor failures. Hardware-based MIR designs eliminate rollback data hazards by providing data redundancy implemented in hardware. Compiler-based MIR designs have also been developed which remove rollback data hazards directly with data-flow transformations. This paper describes compiler-assisted techniques to achieve multiple instruction rollback recovery. We observe that some data hazards resulting from instruction rollback can be resolved efficiently by providing an operand read buffer while others are resolved more efficiently with compiler transformations. The compiler-assisted scheme presented consists of hardware that is less complex than shadow files, history files, history buffers, or delayed write buffers, while experimental evaluation indicates performance improvement over compiler-based schemes.

Alewine, Neal J.↗

Accurate upwind methods for the Euler equations

A new class of piecewise linear methods for the numerical solution of the one-dimensional Euler equations of gas dynamics is presented. These methods are uniformly second-order accurate, and can be considered as extensions of Godunov's scheme. With an appropriate definition of monotonicity preservation for the case of linear convection, it can be shown that they preserve monotonicity. Similar to Van Leer's MUSCL scheme, they consist of two key steps: a reconstruction step followed by an upwind step. For the reconstruction step, a monotonicity constraint that preserves uniform second-order accuracy is introduced. Computational efficiency is enhanced by devising a criterion that detects the 'smooth' part of the data where the constraint is redundant. The concept and coding of the constraint are simplified by the use of the median function. A slope steepening technique, which has no effect at smooth regions and can resolve a contact discontinuity in four cells, is described. As for the upwind step, existing and new methods are applied in a manner slightly different from those in the literature. These methods are derived by approximating the Euler equations via linearization and diagonalization. At a 'smooth' interface, Harten, Lax, and Van Leer's one intermediate state model is employed. A modification for this model that can resolve contact discontinuities is presented. Near a discontinuity, either this modified model or a more accurate one, namely, Roe's flux-difference splitting. is used. The current presentation of Roe's method, via the conceptually simple flux-vector splitting, not only establishes a connection between the two splittings, but also leads to an admissibility correction with no conditional statement, and an efficient approximation to Osher's approximate Riemann solver. These reconstruction and upwind steps result in schemes that are uniformly second-order accurate and economical at smooth regions, and yield high resolution at discontinuities.

Huynh, Hung T.↗

Velocity autocorrelation functions in model liquid metals

Starting from interatomic potentials and static radial distribution functions, a self-consistent iteration scheme has been used to calculate velocity autocorrelation functions in liquid metals. The interatomic forces are treated directly. The calculation bypasses the details of the many-body dynamics and it is not necessary to introduce any additional parameters. Several simplifications may be used without introducing appreciable deviations. The results are in good agreement with computer experiments on liquid sodium at 383 K, suggesting that the velocity autocorrelation function may be a simpler quantity than previously supposed.

Tsang, T.↗

TAP 2: A finite element program for thermal analysis of convectively cooled structures

A finite element computer program (TAP 2) for steady-state and transient thermal analyses of convectively cooled structures is presented. The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature-dependent thermal parameters is performed using the Newton-Raphson iteration method. Transient analyses are performed using an implicit Crank-Nicolson time integration scheme with consistent or lumped capacitance matrices as an option. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

Detection of interstellar CCD

The first detection of interstellar CCD has been obtained in the N = 3-2, J = 7/2-5/2 line at 216.3732 GHz towards the Kleinman-Low nebula. An abundance ratio CCD/CCH of 0.045 is found, indicating in this molecule a deuterium enhancement similar to that found for HNC, but an order of magnitude higher than for HCN. Negative results towards DR 21(OH), NGC 2264, L 134N, TM Cl, Rho OphB2, and IRC 10216 rule out a stronger D enhancement for CCD than for other deuterated molecules in these sources. Predictions of currently developed ion-molecule reaction schemes are consistent with the CCD line detected and the negative results.

Combes, F.↗

Building a generalized distributed system model

The key elements in the 1992-93 period of the project are the following: (1) extensive use of the simulator to implement and test - concurrency control algorithms, interactive user interface, and replica control algorithms; and (2) investigations into the applicability of data and process replication in real-time systems. In the 1993-94 period of the project, we intend to accomplish the following: (1) concentrate on efforts to investigate the effects of data and process replication on hard and soft real-time systems - especially we will concentrate on the impact of semantic-based consistency control schemes on a distributed real-time system in terms of improved reliability, improved availability, better resource utilization, and reduced missed task deadlines; and (2) use the prototype to verify the theoretically predicted performance of locking protocols, etc.

Mukkamala, R.↗

ERS-1 scatterometer measurements over the Southern Ocean

Backscatter cross section measurements from the ERS-1 Active Microwave Instrument (AMI) scatterometer were reprocessed to vector winds using the Freilich-Dunbar model function and a meteorologically aided ambiguity removal scheme. This consistent data set was used to examine the wind field over the Southern Ocean from 20 to 60 deg South. The large number of ERS-1 measurements allows relatively accurate calculation of annual mean wind, stress, and curl fields as well as overall statistics of the winds at mid to high southern latitudes. The long duration of the data time series allows preliminary examination of low frequency (semi annual) wind variability.

Freilich, M. H.↗

Impact of the 1997-1998 El Nino on Atmospheric Circulation

The 1997-1998 E1 Nino created significant anomalies in global circulation patterns. We monitored the E1 Nino event, using TOVS (TIROS Operational Vertical Sounder) data, as part of the Laboratory for Atmospheres TOVS Pathfinder Data Set. TOVS has flown on the NOAA operational polar orbiting satellites from November 1978 to present. We have analyzed data from NOAA 9, 10, 11, 12, and 14, covering January 1985 to the present, using a consistent processing scheme. The data set contains 2-4 times daily global fields of land/ocean surface skin temperature, atmospheric temperature-moisture profiles, cloud top pressure and fractional cloud cover, OLR and clear sky Outgoing Longwave Radiation (OLR), and precipitation estimate. We have generated anomalies of all fields compared to a 12 year climatology of our data set, covering the period 1985-1996. Anomalies during the current E1 Nino were compared to those of other E1 Nino-La Nina episodes since 1985, with particular attention paid to tropical anomalies of surface skin and air temperature, precipitation, upper tropospheric water vapor, OLR, and clear sky OLR. All the El Nino-La Nina events produce similar patterns, but the magnitude of the current anomalies in all fields is considerably greater than those during the previous weaker surface skin temperature anomaly episodes. Significant extra-tropical anomalies were found during the current E1 Nino as well.

Susskind, Joel↗

Ultraspectral Sounding Retrieval Error Budget and Estimation

The ultraspectral infrared radiances obtained from satellite observations provide atmospheric, surface, and/or cloud information. The intent of the measurement of the thermodynamic state is the initialization of weather and climate models. Great effort has been given to retrieving and validating these atmospheric, surface, and/or cloud properties. Error Consistency Analysis Scheme (ECAS), through fast radiative transfer model (RTM) forward and inverse calculations, has been developed to estimate the error budget in terms of absolute and standard deviation of differences in both spectral radiance and retrieved geophysical parameter domains. The retrieval error is assessed through ECAS without assistance of other independent measurements such as radiosonde data. ECAS re-evaluates instrument random noise, and establishes the link between radiometric accuracy and retrieved geophysical parameter accuracy. ECAS can be applied to measurements of any ultraspectral instrument and any retrieval scheme with associated RTM. In this paper, ECAS is described and demonstration is made with the measurements of the METOP-A satellite Infrared Atmospheric Sounding Interferometer (IASI)..

Zhou, Daniel K.↗

Bit-Serial Adder Based on Quantum Dots

A proposed integrated circuit based on quantum-dot cellular automata (QCA) would function as a bit-serial adder. This circuit would serve as a prototype building block for demonstrating the feasibility of quantum-dots computing and for the further development of increasingly complex and increasingly capable quantum-dots computing circuits. QCA-based bit-serial adders would be especially useful in that they would enable the development of highly parallel and systolic processors for implementing fast Fourier, cosine, Hartley, and wavelet transforms. The proposed circuit would complement the QCA-based circuits described in "Implementing Permutation Matrices by Use of Quantum Dots" (NPO-20801), NASA Tech Briefs, Vol. 25, No. 10 (October 2001), page 42 and "Compact Interconnection Networks Based on Quantum Dots" (NPO-20855), which appears elsewhere in this issue. Those articles described the limitations of very-large-scale-integrated (VLSI) circuitry and the major potential advantage afforded by QCA. To recapitulate: In a VLSI circuit, signal paths that are required not to interact with each other must not cross in the same plane. In contrast, for reasons too complex to describe in the limited space available for this article, suitably designed and operated QCA-based signal paths that are required not to interact with each other can nevertheless be allowed to cross each other in the same plane without adverse effect. In principle, this characteristic could be exploited to design compact, coplanar, simple (relative to VLSI) QCA-based networks to implement complex, advanced interconnection schemes. To enable a meaningful description of the proposed bit-serial adder, it is necessary to further recapitulate the description of a quantum-dot cellular automation from the first-mentioned prior article: A quantum-dot cellular automaton contains four quantum dots positioned at the corners of a square cell. The cell contains two extra mobile electrons that can tunnel (in the quantum-mechanical sense) between neighboring dots within the cell. The Coulomb repulsion between the two electrons tends to make them occupy antipodal dots in the cell. For an isolated cell, there are two energetically equivalent arrangements (denoted polarization states) of the extra electrons. The cell polarization is used to encode binary information. Because the polarization of a nonisolated cell depends on Coulomb-repulsion interactions with neighboring cells, universal logic gates and binary wires could be constructed, in principle, by arraying QCA of suitable design in suitable patterns. Again, for reasons too complex to describe here, in order to ensure accuracy and timeliness of the output of a QCA array, it is necessary to resort to an adiabatic switching scheme in which the QCA array is divided into subarrays, each controlled by a different phase of a multiphase clock signal. In this scheme, each subarray is given time to perform its computation, then its state is frozen by raising its inter-dot potential barriers and its output is fed as the input to the successor subarray. The successor subarray is kept in an unpolarized state so it does not influence the calculation of preceding subarray. Such a clocking scheme is consistent with pipeline computation in the sense that each different subarray can perform a different part of an overall computation. In other words, QCA arrays are inherently suitable for pipeline and, moreover, systolic computations. This sequential or pipeline aspect of QCA would be utilized in the proposed bit-serial adders.

Fijany, Amir↗

Adaptive Aft Signature Shaping of a Low-Boom Supersonic Aircraft Using Off-Body Pressures

The design and optimization of a low-boom supersonic aircraft using the state-of-the- art o -body aerodynamics and sonic boom analysis has long been a challenging problem. The focus of this paper is to demonstrate an e ective geometry parameterization scheme and a numerical optimization approach for the aft shaping of a low-boom supersonic aircraft using o -body pressure calculations. A gradient-based numerical optimization algorithm that models the objective and constraints as response surface equations is used to drive the aft ground signature toward a ramp shape. The design objective is the minimization of the variation between the ground signature and the target signature subject to several geometric and signature constraints. The target signature is computed by using a least-squares regression of the aft portion of the ground signature. The parameterization and the deformation of the geometry is performed with a NASA in- house shaping tool. The optimization algorithm uses the shaping tool to drive the geometric deformation of a horizontal tail with a parameterization scheme that consists of seven camber design variables and an additional design variable that describes the spanwise location of the midspan section. The demonstration cases show that numerical optimization using the state-of-the-art o -body aerodynamic calculations is not only feasible and repeatable but also allows the exploration of complex design spaces for which a knowledge-based design method becomes less effective.

Ordaz, Irian↗

Real-Time Distributed Embedded Oscillator Operating Frequency Monitoring

A document discusses the utilization of embedded clocks inside of operating network data links as an auxiliary clock source to satisfy local oscillator monitoring requirements. Modem network interfaces, typically serial network links, often contain embedded clocking information of very tight precision to recover data from the link. This embedded clocking data can be utilized by the receiving device to monitor the local oscillator for tolerance to required specifications, often important in high-integrity fault-tolerant applications. A device can utilize a received embedded clock to determine if the local or the remote device is out of tolerance by using a single link. The local device can determine if it is failing, assuming a single fault model, with two or more active links. Network fabric components, containing many operational links, can potentially determine faulty remote or local devices in the presence of multiple faults. Two methods of implementation are described. In one method, a recovered clock can be directly used to monitor the local clock as a direct replacement of an external local oscillator. This scheme is consistent with a general clock monitoring function whereby clock sources are clocking two counters and compared over a fixed interval of time. In another method, overflow/underflow conditions can be used to detect clock relationships for monitoring. These network interfaces often provide clock compensation circuitry to allow data to be transferred from the received (network) clock domain to the internal clock domain. This circuit could be modified to detect overflow/underflow conditions of the buffering required and report a fast or slow receive clock, respectively.

Pollock, Julie↗

Enhancing GPM Passive and Combined Microwave Algorithms with Dynamic Surface Information for Drizzle Retrieval and Improved Precipitation Detection Over Land

Following the 2014 launch of the Global Precipitation Measurement Mission (GPM), an unprecedented combination of coincident active and passive microwave observations are available for state of the art precipitation retrieval. The GPM Combined Algorithm forms the backbone of this effort, optimizing geophysical variables for agreement with the full suite of multi-spectral information content. These combined retrievals are then utilized, along with a radiative transfer model, as a database applied for retrievals across a constellation of passive microwave radiometers of varying frequencies. By keeping such retrievals related through the transfer standard of the combined algorithm, level 3 products such as the Integrated Multi-satellitE Retrievals for GPM (IMERG) are able to provide consistent global products for users at the higher temporal resolution required for hydrological applications. In initial versions of the combined product, precipitation retrievals are carried out only in the presence of a signal from the active radar. As a result, light precipitation and drizzle below the threshold of DPR sensitivity are not included in any of the products down the chain from the constellation to IMERG. In this work, the effects of enhancing the retrievals with a surface emissivity and non-raining water vapor retrieval using the passive observations are explored. Over both ocean and land, the surface retrieval is used to identify areas with high probability of light precipitation and drizzle which is then quantified using techniques derived from the higher sensitivity CloudSat mission. Results indicate successful inclusion of drizzle in the retrievals that can then be included in the constellation databases, as well as improvement in passive microwave false positive precipitation signals over land in cases where surface scattering was misinterpreted as precipitation signal. The inclusion of the dynamic surface information also creates a more robust, radiometrically consistent retrieval scheme for process studies and hydrologic applications.

Ringerud, Sarah↗

Deriving Climate Change Signal from Hyperspectral Sounders Using Spectral Fingerprinting Method

Hyperspectral observations from satellite-based sensors provide high information content for the Earth’s atmospheric temperature, water vapor and trace gas vertical profiles. We have developed a radiometrically consistent spectral fingerprinting method to derive climate change signals from Aqua AIRS/AMSU and S-NPP CrIS/ATMS data. The climate variables include temperature and water vapor profiles, cloud, trace gases, and surface skin temperature. The radiative kernels obtained via a single field of view physical retrieval algorithm under all-sky conditions. A key component to this work is a Principal Component-based Radiative Transfer Model (PCRTM). It is 4 orders of magnitude faster than a line-by-line radiative transfer model while keeping a similar accuracy (0.03 K RMS errors with close to zero bias). The PCRTM includes multiple scattering of clouds and non-thermodynamics equilibrium of CO2 in the RT calculations. Instead of quantifying the radiometric differences between AIRS/AMSU and CrIS/ATMS measurements directly using Simultaneous Nadir Overpass (SNO) or Double Difference Technique (DDT), we use the radiometric consistent fingerprinting scheme to derive two sets of space-time averaged anomalies from the Level 1 data of AIRS/AMSU and CrIS/ATMS. The derived anomalies in geophysical space will form a long-term, stable, and continuous climate data record. We can further infer the causes of any offset or drift by studying the differences between two overlapping data sets. For example, the offset in surface skin temperature anomaly time series will most likely caused by the Blackbody temperature calibration errors of the sounder instruments.

climate↗

An energy and potential enstrophy conserving scheme for the shallow water equations

A three-dimensional finite difference scheme for the solution of the shallow water momentum equations which accounts for the conservation of potential enstrophy in the flow of a homogeneous incompressible shallow atmosphere over steep topography as well as for total energy conservation is presented. The scheme is derived to be consistent with a reasonable scheme for potential vorticity advection in a long-term integration for a general flow with divergent mass flux. Numerical comparisons of the characteristics of the present potential enstrophy-conserving scheme with those of a scheme that conserves potential enstrophy only for purely horizontal nondivergent flow are presented which demonstrate the reduction of computational noise in the wind field with the enstrophy-conserving scheme and its convergence even in relatively coarse grids.

Arakawa, A.↗

A PDE Sensitivity Equation Method for Optimal Aerodynamic Design

The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.

Borggaard, Jeff↗