Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

Interpretation of multiple solutions in fully iterative GF2 and GW schemes using local analysis of two-particle density matrices

Due to the presence of non-linear equations, iterative Green’s function methods can result in multiple different solutions even for simple molecular systems. In contrast to the wave-function methods, a detailed and careful analysis of such molecular solutions was not performed before. In this work, we use two-particle density matrices to investigate local spin and charge correlators that quantify the charge resonance and covalent characters of these solutions. When applied within the unrestricted orbital set, spin correlators elucidate the broken symmetry of the solutions, containing necessary information for building effective magnetic Hamiltonians. Based on GW and GF2 calculations of simple molecules and transition metal complexes, we construct Heisenberg Hamiltonians, four-spin-four-center corrections, and biquadratic spin–spin interactions. These Hamiltonian parameterizations are compared to previous wave-function calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geometry optimizations with spinor-based relativistic coupled-cluster theory

Development of analytic gradients for relativistic coupled-cluster singles and doubles augmented with a non-iterative triples [CCSD(T)] method using an all-electron exact two-component Hamiltonian with atomic mean-field spin–orbit integrals (X2CAMF) is reported. This enables efficient CC geometry optimizations with spin–orbit coupling included in orbitals. Here, the applicability of the implementation is demonstrated using benchmark X2CAMF-CCSD(T) calculations of equilibrium structures and harmonic vibrational frequencies for methyl halides, CH 3 X (X = Br, I, and At), as well as calculations of rotational constants and infrared spectrum for RaSH+, a radioactive molecular ion of interest to spectroscopic study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural mode coupling in perovskite oxides using hypothesis-driven active learning

Abstract Finding the ground-state structure with minimum energy is paramount to designing any material. In ABO 3 -type perovskite oxides with Pnma symmetry, the lowest energy phase is driven by an inherent trilinear coupling between the two primary order parameters such as rotation and tilt with antiferroelectric displacement of the A-site cations as established via hybrid improper ferroelectric mechanism. Conventionally, finding the relevant mode coupling driving phase transition requires performing first-principles calculations which is computationally time-consuming as well as expensive. It involves following an intuitive iterative hit and trial method of (a) adding two or multiple mode vectors, followed by (b) evaluating which combination would lead to the ground-state energy. In this study, we show how a hypothesis-driven active learning framework can identify suitable mode couplings within the Landau free energy expansion with minimal information on amplitudes of modes for a series of double perovskite oxides with A-site layered, columnar and rocksalt ordering. This scheme is expected to be applicable universally for understanding atomistic mechanisms derived from various structural mode couplings behind functionalities, for e.g. polarization, magnetization and metal–insulator transitions.

36 MATERIALS SCIENCE↗

Phase retrieval and reconstruction of coherent synthesis by genetic algorithm

Abstract In the context of diffractive optics, phase retrieval is a heavily investigated process of recreating an entire complex electric field from partial amplitude-only information through iterative algorithms. However, existing methods can fall into local minima during reconstructions or struggle to recover unusual and novel electric field distributions. We present a numerical method based on a global-optimization genetic algorithm that reconstructs non-trivial electric field distributions from single diffracted intensity distributions. Diffraction and propagation of the optical fields over arbitrary distances is modeled through implementation of the angular spectrum technique. Additionally, a coherently-locked laser array system is used as an experimental case-study demonstrating 0.09 π phase reconstruction accuracy of initial laser parameters from single intensity images.

59 BASIC BIOLOGICAL SCIENCES↗

An Edge Alignment-Based Orientation Selection Method for Neutron Tomography

Neutron computed tomography (nCT) is a 3D char-acterization technique used to image the internal morphology or chemical composition of samples in biology and materials sciences. A typical workflow involves placing the sample in the path of a neutron beam, acquiring projection data at a predefined set of orientations, and processing the resulting data using an analytic reconstruction algorithm. Typical nCT scans require hours to days to complete and are then processed using conventional filtered back-projection (FBP), which performs poorly with sparse views or noisy data. Hence, the main methods in order to reduce overall acquisition time are the use of an improved sampling strategy combined with the use of advanced reconstruction methods such as model-based iterative reconstruction (MBIR). In this paper, we propose an adaptive orientation selection method in which an MBIR reconstruction on previously-acquired measurements is used to define an objective function on orientations that balances a data-fitting term promoting edge alignment and a regularization term promoting orientation diversity. Using simulated and experimental data, we demonstrate that our method produces high-quality reconstructions using significantly fewer total measurements than the conventional approach.

Yang, Diyu↗

Transactive Coordination of Flexible Loads with Energy Storage through Day-ahead Scheduling

We are motivated by the problem faced by the aggregator (market coordinator) to optimally schedule the day-ahead operation of flexible loads and energy storages through transactive coordination. The aggregator, which may be the ISO/DSO/LSE/Utility company/Microgrid controller depending on the system considered, aims to maximize social welfare of interconnected distributed energy resources (DERs) by determining the optimal hourly prices for the next operating day. Arbitrage and private objectives are incorporated to ensure participation of third-party energy storages. We show that the coordination problem can be mathematically formulated as a multi-period optimal control problem with global power flow constraints and local DER-level constraints, and propose an iterative hierarchical decentralized solution method which respects the privacy of DERs. Under proper convexity assumptions, the algorithm will converge to the optimal day-ahead market clearing prices. We illustrate the algorithm using a four-bus power system with DC power flow constraints. Results show that incorporating degradation cost decreases the energy storage's ability of peak shaving and that social welfare of the system increases when discharging and charging efficiencies increase.

Ma, Ke↗

Endogenizing Probabilistic Resource Adequacy Risks in Deterministic Capacity Expansion Models

In this work, we demonstrate how power system capacity expansion models can understate the stochastic effects of thermal outages when considering resource availabilities on an hourly expected value basis, yielding system designs with multiple orders of magnitude more shortfall risk than stated adequacy targets. We develop a novel approximation approach to efficiently endogenize awareness of this risk in a deterministic, linear capacity expansion framework. We compare this approach to exogenous tuning of an energy reserve margin, the leading alternative method to compensate for unmodeled probabilistic shortfall risk. Empirical results from a test system show that the new endogenous method cost-effectively meets all regional reliability targets with a single optimization solve, and produces a near-identical system design as the incumbent method without the need for repeated re-optimizations to find an appropriate reserve level. The endogenous method may also use iterative re-optimizations to further improve solution quality, although these incremental benefits were modest in the system studied.

capacity expansion modeling↗

Control optimization for parametric Hamiltonians by pulse reconstruction

Optimal control techniques provide a means to tailor the control pulses required to generate customized quantum gates, which helps to improve the resilience of quantum simulations to gate errors and device noise. However, the significant amount of (classical) computation required to generate customized gates can quickly undermine the effectiveness of this approach, especially when pulse optimization needs to be iterated. We propose a method to reduce the computational time required to generate the control pulse for a Hamiltonian that is parametrically dependent on a time-varying quantity. We use simple interpolation schemes to accurately reconstruct the control pulses from a set of pulses obtained in advance for a discrete set of predetermined parameter values. We obtain a reconstruction with very high fidelity and a significant reduction in computational effort. We report the results of the application of the proposed method to device-level quantum simulations of the unitary (real) time evolution of two interacting neutrons based on superconducting qubits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Parametric study of optical transmission through plasmonic hole arrays modulated by the phase transition of vanadium dioxide

We have performed comprehensive electromagnetic simulations and preliminary experiments to explore the effects of geometrical and material parameters on the extraordinary optical transmission (EOT) through periodic arrays of subwavelength holes in a bilayer stack consisting of a gold or silver film atop a vanadium dioxide film (Au/Ag + VO 2 ), where the latter undergoes a semiconductor-to-metal phase transition. Using the finite-difference time-domain (FDTD) and finite-element methods (FEM), we vary iteratively the array periodicity, VO 2 film thickness and hole diameters, as well as the refractive index inside the VO 2 -layer holes and the VO 2 optical constants. For each variation, we compare the metallic-to-semiconducting ratios of the zero-order transmission ( T 00 ) peaks and find sharp maxima in these ratios within narrow parameter ranges. The maxima arise from Fabry-Perot and Fano-type resonances that minimize T 00 in the semiconducting phase of the perforated bilayers. At a fixed array period, the primary factors controlling the VO 2 -enabled EOT modulation are the VO 2 thickness, diameter of the VO 2 -layer holes, and absorption in the two VO 2 phases. Besides uncovering the origins of the higher metallic-phase T 00 , this study provides a protocol for optimizing the performance of the bilayer hole arrays for potential uses as dynamically tunable nano-optical devices.

Donev, Eugenii U. (ORCID:000000030555219X)↗

12th IAEA Technical Meeting on Energetic Particles in Magnetic Confinement Systems (7–11 September 2011) and 5th IAEA Technical Meeting on Theory of Plasma Instabilities (5–7 September 2011)

The 12th IAEA Technical Meeting on Energetic Particles in Magnetic Confinement Systems took place in Austin, Texas (7–11 September 2011). This meeting was organized jointly with the 5th IAEA Technical Meeting on Theory of Plasma Instabilities (5–7 September 2011). The two meetings shared one day (7 September 2011) with presentations relevant to both groups. Some of the work reported at these meetings was then published in a special issue of Nuclear Fusion [Nucl. Fusion 52 (2012)]. Summaries of the Energetic Particle Conference presentations were given by Kazuo Toi and Boris Breizman. They respectively discussed the experimental and theoretical progress presented at the meeting. Highlights of this meeting include the tremendous progress that has been achieved in the development of diagnostics that enables the ‘viewing’ of internal fluctuations and allows comparison with theoretical predictions, as demonstrated, for example, in the talks of P. Lauber and M. Osakabe. The need and development of hardened diagnostics in the severe radiation environment, such as those that will exist in ITER, was discussed in the talks of V. Kiptily and V.A. Kazakhov. In theoretical studies, much of the effort is focused on nonlinear phenomena. For example, detailed comparison of theory and experiment on D-III-D on the n = 0 geodesic mode was reported in separate papers by R. Nazikian and G. Fu. A large number of theoretical papers were presented on wave chirping including a paper by B.N. Breizman, which notes that wave chirping from a single frequency may emanate continuously once marginal stability conditions have been established. Another area of wide interest was the detailed study of alpha orbits in a burning plasma, where losses can come from symmetry breaking due to finite coil number or magnetic field imperfections introduced by diagnostic or test modules. An important area of development, covered by M.A. Hole and D.A. Spong, is concerned with the self-consistent treatment of the induced fields that accounts for toroidally asymmetric MHD response. In addition, a significant number of studies focused on understanding nonlinear behavior by means of computer simulation of energetic particle driven instability. An under-represented area of investigation was the study of electron runaway formation during major tokamak disruptions. It was noted in an overview by S. Putvinski that electron energies in the 10–20 MeV range is to be expected during projected major disruptions in ITER and that reliable methods for mitigation of the runaway process needs to be developed. Significant recent work in the field of the disruption induced electron runaway, which was reported by J. Riemann, had been submitted to Physics of Plasmas. Overall it is clear that reliable mitigation of electron runaway is an extremely important topic that is in need of better understanding and solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extending Measurements to En=30 MeV and Beyond

The majority of energy release in the fission process is due to the kinetic energy of the fission fragments. Average Total Kinetic Energy measurements for the major actinides over a wide range of incident neutron energies were performed at LANSCE using a Frisch-gridded ionization chamber. The experiments and results of the 238U(n,f) and 235U(n,f) will be presented, including (En), (A), and mass yield distributions as a function of neutron energy. A preliminary (En) for 239Pu(n,f) will also be shown. The (En) shows a clear structure at multichance fission thresholds for all the reactions that we studied. The fragment masses are determined using the iterative double energy (2E) method, with a resolution of A = 4 - 5 amu. The correction for the prompt fission neutrons is the main source of uncertainty, especially at high incident neutron energies, since the behavior of nubar(A,En) is largely unknown. Different correction methods will be discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Atmospheric Emitted Radiance Interferometer Optimal Estimation (AERIoe) Value-Added Product Report

The Atmospheric Emitted Radiance Interferometer (AERI) Optimal Estimation (AERIoe) algorithm (Turner and Löhnert 2014, Turner and Blumberg 2019) retrieves boundary-layer profiles of temperature and water vapor mixing ratio, together with cloud properties for a single-layer cloud (i.e., LWP, effective radius), from AERI-observed infrared radiance spectra and other measurements. The current AERIoe Value-Added Product (VAP) uses version 2.8 of David Turner's code, which uses additional inputs to constrain the retrieval including surface meteorology, microwave radiometer brightness temperatures at any number of frequencies (the U.S. Department of Energy Atmospheric Radiation Measurement [ARM] user facility typically used observations at 23.8 and 30.0 GHz), and Rapid Refresh numerical weather prediction model output. The method is a physical iterative retrieval that uses the optimal estimation framework (Rodgers 2000), so a full error covariance of the solution and the information content (in the form of degrees of freedom for signal) is provided for each retrieval.

54 ENVIRONMENTAL SCIENCES↗

Impact of Glass Irradiation on Laser-induced Breakdown Spectroscopy Diagnostics in the Visible and NIR Range

Increased absorption of optical materials arising from the exposure to ionizing radiation must be accounted for to accurately analyze the laser-induced breakdown spectroscopy (LIBS) data retrieved from high-radiation environments. We evaluate this effect on a specifc example that mimics the diagnostics placed within a gas-cooled fast reactor coolant stream. Analysis is performed on the LIBS data measured with 1% Xe gas in an ambient He environment and the measured optical absorption from the gamma- and neutron-irradiated low-OH fused silica and sapphire glasses. No significant change in the number of shots required to reach a 3σ detection level was observed up to 10 Mrad gamma dose and exposure to a 1.7×10 17 n/cm 2 neutron fluence. In contrast, the spectral dependence of absorption results in either a systematic overestimation or underestimation of line intensity ratios, depending on the line of interest selected for analysis. Moreover, if lines from different spectral regions are used to create Boltzmann plots, this attenuation also leads to statistically significant changes in the calculated temperatures for calculations using Xe II lines only, lowering them from 8000±610 K to 6800±810 K for the case of exposure to the 1.7×10 17 n/cm 2 flux, and broadening the temperature range required for a 95% confidence interval. In the case of measuring the Xe spectrum, these effects may be mitigated using only the longer-wavelength spectral region, where radiation-attenuation is relatively small, or through analysis using the iterative Saha-Boltzmann method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part II - Clustering Extremely High-Dimensional Grid-Based Data

Building an accurate surrogate model for the spatio-temporal outputs of a computer simulation is a challenging task. A simple approach to improve the accuracy of the surrogate is to cluster the outputs based on similarity and build a separate surrogate model for each cluster. This clustering is relatively straightforward when the output at each time step is of moderate size. However, when the spatial domain is represented by a large number of grid points, numbering in the millions, the clustering of the data becomes more challenging. In this report, we consider output data from simulations of a jet interacting with high explosives. These data are available on spatial domains of different sizes, at grid points that vary in their spatial coordinates, and in a format that distributes the output across multiple files at each time step of the simulation. We first describe how we bring these data into a consistent format prior to clustering. Borrowing the idea of random projections from data mining, we reduce the dimension of our data by a factor of thousand, making it possible to use the iterative k-means method for clustering. We show how we can use the randomness of both the random projections, and the choice of initial centroids in k-means clustering, to determine the number of clusters in our data set. Our approach makes clustering of extremely high dimensional data tractable, generating meaningful cluster assignments for our problem, despite the approximation introduced in the random projections.

97 MATHEMATICS AND COMPUTING↗

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of in situ thickness sensor for vapor diffused Nb3Sn films

Fermilab is one of the leaders in development of vapor diffused Nb3Sn films inside niobium cavities. This material has a higher critical temperature (Tc) than niobium, enabling cavity operation at 4.2 K. This higher operational temperature significantly reduces the infrastructure required for cooling compared to 2 K systems, making superconducting radio-frequency (SRF) technology more accessible. Current deposition methods have relied on iterative testing to determine nominal film thickness, a process that can be time-consuming and imprecise. To address this, we are developing a sensor to measure the thickness of Nb3Sn thin film in situ during vapor diffusion. Our design involves the resistance measurement of a thin film of niobium. During coating, the change in resistance reflects the conversion of the film to Nb3Sn, which allows simple integration with the current furnace infrastructure. This sensor would allow real time measurement of the thickness, allowing for increased precision in future depositions.

Howard, K. [Unlisted] (ORCID:0000000325613308)↗

Development of in situ thickness sensor for vapor diffused Nb3Sn films

Fermilab is one of the leaders in development of vapor diffused Nb3Sn films inside niobium cavities. This material has a higher critical temperature (Tc) than niobium, enabling cavity operation at 4.2 K. This higher operational temperature significantly reduces the infrastructure required for cooling compared to 2 K systems, making superconducting radio-frequency (SRF) technology more accessible. Current deposition methods have relied on iterative testing to determine nominal film thickness, a process that can be time-consuming and imprecise. To address this, we are developing a sensor to measure the thickness of Nb3Sn thin film in situ during vapor diffusion. Our design involves a four-point resistance measurement of a thin film of niobium, inside the coating region. During coating, the change in resistance reflects the conversion of the film from Nb to Nb3Sn, which allows simple integration with the current furnace infrastructure. This sensor would allow real time measurement of the Nb3Sn film thickness, allowing for increased precision in future depositions for cavity applications.

Howard, K. [Chicago U.] (ORCID:0000000325613308)↗

Noninvasive cross section reconstruction with transport theory constraints

We consider the inverse problem of identifying the spatially variable absorption and scattering properties of a medium by measuring the exiting radiation when the body is actively interrogated. We formulate this inverse problem as a PDE-constrained optimization problem and solve it iteratively with Newton's method. The constraint is given by the radiative transport equation for neutral particles. Two examples are considered. The first is a dual inclusion domain with no scattering. This problem explores the convergence patterns of the method. The second problem is a central inclusion problem with scattering. This problem explores the optical thickness limit of the method. This optical thickness was determined to be 2-3 mean free paths. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗