Engineering Papers⌕ Search

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.

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

At least 721 records · Page 40

Systems and methods for modeling water quality

A system, method, device and computer-readable medium for creating an ensemble model of water quality. The ensemble model is generated by determining a set of optimal component models for spectral regions of a body of water, and combining the optimal models. The optimal models can be based on remote sensing data, including satellite imagery. A K-fold partition approach or a global approach can be used to determine the optimal component models, and the optimal component models can be combined through spectral space partition rules to generate an ensemble model of water quality. The ensemble model not only has improved water quality prediction ability, but also has strong spatial and temporal extensibility. The spatial and temporal extensibility of the ensemble model is fundamentally important and desirable for long-term and large-scale remote sensing monitoring and assessment of water quality.

Liu, Hongxing↗

Estimation of Supersonic Stage Separation Aerodynamics of Winged-Body Launch Vehicles Using Response Surface Methods

A collection of statistical and mathematical techniques referred to as response surface methodology was used to estimate the longitudinal stage separation aerodynamic characteristics of a generic, bimese, winged multi-stage launch vehicle configuration using data obtained on small-scale models at supersonic speeds in the NASA Langley Research Center Unitary Plan Wind Tunnel. The simulated Mach 3 staging was dominated by multiple shock wave interactions between the orbiter and booster vehicles throughout the relative spatial locations of interest. This motivated a partitioning of the overall inference space into several contiguous regions within which the separation aerodynamics were presumed to be well-behaved and estimable using cuboidal and spherical central composite designs capable of fitting full second-order response functions. The primary goal was to approximate the underlying overall aerodynamic response surfaces of the booster vehicle in belly-to-belly proximity to the orbiter vehicle using relatively simple, lower-order polynomial functions that were piecewise-continuous across the full independent variable ranges of interest. The quality of fit and prediction capabilities of the empirical models were assessed in detail, and the issue of subspace boundary discontinuities was addressed. The potential benefits of augmenting the central composite designs to full third order using computer-generated D-optimality criteria were also evaluated. The usefulness of central composite designs, the subspace sizing, and the practicality of fitting low-order response functions over a partitioned inference space dominated by highly nonlinear and possibly discontinuous shock-induced aerodynamics are discussed.

Erickson, Gary E.↗

Quantum Simulation of the First-Quantized Pauli-Fierz Hamiltonian

We provide an explicit recursive divide-and-conquer approach for simulating quantum dynamics and derive a discrete first-quantized nonrelativistic QED Hamiltonian based on the many-particle Pauli-Fierz Hamiltonian. We apply this recursive divide-and-conquer algorithm to this Hamiltonian and compare it to a concrete simulation algorithm that uses qubitization. Our divide-and-conquer algorithm, using lowest-order Trotterization, scales for fixed grid spacing as O ~ ( Λ N 2 η 2 t 2 / ϵ ) for grid size N , η particles, simulation time t , field cutoff Λ , and error ϵ . Our qubitization algorithm scales as O ~ ( N ( η + N ) ( η + Λ 2 ) t log ( 1 / ϵ ) ) . This shows that even a naive partitioning and low-order splitting formula can yield, through our divide-and-conquer formalism, superior scaling to qubitization for large Λ . We compare the relative costs of these two algorithms on systems that are relevant for applications such as the spontaneous emission of photons and the photoionization of electrons. We observe that for different parameter regimes, one method can be favored over the other. Finally, we give new algorithmic and circuit-level techniques for gate optimization, including a new way of implementing a group of multicontrolled- X gates that can be used for better analysis of circuit cost. Published by the American Physical Society 2024

Mukhopadhyay, Priyanka (ORCID:0000000164639100)↗

Exploratory Thermal-mechanical Fatigue Results for Rene' 80 in Ultrahigh Vacuum

A limited study was conducted of the use of strainage partitioning for predicting the thermalmechanical fatigue life of cast nickel-base superalloy Rene' 80. The fatigue lives obtained by combined inphase thermal and mechanical strain cycling between 400 C (752 F) and 1000 C (1802 F) in an ultrahigh vacuum were considerably shorter than those represented by the four basic partitioned inelastic strainrange fatigue life relationships established previously for this alloy at 871 C (1600 F) and 1000 C (1832 F) in an ultrahigh vacuum. This behavior was attributed to the drastic decrease in ductility with decreasing temperature for this alloy. These results indicated that the prediction of the thermal-mechanical fatigue life of Rene' 80 by the method of strainrange partioning may be improved if based on the four basic fatigue life relationships determined at a lower temperature in the thermal-mechanical strain cycle.

Sheinker, A. A.↗

Problem size, parallel architecture and optimal speedup

The communication and synchronization overhead inherent in parallel processing can lead to situations where adding processors to the solution method actually increases execution time. Problem type, problem size, and architecture type all affect the optimal number of processors to employ. The numerical solution of an elliptic partial differential equation is examined in order to study the relationship between problem size and architecture. The equation's domain is discretized into n sup 2 grid points which are divided into partitions and mapped onto the individual processor memories. The relationships between grid size, stencil type, partitioning strategy, processor execution time, and communication network type are analytically quantified. In so doing, the optimal number of processors was determined to assign to the solution, and identified (1) the smallest grid size which fully benefits from using all available processors, (2) the leverage on performance given by increasing processor speed or communication network speed, and (3) the suitability of various architectures for large numerical problems.

Nicol, David M.↗

Relationship between fatigue life in the creep-fatigue region and stress-strain response

On the basis of mechanical tests and metallographic studies, strainrange partitioned lives were predicted by introducing stress-strain materials parameters into the Universal Slopes Equation. This was the result of correlating fatigue damage mechanisms and deformation mechanisms operating at elevated temperatures on the basis of observed mechanical and microstructural behavior. Correlation between high temperature fatigue and stress strain properties for nickel base superalloys and stainless steel substantiated the method. Parameters which must be evaluated for PP- and CC- life are the maximum stress achievable under entirely plastic and creep conditions respectively and corresponding inelastic strains, and the elastic modulus. For plasticity/creep interaction conditions (PC and CP) two more pairs of stress strain parameters must be ascertained.

Berkovits, A.↗

Problem size, parallel architecture, and optimal speedup

The communication and synchronization overhead inherent in parallel processing can lead to situations where adding processors to the solution method actually increases execution time. Problem type, problem size, and architecture type all affect the optimal number of processors to employ. The numerical solution of an elliptic partial differential equation is examined in order to study the relationship between problem size and architecture. The equation's domain is discretized into n sup 2 grid points which are divided into partitions and mapped onto the individual processor memories. The relationships between grid size, stencil type, partitioning strategy, processor execution time, and communication network type are analytically quantified. In so doing, the optimal number of processors was determined to assign to the solution, and identified (1) the smallest grid size which fully benefits from using all available processors, (2) the leverage on performance given by increasing processor speed or communication network speed, and (3) the suitability of various architectures for large numerical problems.

Nicol, David M.↗

Relationship between fatigue life in the creep-fatigue region and stress-strain response

On the basis of mechanical tests and metallographic studies, strainrange partitioned lives were predicted by introducing stress-strain materials parameters into the Universal Slopes Equation. This was the result of correlating fatigue damage mechanisms and deformation mechanisms operating at elevated temperatures on the basis of observed mechanical and microstructural behavior. Correlation between high temperature fatigue and stress strain properties for nickel base superalloys and stainless steel substantiated the method. Parameters which must be evaluated for PP- and CC- life are the maximum stress achievable under entirely plastic and creep conditions respectively and corresponding inelastic strains, and the two more pairs of stress strain parameters must be ascertained.

Berkovits, A.↗

Preconditioned implicit solvers for the Navier-Stokes equations on distributed-memory machines

The GMRES method is parallelized, and combined with local preconditioning to construct an implicit parallel solver to obtain steady-state solutions for the Navier-Stokes equations of fluid flow on distributed-memory machines. The new implicit parallel solver is designed to preserve the convergence rate of the equivalent 'serial' solver. A static domain-decomposition is used to partition the computational domain amongst the available processing nodes of the parallel machine. The SPMD (Single-Program Multiple-Data) programming model is combined with message-passing tools to develop the parallel code on a 32-node Intel Hypercube and a 512-node Intel Delta machine. The implicit parallel solver is validated for internal and external flow problems, and is found to compare identically with flow solutions obtained on a Cray Y-MP/8. A peak computational speed of 2300 MFlops/sec has been achieved on 512 nodes of the Intel Delta machine,k for a problem size of 1024 K equations (256 K grid points).

Ajmani, Kumud↗

Compositional and Microstructural Evolution of Olivine During Pulsed Laser Irradiation: Insights Based on a FIB/Field-Emission TEM Study

Introduction: The use of pulsed laser irradiation to simulate the short duration, high-energy conditions characteristic of micrometeorite impacts is now an established approach in experimental space weathering studies. The laser generates both melt and vapor deposits that contain nanophase metallic Fe (npFe(sup 0)) grains with size distributions and optical properties similar to those in natural impact-generated melt and vapor deposits. There remains uncertainty, however, about how well lasers simulate the mechanical work and internal (thermal) energy partitioning that occurs in actual impacts. We are currently engaged in making a direct comparison between the products of laser irradiation and experimental/natural hypervelocity impacts. An initial step reported here is to use analytical TEM is to attain a better understanding of how the microstructure and composition of laser deposits evolve over multiple cycles of pulsed laser irradiation. Experimental Methods: We irradiated pressed-powder pellets of San Carlos olivine (Fo(sub 90)) with up to 99 rastered pulses of a GAM ArF excimer laser. The irradiated surface of the sample were characterized by SEM imaging and areas were selected for FIB cross sectioning for TEM study using an FEI Quanta dual-beam electron/focused ion beam instrument. FIB sections were characterized using a JEOL2500SE analytical field-emission scanning transmission electron microscope (FE-STEM) optimized for quantitative element mapping at less than 10 nm spatial resolutions. Results: In the SEM the 99 pulse pressed pellet sample shows a complex, inhomogeneous, distribution of laser-generated material, largely concentrated in narrow gaps and larger depressions between grains. Local concentrations of npFe0 spherules 0.1 to 1 micrometers in size are visible within these deposits in SEM back-scatter images. Fig. 1 shows bright-field STEM images of a FIB cross-section of a one of these deposits that continuously covers the top and sloping side of an olivine grain. The deposit has 3 microstructurally distinct sub-layers composed of silicate glass with varying modal fractions and size distributions of npFe( sup 0) spherules, along with nanocrystalline silicate material. A relatively thin (50-300 nm) topmost surface layer has a high-concentration of npFe0 spherules 5-20 nm in size. Element mapping shows the layer to be enriched in Fe by a factor of 2.5 relative to the olivine substrate, with Mg and Si depleted by 20% and 10% respectively. This is compositionally complementary to the underlying, middle layer of the deposit that is depleted in Fe, enriched in Mg and has a much lower npFe0 concentration. A third layer of nanocrystalline olivine occurs at the substrate interface. Discussion: The FE-STEM results suggest the topmost layer is a vapor deposit, underlain by a thicker microstructurally complex melt-generated layer. The compositional relations suggest the melt layer was partially vaporized, preferentially losing more volatile elements (e.g., Fe). The vaporized material re-condensed to form the thin, npFe(sup 0)-rich surface deposit during or immediately after the scan cycle. Nanocrystalline olivine that grew within the melt layer as it formed and cooled is similar in volume and microstructure to what we have observed in the impact melt lining of a micrometeorite impact crater in olivine. This suggest the time-temperature relations attained in the laser sample may not be too different from a micrometeorite impact. Our TEM observations, however, do not show evidence for the same level of mechanical dam-age (e.g., fracturing) seen around the natural micrometeorite crater.

Christoffersen, R.↗

Continuous Counter‐Current Microfluidic Liquid–Liquid Extraction Achieved Using a Pair of Wettable Screen Meshes

Continuous counter‐current microfluidic liquid–liquid extraction performs separations by flowing immiscible liquids in opposing directions within a single flow channel. In principle, this flow arrangement enables a large number of theoretical separation units in a small footprint, without using interstage valving, pumping, and phase separation. Despite its potential for excellent separation performance, this microfluidic scheme rarely appears in literature due to the requirement for capillary forces to be greater than hydrodynamic forces for stable flow. We present a novel microfluidic device and flow approaches that overcome this force‐balance challenge, enabling stable, long‐duration continuous counter‐current flow. Additionally, we cover a suite of methodologies for quantifying the performance of the microfluidic device, revealing the number of theoretical equilibrium stages achieved. The enabling technologies include a woven mesh screen‐based microfluidic device architecture that is easily fabricated outside of a clean room, surface functionalization strategies to promote conjugate (organic/aqueous) wettability, flow approaches to eliminate bubbles and carryover, and computer‐aided flow automation with optical measurement of extraction performance. The reported experiments lasted for over 36 h, terminated only at experiment conclusion, where the device still exhibited good performance. Automated Raman spectroscopy was used for solute quantitation of the ternary system tert‐butanol in a toluene/water matrix, a ternary system that was specifically chosen to analyze the device's performance with a small solute partition ratio and to enable in‐line Raman measurements of solute concentrations in both phases. The microfluidic device possessed a 55 mm contact length and a 38.5 µL internal volume. During counter‐current flow, we observed approximately 37 equilibrium stages (37 ± 13) based on a best‐fit of the solute fraction remaining in the aqueous phase using a Kremser Group Method analysis.

36 MATERIALS SCIENCE↗

Attention Enabled Multi-Agent DRL for Decentralized Volt-VAR Control of Active Distribution System Using PV Inverters and SVCs

This paper proposes attention enabled multi-agent deep reinforcement learning (MADRL) framework for active distribution network decentralized Volt-VAR control. Using the unsupervised clustering, the whole distribution system can be decomposed into several sub-networks according to the voltage and reactive power sensitivity relationships. Then, the distributed control problem of each sub-network is modeled as Markov games and solved by the improved MADRL algorithm, where each sub-network is modeled as an adaptive agent. An attention mechanism is developed to help each agent focus on specific information that is mostly related to the reward. All agents are centrally trained offline to learn the optimal coordinated Volt-VAR control strategy and executed in a decentralized manner to make online decisions with only local information. Compared with other distributed control approaches, the proposed method can effectively deal with uncertainties, achieve fast decision makings, and significantly reduce the communication requirements. Comparison results with model-based and other data-driven methods on IEEE 33-bus and 123-bus systems demonstrate the benefits of the proposed approach.

distribution network↗

Computations of Viscous Flows in Complex Geometries Using Multiblock Grid Systems

Generating high quality, structured, continuous, body-fitted grid systems (multiblock grid systems) for complicated geometries has long been a most labor-intensive and frustrating part of simulating flows in complicated geometries. Recently, new methodologies and software have emerged that greatly reduce the human effort required to generate high quality multiblock grid systems for complicated geometries. These methods and software require minimal input form the user-typically, only information about the topology of the block structure and number of grid points. This paper demonstrates the use of the new breed of multiblock grid systems in simulations of internal flows in complicated geometries. The geometry used in this study is a duct with a sudden expansion, a partition, and an array of cylindrical pins. This geometry has many of the features typical of internal coolant passages in turbine blades. The grid system used in this study was generated using a commercially available grid generator. The simulations were done using a recently developed flow solver, TRAF3D.MB, that was specially designed to use multiblock grid systems.

Steinthorsson, Erlendur↗

SpF: Enabling Petascale Performance for Pseudospectral Dynamo Models

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

Pseudospectral (PS)↗

A new use of high resolution magnetograms

Ground-based solar magnetograms are frequently in error by as much as twenty percent and contribute to the poor correlation between magnetic changes and solar flares. High resolution measurement of the magnetic field component, which is normal to the photosphere and measured at photospheric height, can be used to construct a magnetic flux partition function F. Therefore, dF/dt is an EMF which drives atmospheric currents in reconnecting solar active regions. With a high quality magnetograph, the solar probe can be used to obtain good estimates of F and dF/dt and thereby the energy stored as induced solar atmospheric currents during quiescent interflare periods. Should a flare occur during a favorable observing period, the present method of analysis should show characteristic signatures in F, DF/dt, and especially, in the stored flux computed from dF/dt.

Baum, P. J.↗

The equation of state for stellar envelopes. II - Algorithm and selected results

A free-energy-minimization method for computing the dissociation and ionization equilibrium of a multicomponent gas is discussed. The adopted free energy includes terms representing the translational free energy of atoms, ions, and molecules; the internal free energy of particles with excited states; the free energy of a partially degenerate electron gas; and the configurational free energy from shielded Coulomb interactions among charged particles. Internal partition functions are truncated using an occupation probability formalism that accounts for perturbations of bound states by both neutral and charged perturbers. The entire theory is analytical and differentiable to all orders, so it is possible to write explicit analytical formulas for all derivatives required in a Newton-Raphson iteration; these are presented to facilitate future work. Some representative results for both Saha and free-energy-minimization equilibria are presented for a hydrogen-helium plasma with N(He)/N(H) = 0.10. These illustrate nicely the phenomena of pressure dissociation and ionization, and also demonstrate vividly the importance of choosing a reliable cutoff procedure for internal partition functions.

Mihalas, Dimitri↗

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↗

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.↗