Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equivalence Method”

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 451 records · Page 25

Integrating Terrain Maps Into a Reactive Navigation Strategy

An improved method of processing information for autonomous navigation of a robotic vehicle across rough terrain involves the integration of terrain maps into a reactive navigation strategy. Somewhat more precisely, the method involves the incorporation, into navigation logic, of data equivalent to regional traversability maps. The terrain characteristic is mapped using a fuzzy-logic representation of the difficulty of traversing the terrain. The method is robust in that it integrates a global path-planning strategy with sensor-based regional and local navigation strategies to ensure a high probability of success in reaching a destination and avoiding obstacles along the way. The sensor-based strategies use cameras aboard the vehicle to observe the regional terrain, defined as the area of the terrain that covers the immediate vicinity near the vehicle to a specified distance a few meters away.

Howard, Ayanna↗

M-MRAC Backstepping for Systems with Unknown Virtual Control Coefficients

The paper presents an over-parametrization free certainty equivalence state feedback backstepping adaptive control design method for systems of any relative degree with unmatched uncertainties and unknown virtual control coefficients. It uses a fast prediction model to estimate the unknown parameters, which is independent of the control design. It is shown that the system's input and output tracking errors can be systematically decreased by the proper choice of the design parameters. The benefits of the approach are demonstrated in numerical simulations.

Certainty equivalence control↗

A comparison between protein crystals grown with vapor diffusion methods in microgravity and protein crystals using a gel liquid-liquid diffusion ground-based method

Crystals of human serum albumin have been successfully grown in a variety of gels using crystallization conditions otherwise equivalent to those utilized in the popular hanging-drop vapor-equilibrium method. Preliminary comparisons of gel grown crystals with crystals grown by the vapor diffusion method via both ground-based and microgravity methods indicate that crystals superior in size and quality may be grown by limiting solutal convection. Preliminary X-ray diffraction statistics are presented.

Miller, Teresa Y.↗

Method for calculating alloy energetics

A semiempirical method for the computation of alloy energies is introduced. It is based on the equivalent-crystal theory of defect-formation energies in elemental solids. The method is both simple and accurate. Heats of formation as a function of composition are computed for some binary alloys of Cu, Ni, Al, Ag, Pd, Pt, and Au using the heats of solution in the dilute limit as experimental input. The separation of heats into strain and chemical components helps in understanding the energetics. In addition, lattice-parameter contractions seen in solid solutions of Ag and Au are accurately predicted. Good agreement with experiment is obtained in all cases.

Bozzolo, Guillermo↗

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗

Models of the hard X-ray spectrum of AM Herculis and implications for the accretion rate

Phenomenological fits to the hard X-ray spectrum of AM Hercules have failed to explain the high equivalent width (0.8 + or - 0.1 keV) of Fe K-alpha emission. The scattering and fluorescent line production in a cold or partially ionized accretion column of hard X-rays emitted at the base are investigated using the Monte Carlo method. It is shown that the Fe emission equivalent width observed in the spectrum of AM Her can be explained without resorting to an overabundance of Fe in the accreting material. For sufficient optical depth across the accretion column, equivalent widths of K-alpha emission even larger than the observed values can be obtained. The accretion rate and gravitational power can be deduced from the optical depth across the column, if the column size is known, and, together with the observed hard X-ray and polarized light luminosities, imply a lower limit for the luminosity in the UV to soft X-ray range, for which the observations give model-dependent values.

Swank, J. H.↗

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗

Generative Representations for Automated Design of Robots

A method of automated design of complex, modular robots involves an evolutionary process in which generative representations of designs are used. The term generative representations as used here signifies, loosely, representations that consist of or include algorithms, computer programs, and the like, wherein encoded designs can reuse elements of their encoding and thereby evolve toward greater complexity. Automated design of robots through synthetic evolutionary processes has already been demonstrated, but it is not clear whether genetically inspired search algorithms can yield designs that are sufficiently complex for practical engineering. The ultimate success of such algorithms as tools for automation of design depends on the scaling properties of representations of designs. A nongenerative representation (one in which each element of the encoded design is used at most once in translating to the design) scales linearly with the number of elements. Search algorithms that use nongenerative representations quickly become intractable (search times vary approximately exponentially with numbers of design elements), and thus are not amenable to scaling to complex designs. Generative representations are compact representations and were devised as means to circumvent the above-mentioned fundamental restriction on scalability. In the present method, a robot is defined by a compact programmatic form (its generative representation) and the evolutionary variation takes place on this form. The evolutionary process is an iterative one, wherein each cycle consists of the following steps: 1. Generative representations are generated in an evolutionary subprocess. 2. Each generative representation is a program that, when compiled, produces an assembly procedure. 3. In a computational simulation, a constructor executes an assembly procedure to generate a robot. 4. A physical-simulation program tests the performance of a simulated constructed robot, evaluating the performance according to a fitness criterion to yield a figure of merit that is fed back into the evolutionary subprocess of the next iteration. In comparison with prior approaches to automated evolutionary design of robots, the use of generative representations offers two advantages: First, a generative representation enables the reuse of components in regular and hierarchical ways and thereby serves a systematic means of creating more complex modules out of simpler ones. Second, the evolved generative representation may capture intrinsic properties of the design problem, so that variations in the representations move through the design space more effectively than do equivalent variations in a nongenerative representation. This method has been demonstrated by using it to design some robots that move, variously, by walking, rolling, or sliding. Some of the robots were built (see figure). Although these robots are very simple, in comparison with robots designed by humans, their structures are more regular, modular, hierarchical, and complex than are those of evolved designs of comparable functionality synthesized by use of nongenerative representations.

Homby, Gregory S.↗

Deep compressed seismic learning for fast location and moment tensor inferences with natural and induced seismicity

Fast detection and characterization of seismic sources is crucial for decision-making and warning systems that monitor natural and induced seismicity. However, besides the laying out of ever denser monitoring networks of seismic instruments, the incorporation of new sensor technologies such as Distributed Acoustic Sensing (DAS) further challenges our processing capabilities to deliver short turnaround answers from seismic monitoring. In response, this work describes a methodology for the learning of the seismological parameters: location and moment tensor from compressed seismic records. In this method, data dimensionality is reduced by applying a general encoding protocol derived from the principles of compressive sensing. The data in compressed form is then fed directly to a convolutional neural network that outputs fast predictions of the seismic source parameters. Thus, the proposed methodology can not only expedite data transmission from the field to the processing center, but also remove the decompression overhead that would be required for the application of traditional processing methods. An autoencoder is also explored as an equivalent alternative to perform the same job. We observe that the CS-based compression requires only a fraction of the computing power, time, data and expertise required to design and train an autoencoder to perform the same task. Implementation of the CS-method with a continuous flow of data together with generalization of the principles to other applications such as classification are also discussed.

54 ENVIRONMENTAL SCIENCES↗

The Precipitation Imaging Package: Assessment of Microphysical and Bulk Characteristics of Snow

Remote-sensing observations are needed to estimate the regional and global impacts of snow. However, to retrieve accurate estimates of snow mass and rate, these observations require augmentation through additional information and assumptions about hydrometeor properties. The Precipitation Imaging Package (PIP) provides information about precipitation characteristics and can be utilized to improve estimates of snowfall rate and accumulation. Here, the goal is to demonstrate the quality and utility of two higher-order PIP-derived products: liquid water equivalent snow rate and an approximation of volume-weighted density called equivalent density. Accuracy of the PIP snow rate and equivalent density is obtained through intercomparison with established retrieval methods and through evaluation with colocated ground-based observations. The results confirm the ability of the PIP-derived products to quantify properties of snow rate and equivalent density, and demonstrate that the PIP produces physically realistic snow characteristics. When compared to the National Weather Service (NWS) snow field measurements of six-hourly accumulation, the PIP-derived accumulations were biased only +2.48% higher. Additionally, this work illustrates fundamentally different microphysical and bulk features of low and high snow-to-liquid ratio events, through assessment of observed particle size distributions, retrieved mass coefficients, and bulk properties. Importantly, this research establishes the role that PIP observations and higher-order products can serve for constraining microphysical assumptions in ground-based and spaceborne remotely sensed snowfall retrievals.

Claire Pettersen↗

Further comparison of two synoptic surface wind and pressure analysis methods

The geostrophic stream and potential functions of Sangster and the derived flat pressure field of Pielke and Cram are further compared. A simple numerical experiment with an idealized mountain-atmosphere system compares the two methods and their sensitivity to lateral boundary conditions. The flat pressure field of Pielke and Cram is essentially equivalent to the stream-function of Sangster, although the two methods have different boundary condition formulations and resulting different sensitivities to these boundary conditions.

Cram, Jennifer M.↗

Nonlinear Rescaling and Proximal-Like Methods in Convex Optimization

The nonlinear rescaling principle (NRP) consists of transforming the objective function and/or the constraints of a given constrained optimization problem into another problem which is equivalent to the original one in the sense that their optimal set of solutions coincides. A nonlinear transformation parameterized by a positive scalar parameter and based on a smooth scaling function is used to transform the constraints. The methods based on NRP consist of sequential unconstrained minimization of the classical Lagrangian for the equivalent problem, followed by an explicit formula updating the Lagrange multipliers. We first show that the NRP leads naturally to proximal methods with an entropy-like kernel, which is defined by the conjugate of the scaling function, and establish that the two methods are dually equivalent for convex constrained minimization problems. We then study the convergence properties of the nonlinear rescaling algorithm and the corresponding entropy-like proximal methods for convex constrained optimization problems. Special cases of the nonlinear resealing algorithm are presented. In particular a new class of exponential penalty-modified barrier functions methods is introduced.

Polyak, Roman↗

Spectra-to-exposure conversion using polynomial response models for gamma-ray field characterization

Accurate measurement of exposure rate from gamma-ray spectral data remains a critical challenge during radiological emergency response operations. Conventional methods rely on pre-defined static conversion factors derived from fixed geometries and isotopic compositions, which often fail to capture real-world environmental variability. This study presents a generalized approach as a "next-step" for converting gamma-ray spectral data into exposure rate using polynomial response models. The method introduces a flexible weighting scheme based on the in-situ detector response to distributed sources, enabling a pathway towards improved correspondence between measured spectra and "ground-truth" exposure rates. Experimental data from sodium iodide NaI(Tl) detectors were used to validate the approach as, at least equivalent to the current count-to-exposure method employed in emergency response CONOPS. Results show that the polynomial weighting model is sufficiently equal to the count-to-exposure method and may help improve accuracy given its adaptability to real-world conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗

Optimal processing of satellite-derived magnetic anomaly data

It is shown how the concept of the power spectrum can be extended to two-dimensional spatial power spectra and how it can be used in determining optimal data processing methods for satellite-derived magnetic anomaly data and planning missions to obtain such data. The analysis techniques are applied to the data set and data-processing procedure described by Mayhew et al. (1980), a study that treats magnetic anomaly data for Australia and the surrounding ocean obtained by the polar orbit POGO series satellites. It is shown that the data-processing method used by Mayhew et al. is approximately equivalent to an invariant two-dimensional linear filter and that it is reasonably close to optimal with respect to accuracy, although some possible improvements are suggested. However, as is usual when filtering data, some real 'signal' is unavoidably removed along with the 'noise' resulting in errors that can be quite large. A method for reducing these errors by using additional data from a medium inclination orbit satellite (for example, 60 deg inclination) is proposed.

Mcleod, M. G.↗

Structural response reconstruction using a system-equivalent singular vector basis

Here, this paper develops a novel method for reconstructing the full-field response of structural dynamic systems using sparse measurements. The singular value decomposition is applied to a frequency response matrix relating the structural response to physical loads, base motion, or modal loads. The left singular vectors form a non-physical reduced basis that can be used for response reconstruction with far fewer sensors than existing methods. The contributions of the singular vectors to measured response are termed singular-vector loads (SVLs) and are used in a regularized Bayesian framework to generate full-field response estimates and confidence intervals. The reconstruction framework is applicable to the estimation of single data records and power spectral densities from multiple records. Reconstruction is successfully performed in configurations where the number of SVLs to identify is less than, equal to, and greater than the number of sensors used for reconstruction. In a simulation featuring a seismically excited shear structure, SVL reconstruction significantly outperforms modal FRF-based reconstruction and successfully estimates full-field responses with as few as two uniaxial accelerometers. SVL reconstruction is further verified in a simulation featuring an acoustically excited cylinder. Finally, response reconstruction and uncertainty quantification are performed on an experimental structure with three shaker inputs and 27 triaxial accelerometer outputs.

42 ENGINEERING↗

Nitrate Removal Across Ecogeomorphic Zones in Wax Lake Delta, Louisiana (USA)

Human activities have increased nitrate export from rivers, degrading coastal water quality. At deltaic river mouths, the flow of water through wetlands increases nitrate removal, and the spatial organization of removal rates influences coastal water quality. To understand the spatial distribution of nitrate removal in a river-dominated delta, in this study we deployed 23 benthic chambers across ecogeomorphic zones with varying elevation, vegetation, and sediment properties in Wax Lake Delta (Louisiana, USA) in June 2018. Regression analyses indicate that normalized difference vegetation index is a useful predictor of summertime nitrate removal. Mass transfer velocity were approximately three times greater on a vegetated submerged levee (13 mm hr –1 ), where normalized difference vegetation index was greatest, compared to other locations (4.6 mm hr –1 ). Two methods were developed to upscale nitrate removal across the delta. The flooded-delta method integrates spatially explicit potential removal rates across submerged portions of the delta and suggests that intermediate elevations on the delta—including submerged levees—are responsible for 70% of potential nitrate removal despite covering only 33% of the flooded area. The channel network method treats the delta as a network of river channels and suggests that although secondary channels are more efficient than primary channels at removing received nitrate, primary channels collectively contribute more to overall removal because they convey more of the total nitrate load. The two upscaling methods predict similar rates of nitrate removal, equivalent to less than 4% of nitrate entering the delta. To protect coastal waters against high nitrate loads, management policies should aim to reduce upstream nutrient loads.

54 ENVIRONMENTAL SCIENCES↗

Investigation of experimental observables in search of the chiral magnetic effect in heavy-ion collisions in the STAR experiment *

The chiral magnetic effect (CME) is a novel transport phenomenon, arising from the interplay between quantum anomalies and strong magnetic fields in chiral systems. In high-energy nuclear collisions, the CME may survive the expansion of the quark-gluon plasma fireball and be detected in experiments. Over the past two decades, experimental searches for the CME have attracted extensive interest at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). The main goal of this study is to investigate three pertinent experimental approaches: the $γ$ correlator, the R correlator, and the signed balance functions. We exploit simple Monte Carlo simulations and a realistic event generator (EBE-AVFD) to verify the equivalence of the core components among these methods and to ascertain their sensitivities to the CME signal and the background contributions for the isobar collisions at the RHIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗