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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hydrodynamic and Radiographic Toolbox (HART)

With a multi-lab and university team, we propose to develop new methods for a Hydrodynamic and Radiographic Toolbox (HART) that will enable a fuller and more extensive use of experimental radiographic data towards better characterizing and reducing uncertainties in predictive modeling of weapons performance. We will achieve this by leveraging recent developments in areas of computational imaging, statistical and machine learning, and reduced order modeling of hydrodynamics. In terms of software practices, by partnering with XCP (RISTRA Project) we will conform to recent XCP standards that are inline with modern software practices and standards and ensure compatibility and ease of inter-operability with existing codes. Three new activities under this proposal include (a) the development and use of deep learning-based surrogates to accelerate reconstruction and variational inference of density fields from radiographs of hydrotests, (b) a model-data fusion strategy that couples deep learning-based density reconstructions with fast hydrodynamics simulators to better constrain the reconstruction, and (c) a method for treating asymmetries using techniques adopted from limited view tomography. All three new activities will be based on improved treatment of scatter, noise, beam spot movement, detector blur, and flat fielding in the forward model, and a use of sophisticated priors to aid in the re construction. The improvements to the forward model and improved algorithmic design of the reconstruction when complete will be contained in the iterative reconstruction code SHIVA—a code project that we have recently initiated. The many ways in which machine learning can be used in the reconstruction work will be contained in a code HERMES that has been initiated with DTRA support. For example, the significant levels of acceleration that will likely be achieved by the use of machine learning techniques will permit us (and are required) to quantify uncertainties in density retrievals. Next, the two-way coupling between density reconstruction and model-based simulation of the hydrodynamics will be contained in code EREBUS, and will permit a fuller realization of the potential of the data to constrain the hydrodynamic model and better address issues related to asymmetries in the problem. Finally, we anticipate that the better consistency with physics achieved in our reconstructions will allow them to be used by X-Division more so than today.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Feasibility Assessment of Using Electrical Impedance Tomography for Damage Localization in Graphite Microreactor Components

Microreactors have great potential to decrease capital costs and construction timelines, reducing the barriers to implementing advanced nuclear reactor technologies. However, the lower power output of these microreactors introduces economic challenges if the operation and maintenance costs cannot be reduced sufficiently. Many microreactor concepts use graphite materials for in-core moderator and structural purposes, which will require periodic inspection or, ideally, in situ structural health monitoring. This work describes an initial evaluation of the feasibility of leveraging graphite’s semiconducting properties to perform electrical impedance tomography (EIT) for defect localization. First, a study was performed to identify the best methods for bonding electrical contacts to minimize the effects of contact resistance that interfere with impedance measurements of the graphite. After determining the best bonding approach, electrodes were attached to multiple graphite components with varying geometries (e.g., simple block and more representative microreactor hexagonal block). In parallel, finite element analysis approaches were developed and implemented to analyze the impact of defects and contact resistances and to inform an iterative inverse model for reconstructing the conductivity distribution. The results of this study show that defect localization in graphite components is possible using EIT if a sufficient number of electrodes (to improve spatial resolution) can be bonded using a technique with very low contact resistance and those contact resistances remain stable during reactor operation. If the reactor core and vessel design do not allow electrodes to be bonded during reactor operation, then it may be possible to detach/reattach electrodes between reactor operation and outages if a suitable mechanical connection that offers repeatable contact resistance can be identified. Alternatively, future work could focus on identifying contacts that can survive the harsh reactor operating environment so that contacts do not need to be removed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Sequential linear integer programming for integer optimal control with total variation regularization

We propose a trust-region method that solves a sequence of linear integer programs to tackle integer optimal control problems regularized with a total variation penalty. The total variation penalty implies that the considered integer control problems admit minimizers. We introduce a local optimality concept for the problem, which arises from the infinite-dimensional perspective. In the case of a one-dimensional domain of the control function, we prove convergence of the iterates produced by our algorithm to points that satisfy first-order stationarity conditions for local optimality. We demonstrate the theoretical findings on a computational example.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Reconstructing Quasar Spectra and Measuring the Lyα Forest with SpenderQ

Quasar spectra carry the imprint of foreground intergalactic medium (IGM) through absorption features. In particular, absorption caused by neutral hydrogen gas, the "Lyα forest," is a key spectroscopic tracer for cosmological analyses used to measure cosmic expansion and test physics beyond the standard model. Despite their importance, current methods for measuring Lyα absorption cannot directly derive the intrinsic quasar continuum and make strong assumptions on its shape, thus distorting the measured Lyα clustering. We present SpenderQ , a ML-based approach for directly reconstructing the intrinsic quasar spectra and measuring the Lyα forest from observations. SpenderQ uses the Spender spectrum autoencoder to learn a compact and redshift-invariant latent encoding of quasar spectra, combined with an iterative procedure to identify and mask absorption regions. To demonstrate its performance, we apply SpenderQ to 400,000 synthetic quasar spectra created to validate the Dark Energy Spectroscopic Instrument Year 1 Lyα cosmological analyses. SpenderQ accurately reconstructs the true intrinsic quasar spectra, including the broad Lyβ, Lyα, SiIV, CIV, and CIII emission lines. Redward of Lyα, SpenderQ provides percent-level reconstructions of the true quasar spectra. Blueward of Lyα, SpenderQ reconstructs the true spectra to < 5%. SpenderQ reproduces the shapes of individual quasar spectra more robustly than the current state-of-the-art. We, thus, expect it will significantly reduce biases in Lyα clustering measurements and enable studies of quasars and their physical properties. SpenderQ also provides informative latent variable encodings that can be used to, e.g., classify quasars with Broad Absorption Lines. Overall, SpenderQ provides a new data-driven approach for unbiased Lyα forest measurements in cosmological, quasar, and IGM studies.

Hahn, ChangHoon [Arizona U., Astron. Dept. - Stewa↗
DOE OSTI2021

Towards a reliable assessment of charging effects during surface analysis: Accurate spectral shapes of ZrO 2 and Pd/ZrO 2 via X-ray Photoelectron Spectroscopy

X-ray Photoelectron Spectroscopy of large bandgap or insulating material surfaces relies on an effective mechanism that compensates for the emission (loss) of electrons by maintaining the material surface at a steady-state uniform potential. While a steady-state may be attained by utilizing an active compensation, such as low power electron emitting filament, there is the possibility that the surface potential is not uniform over the area analyzed, leading to peak shifts and incorrect spectral interpretation. Here, in this work, a spectral data processing method based on mapping the ZrO 2 and Pd/ZrO 2 surfaces utilizing photoemission peak binding energy is proposed, which provides information about the response of specific material surfaces to charge compensation. Spectromicroscopy of ZrO 2 and Pd/ZrO 2 surfaces without spatial information is used to monitor the efficacy of charge compensation. Exploiting counts distributed over many bins require the use of procedures and algorithms essential to practical mapping peak positions. Iterative singular value decomposition is therefore introduced and utilized as a means of efficiently delivering spatially resolved spectra from which binding energy for peaks is computed. The concepts developed in this work result in robust and accurate peak models of ZrO 2 and Pd/ZrO 2 that can be applied in XPS analysis of not only ZrO 2 but other large bandgap or insulating material surfaces. Supporting arguments for a peak model representing signal from Zr 3p and Pd 3d are developed within this work are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2020

N 5 -Scaling Excited-State-Specific Perturbation Theory

In this paper, we show that by working in a basis similar to that of the natural transition orbitals and using a modified zeroth-order Hamiltonian, the cost of a recently introduced perturbative correction to excited-state mean field theory can be reduced from seventh to fifth order in the system size. The (occupied) 2 (virtual) 3 asymptotic scaling matches that of ground-state second-order Møller-Plesset theory but with a significantly higher prefactor because the bottleneck is iterative: it appears in the Krylov-subspace-based solution of the linear equation that yields the first-order wave function. Here, we discuss the details of the modified zeroth-order Hamiltonian we use to reduce the cost and the automatic code generation process we used to derive and verify the cost scaling of the different terms. Overall, we find that our modifications have little impact on the method's accuracy, which remains competitive with singles and doubles equation-of-motion coupled cluster.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Electromagnetic Analysis of ITER Electron Cyclotron Emission Components and Support Structures for the Preliminary Design Review

The ITER electron cyclotron emission (ECE) diagnostic system is located at Design Shield Module (DSM) 2, Equatorial Port (EP) 9 to measure electron temperature profile and electron temperature fluctuations and also to assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two couples of mirrors for two different optical views. This ECE diagnostic system shall be exposed to significant power due to unabsorbed electron cyclotron heating power in the plasma. It shall also receive large electromagnetic (EM) loads up to 100 MN/m 3 force density due to the eddy currents generated by the short 16-ms transient plasma disruption. The global EM models of EP12 and EP11 with the worst case of plasma disruption Major Disruption Downward Exponential 16 ms Category III (MD_DW_EXP16MS_CATIII) have been provided by ITER. However, the support structure for these critical components, including four mirrors, is still under the design and development phase, especially during the period of preliminary design reviews (PDRs). For those small in-vessel components that are not modeled explicitly for the PDRs, we can extract the magnetic field (B) data and flux variations (dB/dt) at those locations from the global model and then use the B and dB/dt method to calculate the eddy currents, EM forces, and torques on those components in the local submodels. In conclusion, this paper will present the detailed PDR EM analysis results of the ECE components and support structures.

Electromagnetic↗
DOE OSTI2024

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure of the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

Here we present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. The new model enables fully 3D simulations of complete three-component waveforms and provides improved fits that were not possible with previous models. We inverted over 32,000 waveforms from 192 earthquakes recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region with highly uneven coverage. Inversion iterations proceeded from the period band 50–100 s in six stages and 54 total iterations reducing the minimum period to 30 s. Our final model, MESWA, improves waveform fits compared to the starting and other models for both the data used in the inversion and an independent validation set of 66 events. Restitution tests indicate that the model resolves features in the central part of the model to depths of about 150 km. The new model reveals tectonic features imaged by other studies and methods but in a new holistic model of anisotropic shear and compressional wavespeeds (V S and V P , respectively) covering a larger domain with smaller scale length and amplified features. Examples include low crustal V S in the Tethyan belt and low mantle VS following divergent (Gulf of Aden, Red Sea) and transform (Dead Sea fault) margins of the Arabian plate. Low V S is imaged below Cenozoic volcanic centers of the Mecca–Madina–Nafud Line, Arabian Peninsula, and the Türkiye–Iran border region. Elevated V S tracks Makran subduction under southeast Iran with near vertical dip. MESWA could be used as a starting model for further improvements, say, using waveforms from in-country seismic networks that are not currently openly available and/or smaller-scale studies targeting a shorter period. The model could be used to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A dynamic Bayesian optimized active recommender system for curiosity-driven partially Human-in-the-loop automated experiments

Optimization of experimental materials synthesis and characterization through active learning methods has been growing over the last decade, with examples ranging from measurements of diffraction on combinatorial alloys at synchrotrons, to searches through chemical space with automated synthesis robots for perovskites. In virtually all cases, the target property of interest for optimization is defined a priori with the ability to shift the trajectory of the optimization based on human-identified findings during the experiment is lacking. Thus, to highlight the best of both human operators and AI-driven experiments, here we present the development of a human–AI collaborated experimental workflow, via a Bayesian optimized active recommender system (BOARS), to shape targets on the fly with human real-time feedback. Here, the human guidance overpowers AI at early iteration when prior knowledge (uncertainty) is minimal (higher), while the AI overpowers the human during later iterations to accelerate the process with the human-assessed goal. We showcase examples of this framework applied to pre-acquired piezoresponse force spectroscopy of a ferroelectric thin film, and in real-time on an atomic force microscope, with human assessment to find symmetric hysteresis loops. It is found that such features appear more affected by subsurface defects than the local domain structure. This work shows the utility of human–AI approaches for curiosity driven exploration of systems across experimental domains.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling the co-assembly of binary nanoparticles

Abstract In this work, we present a binary assembly model that can predict the co-assembly structure and spatial frequency spectra of monodispersed nanoparticles with two different particle sizes. The approach relies on an iterative algorithm based on geometric constraints, which can simulate the assembly patterns of particles with two distinct diameters, size distributions, and at various mixture ratios on a planar surface. The two-dimensional spatial-frequency spectra of the modeled assembles can be analyzed using fast Fourier transform analysis to examine their frequency content. The simulated co-assembly structures and spectra are compared with assembled nanoparticles fabricated using transfer coating method are in qualitative agreement with the experimental results. The co-assembly model can also be used to predict the peak spatial frequency and the full-width at half-maximum bandwidth, which can lead to the design of the structure spectra by selection of different monodispersed particles. This work can find applications in fabrication of non-periodic nanostructures for functional surfaces, light extraction structures, and broadband nanophotonics.

Mohanty, Saurav (ORCID:0009000302592855)↗
DOE OSTI2024

FEM Analysis of Hybrid LTS/HTS Cos-Theta Dipole Magnet

After the recent results on Bi-2212 superconductive magnets realized and tested in canted cosine-theta and solenoid designs, respectively, at Lawrence Berkeley National Laboratory (LBNL) and NHMFL, the first Bi-2212 stress-managed cosine-theta insert magnet is in the assembly phase at Fermilab. This insert will be part of the first hybrid cosine-theta magnet made of Nb3Sn outer layers within the US-MDP effort to reach 20 T bore field. This paper reports the analytical analysis of the cosine-theta Nb3Sn/Bi-2212 hybrid magnet. As an introduction, it shows the validation of the finalized HTS coils' design optimized using the ROXIE code. Subsequently, it reports the parameters, logic, and implementation method of the 2D electromagnetic and mechanical FEM analysis of the LTS/HTS hybrid magnet. Results from the ANSYS detailed model made with sub-modeling geometry are compared with the homogeneous model implemented in the past. Moreover, a Python code was implemented to simulate the current degradation due to stresses in the detail-modeled conductor areas. The current degradation is introduced in the simulation dynamics for both conductors as an iteration process during the energization load step only after applying pre-stress and thermal loads. The numerical and graphic results of the 2D cosine-theta LTS/HTS magnet cross-section will be described and analyzed.

D'Agliano, A.↗
DOE OSTI2020

Development of the finite slab geometry in the NORA code

In this work we develop the complete linearization formalism for the radiation transfer equation in slab geometry. This development was motivated by an observed convergence instability arising from the NORA code’s inner boundary while modeling a suite of simple test problems in slab geometry presented at the Line Transfer Workshop held in LLNL in May of 2019. The convergence issue was seen to emerge early on in Newton-Raphson iterations. These test problems required parameters in the former NORA inner boundary to achieve the asymptotic limit of zero incoming intensity. In this work we derive the limiting case directly to gain insight into the former algorithm’s pathology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A note on higher-order and nonlinear limiting approaches for continuously bounds-preserving discontinuous Galerkin methods

In Dzanic (2024), a limiting approach for high-order discontinuous Galerkin schemes was introduced which allowed for imposing constraints on the solution continuously (i.e., everywhere within the element). While exact for linear constraint functionals, this approach only imposed a sufficient (but not the minimum necessary) amount of limiting for nonlinear constraint functionals. This short note shows how this limiting approach can be extended to allow exactness for general nonlinear quasiconcave constraint functionals through a nonlinear limiting procedure, reducing unnecessary numerical dissipation. Finally, some examples are shown for nonlinear pressure and entropy constraints in the compressible gas dynamics equations, where both analytic and iterative approaches are used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

MAPPRAISER: A massively parallel map-making framework for multi-kilo pixel CMB experiments

Forthcoming cosmic microwave background (CMB) polarized anisotropy experiments have the potential to revolutionize our understanding of the Universe and fundamental physics. The sought-after, tale-telling signatures will be however distributed over voluminous data sets which these experiments will collect. These data sets will need to be efficiently processed and unwanted contributions due to astrophysical, environmental, and instrumental effects characterized and efficiently mitigated in order to uncover the signatures. This poses a significant challenge to data analysis methods, techniques, and software tools which will not only have to be able to cope with huge volumes of data but to do so with unprecedented precision driven by the demanding science goals posed for the new experiments. A keystone of efficient CMB data analysis is solvers of very large linear systems of equations. Such systems appear in very diverse contexts throughout CMB data analysis pipelines, however they typically display similar algebraic structures and can therefore be solved using similar numerical techniques. Linear systems arising in the so-called map-making problem are one of the most prominent and common ones. In this work we present a massively parallel, flexible and extensible framework, comprised of a numerical library, MIDAPACK, and a high level code, MAPPRAISER, which provide tools for solving efficiently such systems. Here, the framework implements iterative solvers based on conjugate gradient techniques: enlarged and preconditioned using different preconditioners. We demonstrate the framework on simulated examples reflecting basic characteristics of the forthcoming data sets issued by ground-based and satellite-borne instruments, executing it on as many as 16,384 compute cores. The software is developed as an open source project freely available to the community at: https://github.com/B3Dcmb/midapack.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Physics basis for the Wisconsin HTS Axisymmetric Mirror (WHAM)

The Wisconsin high-temperature superconductor axisymmetric mirror experiment (WHAM) will be a high-field platform for prototyping technologies, validating interchange stabilization techniques and benchmarking numerical code performance, enabling the next step up to reactor parameters. A detailed overview of the experimental apparatus and its various subsystems is presented. WHAM will use electron cyclotron heating to ionize and build a dense target plasma for neutral beam injection of fast ions, stabilized by edge-biased sheared flow. At 25 keV injection energies, charge exchange dominates over impact ionization and limits the effectiveness of neutral beam injection fuelling. This paper outlines an iterative technique for self-consistently predicting the neutral beam driven anisotropic ion distribution and its role in the finite beta equilibrium. Beginning with recent work by Egedal et al. ( Nucl. Fusion , vol. 62, no. 12, 2022, p. 126053) on the WHAM geometry, we detail how the FIDASIM code is used to model the charge exchange sources and sinks in the distribution function, and both are combined with an anisotropic magnetohydrodynamic equilibrium solver method to self-consistently reach an equilibrium. We compare this with recent results using the CQL3D code adapted for the mirror geometry, which includes the high-harmonic fast wave heating of fast ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accelerating the design of lattice structures using machine learning

Lattices remain an attractive class of structures due to their design versatility; however, rapidly designing lattice structures with tailored or optimal mechanical properties remains a significant challenge. With each added design variable, the design space quickly becomes intractable. To address this challenge, research efforts have sought to combine computational approaches with machine learning (ML)-based approaches to reduce the computational cost of the design process and accelerate mechanical design. While these efforts have made substantial progress, significant challenges remain in (1) building and interpreting the ML-based surrogate models and (2) iteratively and efficiently curating training datasets for optimization tasks. Here, we address the first challenge by combining ML-based surrogate modeling and Shapley additive explanation (SHAP) analysis to interpret the impact of each design variable. We find that our ML-based surrogate models achieve excellent prediction capabilities (R 2 > 0.95) and SHAP values aid in uncovering design variables influencing performance. We address the second challenge by utilizing active learning-based methods, such as Bayesian optimization, to explore the design space and report a 5 × reduction in simulations relative to grid-based search. Collectively, these results underscore the value of building intelligent design systems that leverage ML-based methods for uncovering key design variables and accelerating design.

36 MATERIALS SCIENCE↗