Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced basis method”

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 415 records · Page 23

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simulating scalar field theories on quantum computers with limited resources

Here we present a quantum algorithm for implementing φ 4 lattice scalar field theory on qubit computers. The field is represented in the discretized field amplitude basis. The number of qubits and elementary gates required by the implementation of the evolution operator is proportional to the lattice size. The algorithm allows efficient φ 4 state preparation for a large range of input parameters in both the normal and broken-symmetry phases. The states are prepared using a combination of variational and adiabatic evolution methods. First, the ground state of a local Hamiltonian, which includes the φ 4 self-interaction, is prepared using short variational circuits. Next, this state is evolved by switching on the coupling between the lattice sites adiabatically. The parameters defining the local Hamiltonian are adjustable and constitute the input of our algorithm. We present a method to optimize these parameters in order to reduce the adiabatic time required for state preparation. For preparing broken-symmetry states, the adiabatic evolution problems caused by crossing the phase transition critical line and by the degeneracy of the broken-symmetry ground state can be addressed using an auxiliary external field which gradually turns off during the adiabatic process. We show that the time dependence of the external field during the adiabatic evolution is important for addressing the broken-symmetry ground state degeneracy. The adiabatic time dependence on the inverse error tolerance can be reduced from quadratic to linear by using a field strength that decreases exponentially in time relative to one that decreases linearly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An Applied Strategy for Using Empirical and Hybrid Models in Online Monitoring

The monitoring of plant equipment for failure prediction is one of the key contributors to operation and maintenance (O&M) costs for a nuclear power plant (NPP) because O&M monitoring depends on labor-intensive activities that are required to meet high equipment reliability standards. These activities rely primarily on humans for information gathering, condition diagnosis, and predictive analysis. Online monitoring aims to automate these activities by relying on sensors to replace human information gathering and machine learning to replace human analysis and decision making. To facilitate automated monitoring, a systematic strategy for anomaly detection is needed to optimally use the available sensor data, empirical models, and physics-supported models. This strategy is essential to provide credible reasoning on why and when an empirical (i.e., purely data-driven) versus hybrid (i.e., physics-supported) approach should be used and to determine the ideal mix of these two approaches for a defined anomaly detection scope. The extant methods usually adopt an ad hoc trial-and-error approach that, in addition to being time-consuming and costly, is also highly subjective; it is impacted by the background and the skill set of the personnel making the decisions. Thus, such an approach cannot guarantee an optimum outcome. This represents the motivation of the current research effort, which is focused on devising a scientifically supported strategy for the optimum selection of anomaly detection methods. This report presents a detailed assessment of the main anomaly detection techniques within the empirical or hybrid method streams. Empirical methods include pattern, statistical, and causal inference. Hybrid methods include the use of physics models to train and test data methods, reduce data dimensionality, reduce data-model complexity, augment data, and reduce empirical uncertainty; hybrid methods also include the use of data to tune physics models. The listed techniques within these two streams represent the vast majority of techniques performed for anomaly detection. Using the techniques as outcomes, a strategy was developed to enable a systematic decision-making process to lead to one of these techniques. The strategy is driven by key decision points related to data relevance, simple modeling feasibility, data inference, physics-modeling value, data dimensionality, physics knowledge, method of validation, performance, data availability and suitability for training and testing, cause-effect, entropy inference, and model fitting. Each of these decision points in the strategy is explained in detail in this report with examples, along with the scientific basis behind the decisions and outcomes in common and simplified terminology. The strategy is developed for use by any NPP staff with basic engineering or science knowledge. A user-friendly graphical state flow diagram was also developed as a visual presentation of the strategy. The strategy was tested and demonstrated through two pilot projects for the application of anomaly detection at an NPP. Each pilot had two use cases: an initial case in which certain decisions were made that resulted in one or more empirical techniques and a revised use case where one or more key decisions were modified resulting in using a set of hybrid methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗

Automated pixel screening and selection technique

One approach for estimating crop acreage on the basis of Landsat data involves the sampling of portions of Landsat imagery which represent 5 x 6-n mi areas on the ground, taking into account acreage estimates for each of these areas. Difficulties concerning the estimating procedures are related to labeling. These difficulties arise in particular in connection with mixed pixels. The present investigation is concerned with the employment of a suitable automated technique for overcoming these difficulties. A description of the selected technique is presented, and the test conducted for the evaluation of the technique is discussed. It is found that the employed automated pixel screening and selection procedure provides representative samples of the scene. The amount of time spent on pure pixel selection is significantly reduced by using the automated method.

Dailey, C. L.↗

Investigation of UV, ns-laser damage resistance of hafnia films produced by electron beam evaporation and ion beam sputtering deposition methods

Laser-induced damage in coating materials with a high index of refraction, such as hafnia, limits the performance of high power and high energy laser systems. Understanding the underlying physics responsible for laser damage holds the key for developing damage-resistant optical films. Previous studies have reported a substantial difference in laser damage onset for hafnia films produced by different deposition methods, yet the underlying mechanisms for the observed difference remain elusive. We combined laser damage testing with analytical characterizations and theoretical simulations to investigate the response of hafnia films produced by electron (e-) beam evaporation vs ion beam sputtering (IBS) methods upon UV ns-laser exposure. We found that e-beam produced hafnia films were overall more damage resistant; in addition, we observed a polarization anisotropy associated with the onset of damage in the e-beam films, while this effect was absent in the latter films. The observed differences can be attributed to the stark contrast in the pressure inside the pores inherent in both films. The high pressure inside the IBS-induced nanobubbles has been shown to reduce the threshold for laser-induced plasma breakdown leading to film damage. The polarization effects in the e-beam coatings can be related to the asymmetric electric field intensification induced by the columnar void structure. Our findings provide a fundamental basis for developing strategies to produce laser damage-resistant coatings for UV pulsed laser applications.

36 MATERIALS SCIENCE↗

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗

Light Water Reactor Sustainability: Program Materials Research Pathway (FY22 Technical Program Plan)

The Materials Research (MR) Pathway within the Light Water Reactor Sustainability (LWRS) program is charged with performing R&D to develop the scientific basis for understanding and predicting the long-term environmental degradation behavior of materials in nuclear reactors. Furthermore, the mechanistic understanding of degradation phenomena in materials must be leveraged to develop mitigation and repair strategies and new material alternatives for existing components. This research will provide data and methods to assess the performance of systems, structures, and components essential to safe and sustained reactor operations. The R&D products developed from the LWRS program will be used by stakeholders—including utilities, industry groups, and regulators—to inform operational and regulatory requirements for materials in reactor systems, structures, and components subjected to long-term operation conditions, providing key inputs to regulators and industry. Therefore, the intent of this research is to provide options to reduce the operating costs, which may be in the form of offset maintenance costs due to better predictive models for component lifetimes, improved analyses of materials through nondestructive evaluation, reduced costs for repairs, or extended performance of plants through the selection of improved replacement materials. To provide the best options, industry experience and guidance are important because of their role in coordinated or collaborative research projects. The objectives of this report are to describe the motivation and organization of the MR Pathway within the LWRS program; provide details on the individual research tasks within the MR Pathway; describe the outcomes and deliverables of the MR Pathway, including recent technical highlights and progress; and describe the requirements for performing this critically important research.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Concentrator hot-spot testing, phase 1

Results of a study to determine the hot-spot susceptibility of concentrator cells, to provide a hot-spot qualification test for concentrator modules, and to provide guidelines for reducing hot-spot susceptibility are presented. Hot-spot heating occurs in a photovoltaic module when the short-circuit current of a cell is lower than the string operating current forcing the cell into reverse bias with a concurrent power dissipation. Although the basis for the concentrator module hot-spot qualification test is the test developed for flat-plate modules, issues, such as providing cell illumination, introduce additional complexities into the testing procedure. The same general guidelines apply for protecting concentrator modules from hot-spot stressing as apply to flat-plate modules. Therefore, recommendations are made on the number of bypass diodes required per given number of series cells per module or source circuit. In addition, a new method for determining the cell temperature in the laboratory or in the field is discussed.

Gonzalez, C. C.↗

Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters

The double ionization potential (DIP) equation-ofmotion (EOM) coupled-cluster (CC) method with 4-hole−2- particle (4h-2p) excitations on top of the CC with singles, doubles, and triples calculation, abbreviated as DIP-EOMCCSDT(4h-2p), along with its perturbative DIP-EOMCCSD(T)(a)(4h-2p) approximation, are extended to a relativistic four-component (4c) framework. In addition, we introduce and test a new computationally practical DIP-EOMCC approach, which we call DIPEOMCCSD( T)(ã)(4h-2p), that approximates the treatment of 4h- 2p correlations within the DIP-EOMCCSD(T)(a)(4h-2p) method and reduces the $\mathcal{N}$ 8 scaling characterizing DIP-EOMCCSDT(4h- 2p) and DIP-EOMCCSD(T)(a)(4h-2p) to $\mathcal{N}$ 7 with the system size $\mathcal{N}$. Further improvements in computational efficiency are obtained using the frozen natural spinor (FNS) approximation to reduce the numbers of unoccupied spinors entering the correlated steps of the DIP-EOMCC calculations according to a well-defined occupation-number-based threshold. The resulting 4c-FNS-DIPEOMCC approaches are used to compute DIPs for the series of inert gas atoms from argon to radon as well as the vertical DIPs in Cl 2 , Br 2 , HBr, and HI, which have been experimentally examined in the past. We demonstrate that, when using complete basis set extrapolations and FNS truncation threshold of 10 −4.5 , the 4c-FNS-DIP-EOMCCSD(T)(ã)(4h-2p) calculations are capable of predicting DIPs in agreement with experimental data, improving upon their nonrelativistic and spin-free scalar-relativistic counterparts, particularly when examining DIPs characterized by stronger spin−orbit coupling effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Utilizing Physics-Informed Synthetic Data to Train a Digital Twin for Predicting Reactor Operations

Understanding techniques to strengthen the nuclear nonproliferation regime is crucial in reducing the creation of nuclear weaponry on the basis of advancements in the nuclear energy industry. Prior to construction of a nuclear power plant, it is necessary to understand the proliferation potential of the plant’s reactor. Digital twins serve as a unique solution to recognizing reactor behavior indicative of nuclear proliferation. The following research conducted serves as a validation to training a digital twin on synthetic data fabricated via means of reactor physics simulations based on parameters of Idaho State University’s AGN-201 reactor. The synthetic data is utilized to train a long short-term memory (LSTM) recurrent neural network model. The accuracy of the predicted data is measured against real operational data to verify the reliability of the synthetic data creation methods and if these methods should be used in the future to inform inspectors of a reactor’s proliferation capabilities.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Correlation between nonequilibrium thermodynamic theory and V-I curves of a thermionic converter

The methods of irreversible thermodynamics are used as a basis for developing a transport theory applicable to thermionic cesium diodes. This theory is applied in describing certain aspects of experimental voltage-current (V-I) characteristic curves of a thermionic diode. Emphasis is placed on the open-circuit voltage and short-circuit current of the V-I curves. Experimental results show that the open-circuit voltage is essentially independent of the interelectrode spacing, whereas the short-circuit current varies inversely with the interelectrode separation. The inverse dependence is in agreement with the results of the present analysis and also with results expected from diffusion theory, which has been used in other analyses. In fact, the theory used in this investigation reduces to diffusion theory under special assumptions. Deviations from diffusion theory, which were found in experimental results, can be accounted for using the theory developed in this paper.

Robinson, L. B.↗

Investigation of articulated panel dynamics

An articulated panel dynamics experiment to evaluate present analysis and ground test methods as well as damping effectiveness is proposed. The experiment uses an existing panel design which was extensively analyzed and tested. These data provide a firm basis for evaluating the adequacy of panel analysis and ground test methodology. The key issues for future large space structure panel designs are addressed: the critical launch transient and vibroacoustic loading; the deployment analysis adequacy including air and gravitational effects; and the orbital resonant frequencies and mode shapes of deployed panel assemblies. By using an existing mature design that was thoroughly tested, the effort can focus on correlation of actual flight results with existing predictions. A second panel assembly incorporating passive damping is proposed to provide a direct measure of damped panel benefits that can be obtained. These benefits include: reduced launch loads and responses; and highly damped deployed modes. Existing space-qualified viscoelastic epoxy, in combination with composite materials or an alternate more effective damping material, will be used.

Stahle, C. V., Jr.↗

Hot-spot qualification testing of concentrator modules

Results of a study to determine the hot-spot susceptibility of concentrator cells, to provide a hot-spot qualification test for concentrator modules, and to provide guidelines for reducing hot-spot susceptibility are presented. Hot-spot heating occurs in a photovoltaic module when the short-circuit current of a cell is lower than the string operating current, forcing the cell into reverse bias with a concurrent power dissipation. Although the basis for the concentrator-module hot-spot qualification test is the test developed for flat-plate modules, issues such as providing cell illumination introduce additional complexities into the testing procedure. The results indicate that the same general guidelines apply to protecting concentrator modules from hot-spot stressing as apply to flat-plate modules, and recommendations are made on the number of bypass diodes required per given number of series cells per module or source circuit. A method for determining the cell temperature in the laboratory or in the field is discussed.

Gonzalez, C. C.↗

Quantitatively controlled electrophoretic deposition of nanocrystal films from non-aqueous suspensions

This study presents a novel method to correlate the mass and charge transfer kinetics during the electrophoretic deposition of nanocrystal films by using a purpose-built double quartz crystal microbalance combined with simultaneous current-measurement. Our data support a multistep process for film formation: generation of charged nanocrystal flux, charge transfer at the electrode, and polarization of neutral nanocrystals near the electrode surface. The polarized particles are then subject to dielectrophoretic forces that reduce diffusion away from the interface, generating a sufficiently high neutral particle concentration at the interface to form a film. The correlation of mass and charge transfer enables quantification of the nanocrystal charge, the fraction of charged nanocrystals, and the initial sticking coefficient of the particles. These quantities permit calculation of the film thickness, providing a theoretical basis for using concentration and voltage as process parameters to grow films of targeted thicknesses.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Weak entanglement approximation for nuclear structure

The interacting shell model, a configuration-interaction method, is a venerable approach for low-lying nuclear structure calculations, but it is hampered by the exponential growth of its basis dimension as one increases the single-particle space and/or the number of active particles. Recent, quantum-information-inspired work has demonstrated that the proton and neutron sectors of a nuclear wave function are weakly entangled. Furthermore, the entanglement is smaller for nuclides away from N = Z, such as heavy, neutron-rich nuclides. Here, in this study, we implement a weak entanglement approximation to bipartite configuration-interaction wave functions, approximating low-lying levels by coupling a relatively small number of many-proton and many-neutron states. This truncation scheme, which we present in the context of past approaches, reduces the basis dimension by many orders of magnitude while preserving essential features of nuclear spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Determination of in-flight AVIRIS spectral, radiometric, spatial and signal-to-noise characteristics using atmospheric and surface measurements from the vicinity of the rare-earth-bearing carbonatite at Mountain Pass, California

An assessment of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) performance was made for a flight over Mountain Pass, California, July 30, 1987. The flight data were reduced to reflectance using an empirical algorithm which compensates for solar, atmospheric and instrument factors. AVIRIS data in conjunction with surface and atmospheric measurements acquired concurrently were used to develop an improved spectral calibration. An accurate in-flight radiometric calibration was also performed using the LOWTRAN 7 radiative transfer code together with measured surface reflectance and atmospheric optical depths. A direct comparison with coincident Thematic Mapper imagery of Mountain Pass was used to demonstrate the high spatial resolution and good geometric performance of AVIRIS. The in-flight instrument noise was independently determined with two methods which showed good agreement. A signal-to-noise ratio was calculated using data from a uniform playa. This ratio was scaled to the AVIRIS reference radiance model, which provided a basis for comparison with laboratory and other in-flight signal-to-noise determinations.

Green, Robert O.↗

Modeling the Proton Radiation Belt With Van Allen Probes Relativistic Electron-Proton Telescope Data

An empirical model of the proton radiation belt is constructed from data taken during 2013-2017 by the Relativistic Electron-Proton Telescopes on the Van Allen Probes satellites. The model intensity is a function of time, kinetic energy in the range 18-600 megaelectronvolts, equatorial pitch angle, and L shell of proton guiding centers. Data are selected, on the basis of energy deposits in each of the nine silicon detectors, to reduce background caused by hard proton energy spectra at low L. Instrument response functions are computed by Monte Carlo integration, using simulated proton paths through a simplified structural model, to account for energy loss in shielding material for protons outside the nominal field of view. Overlap of energy channels, their wide angular response, and changing satellite orientation require the model dependencies on all three independent variables be determined simultaneously. This is done by least squares minimization with a customized steepest descent algorithm. Model uncertainty accounts for statistical data error and systematic error in the simulated instrument response. A proton energy spectrum is also computed from data taken during the 8 January 2014 solar event, to illustrate methods for the simpler case of an isotropic and homogeneous model distribution. Radiation belt and solar proton results are compared to intensities computed with a simplified, on-axis response that can provide a good approximation under limited circumstances.

Selesnick, R. S.↗