Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximate computing”

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 1,063 records · Page 59

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interface learning in fluid dynamics: Statistical inference of closures within micro–macro-coupling models

Many complex multiphysics systems in fluid dynamics involve using solvers with varied levels of approximations in different regions of the computational domain to resolve multiple spatiotemporal scales present in the flow. The accuracy of the solution is governed by how the information is exchanged between these solvers at the interface and several methods have been devised for such coupling problems. In this article, we construct a data-driven model by spatially coupling a microscale lattice Boltzmann method (LBM) solver and macroscale finite difference method (FDM) solver for reaction-diffusion systems. The coupling between the micro-macro solvers has one to many mapping at the interface leading to the interface closure problem, and we propose a statistical inference method based on neural networks to learn this closure relation. The performance of the proposed framework in a bifidelity setting partitioned between the FDM and LBM domain shows its promise for complex systems where analytical relations between micro-macro solvers are not available.

42 ENGINEERING↗

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Distribution of blackouts in the power grid and the Motter and Lai model

Carreras, Dobson, and colleagues have studied empirical data on the sizes of the blackouts in real grids and modeled them with computer simulations using the direct current approximation. They have found that the resulting blackout sizes are distributed as a power law and suggested that this is because the grids are driven to the self-organized critical state. In contrast, more recent studies found that the distribution of cascades is bimodal resulting in either a very small blackout or a very large blackout, engulfing a finite fraction of the system. Here we reconcile the two approaches and investigate how the distribution of the blackouts changes with model parameters, including the tolerance criteria and the dynamic rules of failure of the overloaded lines during the cascade. Finally, we study the same problem for the Motter and Lai model and find similar results, suggesting that the physical laws of flow on the network are not as important as network topology, overload conditions, and dynamic rules of failure.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multigluon Correlations and Evidence of Saturation from Dijet Measurements at an Electron-Ion Collider

In this work, we study inclusive and diffractive dijet production in electron-proton and electron-nucleus collisions within the color glass condensate effective field theory. We compute dijet cross sections differentially in both mean dijet transverse momentum P and recoil momentum Δ , as well as the anisotropy in the relative angle between P and Δ . Our results cover a much larger kinematic range than accessible in previous computations performed in the correlation limit approximation, where it is assumed that | P | &raquo | Δ | . We validate this approximation in its range of applicability and quantify its failure for | P | ≲ | Δ | . We also predict significant target-dependent deviations from the correlation limit approximation for | P | > | Δ | and | P | ≲ Q s , which offers a straightforward test of gluon saturation and access to multigluon distributions at a future Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Allocation and Operation of A Hydropneumatic Energy Storage with Building Microgrid

Behind-meter application for a recently invented hydropneumatic Ground-Level Integrated Diverse Energy Storage (GLIDES) is focused in this research. Considering possible uncertainties from solar radiation and electricity load in a building microgrid, two-stage stochastic programming is adopted for 15-min operation and a mixed-integer nonlinear stochastic model is built with integrated energy systems of GLIDES, solar panel and power generating unit. In order to make the model computationally tractable, we use sample average approximate method to decrease associated complexity with sample size and improve solving efficiency. The experimental results in deterministic and stochastic operation indicate that best candidate solution from SAA are superior to the solution from expected value problem in terms of optimality gap and variance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Selecting Critical Scenarios of DER Adoption in Distribution Grids Using Bayesian Optimization

We develop a new methodology to select scenarios of DER adoption most critical for distribution grids. Anticipating risks of future voltage and line flow violations due to additional PV adopters is central for utility investment planning but continues to rely on deterministic or ad hoc scenario selection. We propose a highly efficient search framework based on multi-objective Bayesian Optimization. We treat underlying grid stress metrics as computationally expensive black-box functions, approximated via Gaussian Process surrogates and design an acquisition function based on probability of scenarios being Pareto-critical across a collection of line- and bus-based violation objectives. Our approach provides a statistical guarantee and offers an order of magnitude speed-up relative to a conservative exhaustive search. Case studies on realistic feeders with 200-400 buses demonstrate the effectiveness and accuracy of our approach.

Mulkin, Olivier↗

pnnl/NWHypergraph

NWHypergraph is a C++ hypergraph processing framework for shared-memory architecture. NWHypergraph provides efficient algorithms to construct s-line graphs, a lower-order approximation of a given hypergraph, and computes different graph metrics of a s-line graph such as s-connected components, s-betweenness centrality, s-closeness centrality, etc. It also provides Python APIs for s-line graph computation. The Python APIs are provided using Pybind11

Lumsdaine, Andrew↗

On the modeling of $t\bar{t}W^\pm$ signatures at the LHC

We discuss our recent progress on improving the theoretical description of the pp \to t\bar{t}W^\pm p p → t t ‾ W ± process with particular focus on fiducial signatures that include top-quark decays. We employ a wide range of techniques from parton shower matching, narrow width approximation to full off-shell computations to assess the importance of the different modeling approaches.

Kraus, Manfred↗

Solving Coupled Cluster Equations by the Newton Krylov Method

We describe using the Newton Krylov method to solve the coupled cluster equation. The method uses a Krylov iterative method to compute the Newton correction to the approximate coupled cluster amplitude. The multiplication of the Jacobian with a vector, which is required in each step of a Krylov iterative method such as the Generalized Minimum Residual (GMRES) method, is carried out through a finite difference approximation, and requires an additional residual evaluation. The overall cost of the method is determined by the sum of the inner Krylov and outer Newton iterations. We discuss the termination criterion used for the inner iteration and show how to apply pre-conditioners to accelerate convergence. We will also examine the use of regularization technique to improve the stability of convergence and compare the method with the widely used direct inversion of iterative subspace (DIIS) methods through numerical examples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A note on mass loss during collisions between galaxies and the formation of giant systems.

Stars will be liberated from the outer portions of galaxies which interact tidally. Using models based on the density distribution observed in the E1 galaxy NGC 3379, the mass loss is computed for hyperbolic collisions on an approximate but conservative basis. Binary collisions are then summed to show that large numbers of stars will be evaporated from elliptical members of dense clusters of galaxies. An estimate of one trillion solar masses is obtained for the amount of material dispersed in the Coma cluster, and in rich very dense systems (Morgan's Type I) the large amount of tidally dispersed debris may constitute the centrally located galaxies designated cD systems which often dominate these clusters.

Gallagher, J. S., III↗

The effect of thermal neutron field slagging caused by cylindrical BF3 counters in diffusion media

Computations are carried out in transport approximation (first collision method) for the attenuation of the field of thermal neutrons formed in counters of the CHM-8 and CHMO-5 type. The deflection of the thermal neutron field is also obtained near the counters and in the air (shade effect) and in various decelerating media (water, paraffin, plexiglas) for which the calculations are carried out on the basis of diffusion theory. To verify the calculations, the distribution of the density of the thermal neutrons at various distances from the counter in the water is measured.

Gorshkov, G. V.↗

Rapid methods for calculating radiation transport in the entry environment

A procedure is developed for the prediction of radiation transport events in the context of the entry heating environment. The equivalent-width and exponential approximations are employed to reduce the computational requirements of the governing equations. Novel features are introduced in the use of an Elsasser band model to estimate line overlapping events and in the formulation of the wall reflection events. A matrix of calculations is presented to allow an assessment of the trade-offs between accuracy and computational efforts. It is concluded that the equivalent width model allows significant savings in computational effort with only modest penalties in accuracy. Consequently, the model should be viewed as an attractive candidate for use in radiation-coupled flow field prediction procedures.

Nicolet, W. E.↗

Preliminary design of composite wings for buckling, strength and displacement constraints

An unstiffened panel buckling constraint for balanced, symmetric laminated composites is included on the global design level in a mathematical programming structural optimization procedure for designing wing structures. Constraints are introduced by penalty functions, and Newton's method based on approximate second derivatives of the penalty terms is used as the search algorithm to obtain minimum-mass designs. Constraint approximations used during the optimization process contribute to the computational efficiency of the procedure. A criterion is developed that identifies the appropriate conservative form of the constraint approximations that are used with the optimization procedure. Minimum-mass design results are obtained for a multispar high-aspect-ratio wing subjected to material strength, minimum-gage, displacement, panel buckling and twist constraints. The material systems considered for the examples are all graphite-epoxy, graphite-epoxy with boron-epoxy spar caps, and all aluminum. The composite material designs are shown to have an advantage over the aluminum designs since they can often satisfy additional constraints with only small mass increases.

Starnes, J. H., Jr.↗

Dynamic response and input identification of MDOF structures subjected to coupled random vector inputs

A method of random dynamic analysis is presented which is based on the classical approach applied to a discretized structure. The method assumes that the system identification is available in the form of natural modes and frequencies. These modes and frequencies can be found from either available solutions or approximately from finite element programs. A computer program has been developed to perform the computations required in the analysis. All computations are performed with transformed modal variables, which results in significant economy since the number of modal degrees of freedom is almost always less than the number of physical degrees of freedom. The program computes the response to force and base inputs which are statistically coupled. A method is also presented for predicting one of the inputs if the second input and the acceleration response at a point on the structure are known. Finally, results are presented for the random response of a rectangular plate subjected to a random pressure and a random base input. The inputs are considered individually and with various degrees of statistical coupling.

Ocallahan, J. C.↗

Reducing the error of geoid undulation computations by modifying Stokes' function

The truncation theory as it pertains to the calculation of geoid undulations based on Stokes' integral, but from limited gravity data, is reexamined. Specifically, the improved procedures of Molodenskii et al. are shown through numerical investigations to yield substantially smaller errors than the conventional method that is often applied in practice. In this improved method, as well as in a simpler alternative to the conventional approach, the Stokes' kernel is suitably modified in order to accelerate the rate of convergence of the error series. These modified methods, however, effect a reduction in the error only if a set of low-degree potential harmonic coefficients is utilized in the computation. Consider, for example, the situation in which gravity anomalies are given in a cap of radius 10 deg and the GEM 9 (20,20) potential field is used. Then, typically, the error in the computed undulation (aside from the spherical approximation and errors in the gravity anomaly data) according to the conventional truncation theory is 1.09 m; with Meissl's modification it reduces to 0.41m, while Molodenskii's improved method gives 0.45 m. A further alteration of Molodenskii's method is developed and yields an RMS error of 0.33 m. These values reflect the effect of the truncation, as well as the errors in the GEM 9 harmonic coefficients. The considerable improvement, suggested by these results, of the modified methods over the conventional procedure is verified with actual gravity anomaly data in two oceanic regions, where the GEOS-3 altimeter geoid serves as the basis for comparison. The optimal method of truncation, investigated by Colombo, is extremely ill-conditioned. It is shown that with no corresponding regularization, this procedure is inapplicable.

Jekeli, C.↗

Cyclic plasticity and failure of structural components

An analytical approach for low-cycle fatigue prediction is presented. The approach combines a cyclic plasticity model with the finite element method and a damage accumulation criterion for ductile metals. The cyclic plasticity model is based on the concept of the combination of several yield surfaces. The surfaces are related to the material uniaxial stress-strain curve idealized by piecewise linear segments. The damage criterion is based on the Coffin-Manson formulae modified for the mean stress variation effect. It is extended to the multiaxial varying stress-strain field and applied for both the crack initiation and the crack growth processes. The stable slow crack growth rate is approximated by the damage accumulation gradient computed from the cracked finite element models. This procedure requires fatigue testing data of only smooth specimens under constant strain amplitudes. The present approach is illustrated by numerical examples of an aircraft wing stiffened panel subjected to compression, which causes material yielding and residual tension.

Kalev, I.↗