Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “basis sets”

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 901 records · Page 50

Quantitative tests of a steady state theory of solar wind electrons

A comparison is made of IMP 6, 7 and 8 electron data with the predictions of a solar wind electron steady state theory in which the control of transport by the macroscopic interplanetary electric and magnetic fields, as well as elastic Coulomb collisions with solar wind protons and thermal electrons, is assumed. While a ratio of forward to backward phase density for field-aligned extrathermal electrons of 6:1 is predicted, electron distribution measurements within the high speed solar wind show this ratio to be typically about an order of magnitude larger. A set of solar wind bulk speed anticorrelations predicted by the theory on the basis of a larger set of assumptions cannot be found in the IMP electron data set, so that improved agreement may require such modifications of the theory's assumptions as the inclusion of inelastic Coulomb and/or wave electron collisions.

Feldman, W. C.↗

Per-Phase and 3-Phase Optimal Coordination of Directional Overcurrent Relays Using Genetic Algorithm

Penetration of the power grid by renewable energy sources, distributed storage, and distributed generators is becoming increasingly common. Increased utilization of these distributed energy resources (DERs) has given rise to additional protection coordination concerns, particularly where they are utilized in an unbalanced manner or where loading among phases is unbalanced. Digital relays such as the SEL-751 (produced by Schweitzer Engineering Laboratories, Pullman, WA, USA) series have the capability of being set on a per-phase basis. This capability is underutilized in common practice. Additionally, in optimization algorithms for determining relay settings, the time-overcurrent characteristics (TOCs) of relays are generally not treated as variables and are assigned before running the optimization algorithm. In this paper, TOC options themselves are treated as discrete variables to be considered in the optimization algorithm. A mixed integer nonlinear programming problem (MINLP) is set up where the goal is to minimize relay operating times. A genetic algorithm (GA) approach is implemented in MATLAB where two cases are considered. In the first case, the TOC and Time dial setting (TDS) of each relay is set on a three-phase basis. In the second case, per-phase settings are considered. Relay TDSs and TOCs are both considered as simultaneous discrete control variables. Despite the stochastic nature of using per-phase settings for unbalanced systems is found to generally allow for shorter operating times. However, for relatively balanced systems, it is best to use three-phase settings if computation time is of importance.

Matthews, Ronald C.↗

MSG-7: Molecular absorption processes related to the penetration of ultraviolet solar radiation into the middle atmosphere

The information presently available on the absorption cross sections of O2 and O3 with attention to the application of these data in middle atmospheric science is reviewed. The cross sections values reported by different groups are intercompared in tabular form where feasible, and specific values are recommended when there is a basis for preferring a particular set of results over other available data. When no such basis exists, the differences among published cross sections then serve to indicate a range of uncertainty. In these cases the need for additional work is indicated. Specific topics addressed are the absorption of molecular oxygen at Lyman alpha, in the Schumann-Runge continuum, in the Schumann-Runge bands, and in the Herzberg continuum. For ozone, the Hartley and Huggins bands are considered.

Frederick, J. E.↗

Calculation of Coupled Vibroacoustics Response Estimates from a Library of Available Uncoupled Transfer Function Sets

The design and theoretical basis of a new database tool that quickly generates vibroacoustic response estimates using a library of transfer functions (TFs) is discussed. During the early stages of a launch vehicle development program, these response estimates can be used to provide vibration environment specification to hardware vendors. The tool accesses TFs from a database, combines the TFs, and multiplies these by input excitations to estimate vibration responses. The database is populated with two sets of uncoupled TFs; the first set representing vibration response of a bare panel, designated as H(sup s), and the second set representing the response of the free-free component equipment by itself, designated as H(sup c). For a particular configuration undergoing analysis, the appropriate H(sup s) and H(sup c) are selected and coupled to generate an integrated TF, designated as H(sup s +c). This integrated TF is then used with the appropriate input excitations to estimate vibration responses. This simple yet powerful tool enables a user to estimate vibration responses without directly using finite element models, so long as suitable H(sup s) and H(sup c) sets are defined in the database libraries. The paper discusses the preparation of the database tool and provides the assumptions and methodologies necessary to combine H(sup s) and H(sup c) sets into an integrated H(sup s + c). An experimental validation of the approach is also presented.

Smith, Andrew↗

On data set tensions and signatures of new cosmological physics

ABSTRACT Can new cosmic physics be uncovered through tensions amongst data sets? Tensions in parameter determinations amongst different types of cosmological observation, especially the ‘Hubble tension’ between probes of the expansion rate, have been invoked as possible indicators of new physics, requiring extension of the ΛCDM paradigm to resolve. Within a fully Bayesian framework, we show that the standard tension metric gives only part of the updating of model probabilities, supplying a data co-dependence term that must be combined with the Bayes factors of individual data sets. This shows that, on its own, a reduction of data set tension under an extension to ΛCDM is insufficient to demonstrate that the extended model is favoured. Any analysis that claims evidence for new physics solely on the basis of alleviating data set tensions should be considered incomplete and suspect. We describe the implications of our results for the interpretation of the Hubble tension.

Cortês, Marina (ORCID:0000000304853767)↗

3D spherical functional expansion tallies in Serpent 2 Monte Carlo code

This work extends the application of functional expansion tallies to 3D spherical geometries. The 3D Zernike polynomials are set as an orthonormal polynomials basis for the functional reconstruction. The study describes the construction of the complete set of polynomials, a natural expansion of the spherical harmonics polynomials where 3D Zernike moments can be evaluated as a linear combination of the geometrical moments. The 3D Zernike polynomials formulation and the computational approach implemented in Serpent 2 are presented and tested through the Godiva model from the ICSBEP criticality benchmark test cases. The implementation results are in agreement with a reference solution described in a fine-resolution mesh, enhancing also the performance and memory demand. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A proposal to extend our understanding of the global economy

Satellites acquire information on a global and repetitive basis. They are thus ideal tools for use when global scale and analysis over time is required. Data from satellites comes in digital form which means that it is ideally suited for incorporation in digital data bases and that it can be evaluated using automated techniques. The development of a global multi-source data set which integrates digital information is proposed regarding some 15,000 major industrial sites worldwide with remotely sensed images of the sites. The resulting data set would provide the basis for a wide variety of studies of the global economy. The preliminary results give promise of a new class of global policy model which is far more detailed and helpful to local policy makers than its predecessors. The central thesis of this proposal is that major industrial sites can be identified and their utilization can be tracked with the aid of satellite images.

Hough, Robbin R.↗

Composing preconditioners for multiphysics PDE systems with applications to Generalized MHD

New patch smoothers or relaxation techniques are developed for solving linear matrix equations coming from systems of discretized partial differential equations (PDEs). One key linear solver challenge for many PDE systems arises when the resulting discretization matrix has a near null space that has a large dimension, which can occur in generalized magnetohydrodynamic (GMHD) systems. Patch-based relaxation is highly effective for problems when the null space can be spanned by a basis of locally supported vectors. The patch-based relaxation methods that we develop can be used either within an algebraic multigrid (AMG) hierarchy or as stand-alone preconditioners. These patch-based relaxation techniques are a form of well-known overlapping Schwarz methods where the computational domain is covered with a series of overlapping sub-domains (or patches). Patch relaxation then corresponds to solving a set of independent linear systems associated with each patch. In the context of GMHD, we also reformulate the underlying discrete representation used to generate a suitable set of matrix equations. In general, deriving a discretization that accurately approximates the curl operator and the Hall term while also producing linear systems with physically meaningful near null space properties can be challenging. Unfortunately, many natural discretization choices lead to a near null space that includes non-physical oscillatory modes and where it is not possible to span the near null space with a minimal set of locally supported basis vectors. Further discretization research is needed to understand the resulting trade-offs between accuracy, stability, and ease in solving the associated linear systems.

97 MATHEMATICS AND COMPUTING↗

Planetary Transmission Diagnostics

This report presents a methodology for detecting and diagnosing gear faults in the planetary stage of a helicopter transmission. This diagnostic technique is based on the constrained adaptive lifting algorithm. The lifting scheme, developed by Wim Sweldens of Bell Labs, is a time domain, prediction-error realization of the wavelet transform that allows for greater flexibility in the construction of wavelet bases. Classic lifting analyzes a given signal using wavelets derived from a single fundamental basis function. A number of researchers have proposed techniques for adding adaptivity to the lifting scheme, allowing the transform to choose from a set of fundamental bases the basis that best fits the signal. This characteristic is desirable for gear diagnostics as it allows the technique to tailor itself to a specific transmission by selecting a set of wavelets that best represent vibration signals obtained while the gearbox is operating under healthy-state conditions. However, constraints on certain basis characteristics are necessary to enhance the detection of local wave-form changes caused by certain types of gear damage. The proposed methodology analyzes individual tooth-mesh waveforms from a healthy-state gearbox vibration signal that was generated using the vibration separation (synchronous signal-averaging) algorithm. Each waveform is separated into analysis domains using zeros of its slope and curvature. The bases selected in each analysis domain are chosen to minimize the prediction error, and constrained to have the same-sign local slope and curvature as the original signal. The resulting set of bases is used to analyze future-state vibration signals and the lifting prediction error is inspected. The constraints allow the transform to effectively adapt to global amplitude changes, yielding small prediction errors. However, local wave-form changes associated with certain types of gear damage are poorly adapted, causing a significant change in the prediction error. The constrained adaptive lifting diagnostic algorithm is validated using data collected from the University of Maryland Transmission Test Rig and the results are discussed.

Lewicki, David G.↗

The Global Precipitation Climatology Project (GPCP): Results, Status and Future

The Global Precipitation Climatology Project (GPCP) is one of a number of long-term, satellite-based, global analyses routinely produced under the auspices of the World Climate Research Program (WCRP) and its Global Energy and Watercycle EXperiment (GEWEX) program. The research quality analyses are produced a few months after real-time through the efforts of scientists at various national agencies and universities in the U.S., Europe and Japan. The primary product is a monthly analysis of surface precipitation that is globally complete and spans the period 1979-present. There are also pentad analyses for the same period and a daily analysis for the 1997-present period. Although generated with somewhat different data sets and analysis schemes, the pentad and daily data sets are forced to agree with the primary monthly analysis on a grid box by grid box basis. The primary input data sets are from low-orbit passive microwave observations, geostationary infrared observations and surface raingauge information. Examples of research with the data sets are discussed, focusing on tropical (25N-25s) rainfall variations and possible long-term changes in the 28-year (1979-2006) monthly dataset. Techniques are used to discriminate among the variations due to ENSO, volcanic events and possible long-term changes for rainfall over both land and ocean. The impact of the two major volcanic eruptions over the past 25 years is estimated to be about a 5% maximum reduction in tropical rainfall during each event. Although the global change of precipitation in the data set is near zero, a small upward linear change over tropical ocean (0.06 mm/day/l0yr) and a slight downward linear change over tropical land (-0.03 mm/day/l0yr) are examined to understand the impact of the inhomogeneity in the data record and the length of the data set. These positive changes correspond to about a 5% increase (ocean) and 3% increase (ocean plus land) during this time period. Relations between variations in surface temperature and precipitation are analyzed on seasonal to inter-decadal time scales. A new, version 3 of GPCP is being planned to incorporate new satellite information (e.g., TRMM) and provide higher spatial and temporal resolution for at least part of the data record. The goals and plans for that GPCP re-processing will be outlined.

Adler, Robert F.↗

Observing Tropospheric Ozone From Space

The importance of tropospheric ozone embraces a spectrum of relevant scientific issues ranging from local environmental concerns, such as damage to the biosphere and human health, to those that impact global change questions, Such is climate warming. From an observational perspective, the challenge is to determine the tropospheric ozone global distribution. Because its lifetime is short compared with other important greenhouse gases that have been monitored over the past several decades, the distribution of tropospheric ozone cannot be inferred from a relatively small set of monitoring stations. Therefore, the best way to obtain a true global picture is from the use of space-based instrumentation where important spatial gradients over vast ocean expanses and other uninhabited areas can be properly characterized. In this paper, the development of the capability to measure tropospheric ozone from space over the past 15 years is summarized. Research in the late 1980s successfully led to the determination of the climatology of tropospheric ozone as a function of season; more recently, the methodology has improved to the extent where regional air pollution episodes can be characterized. The most recent modifications now provide quasi-global (50 N) to 50 S) maps on a daily basis. Such a data set would allow for the study of long-range (intercontinental) transport of air pollution and the quantification of how regional emissions feed into the global tropospheric ozone budget. Future measurement capabilities within this decade promise to offer the ability to provide Concurrent maps of the precursors to the in situ formation of tropospheric ozone from which the scientific community will gain unprecedented insight into the processes that control global tropospheric chemistry

Fishman, Jack↗

Bayes goes fast: Uncertainty quantification for a covariant energy density functional emulated by the reduced basis method

A covariant energy density functional is calibrated using a principled Bayesian statistical framework informed by experimental binding energies and charge radii of several magic and semi-magic nuclei. The Bayesian sampling required for the calibration is enabled by the emulation of the high-fidelity model through the implementation of a reduced basis method (RBM)—a set of dimensionality reduction techniques that can speed up demanding calculations involving partial differential equations by several orders of magnitude. The RBM emulator we build—using only 100 evaluations of the high-fidelity model—is able to accurately reproduce the model calculations in tens of milliseconds on a personal computer, an increase in speed of nearly a factor of 3,300 when compared to the original solver. Besides the analysis of the posterior distribution of parameters, we present model calculations for masses and radii with properly estimated uncertainties. We also analyze the model correlation between the slope of the symmetry energy L and the neutron skin of 48 Ca and 208 Pb. The straightforward implementation and outstanding performance of the RBM makes it an ideal tool for assisting the nuclear theory community in providing reliable estimates with properly quantified uncertainties of physical observables. Such uncertainty quantification tools will become essential given the expected abundance of data from the recently inaugurated and future experimental and observational facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Training and projecting: A reduced basis method emulator for many-body physics

Here, we present the reduced basis method as a tool for developing emulators for equations with tun able parameters within the context of the nuclear many-body problem. The method uses a basis expansion informed by a set of solutions for a few values of the model parameters and then projects the equations over a well-chosen low-dimensional subspace. We connect some of the results in the eigenvector continuation literature to the formalism of reduced basis methods and show how these methods can be applied to a broad set of problems. As we illustrate, the possible success of the formalism on such problems can be diagnosed beforehand by a principal component analysis. We apply the reduced basis method to the one-dimensional Gross-Pitaevskii equation with a harmonic trap ping potential and to nuclear density functional theory for 48 Ca, achieving speed-ups of more than x150 and x250, respectively, when compared to traditional solvers. The outstanding performance of the approach, together with its straightforward implementation, show promise for its application to the emulation of computationally demanding calculations, including uncertainty quantification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A resolution independent neural operator

The Deep operator network (DeepONet) is a powerful yet simple neural operator architecture that utilizes two deep neural networks to learn mappings between infinite-dimensional function spaces. This architecture is highly flexible, allowing the evaluation of the solution field at any location within the desired domain. However, it imposes a strict constraint on the input space, requiring all input functions to be discretized at the same locations; this limits its practical applications. Here, in this work, we introduce a general framework for operator learning from input–output data with arbitrary number and locations of sensors. This begins by introducing a resolution-independent DeepONet (RI-DeepONet), enabling it to handle input functions that are arbitrarily, but sufficiently finely, discretized. To this end, we propose two dictionary learning algorithms to adaptively learn a set of appropriate continuous basis functions, parameterized as implicit neural representations (INRs), from correlated signals defined on arbitrary point cloud data. These basis functions are then used to project arbitrary input function data as a point cloud onto an embedding space (i.e., a vector space of finite dimensions) with dimensionality equal to the dictionary size, which can be directly used by DeepONet without any architectural changes. In particular, we utilize sinusoidal representation networks (SIRENs) as trainable INR basis functions. The introduced dictionary learning algorithms are then used in a similar way to learn an appropriate dictionary of basis functions for the output function data, which defines a new neural operator architecture referred to as the R esolution I ndependent N eural O perator (RINO). In the RINO, the operator learning task simplifies to learning a mapping from the coefficients of input basis functions to the coefficients of output basis functions. We demonstrate the robustness and applicability of RINO in handling arbitrarily (but sufficiently richly) sampled input and output functions during both training and inference through several numerical examples.

Deep operator network (DeepONet)↗

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The high-power helicon program at DIII-D: gearing up for first experiments

Helicon current drive, also called fast wave current drive in the lower hybrid range of frequencies, has long been regarded as a promising current drive tool for reactor grade plasmas. A newly installed MW-level system at DIII-D will be the first test of this technology in reactor-relevant plasmas, in the sense that full single-pass absorption is expected. A 30-module traveling wave antenna has been installed and optimized in-vessel in early 2020. The linear electromagnetic characteristics of the unloaded module array have been extensively tested both on the bench and in the vessel at instrumentation power levels. Excellent performance has been achieved, ~2% reflected power and ~1.5% dissipated power per module in air, in a 10 MHz band around 476 MHz. Stripline feeds on both ends of the antenna allow either co or counter current drive. The installation of a 1.2 MW klystron and associated high-power electronics was completed in Fall 2020. Commissioning of the antenna is ongoing. An important goal of this experiment is to validate the helicon current drive physics basis using an extensive set of new and upgraded diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cheaper and more noise-resilient quantum state preparation using eigenvector continuation

Subspace methods are powerful, noise-resilient methods that can effectively prepare ground states on quantum computers. The challenge is to get a subspace with a small condition number that spans the states of interest using minimal quantum resources. In this work, we will use eigenvector continuation to build a subspace from the low-lying states of a set of Hamiltonians. The basis vectors are prepared using truncated versions of standard state preparation methods such as imaginary-time evolution (ITE), adiabatic state preparation (ASP), and variational quantum eigensolver. By using these truncated methods combined with eigenvector continuation, we can directly improve upon them, obtaining more accurate ground-state energies at a reduced cost. We use several spin systems to demonstrate convergence even when methods like ITE and ASP fail, such as ASP in the presence of level crossings and ITE with vanishing energy gaps. We also showcase the noise resilience of this approach beyond the gains already made by having shallower quantum circuits. Furthermore, our findings suggest that eigenvector continuation can be used to improve existing state preparation methods in the near term.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of PV Module Packaging Strategies of Monofacial and Bifacial PERC Using Degradation Pathway Network Modeling

As the PV industry is rapidly expanding, it is important to thoroughly investigate the long-term impact of packaging strategies on the performance of PV modules. In this study, the variants in sets differ on the basis of manufacturer (A/B), encapsulant (EVA/POE), rear encapsulant (UV-cutoff/opaque/transparent), module architecture (GB/DG) and cell type (monofacial/bifacial). The minimodules were exposed for 2520 hours in modified damp heat, with or without full spectrum light. Every 504 hours, stepwise electrical characterization techniques were employed to track changes in minimodules. Degradation pathway modeling using network structural equation modeling was employed to study pairwise relationships between variables and service lifetime prediction in minimodules. Through this study, differences in quality control are identified in minimodules made by different manufacturers. Minimodules with UV-cutoff rear encapsulant show relatively better stability, whereas the ones with opaque rear encapsulant show greater power loss. In addition, GB having UV-cutoff rear encapsulation and GB with POE having opaque rear encapsulation were identified to be stable as they lack a best model fit. Here, the primary power loss mechanism in degrading variants is interconnect corrosion.

14 SOLAR ENERGY↗