Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convergence analysis”

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 253 records · Page 14

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗

Evaluation of top-down and bottom-up global terrestrial respiration estimates and their mismatch with model simulations

Terrestrial respiration is one of the most poorly understood processes in the global carbon cycle, making respiration predictions uncertain. However, expanding observations and machine learning approaches have led to a proliferation of estimates. We compiled total ecosystem and heterotrophic respiration estimates derived from top-down atmospheric inversions and bottom-up upscaling of ecosystem observations and compared them with dynamic vegetation models (DGVM) simulations over the 1980-2020 period. Our analysis revealed a convergence in mean annual global total ecosystem respiration estimates between top-down 97.1 (± SD 6.8) PgC yr-1 and bottom-up 98.5 (+/-13.4) PgC yr-1, which were both significantly lower than the ensemble mean from DGVMs estimates 133.7 (±4.7) PgC yr-1. We also found similar temporal trends between top-down estimates with a mean of 0.075 (±0.05) PgC yr-2, and bottom-up estimates of 0.05 (±0.05) PgC yr-2, which were 5 to 7 times smaller than the ensemble mean trend of 0.34 PgC yr-2 simulated by DGVMs. Global heterotrophic respiration showed much less agreement, ranging from top-down estimates of 42.7 (±4.0) PgC yr-1 to bottom-up estimates of 51.5 (±4.0) PgC yr-1 and a significantly larger ensemble model mean estimate of 60.8 PgC yr-1 (±1.9). The temporal trends in observation-based bottom-up estimates of heterotrophic respiration of 0.03 PgC yr-2 were five times lower than the model ensemble mean 0.15 PgC yr-2. Large regional disagreements in heterotrophic respiration estimates and simulations were evident in tropical and boreal latitudes. Therefore, improved regional and heterotrophic respiration estimates are necessary to reduce uncertainties regarding the future vulnerability of soil carbon.

Ballantyne, Ashley↗

Dynamical evidence for causality between galactic cosmic rays and interannual variation in global temperature

As early as 1959, it was hypothesized that an indirect link between solar activity and climate could be mediated by mechanisms controlling the flux of galactic cosmic rays (CR). Although the connection between CR and climate remains controversial, a significant body of laboratory evidence has emerged at the European Organization for Nuclear Research and elsewhere, demonstrating the theoretical mechanism of this link. In this article, we present an analysis based on convergent cross mapping, which uses observational time series data to directly examine the causal link between CR and year-to-year changes in global temperature. Despite a gross correlation, we find no measurable evidence of a causal effect linking CR to the overall 20th-century warming trend. Furthermore, on short interannual timescales, we find a significant, although modest, causal effect between CR and short-term, year-to-year variability in global temperature that is consistent with the presence of nonlinearities internal to the system. Thus, although CR do not contribute measurably to the 20th-century global warming trend, they do appear as a nontraditional forcing in the climate system on short interannual timescales.

54 ENVIRONMENTAL SCIENCES↗

7 Li( d , p ) 8 Li transfer reaction in the NCSM/RGM approach

Recently, we applied an ab initio method, the no-core shell model combined with the resonating group method, to the transfer reactions with light p-shell nuclei as targets and deuteron as the projectile. In particular, we studied the elastic scattering of deuterium on 7Li and the 7Li(d,p)8Li transfer reaction starting from a realistic two-nucleon interaction. In this contribution, we review of our main results on the 7Li(d,p)8Li transfer reaction, and we extend the study of the relevant reaction channels, by showing the dominant resonant phase shifts of the scattering matrix. We assess also the impact of the polarization effects of the deuteron below the breakup on the positive-parity resonant states in the reaction. For this purpose, we perform an analysis of the convergence trend of the phase and eigenphase shifts, with respect to the number of deuteron pseudostates included in the model space

Physics↗

Magnetic structure of few-nucleon systems at high momentum transfers in a chiral effective field theory approach

The five low-energy constants (LECs) in the electromagnetic current derived in chiral effective field theory ( χ EFT ) up to one loop are determined by a simultaneous fit to the A = 2 – 3 nuclei magnetic moments and to the deuteron magnetic form factor and threshold electrodisintegration at backward angles over a wide range of momentum transfers. The resulting parametrization then yields predictions for the 3 He/ 3 H magnetic form factors in excellent accord with the experimental values for momentum transfers ranging up to ≈0.8 GeV/c, beyond the expected regime of validity of the χEFT approach. The calculations are based on last-generation two-nucleon interactions including high orders in the chiral expansion and derived by Entem, Machleidt, and Nosyk [Phys. Rev. C 96, 024004 (2017)] and by Piarulli et al. [Phys. Rev. C 94, 054007 (2016)], using different χEFT formulations. In the A = 3 calculations, (chiral) three-nucleon interactions are also accounted for. The model dependence resulting from these different formulations of the interactions is found to be mild for momentum transfer below ≈0.8 GeV/c . Further, an analysis of the convergence of the chiral expansion is also provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Privacy-Preserving Robust Consensus for Distributed Microgrid Control Applications

Consensus-based distributed control has been proposed for coordinating distributed energy resources (DERs) in microgrids (MGs). As one key component, distributed average observers are used to estimate the average of a group of reference signals (e.g., voltage, current, or power). State-of-the-art distributed average observers could lead to loss of privacy due to information exchange on the communication channels. The DERs' reference signals, which contain private information, could be inferred by an eavesdropper. In this article, a privacy-preserving distributed average observer is proposed that is based on robust consensus and uses the state decomposition method to preserve privacy. Compared to the existing methods, the proposed observer does not require the knowledge of the reference signal's derivative and gives accurate and smooth estimation, and is thus applicable for MG distributed control applications. A detailed analysis regarding the convergence and privacy properties of the proposed observer is presented. Here, the proposed observer is implemented on hardware controllers and validated in the context of distributed MG control applications through hardware-in-the-loop (HIL) tests.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Final Scientific/Technical Report

In this work, we have built upon our results from previous DOE funding (DEFG 0204ER25655), where we developed new and more efficient methods for solving certain optimization problems with many inequality constraints. This past work resulted in efficient algorithms (and analysis of their convergence) for linear programming, convex quadratic programming, and the training of support vector machines. The algorithms are based on using constraint reduction in interior point methods: at each iteration we consider only a smaller subset of the inequality constraints, focusing on the constraints that are close enough to be relevant. Surprisingly, we have been able to show theoretically that such algorithms are globally convergent and to demonstrate experimentally that they are much more efficient than standard interior point methods.

97 MATHEMATICS AND COMPUTING↗

A Low Order Mimetic Finite Difference Method for Resistive Magnetohydrodynamics in 2D

We developed a low order Mimetic Finite Difference (MFD) discretization for the equations of magnetohydrodynamics (MHD) in two space dimensions. These equations describe the evolution of the electric and magnetic fields in the presence of prescribed velocity field.The method is designed to work on general polygonal meshes and preserves the divergence-free condition on the magnetic field. The electric field is discretized at the vertices/nodes and the magnetic field uses edge-based discretization. The method reconstructs the magnetic field to extract nodal values necessary to approximate some terms present in Ohm’s law. We test the robustness of our numerical scheme on three different types of meshes: with triangular elements, quadrilateral and unstructured polyhedrons obtain from a Voronoi tesselation. Analysis of the convergence for each of the aforementioned mesh types is presented. We finish with a test problem that shows the method is capable of modelling magnetic reconnection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Variational Quantum Algorithms for Semidefinite Programming

A semidefinite program (SDP) is a particular kind of convex optimization problem with applications in operations research, combinatorial optimization, quantum information science, and beyond. In this work, we propose variational quantum algorithms for approximately solving SDPs. For one class of SDPs, we provide a rigorous analysis of their convergence to approximate locally optimal solutions, under the assumption that they are weakly constrained (i.e., N " M, where N is the dimension of the input matrices and M is the number of constraints). We also provide algorithms for a more general class of SDPs that requires fewer assumptions. Finally, we numerically simulate our quantum algorithms for applications such as MaxCut, and the results of these simulations provide evidence that convergence still occurs in noisy settings.

97 MATHEMATICS AND COMPUTING↗

Variational Quantum Algorithms for Semidefinite Programming

A semidefinite program (SDP) is a particular kind of convex optimization problem with applications in operations research, combinatorial optimization, quantum information science, and beyond. In this work, we propose variational quantum algorithms for approximately solving SDPs. For one class of SDPs, we provide a rigorous analysis of their convergence to approximate locally optimal solutions, under the assumption that they are weakly constrained (i.e., N$\gg$M, where N is the dimension of the input matrices and M is the number of constraints). We also provide algorithms for a more general class of SDPs that requires fewer assumptions. Finally, we numerically simulate our quantum algorithms for applications such as MaxCut, and the results of these simulations provide evidence that convergence still occurs in noisy settings.

97 MATHEMATICS AND COMPUTING↗

Analysis and testing of two-dimensional slot nozzle ejectors with variable area mixing sections

Finite difference computer techniques have been used to calculate the detailed performance of air to air two-dimensional ejectors with symmetric variable area mixing sections and coaxial coverging primary nozzles. The analysis of the primary nozzle assumed correct expansion of the flow and is suitable for subsonic and slightly supersonic velocity levels. The variation of the mixing section channel walls is assumed to be gradual so that the static pressure can be assumed uniform on planes perpendicular to the axis. A test program was run to provide two-dimensional ejector test data for verification of the computer analysis. A primary converging nozzle with a discharge geometry of 0.125 inch x 8.0 inch was supplied with 600 SCFM of air at about 35 psia and 180 F. This nozzle was combined with two mixing section geometries with throat sizes of 1.25 inch x 8.0 inch and 1.875 inch x 8.0 inch and was tested at a total of 11 operating points.

Gilbert, G. B.↗

Analysis and testing of two-dimensional vented Coanda ejectors with asymmetric variable area mixing sections

The analysis of asymmetric, curved (Coanda) ejector flow has been completed using a finite difference technique and a quasi-orthogonal streamline coordinate system. The boundary layer type jet mixing analysis accounts for the effect of streamline curvature in pressure gradients normal to the streamlines and on eddy viscosities. The analysis assured perfect gases, free of pressure discontinuities and flow separation and treated three compound flows of supersonic and subsonic streams. Flow parameters and ejector performance were measured in a vented Coanda flow geometry for the verification of the computer analysis. A primary converging nozzle with a discharge geometry of 0.003175 m x 0.2032 m was supplied with 0.283 cu m/sec of air at about 241.3 KPa absolute stagnation pressure and 82 C stagnation temperature. One mixing section geometry was used with a 0.127 m constant radius Coanda surface. Eight tests were run at spacing between the Coanda surface and primary nozzle 0.01915 m and 0.318 m and at three angles of Coanda turning: 22.5 deg, 45.0 deg, and 75.0 deg. The wall static pressures, the loci of maximum stagnation pressures, and the stagnation pressure profiles agree well between analytical and experimental results.

Maroti, L. A.↗

Single scattering from nonspherical Chebyshev particles: A compendium of calculations

A large set of exact calculations of the scattering from a class of nonspherical particles known as Chebyshev particles' has been performed. Phase function and degree of polarization in random orientation, and parallel and perpendicular intensities in fixed orientations, are plotted for a variety of particles shapes and sizes. The intention is to furnish a data base against which both experimental data, and the predictions of approximate methods, can be tested. The calculations are performed with the widely-used Extended Boundary Condition Method. An extensive discussion of this method is given, including much material that is not easily available elsewhere (especially the analysis of its convergence properties). An extensive review is also given of all extant methods for nonspherical scattering calculations, as well as of the available pool of experimental data.

Wiscombe, W. J.↗

Data integration in multi-sensor based robotic workstations

A relaxation labeling algorithm is developed. The major advantage of this algorithm over the existing ones is that the mathematic operation is simplified. The simplification eases the analysis of the convergence properties. Both the theoretical and application aspects of the proposed algorithm are investigated. The local convergence properties of a labeling process with n labels and m labels are established. The investigation of the interaction among the nodes in a multinode labeling process reveals some insight into the mathematical issues involved in the relaxation operations.

Chen, Qin↗

Fatigue Characterization of Alloy 10: a 1300F Disk Alloy for Small Gas Turbine Engines

A detailed fatigue characterization of Alloy 10, a high strength nickel-based disk alloy, was conducted on test coupons machined from a 'pancake' forging. Smooth bar, strain controlled fatigue testing at various R-ratios was run at representative bore, 750 F, and rim, 1300 F, temperatures. This was followed by notch fatigue testing (Kt=2) run under load control. Analysis of the fatigue data using a Smith-Watson-Topper approach and finite element analysis of the notch root was employed to understand material behavior in these tests. Smooth bar fatigue data showed a significant R-ratio dependence at either test temperature which could be accounted for using a Smith-Watson-Topper parameter (SWT). In general, fatigue life was longer at 750 F than 1300 F for a given SWT. For notch fatigue tests, life was longer at 750 F than 1300 F but only at higher stresses. This was attributed to differences in alloy strength. At lower stresses, finite element analysis suggested that convergence of fatigue life at both temperatures resulted from relaxation of stresses at the notch root in the 1300 F tests.

Gayda, John↗

Validation of Vehicle Panel/Equipment Response from Diffuse Acoustic Field Excitation Using Spatially Correlated Transfer Function Approach

An approach for predicting the vibration, strain, and force responses of a flight-like vehicle panel assembly to acoustic pressures is presented. Important validation for the approach is provided by comparison to ground test measurements in a reverberant chamber. The test article and the corresponding analytical model were assembled in several configurations to demonstrate the suitability of the approach for response predictions when the vehicle panel is integrated with equipment. Critical choices in the analysis necessary for convergence of the predicted and measured responses are illustrated through sensitivity studies. The methodology includes representation of spatial correlation of the pressure field over the panel surface. Therefore, it is possible to demonstrate the effects of hydrodynamic coincidence in the response. The sensitivity to pressure patch density clearly illustrates the onset of coincidence effects on the panel response predictions.

Smith, Andrew↗