Search NASA⌕ Search

SEARCH · Search NASA

Results for “under-determined inverse problem”

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.

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Decentralizing Network Inference Problems With Multiple-Description Fusion Estimation (MDFE)

Network inference (or tomography) problems, such as traffic matrix estimation or completion and link loss inference, have been studied rigorously in different networking applications. These problems are often posed as under-determined linear inverse (UDLI) problems and solved in a centralized manner, where all the measurements are collected at a central node, which then applies a variety of inference techniques to estimate the attributes of interest. This paper proposes a novel framework for decentralizing these large-scale under-determined network inference problems by intelligently partitioning it into smaller subproblems and solving them independently and in parallel. The resulting estimates, referred to as multiple descriptions, can then be fused together to compute the global estimate. We apply this Multiple Description and Fusion Estimation (MDFE) framework to three classical problems: traffic matrix estimation, traffic matrix completion, and loss inference. Using real topologies and traces, we demonstrate how MDFE can speed up computation while maintaining (even improving) the estimation accuracy and how it enhances robustness against noise and failures. In conclusion, we also show that our MDFE framework is compatible with a variety of existing inference techniques used to solve the UDLI problems.

97 MATHEMATICS AND COMPUTING↗

Accurate Holdup Calculations with Predictive Modeling & Data Integration

In facilities that process special nuclear material (SNM) it is important to account accurately for the fissile material that enters and leaves the plant. Although there are many stages and processes through which materials must be traced and measured, the focus of this project is material that is “held-up” in equipment, pipes, and ducts during normal operation and that can accumulate over time into significant quantities. Accurately estimating the holdup is essential for proper SNM accounting (vis-à-vis nuclear non-proliferation), criticality and radiation safety, waste management, and efficient plant operation. Usually it is not possible to directly measure the holdup quantity and location, so these must be inferred from measured radiation fields, primarily gamma and less frequently neutrons. Current methods to quantify holdup, i.e. Generalized Geometry Holdup (GGH), primarily rely on simple source configurations and crude radiation transport models aided by ad hoc correction factors. This project seeks an alternate method of performing measurement-based holdup calculations using a predictive model that employs state-of-the-art radiation transport codes capable of accurately simulating such situations. Inverse and data assimilation methods use the forward transport model to search for a source configuration that best matches the measured data and simultaneously provide an estimate of the level of confidence in the correctness of such configuration. In this work the holdup problem is re-interpreted as an inverse problem that is under-determined, hence may permit multiple solutions. A probabilistic approach is applied to solving the resulting inverse problem. This approach rates possible solutions according to their plausibility given the measurements and initial information. This is accomplished through the use of Bayes’ Theorem that resolves the issue of multiple solutions by giving an estimate of the probability of observing each possible solution. To use Bayes’ Theorem, one must have a model y(x) that maps the state variables x (the solution in this case) to the measurements y. In this case, the unknown state variables are the configuration and composition of the heldup SNM. The measurements are the detector readings. Thus, the natural model is neutral-particle radiation transport where a wealth of computational tools exists for performing these simulations accurately and efficiently. The combination of predictive model and Bayesian inference forms the Data Integration with Modeled Predictions (DIMP) method that serves as foundation for this project. The cost functional describing the model-to-data misfit is computed via a norm created by the inverse of the covariance matrix of the model parameters and responses. Since the model y(x) for the holdup problem is nonlinear, a nonlinear optimization on Q is conducted via Newton-type iterative methods to find the optimal values of the model parameters x. This project comprised a collaboration between NC State University (NCSU), the University of South Carolina (USC), and Oak Ridge National Laboratory (ORNL). The project was originally proposed in seven main tasks with an eighth contingency task to be performed if time and funding permitted; in fact time did not permit commencement of the contingency task and it was not performed. The remaining tasks involved holdup analysis with gamma detection strategies and separately with neutrons based on coincidence counting. Early in the project, and upon consultation with experts in coincidence counting it became evident that this approach is not viable for holdup applications and this task was replaced with an alternative, but valuable investigation that was carried out by the USC partner. Nevertheless, the experimental 4 measurements at ORNL of both gamma and neutron sources for the purpose of constructing Detector Response Functions (DRFs) with the associated uncertainties were indeed completed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Sparse measurement medical CT reconstruction using multi-fused block matching denoising priors

A major challenge for medical X-ray CT imaging is reducing the number of X-ray projections to lower radiation dosage and reduce scan times without compromising image quality. However these under-determined inverse imaging problems rely on the formulation of an expressive prior model to constrain the solution space while remaining computationally tractable. Traditional analytical reconstruction methods like Filtered Back Projection (FBP) often fail with sparse measurements, producing artifacts due to their reliance on the Shannon-Nyquist Sampling Theorem. Consensus Equilibrium, which is a generalization of Plug and Play, is a recent advancement in Model-Based Iterative Reconstruction (MBIR), has facilitated the use of multiple denoisers are prior models in an optimization free framework to capture complex, non-linear prior information. However, 3D prior modelling in a Plug and Play approach for volumetric image reconstruction requires long processing time due to high computing requirement. Instead of directly using a 3D prior, this work proposes a BM3D Multi Slice Fusion (BM3D-MSF) prior that uses multiple 2D image denoisers fused to act as a fully 3D prior model in Plug and Play reconstruction approach. Our approach does not require training and are thus able to circumvent ethical issues related with patient training data and are readily deployable in varying noise and measurement sparsity levels. In addition, reconstruction with the BM3D-MSF prior achieves similar reconstruction image quality as fully 3D image priors, but with significantly reduced computational complexity. We test our method on clinical CT data and demonstrate that our approach improves reconstructed image quality.

Hossain, Maliha [ORNL]↗

Improving Limited Angle CT Reconstruction with a Robust GAN Prior

Limited angle CT reconstruction is an under-determined linear inverse problem that requires appropriate regularization techniques to be solved. In this work we study how pre-trained generative adversarial networks (GANs) can be used to clean noisy, highly artifact laden reconstructions from conventional techniques, by effectively projecting onto the inferred image manifold. In particular, we use a robust version of the popularly used GAN prior for inverse problems, based on a recent technique called corruption mimicking, that significantly improves the reconstruction quality. The proposed approach operates in the image space directly, as a result of which it does not need to be trained or require access to the measurement model, is scanner agnostic, and can work over a wide range of sensing scenarios.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Regularization of soft-X-ray data in a tokamak application

An image inversion scheme for soft X-ray (SXR) diagnostics at the DIII-D tokamak is developed to obtain the local soft X-ray emission at a poloidal cross-section from the spatially line-integrated image taken by either a set of SXR arrays or a tangential camera. The scheme uses the Tikhonov regularization method since the inversion problem is generally ill-posed. The regularization technique uses the generalized singular value decomposition (GSVD) to determine a solution that depends on a free regularization parameter. The latter has to be chosen carefully, and the so called L -curve method to find the optimum regularization parameter is outlined. A representative test image is used to study the properties of the inversion scheme with respect to inversion accuracy, amount/strength of regularization and image resolution. The SXR arrays thereby represent an under-determined case while the tangential camera is an over-determined case of the linear inversion problem. In conclusion, it is found that the under-determined case suffers from lack of information but can still return a coarse, but fair image, while the camera data can be well inverted with high resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data Augmentation for Neutron Spectrum Unfolding with Neural Networks

Neural networks require a large quantity of training spectra and detector responses in order to learn to solve the inverse problem of neutron spectrum unfolding. In addition, due to the under-determined nature of unfolding, non-physical spectra which would not be encountered in usage should not be included in the training set. While physically realistic training spectra are commonly determined experimentally or generated through Monte Carlo simulation, this can become prohibitively expensive when considering the quantity of spectra needed to effectively train an unfolding network. In this paper, we present three algorithms for the generation of large quantities of realistic and physically motivated neutron energy spectra. Using an IAEA compendium of 251 spectra, we compare the unfolding performance of neural networks trained on spectra from these algorithms, when unfolding real-world spectra, to two baselines. We also investigate general methods for evaluating the performance of and optimizing feature engineering algorithms.

McGreivy, James (ORCID:0000000321723411)↗

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improved Bayesian regularization of inverse problems in vibrations and acoustics using noise-only measurements

Here, this paper studies Tikhonov regularization (ridge regression) parameter selection for problems in vibrations and acoustics. The selection method is based on a popular Bayesian method, but it incorporates measurements of sensor noise. The regularization parameter is closely related to the ratio of system input energy to noise energy, so noise measurements inform the inference procedure and improve parameter identification. In cases where standard Bayesian regularization identifies zero as the optimal regularization parameter, noise measurements guarantee a unique nonzero optimum. Sufficient theoretical criteria are developed for this guarantee. The method is verified in even-determined and under-determined configurations in an acoustic source localization simulation and a vibration load identification experiment. It is shown to yield significant improvements over existing empirical Bayesian regularization. Improvements are larger in the even-determined case and smaller in the under-determined case, wherein the inverse solution is less sensitive to the regularization parameter.

42 ENGINEERING↗