Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Summing”

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 91 records · Page 5

Matrix elements in the coupled-cluster approach - With application to low-lying states in Li

A procedure is suggested for evaluating matrix elements of an operator between wavefunctions in the coupled-cluster form. The use of the exponential ansatz leads to compact exponential expressions also for matrix elements. Algorithms are developed for summing all effects of one-particle clusters and certain chains of two-particle clusters (containing the well-known random-phase approximation as a subset). The treatment of one-particle perturbations in single valence states is investigated in detail. As examples the oscillator strength for the 2s-2p transition in Li as well as the hyperfine structure for the two states are studied and compared to earlier work.

Martensson-Pendrill, Ann-Marie↗

Using Controlled Impact Experiments to Understand the Effects of Cohesive Blocks on the Cratering Process

When considering large-scale impact events that occur in the gravity-regime, smallscale target heterogeneities (on order of the projectile diameter) are anticipated to have a minimal effect on the cratering process and are therefore not considered in scaling relations [1-3]. Such heterogeneities and other details of the target's strength and structure, however, become increasingly influential as the magnitude of the impact decreases [4]. At the scale of craters forming mostly in regolith on the Moon (< a few meters), for example, target properties such as cohesive strength, density, porosity, and friction have more leverage than gravity to influence the resultant crater dimensions [5]. Further complications beyond just target properties can manifest in many forms, such as subsurface layers [6]; the presence of random, competent blocks [7]; buried lenses of impact melt; and/or pockets of regolith compacted or indurated to varying degrees by previous impacts. Summing the total effect of target properties and any additional heterogeneities for integration into scaling relations is not straightforward, as each individual parameter will have unequal, and so far, unquantified, effects on the final crater. Even when craters are produced under a best-case scenario of controlled laboratory conditions, understanding the entire relevant parameter space can be a formidable task, and deconvolution of all the contributing factors can be a complex undertaking. To begin the process of quantifying such effects, we are conducting an extensive series of impact experiments to build a framework for the interpretation of craters observed on the characteristically complex surfaces throughout the solar system. Here we present the first, and simplest, configuration to be examined in this multi-year experimental campaign: one in which solitary, cohesive blocks are placed in the path of the projectile’s trajectory at different burial depths in granular, sand targets. Among the variables to be investigated in this configuration are the influence of block size, compressive strength, and burial depth.

impact crater↗

ESTIMATION OF THE SPACE SHUTTLE ROLLOUT FORCING FUNCTION

The Space Shuttle Vehicle is assembled in the Vertical Assembly Building (VAB) at Kennedy Space Flight Center in Florida. The Vehicle is stacked on a Mobile Launch Platform (MLP) that weighs eight million pounds. A Crawler Transporter (CT) then carries the MLP and the stacked vehicle (12 million pounds total weight) to the launch complex located 5 miles away. This operation is performed at 0.9 mph resulting in a 4.5-hour transport. A recent test was performed to monitor the dynamic environment that was produced during rollout. It was found that the rollout is a harmonic-rich dynamic environment that was previously not understood. This paper will describe work that has been performed to estimate the forcing function that is produced in the transportation process. The rollout analysis team has determined that there are two families of harmonics of the drive train, which excite the system as a function of CT speed. There are also excitation sources, which are random or narrow-band in frequency and are not a function of CT speed. This presentation will discuss the application of the Sum of Weighted Accelerations Technique (SWAT) to further refine this understanding by estimating the forces and moments at the center-of-mass.

James, George H., III↗

Model of Image Artifacts from Dust Particles

A mathematical model of image artifacts produced by dust particles on lenses has been derived. Machine-vision systems often have to work with camera lenses that become dusty during use. Dust particles on the front surface of a lens produce image artifacts that can potentially affect the performance of a machine-vision algorithm. The present model satisfies a need for a means of synthesizing dust image artifacts for testing machine-vision algorithms for robustness (or the lack thereof) in the presence of dust on lenses. A dust particle can absorb light or scatter light out of some pixels, thereby giving rise to a dark dust artifact. It can also scatter light into other pixels, thereby giving rise to a bright dust artifact. For the sake of simplicity, this model deals only with dark dust artifacts. The model effectively represents dark dust artifacts as an attenuation image consisting of an array of diffuse darkened spots centered at image locations corresponding to the locations of dust particles. The dust artifacts are computationally incorporated into a given test image by simply multiplying the brightness value of each pixel by a transmission factor that incorporates the factor of attenuation, by dust particles, of the light incident on that pixel. With respect to computation of the attenuation and transmission factors, the model is based on a first-order geometric (ray)-optics treatment of the shadows cast by dust particles on the image detector. In this model, the light collected by a pixel is deemed to be confined to a pair of cones defined by the location of the pixel s image in object space, the entrance pupil of the lens, and the location of the pixel in the image plane (see Figure 1). For simplicity, it is assumed that the size of a dust particle is somewhat less than the diameter, at the front surface of the lens, of any collection cone containing all or part of that dust particle. Under this assumption, the shape of any individual dust particle artifact is the shape (typically, circular) of the aperture, and the contribution of the particle to the attenuation factor for a given pixel is the fraction of the cross-sectional area of the collection cone occupied by the particle. Assuming that dust particles do not overlap, the net transmission factor for a given pixel is calculated as one minus the sum of attenuation factors contributed by all dust particles affecting that pixel. In a test, the model was used to synthesize attenuation images for random distributions of dust particles on the front surface of a lens at various relative aperture (F-number) settings. As shown in Figure 2, the attenuation images resembled dust artifacts in real test images recorded while the lens was aimed at a white target.

Willson, Reg↗

Quantum interferometry and pathway selectivity in the nonlinear response of photosynthetic excitons

We propose a time–frequency resolved spectroscopic technique which employs nonlinear interferometers to study exciton–exciton scattering in molecular aggregates. A higher degree of control over the contributing Liouville pathways is obtained as compared to classical light. We show how the nonlinear response can be isolated from the orders-of-magnitude stronger linear background by either phase matching or polarization filtering. Both arise due to averaging the signal over a large number of noninteracting, randomly oriented molecules. We apply our technique to the Frenkel exciton model which excludes charge separation for the photosystem II reaction center. We show how the sum of the entangled photon frequencies can be used to select two-exciton resonances, while their delay times reveal the single-exciton levels involved in the optical process.

Science & Technology - Other Topics↗

Distributed memory approaches for robotic neural controllers

The suitability is explored of two varieties of distributed memory neutral networks as trainable controllers for a simulated robotics task. The task requires that two cameras observe an arbitrary target point in space. Coordinates of the target on the camera image planes are passed to a neural controller which must learn to solve the inverse kinematics of a manipulator with one revolute and two prismatic joints. Two new network designs are evaluated. The first, radial basis sparse distributed memory (RBSDM), approximates functional mappings as sums of multivariate gaussians centered around previously learned patterns. The second network types involved variations of Adaptive Vector Quantizers or Self Organizing Maps. In these networks, random N dimensional points are given local connectivities. They are then exposed to training patterns and readjust their locations based on a nearest neighbor rule. Both approaches are tested based on their ability to interpolate manipulator joint coordinates for simulated arm movement while simultaneously performing stereo fusion of the camera data. Comparisons are made with classical k-nearest neighbor pattern recognition techniques.

Jorgensen, Charles C.↗

Performance of a Strapdown Ring Laser Gyro Tetrad Inertial Navigation System in a helicopter flight environment

Results are presented for the first flight demonstration test of a strapdown laser inertial navigation system in a helicopter environment. It is noted that acceptable navigation accuracies were obtained with system alignment times as short as 5 min or less, and that the use of all four available sensors in the navigation algorithm did not significantly improve performance over that achievable by means of only three orthogonal sensors. Failure detection by means of the parity residual from four axes proved successful for failure levels far below the noise threshold of flight control sensing requirements. Summing the gyro parity residual for a period of 13-20 min can yield step-shifted laser gyro sensor bias error detection comparable to the gyro random shift specification.

Carestia, R. A.↗

Ab initio calculations of monopole sum rules: From finite nuclei to infinite nuclear matter

We compute moments of the isoscalar monopole response of 𝑁 = 𝑍 closed-shell nuclei based on chiral nucleon-nucleon plus three-nucleon interactions. We employ the random-phase approximation (RPA) and two ab initio many-body approaches, the in-medium similarity renormalization group (IMSRG) and coupled-cluster theory (CC). In the IMSRG framework, the moments are obtained as ground-state expectation values, whereas in the CC approach, they are evaluated through excited-state calculations. We find good agreement between the IMSRG and CC results across all nuclei studied. RPA provides a reasonable approximation to the correlated methods if the interaction is soft. From the calculated moments, we extract average energies of the monopole response, compute finite-nucleus incompressibilities, and estimate the incompressibility of symmetric nuclear matter by a fit to a leptodermous expansion. Our extrapolated values are lower than those obtained in nuclear-matter calculations with the same interactions, but the values are consistent with phenomenological ranges.

Bonaiti, Francesca [Michigan State Univ., East Lan↗

Convergence acceleration of Monte Carlo many-body perturbation methods by direct sampling

In the Monte Carlo many-body perturbation (MC-MP) method, the conventional correlation-correction formula, which is a long sum of products of low-dimensional integrals, is first recast into a short sum of high-dimensional integrals over electron-pair and imaginary-time coordinates. These high-dimensional integrals are then evaluated by the Monte Carlo method with random coordinates generated by the Metropolis–Hasting algorithm according to a suitable distribution. The latter algorithm, while advantageous in its ability to sample nearly any distribution, introduces autocorrelation in sampled coordinates, which in turn increases the statistical uncertainty of the integrals and thus the computational cost. It also involves wasteful rejected moves and an initial “burn-in” step as well as displays hysteresis. Here, an algorithm is proposed that directly produces a random sequence of electron-pair coordinates for the same distribution used in the MC-MP method, which is free from autocorrelation, rejected moves, a burn-in step, or hysteresis. Furthermore, this direct-sampling algorithm is shown to accelerate second- (MC-MP2) and third-order Monte Carlo many-body perturbation (MC-MP3) calculations by up to 222% and 38%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Television noise-reduction device

System greatly improves signal-to-noise ratio with little or no loss in picture resolution. By storage of luminance component, which is summed with chrominance component, system performs mathematical integration of basically-repetitive television signals. Integration of signals over interval of their repetition causes little change in original signals and eliminates random noise.

Stamps, J. C.↗

3D Deep Learning Joint Inversion of Active Seismic Full Waveform and Passive Seismic Traveltime Data for Reservoir Imaging and Uncertainty Quantification

Here, we present deep learning (DL) networks for three-dimensional (3D) joint inversion of active seismic full waveform and passive seismic traveltime data to image reservoirs and their properties and quantify imaging uncertainties. Active seismic full-waveform data can provide high-resolution monitoring images but are collected only intermittently because of their high acquisition cost. In contrast, passive seismic data can be gathered at relatively low cost between regular active surveys, although their imaging quality can be compromised by factors such as low signal-to-noise ratios and limited ray coverage of the target. Although these datasets are routinely acquired together at CO 2 storage sites, their combined inversion within a 3D DL framework has not been previously demonstrated. To our knowledge, this is the first study to address this gap, combining the strength of both data types. For efficient data storage and DL training with large 3D seismic datasets, we use a 3D data matrix in which a random number of passive seismic traveltime data are stored as parabolic envelopes using one-hot encoding and a 3D full-waveform data matrix in which multiple shot gathers are summed. Two network architectures are evaluated: a single-encoder U-Net for single-data type inversion and a dual-encoder U-Net for joint inversion of active and passive seismic data. We also evaluate the single-encoder U-Net for joint inversion by concatenating full-waveform data and traveltime data. We propose a systematic approach for selecting an optimal dropout rate that balances regularization during training and Monte Carlo dropout-based uncertainty quantification during prediction by examining the correlation coefficient between standard deviation and prediction error, along with the training misfit, across a range of dropout rates. 3D DL inversion experiments include five different network configurations, with evaluations under ideal, noisy and dropout-enabled conditions. Both model and data uncertainties are assessed, as well as their combined effects. Across all conditions, the networks consistently predict accurate CO 2 saturation models with low prediction errors, such as a structural similarity index of 0.993 and CO 2 difference of 1.1%. Uncertainty estimates show strong spatial correlation with prediction errors, confirming the effectiveness of the proposed dropout selection approach. The results demonstrate that our DL approach, utilizing compact data representations and appropriate uncertainty quantification, yields accurate subsurface images under various inversion conditions and provides valuable insights into the reliability of predictions.

Um, Evan Schankee [Lawrence Berkeley National Labo↗

Explaining persistent incomplete mixing in multicomponent reactive transport with Eulerian stochastic model

We present an Eulerian stochastic advection–diffusion–reaction (SADR) model and use it to explain incomplete mixing typically observed in transport experiments with bimolecular reactions. Unlike traditional advection– dispersion–reaction (ADR) models, the SADR model describes mechanical and diffusive mixing as two separate processes. In the SADR model, mechanical mixing is driven by random advective velocity whose variance is given by the coefficient of mechanical dispersion. The diffusive mixing is modeled as a Fickian diffusion process with the effective diffusion coefficient. We demonstrate that the sum of the two coefficients is equal to the dispersion coefficient, but only the effective diffusion coefficient contributes to the mixing-controlled reactions. We use experimental results of Gramling et al. * * to show that for transport and bimolecular reactions in porous media, the SADR model is significantly more accurate than the ADR model, which overestimates the concentration of the reaction product by as much as 60%. We also show that the SADR model predicts an experimentally observed bell- shaped spatial distribution of the reactive product concentration, while the ADR model results in a concentration distribution with an unphysical kink.

Bimolecular reactions, incomplete mixing, Stochast↗

Minimizing Timing Jitter’s Impact on Ground-Penetrating Radar Array Coupling Signals

This article presents a novel investigation into the effect of receiver timing jitter on the quality of imaging of ground-penetrating radar systems. We explicitly show the impact of timing jitter throughout an entire example processing pipeline. The process of coupling removal in systems whose receivers have a random Gaussian-distributed jitter with a high standard deviation leaves randomly distributed residue in the radar image. This residue, called jitter-induced coupling noise, extends even into the expected target domain which may cause false alarms. This residue approximately follows a sum of Gaussian and Gamma distributions based on the higher order derivatives of the coupling signal. Once migrated to create an image, the coupling residue is best described by a Weibull distribution. We show that simple filtering does not remove enough residue and propose two preventative design specifications called the linear slope and migrated probability specifications. Antenna and system designers should use these specifications to assess the true severity of coupling signals in a given jitter environment.

47 OTHER INSTRUMENTATION↗

An investigation of Newton-Sketch and subsampled Newton methods

Sketching, a dimensionality reduction technique, has received much attention in the statistics community. In this paper, we study sketching in the context of Newton's method for solving finite-sum optimization problems in which the number of variables and data points are both large. In this work, we study two forms of sketching that perform dimensionality reduction in data space: Hessian subsampling and randomized Hadamard transformations. Each has its own advantages, and their relative tradeoffs have not been investigated in the optimization literature. Additionally, our study focuses on practical versions of the two methods in which the resulting linear systems of equations are solved approximately, at every iteration, using an iterative solver. The advantages of using the conjugate gradient method vs. a stochastic gradient iteration are revealed through a set of numerical experiments, and a complexity analysis of the Hessian subsampling method is presented.

97 MATHEMATICS AND COMPUTING↗

Hybrid analog-digital associative neural network

Random access memory is used to store synaptic information in the form of a matrix of rows and columns of binary digits. N rows read in sequence are processed through switches and resistors, and a summing amplifier to N neural amplifiers in sequence, one row for each amplifier, using a first array of sample-and-hold devices S/H1 for commutation. The outputs of the neural amplifiers are stored in a second array of sample-and-hold devices S/H2 so that after N rows are processed, all of the second array of sample-and-hold devices are updated. A second memory may be added for binary values of 0 and -1, and processed simultaneously with the first to provide for values of 1, 0, and -1, the results of which are combined in a difference amplifier.

Moopenn, Alexander W.↗

Hybrid analog-digital associative neural network

Random access memory is used to store synaptic information in the form of a matrix of rows and columns of binary digits. N rows read in sequence are processed through switches and resistors, and a summing amplifier to N neural amplifiers in sequence, one row for each amplifier, using a first array of sample-and-hold devices S/H1 for commutation. The outputs of the neural amplifiers are stored in a second array of sample-and-hold devices S/H2 so that after N rows are processed, all of said second array of sample-and-hold devices are updated. A second memory may be added for binary values of 0 and -1, and processed simultaneously with the first to provide for values of 1, 0, and -1, the results of which are combined in a difference amplifier.

Moopenn, Alexander W.↗

Square Root Statistics of Density Matrices and Their Applications

To estimate the degree of quantum entanglement of random pure states, it is crucial to understand the statistical behavior of entanglement indicators such as the von Neumann entropy, quantum purity, and entanglement capacity. These entanglement metrics are functions of the spectrum of density matrices, and their statistical behavior over different generic state ensembles have been intensively studied in the literature. As an alternative metric, in this work, we study the sum of the square root spectrum of density matrices, which is relevant to negativity and fidelity in quantum information processing. In particular, we derive the finite-size mean and variance formulas of the sum of the square root spectrum over the Bures–Hall ensemble, extending known results obtained recently over the Hilbert–Schmidt ensemble.

97 MATHEMATICS AND COMPUTING↗