Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tracking errors”

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 595 records · Page 33

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗

WE-AB-303-01: FEATURED PRESENTATION: A Dual-Detector Phase-Matched Digital Tomosynthesis (DTS) Imaging Scheme Using Aggregated KV and MV Projections for Intra-Treatment Lung Tumor Tracking

Purpose: To develop a dual-detector phase-matched DTS technique for continuous and fast intra-treatment lung tumor localization. Methods: Tumor localization accuracy of limited-angle DTS imaging is affected by low inter-slice resolution. The dual-detector DTS technique aims to overcome this limitation through combining orthogonally acquired beam’s eye view MV projections and kV projections for intra-treatment DTS reconstruction and localization. To aggregate the kV and MV projections for reconstruction, the MV projections were linearly converted to synthesize corresponding kV projections. To further address the lung motion induced localization errors, this technique uses respiratory phase-matching to match the motion information between on-board DTS and reference DTS to offset the adverse effects of motion blurriness in tumor localization.A study was performed using the CIRS008A lung phantom to simulate different on-board target variation scenarios for localization. The intra-treatment kV and MV acquisition was achieved through the Varian TrueBeam Developer Mode. Four methods were compared for their localization accuracy: 1. the proposed dual-detector phase-matched DTS technique; 2. the single-detector phase-matched DTS technique; 3. the dual-detector 3D-DTS technique without phase-matching; and 4. the single-detector 3D-DTS technique without phase-matching. Results: For scan angles of 2.5°, 5°, 10°, 20° and 30°, the dual-detector phase-matched DTS technique localized the tumor with average(±standard deviations) errors of 0.4±0.3 mm, 0.5±0.3 mm, 0.6±0.2 mm, 0.9±0.4 mm and 1.0±0.3 mm, respectively. The corresponding values of single-detector phase-matched DTS technique were 4.0±2.5 mm, 2.7±1.1 mm, 1.7±1.2 mm, 2.2±0.9 mm and 1.5±0.8 mm, respectively. The values of dual-detector 3D-DTS technique were 6.2±1.7 mm, 6.3±1.2 mm, 5.3±1.3 mm, 2.0±2.2 mm and 1.5±0.5 mm, respectively. And the values of single-detector 3D-DTS technique were 9.7±8.9 mm, 9.8±8.8 mm, 10.0±9.7 mm, 3.9±2.7 mm and 2.2±1.3 mm, respectively. Conclusion: The dual-detector phase-matched DTS technique substantially improves the tumor localization accuracy, which can be applied to real-time intra-treatment lung tumor localization. The research was funded by the National Institutes of Health Grant No. R01-CA184173 and a grant from Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

S2M: A Stochastic Simulation Model of Poliovirus Genetic State Transition

Modeling the molecular mechanisms that govern genetic variation can be useful in understanding the dynamics that drive genetic state transition in quasispecies viruses. For example, there is considerable interest in understanding how the relatively benign vaccine strains of poliovirus eventually revert to forms that confer neurovirulence and cause disease (ie, vaccine-derived poliovirus). This report describes a stochastic simulation model, S2M, which can be used to generate hypothetical outcomes based on known mechanisms of genetic diversity. S2M begins with predefined genotypes based on the Sabin-1 and Mahoney wild-type sequences, constructs a set of independent cell-based populations, and performs in-cell replication and cell-to-cell infection cycles while quantifying genetic changes that track the transition from Sabin-1 toward Mahoney. Realism is incorporated into the model by assigning defaults for variables that constrain mechanisms of genetic variability based roughly on metrics reported in the literature, yet these values can be modified at the command line in order to generate hypothetical outcomes driven by these parameters. To demonstrate the utility of S2M, simulations were performed to examine the effects of the rates of replication error and recombination and the presence or absence of defective interfering particles, upon reaching the end states of Mahoney resemblance (semblance of a vaccine-derived state), neurovirulence, genome fitness, and cloud diversity. Simulations provide insight into how modeled biological features may drive hypothetical outcomes, independently or in combination, in ways that are not always intuitively obvious.

59 BASIC BIOLOGICAL SCIENCES↗

Learning turbulent flows with generative models for super resolution and sparse flow reconstruction

Neural operators are promising surrogates for dynamical systems but when trained with standard L 2 losses they tend to oversmooth fine-scale turbulent structures. Here, we show that combining operator learning with generative modeling overcomes this limitation. We consider three practical turbulent-flow challenges where conventional neural operators fail: spatio-temporal super-resolution, forecasting, and sparse flow reconstruction. For Schlieren jet super-resolution, an adversarially trained neural operator (adv-NO) reduces the energy-spectrum error by 15 × while preserving sharp gradients at neural operator-like inference cost. For 3D homogeneous isotropic turbulence, adv-NO trained on only 160 timesteps from a single trajectory forecasts accurately for five eddy-turnover times and offers 114 × wall-clock speed-up at inference than the baseline diffusion-based forecasters, enabling near-real-time rollouts. For reconstructing cylinder wake flows from highly sparse Particle Tracking Velocimetry-like inputs, a conditional generative model infers full 3D velocity and pressure fields with correct phase alignment and statistics. These advances enable accurate reconstruction and forecasting at low compute cost, bringing near-real-time analysis and control within reach in experimental and computational fluid mechanics.

Fluid dynamics↗

Systematic non-LTE study of the −2.6 ≤ [Fe/H] ≤ 0.2 F and G dwarfs in the solar neighborhood. I. Stellar atmosphere parameters

We present atmospheric parameters for 51 nearby F and G dwarf and subgiant stars uniformly distributed over the −2.60<[Fe/H]<+0.20 metallicity range that is suitable for the Galactic chemical evolution research. Lines of iron in the two ionization stages, Fe i and Fe ii, were used to derive a homogeneous set of effective temperatures, surface gravities, iron abundances, and microturbulence velocities. Our spectroscopic analyses took advantage of employing high-resolution (R ≥ 60,000) Shane/Hamilton and Canada–France–Hawaii Telescope/ESPaDOnS observed spectra and non-LTE (NLTE) line formation for Fe i and Fe ii in the classical one-dimensional model atmospheres. The spectroscopic method was tested in advance with the 20 benchmark stars, for which there are multiple measurements of the infrared flux method effective temperature and their Hipparcos parallax error is less than 10%. We found NLTE abundances from lines of Fe i and Fe ii to be consistent within 0.06 dex for every benchmark star, when applying a scaling factor of S{sub H} = 0.5 to the Drawinian rates of inelastic Fe+H collisions. The obtained atmospheric parameters were checked for each program star by comparing its position in the log g–T{sub eff} plane with the theoretical evolutionary track of given metallicity and α-enhancement in the Yi et al. grid. Our final effective temperatures lie exactly in between the T{sub IRFM} scales of Alonso et al. and Casagrande et al., with a mean difference of +46 and −51 K, respectively. NLTE leads to higher surface gravity compared with that for LTE. The shift in log g is smaller than 0.1 dex for stars with [Fe/H] ⩾−0.75, T{sub eff} ≤ 5750 K, or log g ≥ 4.20. NLTE analysis is crucial for the very metal-poor turnoff and subgiant stars, for which the shift in log g between NLTE and LTE can be up to 0.5 dex. The obtained accurate atmospheric parameters will be used in the forthcoming papers to determine NLTE abundances of important astrophysical elements from lithium to europium and to improve observational constraints on the chemodynamical models of the Galaxy evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Modelling Curvilinear Beamline Effects on Beam- and Spin-Dynamics in the Fermilab Muon g-2 Storage Ring.

A model of the Fermilab Muon g-2 Experiment storage ring is developed with the BMAD simulation library for beam- and spin-dynamics studies. Particular focus is given to the electrostatic quadrupoles, the plates of which follow the curvature of the particle design trajectory of the magnetically uniform storage ring. To model these, standard electric multipoles and a novel method based on cylindrical harmonics are investigated and found to accurately match the electric field maps produced for the bulk and fringe regions of the curved quadrupoles respectively, despite these methods being specifically for straight beamlines. The toroidal multipole description is then developed, to include curvature and longitudinal field variations (to model the fringe field regions). Various combinations of the different descriptions (including use of the field maps) are used for the bulk and fringe fields of the quadrupoles, and the toroidal multipoles are found to provide the best agreement with the measured tunes of the storage ring. Two integrators are also developed, based on the method of Wu, Forest and Robin. The first of these describes motion in a uniform magnetic dipole field with arbitrary electric field, and produces results comparable to those of the adaptive-step fourth-order Runge–Kutta integrator implemented in BMAD, suitable for tracking over thousands of turns. The second is a generalisation of the first to arbitrary magnetic vector potentials superimposed on that of the uniform dipole, with arbitrary electric fields, and was also found to agree with the Runge–Kutta integrator. Simulations of antimuon storage for 2000 turns of the storage ring were performed, where losses were studied along with the evolution of beam polarisation. The storage losses, important for reducing the systematic error in the measurement of the muon g-2, were found to be largely due to collimator shape. The correlation between momentum and spin precession frequency reflect the effects expected from the electric field correction in the T-BMT equation, on the order of 10 -6 . Particle injection was also studied, for various injection angles, quadrupole strengths and strengths of injection kick, and in the cases where experimental data was provided, the model was found to be consistent with the experiment.

43 PARTICLE ACCELERATORS↗

Filtering adaptive output feedback control for multivariable nonlinear systems with mismatched uncertainties and unmodeled dynamics

This article synthesizes a filtering adaptive output feedback controller for multivariable nonlinear systems with mismatched uncertainties and unmodeled dynamics. The multivariable nonlinear systems under consideration have both matched and mismatched uncertainties, which satisfy the semiglobal Lipschitz condition. The unmodeled dynamics are bounded-input bounded-output stable. By adopting an estimation/cancellation strategy, a piecewise constant adaptive law drives the estimation error to zero at every time instant, which yields the adaptive parameters; a disturbance rejection control law is designed to compensate the nonlinear uncertainties within the bandwidth of low-pass filters. The matched uncertainties are cancelled directly by adopting their opposite in the control signal, while a dynamic inversion of the system is required to eliminate the effect of the mismatched uncertainties on the output. A feedforward control law is designed to track a given command. Since the virtual reference system defines the best performance that can be achieved by the closed-loop system, the uniform performance bounds are derived for the states and control signals via comparison. Both numerical and practical examples are provided to illustrate the effectiveness of the proposed filtering adaptive output feedback control architecture, comparisons with the model reference adaptive control demonstrates the superiority of the proposed control method.

filtering adaptive control, mismatched uncertainit↗

Fault isolation and fault-tolerant control for nonlinear stochastic distribution control systems with multiplicative faults

Here, in this paper, a fault isolation, diagnosis and fault tolerant control algorithm is proposed for nonlinear multiple multiplicative faults stochastic distribution control systems employing Takagi–Sugeno fuzzy system. To obtain the detailed fault information, a fault detection algorithm is introduced to discover the fault occurrence time. Then a fault isolation observer is built to produce the residual, and the error system is separated to subsystems affected only by disturbance and multiplicative faults. Moreover, a fault estimation scheme is presented to obtain the fault magnitude information. When faults occur, the system output probability density function will deviate from the desired distribution. So the model predictive control fault tolerant control scheme is needed to minimize the impact of faults as much as possible to make sure that the post fault output probability density function track the desired probability density function. The validity of the designed algorithm is demonstrated through a simulation example, where the fault tolerant control algorithm ensures that the system output probability density function still track the given output probability density function despite the complex case of multiple multiplicative faults occurring simultaneously.

42 ENGINEERING↗

A broadly-applicable unified closure relation for Taylor bubble rise velocity in pipes with stagnant liquid

Taylor bubble velocity in slug flow is a closure relation which significantly affects the prediction of liquid holdup (or void fraction) and pressure gradient in mechanistic models of slug flow for oil and gas pipe applications. In this work, we use a validated Computational Fluid Dynamics (CFD) approach to simulate the motion of Taylor bubbles in pipes; the interface is tracked with a Level-Set method implemented in a commercial code. A large numerical database is generated covering the most ample range of fluid properties and pipe inclination angles explored to date (Eo ∊ [10, 700],Mo ∊ [1 x 10 -6 , 5 x 10 3 , and θ ∊ [0°, 90°]). A unified Taylor bubble rise velocity correlation is extracted from that database. The new correlation predicts the numerical database with 8.6% absolute average relative error and a coefficient of determination R 2 = 0.97 and other available experimental data with 13.0% absolute average relative error and R 2 = 0.84 outperforming existing correlations and models.

Mechanics↗

LDRD FY25 Program Overview

As Lawrence Livermore National Laboratory’s (LLNL’s) Laboratory Directed Research and Development (LDRD) program enters its fifth decade of leading-edge research and development, its impact and importance have never been stronger. The program continues to advance strategic investments in pioneering science, technology, and engineering, ensuring LLNL will be ready to deliver on our mission as it evolves over the coming decades. Investing in LDRD research, and the people who perform this critical work, gives LLNL the ability to sustain our role as a leader in the Department of Energy and National Nuclear Security Administration enterprise. The LDRD program enables high-risk, high-payoff research that anticipates emerging threats and future mission needs. By nurturing the ingenuity of the Lab’s greatest asset, its people, LDRD funding advances not only our research but also grows and nurtures our workforce: engaging future innovators with student mentoring, challenging postdoctoral researchers to apply their skills to support national security, and strengthening the leadership skills of early career staff. This annual report documents how LDRD investments advance LLNL’s science, technology, and engineering across our mission space. To assess LDRD’s impact we track both short and long-term metrics such as peer-reviewed publications, number of students, or professional fellows. In addition to reviewing these metrics, I encourage you to delve deeper into the breadth of science and technology that illustrate the strategic value of this research portfolio. For instance, a recent exploratory research project used advanced manufacturing to construct miniaturized three-dimensional ion traps for a quantum computer with reduced quantum error rates to enable applications that address national security missions and support basic science. Another project has delved into studying detonation by examining deflagration to enhance the safety and security of the nuclear weapons stockpile. LDRD researchers are also deploying AI agents on two of the world’s most powerful supercomputers to automate and accelerate inertial confinement fusion experiments. Other teams are delivering more accurate optical constants to enable improved validation for aluminum to advance atomic and molecular physics models. LDRD-driven discoveries of how metals deform under extreme conditions strengthen our ability to model and design materials for demanding national security environments. National security challenges are increasingly complex and continuously evolving. LDRD focuses our most innovative science and technology on these challenges, ensuring the Laboratory is developing creative, forward-leaning solutions for our nation and the world. The following pages feature highlights of published scientific advances, patents, and honors that stem from LDRD investments. As you read this report, I hope you will understand how these investments position the Laboratory, and our partners, to meet the demands of the decades ahead.

36 MATERIALS SCIENCE↗

On bias of kinetic temperature measurements in complex plasmas

The kinetic temperature in complex plasmas is often measured using particle tracking velocimetry. Here, we introduce a criterion which minimizes the probability of faulty tracking of particles with normally distributed random displacements in consecutive frames. Faulty particle tracking results in a measurement bias of the deduced velocity distribution function and hence the deduced kinetic temperature. For particles with a normal velocity distribution function, mistracking biases the obtained velocity distribution function towards small velocities at the expense of large velocities, i.e., the inferred velocity distribution is more peaked and its tail is less pronounced. The kinetic temperature is therefore systematically underestimated in measurements. We give a prescription to mitigate this type of error.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impacts on cloud radiative effects induced by coexisting aerosols converted from international shipping and maritime DMS emissions

International shipping emissions (ISE), particularly sulfur dioxide, can influence the global radiation budget by interacting with clouds and radiation after being oxidized into sulfate aerosols. A better understanding of the uncertainties in estimating the cloud radiative effects (CREs) of ISE is of great importance in climate science. Many international shipping tracks cover oceans with substantial natural dimethyl sulfide (DMS) emissions. The interplay between these two major aerosol sources on CREs over vast oceanic regions with a relatively low aerosol concentration is an intriguing yet poorly addressed issue confounding estimation of the CREs of ISE. Using an Earth system model including two aerosol modules with different aerosol mixing configurations, we derive a significant global net CRE of ISE (-0.153 W m -2 with a standard error of ±0.004 W m -2 ) when using emissions consistent with current ship emission regulations. This global net CRE would become much weaker and actually insignificant (-0.001 W m -2 standard error of ±0.007 W m -2 ) if a more stringent regulation were adopted. We then reveal that the ISE-induced CRE would achieve a significant enhancement when a lower DMS emission is prescribed in the simulations, owing to the sublinear relationship between aerosol concentration and cloud response. In addition, this study also demonstrates that the representation of certain aerosol processes, such as mixing states, can influence the magnitude and pattern of the ISE-induced CRE. These findings suggest a reevaluation of the ISE-induced CRE with consideration of DMS variability.

54 ENVIRONMENTAL SCIENCES↗

TH-C-12A-08: New Compact 10 MV S-Band Linear Accelerator: 3D Finite-Element Design and Monte Carlo Dose Simulations

Purpose: To design a new compact S-band linac waveguide capable of producing a 10 MV x-ray beam, while maintaining the length (27.5 cm) of current 6 MV waveguides. This will allow higher x-ray energies to be used in our linac-MRI systems with the same footprint. Methods: Finite element software COMSOL Multiphysics was used to design an accelerator cavity matching one published in an experiment breakdown study, to ensure that our modeled cavities do not exceed the threshold electric fields published. This cavity was used as the basis for designing an accelerator waveguide, where each cavity of the full waveguide was tuned to resonate at 2.997 GHz by adjusting the cavity diameter. The RF field solution within the waveguide was calculated, and together with an electron-gun phase space generated using Opera3D/SCALA, were input into electron tracking software PARMELA to compute the electron phase space striking the x-ray target. This target phase space was then used in BEAM Monte Carlo simulations to generate percent depth doses curves for this new linac, which were then used to re-optimize the waveguide geometry. Results: The shunt impedance, Q-factor, and peak-to-mean electric field ratio were matched to those published for the breakdown study to within 0.1% error. After tuning the full waveguide, the peak surface fields are calculated to be 207 MV/m, 13% below the breakdown threshold, and a d-max depth of 2.42 cm, a D10/20 value of 1.59, compared to 2.45 cm and 1.59, respectively, for the simulated Varian 10 MV linac and brehmsstrahlung production efficiency 20% lower than a simulated Varian 10 MV linac. Conclusion: This work demonstrates the design of a functional 27.5 cm waveguide producing 10 MV photons with characteristics similar to a Varian 10 MV linac.

07 ISOTOPE AND RADIATION SOURCES↗

SU-F-J-55: Feasibility of Supraclavicular Field Treatment by Investigating Variation of Junction Position Between Breast Tangential and Supraclavicular Fields for Deep Inspiration Breath Hold (DIBH) Left Breast Radiation

Purpose: To explore the feasibility of supraclavicular field treatment by investigating the variation of junction position between tangential and supraclavicular fields during left breast radiation using DIBH technique. Methods: Six patients with left breast cancer treated using DIBH technique were included in this study. AlignRT system was used to track patient’s breast surface. During daily treatment, when the patient’s DIBH reached preset AlignRT tolerance of ±3mm for all principle directions (vertical, longitudinal, and lateral), the remaining longitudinal offset was recorded. The average with standard-deviation and the range of daily longitudinal offset for the entire treatment course were calculated for all six patients (93 fractions totally). The ranges of average ± 1σ and 2σ were calculated, and they represent longitudinal field edge error with the confidence level of 68% and 95%. Based on these longitudinal errors, dose at junction between breast tangential and supraclavicular fields with variable gap/overlap sizes was calculated as a percentage of prescription (on a representative patient treatment plan). Results: The average of longitudinal offset for all patients is 0.16±1.32mm, and the range of longitudinal offset is −2.6 to 2.6mm. The range of longitudinal field edge error at 68% confidence level is −1.48 to 1.16mm, and at 95% confidence level is −2.80 to 2.48mm. With a 5mm and 1mm gap, the junction dose could be as low as 37.5% and 84.9% of prescription dose; with a 5mm and 1mm overlap, the junction dose could be as high as 169.3% and 117.6%. Conclusion: We observed longitudinal field edge error at 95% confidence level is about ±2.5mm, and the junction dose could reach 70% hot/cold between different DIBH. However, over the entire course of treatment, the average junction variation for all patients is within 0.2mm. The results from our study shows it is potentially feasible to treat supraclavicular field with breast tangents.

60 APPLIED LIFE SCIENCES↗

A Data-Driven Framework for Direct Local Tensile Property Prediction of Laser Powder Bed Fusion Parts

This article proposes a generalizable, data-driven framework for qualifying laser powder bed fusion additively manufactured parts using part-specific in situ data, including powder bed imaging, machine health sensors, and laser scan paths. To achieve part qualification without relying solely on statistical processes or feedstock control, a sequence of machine learning models was trained on 6299 tensile specimens to locally predict the tensile properties of stainless-steel parts based on fused multi-modal in situ sensor data and a priori information. A cyberphysical infrastructure enabled the robust spatial tracking of individual specimens, and computer vision techniques registered the ground truth tensile measurements to the in situ data. The co-registered 230 GB dataset used in this work has been publicly released and is available as a set of HDF5 files. The extensive training data requirements and wide range of size scales were addressed by combining deep learning, machine learning, and feature engineering algorithms in a relay. The trained models demonstrated a 61% error reduction in ultimate tensile strength predictions relative to estimates made without any in situ information. Lessons learned and potential improvements to the sensors and mechanical testing procedure are discussed.

36 MATERIALS SCIENCE↗

The Road to Useful Quantum Computers

Building a useful quantum computer is a grand science and engineering challenge, currently pursued intensely by teams around the world. In the 1980s, Richard Feynman and Yuri Manin observed independently that computers based on quantum mechanics might enable better simulations of quantum phenomena. Their vision remained an intellectual curiosity until Peter Shor published his famous quantum algorithm for integer factoring, and shortly thereafter a proof that errors in quantum computations can be corrected. Since then, quantum computing R&D has progressed rapidly, from small-scale experiments in university physics laboratories to well-funded industrial efforts and prototypes. Hype notwithstanding, quantum computers have yet to solve scientifically or practically important problems -- a target often called quantum utility. In this article, we describe the capabilities of contemporary quantum computers, compare them to the requirements of quantum utility, and illustrate how to track progress from today to utility. We highlight key science and engineering challenges on the road to quantum utility, touching on relevant aspects of our own research.

Emerging Technologies (cs.ET)↗

On the Convergence of the Inexact Running Krasnosel'skii-Mann Method

This letter leverages a framework based on averaged operators to tackle the problem of tracking fixed points associated with maps that evolve over time. In particular, this letter considers the Krasnosel'skii-Mann (KM) method in a settings where: 1) the underlying map may change at each step of the algorithm, thus leading to a 'running' implementation of the KM method, and 2) an imperfect information of the map may be available. An imperfect knowledge of the maps can capture cases where processors feature a finite precision or quantization errors, or the case where (part of) the map is obtained from measurements. The analytical results are applicable to inexact running algorithms for solving optimization problems, whenever the algorithmic steps can be written in the form of (a composition of) averaged operators; examples are provided for inexact running gradient methods and the forward-backward splitting method. Convergence of the average fixed-point residual is investigated for the non-expansive case; linear convergence to a unique fixed-point trajectory is showen in the case of inexact running algorithms emerging from contractive operators.

97 MATHEMATICS AND COMPUTING↗

SU-E-J-58: Comparison of Conformal Tracking Methods Using Initial, Adaptive and Preceding Image Frames for Image Registration

Purpose: Accounting for tumor motion during radiation therapy is important to ensure that the tumor receives the prescribed dose. Increasing the field size to account for this motion exposes the surrounding healthy tissues to unnecessary radiation. In contrast to using motion-encompassing techniques to treat moving tumors, conformal radiation therapy (RT) uses a smaller field to track the tumor and adapts the beam aperture according to the motion detected. This work investigates and compares the performance of three markerless, EPID based, optical flow methods to track tumor motion with conformal RT. Methods: Three techniques were used to track the motions of a 3D printed lung tumor programmed to move according to the tumor of seven lung cancer patients. These techniques utilized a multi-resolution optical flow algorithm as the core computation for image registration. The first method (DIR) registers the incoming images with an initial reference frame, while the second method (RFSF) uses an adaptive reference frame and the third method (CU) uses preceding image frames for registration. The patient traces and errors were evaluated for the seven patients. Results: The average position errors for all patient traces were 0.12 ± 0.33 mm, −0.05 ± 0.04 mm and −0.28 ± 0.44 mm for CU, DIR and RFSF method respectively. The position errors distributed within 1 standard deviation are 0.74 mm, 0.37 mm and 0.96 mm respectively. The CU and RFSF algorithms are sensitive to the characteristics of the patient trace and produce a wider distribution of errors amongst patients. Although the mean error for the DIR method is negatively biased (−0.05 mm) for all patients, it has the narrowest distribution of position error, which can be corrected using an offset calibration. Conclusion: Three techniques of image registration and position update were studied. Using direct comparison with an initial frame yields the best performance. The authors would like to thank Dr.YeLin Suh for making the Cyberknife dataset available to us. Scholarship funding from the Natural Sciences and Engineering Research Council of Canada (NSERC) and CancerCare Manitoba Foundation is acknowledged.

60 APPLIED LIFE SCIENCES↗