Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Trajectory optimization”

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 487 records · Page 27

HARMONY: Large-Scale Architecture Search for Efficient Hybrid Language Models

As large language models scale to trillions of parameters, their computational and memory requirements present critical challenges for efficient training and deployment. While Mixture of Experts (MoE) architectures enable efficient scaling through sparse parameter activation, and state-space models like Mamba offer linear-time complexity, principled methods for combining these paradigms remain undeveloped. We introduce HARMONY (Hybrid Architecture Research for Mamba, Optimized with Neural efficiencY), a multi-objective evolutionary neural architecture search framework for discovering efficient hybrid language models that integrate Transformer attention mechanisms, Mixture-of-Experts routing, and Mamba state-space components. Through large-scale distributed search using 16,384 MI250X GPUs on the Frontier supercomputer, HARMONY explores a comprehensive design space encompassing six attention variants (MHA, MQA, GQA, MLA, SWA, and Mamba-2), variable MoE configurations with both routed and shared experts, and extensive Mamba hyperparameters. Our framework discovers heterogeneous architectures that balance training performance with computational efficiency through multi-objective optimization incorporating latency penalties and fitness-based selection. Analysis of discovered architectures reveals that optimal hybrid designs favor heterogeneous component mixing rather than homogeneous patterns, with Mamba-2 and Multi-Head Latent Attention (MLA) emerging as preferred mechanisms. Discovered architectures demonstrate superior training efficiency: our best configuration achieves a final perplexity of 1.0874 with 2.38B parameters while processing 4,320 tokens/second, outperforming significantly larger manually designed models. Full-scale evaluation shows HARMONY's top architectures achieve better loss trajectories than equivalently-sized models using state-of-the-art configurations including Mixtral, Jamba, and Samba. Additionally, we demonstrate 91% weak scaling efficiency when training discovered 36B-parameter models across 1,024 GPUs. HARMONY is released as an open framework with comprehensive tools for building and training hybrid models using expert-data-pipeline parallelism, democratizing access to automated architecture design for next-generation language models.

Herron, Emily [ORNL] (ORCID:0000000273008172)↗

Measurement of helicon waves with phase contrast imaging on DIII-D – A theoretical feasibility study

A DIII-D high-beta H-mode discharge, with I p = 850 kA, B t = 2.1 T n e = 4 10 19 m –3 , has been designed to validate full wave modeling of helicon waves by optimizing the expected response of the Phase Contrast Imaging diagnostic. Helicon waves have been predicted to have high current drive efficiency off-axis without facing the accessibility issues of lower-hybrid waves. To test these predictions experimentally, DIII-D has recently commissioned a high-power helicon antenna. To confidently predict the behavior of helicon waves in future devices, measurements of their fundamental properties and validation against models will be essential. Phase contrast imaging (PCI) is an absolutely calibrated internal reference interferometer able to measure density fluctuations with radial wavenumbers k R between 1.5 cm –1 and 20 cm –1 . For helicon waves 2 cm –1 < k R < 10 cm –1 is expected, allowing PCI to measure their envelope and wavenumber spectrum. This makes PCI a powerful tool for the validation of state-of-the-art models, like the AORSA full wave code. AORSA is used to compute the density perturbations measured by the PCI with 2D calculations corresponding to 11 different toroidal mode numbers combined to resolve the trajectory of the helicon wave in 3D. This is necessary because the waves travel 120 degrees toroidally from the antenna to the PCI. Here, a cold plasma finite element model (CPFEM) [4] is used to predict propagation through the scrape-off layer. The result of the CPFEM model is connected to AORSA by creating an artificial Gaussian antenna on the last closed flux surface. PCI shows best results for waves with small vertical wavenumbers k z . Modeling the helicon waves for several past DIII-D experiments shows that k z is minimized if the intersection of the helicon and the PCI laser beams occurs in the midplane. For such an optimized scenario the predicted signal level is two orders of magnitude larger than the background density fluctuations arising from broadband turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Unraveling the impact of initial choices and in-loop interventions on learning dynamics in autonomous scanning probe microscopy

The current focus in Autonomous Experimentation (AE) is on developing robust workflows to conduct the AE effectively. This entails the need for well-defined approaches to guide the AE process, including strategies for hyperparameter tuning and high-level human interventions within the workflow loop. This paper presents a comprehensive analysis of the influence of initial experimental conditions and in-loop interventions on the learning dynamics of Deep Kernel Learning (DKL) within the realm of AE in scanning probe microscopy. We explore the concept of the “seed effect,” where the initial experiment setup has a substantial impact on the subsequent learning trajectory. Additionally, we introduce an approach of the seed point interventions in AE allowing the operator to influence the exploration process. Using a dataset from Piezoresponse Force Microscopy on PbTiO 3 thin films, we illustrate the impact of the “seed effect” and in-loop seed interventions on the effectiveness of DKL in predicting material properties. The study highlights the importance of initial choices and adaptive interventions in optimizing learning rates and enhancing the efficiency of automated material characterization. This work offers valuable insights into designing more robust and effective AE workflows in microscopy with potential applications across various characterization techniques.

47 OTHER INSTRUMENTATION↗

AEcroscopy: A Software–Hardware Framework Empowering Microscopy Toward Automated and Autonomous Experimentation

Microscopy has been pivotal in improving the understanding of structure-function relationships at the nanoscale and is by now ubiquitous in most characterization labs. However, traditional microscopy operations are still limited largely by a human-centric click-and-go paradigm utilizing vendor-provided software, which limits the scope, utility, efficiency, effectiveness, and at times reproducibility of microscopy experiments. Here, in this work, a coupled software–hardware platform is developed that consists of a software package termed AEcroscopy (short for Automated Experiments in Microscopy), along with a field-programmable-gate-array device with LabView-built customized acquisition scripts, which overcome these limitations and provide the necessary abstractions toward full automation of microscopy platforms. The platform works across multiple vendor devices on scanning probe microscopes and electron microscopes. It enables customized scan trajectories, processing functions that can be triggered locally or remotely on processing servers, user-defined excitation waveforms, standardization of data models, and completely seamless operation through simple Python commands to enable a plethora of microscopy experiments to be performed in a reproducible, automated manner. This platform can be readily coupled with existing machine-learning libraries and simulations, to provide automated decision-making and active theory-experiment optimization to turn microscopes from characterization tools to instruments capable of autonomous model refinement and physics discovery.

47 OTHER INSTRUMENTATION↗

Effects of Magnetic and Electric Field Uniformity on Coded Aperture Imaging Quality in a Cycloidal Mass Analyzer

Cycloidal mass analyzers are unique sector mass analyzers as they exhibit perfect double focusing, making them ideal for incorporating spatial aperture coding, which can increase the throughput of a mass analyzer without affecting the resolving power. However, the focusing properties of the cycloidal mass analyzer depend on the uniformity of the electric and magnetic fields. In this paper, finite element simulation and charged particle tracing were used to investigate the effect of field uniformity on imaging performance of a cycloidal mass analyzer. For the magnetic field, we evaluate a new permanent magnet geometry by comparing it to a traditional geometry. Results indicate that creating an aperture image in a cycloidal mass spectrometer with the same FWHM as the slit requires less than 1% variation in magnetic field strength along the ion trajectories. The new magnet design, called the opposed dipole magnet, has less than 1% field variation over an area approximately 62 × 65 mm; nearly twice the area available in a traditional design of similar size and weight. This allows ion imaging across larger detector arrays without loss of resolving power. In addition, we compare the aperture imaging quality of a traditionally used cycloidal mass spectrometer electric design with a new optimized design with improved field uniformity. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-J-61: Electrodynamics and Nano-Scale Fluid Dynamics in Protein Localization of Nuclear Pore Complexes

Purpose: To develop a simulation to catalyze a reevaluation of common assumptions about 3 dimensional diffusive processes and help cell biologists gain a more nuanced, intuitive understanding of the true physical hurdles of protein signaling cascades. Furthermore, to discuss the possibility of intracellular electrodynamics as a critical, unrecognized component of cellular biology and protein dynamics that is necessary for optimal information flow from the cell membrane to the nucleus. Methods: The Unity 3D gaming physics engine was used to build an accurate virtual scale model of the cytoplasm within a few hundred nanometers of the nuclear membrane. A cloud of simulated pERK proteins is controlled by the physics simulation, where diffusion is based on experimentally measured values and the electrodynamics are based on theoretical nano-fluid dynamics. The trajectories of pERK within the cytoplasm and through the 1250 nuclear pores on the nuclear surface is recorded and analyzed. Results: The simulation quickly demonstrates that pERKs moving solely by diffusion will rarely locate and come within capture distance of a nuclear pore. The addition of intracellular electrodynamics between charges on the nuclear pore complexes and on pERKs increases the number of successful translocations by allowing the electro-physical attractive effects to draw in pERKs from the cytoplasm. The effects of changes in intracellular shielding ion concentrations allowed for estimation of the “capture radius” under varying conditions. Conclusion: The simulation allows a shift in perspective that is paramount in attempting to communicate the scale and dynamics of intracellular protein cascade mechanics. This work has allowed researchers to more fully understand the parameters involved in intracellular electrodynamics, such as shielding anion concentration and protein charge. As these effects are still far below the spatial resolution of currently available measurement technology this simulation is invaluable to researchers to visualize and study current hypotheses concerning protein dynamics.

60 APPLIED LIFE SCIENCES↗

CBETA 4-pass Orbit Correction and Tolerance Study (CBETA Machine Note #037)

The main goal of this study is to find the tolerance of the CBETA FFAG lattice subjected to various errors. Common error sources include non-ideal magnet positioning and alignment, undesired multipole field strengths, and BPM offsets. In order to fix the beam trajectories (orbits) affected by the errors, dipole kickers have been allocated throughout the FFAG beamline to serve as correctors. Since there are four coexisting orbits in the FFAG, the correction scheme must aim to correct all four of them simultaneously. The simulation is performed using Bmad, developed by Cornell University to model relativistic beam dynamics in customized accelerator lattices. The core of orbit correction is the singular value decomposition (SVD) optimization. As the error magnitudes increase, the beam size and emittance after correction could undesirably increase, and our goal is to find the "individual limit" of each error type beyond which the increases become unacceptable. Typically this limit is reached before the orbits become un-correctable. Similarly we can also find the "combined limit" of various errors. Both limits will be properly defined in the later section.

43 PARTICLE ACCELERATORS↗

Ni and Co Incorporation in Forsterite: A Density Functional Theory Study with Hubbard Correction

Ni and Co are critical elements for the world economy and modern technologies. Mafic and ultramafic deposits represent low-grade yet abundant alternatives to traditional Ni and Co ores. Here, in this work, density functional theory (DFT) with the Hubbard U correction (DFT+U) was used to simulate the incorporation of Ni and Co in forsterite (Mg 2 SiO 4 ), the Mg endmember of olivine, a common mineral in mafic and ultramafic rocks. Hubbard U terms for Ni and Co were parametrized using a series of oxide, hydroxide, carbonate, silicate, and sulfide minerals relevant to extraction and recovery of Ni and Co from mafic and ultramafic deposits. Electronic, energetic, magnetic, and structural properties were considered in the parametrization. For each of Ni and Co, an effective Hubbard correction (U eff ) value that optimized agreement with either experimental data or a hybrid exchange-correlation functional for all of the minerals considered is reported. DFT+U ab initio molecular dynamics (AIMD) simulations of Ni and Co incorporated into the M1 and M2 octahedral sites of forsterite were then performed. Ni and Co substitution in the M1 site was more energetically favorable than substitution in the M2 site, in agreement with published partition coefficients. AIMD trajectories were used to compute extended X-ray absorption fine structure (EXAFS) spectra of Ni in the M1 and M2 sites for direct fitting to a published experimental spectrum of Ni in a natural San Carlos olivine sample. The results of the fit indicated that ordering of Ni in the M1 site was not as strong at the low Ni concentrations relevant to mafic and ultramafic silicate minerals as that at the higher concentrations of the Ni-Mg olivine solid solutions studied to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implosion dynamics measurements at the National Ignition Facility

Measurements have been made of the in-flight dynamics of imploding capsules indirectly driven by laser energies of 1-1.7 MJ at the National Ignition Facility [Miller et al., Nucl. Fusion 44, 228 (2004)]. These experiments were part of the National Ignition Campaign [Landen et al., Phys. Plasmas 18, 051002 (2011)] to iteratively optimize the inputs required to achieve thermonuclear ignition in the laboratory. Using gated or streaked hard x-ray radiography, a suite of ablator performance parameters, including the time-resolved radius, velocity, mass, and thickness, have been determined throughout the acceleration history of surrogate gas-filled implosions. These measurements have been used to establish a dynamically consistent model of the ablative drive history and shell compressibility throughout the implosion trajectory. First results showed that the peak velocity of the original 1.3-MJ Ge-doped polymer (CH) point design using Au hohlraums reached only 75% of the required ignition velocity. Several capsule, hohlraum, and laser pulse changes were then implemented to improve this and other aspects of implosion performance and a dedicated effort was undertaken to test the sensitivity of the ablative drive to the rise time and length of the main laser pulse. Changing to Si rather than Ge-doped inner ablator layers and increasing the pulse length together raised peak velocity to 93% {+-} 5% of the ignition goal using a 1.5 MJ, 420 TW pulse. Further lengthening the pulse so that the laser remained on until the capsule reached 30% (rather than 60%-70%) of its initial radius, reduced the shell thickness and improved the final fuel {rho}R on companion shots with a cryogenic hydrogen fuel layer. Improved drive efficiency was observed using U rather than Au hohlraums, which was expected, and by slowing the rise time of laser pulse, which was not. The effect of changing the Si-dopant concentration and distribution, as well as the effect of using a larger initial shell thickness were also examined, both of which indicated that instabilities seeded at the ablation front are a significant source of hydrodynamic mix into the central hot spot. Additionally, a direct test of the surrogacy of cryogenic fuel layered versus gas-filled targets was performed. Together all these measurements have established the fundamental ablative-rocket relationship describing the dependence of implosion velocity on fractional ablator mass remaining. This curve shows a lower-than-expected ablator mass at a given velocity, making the capsule more susceptible to feedthrough of instabilities from the ablation front into the fuel and hot spot. This combination of low velocity and low ablator mass indicates that reaching ignition on the NIF will require >20 {mu}m ({approx}10%) thicker targets and laser powers at or beyond facility limits.

36 MATERIALS SCIENCE↗

Future Aviation Biofuel Analysis Using the Biomass Scenario Model

This presentation summarizes the findings from two articles that use the Biomass Scenario Model (BSM) to explore potential future aviation biofuels scenarios. Newes et al. (2017) concluded that construction/build out capabilities and development of the feedstock market are key bottlenecks in the initial years. Displacement of jet fuel by 30% with biofuels by 2030 is possible, but several factors related to policy design-in the absence of high oil prices or policy uncertainty-contribute to the timing and magnitude of aviation biofuels production. First, incentives targeted toward jet fuel production, such as financial incentives (e.g., a producer tax credit or a CO2 tax) could be sufficient to reach six billion gallons. Second, R&D investment in pre-commercial technologies is needed to reduce the cost of production through learning-by-doing. Third, reduction of investment risk through loan guarantees and offtake agreements may allow production to ramp up more quickly through accelerating industry learning. Lewis et al. (2018) showed that both BSM and the Freight and Fuel Transportation Optimization Tool suggest that 200 million to 1 billion gallons per year of alternative jet fuel production are possible by 2030, given multiple incentives and a favorable investment climate. However, different capital costs and technology maturation rates in the two models will affect deployment of different fuel production technologies and therefore the feedstocks needed. Further collaboration on these modeling approaches would reduce methodological blind spots while providing insights into future industry trajectories.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Target tracking via recursive Bayesian state estimation in cognitive radar networks

To cope with complicated environments and stealthier targets, incorporating intelligence and cognition cycles into target tracking is of great importance in modern sensor network management. With remarkable advances in sensor techniques and deployable platforms, a sensing system has freedom to select a subset of available radars, plan their trajectories, and transmit designed waveforms. In this paper, we propose a general framework for single target tracking in cognitive networks of radars, including consideration of waveform design, path planning, and radar selection, which are separately but not jointly taken into account in existing work. The tracking procedure, built on the theories of dynamic graphical models (DGM) and recursive Bayesian state estimation (RBSE), is formulated as two iterative steps: (i) solving a combinatorial optimization problem to select the optimal subset of radars, waveforms, and locations for the next tracking instant, and (ii) acquiring the recursive Bayesian state estimation to accurately track the target. Further, an illustrative example introduces a specific scenario in 2-D space. Simulation results based on the scenario demonstrate that the proposed framework can accurately track the target under the management of the network of radars.

47 OTHER INSTRUMENTATION↗

Data-Enabled Predictive Control for Building HVAC Systems

Model predictive control is widely used as a control technology for the computation of optimal control inputs of building heating, ventilating, and air conditioning (HVAC) systems. However, both the benefits and widespread adoption of model predictive control (MPC) are hindered by the effort of model creation, calibration, and accuracy of the predictions. In this paper, we apply the data-enabled predictive control (DeePC) algorithm for designing controls for building HVAC systems. The algorithm solely depends on input/output data from the system to predict future state trajectories without the need for system identification. The algorithm relies on the idea that a vector space of all input–output trajectories of a discrete-time linear time-invariant (LTI) system is spanned by time-shifts of a single measured trajectory, given the input signal is persistently exciting. Closed-loop simulations using EnergyPlus are performed to demonstrate the approach. The simulated building modeled in EnergyPlus is a modified commercial large office prototype building served by an air handling unit-variable air volume HVAC system. Temperature setpoints of zones are used as control variables to minimize the HVAC energy cost of the building considering a time-of-use electricity rate structure. Furthermore, sensitivity analysis is conducted to gain insights into the effect of parameter tuning on DeePC performance. Simulation results are used to illustrate the performance of the algorithm and compare the algorithm with model-based MPC and occupancy-based setpoint controller. Overall, DeePC achieves similar performance compared to MPC for lower engineering effort.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Three-dimensional anisotropic adaptive filtering of projection data for noise reduction in cone beam CT

Purpose: The combination of quickly rotating C-arm gantry with digital flat panel has enabled the acquisition of three-dimensional data (3D) in the interventional suite. However, image quality is still somewhat limited since the hardware has not been optimized for CT imaging. Adaptive anisotropic filtering has the ability to improve image quality by reducing the noise level and therewith the radiation dose without introducing noticeable blurring. By applying the filtering prior to 3D reconstruction, noise-induced streak artifacts are reduced as compared to processing in the image domain. Methods: 3D anisotropic adaptive filtering was used to process an ensemble of 2D x-ray views acquired along a circular trajectory around an object. After arranging the input data into a 3D space (2D projections + angle), the orientation of structures was estimated using a set of differently oriented filters. The resulting tensor representation of local orientation was utilized to control the anisotropic filtering. Low-pass filtering is applied only along structures to maintain high spatial frequency components perpendicular to these. The evaluation of the proposed algorithm includes numerical simulations, phantom experiments, and in-vivo data which were acquired using an AXIOM Artis dTA C-arm system (Siemens AG, Healthcare Sector, Forchheim, Germany). Spatial resolution and noise levels were compared with and without adaptive filtering. A human observer study was carried out to evaluate low-contrast detectability. Results: The adaptive anisotropic filtering algorithm was found to significantly improve low-contrast detectability by reducing the noise level by half (reduction of the standard deviation in certain areas from 74 to 30 HU). Virtually no degradation of high contrast spatial resolution was observed in the modulation transfer function (MTF) analysis. Although the algorithm is computationally intensive, hardware acceleration using Nvidia's CUDA Interface provided an 8.9-fold speed-up of the processing (from 1336 to 150 s). Conclusions: Adaptive anisotropic filtering has the potential to substantially improve image quality and/or reduce the radiation dose required for obtaining 3D image data using cone beam CT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

High precision magnetic field measurement and mapping of the LEReC 180 degree bending magnet using very low field NMR with hall combined probe (140 - 350 gauss)

The Relativistic Heavy Ion Collider (RHIC) at BNL are using the Low Energy RHIC Electron Cooling (LEReC) to conduct experiments that search for the quantum chromodynamic (QCD) critical point. The first ever electron cooling based on the RF acceleration of electron beams was experimentally demonstrated on April 5, 2019 at LEReC at BNL. The first critical step in obtaining successful 3D non-magnetized cooling of the Au ion bunches in the RHIC cooling section was matching the electron beam energy with a relative error less than 5*10 -4 to the ion beam energy. Part of the LEReC beamline is a dipole magnet that bends the electron beam 180 degree. One of the most outstanding measurement challenges is that the dipole field is so low (≈200 G). Most of the existing NMR probes can only measure fields >400 G. Lower signal-to-noise ratio at low fields is requires the use of larger sample volumes. Working with CAYLAR, a NMR probe has been redesigned and optimized for these low field measurements with high resolution. We report the methods, challenges, and results for extensive magnetic field mappings of the 180 dipole magnet. A combination of NMR and Hall sensors has been successfully implemented to measure uniform field regimes inside the magnet center area and non-uniform field regimes at the magnet ends. Detailed measurement and mapping have been performed at five radii and five heights along the beam trajectory. Meanwhile a finite element magnetic modeling simulation of the magnet using Opera software has been performed. The calculated and measured data are compared, and the calculated data is good reference for measured data over long length mapping from magnet edge to center. The measured magnetic measurement data is directly useful for beam instrumentation, diagnostics and operation.

Magnetic Field Measurement↗

Speeding up particle track reconstruction using a parallel Kalman filter algorithm

One of the most computationally difficult problems expected for the High-Luminosity Large Hadron Collider (HL-LHC) is determining the trajectory of charged particles during event reconstruction. Algorithms used at the LHC today rely on Kalman filtering, which builds physical trajectories incrementally while incorporating material effects and error estimation. Recognizing the need for faster computational throughput, we have adapted Kalman-filter-based methods for highly parallel, many-core SIMD architectures that are now prevalent in high-performance hardware. In this paper, we discuss the design and performance of the improved tracking algorithm, referred to as mkFit. A key piece of the algorithm is the Matriplex library, containing dedicated code to optimally vectorize operations on small matrices. The physics performance of the mkFit algorithm is comparable to the nominal CMS tracking algorithm when reconstructing tracks from simulated proton-proton collisions within the CMS detector. We study the scaling of the algorithm as a function of the parallel resources utilized and find large speedups both from vectorization and multi-threading. mkFit achieves a speedup of a factor of 6 compared to the nominal algorithm when run in a single-threaded application within the CMS software framework.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-F-J-80: Deformable Image Registration for Residual Organ Motion Estimation in Respiratory Gated Treatments with Scanned Carbon Ion Beams

Purpose: At the Centro Nazionale di Adroterapia Oncologica (CNAO, Pavia, Italy) C-ions respiratory gated treatments of patients with abdominal tumours started in 2014. In these cases, the therapeutic dose is delivered around end-exhale. We propose the use of a respiratory motion model to evaluate residual tumour motion. Such a model requires motion fields obtained from deformable image registration (DIR) between 4DCT phases, estimating anatomical motion through interpolation. The aim of this work is to identify the optimal DIR technique to be integrated in the modeling pipeline. Methods: We used 4DCT datasets from 4 patients to test 4 DIR algorithms: Bspline, demons, log-domain and symmetric log domain diffeomorphic demons. We evaluate DIR performance in terms of registration accuracy (RMSE between registered images) and anatomical consistency of the motion field (Jacobian) when registering end-inhale to end-exhale. We subsequently employed the model to estimate the tumour trajectory within the ideal gating window. Results: Within the liver contour, the RMSE is in the range 31–46 HU for the best performing algorithm (Bspline) and 43–145 HU for the worst one (demons). The Jacobians featured zero negative voxels (which indicate singularities in the motion field) for the Bspline fields in 3 of 4 patients, whereas diffeomorphic demons fields showed a non-null number of negative voxels in every case. GTV motion in the gating window measured less than 7 mm for every patient, displaying a predominant superior-inferior (SI) component. Conclusion: The Bspline algorithm allows for acceptable DIR results in the abdominal region, exhibiting the property of anatomical consistency of the computed field. Computed trajectories are in agreement with clinical expectations (small and prevalent SI displacements), since patients lie wearing semi-rigid immobilizing masks. In future, the model could be used for retrospective estimation of organ motion during treatment, as guided by the breathing surrogate signal.

60 APPLIED LIFE SCIENCES↗

Network Algorithms for Detection of Radiation Sources

In support of national defense, Domestic Nuclear Detection Office’s (DNDO) Intelligent Radiation Sensor Systems (IRSS) program supported the development of networks of radiation counters for detecting, localizing and identifying low-level, hazardous radiation sources. Industry teams developed the first generation of such networks with tens of counters, and demonstrated several of their capabilities in indoor and outdoor characterization tests. Subsequently, these test measurements have been used in algorithm replays using various sub-networks of counters. Test measurements combined with algorithm outputs are used to extract Key Measurements and Benchmark (KMB) datasets. We present two selective analyses of these datasets: (a) a notional border monitoring scenario that highlights the benefits of a network of counters compared to individual detectors, and (b) new insights into the Sequential Probability Ratio Test (SPRT) detection method, which lead to its adaptations for improved detection. Using KMB datasets from an outdoor test, we construct a notional border monitoring scenario, wherein twelve 2”*2” NaI detectors are deployed on the periphery of 21*21meter square region. A Cs-137 (175 uCi) source is moved across this region, starting several meters from outside and finally moving away. The measurements from individual counters and the network were processed using replays of a particle filter algorithm developed under IRSS program. The algorithm outputs from KMB datasets clearly illustrate the benefits of combining measurements from all networked counters: the source was detected before it entered the region, during its trajectory inside, and until it moved several meters away. When individual counters are used for detection, the source was detected for much shorter durations, and sometimes was missed in the interior region. The application of SPRT for detecting radiation sources requires choosing the detection threshold, which in turn requires a source strength estimate, typically specified as a multiplier of the background radiation level. A judicious selection of this source multiplier is essential to achieve optimal detection probability at a specified false alarm rate. Typically, this threshold is chosen from the Receiver Operating Characteristic (ROC) by varying the source multiplier estimate. ROC is expected to have a monotonically increasing profile between the detection probability and false alarm rate. We derived ROCs for multiple indoor tests using KMB datasets, which revealed an unexpected loop shape: as the multiplier increases, detection probability and false alarm rate both increase until a limit, and then both contract. Consequently, two detection probabilities correspond to the same false alarm rate, and the higher is achieved at a lower multiplier, which is the desired operating point. Using the Chebyshev’s inequality we analytically confirm this shape. Then, we present two improved network-SPRT methods by (a) using the threshold off-set as a weighting factor for the binary decisions from individual detectors in a weighted majority voting fusion rule, and (b) applying a composite SPRT derived using measurements from all counters.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Poster — Thur Eve — 28: Enabling trajectory-based radiotherapy on a TrueBeam accelerator with the Eclipse treatment planning system

The TrueBeam linear accelerator platform has a developer's mode which permits the user dynamic control over many of the machine's mechanical and radiation systems. Using this research tool, synchronous couch and gantry motion can be programmed to simulate isocentric treatment with a shortened SAD, with benefits such as smaller projected MLC leaf widths and an increased dose rate. In this work, water tank measurements were used to commission a virtual linear accelerator with an 85 cm SAD in Eclipse, from which several arc-based radiotherapy treatments were generated, including an inverse optimized VMAT delivery. For each plan, the pertinent treatment delivery information was extracted from control points specified in the Eclipse-exported DICOM files using the pydicom package in Python, allowing construction of an XML control file. The dimensions of the jaws and MLC positions, defined for an 85 cm SAD in Eclipse, were scaled for delivery on a conventional SAD linear accelerator, and translational couch motion was added as a function of gantry angle to simulate delivery at 85 cm SAD. Ionization chamber and Gafchromic film measurements were used to compare the radiation delivery to dose calculations in Eclipse. With the exception of the VMAT delivery, ionization chamber measurements agreed within 3.3% of the Eclipse calculations. For the VMAT delivery, the ionization chamber was located in an inhomogeneous region, but gamma evaluation of the Gafchromic film plane resulted in a 94.5% passing rate using criteria of 3 mm/3%. The results indicate that Eclipse calculation infrastructure can be used.

07 ISOTOPE AND RADIATION SOURCES↗