Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matlab data”

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 505 records · Page 28

CP decomposition for tensors via alternating least squares with QR decomposition

The CP tensor decomposition is used in applications such as machine learning and signal processing to discover latent low-rank structure in multidimensional data. Computing a CP decomposition via an alternating least squares (ALS) method reduces the problem to several linear least squares problems. The standard way to solve these linear least squares subproblems is to use the normal equations, which inherit special tensor structure that can be exploited for computational efficiency. However, the normal equations are sensitive to numerical ill-conditioning, which can compromise the results of the decomposition. In this paper, we develop versions of the CP-ALS algorithm using the QR decomposition and the singular value decomposition, which are more numerically stable than the normal equations, to solve the linear least squares problems. Our algorithms utilize the tensor structure of the CP-ALS subproblems efficiently, have the same complexity as the standard CP-ALS algorithm when the input is dense and the rank is small, and are shown via examples to produce more stable results when ill-conditioning is present. Our MATLAB implementation achieves the same running time as the standard algorithm for small ranks, and we show that the new methods can obtain lower approximation error.

97 MATHEMATICS AND COMPUTING↗

REBUS-PC v1.5 Verification Test Summary

Motivated in part by its usage for analysis of the University of Missouri Research Reactor (MURR), the verification of the REBUS-PC software for θ-R-Z geometries is discussed and displayed in this manuscript. Of particular importance is to show that the eigenvalue, flux, power, and burnup outputs of the software given a user problem is correct. REBUS-PC contains many features and MURR only makes use of some of those aspects. As an example, the use of region specific cross section sets, interpolation of the cross section data with respect to burnup, and automatic movement of control rods during the fuel cycle is discussed and displayed. To demonstrate the accuracy, 22 verification tests based upon 7 physical geometries are discussed. Geometry descriptions are provided in Section 2 and the REBUS-PC results for each verification test are discussed in Section 3. These verification tests consist of both analytic problems with known solutions and a more complicated input consistent with usage of REBUS-PC on the MURR reactor itself called STRZ. In this verification test summary report, the REBUS-PC software is repeatedly demonstrated to reproduce, or closely match, the eigenvalue, flux, and power distributions from analytic test problem solutions, demonstrating that the implemented finite difference approximation in θ-R-Z is working as intended. Several additional REBUS-PC features are tested, such as control rod movement by detailed comparisons of the results against reference solutions obtained using an auxiliary implemented MatLab program. As expected, the REBUS-PC software produces results that are 100% consistent with the stated capabilities and thus the verification exercise is completed. These verification tests should be executed with each re-compilation of the REBUS-PC source and after any source code changes. All output results from any such new work should be rigorously compared against the stored output of the 22 verification tests cited in this report that forms the basis of the verification process.

97 MATHEMATICS AND COMPUTING↗

Application of Crunch-Flow Routines to Constrain Present and Past Carbon Fluxes at Gas-Hydrate Bearing Sites

In November 2012, Oregon State University initiated the project entitled: Application of Crunch-Flow routines to constrain present and past carbon fluxes at gas-hydrate bearing sites. Within this project we developed Crunch-Flow based modeling modules that include important biogeochemical processes that need to be considered in gas hydrate environments. Our modules were applied to quantify carbon cycling in present and past systems, using data collected during several DOE-supported drilling expeditions, which include the Cascadia margin in US, Ulleung Basin in South Korea, and several sites drilled offshore India on the Bay of Bengal and Andaman Sea. Specifically, we completed modeling efforts that: 1) Reproduce the compositional and isotopic profiles observed at the eight drilled sites in the Ulleung Basin that constrain and contrast the carbon cycling pathways at chimney (high methane flux) and non-chimney sites (low methane, advective systems); 2) Simulate the Ba record in the sediments to quantify the past dynamics of methane flux in the southern Hydrate Ridge, Cascadia margin; and 3) Provide quantitative estimates of the thickness of individual mass transport deposits (MTDs), time elapsed after the MTD event, rate of sulfate reduction in the MTD, and time required to reach a new steady state at several sites drilled in the Krishna-Godavari (K-G) Basin off India. In addition we developed a hybrid model scheme by coupling a home-made MATLAB code with CrunchFlow to address the methane transport and chloride enrichment at the Ulleung Basins chimney sites, and contributed the modeling component to a study focusing on pore-scale controls on gas hydrate distribution in sediments from the Andaman Sea. These efforts resulted in two manuscripts currently under review, and contributed the modeling component of another pare, also under review. Lessons learned from these efforts are the basis of a mini-workshop to be held at Oregon State University (Feb 2014) to instruct graduate students (OSU and UW) as well as DOE staff from the NETL lab in Albany on the use of Crunch Flow for geochemical applications.

03 NATURAL GAS↗

Distributed Acoustic Sensing (DAS) of Strain at Earth Tide Frequencies: Laboratory Tests

The solid Earth strains in response to the gravitational pull from the Moon, Sun, and other planetary bodies. Measuring the flexure of geologic material in response to these Earth tides provides information about the geomechanical properties of rock and sediment. Such measurements are particularly useful for understanding dilation of faults and fractures in competent rock. A new approach to measuring earth tides using fiber optic distributed acoustic sensing (DAS) is presented here. DAS was originally designed to record acoustic vibration through the measurement of dynamic strain on a fiber optic cable. Here, laboratory experiments demonstrate that oscillating strain can be measured with DAS in the microHertz frequency range, corresponding to half-day (M2) lunar tidal cycles. Although the magnitude of strain measured in the laboratory is larger than what would be expected due to earth tides, a clear signal at half-day period was extracted from the data. With the increased signal-to-noise expected from quiet field applications and improvements to DAS using engineered fiber, earth tides could potentially be measured in deep boreholes with DAS. Because of the distributed nature of the sensor (0.25 m measurement interval over kilometers), fractures could be simultaneously located and evaluated. Such measurements would provide valuable information regarding the placement and stiffness of open fractures in bedrock. Characterization of bedrock fractures is an important goal for multiple subsurface operations such as petroleum extraction, geothermal energy recovery, and geologic carbon sequestration.

15 GEOTHERMAL ENERGY↗

SU-E-T-657: Quantitative Assessment of Plan Robustness for Helical Tomotherapy for Head and Neck Cancer Radiotherapy

Introduction: Geometric uncertainties in daily patient setup can lead to variations in the planned dose, especially when using highly conformal techniques such as helical Tomotherapy. To account for the potential effect of geometric uncertainty, our clinical practice is to expand critical structures by 3mm expansion into planning risk volumes (PRV). The PRV concept assumes the spatial dose cloud is insensitive to patient positioning. However, no tools currently exist to determine if a Tomotherapy plan is robust to the effects of daily setup variation. We objectively quantified the impact of geometric uncertainties on the 3D doses to critical normal tissues during helical Tomotherapy. Methods: Using a Matlab-based program created and validated by Accuray (Madison, WI), the planned Tomotherapy delivery sinogram recalculated dose on shifted CT datasets. Ten head and neck patients were selected for analysis. To simulate setup uncertainty, the patient anatomy was shifted ±3mm in the longitudinal, lateral and vertical axes. For each potential shift, the recalculated doses to various critical normal tissues were compared to the doses delivered to the PRV in the original plan Results: 18 shifted scenarios created from Tomotherapy plans for three patients with head and neck cancers were analyzed. For all simulated setup errors, the maximum doses to the brainstem, spinal cord, parotids and cochlea were no greater than 0.6Gy of the respective original PRV maximum. Despite 3mm setup shifts, the minimum dose delivered to 95% of the CTVs and PTVs were always within 0.4Gy of the original plan. Conclusions: For head and neck sites treated with Tomotherapy, the use of a 3mm PRV expansion provide a reasonable estimate of the dosimetric effects of 3mm setup uncertainties. Similarly, target coverage appears minimally effected by a 3mm setup uncertainty. Data from a larger number of patients will be presented. Future work will include other anatomical sites.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

TU-F-17A-01: BEST IN PHYSICS (JOINT IMAGING-THERAPY) - An Automatic Toolkit for Efficient and Robust Analysis of 4D Respiratory Motion

Purpose: To provide an automatic image analysis toolkit to process thoracic 4-dimensional computed tomography (4DCT) and extract patient-specific motion information to facilitate investigational or clinical use of 4DCT. Methods: We developed an automatic toolkit in MATLAB to overcome the extra workload from the time dimension in 4DCT. This toolkit employs image/signal processing, computer vision, and machine learning methods to visualize, segment, register, and characterize lung 4DCT automatically or interactively. A fully-automated 3D lung segmentation algorithm was designed and 4D lung segmentation was achieved in batch mode. Voxel counting was used to calculate volume variations of the torso, lung and its air component, and local volume changes at the diaphragm and chest wall to characterize breathing pattern. Segmented lung volumes in 12 patients are compared with those from a treatment planning system (TPS). Voxel conversion was introduced from CT# to other physical parameters, such as gravity-induced pressure, to create a secondary 4D image. A demon algorithm was applied in deformable image registration and motion trajectories were extracted automatically. Calculated motion parameters were plotted with various templates. Machine learning algorithms, such as Naive Bayes and random forests, were implemented to study respiratory motion. This toolkit is complementary to and will be integrated with the Computational Environment for Radiotherapy Research (CERR). Results: The automatic 4D image/data processing toolkit provides a platform for analysis of 4D images and datasets. It processes 4D data automatically in batch mode and provides interactive visual verification for manual adjustments. The discrepancy in lung volume calculation between this and the TPS is <±2% and the time saving is by 1–2 orders of magnitude. Conclusion: A framework of 4D toolkit has been developed to analyze thoracic 4DCT automatically or interactively, facilitating both investigational and clinical use. More software tools are under development for automatic analysis of other 4D imaging modalities. This research is in part supported by NIH (U54CA137788/132378) and City Seed FY13(Grant No. 93348-14 01). AY would like to thank MSKCC summer medical student research program supported by National Cancer Institute and hosted by Department of Medical Physics at MSKCC.

60 APPLIED LIFE SCIENCES↗

MO-F-CAMPUS-T-05: Correct Or Not to Correct for Rotational Patient Set-Up Errors in Stereotactic Radiosurgery

Purpose: To quantify the effect of patient rotation in stereotactic radiation therapy and establish a threshold where rotational patient set-up errors have a significant impact on target coverage. Methods: To simulate rotational patient set-up errors, a Matlab code was created to rotate the patient dose distribution around the treatment isocentre, located centrally in the lesion, while keeping the structure contours in the original locations on the CT and MRI. Rotations of 1°, 3°, and 5° for each of the pitch, roll, and yaw, as well as simultaneous rotations of 1°, 3°, and 5° around all three axes were applied to two types of brain lesions: brain metastasis and acoustic neuroma. In order to analyze multiple tumour shapes, these plans included small spherical (metastasis), elliptical (acoustic neuroma), and large irregular (metastasis) tumour structures. Dose-volume histograms and planning target volumes were compared between the planned patient positions and those with simulated rotational set-up errors. The RTOG conformity index for patient rotation was also investigated. Results: Examining the tumour volumes that received 80% of the prescription dose in the planned and rotated patient positions showed decreases in prescription dose coverage of up to 2.3%. Conformity indices for treatments with simulated rotational errors showed decreases of up to 3% compared to the original plan. For irregular lesions, degradation of 1% of the target coverage can be seen for rotations as low as 3°. Conclusions: This data shows that for elliptical or spherical targets, rotational patient set-up errors less than 3° around any or all axes do not have a significant impact on the dose delivered to the target volume or the conformity index of the plan. However the same rotational errors would have an impact on plans for irregular tumours.

60 APPLIED LIFE SCIENCES↗

Parametric and Sensitivity Analysis of a Steam Generator Model Using Python and Machine-Learning Tools

For this study, we used Python and machine-learning tools to perform a comprehensive parametric and sensitivity analysis on a steam generator (SG) model. (The Python model was based on a previously completed MATLAB framework for the Holtec SMR-160 SG.) We investigated the influence of various input parameters (e.g., heat transfer coefficient [HTC], Nusselt number, and heat exchanger effectiveness) on the system’s output. With machine-learning tools such as the Risk Analysis Virtual Environment (RAVEN), which was developed at Idaho National Laboratory, we were then able to perform an automated analysis of the SG inputs’ effect on the HTC. The analysis results give valuable insights into the performance and optimization of SG systems. We found the inlet mass flow rate (MFR) to have the greatest impact on the HTC, followed closely by the inlet temperature, and then pressure. Shifting of the input parameters causes the location of the maximum HTC along the SG length to change incrementally. The cold leg (CL) MFR was also found to impact the HTC magnitude as well as the location of the maximum HTC. At between 0.4–0.9 of the total SG length, the input parameters experience maximum impact on the HTC, leading us to suggest that sensors be efficiently placed on the SG so as to closely and effectively monitor thermal-hydraulic properties during reactor operation. We also found that the sensitivity data calculated manually agrees with the RAVEN – based data, confirming the same range of maximum sensitivity. However, the RAVEN-based analysis showed that cold leg pressure and hot leg temperature have a greater impact on the heat transfer coefficient than the mass flow rate, implying that a manual sensitivity study taking only two samples is not accurate.

20 FOSSIL-FUELED POWER PLANTS↗

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

SU-E-T-294: Simulations to Investigate the Feasibility of ‘dry’ Optical-CT Imaging for 3D Dosimetry

Purpose: To perform simulations investigating the feasibility of “dry” optical-CT, and determine optimal design and scanning parameters for a novel dry tank telecentric optical-CT 3D dosimetry system. Such a system would have important advantages in terms of practical convenience and reduced cost. Methods: A Matlab based ray-tracing simulation platform, ScanSim, was used to model a telecentric system with a polyurethane dry tank, cylindrical dosimeter, and surrounding fluid. This program's capabilities were expanded for the geometry and physics of dry scanning. To categorize the effects of refractive index (RI) mismatches, simulations were run for several dosimeter (RI = 1.5−1.48) and fluid (RI = 1.55−1.33) combinations. Additional simulations examined the effect of increasing gap size (1–5mm) between the dosimeter and tank wall, and of changing the telecentric lens tolerance (0.5°−5°). The evaluation metric is the usable radius; the distance from the dosimeter center where the measured and true doses differ by less than 2%. Results: As the tank/dosimeter RI mismatch increases from 0–0.02, the usable radius decreases from 97.6% to 50.2%. The fluid RI for matching is lower than either the tank or dosimeter RI. Changing gap sizes has drastic effects on the usable radius, requiring more closely matched fluid at large gap sizes. Increasing the telecentric tolerance through a range from 0.5°–5.0° improved the usable radius for every combination of media. Conclusion: Dry optical-CT with telecentric lenses is feasible when the dosimeter and tank RIs are closely matched (<0.01 difference), or when data in the periphery is not required. The ScanSim tool proved very useful in situations when the tank and dosimeter have slight differences in RI by enabling estimation of the optimal choice of RI of the small amount of fluid still required. Some spoiling of the telecentric beam and increasing the tolerance helps recover the usable radius.

07 ISOTOPE AND RADIATION SOURCES↗

A Comparison of Image Quality Evaluation Techniques for Transmission X-Ray Microscopy

Beamline 6-2c at Stanford Synchrotron Radiation Lightsource (SSRL) is capable of Transmission X-ray Microscopy (TXM) at 30 nm resolution. Raw images from the microscope must undergo extensive image processing before publication. Since typical data sets normally contain thousands of images, it is necessary to automate the image processing workflow as much as possible, particularly for the aligning and averaging of similar images. Currently we align images using the 'phase correlation' algorithm, which calculates the relative offset of two images by multiplying them in the frequency domain. For images containing high frequency noise, this algorithm will align noise with noise, resulting in a blurry average. To remedy this we multiply the images by a Gaussian function in the frequency domain, so that the algorithm ignores the high frequency noise while properly aligning the features of interest (FOI). The shape of the Gaussian is manually tuned by the user until the resulting average image is sharpest. To automatically optimize this process, it is necessary for the computer to evaluate the quality of the average image by quantifying its sharpness. In our research we explored two image sharpness metrics, the variance method and the frequency threshold method. The variance method uses the variance of the image as an indicator of sharpness while the frequency threshold method sums up the power in a specific frequency band. These metrics were tested on a variety of test images, containing both real and artificial noise. To apply these sharpness metrics, we designed and built a MATLAB graphical user interface (GUI) called 'Blur Master.' We found that it is possible for blurry images to have a large variance if they contain high amounts of noise. On the other hand, we found the frequency method to be quite reliable, although it is necessary to manually choose suitable limits for the frequency band. Further research must be performed to design an algorithm which automatically selects these parameters.

99 GENERAL AND MISCELLANEOUS↗

Boil-off losses along LH2 pathway

Losses along the LH 2 pathway are intrinsic to the utilization of a cryogenic fluid. They occur when the molecule is transferred between 2 vessels (liquefaction plant to trailer, trailer to station storage, station storage to pump or compressor, then fuel cell electric vehicles …) and when the fluid is warmed up due to heat transfer with the environment. Those losses can be estimated with good accuracy using thermodynamic models based on conservation of mass and energy, providing the thermodynamic states are correctly described. Indeed, the fluid undergoes various changes as it moves along the entire pathway (2 phase transition, super-heated warming, non-uniform temperature distributions across the saturation film) and accurate equations of state and 2 phase behavior implementations are essential. The balances of mass and energy during the various dynamics processes then enable to quantify the boil-off losses. In this work, a MATLAB code previously developed by NASA to simulate rocket loading is used as the basis for the LH 2 transfer model. This code implements complex physical phenomena such as the competition between condensation and evaporation and the convection vs. conduction heat transfer as a function of the relative temperatures on both sides of the saturated film. The original code was modified to consider real gas equations of state, and some semi-empirical relationships, such as between the heat of vaporization and the critical temperature, were also replaced by a REFPROP equivalent expression, assumed to be more accurate. Non-constant liquid temperature equations were added to simulate sub-cooled conditions. It is shown that although trailer depressurization and transfer losses from the receiving vessel during a LH 2 delivery may constitute the largest sources of losses, those can be greatly reduced and even eliminated if CFR/DOT regulations are properly applied and if no-vent fill methods are implemented, particularly for a stationary LH2 storage operating at reasonably low pressure (45 psia or so). As such, it is expected that the only remaining boil-off losses for a refueling station would come from the LH 2 pump (utilization and idling), pump vessel cool down/warm up, and environment heat transfer. Based on experimental data measured at LLNL on the Linde 875 bar LH 2 cryo-pump and results from the model, and extrapolating for refueling stations of various sizes, it can be shown that boil-off losses can vary from 15% of delivered LH 2 for a 100 kg/day, 5% at 400 kg/day, and down to less than 2% for stations above 1,800 kg/day (for a system similar to the one at LLNL). Less boil-off is to be expected for LH 2 pumps dispensing at 350 bar, so that less than 0.7% can be expected above 1,800 kg/day station capacities. At last, boil-off recovery solutions are briefly analyzed, including compressors, cryo-coolers and stationary fuel cell. Their economic value is strongly related to the costs of H2 and of the electricity needed to perform the recovery, on top of the initial capital expenditures. Effective boil-off could be reduced to 0.6-1% for a 2,000 kg/day station if adequate recovery solutions are used.

08 HYDROGEN↗

SU-E-T-194: From Dicom-RT to Radiobiological Dose Metrics in 5 Minutes

Purpose: To develop a flexible and standalone framework for batch calculation of radiobiological dose metrics from Dicom-RT. Methods: Software has been developed which allows (1) The calculation of DVH data from DICOM dose and structure files (DVHgenerator), (2) Calculation of a wide range of radiobiological metrics from this data (CompPlanGui). Both these tools are run via graphical user interface (GUI), making them fast and simple. Part 1 is a new tool which has not previously been published, whilst part 2 is a GUI overlay for the previously published software ‘Comp-Plan’ (Holloway et. al., Medical Dosimetry, 2012), previously reliant on command line interface. The time taken for an experienced user to evaluate a test case of 6 plans with and without CompPlanGUI was quantified. Results: The DVH-generator has been found to be faster, more robust and require far less physical memory then using alternative software solutions for the same purpose. The Comp Plan GUI significantly reduces the amount of time required to set up a base directory, eliminates code crashes arising from typographical errors, and renders the code far more accessible to non-expert users. It took an experienced user of the code around 3 minutes to set up a base directory of 6 plans compared around 8 minutes without, indicating that using CompPlanGUI reduced setup time by over 50%. Conclusion: A standalone GUI based framework has developed which allows for the batch calculation of radiobiological dose metrics directly from Dicom-RT files. As with the original code, this work will be made freely available on request, as well as via matlab file exchange.

07 ISOTOPE AND RADIATION SOURCES↗

SU-E-T-664: Radiobiological Modeling of Prophylactic Cranial Irradiation in Mice

Purpose: Prophylactic cranial irradiation (PCI) is a clinical technique used to reduce the incidence of brain metastasis and improve overall survival in select patients with ALL and SCLC, and we have shown the potential of PCI in select breast cancer patients through a mouse model (manuscript in preparation). We developed a computational model using our experimental results to demonstrate the advantage of treating brain micro-metastases early. Methods: MATLAB was used to develop the computational model of brain metastasis and PCI in mice. The number of metastases per mouse and the volume of metastases from four- and eight-week endpoints were fit to normal and log-normal distributions, respectively. Model input parameters were optimized so that model output would match the experimental number of metastases per mouse. A limiting dilution assay was performed to validate the model. The effect of radiation at different time points was computationally evaluated through the endpoints of incidence, number of metastases, and tumor burden. Results: The correlation between experimental number of metastases per mouse and the Gaussian fit was 87% and 66% at the two endpoints. The experimental volumes and the log-normal fit had correlations of 99% and 97%. In the optimized model, the correlation between number of metastases per mouse and the Gaussian fit was 96% and 98%. The log-normal volume fit and the model agree 100%. The model was validated by a limiting dilution assay, where the correlation was 100%. The model demonstrates that cells are very sensitive to radiation at early time points, and delaying treatment introduces a threshold dose at which point the incidence and number of metastases decline. Conclusion: We have developed a computational model of brain metastasis and PCI in mice that is highly correlated to our experimental data. The model shows that early treatment of subclinical disease is highly advantageous.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Simulink models for performance analysis of high speed DQPSK modulated optical link

This paper attempts to present the design approach for development of simulation models to study and analyze the transmission of 10 Gbps DQPSK signal over a single channel Peer to Peer link using Matlab Simulink. The simulation model considers the different optical components used in link design with their behavior represented initially by theoretical interpretation, including the transmitter topology, Mach Zehnder Modulator(MZM) module and, the propagation model for optical fibers etc. thus allowing scope for direct realization in experimental configurations. It provides the flexibility to incorporate the various photonic components as either user-defined or fixed and, can also be enhanced or removed from the model as per the design requirements. We describe the detailed operation and need of every component model and its representation in Simulink blocksets. Moreover the developed model can be extended in future to support Dense Wavelength Division Multiplexing (DWDM) system, thereby allowing high speed transmission with N × 40 Gbps systems. The various compensation techniques and their influence on system performance can be easily investigated by using such models.

42 ENGINEERING↗

Spatiotemporal measurements of striations in a glow discharge’s positive column using laser-collisional induced fluorescence

Here we have observed the behavior of striations caused by ionization waves propagating in low-pressure helium DC discharges using the non-invasive laser-collision induced fluorescence (LCIF) diagnostic. To achieve this, we developed an analytic fit of collisional radiative model (CRM) predictions to interpret the LCIF data and recover quantitative two-dimensional spatial maps of the electron density, n e , and the ratios of LCIF emission states that can be correlated with T e with the use of accurate distribution functions at localized positions within striated helium discharges at 500 mTorr, 750 mTorr, and 1 Torr. To our knowledge, these are the first spatiotemporal, laser-based, experimental measurements of n e in DC striations. The n e and 447:588 ratio distributions align closely with striation theory. Constriction of the positive column appears to occur with decreased gas pressure, as shown by the radial n e distribution. We identify a transition from a slow ionization wave to a fast ionization wave between 750 mTorr and 1 Torr. These experiments validate our analytic fit of n e , allowing the implementation of an LCIF diagnostic in helium without the need to develop a CRM.

42 ENGINEERING↗

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Solvation of the E. coli CheY Phosphorylation Site Mapped by XFMS

The Escherichia coli CheY protein belongs to a large bacterial response regulator superfamily. X-ray hydroxy radical foot-printing with mass spectroscopy (XFMS) has shown that allosteric activation of CheY by its motor target triggers a concerted internalization of aromatic sidechains. We reanalyzed the XFMS data to compare polar versus non-polar CheY residue positions. The polar residues around and including the 57D phosphorylated site had an elevated hydroxy radical reactivity. Bioinformatic measures revealed that a water-mediated hydrogen bond network connected this ring of residues with the central 57D. These residues solvated 57D to energetically stabilize the apo-CheY fold. The abundance of these reactive residues was reduced upon activation. This result was supported by the bioinformatics and consistent with the previously reported activation-induced increase in core hydrophobicity. It further illustrated XFMS detection of structural waters. Direct contacts between the ring residues and the phosphorylation site would stabilize the aspartyl phosphate. In addition, we report that the ring residue, 18R, is a constant central node in the 57D solvation network and that 18R non-polar substitutions determine CheY diversity as assessed by its evolutionary trace in bacteria with well-studied chemotaxis. These results showcase the importance of structured water dynamics for phosphorylation-mediated signal transduction.

59 BASIC BIOLOGICAL SCIENCES↗