Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iteration 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 991 records · Page 55

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)↗

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)↗

A comparison of several computational techniques for solving some common aeronomic problems

Several numerical integration techniques for solving common aeronomic problems involving species rate equations are compared for speed and accuracy. A newer technique that defines families of species that are nearly conserved is found to be superior to an iterative technique when both methods are applied to simple test problems. The 'conservation' technique is also found to be more economical than the more complex Gear (1969) integration scheme for comparable accuracy.

Turco, R. P.↗

Program for analysis of nonlinear equilibrium and stability (PANES)

PANES utilizes improved techniques for analysis of structures with material and geometric nonlinearities, including limit point and bifurcations behavior which occurs in buckling and collapse problems. Incremental loading, Newton-Raphson iteration, and higher order methods are used in program.

Vos, R. G.↗

Sequential filtering applied to the determination of tracking station locations

The extended sequential filter has been applied to the problem of dynamically determining the geocentric coordinates of two laser satellite tracking stations. This filter provides significant advantages over the classical batch methods through (1) fewer iterations required for convergence, (2) wider radius of convergence, and (3) availability of the parameter estimate evolution. Processing the data sequentially readily identifies the data arcs required to minimize the effects of geopotential model errors. By means of the Smithsonian standard earth 2 and the Goddard earth model 1 geopotentials to reduce laser range observations of the Beacon Explorer-C satellite, it is demonstrated that a two-pass arc is optimal for estimating the height of one station and all coordinates of the second station while minimizing the effect of geopotential model error. These two-pass estimates are in good agreement with other determinations that utilize considerably more data as well as different satellites.

Schutz, B. E.↗

Ocean tides from Seasat-A

Two procedures for the evaluation of global tides from SEASAT-A altimetry data are elaborated: an empirical method leading to the response functions for a grid of about 500 points from which the tide can be predicted for any point in the oceans, and a dynamic method which consists of iteratively modifying the parameters in a numerical solution to Laplace tide equations. It is assumed that the shape of the received altimeter signal can be interpreted for sea state and that orbit calculations are available so that absolute sea levels can be obtained.

Hendershott, M. C.↗

Correlation between multispectral photography and near-surface turbidities

Four-band multispectral photography obtained from an aerial platform at an altitude of about 10,000 feet has been utilized to measure near-surface turbidity at numerous sampling sites in the Ross Barnett Reservoir, Mississippi. Correlation of the photographs with turbidity measurements has been accomplished via an empirical mathematical model which depends upon visual color recognition when the composited photographs are examined on either an I squared S model 600 or a Spectral Data model 65 color-additive viewer. The mathematical model was developed utilizing least-squares, iterative, and standard statistical methods and includes a time-dependent term related to sun angle. This model is consistent with information obtained from two overflights of the target area - July 30, 1973 and October 30, 1973 - and now is being evaluated with regard to information obtained from a third overflight on November 8, 1974.

Wertz, D. L.↗

CLASSY: An adaptive maximum likelihood clustering algorithm

The CLASSY clustering method alternates maximum likelihood iterative techniques for estimating the parameters of a mixture distribution with an adaptive procedure for splitting, combining, and eliminating the resultant components of the mixture. The adaptive procedure is based on maximizing the fit of a mixture of multivariate normal distributions to the observed data using its first through fourth central moments. It generates estimates of the number of multivariate normal components in the mixture as well as the proportion, mean vector, and covariance matrix for each component. The basic mathematical model for CLASSY and the actual operation of the algorithm as currently implemented are described. Results of applying CLASSY to real and simulated LANDSAT data are presented and compared with those generated by the iterative self-organizing clustering system algorithm on the same data sets.

Lennington, R. K.↗

A computer simulation of the midlatitude plasmasphere and ionosphere

A computer model has been developed to simulate species density, temperature, and plasma flow in the ionosphere and plasmasphere. A new approach, the flux preserving method involving a Newton iteration, is used to solve a system of governing equations comprising four second-order partial differentials. Tests are performed to demonstrate that the simulation converges to a stable steady state solution, then steady state ion fluxes are analyzed. Finally, simulations of the collapse of the sunset ionosphere are presented. Comparisons with satellite and radar data show good agreement in a number of cases.

Young, E. R.↗

A finite element model with nonviscous damping

A constitutive law by which structural damping is modeled as a relationship between stress, strain, and strain rate in a material is used in conjunction with the finite element method to develop general integral expressions for viscous and nonviscous damping matrices. To solve the set of nonlinear equations resulting from the presence of nonviscous damping, a solution technique is developed by modifying the Newmark method to accommodate an iterative solution and treat the nonviscous damping as a pseudo-force. The technique is then checked for accuracy and convergence in single- and multi-degree-of-freedom problems, and is found to be accurate and efficient for initial-condition problems with small nonviscous damping.

Roussos, L. A.↗