Search NASA⌕ Search

SEARCH · Search NASA

Results for “Iterative reconstruction”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Real-time sparse-sampled Ptychographic imaging through deep neural networks

Ptychography has rapidly grown in the fields of X-ray and electron imaging for its unprecedented ability to achieve nano or atomic scale resolution while simultaneously retrieving chemical or magnetic information from a sample. A ptychographic reconstruction is achieved by means of solving a complex inverse problem that imposes constraints both on the acquisition and on the analysis of the data, which typically precludes real-time imaging due to computational cost involved in solving this inverse problem. In this work we propose PtychoNN, a novel approach to solve the ptychography reconstruction problem based on deep convolutional neural networks. We demonstrate how the proposed method can be used to predict real-space structure and phase at each scan point solely from the corresponding far-field diffraction data. The presented results demonstrate how PtychoNN can effectively be used on experimental data, being able to generate high quality reconstructions of a sample up to hundreds of times faster than state-of-the-art ptychography reconstruction solutions once trained. By surpassing the typical constraints of iterative model-based methods, we can significantly relax the data acquisition sampling conditions and produce equally satisfactory reconstructions. Besides drastically accelerating acquisition and analysis, this capability can enable new imaging scenarios that were not possible before, in cases of dose sensitive, dynamic and extremely voluminous samples.

Cherukara, Mathew J↗
DOE OSTI2022

Iterative redox activation promotes interfacial synergy in an Ag/Cu x O catalyst for oxygen reduction

In this study, an iterative electrochemical activation process is demonstrated to maximize the interface between Ag and Cu x O. Morphological reconstruction and phase transformation of Cu x O occurs through the iterative electrochemical redox reaction, as indicated by energy-dispersive spectroscopy and X-ray photoelectron spectroscopy. The activated Ag/Cu x O/C catalysts exhibits enhanced oxygen reduction reaction performance, with an onset potential of 0.86 V vs. the reversible hydrogen electrode, a Tafel slope of 46mV dec -1 , and a high stability as compared to Ag/C. A mechanistic study using density functional theory shows that the weakened binding energy of the OH intermediate, originating from the charge transfer from Ag to Cu x O, improves the ORR activity at the interface sites of the Ag/Cu x O electrocatalysts. Ag/Cu x O shows a higher limiting potential of 0.14V for interface sites than isolated Ag nanoparticles. The interfacial charge transfer between Ag and CuO is verified experimentally.

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

Reconstructing magnetic deflections from sets of proton images using differential evolution

Proton imaging is a powerful technique for imaging electromagnetic fields within an experimental volume, in which spatial variations in proton fluence are a result of deflections to proton trajectories due to interaction with the fields. When deflections are large, proton trajectories can overlap, and this nonlinearity creates regions of greatly increased proton fluence on the image, known as caustics. The formation of caustics has been a persistent barrier to reconstructing the underlying fields from proton images. We have developed a new method for reconstructing the path-integrated magnetic fields, which begins to address the problem posed by caustics. Our method uses multiple proton images of the same object, each image at a different energy, to fill in the information gaps and provide some uniqueness when reconstructing caustic features. We use a differential evolution algorithm to iteratively estimate the underlying deflection function, which accurately reproduces the observed proton fluence at multiple proton energies simultaneously. We test this reconstruction method using synthetic proton images generated for three different, cylindrically symmetric field geometries at various field amplitudes and levels of proton statistics and present reconstruction results from a set of experimental images. Here, the method we propose requires no assumption of deflection linearity and can reliably solve for fields underlying linear, nonlinear, and caustic proton image features for the selected geometries and is shown to be fairly robust to noise in the input proton intensity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Learning neural representations for X-ray ptychography reconstruction with unknown probes

X-ray ptychography provides exceptional nanoscale resolution and is widely applied in materials science, biology, and nanotechnology. However, its full potential is constrained by the critical challenge of accurately reconstructing images when the illuminating probe is unknown. Conventional iterative methods and deep learning approaches are often suboptimal, particularly under the low-signal conditions inherent to low-dose and high-speed experiments. These limitations compromise reconstruction fidelity and restrict the broader adoption of the technique. In this work, we introduce the Ptychographic Implicit Neural Representation (PtyINR), a self-supervised framework that simultaneously addresses the object- and probe-recovery problem. By parameterizing both as continuous neural representations, PtyINR performs end-to-end reconstruction directly from raw diffraction patterns without requiring any pre-characterization of the probe. Extensive evaluations demonstrate that PtyINR achieves superior reconstruction quality on both simulated and experimental data, with remarkable robustness under challenging low-signal conditions. Furthermore, PtyINR offers a generalizable, physics-informed framework for addressing probe-dependent inverse problems, making it applicable to a wide range of computational microscopy problems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Watertight Tensor-Product Spline Reconstruction (KCNSC PDRD Technical Report)

The standard for representing geometric models in Computer Aided Design (CAD) systems is boundary representation with trimmed surfaces. However, the approximate nature of surface-surface intersections in CAD systems leaves many small gaps between mating surfaces. These gaps lead to problems for downstream analysis and manufacturing applications, requiring time consuming repair operations and further approximation of the geometry as polygonal meshes, which do not support design iteration. In this study, nVariate, Inc. evaluated its watertight reconstruction technology on a test part. The resulting watertight model has the same absolute model accuracy of the original model zero gaps, providing a representation that requires no repair and supports further design iterations in the CAD system.

42 ENGINEERING↗
DOE OSTI2022

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↗
DOE OSTI2021

Feasibility of DEIM for retrieving the initial field via dimensionality reduction

When parameter estimation is solved in a high-dimensional space, the dimensionality reduction strategy becomes the primary consideration for alleviating the tremendous computational cost. Here, the discrete empirical interpolation method (DEIM) is explored to retrieve the initial condition (IC) by combining the polynomial chaos (PC) based ensemble Kalman filter (i.e. PC-EnKF), where a non-intrusive PC expansion is considered as a surrogate model in place of the forward model in the prediction step of the ensemble Kalman filter, resulting in fewer forward model integrations but with a comparable accuracy as Monte Carlo-based approaches. The DEIM acts as a hyper-reduction tool to provide the low-dimensional input for the high-dimensional initial field, which can be reconstructed using the information on the sparse interpolation grid points that is adaptively obtained through PC-EnKF data assimilation method. Thus an innovative framework to reconstruct the IC is developed. The detailed procedure at each assimilation iteration includes: the determination of the spatial interpolation points, the estimation of the initial values on the interpolation locations using the optimal observations, and the reconstruction of IC in the full space. The current study uses the reconstruction field of initial conditions of the Navier-Stokes equations as an example to illustrate the efficacy of our method. The experimental results demonstrate the proposed algorithm achieves a satisfactory reconstruction for the initial field. The proposed method helps to extend the applicable area of DEIM in solving inverse problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Lagrangian displacement field estimators in cosmology

The late-time nonlinear Lagrangian displacement field is highly correlated with the initial field, so reconstructing it could enable us to extract primordial cosmological information. Our previous work carefully studied the displacement field reconstructed from the late-time density field using the iterative method proposed and found that it does not fully converge to the true, underlying displacement field (e.g., ~8 % offset at k ~ 0.2 h Mpc –1 at z = 0.6 ). We also constructed the Lagrangian perturbation theory model for the reconstructed field, but the model could not explain the discrepancy between the true and the reconstructed fields in the previous work. Here, the main sources of the discrepancy were speculated to be a numerical artifact in the displacement estimator due to the discreteness of the sample. In this paper, we develop two new estimators of the displacement fields to reduce such a numerical discreteness effect, the normalized momentum estimator and the rescaled resumed estimator. We show that the discrepancy reported is not due to the numerical artifacts. We conclude that the method cannot fully reconstruct the shape of the nonlinear displacement field at the redshift we studied, while it is still an efficient baryon acoustic oscillation reconstruction method. In parallel, by properly accounting for the UV-sensitive term in a reconstruction procedure with an effective field theory approach, we improve the theoretical model for the reconstructed displacement field, by almost 5 times, from ~15% to the level of a few percent at k ~ 0.2 h Mpc –1 at the redshift z = 0.6.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Lacustrine leaf wax hydrogen isotopes indicate strong regional climate feedbacks in Beringia since the last ice age

The Late-Quaternary climate of Beringia remains unresolved despite the region's role in modulating glacial-interglacial climate and as the likely conduit for human dispersal into the Americas. In this work, we investigate Beringian temperature change using an ~32,000-year lacustrine record of leaf wax hydrogen isotope ratios (δ2Hwax) from Arctic Alaska. Based on Monte Carlo iterations accounting for multiple sources of uncertainty, the reconstructed summertime temperatures were ~3 °C colder (range: –8 to +3 °C) during the Last Glacial Maximum (LGM; 21-25 ka) than the pre-industrial era (PI; 2–0.1 ka). This ice-age summer cooling is substantially smaller than in other parts of the Arctic, reflecting altered atmospheric circulation and increased continentality which weakened glacial cooling in the region. Deglacial warming was punctuated by abrupt events that are largely synchronous with events seen in Greenland ice cores that originate in the North Atlantic but which are also controlled locally, such as by the opening of the Bering Strait between 13.4 and 11 ka. Our reconstruction, together with climate modeling experiments, indicates that Beringia responds more strongly to North Atlantic freshwater forcing under modern-day, open-Bering Strait conditions than under glacial conditions. Furthermore, a 2 °C increase (Monte Carlo range: –1 to +5 °C) over the anthropogenic era reverses a 6 °C decline (Monte Carlo range: –10 to 0 °C) through the Holocene, indicating that recent warming in Arctic Alaska has not surpassed peak Holocene summer warmth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

AI-enabled high-resolution scanning coherent diffraction imaging

Ptychographic imaging is a powerful means of imaging beyond the resolution limits of typical x-ray optics. Capturing recovered images from raw ptychographic data, however, requires the solution of an inverse problem, namely, phase retrieval. Phase retrieval algorithms are computationally expensive, which precludes real-time imaging. In this work, we propose PtychoNN, an approach to solve the ptychography data inversion problem based on a deep convolutional neural network. We demonstrate how the proposed method can be used to predict real-space structure and phase at each scan point solely from the corresponding far-field diffraction data. Our results demonstrate the practical application of machine learning to recover high fidelity amplitude and phase contrast images of a real sample hundreds of times faster than current ptychography reconstruction packages. Furthermore, by overcoming the constraints of iterative model-based methods, we can significantly relax sampling constraints on data acquisition while still producing an excellent image of the sample. Besides drastically accelerating acquisition and analysis, this capability has profound implications for the imaging of dose sensitive, dynamic, and extremely voluminous samples.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Inversion of Many-Beam Bragg Intensities for Phasing by Iterated Projections: Removal of Multiple Scattering Artifacts from Diffraction Data

An iterated projection algorithm (N-Phaser) is developed that reconstructs a scattering potential from N-beam multiple Bragg scattered intensities. The method may be used to eliminate multiple scattering artifacts from electron diffraction data, solving the phase problem and increasing the thicknesses of samples used in materials science, solid-state chemistry, and small molecule crystallography. For high-energy transmission electron diffraction, we show that the algorithm recovers accurate complex structure factors from a wide range of thicknesses, orientations, and relativistic beam energies, and does not require known thickness or atomic-resolution data if sufficient multiple scattering occurs. Furthermore, extensions to Cryo-electron microscopy and Micro-electron diffraction are suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An empirical method for geometric calibration of a photon counting detector-based cone beam CT system

BACKGROUND: Geometric calibration is essential in developing a reliable computed tomography (CT) system. It involves estimating the geometry under which the angular projections are acquired. Geometric calibration of cone beam CTs employing small area detectors, such as currently available photon counting detectors (PCDs), is challenging when using traditional-based methods due to detectors’ limited areas. OBJECTIVE: This study presented an empirical method for the geometric calibration of small area PCD-based cone beam CT systems. METHODS: Unlike the traditional methods, we developed an iterative optimization procedure to determine geometric parameters using the reconstructed images of small metal ball bearings (BBs) embedded in a custom-built phantom. An objective function incorporating the sphericities and symmetries of the embedded BBs was defined to assess performance of the reconstruction algorithm with the given initial estimated set of geometric parameters. The optimal parameter values were those which minimized the objective function. The TIGRE toolbox was employed for fast tomographic reconstruction. To evaluate the proposed method, computer simulations were carried out using various numbers of spheres placed in various locations. Furthermore, efficacy of the method was experimentally assessed using a custom-made benchtop PCD-based cone beam CT. RESULTS: Computer simulations validated the accuracy and reproducibility of the proposed method. The precise estimation of the geometric parameters of the benchtop revealed high-quality imaging in CT reconstruction of a breast phantom. Within the phantom, the cylindrical holes, fibers, and speck groups were imaged in high fidelity. The CNR analysis further revealed the quantitative improvements of the reconstruction performed with the estimated parameters using the proposed method. CONCLUSION: Apart from the computational cost, we concluded that the method was easy to implement and robust.

Instruments & Instrumentation↗
DOE OSTI2023

Physics constrained unsupervised deep learning for rapid, high resolution scanning coherent diffraction reconstruction

By circumventing the resolution limitations of optics, coherent diffractive imaging (CDI) and ptychography are making their way into scientific fields ranging from X-ray imaging to astronomy. Yet, the need for time consuming iterative phase recovery hampers real-time imaging. While supervised deep learning strategies have increased reconstruction speed, they sacrifice image quality. Furthermore, these methods’ demand for extensive labeled training data is experimentally burdensome. Here, we propose an unsupervised physics-informed neural network reconstruction method, PtychoPINN, that retains the factor of 100-to-1000 speedup of deep learning-based reconstruction while improving reconstruction quality by combining the diffraction forward map with real-space constraints from overlapping measurements. In particular, PtychoPINN gains a factor of 4 in linear resolution and an 8 dB improvement in PSNR while also accruing improvements in generalizability and robustness. This blend of performance and computational efficiency offers exciting prospects for high-resolution real-time imaging in high-throughput environments such as X-ray free electron lasers (XFELs) and diffraction-limited light sources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Quantitative imaging and automated fuel pin identification for passive gamma emission tomography

Compliance of member States to the Treaty on the Non-Proliferation of Nuclear Weapons is monitored through nuclear safeguards. The Passive Gamma Emission Tomography (PGET) system is a novel instrument developed within the framework of the International Atomic Energy Agency (IAEA) project JNT 1510, which included the European Commission, Finland, Hungary and Sweden. The PGET is used for the verification of spent nuclear fuel stored in water pools. Advanced image reconstruction techniques are crucial for obtaining high-quality cross-sectional images of the spent-fuel bundle to allow inspectors of the IAEA to monitor nuclear material and promptly identify its diversion. In this work, we have developed a software suite to accurately reconstruct the spent-fuel cross sectional image, automatically identify present fuel rods, and estimate their activity. Unique image reconstruction challenges are posed by the measurement of spent fuel, due to its high activity and the self-attenuation. While the former is mitigated by detector physical collimation, we implemented a linear forward model to model the detector responses to the fuel rods inside the PGET, to account for the latter. The image reconstruction is performed by solving a regularized linear inverse problem using the fast-iterative shrinkage-thresholding algorithm. We have also implemented the traditional filtered back projection (FBP) method based on the inverse Radon transform for comparison and applied both methods to reconstruct images of simulated mockup fuel assemblies. Higher image resolution and fewer reconstruction artifacts were obtained with the inverse-problem approach, with the mean-square-error reduced by 50%, and the structural-similarity improved by 200%. We then used a convolutional neural network (CNN) to automatically identify the bundle type and extract the pin locations from the images; the estimated activity levels finally being compared with the ground truth. The proposed computational methods accurately estimated the activity levels of the present pins, with an associated uncertainty of approximately 5%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Numerical simulation study of impurity B transport during real-time B powder injection in EAST

To understand the beneficial effects of impurity boron (B) on the plasma, the transport of B species in both atomic and ionic states in the edge plasma and the ensuing evolution of the background plasma are simulated during real-time B powder injection by using two-dimensional fluid code SOLPS-ITER. The initial background plasma profiles for the simulation are reconstructed based on an EAST upper-single null discharge (shot No. 70601). With the injection rate of B atoms, 4.0 × 10 21 atoms/s, the simulation shows that effective charge number is less than 1.1 at the inner boundary R - R sep = -5.2 cm (background plasma density, 5.1 × 10 19 m -3 ) at the outer mid-plane indicating that the B impurity does not have much influence on the energy confinement. The present work finds that the impurity B exhibits very specific features in the edge plasma, mostly present in the narrow region just outside the separatrix. The particular distribution pattern can reduce the heat flux greatly, leading to the divertor detachment. The divertor detachment can be achieved easily under an injection rate of B powder in a reasonable range. This study suggests that the impurity B has a promising potential application in future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Compression of tokamak boundary plasma simulation data using a maximum volume algorithm for matrix skeleton decomposition

This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗