Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized simulation”

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 127 records · Page 7

SU-F-T-62: Three-Dimensional Dosimetric Gamma Analysis for Impacts of Tissue Inhomogeneity Using Monte Carlo Simulation in Intracavitary Brachytheray for Cervix Carcinoma

Purpose: The aim of this study was to investigate the impacts of tissue inhomogeneity on dose distributions using a three-dimensional (3D) gamma analysis in cervical intracavitary brachytherapy using Monte Carlo (MC) simulations. Methods: MC simulations for comparison of dose calculations were performed in a water phantom and a series of CT images of a cervical cancer patient (stage: Ib; age: 27) by employing a MC code, Particle and Heavy Ion Transport Code System (PHIT) version 2.73. The {sup 192}Ir source was set at fifteen dwell positions, according to clinical practice, in an applicator consisting of a tandem and two ovoids. Dosimetric comparisons were performed for the dose distributions in the water phantom and CT images by using gamma index image and gamma pass rate (%). The gamma index is the minimum Euclidean distance between two 3D spatial dose distributions of the water phantom and CT images in a same space. The gamma pass rates (%) indicate the percentage of agreement points, which mean that two dose distributions are similar, within an acceptance criteria (3 mm/3%). The volumes of physical and clinical interests for the gamma analysis were a whole calculated volume and a region larger than t% of a dose (close to a target), respectively. Results: The gamma pass rates were 77.1% for a whole calculated volume and 92.1% for a region within 1% dose region. The differences of 7.7% to 22.9 % between two dose distributions in the water phantom and CT images were found around the applicator region and near the target. Conclusion: This work revealed the large difference on the dose distributions near the target in the presence of the tissue inhomogeneity. Therefore, the tissue inhomogeneity should be corrected in the dose calculation for clinical treatment.

60 APPLIED LIFE SCIENCES↗

WE-AB-BRB-05: Toward a 2D Water-Equivalent Dosimetry Panel Using KCl:Eu2+

Purpose: KCl:Eu2+ storage phosphor shows promise for radiation therapy dosimetry. The purpose of this work is to investigate several important aspects of this material for potential commercial use. Methods: KCl:Eu2+ chips were fabricated and a conformal coating using Parylene was applied. Material’s dose response in a 6 MV beam was investigated using Monte-Carlo simulations. We attempted to micronize the materials using a spiral jet mill. As we did not have a water-free glovebox, we used commercially available non-hygroscopic BaFBr0.85I0.15:Eu2+ computed radiography material to test if a homogeneous panel can be made using micron-sized phosphors. Results: Dosimeters remained intact and showed no change in PSL intensity after eight hrs of submersion in water. We then optically bleached the samples for reuse, irradiated and immersed for another 24 hrs. We observed marginal worsening of the PSL signal for both the soaked and un-soaked chips. By contrast, we were unable to measure PSL intensity of the un-coated pellets due to these pellets dissolving within minutes of being immersed in water. MC data indicate that the micron-sized KCl:Eu2+ is predicted to have a nearly water-equivalent response. KCl:Eu2+ particles with a median size of 3 microns can be produced using a jet mill, which could be reduced further if necessary. While the particles tend to agglomerate over time when stored in a desiccator, they still possess favorable d50’s and d99’s even after 100 minutes, providing an adequate time window for making a panel via tape casting. A panel cast using optimized methods exhibits nearly perfect particle arrangement. Conclusions: Data shown here support ongoing efforts in fabricating a reusable, high resolution dosimetry panel in a water-free glovebox using micron-sized KCl:Eu2+ particles separated by water-equivalent polymers. The conformal coating thereafter will provide good humidity resistance. HL is the founder of DoseImaging, LLC that is exclusively dedicated to commercializing this technology.

60 APPLIED LIFE SCIENCES↗

SU-F-I-38: Patient Organ Specific Dose Assessment in Coronary CT Angiograph Using Voxellaized Volume Dose Index in Monte Carlo Simulation

Purpose: Clinical use of multi detector computed tomography(MDCT) in diagnosis of diseases due to high speed in data acquisition and high spatial resolution is significantly increased. Regarding to the high radiation dose in CT and necessity of patient specific radiation risk assessment, the adoption of new method in the calculation of organ dose is completely required and necessary. In this study by introducing a conversion factor, patient organ dose in thorax region based on CT image data using MC system was calculated. Methods: The geometry of x-ray tube, inherent filter, bow tie filter and collimator were designed using EGSnrc/BEAMnrc MC-system component modules according to GE-Light-speed 64-slices CT-scanner geometry. CT-scan image of patient thorax as a specific phantom was voxellised with 6.25mm3 in voxel and 64×64×20 matrix size. Dose to thorax organ include esophagus, lung, heart, breast, ribs, muscle, spine, spinal cord with imaging technical condition of prospectively-gated-coronary CT-Angiography(PGT) as a step and shoot method, were calculated. Irradiation of patient specific phantom was performed using a dedicated MC-code as DOSXYZnrc with PGT-irradiation model. The ratio of organ dose value calculated in MC-method to the volume CT dose index(CTDIvol) reported by CT-scanner machine according to PGT radiation technique has been introduced as conversion factor. Results: In PGT method, CTDIvol was 10.6mGy and Organ Dose/CTDIvol conversion factor for esophagus, lung, heart, breast, ribs, muscle, spine and spinal cord were obtained as; 0.96, 1.46, 1.2, 3.28. 6.68. 1.35, 3.41 and 0.93 respectively. Conclusion: The results showed while, underestimation of patient dose was found in dose calculation based on CTDIvol, also dose to breast is higher than the other studies. Therefore, the method in this study can be used to provide the actual patient organ dose in CT imaging based on CTDIvol in order to calculation of real effective dose(ED) based on organ dose. This work has been supported by the research chancellor of tehran university of medical sciences(tums), school of medicine, Tehran, Iran.

60 APPLIED LIFE SCIENCES↗

Optimization of an on-board imaging system for extremely rapid radiation therapy

Purpose: Next-generation extremely rapid radiation therapy systems could mitigate the need for motion management, improve patient comfort during the treatment, and increase patient throughput for cost effectiveness. Such systems require an on-board imaging system that is competitively priced, fast, and of sufficiently high quality to allow good registration between the image taken on the day of treatment and the image taken the day of treatment planning. In this study, three different detectors for a custom on-board CT system were investigated to select the best design for integration with an extremely rapid radiation therapy system. Methods: Three different CT detectors are proposed: low-resolution (all 4 × 4 mm pixels), medium-resolution (a combination of 4 × 4 mm pixels and 2 × 2 mm pixels), and high-resolution (all 1 × 1 mm pixels). An in-house program was used to generate projection images of a numerical anthropomorphic phantom and to reconstruct the projections into CT datasets, henceforth called “realistic” images. Scatter was calculated using a separate Monte Carlo simulation, and the model included an antiscatter grid and bowtie filter. Diagnostic-quality images of the phantom were generated to represent the patient scan at the time of treatment planning. Commercial deformable registration software was used to register the diagnostic-quality scan to images produced by the various on-board detector configurations. The deformation fields were compared against a “gold standard” deformation field generated by registering initial and deformed images of the numerical phantoms that were used to make the diagnostic and treatment-day images. Registrations of on-board imaging system data were judged by the amount their deformation fields differed from the corresponding gold standard deformation fields—the smaller the difference, the better the system. To evaluate the registrations, the pointwise distance between gold standard and realistic registration deformation fields was computed. Results: By most global metrics (e.g., mean, median, and maximum pointwise distance), the high-resolution detector had the best performance but the medium-resolution detector was comparable. For all medium- and high-resolution detector registrations, mean error between the realistic and gold standard deformation fields was less than 4 mm. By pointwise metrics (e.g., tracking a small lesion), the high- and medium-resolution detectors performed similarly. For these detectors, the smallest error between the realistic and gold standard registrations was 0.6 mm and the largest error was 3.6 mm. Conclusions: The medium-resolution CT detector was selected as the best for an extremely rapid radiation therapy system. In essentially all test cases, data from this detector produced a significantly better registration than data from the low-resolution detector and a comparable registration to data from the high-resolution detector. The medium-resolution detector provides an appropriate compromise between registration accuracy and system cost.

60 APPLIED LIFE SCIENCES↗

Computed Tomography of Dry Cask Storage for Used Fuel

The demand for dry cask storage for used nuclear fuel is on the rise due to the continuous operation of currently existing nuclear plants which are reaching or have reached the capacity of their used fuel pools. Dry cask storage allows used fuel that has already been cooled in the used fuel pool for at least five years to be surrounded by an inert gas inside a specially designed and licensed cask. Currently, 75 of the 78 U.S.A reactor sites either have or are pursuing an independent spent fuel storage installation (ISFSI) for dry cask fuel storage with over 2,000 casks already loaded. Cask loading will only continue to increase until a final repository solution is operational or the U.S. begins to recycle fuel. The challenge of monitoring the condition of the used fuel inside the cask without exposing the used fuel has been apparent since the 1980's. The objective of the University of Florida (UF) project is ultimately to investigate the feasibility of imaging the used fuel assemblies' location and evaluating the inner structural condition using emission source computed tomography (CT). A combination of radiation transport models will be used to simulate gamma-rays and neutrons traversing through the used fuel, structure, and cask body materials. In the simulations, multiple linear detector arrays will record the radiation outside the cask wall and will be placed at various angles in order to model rotating a linear detector array around the entire cask. The combination of data from all the detectors will give the planar radiation projection onto the detectors, more commonly called the Radon transform, of the structure of used fuel assemblies. This project can be processed with filtered back projection or other more advanced reconstructive imaging algorithms to recreate the fuel structures. A majority of the historical work done for used fuel has been on verification, diversion detection, and burnup credit validation in order to meet regulatory guidelines. One of the more commonly used systems for single assembly verification is the Fork Detector developed at Los Alamos National Laboratory in 1988. A group at the Idaho National Laboratory developed a Compton Dry-Cask Imaging Scanner that has the ability to detect empty vs full storage positions but it suffered a setback when transported to the Doel Nuclear Power Plant near Antwerp, Belgium; the system was unable to detect full energy peaks because of the additional scattering caused by a thick ballistic steel shield. As detectors, computational power, and modeling improves, a few groups have leveraged these advancements and started to investigate CT hoping, that despite the difficulties, the technique can give information on fuel condition. So far the groups have mostly limited the tomographic work to a per assembly basis , rather than the full dry cask as the University of Florida plans to investigate. The downsides of CT versus simpler verification and validation are the amount of data required and the time it takes to acquire that data. Given that dry casks are intended to sit undistributed for up to 60 years, a scanning time of several weeks to several months is not unreasonable if such systems do not require human operators beyond the initial setup. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Decontamination Experiments and Simulations Supporting Recovery After Widespread Contamination - 18685

Clean-up after an uncontrolled nuclear material release differs from common decontamination practices in the nuclear industry and military combat scenarios. Uncontrolled releases, particularly acts of terror, require specific protective actions guides (PAGs) to optimize short and long-term repercussions. Ideally, non-destructive decontamination methods would be employed. To help fulfill this need, the IWATERS (Integrated Wash Aid Treatment Emergency Reuse System) was proposed in 2014 for remediation outside the immediate blast area. 'Wash-Aid' is hosed onto affected areas, and the used wash solution is collected and recycled through sand/clay beds. To help characterize the IWATERS technology, the below was completed: - Decontamination experiments with a pressurized wash solution using Cs-137, Sr- 85, and Eu-152 contaminants on porous urban materials. - Sensitivity analyses on recycling bed variables using GoldSim{sup TM} CT. - A proposed framework that expands the recycling bed model to include probabilistic simulation of immediate recovery efforts. Removal efficacy of tap water and 0.5 M KCl were compared for Cs, Sr, and Eu across asphalt, brick, and concrete. The 0.5 M KCl solution significantly outperformed tap water in every case except for Eu on concrete and Cs on brick. Trends in the high-pressure decontamination with salt additive compared to low pressure decontamination show evidence of physical and chemical removal. Variations in cleaning rate yielded no significant trend. Pressure changes (by adjusting nozzle-to-coupon distance) were only significant for the highest pressure. Hypothetical Contamination Scenario: To tie the decontamination experiments and recycling bed model together, a small contamination scenario was created. Imagine the sides of a small concrete building uniformly contaminated with Cs-137. Power washing the building with 0.5 M KCl solution would result in 857 L of contaminated solution generated in 185 minutes. Two sand clay beds were designed to clean the used solution. With a sand mass fraction of 0.9, a 10 ft{sup 2} bed 4 ft tall could recycle the solution in 1.05 hrs. Changing the sand mass fraction to 0.5 yielded a 2.5 ft{sup 2} bed 4 ft tall (∼55-gallon drum) that recycled the solution in 7.43 hrs. This exercise showed how the sand/clay bed design can easily change to match the primary remediation goals. The second sand/clay bed would reduce the waste generation while the first bed reduces the remediation time. Conclusions: Variations on speed and nozzle distance (pressure) yielded insignificant differences in efficacy except for the closest distance. High-pressure washing showed increased decontamination efficacy with KCl additive in the majority of cases. Except for bed composition, key recycling bed parameters can be varied on-the-fly. This allows for optimization of, for example, time and waste generation. The hypothetical contamination scenario showed pairing decontamination experiment results with modeling tools can lead to increased preparedness for a large-scale contamination event. With collaboration and time, an expanded recovery/remediation model could aid in development of effective protective action guides.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

TU-E-201-03: Eye Lens Dosimetry in Radiotherapy Using Contact Lens-Shaped Applicator

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗

TU-E-201-02: Eye Lens Dosimetry From CT Perfusion Studies

Madan M. Rehani, Massachusetts General Hospital and Harvard Medical School, Boston Methods for Eye Lens Dosimetry and Studies On Lens Opacities with Interventionalists Radiation induced cataract is a major threat among staff working in interventional suites. Nearly 16 million interventional procedures are performed annually in USA. Recent studies by the principal investigator’s group, primarily among interventional cardiologists, on behalf of the International Atomic Energy Agency, show posterior subcapsular (PSC) changes in the eye lens in 38–53% of main operators and 21–45% of support staff. These changes have potential to lead to cataract in future years, as per information from A-Bomb survivors. The International Commission on Radiological Protection has reduced dose limit for staff by a factor of 7.5 (from 150 mSv/y to 20 mSv/y). With increasing emphasis on radiation induced cataracts and reduction in threshold dose for eye lens, there is a need to implement strategies for estimating eye lens dose. Unfortunately eye lens dosimetry is at infancy when it comes to routine application. Various approaches are being tried namely direct measurement using active or passive dosimeters kept close to eyes, retrospective estimations and lastly correlating patient dose in interventional procedures with staff eye dose. The talk will review all approaches available and ongoing active research in this area, as well as data from surveys done in Europe on status of eye dose monitoring in interventional radiology and nuclear medicine. The talk will provide update on how good is Hp(10) against Hp(3), estimations from CTDI values, Monte Carlo based simulations and current status of eye lens dosimetry in USA and Europe. The cataract risk among patients is in CT examinations of the head. Since radiation induced cataract predominantly occurs in posterior sub-capsular (PSC) region and is thus distinguishable from age or drug related cataracts and is also preventable, actions on awareness can lead to avoidance or even prevention. Learning Objectives: To understand recent changes in eye lens dose limits and thresholds for tissue reactions To understand different approaches to dose estimation for eye lens To learn about challenges in eye lens opacities among staff in interventional fluoroscopy Di Zhang, Toshiba America Medical Systems, Tustin, CA, USA Eye lens radiation dose from brain perfusion CT exams CT perfusion imaging requires repeatedly exposing one location of the head to monitor the uptake and washout of iodinated contrast. The accumulated radiation dose to the eye lens can be high, leading to concerns about potential radiation injury from these scans. CTDIvol assumes continuous z coverage and can overestimate eye lens dose in CT perfusion scans where the table do not increment. The radiation dose to the eye lens from clinical CT brain perfusion studies can be estimated using Monte Carlo simulation methods on voxelized patient models. MDCT scanners from four major manufacturers were simulated and the eye lens doses were estimated using the AAPM posted clinical protocols. They were also compared to CTDIvol values to evaluate the overestimation from CTDIvol. The efficacy of eye lens dose reduction techniques such as tilting the gantry and moving the scan location away from the eyelens were also investigated. Eye lens dose ranged from 81 mGy to 279 mGy, depending on the scanner and protocol used. It is between 59% and 63% of the CTDIvol values reported by the scanners. The eye lens dose is significantly reduced when the eye lenses were not directly irradiated. CTDIvol should not be interpreted as patient dose; this study has shown it to overestimate dose to the eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy. Learning Objectives: To become familiar with method of eye dose estimation for patient in specific situation of brain perfusion CT To become familiar with level of eye lens radiation doses in patients undergoing brain perfusion MDCT To understand methods for reducing eye lens dose to patient Jong Min Park, Department of Radiation Oncology, Seoul National University Hospital, Seoul, Korea Eye lens dosimetry in radiotherapy using contact lens-shaped applicator Dose calculation accuracy of commercial treatment planning systems is relatively low at shallow depths. Therefore, in-vivo measurements are often performed in the clinic to verify delivered doses to eye lens which are located at shallow depth. Current in-vivo dosimetry for eye lens during radiotherapy is generally performed with small in-vivo dosimeters on the surface of patient eyelid. Since this procedure potentially contains considerable uncertainty, a contact lens-shaped applicator made of acrylic (lens applicator) was developed for in-vivo measurements of eye lens dose during radiotherapy to reduce uncertainty. The lens applicator allows the insertion of commercially available metal oxide semiconductor field effect transistor (MOSFET) dosimeters. Computed tomography (CT) images of an anthropomorphic phantom with and without the lens applicator were acquired. A total of 20 VMAT plans were delivered to an anthropomorphic phantom and the doses with the lens applicator and the doses at the surface of the eyelid were measured using both micro and standard MOSFET dosimeters. The differences in measured dose at the surface of the eyelid from the calculated lens dose were acquired. The differences between the measured and the calculated doses at the lens applicator, as well as the differences between the measured and the calculated doses at the surface of the eyelid were acquired. The statistical significance of the differences was analyzed. The average difference between the measured and the calculated dose with the lens applicator was 16.8 % ± 10.4 % with a micro MOSFET dosimeter and 16.6 % ± 10.9% with a standard MOSFET dosimeter. The average difference without the lens applicator was 35.9% ± 41.5% with micro MOSFET dosimeter and 42.9% ± 52.2% with standard MOSFET dosimeter. The maximum difference with micro MOSFET dosimeter was 46% with the applicator and 188.4% without the applicator. For the standard MOSFET dosimeter, the maximum difference was 44.4% with the applicator and 246.4% without the applicator. The lens applicator allowed reduction of the differences between the calculated and the measured dose during in-vivo measurement for the eye lens as compared to in-vivo measurement at the surface of the eyelid. Learning Objectives: To understand limitations of dose calculation with commercial treatment planning system for eye lens during radiotherapy To learn about current in-vivo dosimetry methods for eye lens in the clinic To understand limitations of in-vivo dosimetry for eye lens during radiotherapy Di Zhang is an employee of Toshiba America Medical Systems.

60 APPLIED LIFE SCIENCES↗

How the “main condition” of phase stability can explain the effect of the velocity deviation of secondary electrons in DC-biased single-sided multipactors

In this work, a “main condition” for phase stability has been employed to investigate the effects of the velocity deviation of the electrons in DC-biased single-sided multipactors (MPs). In a previous study [M. Mostajeran, Phys. Plasmas 21, 053108 (2014)], a stability equation was derived, where the secondary electron was assumed to have zero initial velocity and the phase deviation from the resonant phase was considered. In this work, both deviations in phase and velocity from the resonant condition are taken into account, assuming nonzero initial velocity for the secondary electrons. Using the main condition for stability, it is shown that MP discharge can rise in situations, where large velocity deviations from initial velocity and large phase deviations from resonant phase exist. This is contrary to what can be predicted on the basis of the “simple stability condition.” This result is further confirmed by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of short electron bunches by a laser pulse crossing a sharp boundary of inhomogeneous plasma

The formation of short electron bunches during the passage of a laser pulse of relativistic intensity through a sharp boundary of semi-bounded plasma has been analytically studied. It is shown in one-dimensional geometry that one physical mechanism that is responsible for the generation of electron bunches is their self-injection into the wake field of a laser pulse, which occurs due to the mixing of electrons during the action of the laser pulse on plasma. Simple analytic relationships are obtained that can be used for estimating the length and charge of an electron bunch and the spread of electron energies in the bunch. The results of the analytical investigation are confirmed by data from numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stability and permanence of an eco-epidemiological SEIN model with impulsive biological control

Natural enemies of insects are extremely important in preventing pest outbreaks in crop fields. Therefore, we investigate the dynamical behavior of a pest-dependent consumption pest–natural enemy (i.e., prey–predator) SEIN model concerning diseases in pest population with three classes (susceptible–exposed–infectious) and impulsive releasing of infectious pests and natural enemies at fixed moments of time. We prove that all solutions of the system are uniformly ultimately bounded. In first part of the main results, the sufficient conditions for local as well as global asymptotic stability of the susceptible and exposed pest extinction periodic solution are determined using a Floquet’s theorem of impulsive differential equations, small-amplitude perturbation skills and comparison theorem. In second part, the sufficient condition for the permanence of a system is determined. These dynamics imply that susceptible and exposed pest populations become extinct when impulse period is less than some critical value and pests coexist with natural enemies at low level when impulse period crosses the critical value. Thus, our results provide some reliable theoretical tactics for pest management and finally these are verified by performing some numerical simulations.

97 MATHEMATICS AND COMPUTING↗

Prediction of plasma simulation data with the Gaussian process method

The simulation of plasma-wall interactions of fusion plasmas is extremely costly in computer power and time - the running time for a single parameter setting is easily in the order of weeks or months. We propose to exploit the already gathered results in order to predict the outcome for parametric studies within the high dimensional parameter space. For this we utilize Gaussian processes within the Bayesian framework and perform validation with one and two dimensional test cases from which we learn how to assess the outcome. Finally, the newly implemented method is applied to simulated data from the scrape-off layer of a fusion plasma. Uncertainties of the predictions are provided which point the way to parameter settings of further (expensive) simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetics of subdiffusive growth of new phase particles in supersaturated solid solutions

The kinetics of the subdiffusion-limited growth of spherical precipitates is studied. The process is described by the equation of anomalous diffusion with a fractional derivative with respect to time. It is shown that a decrease in the concentration of monomers is described by the law exp(−kt{sup 3α/2}) at the initial stage and the power law t{sup −α} at large times, where 0 < α ≤ 1 is the dispersion parameter coinciding with the order of time derivative in the subdiffusion equation. The time dependence of the size of a spherical precipitate is obtained. The results generalize the Ham diffusion theory and are in agreement with the Monte Carlo simulation data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Increasing ROI Using Your Digital Twin for Model Based Engineering

The purpose of this work is to identify innovative ways of using the plant reference simulator or digital twin differently than it has traditionally been used in the past. This paper will identify reasons that the plant reference simulator will support varying levels of model based design and engineering. I will point out real world examples in which the plant reference simulator was used to identify problems with engineering designs, digital control systems and plant modifications. The earlier these issues can be identified and resolved, the less impact they have on the plant modification or upgrade. By using the plant reference simulator earlier in the engineering process, utilities and engineering organizations can minimize costly outage extensions and increase the return on the investment made in the simulator.

42 ENGINEERING↗

On the Validity of Linear Response Theory in High-Dimensional Deterministic Dynamical Systems

This theoretical work considers the following conundrum: linear response theory is successfully used by scientists in numerous fields, but mathematicians have shown that typical low-dimensional dynamical systems violate the theory’s assumptions. Here we provide a proof of concept for the validity of linear response theory in high-dimensional deterministic systems for large-scale observables. We introduce an exemplary model in which observables of resolved degrees of freedom are weakly coupled to a large, inhomogeneous collection of unresolved chaotic degrees of freedom. By employing statistical limit laws we give conditions under which such systems obey linear response theory even if all the degrees of freedom individually violate linear response. We corroborate our result with numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dependence of Strain Rate Sensitivity on the Slip System: A Molecular Dynamics Simulation

The strain rate of deformation might affect mechanical properties of metallic materials, especially at elevated temperatures. Due to the nature of dislocation slip, it is anticipated that strain rate sensitivity (SRS) depends on the slip system. While the dependency of SRS on the temperature and strain rate of the deformation is well recognized, its dependence on the slip system is not well understood. Accordingly, the molecular dynamics simulation was used to investigate the dependence of strain rate sensitivity of pure Al single crystals on the slip system. In this study, the embedded atom method (EAM) potential for Al was employed. SRS and shear strength of the material were studied in four different slip systems and at two temperatures of 300 and 500 K. It was found that SRS of the material depends on the slip system in addition to the temperature, and SRS was higher in less compact systems with more difficult slip. The dislocation theories were used to rationalize the simulation results.

36 MATERIALS SCIENCE↗

Combination of the pair density approximation and the Takahashi–Imada approximation for path integral Monte Carlo simulations

We construct an accurate imaginary time propagator for path integral Monte Carlo simulations for heterogeneous systems consisting of a mixture of atoms and molecules. We combine the pair density approximation, which is highly accurate but feasible only for the isotropic interactions between atoms, with the Takahashi–Imada approximation for general interactions. We present finite temperature simulations results for energy and structure of molecules–helium clusters X{sup 4}He{sub 20} (X=HCCH and LiH) which show a marked improvement over the Trotter approximation which has a 2nd-order time step bias. We show that the 4th-order corrections of the Takahashi–Imada approximation can also be applied perturbatively to a 2nd-order simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rossby wave radiation by an eddy on a beta-plane: Experiments with laboratory altimetry

Results from the laboratory experiments on the evolution of vortices (eddies) generated in a rotating tank with topographic β-effect are presented. The focus of the experiments is on the far-field flow which contains Rossby waves emitted by travelling vortices. The surface elevation and velocity fields are measured by the altimetric imaging velocimetry. The experiments are supplemented by shallow water numerical simulations as well as a linear theory which describes the Rossby wave radiation by travelling vortices. The cyclonic vortices observed in the experiments travel to the northwest and continuously radiate Rossby waves. Measurements show that initially axisymmetric vortices develop a dipolar component which enables them to perform translational motion. A pattern of alternating zonal jets to the west of the vortex is created by Rossby waves with approximately zonal crests. Energy spectra of the flows in the wavenumber space indicate that a wavenumber similar to that introduced by Rhines for turbulent flows on the β-plane can be introduced here. The wavenumber is based on the translational speed of a vortex rather than on the root-mean-square velocity of a turbulent flow. The comparison between the experiments and numerical simulations demonstrates that evolving vortices also emit inertial waves. While these essentially three-dimensional non-hydrostatic waves can be observed in the altimetric data, they are not accounted for in the shallow water simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗