SEARCH · Engineering Papers

Results for “partitioned methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS2013

Regression Verification Using Impact Summaries

Regression verification techniques are used to prove equivalence of syntactically similar programs. Checking equivalence of large programs, however, can be computationally expensive. Existing regression verification techniques rely on abstraction and decomposition techniques to reduce the computational effort of checking equivalence of the entire program. These techniques are sound but not complete. In this work, we propose a novel approach to improve scalability of regression verification by classifying the program behaviors generated during symbolic execution as either impacted or unimpacted. Our technique uses a combination of static analysis and symbolic execution to generate summaries of impacted program behaviors. The impact summaries are then checked for equivalence using an o-the-shelf decision procedure. We prove that our approach is both sound and complete for sequential programs, with respect to the depth bound of symbolic execution. Our evaluation on a set of sequential C artifacts shows that reducing the size of the summaries can help reduce the cost of software equivalence checking. Various reduction, abstraction, and compositional techniques have been developed to help scale software verification techniques to industrial-sized systems. Although such techniques have greatly increased the size and complexity of systems that can be checked, analysis of large software systems remains costly. Regression analysis techniques, e.g., regression testing [16], regression model checking [22], and regression verification [19], restrict the scope of the analysis by leveraging the differences between program versions. These techniques are based on the idea that if code is checked early in development, then subsequent versions can be checked against a prior (checked) version, leveraging the results of the previous analysis to reduce analysis cost of the current version. Regression verification addresses the problem of proving equivalence of closely related program versions [19]. These techniques compare two programs with a large degree of syntactic similarity to prove that portions of one program version are equivalent to the other. Regression verification can be used for guaranteeing backward compatibility, and for showing behavioral equivalence in programs with syntactic differences, e.g., when a program is refactored to improve its performance, maintainability, or readability. Existing regression verification techniques leverage similarities between program versions by using abstraction and decomposition techniques to improve scalability of the analysis [10, 12, 19]. The abstractions and decomposition in the these techniques, e.g., summaries of unchanged code [12] or semantically equivalent methods [19], compute an over-approximation of the program behaviors. The equivalence checking results of these techniques are sound but not complete-they may characterize programs as not functionally equivalent when, in fact, they are equivalent. In this work we describe a novel approach that leverages the impact of the differences between two programs for scaling regression verification. We partition program behaviors of each version into (a) behaviors impacted by the changes and (b) behaviors not impacted (unimpacted) by the changes. Only the impacted program behaviors are used during equivalence checking. We then prove that checking equivalence of the impacted program behaviors is equivalent to checking equivalence of all program behaviors for a given depth bound. In this work we use symbolic execution to generate the program behaviors and leverage control- and data-dependence information to facilitate the partitioning of program behaviors. The impacted program behaviors are termed as impact summaries. The dependence analyses that facilitate the generation of the impact summaries, we believe, could be used in conjunction with other abstraction and decomposition based approaches, [10, 12], as a complementary reduction technique. An evaluation of our regression verification technique shows that our approach is capable of leveraging similarities between program versions to reduce the size of the queries and the time required to check for logical equivalence. The main contributions of this work are: - A regression verification technique to generate impact summaries that can be checked for functional equivalence using an off-the-shelf decision procedure. - A proof that our approach is sound and complete with respect to the depth bound of symbolic execution. - An implementation of our technique using the LLVMcompiler infrastructure, the klee Symbolic Virtual Machine [4], and a variety of Satisfiability Modulo Theory (SMT) solvers, e.g., STP [7] and Z3 [6]. - An empirical evaluation on a set of C artifacts which shows that the use of impact summaries can reduce the cost of regression verification.

Backes, John↗

NASA NTRS2019

Definition of a Technology Validation Mission for P-band Reflectometry using Signals of Opportunity

Root-Zone Soil Moisture (RZSM) (moisture profile in the top meter of soil) and Snow Water Equivalent (SWE) (total snow pack water content) are identified as priority target variables in the ESAS 2017 decadal survey [1] with critical roles in hydrology and water management. RZSM estimates are vital for understanding multiple Earth system processes and forecasting (for example, droughts [2]). Simultaneous knowledge of surface and RZSM could enable a breakthrough in estimating key unobserved hydrologic fluxes and reduce uncertainty in net ecosystem exchange (NEE), carbon balance [3] discharge estimates, and crop yield forecasts [4] .With the high albedo and insulating properties of snow, monitoring, SWE accumulation would provide a key constraint on the potential runoff during spring ablation while monitoring SWE disappearance rates would provide a key constraint on SWE partition into runoff vs. infiltration/recharge. [5] demonstrated that knowledge of early-spring SWE generally contributes most to streamflow forecast skill in the Western U.S. SWE is also a source of water storage that provides the water resources during spring snowmelt. Despite such potentially transformative contributions, accurate RZSM and SWE measurements are unattainable with current technology. While active/passive L-band methods (e.g. SMAP, SMOS) can reliably retrieve surface soil moisture in the top 5 cm of soil [6], [7]. RZSM estimates are only available through model assimilation of brightness temperatures with a radiative transfer and land surface models [8]. SWE estimation uses multi-frequency passive microwave techniques (e.g. [9]-[11]), which have significant problems with deeper snow and in forested and mountainous environments [12]. Signals of opportunity (SoOp) in P-band (200-400 MHz) is a new remote sensing technique with the capability of estimating both essential hydrologic variables, RZSM and SWE, circumventing many of the aforementioned limitations under all weather conditions day and night. SoOp is the re-utilization of existing powerful satellite transmissions within bands allocated for communications or navigation. P-band SoOp sensitivity to soil moisture has been demonstrated in an airborne experiment over Oklahoma in 2016 [13]. Recent theory [14] and experiments [15] have also confirmed that the reflection coefficient phase is proportional to SWE.

Garrison, J. L.↗

NASA NTRS1990

Localized basis functions and other computational improvements in variational nonorthogonal basis function methods for quantum mechanical scattering problems involving chemical reactions

The Generalized Newton Variational Principle for 3D quantum mechanical reactive scattering is briefly reviewed. Then three techniques are described which improve the efficiency of the computations. First, the fact that the Hamiltonian is Hermitian is used to reduce the number of integrals computed, and then the properties of localized basis functions are exploited in order to eliminate redundant work in the integral evaluation. A new type of localized basis function with desirable properties is suggested. It is shown how partitioned matrices can be used with localized basis functions to reduce the amount of work required to handle the complex boundary conditions. The new techniques do not introduce any approximations into the calculations, so they may be used to obtain converged solutions of the Schroedinger equation.

Schwenke, David W.↗

NASA NTRS1977

Trace element geochemistry and the early lunar differentiation

The results are considered of experiments in which thermal neutron activation analysis was used to measure the halogens, U, Te, Ru, Os, and Hg in lunar samples. A calorimetric method was used for P2O5. Three Apollo 17 basalts, one Apollo 11 basalt, and one Apollo 15 basalt were used for the Ru-Os study. Os and Ru were measured in seven samples of breccia 14305 in a series of experiments to check an apparent anomaly in the isotopic composition of the Hg which it contains. A lunar model for Ru-Os partition is discussed. It is suggested that Ru and Os were partitioned and fractionated in the original differentiation of the lunar surface ocean. Fractionation among near-surface materials was less than a factor of 2 according to the lithic breccia data. Fractionations among basaltic sources were very large, exceeding 100, and primarily due to Ru. A clue to Ru and Os distributions at the depths from which mare basalts were derived may be provided by the relative ages of the basalts. These are interpreted as being related to partial melting in response to temperature increases at depth with time.

Jovanovic, S.↗

NASA NTRS1998

Stable Inversion for Nonlinear Nonminimum-Phase Time Varying Systems

In this paper, we extend stable inversion to nonlinear time-varying systems and study computational issues; the technique is applicable to minimum-phase as well as nonminimum-phase systems. The inversion technique is new, even in the linear time-varying case, and relies on partitioning (the dichotomic split of) the linearized system dynamics into time-varying, stable, and unstable, submanifolds. This dichotomic split is used to build time-varying filters which are, in turn, the basis of a contraction used to find a bounded inverse input-state trajectory. Finding the inverse input-state trajectory allows the development or exact-output tracking controllers. The method is local to the time-varying trajectory and requires that the internal dynamics vary slowly; however, the method represents a significant advance relative to presently available tracking controllers. Present techniques are restricted to time-invariant nonlinear systems and, in the general case, track only asymptotically.

Devasia, S.↗

NASA NTRS1989

Woven-grid sealed quasi-bipolar lead-acid battery construction and fabricating method

A quasi-bipolar lead-acid battery construction includes a plurality of bipolar cells disposed in side-by-side relation to form a stack, and a pair of monoplanar plates at opposite ends of the stack, the cell stack and monopolar plates being contained within a housing of the battery. Each bipolar cell is loaded with an electrolyte and composed of a bipolar electrode plate and a pair of separator plates disposed on opposite sides of the electrode plate and peripherally sealed thereto. Each bipolar electrode plate is composed of a partition sheet and two bipolar electrode elements folded into a hairpin configuration and applied over opposite edges of the partition sheet so as to cover the opposite surfaces of the opposite halves thereof. Each bipolar electrode element is comprised of a woven grid with a hot-melt strip applied to a central longitudinal region of the grid along which the grid is folded into the hairpin configuration, and layers of negative and positive active material pastes applied to opposite halves of the grid on opposite sides of the central hot-melt strip. The grid is made up of strands of conductive and non-conductive yarns composing the respective transverse and longitudinal weaves of the grid. The conductive yarn has a multi-stranded glass core surrounded and covered by a lead sheath, whereas the non-conductive yarn has a multi-stranded glass core surrounded and covered by a thermally activated sizing.

Rippel, Wally E.↗

NASA NTRS1997

An Adaptive Flow Solver for Air-Borne Vehicles Undergoing Time-Dependent Motions/Deformations

This report describes a concurrent Euler flow solver for flows around complex 3-D bodies. The solver is based on a cell-centered finite volume methodology on 3-D unstructured tetrahedral grids. In this algorithm, spatial discretization for the inviscid convective term is accomplished using an upwind scheme. A localized reconstruction is done for flow variables which is second order accurate. Evolution in time is accomplished using an explicit three-stage Runge-Kutta method which has second order temporal accuracy. This is adapted for concurrent execution using another proven methodology based on concurrent graph abstraction. This solver operates on heterogeneous network architectures. These architectures may include a broad variety of UNIX workstations and PCs running Windows NT, symmetric multiprocessors and distributed-memory multi-computers. The unstructured grid is generated using commercial grid generation tools. The grid is automatically partitioned using a concurrent algorithm based on heat diffusion. This results in memory requirements that are inversely proportional to the number of processors. The solver uses automatic granularity control and resource management techniques both to balance load and communication requirements, and deal with differing memory constraints. These ideas are again based on heat diffusion. Results are subsequently combined for visualization and analysis using commercial CFD tools. Flow simulation results are demonstrated for a constant section wing at subsonic, transonic, and a supersonic case. These results are compared with experimental data and numerical results of other researchers. Performance results are under way for a variety of network topologies.

Singh, Jatinder↗

NASA NTRS2002

Load Balancing Strategies for Multi-Block Overset Grid Applications

The multi-block overset grid method is a powerful technique for high-fidelity computational fluid dynamics (CFD) simulations about complex aerospace configurations. The solution process uses a grid system that discretizes the problem domain by using separately generated but overlapping structured grids that periodically update and exchange boundary information through interpolation. For efficient high performance computations of large-scale realistic applications using this methodology, the individual grids must be properly partitioned among the parallel processors. Overall performance, therefore, largely depends on the quality of load balancing. In this paper, we present three different load balancing strategies far overset grids and analyze their effects on the parallel efficiency of a Navier-Stokes CFD application running on an SGI Origin2000 machine.

Djomehri, M. Jahed↗

NASA NTRS2023

Thermomechanical Real-Time Hybrid Simulation: Conceptual Framework and Control Requirements

Real-time hybrid simulation (RTHS) is an enabling technology that has transformed engineering experimentation and helped researchers expand modeling capabilities. However, breakthroughs are necessary to expand the range of hybrid simulation methods and, thus, enable experiments with loading conditions representing multiple hazards. This paper discusses the development of a new thermomechanical RTHS framework and a systematic approach to determining RTHS control requirements. First, the framework is established using a representative finite element model of a layered structural system subjected to thermal loading. A complete two-layer system model serves as the reference system, and it is then partitioned into a numerical layer and an experimental layer that share interface conditions. Next, a thermal actuator is introduced to impose dynamic thermal loading on the experimental subsystem, serving as a transfer system. Finally, control and performance metrics are defined to evaluate the realization of interface boundary conditions and map this to the RTHS execution. Through an illustrative example considering the influence of temperature on a lunar habitat, we demonstrate how to establish controller requirements for RTHS and demonstrate that this approach can be used to conduct RTHS on structures with thermomechanical loading.

Herta Montoya↗

NASA NTRS1988

Fatigue life prediction modeling for turbine hot section materials

This paper presents a summary of the life prediction methods developed under the NASA Lewis Research Center's Hot Section Technology (HOST) program. A major objective of the fatigue and fracture efforts under the HOST program was to significantly improve the analytic life prediction tools used by the aeronautical gas turbine engine industry. This has been achieved in the areas of high-temperature thermal and mechanical fatigue of bare and coated high-temperature superalloys. Such technical improvements will eventually reduce life cycle costs. The cyclic crack initiation and propagation resistance of nominally isotropic polycrystalline alloys and highly anisotropic single crystal alloys has been addressed. A sizeable data base has been generated for three alloys [cast PWA 1455 (B-1900 + Hf), wrought Inconel 718, and cast single crystal PWA 1480] in bare and coated conditions. Two coatings systems, diffusion aluminide (PWA 273) and plasma sprayed MCrAIY overlay (PWA 286) were employed. Life prediction modeling efforts were devoted to creep-fatigue interaction, oxidation, coatings interactions, multiaxiality of stress-strain states, mean stress effects, cumulative damage, and thermomechanical fatigue. The fatigue crack initiation life models developed to date include the Cyclic Damage Accumulation (CDA) Model of Pratt & Whitney and the Total Strain Version of Strainrange Partitioning (TS-SRP) of NASA Lewis for nominally isotropic materials, and the Tensile Hysteretic Energy Model of Pratt & Whitney for anisotropic superalloys. The fatigue model being developed by the General Electric Company is based upon the concepts of Path-lndependent Integrals (PII) for describing cyclic crack growth under complex nonlinear response at the crack tip due to thermomechanical loading conditions. A micromechanistic oxidation crack extension model has been derived by researchers at Syracuse University. The models are described and discussed in the paper. Only limited verification has been achieved to-date as several of the technical programs are still in progress and the verification tasks are scheduled, quite naturally, near the conclusion of the program. To-date, efforts have concentrated on development of independent models for cyclic constitutive behavior, cyclic crack initiation, and cyclic crack propagation. The transition between crack initiation and crack propagation has not been thoroughly researched as yet, and the integration of these models into a unified life prediction method has not been addressed.

Metal fatigue↗

NASA NTRS1987

Planning and scheduling for robotic assembly

A system for reasoning about robotic assembly tasks is described. The first element of this system is a facility for itemizing the constraints which determine the admissible orderings over the activities to be sequenced. The second element is a facility which partitions the activities into independent subtasks and produces a set of admissible strategies for each. Finally, the system has facilities for constructing an admissible sequence of activities which is consistent with the given constraints. This can be done off-line, in advance of task execution, or it can be done incrementally, at execution time, according to conditions in the execution environment. The language of temporal constraints and the methods of inference presented in related papers are presented. It is shown how functional and spatial relationships between components impose temporal constraints on the order of assembly and how temporal constraints then imply admissible strategies and feasible sequences.

Fox, Barry R.↗

NASA NTRS1993

Impeller flow field characterization with a laser two-focus velocimeter

Use of Computational Fluid Dynamics (CFD) codes, prevalent in the rocket engine turbomachinery industry, necessitates data of sufficient quality and quantity to benchmark computational codes. Existing data bases for typical rocket engine configurations, in particular impellers, are limited. In addition, traditional data acquisition methods have several limitations: typically transducer uncertainties are 0.5% of transducer full scale and traditional pressure probes are unable to provide flow characteristics in the circumferential (blade-to-blade) direction. Laser velocimetry circumvents these limitations by providing +0.5% uncertainty in flow velocity and +0.5% uncertainty in flow angle. The percent of uncertainty in flow velocity is based on the measured value, not full range capability. The laser electronics multiple partitioning capability allows data acquired between blades as the impeller rotates, to be analyzed separately, thus providing blade-to-blade flow characterization. Unlike some probes, the non-intrusive measurements made with the laser velocimeter does not disturb the flow. To this end,, and under Contract (NAS8-38864) to the National Aeronautics and Space Administration (NASA) at Marshall Space Flight Center (MSFC), an extensive test program was undertaken at Rocketdyne. Impellers from two different generic rocket engine pump configurations were examined. The impellers represent different spectrums of pump design: the Space Shuttle Main Engine (SSME) high pressure fuel turbopump (HPFTP) impeller was designed in the 1 1970's the Consortium for CFD application in Propulsion Technology Pump Stage Technology Team (Pump Consortium) optimized impeller was designed with the aid of modern computing techniques. The tester configuration for each of the impellers consisted of an axial inlet, an inducer, a diffuser, and a crossover discharge. While the tested configurations were carefully chosen to be representative of generic rocket engine pumps, several features of both testers were intentionally atypical. A crossover discharge, downstream of the impeller, rather than a volute discharge was used to minimize asymmetric flow conditions that might be reflected in the impeller discharge flow data. Impeller shroud wear ring radial clearances were purposely close to minimize leakage flow, thus increasing confidence in using the inlet data as an input to CFD programs. The empirical study extensively examined the flow fields of the two impellers via performance of laser two-focus velocimeter surveys in an axial plane upstream of the impellers and in multiple radial planes downstream of the impellers. Both studies were performed at the impeller design flow coefficients. Inlet laser surveys that provide CFD code inlet boundary conditions were performed in one axial plane, with ten radial locations surveyed. Three wall static pressures, positioned circumferentially around the impeller inlet, were used to identify asymmetrical pressure distributions in the inlet survey plane. The impeller discharge flow characterization consisted of three radial planes for the SSME HPFTP impeller and two radial planes for the Pump Consortium optimized impeller. Housing wall static pressures were placed to correspond to the radial locations surveyed with the laser velocimeter. Between five and thirteen axial stations across the discharge channel width were examined in each radial plane during the extensive flow mapping. The largely successful empirical flow characterization of two different impellers resulted in a substantial contribution to the limited existing data base, and yielded accurate data for CFD code benchmarking.

Brozowski, L. A.↗

NASA NTRS1999

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗

NASA NTRS1977

Skylab water balance analysis

The water balance of the Skylab crew was analyzed. Evaporative water loss using a whole body input/output balance equation, water, body tissue, and energy balance was analyzed. The approach utilizes the results of several major Skylab medical experiments. Subsystems were designed for the use of the software necessary for the analysis. A partitional water balance that graphically depicts the changes due to water intake is presented. The energy balance analysis determines the net available energy to the individual crewman during any period. The balances produce a visual description of the total change of a particular body component during the course of the mission. The information is salvaged from metabolic balance data if certain techniques are used to reduce errors inherent in the balance method.

Leonard, J. I.↗

NASA NTRS1987

Parallel processing for digital picture comparison

In picture processing an important problem is to identify two digital pictures of the same scene taken under different lighting conditions. This kind of problem can be found in remote sensing, satellite signal processing and the related areas. The identification can be done by transforming the gray levels so that the gray level histograms of the two pictures are closely matched. The transformation problem can be solved by using the packing method. Researchers propose a VLSI architecture consisting of m x n processing elements with extensive parallel and pipelining computation capabilities to speed up the transformation with the time complexity 0(max(m,n)), where m and n are the numbers of the gray levels of the input picture and the reference picture respectively. If using uniprocessor and a dynamic programming algorithm, the time complexity will be 0(m(3)xn). The algorithm partition problem, as an important issue in VLSI design, is discussed. Verification of the proposed architecture is also given.

Cheng, H. D.↗

NASA NTRS1991

Distributed memory compiler design for sparse problems

A compiler and runtime support mechanism is described and demonstrated. The methods presented are capable of solving a wide range of sparse and unstructured problems in scientific computing. The compiler takes as input a FORTRAN 77 program enhanced with specifications for distributing data, and the compiler outputs a message passing program that runs on a distributed memory computer. The runtime support for this compiler is a library of primitives designed to efficiently support irregular patterns of distributed array accesses and irregular distributed array partitions. A variety of Intel iPSC/860 performance results obtained through the use of this compiler are presented.

Wu, Janet↗

NASA NTRS1989

Efficient detection of a CW signal with a linear frequency drift

An efficient method is presented for the detection of a continuous wave (CW) signal with a frequency drift that is linear in time. Signals of this type occur in transmissions between any two locations that are accelerating relative to one another, e.g., transmissions from the Voyager spacecraft. We assume that both the frequency and the drift are unknown. We also assume that the signal is weak compared to the Gaussian noise. The signal is partitioned into subsequences whose discrete Fourier transforms provide a sequence of instantaneous spectra at equal time intervals. These spectra are then accumulated with a shift that is proportional to time. When the shift is equal to the frequency drift, the signal to noise ratio increases and detection occurs. Here, we show how to compute these accumulations for many shifts in an efficient manner using a variety of Fast Fourier Transformations (FFT). Computing time is proportional to L log L where L is the length of the time series.

Swarztrauber, Paul N.↗

NASA NTRS1987

Mapping a battlefield simulation onto message-passing parallel architectures

Perhaps the most critical problem in distributed simulation is that of mapping: without an effective mapping of workload to processors the speedup potential of parallel processing cannot be realized. Mapping a simulation onto a message-passing architecture is especially difficult when the computational workload dynamically changes as a function of time and space; this is exactly the situation faced by battlefield simulations. This paper studies an approach where the simulated battlefield domain is first partitioned into many regions of equal size; typically there are more regions than processors. The regions are then assigned to processors; a processor is responsible for performing all simulation activity associated with the regions. The assignment algorithm is quite simple and attempts to balance load by exploiting locality of workload intensity. The performance of this technique is studied on a simple battlefield simulation implemented on the Flex/32 multiprocessor. Measurements show that the proposed method achieves reasonable processor efficiencies. Furthermore, the method shows promise for use in dynamic remapping of the simulation.

Nicol, David M.↗