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,315 records · Page 73

GenAI4UQ: A software for forward and inverse uncertainty quantification using conditional generative AI

We introduce GenAI4UQ, a software package for forward and inverse uncertainty quantification in model calibration, parameter estimation, and ensemble forecasting. GenAI4UQ leverages a generative AI-based conditional modeling framework to address limitations of traditional inverse modeling techniques, such as Markov Chain Monte Carlo (MCMC) methods. By replacing computationally intensive iterative processes with a direct, learned mapping, GenAI4UQ enables efficient calibration of input parameters and generation of predictions directly from observations. The software supports rapid ensemble forecasting with robust uncertainty quantification while maintaining computational and storage efficiency. Built-in auto-tuning of hyperparameters simplifies model training, ensuring accessibility for users with varying expertise. Its versatile conditional generative framework is applicable across diverse scientific domains. While GenAI4UQ offers significant advantages in flexibility and efficiency, users should interpret its uncertainty estimates with caution in data-sparse scenarios, as the model may overestimate uncertainty—an effect common to all surrogate-based approaches including MCMC with surrogate models. Despite this, GenAI4UQ transforms inverse modeling by providing a fast, reliable, and user-friendly solution. It empowers researchers and practitioners to quickly estimate parameter distributions and generate model predictions for new observations, facilitating efficient decision-making and advancing the state of uncertainty quantification in computational modeling.

97 MATHEMATICS AND COMPUTING↗

Theoretical studies for the N{sub 2}–N{sub 2}O van der Waals complex: The potential energy surface, intermolecular vibrations, and rotational transition frequencies

Theoretical studies of the potential energy surface (PES) and bound states are performed for the N{sub 2}–N{sub 2}O van der Waals (vdW) complex. A four-dimensional intermolecular PES is constructed at the level of single and double excitation coupled-cluster method with a non-iterative perturbation treatment of triple excitations [CCSD(T)] with aug-cc-pVTZ basis set supplemented with bond functions. Two equivalent T-shaped global minima are located, in which the O atom of N{sub 2}O monomer is near the N{sub 2} monomer. The intermolecular fundamental vibrational states are assigned by inspecting the orientation of the nodal surface of the wavefunctions. The calculated frequency for intermolecular disrotation mode is 23.086 cm{sup −1}, which is in good agreement with the available experimental data of 22.334 cm{sup −1}. A negligible tunneling splitting with the value of 4.2 MHz is determined for the ground vibrational state and the tunneling splitting increases as the increment of the vibrational frequencies. Rotational levels and transition frequencies are calculated for both isotopomers {sup 14}N{sub 2}–N{sub 2}O and {sup 15}N{sub 2}–N{sub 2}O. The accuracy of the PES is validated by the good agreement between theoretical and experimental results for the transition frequencies and spectroscopic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Consistent solutions of the radiation diffusion equation in spherical and cylindrical geometries

We have extended the radiation diffusion model of Hammer and Rosen [Phys. Plasmas 10, 1829 (2003)] to diverging spherical and cylindrical geometries. The effect of curvilinear geometry on the supersonic, expanding wavefront increases as the internal radius of a spherical or cylindrical shell approaches zero. Small spherical geometries are important for modeling systems at the size scale of inertial confinement fusion capsules, at these scales existing quasi-analytic models for planar geometry significantly disagree with the results of simulation. With this method, the benefits of rapid iteration can be applied to common spherical systems at much smaller length scales. We present comparisons between numerical diffusion solutions and the analytic model to give ranges of applicability for the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GPU-BASED MONTE CARLO DUST RADIATIVE TRANSFER SCHEME APPLIED TO ACTIVE GALACTIC NUCLEI

A three-dimensional parallel Monte Carlo (MC) dust radiative transfer code is presented. To overcome the huge computing-time requirements of MC treatments, the computational power of vectorized hardware is used, utilizing either multi-core computer power or graphics processing units. The approach is a self-consistent way to solve the radiative transfer equation in arbitrary dust configurations. The code calculates the equilibrium temperatures of two populations of large grains and stochastic heated polycyclic aromatic hydrocarbons. Anisotropic scattering is treated applying the Heney-Greenstein phase function. The spectral energy distribution (SED) of the object is derived at low spatial resolution by a photon counting procedure and at high spatial resolution by a vectorized ray tracer. The latter allows computation of high signal-to-noise images of the objects at any frequencies and arbitrary viewing angles. We test the robustness of our approach against other radiative transfer codes. The SED and dust temperatures of one- and two-dimensional benchmarks are reproduced at high precision. The parallelization capability of various MC algorithms is analyzed and included in our treatment. We utilize the Lucy algorithm for the optical thin case where the Poisson noise is high, the iteration-free Bjorkman and Wood method to reduce the calculation time, and the Fleck and Canfield diffusion approximation for extreme optical thick cells. The code is applied to model the appearance of active galactic nuclei (AGNs) at optical and infrared wavelengths. The AGN torus is clumpy and includes fluffy composite grains of various sizes made up of silicates and carbon. The dependence of the SED on the number of clumps in the torus and the viewing angle is studied. The appearance of the 10 {mu}m silicate features in absorption or emission is discussed. The SED of the radio-loud quasar 3C 249.1 is fit by the AGN model and a cirrus component to account for the far-infrared emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Internal Delensing of Cosmic Microwave Background Polarization B -Modes with the POLARBEAR Experiment

Using only cosmic microwave background polarization data from the polarbear experiment, we measure B-mode polarization delensing on subdegree scales at more than 5σ significance. We achieve a 14% B-mode power variance reduction, the highest to date for internal delensing, and improve this result to 22% by applying for the first time an iterative maximum a posteriori delensing method. Our analysis demonstrates the capability of internal delensing as a means of improving constraints on inflationary models, paving the way for the optimal analysis of next-generation primordial B-mode experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Data analysis for synchrotron microcrystal native-SAD phasing

De novo structural evaluation of native biomolecules from single-wavelength anomalous diffraction (SAD) is a challenge because of the weakness of the anomalous scattering. For microcrystals with sizes less than 10 μm, native-SAD phasing at synchrotron microdiffraction beamlines is even more challenging due to difficulties in sample manipulation, diffraction data collection and data analysis. Native-SAD analysis from microcrystals by using X-ray free electron lasers has been demonstrated but has required use of thousands of thousands of microcrystals to achieve the necessary accuracy. Here we show that by exploitation of anomalous microdiffraction signals at a synchrotron beamline, by the use of polyimide wellmounts, and by an iterative crystal and frame rejection method, optimized data may be acquired and assembled which allows microcrystal native-SAD phasing from as few as about 1200 crystals. Finally, beyond microcrystals, our data analysis strategy may work equally well for microdiffraction data from larger crystals.

36 MATERIALS SCIENCE↗

Optimal Network Topology for Node-Breaker Representations With AC Power Flow Constraints

It has been demonstrated that network topology optimization (NTO) may change the topology of power system networks, and consequently, provide additional flexibility to reduce network congestion and violations. Most NTO problems are formulated based on the bus-branch model in which it is challenging to represent a realistic picture of all substation configurations. In this paper, we explore advantages of substation reconfiguration modeling based on node-breaker representations for NTO problem with full nonlinear alternating current power flow. It also proposes a tailored solution algorithm to solve this nonconvex mixed-integer nonlinear programming through the outer approximation method. The proposed solution approach iterates between a mixed-integer linear programming and a nonlinear subproblem. Additional enhancements to further accelerate the iteration process are illustrated. Numerical case studies demonstrate the relative economic and operational impact of optimal network topology with node-breaker representations.

42 ENGINEERING↗

Comparative Convergence Analysis of Nonlinear AMLI-Cycle Multigrid

The purpose of our paper is to provide a comprehensive convergence analysis of the nonlinear algebraic multilevel iteration (AMLI)-cycle multigrid (MG) method for symmetric positive definite problems. We show that the nonlinear AMLI-cycle MG method is uniformly convergent, based on classical assumptions for approximation and smoothing properties. Furthermore, under only the assumption that the smoother is convergent, we show that the nonlinear AMLI-cycle method is always better (or not worse) than the respective V-cycle MG method. Finally, numerical experiments are presented to illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗

Matlab Code for Multiscale Topology Optimization

This Matlab script implements a simple topology optimization algorithm that minimizes the compliance of a structure within a 3D prismatic design domain. It uses a simple linear elastic finite element solver to solve the equilibrium equations, and the optimality criterion (00) method to update the design iteratively. This functionality is identical to, or the simplest possible extension from 2D to 3D of, a 99-line Matlab code released in 2001 in a journal article by Sigmund. What is new in this code is our implementation of simple surrogate models for several truss micro-architectures, which converts this code to a multiscale design tool, whereas the original code is single-scale. The development of these surrogate models and their implementation is described in a journal article we intend to submit. We plan to include the full source code in an Appendix in this paper. The code’s primary function is as a pedagogical tool rather than as a finished software product. Sigmund’s 2001 paper describes the implementation line-by-line. We are including this updated version of the code to allow the reader to understand exactly what changes must be made to a “standard” topology optimization code to implement the surrogate models and thus multiscale design optimization.

Watts, Seth↗

Generalized Preconditioned Locally Harmonic Residual Eigensolver (GPLHR) v0.1

The software contains a MATLAB implementation of the Generalized Preconditioned Locally Harmonic Residual (GPLHR) method for solving standard and generalized non-Hermitian eigenproblems. The method is particularly useful for computing a subset of eigenvalues, and their eigen- or Schur vectors, closest to a given shift. The proposed method is based on block iterations and can take advantage of a preconditioner if it is available. It does not need to perform exact shift-and-invert transformation. Standard and generalized eigenproblems are handled in a unified framework.

Vecharynski, Eugene↗

Ultrasonic Model-Based Iterative Reconstruction with Spatially Variant Regularization for One-Sided Non-Destructive Evaluation

One-sided ultrasonic non-destructive evaluation (UNDE) uses ultrasound signals to investigate and inspect structures that are only accessible from one side. A widely used reconstruction technique in UNDE is the synthetic aperture focusing technique (SAFT). SAFT produces fast reconstruction and reasonable images for simple structures. However, for large complex structures, SAFT reconstructions suffer from noise and artifacts. To resolve some of the drawbacks of SAFT, an ultrasonic model-based iterative reconstruction (MBIR) algorithm, a method based on Bayesian estimation, was proposed that showed significant enhancement over SAFT in reducing noise and artifacts. In this paper, we build on previous investigations of the use of MBIR reconstruction on ultrasound data by proposing a spatially varying prior-model to account for artifacts from deeper regions and a 3D regularizer to account for correlations between scans from adjacent regions. We demonstrate that the use of the new prior model in MBIR can significantly improve reconstructions compared to SAFT and the previously proposed MBIR technique.

Almansouri, Hani↗

Denoising of polychromatic CT images based on their own noise properties

Purpose: Because of high diagnostic accuracy and fast scan time, computed tomography (CT) has been widely used in various clinical applications. Since the CT scan introduces radiation exposure to patients, however, dose reduction has recently been recognized as an important issue in CT imaging. However, low-dose CT causes an increase of noise in the image and thereby deteriorates the accuracy of diagnosis. In this paper, the authors develop an efficient denoising algorithm for low-dose CT images obtained using a polychromatic x-ray source. The algorithm is based on two steps: (i) estimation of space variant noise statistics, which are uniquely determined according to the system geometry and scanned object, and (ii) subsequent novel conversion of the estimated noise to Gaussian noise so that an existing high performance Gaussian noise filtering algorithm can be directly applied to CT images with non-Gaussian noise. Methods: For efficient polychromatic CT image denoising, the authors first reconstruct an image with the iterative maximum-likelihood polychromatic algorithm for CT to alleviate the beam-hardening problem. We then estimate the space-variant noise variance distribution on the image domain. Since there are many high performance denoising algorithms available for the Gaussian noise, image denoising can become much more efficient if they can be used. Hence, the authors propose a novel conversion scheme to transform the estimated space-variant noise to near Gaussian noise. In the suggested scheme, the authors first convert the image so that its mean and variance can have a linear relationship, and then produce a Gaussian image via variance stabilizing transform. The authors then apply a block matching 4D algorithm that is optimized for noise reduction of the Gaussian image, and reconvert the result to obtain a final denoised image. To examine the performance of the proposed method, an XCAT phantom simulation and a physical phantom experiment were conducted. Results: Both simulation and experimental results show that, unlike the existing denoising algorithms, the proposed algorithm can effectively reduce the noise over the whole region of CT images while preventing degradation of image resolution. Conclusions: To effectively denoise polychromatic low-dose CT images, a novel denoising algorithm is proposed. Because this algorithm is based on the noise statistics of a reconstructed polychromatic CT image, the spatially varying noise on the image is effectively reduced so that the denoised image will have homogeneous quality over the image domain. Through a simulation and a real experiment, it is verified that the proposed algorithm can deliver considerably better performance compared to the existing denoising algorithms.

60 APPLIED LIFE SCIENCES↗

Predicting objective function weights from patient anatomy in prostate IMRT treatment planning

Purpose: Intensity-modulated radiation therapy (IMRT) treatment planning typically combines multiple criteria into a single objective function by taking a weighted sum. The authors propose a statistical model that predicts objective function weights from patient anatomy for prostate IMRT treatment planning. This study provides a proof of concept for geometry-driven weight determination. Methods: A previously developed inverse optimization method (IOM) was used to generate optimal objective function weights for 24 patients using their historical treatment plans (i.e., dose distributions). These IOM weights were around 1% for each of the femoral heads, while bladder and rectum weights varied greatly between patients. A regression model was developed to predict a patient's rectum weight using the ratio of the overlap volume of the rectum and bladder with the planning target volume at a 1 cm expansion as the independent variable. The femoral head weights were fixed to 1% each and the bladder weight was calculated as one minus the rectum and femoral head weights. The model was validated using leave-one-out cross validation. Objective values and dose distributions generated through inverse planning using the predicted weights were compared to those generated using the original IOM weights, as well as an average of the IOM weights across all patients. Results: The IOM weight vectors were on average six times closer to the predicted weight vectors than to the average weight vector, usingl{sub 2} distance. Likewise, the bladder and rectum objective values achieved by the predicted weights were more similar to the objective values achieved by the IOM weights. The difference in objective value performance between the predicted and average weights was statistically significant according to a one-sided sign test. For all patients, the difference in rectum V54.3 Gy, rectum V70.0 Gy, bladder V54.3 Gy, and bladder V70.0 Gy values between the dose distributions generated by the predicted weights and IOM weights was less than 5 percentage points. Similarly, the difference in femoral head V54.3 Gy values between the two dose distributions was less than 5 percentage points for all but one patient. Conclusions: This study demonstrates a proof of concept that patient anatomy can be used to predict appropriate objective function weights for treatment planning. In the long term, such geometry-driven weights may serve as a starting point for iterative treatment plan design or may provide information about the most clinically relevant region of the Pareto surface to explore.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Robust optimization of intensity modulated proton therapy

Purpose: Intensity modulated proton therapy (IMPT) is highly sensitive to range uncertainties and uncertainties caused by setup variation. The conventional inverse treatment planning of IMPT optimized based on the planning target volume (PTV) is not often sufficient to ensure robustness of treatment plans. In this paper, a method that takes the uncertainties into account during plan optimization is used to mitigate the influence of uncertainties in IMPT. Methods: The authors use the so-called ''worst-case robust optimization'' to render IMPT plans robust in the face of uncertainties. For each iteration, nine different dose distributions are computed--one each for {+-} setup uncertainties along anteroposterior (A-P), lateral (R-L) and superior-inferior (S-I) directions, for {+-} range uncertainty, and the nominal dose distribution. The worst-case dose distribution is obtained by assigning the lowest dose among the nine doses to each voxel in the clinical target volume (CTV) and the highest dose to each voxel outside the CTV. Conceptually, the use of worst-case dose distribution is similar to the dose distribution achieved based on the use of PTV in traditional planning. The objective function value for a given iteration is computed using this worst-case dose distribution. The objective function used has been extended to further constrain the target dose inhomogeneity. Results: The worst-case robust optimization method is applied to a lung case, a skull base case, and a prostate case. Compared with IMPT plans optimized using conventional methods based on the PTV, our method yields plans that are considerably less sensitive to range and setup uncertainties. An interesting finding of the work presented here is that, in addition to reducing sensitivity to uncertainties, robust optimization also leads to improved optimality of treatment plans compared to the PTV-based optimization. This is reflected in reduction in plan scores and in the lower normal tissue doses for the same coverage of the target volume when subjected to uncertainties. Conclusions: The authors find that the worst-case robust optimization provides robust target coverage without sacrificing, and possibly even improving, the sparing of normal tissues. Our results demonstrate the importance of robust optimization. The authors assert that all IMPT plans should be robustly optimized.

61 RADIATION PROTECTION AND DOSIMETRY↗

Iterative X-ray spectroscopic ptychography

Spectroscopic ptychography is a powerful technique to determine the chemical composition of a sample with high spatial resolution. In spectro-ptychography, a sample is rastered through a focused X-ray beam with varying photon energy so that a series of phaseless diffraction data are recorded. Each chemical component in the material under investigation has a characteristic absorption and phase contrast as a function of photon energy. Using a dictionary formed by the set of contrast functions of each energy for each chemical component, it is possible to obtain the chemical composition of the material from high-resolution multi-spectral images. This paper presents SPA (spectroscopic ptychography with alternating direction method of multipliers), a novel algorithm to iteratively solve the spectroscopic blind ptychography problem. First, a nonlinear spectro-ptychography model based on Poisson maximum likelihood is designed, and then the proposed method is constructed on the basis of fast iterative splitting operators. SPA can be used to retrieve spectral contrast when considering either a known or an incomplete (partially known) dictionary of reference spectra. By coupling the redundancy across different spectral measurements, the proposed algorithm can achieve higher reconstruction quality when compared with standard state-of-the-art two-step methods. It is demonstrated how SPA can recover accurate chemical maps from Poisson-noised measurements, and its enhanced robustness when reconstructing reduced-redundancy ptychography data using large scanning step sizes is shown.

47 OTHER INSTRUMENTATION↗

TH-AB-207A-05: A Fully-Automated Pipeline for Generating CT Images Across a Range of Doses and Reconstruction Methods

Purpose: To evaluate the robustness of CAD or Quantitative Imaging methods, they should be tested on a variety of cases and under a variety of image acquisition and reconstruction conditions that represent the heterogeneity encountered in clinical practice. The purpose of this work was to develop a fully-automated pipeline for generating CT images that represent a wide range of dose and reconstruction conditions. Methods: The pipeline consists of three main modules: reduced-dose simulation, image reconstruction, and quantitative analysis. The first two modules of the pipeline can be operated in a completely automated fashion, using configuration files and running the modules in a batch queue. The input to the pipeline is raw projection CT data; this data is used to simulate different levels of dose reduction using a previously-published algorithm. Filtered-backprojection reconstructions are then performed using FreeCT-wFBP, a freely-available reconstruction software for helical CT. We also added support for an in-house, model-based iterative reconstruction algorithm using iterative coordinate-descent optimization, which may be run in tandem with the more conventional recon methods. The reduced-dose simulations and image reconstructions are controlled automatically by a single script, and they can be run in parallel on our research cluster. The pipeline was tested on phantom and lung screening datasets from a clinical scanner (Definition AS, Siemens Healthcare). Results: The images generated from our test datasets appeared to represent a realistic range of acquisition and reconstruction conditions that we would expect to find clinically. The time to generate images was approximately 30 minutes per dose/reconstruction combination on a hybrid CPU/GPU architecture. Conclusion: The automated research pipeline promises to be a useful tool for either training or evaluating performance of quantitative imaging software such as classifiers and CAD algorithms across the range of acquisition and reconstruction parameters present in the clinical environment. Funding support: NIH U01 CA181156; Disclosures (McNitt-Gray): Institutional research agreement, Siemens Healthcare; Past recipient, research grant support, Siemens Healthcare; Consultant, Toshiba America Medical Systems; Consultant, Samsung Electronics.

60 APPLIED LIFE SCIENCES↗

Development of a laser cleaning method for the first mirror surface of the charge exchange recombination spectroscopy diagnostics on ITER

A set of optical diagnostics is expected for measuring the plasma characteristics in ITER. Optical elements located inside discharge chambers are exposed to an intense radiation load, sputtering due to collisions with energetic atoms formed in the charge transfer processes, and contamination due to recondensation of materials sputtered from different parts of the construction of the chamber. Removing the films of the sputtered materials from the mirrors with the aid of pulsed laser radiation is an efficient cleaning method enabling recovery of the optical properties of the mirrors. In this work, we studied the efficiency of removal of metal oxide films by pulsed radiation of a fiber laser. Optimization of the laser cleaning conditions was carried out on samples representing metal substrates polished with optical quality with deposition of films on them imitating the chemical composition and conditions expected in ITER. It is shown that, by a proper selection of modes of radiation exposure to the surface with a deposited film, it is feasible to restore the original high reflection characteristics of optical elements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗