SEARCH · Search NASA

Results for “matched field processing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Stereo Imaging Velocimetry

A system and a method for measuring three-dimensional velocities at a plurality of points in a fluid employing at least two cameras positioned approximately perpendicular to one another. The cameras are calibrated to accurately represent image coordinates in world coordinate system. The two-dimensional views of the cameras are recorded for image processing and centroid coordinate determination. Any overlapping particle clusters are decomposed into constituent centroids. The tracer particles are tracked on a two-dimensional basis and then stereo matched to obtain three-dimensional locations of the particles as a function of time so that velocities can be measured therefrom The stereo imaging velocimetry technique of the present invention provides a full-field. quantitative, three-dimensional map of any optically transparent fluid which is seeded with tracer particles.

McDowell, Mark↗

NASA NTRSDate not supplied

Edge Wear of the Advanced Electric Propulsion System Pole Covers

This work summarizes a wear test focused on characterizing the erosion of the edges and side walls of the 12-kW Advanced Electric Propulsion System (AEPS) pole covers as well as the impact of magnetic field strength on component erosion for operation at a discharge voltage of 600 V. Testing was performed with the AEPS Engineering Test Unit 2 (ETU-2) thruster and accumulated approximately 325 hours of operation at the 600 V/12 kW condition with the magnetic field strength set to 75% of its nominal value. Thruster performance was invariant throughout the wear test for all thruster throttle conditions and matched previous measurements at the nominal magnetic field strength. Contrary to past results, inner front pole cover erosion rates were invariant to magnetic field strength at the 600 V/12 kW condition and no erosion was measured on the outer front pole cover. The erosion profile of the downstream inner front pole cover faces was continuous across the cover, indicating no significant change in erosion processes at the pole cover edges. Side wall erosion was only detected on the cathode-facing surfaces of the inner front pole cover and was equal in magnitude to the rates measured on the adjacent downstream edges. Taken together, these results suggest that erosion of the pole cover edges and sidewalls is driven by cathode-borne ions and will ultimately not impact AEPS life estimates.

Hall thrusters↗

DOE OSTI2014

Few-view cone-beam CT reconstruction with deformed prior image

Purpose: Prior images can be incorporated into the image reconstruction process to improve the quality of subsequent cone-beam CT (CBCT) images from sparse-view or low-dose projections. The purpose of this work is to develop a deformed prior image-based reconstruction (DPIR) strategy to mitigate the deformation between the prior image and the target image. Methods: The deformed prior image is obtained by a projection-based registration approach. Specifically, the deformation vector fields used to deform the prior image are estimated through iteratively matching the forward projection of the deformed prior image and the measured on-treatment projections. The deformed prior image is then used as the prior image in the standard prior image constrained compressed sensing (PICCS) algorithm. A simulation study on an XCAT phantom and a clinical study on a head-and-neck cancer patient were conducted to evaluate the performance of the proposed DPIR strategy. Results: The deformed prior image matches the geometry of the on-treatment CBCT more closely as compared to the original prior image. Consequently, the performance of the DPIR strategy from few-view projections is improved in comparison to the standard PICCS algorithm, based on both visual inspection and quantitative measures. In the XCAT phantom study using 20 projections, the average root mean squared error is reduced from 14% in PICCS to 10% in DPIR, and the average universal quality index increases from 0.88 in PICCS to 0.92 in DPIR. Conclusions: The present DPIR approach provides a practical solution to the mismatch problem between the prior image and target image, which improves the performance of the original PICCS algorithm for CBCT reconstruction from few-view or low-dose projections.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2010

The MOeLLER experiment at Jefferson Lab: search for physics beyond the Standard Model

The MOeLLER experiment at Jefferson Lab will measure the parity-violating analyzing power A{sub z} in the scattering of 11 GeV longitudinally polarized electrons from the atomic electrons in a liquid hydrogen target (Moeller scattering). In the Standard Model a non-zero A{sub z} is due to the interference of the electromagnetic amplitude and the weak neutral current amplitude, the latter mediated by the Z{sup 0} boson. A{sub z} is predicted to be 35.6 parts per billion (ppb) at the kinematics of the experiment. It is the objective of the experiment to measure A{sub z} to a precision of 0.73 ppb. This result would yield a measurement of the weak charge of the electron Q{sub W}{sup e} to a fractional error of 2.3% at an average value Q{sup 2} of 0.0056 (GeV/c){sup 2}. This in turn will yield a determination of the weak mixing angle sin{sup 2{theta}}w with an uncertainty of {+-}0.00026(stat) {+-}0.00013(syst), comparable to the accuracy of the two best determinations at high energy colliders (at the Z{sup 0} pole). Consequently, the result could potentially influence the central value of this fundamental electroweak parameter, which is of critical importance in deciphering any signal of new physics that might be observed at the Large Hadron Collider (LHC). The measurement is sensitive to the interference of the electromagnetic amplitude with new neutral current amplitudes as weak as 10{sup -3} GF from as yet unknown high energy dynamics, a level of sensitivity unlikely to be matched in any experiment measuring a flavor and CP conserving process in the next decade. This provides indirect access to new physics at multi-TeV scales in a manner complementary to direct searches at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Electroweak splitting functions and high energy showering

We derive the electroweak (EW) collinear splitting functions for the Standard Model, including the massive fermions, gauge bosons and the Higgs boson. We first present the splitting functions in the limit of unbroken SU(2) L × U(1) Y and discuss their general features in the collinear and soft-collinear regimes. These are the leading contributions at a splitting scale ( k T ) far above the EW scale (v). We then systematically incorporate EW symmetry breaking (EWSB), which leads to the emergence of additional “ultra-collinear” splitting phenomena and naive violations of the Goldstone-boson Equivalence Theorem. We suggest a particularly convenient choice of non-covariant gauge (dubbed “Goldstone Equivalence Gauge”) that disentangles the effects of Goldstone bosons and gauge fields in the presence of EWSB, and allows trivial book-keeping of leading power corrections in v/ k T . We implement a comprehensive, practical EW showering scheme based on these splitting functions using a Sudakov evolution formalism. Novel features in the implementation include a complete accounting of ultra-collinear effects, matching between shower and decay, kinematic back-reaction corrections in multi-stage showers, and mixed-state evolution of neutral bosons ( γ/Z/h ) using density-matrices. We employ the EW showering formalism to study a number of important physical processes at O (1-10 TeV) energies. They include (a) electroweak partons in the initial state as the basis for vector-boson-fusion; (b) the emergence of “weak jets” such as those initiated by transverse gauge bosons, with individual splitting probabilities as large as O (35%); (c) EW showers initiated by top quarks, including Higgs bosons in the final state; (d) the occurrence of O (1) interference effects within EW showers involving the neutral bosons; and (e) EW corrections to new physics processes, as illustrated by production of a heavy vector boson ( W ') and the subsequent showering of its decay products.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Charged Wellbore Casing Controlled Source Electromagnetics (CWC-CSEM) for Reservoir Imaging and Monitoring (Final report)

This project addresses the needs of the U.S. Department of Energy (DOE) to develop advanced monitoring technologies and protocols to track the fate of subsurface carbon dioxide (CO2) plumes for carbon storage. Specifically, the project seeks to develop and test a unique and novel system of technologies consisting of electromagnetic data acquisition, coupled multiphysics imaging, and reservoir model enhancement to understand the migration and long-term distribution of CO2 in the subsurface. The overarching objective is to develop an integrated approach for long term monitoring of carbon storage. The two main components of the project include the methodology development and the test of the method at a field site. The methodology component consists of 1) developing the field procedure and protocol for collecting time-lapse controlled-source electromagnetic (CSEM) data with source electric current injected into the subsurface through wellbore casings; 2) building of background 3D electrical conductivity utilizing multiple sources of data such as supplemental surface transient EM (TEM) surveys, well-logs, and seismic structural information, for enhancing CSEM signal from reservoir depths; and 3) coupled multiphysics simulations and inversion of CSEM data constrained by production data and by structural information from seismic imaging of the reservoir and overlying formations. The testing component used the field site of Bell Creek Oil Field, which served both as a field laboratory for the method development as well as a test site to evaluate the CSEM signal strengths and the methodology developed in this research project. We have accomplished all the proposed tasks and developed the methodology as planned. These include the procedure for time-lapse CSEM data acquisition, data processing techniques, integration with 3D conductivity model building, fast reservoir simulation for history matching using machine learning, and interpreting CSEM data with coupling to the reservoir modeling. Collectively, the outcome of these tasks form a coherent workflow that can be applied to monitor dedicated carbon storage in saline reservoirs. The testing component evaluated the applicability and limitations of the method, and concluded that the method would be ideal for monitoring dedicated carbon storage sites utilizing saline reservoirs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Synthesis of results for Brine Availability Test in Salt (BATS) DECOVALEX-2023 Task E

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Stochastic Inflation at NNLO

Stochastic Inflation is an important framework for understanding the physics of de Sitter space and the phenomenology of inflation. In the leading approximation, this approach results in a Fokker-Planck equation that calculates the probability distribution for a light scalar field as a function of time. Despite its successes, the quantum field theoretic origins and the range of validity for this equation have remained elusive, and establishing a formalism to systematically incorporate higher order effects has been an area of active study. In this paper, we calculate the next-to-next-to-leading order (NNLO) corrections to Stochastic Inflation using Soft de Sitter Effective Theory (SdSET). In this effective description, Stochastic Inflation manifests as the renormalization group evolution of composite operators. The leading impact of non-Gaussian quantum fluctuations appears at NNLO, which is presented here for the first time; we derive the coefficient of this term from a two-loop anomalous dimension calculation within SdSET. We solve the resulting equation to determine the NNLO equilibrium distribution and the low-lying relaxation eigenvalues. In the process, we must match the UV theory onto SdSET at one-loop order, which provides a non-trivial confirmation that the separation into Wilson-coefficient corrections and contributions to initial conditions persists beyond tree level. Furthermore, these results illustrate how the naive factorization of time and momentum integrals in SdSET no longer holds in the presence of logarithmic divergences. It is these effects that ultimately give rise to the renormalization group flow that yields Stochastic Inflation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2015

Damage imaging in a laminated composite plate using an air-coupled time reversal mirror

We demonstrate the possibility of selectively imaging the features of a barely visible impact damage in a laminated composite plate by using an air-coupled time reversal mirror. The mirror consists of a number of piezoelectric transducers affixed to wedges of power law profiles, which act as unconventional matching layers. The transducers are enclosed in a hollow reverberant cavity with an opening to allow progressive emission of the ultrasonic wave field towards the composite plate. The principle of time reversal is used to focus elastic waves at each point of a scanning grid spanning the surface of the plate, thus allowing localized inspection at each of these points. The proposed device and signal processing removes the need to be in direct contact with the plate and reveals the same features as vibrothermography and more features than a C-scan. More importantly, this device can decouple the features of the defect according to their orientation, by selectively focusing vector components of motion into the object, through air. For instance, a delamination can be imaged in one experiment using out-of-plane focusing, whereas a crack can be imaged in a separate experiment using in-plane focusing. As a result, this capability, inherited from the principle of time reversal, cannot be found in conventional air-coupled transducers.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Numerical investigation of nanosecond pulsed discharge in air at above-atmospheric pressures

Here, this paper focuses on the multi-dimensional simulation of non-equilibrium plasma generated by nanosecond pulsed discharge in air, at pressure values higher than atmospheric. Voltage profiles and electrode geometry closely match those from a complementary experimental study. Simulations highlight the transition between different post-discharge plasma regimes at increasing pressure and tie the characteristics of the streamers to the electric field distribution in the gap between the electrodes. Results from simulations match experimental observations and qualitatively capture the experimental trend in terms of regime transition pressure and structure of the streamers. As a result, this paper validates a numerical tool that captures the physical and chemical properties of the low-temperature plasma and contributes to expand the understanding of low-temperature plasma ignition processes.

42 ENGINEERING↗

NASA NTRS2009

Measurement and Characterization of Space Shuttle Solid Rocket Motor Plume Acoustics

Lift-off acoustic environments generated by the future Ares I launch vehicle are assessed by the NASA Marshall Space Flight Center (MSFC) acoustics team using several prediction tools. This acoustic environment is directly caused by the Ares I First Stage booster, powered by the five-segment Reusable Solid Rocket Motor (RSRMV). The RSRMV is a larger-thrust derivative design from the currently used Space Shuttle solid rocket motor, the Reusable Solid Rocket Motor (RSRM). Lift-off acoustics is an integral part of the composite launch vibration environment affecting the Ares launch vehicle and must be assessed to help generate hardware qualification levels and ensure structural integrity of the vehicle during launch and lift-off. Available prediction tools that use free field noise source spectrums as a starting point for generation of lift-off acoustic environments are described in the monograph NASA SP-8072: "Acoustic Loads Generated by the Propulsion System." This monograph uses a reference database for free field noise source spectrums which consist of subscale rocket motor firings, oriented in horizontal static configurations. The phrase "subscale" is appropriate, since the thrust levels of rockets in the reference database are orders of magnitude lower than the current design thrust for the Ares launch family. Thus, extrapolation is needed to extend the various reference curves to match Ares-scale acoustic levels. This extrapolation process yields a subsequent amount of uncertainty added upon the acoustic environment predictions. As the Ares launch vehicle design schedule progresses, it is important to take every opportunity to lower prediction uncertainty and subsequently increase prediction accuracy. Never before in NASA s history has plume acoustics been measured for large scale solid rocket motors. Approximately twice a year, the RSRM prime vendor, ATK Launch Systems, static fires an assembled RSRM motor in a horizontal configuration at their test facility in Utah. The remaining RSRM static firings will take place on elevated terrain, with the nozzle exit plume being mostly undeflected and the landscape allowing placement of microphones within direct line of sight to the exhaust plume. These measurements will help assess the current extrapolation process by direct comparison between subscale and full scale solid rocket motor data.

Kenny, Jeremy↗

DOE OSTI2016

SU-G-IeP1-01: A Novel MRI Post-Processing Algorithm for Visualization of the Prostate LDR Brachytherapy Seeds and Calcifications Based On B0 Field Inhomogeneity Correction and Hough Transform

Purpose: This study aims at developing an MRI-only workflow for post-implant dosimetry of the prostate LDR brachytherapy seeds. The specific goal here is to develop a post-processing algorithm to produce positive contrast for the seeds and prostatic calcifications and differentiate between them on MR images. Methods: An agar-based phantom incorporating four dummy seeds (I-125) and five calcifications of different sizes (from sheep cortical bone) was constructed. Seeds were placed arbitrarily in the coronal plane. The phantom was scanned with 3T Philips Achieva MR scanner using an 8-channel head coil array. Multi-echo turbo spin echo (ME-TSE) and multi-echo gradient recalled echo (ME-GRE) sequences were acquired. Due to minimal susceptibility artifacts around seeds, ME-GRE sequence (flip angle=15; TR/TE=20/2.3/2.3; resolution=0.7×0.7×2mm3) was further processed.The induced field inhomogeneity due to the presence of titaniumencapsulated seeds was corrected using a B0 field map. B0 map was calculated using the ME-GRE sequence by calculating the phase difference at two different echo times. Initially, the product of the first echo and B0 map was calculated. The features corresponding to the seeds were then extracted in three steps: 1) the edge pixels were isolated using “Prewitt” operator; 2) the Hough transform was employed to detect ellipses approximately matching the dimensions of the seeds and 3) at the position and orientation of the detected ellipses an ellipse was drawn on the B0-corrected image. Results: The proposed B0-correction process produced positive contrast for the seeds and calcifications. The Hough transform based on Prewitt edge operator successfully identified all the seeds according to their ellipsoidal shape and dimensions in the edge image. Conclusion: The proposed post-processing algorithm successfully visualized the seeds and calcifications with positive contrast and differentiates between them according to their shapes. Further assessments on more realistic phantoms and patient study are required to validate the outcome.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Co-Optimization of CO 2 -EOR and Storage Processes in Mature Oil Reservoirs

This article presents an optimization methodology for CO 2 enhanced oil recovery in partially depleted reservoirs. A field-scale compositional reservoir flow model was developed for assessing the performance history of an active CO 2 flood and for optimizing both oil production and CO 2 storage in the Farnsworth Unit (FWU), Ochiltree County, Texas. A geological framework model constructed from geophysical, geological, and engineering data acquired from the FWU was the basis for all reservoir simulations and the optimization method. An equation of state was calibrated with laboratory fluid analyses and subsequently used to predict the thermodynamic minimum miscible pressure (MMP). Initial history calibrations of primary, secondary and tertiary recovery were conducted as the basis for the study. After a good match was achieved, an optimization approach consisting of a proxy or surrogate model was constructed with a polynomial response surface method (PRSM). The PRSM utilized an objective function that maximized both oil recovery and CO 2 storage. Experimental design was used to link uncertain parameters to the objective function. Control variables considered in this study included: water alternating gas cycle and ratio, production rates and bottom-hole pressure of injectors and producers. Other key parameters considered in the modeling process were CO 2 purchase, gas recycle and addition of infill wells and/or patterns. The PRSM proxy model was ‘trained’ or calibrated with a series of training simulations. This involved an iterative process until the surrogate model reached a specific validation criterion. A sensitivity analysis was first conducted to ascertain which of these control variables to retain in the surrogate model. A genetic algorithm with a mixed-integer capability optimization approach was employed to determine the optimum developmental strategy to maximize both oil recovery and CO 2 storage. The proxy model reduced the computational cost significantly. The validation criteria of the reduced order model ensured accuracy in the dynamic modeling results. The prediction outcome suggested robustness and reliability of the genetic algorithm for optimizing both oil recovery and CO 2 storage. The reservoir modeling approach used in this study illustrates an improved approach to optimizing oil production and CO 2 storage within partially depleted oil reservoirs such as FWU. Lastly, this study may serve as a benchmark for potential CO 2 –EOR projects in the Anadarko basin and/or geologically similar basins throughout the world.

02 PETROLEUM↗

DOE OSTI2020

Denoising atomic resolution 4D scanning transmission electron microscopy data with tensor singular value decomposition

Tensor singular value decomposition (SVD) is a method to find a low-dimensional representation of data with meaningful structure in three or more dimensions. Here, tensor SVD has been applied to denoise atomic-resolution 4D scanning transmission electron microscopy (4D STEM) data. On data simulated from a SrTiO 3 [100] perfect crystal and a Si [110] edge dislocation, tensor SVD achieved an average peak signal-to-noise ratio (PSNR) of ~40 dB, which matches or exceeds the performance of other denoising methods, with processing times at least 100 times shorter. On experimental data from SrTiO 3 [100] and LiZnSb [11 2 ¯ 0]/GaSb [110] samples, tensor SVD denoises multiple GB 4D STEM data sets in ten minutes on a typical personal computer. Denoising with tensor SVD improves both convergent beam electron diffraction patterns and virtual-aperture annular dark field images.

36 MATERIALS SCIENCE↗

NASA NTRS2011

First Results from ARTEMIS, a New Two-Spacecraft Lunar Mission: Counter-Streaming Plasma Populations in the Lunar Wake

We present observations from the first passage through the lunar plasma wake by one of two spacecraft comprising ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun), a new lunar mission that re-tasks two of five probes from the THEMIS magnetospheric mission. On Feb 13, 2010, ARTEMIS probe P1 passed through the wake at 3.5 lunar radii downstream from the Moon, in a region between those explored by Wind and the Lunar Prospector, Kaguya, Chandrayaan, and Chang'E missions. ARTEMIS observed interpenetrating proton, alpha particle, and electron populations refilling the wake along magnetic field lines from both flanks. The characteristics of these distributions match expectations from self-similar models of plasma expansion into vacuum, with an asymmetric character likely driven by a combination of a tilted interplanetary magnetic field and an anisotropic incident solar wind electron population. On this flyby, ARTEMIS provided unprecedented measurements of the interpenetrating beams of both electrons and ions naturally produced by the filtration and acceleration effects of electric fields set up during the refilling process. ARTEMIS also measured electrostatic oscillations closely correlated with counter-streaming electron beams in the wake, as previously hypothesized but never before directly measured. These observations demonstrate the capability of the comprehensively instrumented ARTEMIS spacecraft and the potential for new lunar science from this unique two spacecraft constellation.

Halekas, J. S.↗

NASA NTRS2014

First Results from ARTEMIS, A New Two-Spacecraft Lunar Mission: Counter-Streaming Plasma Populations in the Lunar Wake

We present observations from the first passage through the lunar plasma wake by one of two spacecraft comprising ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun), a new lunar mission that re-tasks two of five probes from the THEMIS magnetospheric mission. On Feb 13, 2010, ARTEMIS probe P1 passed through the wake at approximately 3.5 lunar radii downstream from the Moon, in a region between those explored by Wind and the Lunar Prospector, Kaguya, Chandrayaan, and Chang'E missions. ARTEMIS observed interpenetrating proton, alpha particle, and electron populations refilling the wake along magnetic field lines from both flanks. The characteristics of these distributions match expectations from self-similar models of plasma expansion into vacuum, with an asymmetric character likely driven by a combination of a tilted interplanetary magnetic field and an anisotropic incident solar wind electron population. On this flyby, ARTEMIS provided unprecedented measurements of the interpenetrating beams of both electrons and ions naturally produced by the filtration and acceleration effects of electric fields set up during the refilling process. ARTEMIS also measured electrostatic oscillations closely correlated with counter-streaming electron beams in the wake, as previously hypothesized but never before directly measured. These observations demonstrate the capability of the comprehensively instrumented ARTEMIS spacecraft and the potential for new lunar science from this unique two spacecraft constellation.

Plasma↗

DOE OSTI2015

A method for volumetric imaging in radiotherapy using single x-ray projection

Purpose: It is an intriguing problem to generate an instantaneous volumetric image based on the corresponding x-ray projection. The purpose of this study is to develop a new method to achieve this goal via a sparse learning approach. Methods: To extract motion information hidden in projection images, the authors partitioned a projection image into small rectangular patches. The authors utilized a sparse learning method to automatically select patches that have a high correlation with principal component analysis (PCA) coefficients of a lung motion model. A model that maps the patch intensity to the PCA coefficients was built along with the patch selection process. Based on this model, a measured projection can be used to predict the PCA coefficients, which are then further used to generate a motion vector field and hence a volumetric image. The authors have also proposed an intensity baseline correction method based on the partitioned projection, in which the first and the second moments of pixel intensities at a patch in a simulated projection image are matched with those in a measured one via a linear transformation. The proposed method has been validated in both simulated data and real phantom data. Results: The algorithm is able to identify patches that contain relevant motion information such as the diaphragm region. It is found that an intensity baseline correction step is important to remove the systematic error in the motion prediction. For the simulation case, the sparse learning model reduced the prediction error for the first PCA coefficient to 5%, compared to the 10% error when sparse learning was not used, and the 95th percentile error for the predicted motion vector was reduced from 2.40 to 0.92 mm. In the phantom case with a regular tumor motion, the predicted tumor trajectory was successfully reconstructed with a 0.82 mm error for tumor center localization compared to a 1.66 mm error without using the sparse learning method. When the tumor motion was driven by a real patient breathing signal with irregular periods and amplitudes, the average tumor center error was 0.6 mm. The algorithm robustness with respect to sparsity level, patch size, and presence or absence of diaphragm, as well as computation time, has also been studied. Conclusions: The authors have developed a new method that automatically identifies motion information from an x-ray projection, based on which a volumetric image is generated.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

The PFLOTRAN Reaction Sandbox

Abstract. As modern reactive transport simulators evolve to accommodate the demands of a user community, researchers need a platform for prototyping new biogeochemical processes, many of which are niche and specific to laboratory or field experiments. The PFLOTRAN Reaction Sandbox leverages modern, object-oriented Fortran in an attempt to provide such an environment within an existing reactive transport simulator. This work describes the PFLOTRAN Reaction Sandbox concept and implementation through several illustrative examples. Reaction Sandbox Biodegradation Hill customizes the existing microbially mediated biodegradation reaction formulation within PFLOTRAN to better match empirical data. Reaction Sandbox Simple provides an isolated environment for testing numerous preconfigured kinetic rate expressions and developing user intuition. Reaction Sandbox Example serves as a template for creating new sandboxes within PFLOTRAN.

58 GEOSCIENCES↗