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 163 records · Page 9

On the vortical-flow-prediction capability of an unstructured-grid Euler solver

The results from a concentrated computational effort are presented with the primary objective being directed at evaluating the vortical-flow-prediction capability of an unstructured-grid Euler solver. Both viscous and inviscid solutions, obtained from an established structured-grid method, along with experimental wind-tunnel data are used as benchmark measures to assess the validity of the unstructured-grid Euler results. Viscous effects on vortical flows are first identified by comparing the viscous and inviscid solutions obtained from the structured-grid method. Computational data analysis is then presented which reveals excellent correlations between the inviscid structured and unstructured-grid results in terms of off-surface flow structures, surface pressure distribution, and the predicted longitudinal aerodynamic characteristics. The sensitivity of the unstructured-grid inviscid solutions to grid refinement is also discussed along with an analysis of the convergence and performance characteristics for each method.

Unstructured grid↗

On the Vortical-Flow Prediction Capability of an Unstructured-Grid Euler Solver

The results from a concentrated computational effort are presented with the primary objective being directed at evaluating the vortical-flow-prediction capability of an unstructured-grid Euler solver. Both viscous and inviscid solutions, obtained from an established structured-grid method, along with an experimental wind-tunnel data are used as bench-mark measures to assess the validity of the unstructured-grid Euler results. Viscous effects on vortical flows are first identified by comparing the viscous and inviscid solutions obtained form the structured-grid method. Computational data analysis are then presented which reveal excellent correlations between the inviscid structured and unstructured-grid results in terms of off-surface flow structures, surface pressure distribution and the predicted longitudinal aerodynamic characteristics. The sensitivity of the unstructured-grid inviscid solutions to grid refinement is also discussed along with an analysis of the convergence and performance characteristics for each method.

Unstructured grid↗

Semigroup approximation and robust stabilization of distributed parameter systems

Theoretical results that enable rigorous statements of convergence and exponential stability of Galerkin approximations of LQR controls for infinite dimensional, or distributed parameter, systems have proliferated over the past ten years. In addition, extensive progress has been made over the same time period in the derivation of robust control design strategies for finite dimensional systems. However, the study of the convergence of robust finite dimensional controllers to robust controllers for infinite dimensional systems remains an active area of research. We consider a class of soft-constrained differential games evolving in a Hilbert space. Under certain conditions, a saddle point control can be given in feedback form in terms of a solution to a Riccati equation. By considering a related LQR problem, we can show a convergence result for finite dimensional approximations of this differential game. This yields a computational algorithm for the feedback gain that can be derived from similar strategies employed in infinite dimensional LQR control design problems. The approach described in this paper also inherits the additional properties of stability robustness common to game theoretic methods in finite dimensional analysis. These theoretical convergence and stability results are verified in several numerical experiments.

Kurdila, A. J.↗

On the convergence of the coupled-wave approach for lamellar diffraction gratings

Among the many existing rigorous methods for analyzing diffraction of electromagnetic waves by diffraction gratings, the coupled-wave approach stands out because of its versatility and simplicity. It can be applied to volume gratings and surface relief gratings, and its numerical implementation is much simpler than others. In addition, its predictions were experimentally validated in several cases. These facts explain the popularity of the coupled-wave approach among many optical engineers in the field of diffractive optics. However, a comprehensive analysis of the convergence of the model predictions has never been presented, although several authors have recently reported convergence difficulties with the model when it is used for metallic gratings in TM polarization. Herein, three points are made: (1) in the TM case, the coupled-wave approach converges much slower than the modal approach of Botten et al; (2) the slow convergence is caused by the use of Fourier expansions for the permittivity and the fields in the grating region; and (3) is manifested by the slow convergence of the eigenvalues and the associated modal fields. The reader is assumed to be familiar with the mathematical formulations of the coupled-wave approach and the modal approach.

Li, Lifeng↗

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↗

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↗

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↗