Search NASA⌕ Search

SEARCH · Search NASA

Results for “Optimal sizing method”

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 1,441 records · Page 80

Description of reaction and vibrational energetics of CO 2 –NH 3 interaction using quantum computing algorithms

CO 2 capture is critical to solving global warming. Amine-based solvents are extensively used to chemically absorb CO 2 . Thus, it is crucial to study the chemical absorption of CO 2 by amine-based solvents to better understand and optimize CO 2 capture processes. Here, we use quantum computing algorithms to quantify molecular vibrational energies and reaction pathways between CO 2 and a simplified amine-based solvent model—NH 3 . Molecular vibrational properties are important to understanding kinetics of reactions. However, the molecule size correlates with the strength of anharmonicity effect on vibrational properties, which can be challenging to address using classical computing. Quantum computing can help enhance molecular vibrational calculations by including anharmonicity. We implement a variational quantum eigensolver (VQE) algorithm in a quantum simulator to calculate ground state vibrational energies of reactants and products of the CO 2 and NH 3 reaction. The VQE calculations yield ground vibrational energies of CO 2 and NH 3 with similar accuracy to classical computing. In the presence of hardware noise, Compact Heuristic for Chemistry (CHC) ansatz with shallower circuit depth performs better than Unitary Vibrational Coupled Cluster. The “Zero Noise Extrapolation” error-mitigation approach in combination with CHC ansatz improves the vibrational calculation accuracy. Excited vibrational states are accessed with quantum equation of motion method for CO 2 and NH 3 . Using quantum Hartree–Fock (HF) embedding algorithm to calculate electronic energies, the corresponding reaction profile compares favorably with Coupled Cluster Singles and Doubles while being more accurate than HF. Our research showcases quantum computing applications in the study of CO 2 capture reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimization of Spirulina for biomanufacturing and the delivery of protein therapeutics

Lumen Bioscience has developed a novel cyanobacterial platform that enormously decreases the cost of pharmaceuticals used to prevent and treat illnesses. Arthrospira platensis (spirulina), is a photosynthetic microorganism that has been consumed as a dietary supplement for centuries in part because of its high protein content. Commercial cultivation operations have matured over the last 50 years to allow large-scale cultivation in outdoor ponds. Lumen recently discovered methods that transform this commercially important cyanobacterium into a genetically tractable platform for bioengineering. The high protein accumulation in spirulina makes it an ideal chassis for heterologous protein expression, and centuries of safe human consumption suggest its utility as an administration vehicle for biologic drugs. Lumen can express a broad range of recombinant proteins in spirulina and has developed manufacturing processes and strains to mucosally deliver bioactive proteins to treat and/or prevent disease. For example, Lumen has produced an oral cocktail of dried whole-cell spirulina biomass, containing 3 toxin-neutralizing antibodies and 1 endolysin, that is efficacious in a preclinical in vivo models of C. difficile infection. Lumen has also demonstrated that intranasal administration of a spirulina-manufactured neutralizing antibody can prevent disease by a respiratory pathogen, SARS-CoV-2, in a hamster model. Lumen’s current pipeline includes therapeutics to treat or prevent C. difficile infection, COVID-19, inflammatory bowel disease, cardiometabolic disease, and traveler’s diarrhea. Lumen deploys strategies to maximize the expression of these therapeutics for optimal dose sizing, a key limiting factor with past food crop-based expression systems. Spirulina can express exogenous proteins at unusually high levels (>15% of dry weight) but reaching this maximum value requires optimization of multiple interacting factors. We describe our current statistics design of experiments approach that optimizes these factors with minimal replicates. This maximizes exogenous protein expression allowing Lumen to generate spirulina-based therapeutics for a rapidly expanding range of medical uses.

Heinnickel, Mark↗

Seismic Tomographic Modeling of the Crust and Upper Mantle beneath Israel and the Middle East: Improved Resolution through Optimized Model Parameterization

Accurate regional seismic travel-time (RSTT) predictions rely on regional phases (e.g., Pg, Lg, Pn, Sn ) to account for 3D effects in the crust and upper mantle that are not captured by 1D models traditionally used for real-time location. The RSTT prediction model accounts for regional-scale crust and upper mantle structure globally by incorporating regional seismic phases into its travel-time calculations. Previous versions of the RSTT model have used a constant grid cell size of 1°. To improve the tomographic accuracy of recovering velocity structure at regional scales, we, in this study, perform data-driven grid refinement on the RSTT model down to a 0.125° grid (~14 km) in pursuit of two main goals: (1) to test the limits of RSTT capability and accuracy of determined velocity structure through variable grid refinement and (2) to image smaller structures in Israel and the Middle East and illuminate upper mantle dynamics operating in this complex tectonic area. We investigate the effects of model parameterization as grid cell size decreases and the trade-offs between recovered velocity structures. Our final dataset includes 4751 events and 499 stations that recorded 79,344 Pn and 7489 Pg . The variable grid refinement method allows recovery of finer-scale velocity structures and reduces travel-time residuals in areas with the highest data coverage. At smaller grid cell sizes, longer paths need to be upweighted to stabilize the inversion. Results illuminate tectonic features undefined in coarser grid-size models; in particular, we observe mantle perturbations related to the subduction zone around the Cyprian arc and crustal anomalies near the Dead Sea fault and throughout the Anatolian plate.

58 GEOSCIENCES↗

Research and Education of CO{sub 2} Separation from Coal Combustion Flue Gases with Regenerable Magnesium Solutions

A novel method using environment-friendly chemical magnesium hydroxide (Mg(OH) 2 ) solution to capture carbon dioxide from coal-fired power plants flue gas has been studied under this project in the post-combustion control area. The project utilizes the chemistry underlying the CO 2 -Mg(OH) 2 system and proven and well-studied mass transfer devices for high levels of CO 2 removal. The major goals of this research were to select and design an appropriate absorber which can absorb greater than 90% CO 2 gas with low energy costs, and to find and optimize the operating conditions for the regeneration step. During the project period, we studied the physical and chemical characteristics of the scrubbing agent, the reaction taking place in the system, development and evaluation of CO 2 gas absorber, desorption mechanism, and operation and optimization of continuous operation. Both batch and continuous operations were performed to examine the effects of various parameters including liquid-to-gas ratio, residence time, lean solvent concentration, pressure drop, bed height, CO 2 partial pressure, bubble size, pH, and temperature on the absorption.

20 FOSSIL-FUELED POWER PLANTS↗

Transport spatial model for the definition of green routes for city logistics centers

This paper presents a transport spatial decision support model (TSDSM) for carrying out the optimization of green routes for city logistics centers. The TSDSM model is based on the integration of the multi-criteria method of Weighted Linear Combination (WLC) and the modified Dijkstra algorithm within a geographic information system (GIS). The GIS is used for processing spatial data. The proposed model makes it possible to plan routes for green vehicles and maximize the positive effects on the environment, which can be seen in the reduction of harmful gas emissions and an increase in the air quality in highly populated areas. The scheduling of delivery vehicles is given as a problem of optimization in terms of the parameters of: the environment, health, use of space and logistics operating costs. Each of these input parameters was thoroughly examined and broken down in the GIS into criteria which further describe them. The model presented here takes into account the fact that logistics operators have a limited number of environmentally friendly (green) vehicles available. The TSDSM was tested on a network of roads with 127 links for the delivery of goods from the city logistics center to the user. The model supports any number of available environmentally friendly or environmentally unfriendly vehicles consistent with the size of the network and the transportation requirements. - Highlights: • Model for routing light delivery vehicles in urban areas. • Optimization of green routes for city logistics centers. • The proposed model maximizes the positive effects on the environment. • The model was tested on a real network.

54 ENVIRONMENTAL SCIENCES↗

Dosimetric comparison of flattened and unflattened beams for stereotactic ablative radiotherapy of stage I non-small cell lung cancer

Purpose: To compare contribution and accuracy of delivery for two flattening filter free (FFF) beams of the nominal energy 6 and 10 MV and a 6 MV flattened beam for early stage lung cancer. Methods: For each of 11 patients with stage I nonsmall cell lung cancer three volumetric modulated arc therapy plans were prepared utilizing a 6 MV flattened photon beam (X6FF) and two nonflattened beams of nominal energy 6 and 10 MV (X6FFF, X10FFF). Optimization constraints were set to produce dose distributions that meet the criteria of the RTOG-0915 protocol. The radiation schedule used for plan comparison in all patients was 50 Gy in five fractions. Dosimetric parameters of planning target volume (PTV) and organs-at-risk and delivery times were assessed and compared. All plans were subject to verification using Delta{sup 4} unit (Scandidos, Sweden) and absolutely calibrated gafchromic films in a thorax phantom. Results: All plans had a qualitatively comparable outcome. Obtained dose distributions were conformal (CI < 1.17) and exhibited a steep dose fall-off outside the PTV. The ratio of monitor units for FFF versus FF plans in the authors' study ranged from 0.95 to 1.21 and from 0.93 to 1.25 for X6FFF/X6FF and X10FFF/X6FF comparisons, respectively. The ratio systematically increased with increasing size of the PTV (up to +25% for 150 cm{sup 3} PTV). Yet the integral dose to healthy tissue did not follow this trend. Comparison of cumulative dose volume histograms for a patient's body showed that X6FFF plans exhibit improved conformity and reduced the volume of tissue that received more than 50% of the prescription dose. Parameters related to dose gradient showed statistically significant improvement. CI{sub 50%}, CI{sub 60%}, CI{sub 80%}, and CI{sub 100%} were on average reduced by 4.6% (p < 0.001), 4.6% (p = 0.002), 3.1% (p = 0.002), and 1.2% (p = 0.039), respectively. Gradient measure was on average reduced by 4.2% (p < 0.001). Due to dose reduction in the surrounding lung tissue, the V{sub 20} {sub Gy} and V{sub 12.5} {sub Gy} were reduced by 5.5% (p = 0.002) and 4.5% (p < 0.001). These dosimetric improvements in the fall-off were not observed for the X10FFF plans. Differences in sparing of normal tissues were not found to be statistically significant for either of the two FFF beams. Mean beam-on times were 111 s (2SD = 11 s) for X10FFF, 128 s (2SD = 19 s) for X6FFF, and X6FF plans required on average 269 s (2SD = 71 s). While the mean dose rate was 1555 ± 264 and 1368 ± 63 MU/min, for X10FFF and X6FFF, plans using the conventional X6FF were delivered with the constant maximum dose rate of 600 MU/min. Verification of all plans showed acceptable and comparable results for all plans in homogeneous as well as heterogeneous phantoms. Mean GS (3%, 2 mm) using the Delta{sup 4} phantom were 98.9% (2SD = 3.2%), 99.2% (2SD = 2.3%), and 99.2% (2SD = 2.3%) for X6FFF, X6FF, and X10FFF modalities. Verification using a thorax phantom showed GS > 98% in all cases. Conclusions: The use of FFF beams for stereotactic radiation therapy of nonsmall cell lung cancer patients yielded dose distributions qualitatively comparable to flattened beams and significantly reduced treatment delivery time. Utilizing the X6FFF beam improved conformity of dose distribution. On the other hand, X10FFF beam offered a slight improvement in treatment efficiency, and lower skin and peripheral dose. All effects were relatively small.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Sci-Sat AM: Radiation Dosimetry and Practical Therapy Solutions - 02: Dosimetric effects of gold nanoparticle surface coatings

Introduction: Gold nanoparticles (GNPs) can enhance radiation therapy within a tumour, increasing local energy deposition under irradiation, but experimental evidence suggests the enhancement is not as large as predicted by dose enhancement alone. Many studies neglect to account for surface coatings that are frequently used to optimize GNP uptake and biological distribution. This study uses Monte Carlo methods to investigate the consequences on local dose enhancement due to including these surface coatings. Methods: Using the PENELOPE Monte Carlo code system, GNP irradiation was simulated both with and without surface coatings of polyethylene glycol (PEG) of various molecular weights. Dose was scored to the gold, coating, and surrounding water, and the dosimetric differences between these scenarios were examined. Results: The simulated PEG coating absorbs a large portion of the energy that would otherwise be deposited in the medium. The mean dose to water was reduced by up to 2.5, 3.5, and 4.5% for GNPs of diameters 50, 20, and 10 nm, respectively. This effect was more pronounced for smaller GNPs, thicker coatings, and low photon source energies where the enhancement due to GNPs is the greatest. The molecular weight of the coating material did not have a significant impact on the dose. Conclusions: The inclusion of a coating material in GNP enhanced radiation may reduce the dose enhancement due to the nanoparticles. Both the composition and size of the coating play a role in the level of this reduction and should be considered carefully.

60 APPLIED LIFE SCIENCES↗

The EGS Collab Hydroshear Experiment at the Sanford Underground Research Facility – Siting Criteria and Evaluation of Candidate Sites

The objective of the EGS Collab project is to establish a suite of intermediate-scale (~10-20 m) field test beds coupled with stimulation and interwell flow tests to provide a basis to better understand fracture stimulation methods, resulting fracture geometries, and processes that control heat transfer between rock and stimulated fractures. Experiment 1 of the project is being conducted at a depth of ~1.5 km in the Sanford Underground Research Facility (SURF) on the 4850 Level (ft. below the ground surface). The stimulation method planned for Experiment 2 of this project is hydroshearing of an existing natural fracture. In siting this experiment, there are several key geologic criteria that need to be met. These include: 1) the fracture should be least 10 m in length so that it could be intersected by two boreholes that are spaced that far apart; 2) it should be optimally oriented relevant to the stress field so that it is critically stressed; 3) the site should have appropriate stress conditions (not too shallow a depth); 4) the fracture should have sufficient permeability to allow for prestimulation flow testing, but not too high a permeability that would preclude permeability enhancement via shear stimulation; 5) the fracture should not intersect other features (permeable fractures, boreholes, adjacent drifts) that could serve as major leak-off zones; 6) the site should not be complicated by geology (a single lithology would be preferred). In addition, there are logistical criteria that will also influence site selection, including: 1) site availability/access, 2) logistical support, e.g., power, internet, water, operating rail system; 3) drift size and orientation compatible with drilling and site operations; 4) appropriate ground support. There are only two deep levels at SURF that can be accessed by the Yates Shaft: the 4100 and the 4850, so we limited our search to these two levels. Currently, three candidate sites are being evaluated: the existing Experiment 1 site and its surroundings on the West Drift on the 4850 level, and two locations on the 4100 level. This paper will evaluate advantages, disadvantages, and risks associated with these three options at SURF.

EGS Collab, Sanford Underground Research Facility↗

Investigation of an Ultrafast Harmonic Resonant RF Kicker

An Energy Recovery Linac (ERL) based multi-turn electron Circulator Cooler Ring (CCR) is envisaged in the proposed Jefferson Lab Electron Ion Collider (JLEIC) to cool the ion bunches with high energy (55 MeV), high current (1.5 A), high repetition frequency (476.3 MHz), high quality magnetized electron bunches. A critical component in this scheme is a pair of ultrafast kickers for the exchange of electron bunches between the ERL and the CCR. The ultrafast kicker should operate with the rise and fall time in less than 2.1 ns, at the repetition rate of ~10s MHz, and should be able to run continuously during the whole period of cooling. These -and-fall time being combined together, are well beyond the state-of-art of traditional pulsed power supplies and magnet kickers. To solve this technical challenge, an alternative method is to generate this high repetition rate, fast rise-and-fall time short pulse continuous waveform by summing several finite number of (co)sine waves at harmonic frequencies of the kicking repetition frequency, and these harmonic modes can be generated by the Quarter Wave Resonater (QWR) based multifrequency cavities. Assuming the recirculator factor is 10, 10 harmonic modes (from 47.63 MHz to 476.3 MHz) with proper amplitudes and phases, plus a DC offset are combined together, a continuous short pulse waveform with the rise-and-fall time in less than 2.1 ns, repetition rate of 47.63 MHz waveform can be generated. With the compact and matured technology of QWR cavities, the total cost of both hardware development and operation can be reduced to a modest level. Focuse on the technical scheme, three main topics will be discussed in this thesis: the synthetization of the kicking pulse, the design and optimization of the deflecting QWR multi-integer harmonic frequency resonator and the fabrication and bench measurements of a half scale copper prototype. In the kicking pulse synthetization part, we begin with the Fourier Series expansion of an ideal square pulse, and get a Flat-Top waveform which will give a uniform kick over the bunch length of the kicked electron bunches, thus the transverse emittance of these kicked electron bunches can be maintained. By using two identical kickers with the betatron phase advance of 180 degree or its odd multiples, the residual kick voltage wave slopes at the unkicked bunch position will be totally cancelled out. Flat-Top waveform combined with two kicker scheme, the transverse emittance of the cooling electron bunches will be conserved during the whole injection, recirculation, and ejection processes. In the cavity design part, firstly, the cavity geometry is optimized to get high transverse shunt impedance thus less than 100 W of RF losses on the cavity wall can be achieved for all these 10 harmonic modes. To support all these 10 harmonic modes, group of four QWRs are adopted with the mode distribution of 5:3:1:1. In the multi-frequency cavities such as the five-mode-cavity and the three-mode-cavity, tunings are required to achieve the design frequencies for each mode. Slight segments of taper design on the inner conductor help to get the frequencies to be exactly on the odd harmonic modes. Stub tuners equal to the number of resonant modes are inserted to the outer conductor wall to compensate the frequency shifts due manufacturing errors and other perturbations during the operation such as the change of the cavity temperature. Single loop couple is designed for all harmonic modes in each cavity. By adjusting its loop size, position and rotation, it is possible to get the fundamental mode critical coupled and other higher harmonic modes slightly over coupled. A broadband circulator will be considered for absorbing the reflected power. Finally in this part, multipole field components due to the asymmetric cylindrical structure around the beam axis of the cavity as well as the beam-induced higher order mode (HOM) issues will be analyzed and discussed in this thesis. A half-scale copper prototype cavity (resonant frequencies from 95.26 MHz to 857.34 MHz) was fabricated to validate the electromagnetic characteristics. With this half scale prototype, the tuning processes of multiple harmonic frequencies, unloaded quality factor measurements of each mode, and bead-pull measurements are performed. The bench measurement results matched well with the simulation results, which have validated our cavity design and construction methods. Finally, a simple mode combining experiment with five separate signal generators was performed on this prototype cavity and the desired fast rise/fall time (1.2 ns), high repetition rate (95.26 MHz) waveform was captured, which finally proved our design of this ultrafast harmonic kicker.

43 PARTICLE ACCELERATORS↗

Tight-binding calculation studies of vacancy and adatom defects in graphene

Computational studies of complex defects in graphene usually need to deal with a larger number of atoms than the current first-principles methods can handle. We show a recently developed three-center tight-binding potential for carbon is very efficient for large scale atomistic simulations and can accurately describe the structures and energies of various defects in graphene. Using the three-center tight-binding potential, we have systematically studied the stable structures and formation energies of vacancy and embedded-atom defects of various sizes up to 4 vacancies and 4 embedded atoms in graphene. In conclusion, our calculations reveal low-energy defect structures and provide a more comprehensive understanding of the structures and stability of defects in graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scalable Techniques for Stochastic Power Flow Problems (Final Report)

The proposed research focuses on developing scalable algorithms for two-stage security-constrained OPF problems with AC power flow constraints, a class of problems complicated by (i) scale arising from a scenario representation; and (ii) the presence of nonlinearity, nonconvexity, and possibly second-stage discreteness or complementarity. Unfortunately, most existing solvers cannot contend with both challenges simultaneously; accordingly, the proposed research focuses on developing solution techniques that can both scale with the number of scenarios and contend with nonconvexity and second-stage complementarity. We consider three avenues for addressing such problems: (i) Variable sample-size SQP (VS-SQP) methods that combine sparse Quasi-Newton updates with a scalable variance-reduced stochastic gradient scheme for stochastic QP subproblems, allowing for contending with second-stage complementarity via regularization; (ii) Variable sample-size stochastic Interior-point (VS-sIP) schemes that propose a sampling-based regularized (to allow for contending with complementarity) interior-point schemes in which a Schur-complement technique is employed for decomposing the Newton direction computation step; (iii) Variable sample-size tractable ADMM (VS-tADMM) schemes combine variable sample-sizes with carefully designed techniques for resolving each of the nonconvex updates (by leveraging the QCQP structures). We intend to compare the three schemes using performance profiles in terms of solution quality, scalability, etc. and then select one scheme which will then be developed and further refined in Python for purposes of the GO competition.

42 ENGINEERING↗

Automated Quantitative Analysis of Extended Irradiation Defects - Dislocations, Voids and Precipitates in Neutron Irradiated HT-9 Steel

Structural alloys that make up the core components of nuclear reactors are exposed to radiation that can substantially alter the microstructure and mechanical properties of the material. Irradiation can cause the formation of extended defects including dislocation line and loops, voids and various types of precipitates. To understand the irradiation effects, it is critical to characterize these defects individually and collectively to quantify how they evolve throughout irradiation stages. Transmission electron microscopy (TEM) is a standard tool for providing extensive size, number density and spatial distribution of extended irradiation defects, however, manual measurement of defect size and number density requires time-intensive effort to analyze enough defects that provide a statistically meaningful result and are prone to errors of subjectivity. In this work, we developed a set of MATLAB algorithms offering robust and automated identification and quantification of extended irradiation defects like dislocation lines, voids and precipitates commonly observed in neutron irradiated cladding alloys [1]. Using a neutron-irradiated HT-9 martensitic alloy irradiated at 412°C to 111.8 dpa, we demonstrate the capability of these algorithms for quantitative defect analysis. Figure 1 presents defect images of the HT-9 using an advanced diffraction contrast imaging scanning transmission electron microscopy (DCI STEM) that suppresses bend contours and provides better defect contrast compared to conventional TEM [2]. The first step of defect analysis is defect identification, a process that can be digitally achieved by contrast thresholding (for simple defect contrast) and by machine learning computer vision for defect semantic segmentation [3]. Here, labeled binary defect images (Figures 2a and 2c) were used. For dislocation density quantification, we employed an optimized circular-grid intersection method first introduced by Dr. John Steeds [4]. This method is shown in Figures 2a and 2b, where the number of intersections of dislocation lines with circular traces was measured. This produces an unbiased sampling of dislocation lines, regardless of their orientation. This process was optimized with automated digital intersection number and grid length registration. Interestingly, it was found that a relatively large number of circles in a set circular gird is was needed to produce an invariant result (Figure 2b), with a low systematic measurement error. This demonstrates the main advantage of computer-based defect analysis, in performing repetitive operations quickly with a reduced systematic error. For the quantitative analysis of voids and precipitates a challenge encountered with automated defect analysis is feature overlap, particularly when the density of voids/precipitates is high. To mitigate this problem, as shown in Figure 2d, we modified the popular circular Hough transformation and treated ‘circular shape’ as the particle descriptor for voids/precipitates, efficiently resolving the overlapping particles. In summary, this automated digital quantification together with machine learning enabled defect identification [3] opens new opportunities for standardizing the characterization workflow for quantitative analysis of extended irradiation defects.

Sainju, Rajat↗

Joint constraints on galaxy bias and σ{sub 8} through the N-pdf of the galaxy number density

We present a full description of the N-probability density function of the galaxy number density fluctuations. This N-pdf is given in terms, on the one hand, of the cold dark matter correlations and, on the other hand, of the galaxy bias parameter. The method relies on the assumption commonly adopted that the dark matter density fluctuations follow a local non-linear transformation of the initial energy density perturbations. The N-pdf of the galaxy number density fluctuations allows for an optimal estimation of the bias parameter (e.g., via maximum-likelihood estimation, or Bayesian inference if there exists any a priori information on the bias parameter), and of those parameters defining the dark matter correlations, in particular its amplitude (σ{sub 8}). It also provides the proper framework to perform model selection between two competitive hypotheses. The parameters estimation capabilities of the N-pdf are proved by SDSS-like simulations (both, ideal log-normal simulations and mocks obtained from Las Damas simulations), showing that our estimator is unbiased. We apply our formalism to the 7th release of the SDSS main sample (for a volume-limited subset with absolute magnitudes M{sub r} ≤ −20). We obtain b-circumflex = 1.193 ± 0.074 and σ-bar{sub 8} = 0.862 ± 0.080, for galaxy number density fluctuations in cells of the size of 30h{sup −1}Mpc. Different model selection criteria show that galaxy biasing is clearly favoured.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-J-36: Comparison of CBCT Image Quality for Manufacturer Default Imaging Modes

Purpose CBCT is being increasingly used in patient setup for radiotherapy. Often the manufacturer default scan modes are used for performing these CBCT scans with the assumption that they are the best options. To quantitatively assess the image quality of these scan modes, all of the scan modes were tested as well as options with the reconstruction algorithm. Methods A CatPhan 504 phantom was scanned on a TrueBeam Linear Accelerator using the manufacturer scan modes (FSRT Head, Head, Image Gently, Pelvis, Pelvis Obese, Spotlight, & Thorax). The Head mode scan was then reconstructed multiple times with all filter options (Smooth, Standard, Sharp, & Ultra Sharp) and all Ring Suppression options (Disabled, Weak, Medium, & Strong). An open source ImageJ tool was created for analyzing the CatPhan 504 images. Results The MTF curve was primarily dictated by the voxel size and the filter used in the reconstruction algorithm. The filters also impact the image noise. The CNR was worst for the Image Gently mode, followed by FSRT Head and Head. The sharper the filter, the worse the CNR. HU varied significantly between scan modes. Pelvis Obese had lower than expected HU values than most while the Image Gently mode had higher than expected HU values. If a therapist tried to use preset window and level settings, they would not show the desired tissue for some scan modes. Conclusion Knowing the image quality of the set scan modes, will enable users to better optimize their setup CBCT. Evaluation of the scan mode image quality could improve setup efficiency and lead to better treatment outcomes.

60 APPLIED LIFE SCIENCES↗

Electron contamination modeling and reduction in a 1 T open bore inline MRI-linac system

Purpose: A potential side effect of inline MRI-linac systems is electron contamination focusing causing a high skin dose. In this work, the authors reexamine this prediction for an open bore 1 T MRI system being constructed for the Australian MRI-Linac Program. The efficiency of an electron contamination deflector (ECD) in purging electron contamination from the linac head is modeled, as well as the impact of a helium gas region between the deflector and phantom surface for lowering the amount of air-generated contamination. Methods: Magnetic modeling of the 1 T MRI was used to generate 3D magnetic field maps both with and without the presence of an ECD located immediately below the MLC’s. Forty-seven different ECD designs were modeled and for each the magnetic field map was imported into Geant4 Monte Carlo simulations including the linac head, ECD, and a 30 × 30 × 30 cm{sup 3} water phantom located at isocenter. For the first generation system, the x-ray source to isocenter distance (SID) will be 160 cm, resulting in an 81.2 cm long air gap from the base of the ECD to the phantom surface. The first 71.2 cm was modeled as air or helium gas, with the latter encased between two windows of 50 μm thick high density polyethlyene. 2D skin doses (at 70 μm depth) were calculated across the phantom surface at 1 × 1 mm{sup 2} resolution for 6 MV beams of field size of 5 × 5, 10 × 10, and 20 × 20 cm{sup 2}. Results: The skin dose was predicted to be of similar magnitude as the generic systems modeled in previous work, 230% to 1400% ofD {sub max} for 5 × 5 to 20 × 20 cm{sup 2}, respectively. Inclusion of the ECD introduced a nonuniformity to the MRI imaging field that ranged from ∼20 to ∼140 ppm while the net force acting on the ECD ranged from ∼151 N to ∼1773 N. Various ECD designs were 100% efficient at purging the electron contamination into the ECD magnet banks; however, a small percentage were scattered back into the beam and continued to the phantom surface. Replacing a large portion of the extended air-column between the ECD and phantom surface with helium gas is a key element as it significantly minimized the air-generated contamination. When using an optimal ECD and helium gas region, the 70 μm skin dose is predicted to increase moderately inside a small hot spot over that of the case with no magnetic field present for the jaw defined square beams examined here. These increases include from 12% to 40% of D {sub max} for 5 × 5 cm{sup 2}, 18% to 55% of D {sub max} for 10 × 10 cm{sup 2}, and from 23% to 65% of D {sub max} for 20 × 20 cm{sup 2}. Conclusions: Coupling an efficient ECD and helium gas region below the MLCs in the 160 cm isocenter MRI-linac system is predicted to ameliorate the impact electron contamination focusing has on skin dose increases. An ECD is practical as its impact on the MRI imaging distortion is correctable, and the mechanical forces acting on it manageable from an engineering point of view.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Dynamical Large Deviations of Two-Dimensional Kinetically Constrained Models Using a Neural-Network State Ansatz

We use a neural-network ansatz originally designed for the variational optimization of quantum systems to study dynamical large deviations in classical ones. We use recurrent neural networks to describe the large deviations of the dynamical activity of model glasses, kinetically constrained models in two dimensions. We present the first finite size-scaling analysis of the large-deviation functions of the two-dimensional Fredrickson-Andersen model, and explore the spatial structure of the high-activity sector of the South-or-East model. These results provide a new route to the study of dynamical large-deviation functions, and highlight the broad applicability of the neural-network state ansatz across domains in physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-E-I-60: Quality Assurance Testing Methods and Customized Phantom for Magnetic Resonance Imaging and Spectroscopy

Purpose: The objectives of this study are to develop an magnetic resonance imaging and spectroscopy (MRI-MRS) fused phantom along with the inserts for metabolite quantification and to conduct quantitative analysis and evaluation of the layered vials of brain-mimicking solution for quality assurance (QA) performance, according to the localization sequence. Methods: The outer cylindrical phantom body is made of acrylic materials. The section other than where the inner vials are located was filled with copper sulfate and diluted with water so as to reduce the T1 relaxation time. Sodium chloride was included to provide conductivity similar to the human body. All measurements of MRI and MRS were made using a 3.0 T scanner (Achiva Tx 3.0 T; Philips Medical Systems, Netherlands). The MRI scan parameters were as follows: (1) spin echo (SE) T1-weighted image: repetition time (TR), 500ms; echo time (TE), 20ms; matrix, 256×256; field of view (FOV), 250mm; gap, 1mm; number of signal averages (NSA), 1; (2) SE T2-weighted image: TR, 2,500 ms; TE, 80 ms; matrix, 256×256; FOV, 250mm; gap, 1mm; NSA, 1; 23 slice images were obtained with slice thickness of 5mm. The water signal of each volume of interest was suppressed by variable pulse power and optimized relaxation delays (VAPOR) applied before the scan. By applying a point-resolved spectroscopy sequence, the MRS scan parameters were as follows: voxel size, 0.8×0.8×0.8 cm{sup 3}; TR, 2,000ms; TE, 35ms; NSA, 128. Results: Using the fused phantom, the results of measuring MRI factors were: geometric distortion, <2% and ±2 mm; image intensity uniformity, 83.09±1.33%; percent-signal ghosting, 0.025±0.004; low-contrast object detectability, 27.85±0.80. In addition, the signal-to-noise ratio of N-acetyl-aspartate was consistently high (42.00±5.66). Conclusion: The MRI-MRS QA factors obtained simultaneously using the phantom can facilitate evaluation of both images and spectra, and provide guidelines for obtaining MRI and MRS QA factors simultaneously. This study was supported by grant (2012-007883 and 2014R1A2A1A10050270) from the Mid-career Researcher Program through the NRF funded by Ministry of Science. In addition, this study was supported by the Industrial R&D of MOTIE/KEIT (10048997, Development of the core technology for integrated therapy devices based on real-time MRI-guided tumor tracking)

60 APPLIED LIFE SCIENCES↗

Multi-Objective Adaptive Mesh Refinement Using Reinforcement Learning

Finite element methods approximate the solution to a partial differential equation (PDE) on a mesh consisting of many elements. In general, using more, smaller elements results in a lower error in the approximation. However, it is often possible to lower the error substantially by only refining, or decreasing the size of, the few elements in the mesh that have the highest error. Adaptive mesh refinement (AMR) is a process that selectively refines regions of a mesh with high error to achieve a desired accuracy in as few degrees of freedom (DOFs) as possible. AMR is favorable compared to uniform refinement, which refines all elements of the mesh equally, because it can often achieve the same accuracy without wasting extra computation time on refinement of elements that already have low error. However, it is difficult to know which elements to refine. In this report we explore ways to choose which elements to refine such that we minimize both the resulting error and the cumulative DOFs used in computation. In particular, we introduce a Pareto-front learning algorithm that trains a policy to give the optimal refinement actions to minimize the cumulative DOFs used to achieve a given target error. Such a policy is useful because it can be deployed on many different problem types where different accuracy levels are desired. Furthermore, training a single policy for a range of target errors allows us to use transfer learning to reduce the required training time.

97 MATHEMATICS AND COMPUTING↗