Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distance Based Approximation”

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 379 records · Page 21

DESI peculiar velocity survey – Fundamental Plane

The Dark Energy Spectroscopic Instrument (DESI) peculiar velocity survey aims to measure the peculiar velocities of early- and late-type galaxies within the DESI footprint using both the Fundamental Plane and optical Tully–Fisher relations. Direct measurements of peculiar velocities can significantly improve constraints on the growth rate of structure, reducing uncertainty by a factor of approximately 2.5 at redshift 0.1 compared to the DESI Bright Galaxy Survey’s redshift space distortion measurements alone. We assess the quality of stellar velocity dispersion measurements from DESI spectroscopic data. These measurements, along with photometric data from the Legacy Survey, establish the Fundamental Plane relation and determine distances and peculiar velocities of early-type galaxies. During survey validation, we obtain spectra for 6698 unique early-type galaxies, up to a photometric redshift of 0.15. 64 per cent of observed galaxies (4267) have relative velocity dispersion errors below 10 per cent. This percentage increases to 75 per cent if we restrict our sample to galaxies with spectroscopic redshifts below 0.1. We use the measured central velocity dispersion, along with photometry from the DESI Legacy Imaging Surveys, to fit the Fundamental Plane parameters using a 3D Gaussian maximum likelihood algorithm that accounts for measurement uncertainties and selection cuts. In addition, we conduct zero-point calibration using the absolute distance measurements to the Coma cluster, leading to a value of the Hubble constant, H 0 = 76.05 ± 0.35 (statistical) ±0.49 (systematic Fundamental Plane) ±4.86 (statistical due to calibration) km s –1 Mpc –1 ⁠. This H 0 value is within 2σ of Planck cosmic microwave background results and within 1σ of other low-redshift distance indicator-based measurements.

cosmological parameters↗

Location of γ -ray emission and magnetic field strengths in OJ 287

We report the γ-ray BL Lac object OJ 287 is known to exhibit inner-parsec “jet-wobbling”, high degrees of variability at all wavelengths and quasi-stationary features, including an apparent (≈100°) position-angle change in projection on the sky plane. Sub-50 micro-arcsecond resolution 86 GHz observations with the global mm-VLBI array (GMVA) supplement ongoing multi-frequency VLBI blazar monitoring at lower frequencies. Using these maps, together with cm/mm total intensity and γ-ray observations from Fermi-LAT from 2008-2014, we aim to Observations with the GMVA offer approximately double the angular resolution compared with 43 GHz VLBA observations and enable us to observe above the synchrotron self-absorption peak frequency. Fermi-LAT γ-ray data were reduced and analysed. The jet was spectrally decomposed at multiple locations along the jet. From this, we could derive estimates of the magnetic field using equipartition and synchrotron self-absorption arguments. How the field decreases down the jet provided an estimate of the distance to the jet apex and an estimate of the magnetic field strength at the jet apex and in the broad line region. Combined with accurate kinematics, we attempt to locate the site of γ-ray activity, radio flares, and spectral changes. Strong γ-ray flares appeared to originate from either the so-called core region, a downstream stationary feature, or both, with γ-ray activity significantly correlated with radio flaring in the downstream quasi-stationary feature. Magnetic field estimates were determined at multiple locations along the jet, with the magnetic field found to be ≥1.6 G in the core and ≤0.4 G in the downstream quasi-stationary feature. Finally, we therefore found upper limits on the location of the VLBI core as ≲6.0 pc from the jet apex and determined an upper limit on the magnetic field near the jet base of the order of thousands of Gauss.

79 ASTRONOMY AND ASTROPHYSICS↗

Phenomenological Model of Hydrophobic and Hydrophilic Interactions

Hydration forces acting between macroscopic bodies at distances L ≤ 3 nm in pure water are calculated based on the phenomenological model of polar liquids. It is shown that depending on the properties of the bodies, the interacting surfaces polarize the liquid differently, and wetting properties of the surfaces are completely characterized by two parameters. If the surfaces are hydrophilic, liquid molecules are polarized at right angles to the surfaces, and the interaction is the short-range repulsion (the forces of interaction decrease exponentially over the characteristic length λ ≈ 0.2 nm). The interaction between the hydrophobic surfaces is more diversified and has been studied less. For L ≤ 3 nm, the interaction exhibits universal properties, while for L ≤ 3 nm, it considerably depends on the properties of the surfaces and on the distances between them, as well as on the composition of the polar liquid. In full agreement with the available experimental results we find that if the interfaces are mostly hydrophobic, then the interaction is attractive and long-range (the interaction forces diminish exponentially with decay length 1.2 nm). In this case, the resultant polarization of water molecules is parallel to the surface. It is shown that hydration forces are determined by nonlinear effects of polarization of the liquid in the bulk or by analogous nonlinearity of the interaction of water with a submerged body. This means that the forces of interaction cannot be calculated correctly in the linear response approximation. The forces acting between hydrophobic or hydrophilic surfaces are of the entropy type or electrostatic, respectively. It is shown that hydrophobic and hydrophilic surfaces for L ≤ 3 nm repel each other. The calculated intensity of their interaction is in agreement with experimental data. We predict the existence of an intermediate regime in which a body cannot order liquid molecules, which results in a much weaker attraction that decreases in inverse proportion to the squared distance between the surfaces of bodies. The difference between the microscopic structures of liquids confined in nanovolumes from liquids in large volumes is considered. The proposed model is applicable for a quantitative description of the properties of water at temperatures T satisfying the condition 0 < (T–T{sub c})/T{sub c} ≪ 1, where T{sub c} ≈ 230 K is the temperature of the ferroelectric phase transition observed in supercooled water. Under standard conditions, the model can be used for obtaining qualitative estimates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Distance-to-Agreement Investigation of Tomotherapy's Bony Anatomy-Based Autoregistration and Planning Target Volume Contour-Based Optimization

Purpose: To compare Tomotherapy's megavoltage computed tomography bony anatomy autoregistration with the best achievable registration, assuming no deformation and perfect knowledge of planning target volume (PTV) location. Methods and Materials: Distance-to-agreement (DTA) of the PTV was determined by applying a rigid-body shift to the PTV region of interest of the prostate from its reference position, assuming no deformations. Planning target volume region of interest of the prostate was extracted from the patient archives. The reference position was set by the 6 degrees of freedom (dof)—x, y, z, roll, pitch, and yaw—optimization results from the previous study at this institution. The DTA and the compensating parameters were calculated by the shift of the PTV from the reference 6-dof to the 4-dof—x, y, z, and roll—optimization. In this study, the effectiveness of Tomotherapy's 4-dof bony anatomy–based autoregistration was compared with the idealized 4-dof PTV contour-based optimization. Results: The maximum DTA (maxDTA) of the bony anatomy-based autoregistration was 3.2 ± 1.9 mm, with the maximum value of 8.0 mm. The maxDTA of the contour-based optimization was 1.8 ± 1.3 mm, with the maximum value of 5.7 mm. Comparison of Pearson correlation of the compensating parameters between the 2 4-dof optimization algorithms shows that there is a small but statistically significant correlation in y and z (0.236 and 0.300, respectively), whereas there is very weak correlation in x and roll (0.062 and 0.025, respectively). Conclusions: We find that there is an average improvement of approximately 1 mm in terms of maxDTA on the PTV going from 4-dof bony anatomy-based autoregistration to the 4-dof contour-based optimization. Pearson correlation analysis of the 2 4-dof optimizations suggests that uncertainties due to deformation and inadequate resolution account for much of the compensating parameters, but pitch variation also makes a statistically significant contribution.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Simulating the Effect of Fracture Connectivity on Repository Performance with GDSA Framework - 18589

This work examines the possible migration of radionuclides from a deep geologic nuclear waste repository sited in fractured crystalline host rock. The key safety concern is the potential effect on waste isolation of the inter-connectivity of the fracture network, which is primarily established by the temporal evolution of the temperature and stress fields at the time of original rock intrusion and cooling. Two end members are considered in these simulations, one with a high degree of connectivity to the biosphere, such that advective transport through the fracture and fault network controls radionuclide migration, and the other with a low degree of connectivity, such that slow diffusive transport through the crystalline rock matrix is the controlling process for radionuclide migration to the biosphere. Both end members have several areally extensive, high-transmissivity deformation zones, whose distance to the repository is an important factor for waste isolation capability. Uncertainties in fracture properties (transmissivity, orientation, radius) and fracture distribution also give rise to uncertainty in waste isolation capability (i.e., the overall connectivity of the repository with the surface biosphere). These types of natural system uncertainties are generally present in any geologic site-characterization program and thus represent an important factor in assessing repository performance in hard rock environments. The waste considered in these simulations are spent fuel rod assemblies from the U.S. commercial nuclear reactor fleet. Thermal output of the waste must be considered in the simulations and can influence the rate and timing of waste package failure, waste form degradation, and fluid flux due to thermal expansion around the repository horizon. To achieve a high-fidelity representation of radionuclide transport in fractures and rock matrix, combined with thermal energy transport and fluid flow in fractures and matrix, the mathematical model is solved numerically in a parallel high-performance computing (HPC) environment on a finite volume unstructured grid consisting of approximately 5 million cells, using Geologic Disposal Safety Assessment (GDSA) Framework (https://pa.sandia.gov), an open-source performance assessment tool for deep underground disposal of nuclear waste. GDSA Framework uses PFLOTRAN to solve the balance equations on a three-dimensional grid with heterogeneous properties, using multiple processors in a parallel configuration based on domain decomposition. The fracture networks in these simulations are originally generated as discrete fracture networks (DFNs), which are sets of two-dimensional planes distributed in a three-dimensional domain. The method used in GDSA Framework maps the stochastically generated DFN to an equivalent continuous porous medium (ECPM) domain that allows for the simulation of coupled heat flow, fluid flow, and radionuclide transport, including heat conduction through the matrix of the fractured rock, which is a process not easily modeled in a DFN representation. Computational efficiency is also greatly enhanced using the ECPM method, allowing for a realistic representation and analysis of uncertainties in a multi-realization performance assessment of a deep geologic repository. The effect of fracture connectivity on the waste isolation safety function, as brought to light by these GDSA Framework simulations, points to the importance of including a realistic representation of uncertainties in fracture properties and distribution (effectively, uncertainty in spatial heterogeneity) in repository safety assessment simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Third-epoch Magellanic Cloud proper motions. II. The large Magellanic Cloud rotation field in three dimensions

We present the first detailed assessment of the large-scale rotation of any galaxy based on full three-dimensional velocity measurements. We do this for the LMC by combining our Hubble Space Telescope average proper motion (PM) measurements for stars in 22 fields, with existing line-of-sight (LOS) velocity measurements for 6790 individual stars. We interpret these data with a model of circular rotation in a flat disk. The PM and LOS data paint a consistent picture of the LMC rotation, and their combination yields several new insights. The PM data imply a stellar dynamical center that coincides with the H I dynamical center, and a rotation curve amplitude consistent with that inferred from LOS velocity studies. The implied disk viewing angles agree with the range of values found in the literature, but continue to indicate variations with stellar population and/or radius. Young (red supergiant) stars rotate faster than old (red and asymptotic giant branch) stars due to asymmetric drift. Outside the central region, the circular velocity is approximately flat at V {sub circ} = 91.7 ± 18.8 km s{sup –1}. This is consistent with the baryonic Tully-Fisher relation and implies an enclosed mass M(8.7 kpc) = (1.7 ± 0.7) × 10{sup 10} M {sub ☉}. The virial mass is larger, depending on the full extent of the LMC's dark halo. The tidal radius is 22.3 ± 5.2 kpc (24.°0 ± 5.°6). Combination of the PM and LOS data yields kinematic distance estimates for the LMC, but these are not yet competitive with other methods.

79 ASTRONOMY AND ASTROPHYSICS↗

SPiRaL: a multiresolution global tomography model of seismic wave speeds and radial anisotropy variations in the crust and mantle

SUMMARY SPiRaL is a joint global-scale model of wave speeds (P and S) and anisotropy (vertical transverse isotropy, VTI) variations in the crust and mantle. The model is comprised of >2.1 million nodes with five parameters at each node that capture velocity variations for P- and S-waves travelling at arbitrary directions in transversely isotropic media with a vertical symmetry axis (VTI). The crust (including ice, water, sediments and crystalline layers) is directly incorporated into the model. The default node spacing is approximately 2° in the lower mantle and 1° in the crust and upper mantle. The grid is refined with ∼0.25° minimum node spacing in highly sampled regions of the crust and upper mantle throughout North America and Eurasia. The data considered in the construction of SPiRaL includes millions of body wave traveltimes (crustal, regional and teleseismic phases with multiples) and surface wave (Rayleigh and Love) dispersion. A multiresolution inversion approach is employed to capture long-wavelength heterogeneities commonly depicted in global-scale tomography images as well as more localized details that are typically resolved in more focused regional-scale studies. Our previous work has demonstrated that such global-scale models with regional-scale detail can accurately predict both teleseismic and regional body wave traveltimes, which is necessary for more accurate location of small seismic events that may have limited signal at teleseismic distances. SPiRaL was constructed to predict traveltimes for event location and long-period waveform dispersion for seismic source inversion applications in regions without sufficiently tuned models. SPiRaL may also serve as a starting model for full-waveform inversion (FWI) with the goal of fitting waves with periods 10–50 s over multiple broad regions (thousands of kilometres) and potentially the globe. To gain insight to this possibility, we simulated waveforms for a small set of events using SPiRaL and independent waveform-based models for comparison. For the events tested, the performance of the traveltime-based SPiRaL model is shown to be generally on par with regional 3-D waveform-based models in three regions (western United States, Middle East, Korean Peninsula) suggesting SPiRaL may serve as a starting model for FWI over broad regions.

58 GEOSCIENCES↗

Hydraulic Response to Thermal Stimulation Efforts at Raft River Based on Stepped Rate Injection Testing

The injection well stimulation project at the Raft River geothermal field tests the effect of long-term cold water injection and high pressure injection on well injectivity, improvements to which could reduce operating costs. The primary data for analysis and interpretation of the injection test are step-rate flow tests run before each new phase of the injection. These tests were analyzed using a combination of standard pump-test analytical solution methods and methods developed expressly for the observed conditions. The stepped rate injection tests, combined with long-term flow and pressure response data suggest that the well is located within a fractured formation of low transmissivity but high storativity. These calculated parameters appeared to increase with pressure during the first injection test and the higher values were reproduced during the second stepped rate test. Calculated transmissivity and storativity are on the order of 4E-5 m cm 2 and 1E-4 m Pa -1 , respectively. The apparent pressure dependence of fitted hydraulic parameters may reflect near-well fracture compliance that increased the effective radius of the wellbore during the first test. While the type curve fit analysis also suggests that the reservoir behaves as a uniformly fractured reservoir with a radial flow regime, the hydraulic parameters indicate that condition may exist only a very limited distance (<10 m) from the well. Longer-term pressure response suggests that flow in the system effectively reaches steady state in a period of less than a day, which may reflect pressure stabilization resulting from pressure-dependent permeability or a region of much higher permeability located with a few meters of the well. The transmissivity estimates obtained from this analysis, converted to approximate fracture density and aperture, provide useful constraints on the distance to which the thermal front may migrate from the well during the cold water injection phase of the stimulation project. We estimate that the cooling front will migrate less than a tenth of a kilometer over an approximately one-year injection period. Here, the effects of that cooling, however, may be substantial, because increases in permeability have maximum effect nearest the well.

cold water injection↗

Tri-Sectional Approximation of the Shortest Path to Long-Term Voltage Stability Boundary with Distributed Energy Resources

Ensuring long-term voltage stability is critical for reliable operations of power grids. High share of distributed energy resources (DERs) can create complicated system operation modes that may invalidate the traditional long-term voltage stability analysis based on typical operation modes. To address this challenge, this paper investigates how to compute the shortest path to the voltage stability boundary in the DER aggregated load space with large dispersion. Instead of working in the Euclidean space, we establish the analysis and computations on the algebraic power flow manifold to better capture the curvature change of the shortest path along the direction of losing stability. A modified optimal control framework is presented for obtaining the ground-truth of the smooth shortest path on the manifold. To efficiently and accurately solve for the shortest path, we further leverage the geometric features of the power flow manifold and propose a tri-sectional approximation model that is scalable for large-scale systems. Several numerical examples, up to the 1354-bus system, with different DER penetration levels and high dimensional renewable power injection variations are evaluated. The simulation results demonstrate that the tri-sectional approximation achieves high accuracy and efficiency to approximate the shortest path to the voltage stability boundary.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A complete software application for automatic registration of x-ray mammography and magnetic resonance images

Purpose: This work presents a complete and automatic software application to aid radiologists in breast cancer diagnosis. The application is a fully automated method that performs a complete registration of magnetic resonance (MR) images and x-ray (XR) images in both directions (from MR to XR and from XR to MR) and for both x-ray mammograms, craniocaudal (CC), and mediolateral oblique (MLO). This new approximation allows radiologists to mark points in the MR images and, without any manual intervention, it provides their corresponding points in both types of XR mammograms and vice versa. Methods: The application automatically segments magnetic resonance images and x-ray images using the C-Means method and the Otsu method, respectively. It compresses the magnetic resonance images in both directions, CC and MLO, using a biomechanical model of the breast that distinguishes the specific biomechanical behavior of each one of its three tissues (skin, fat, and glandular tissue) separately. It makes a projection of both compressions and registers them with the original XR images using affine transformations and nonrigid registration methods. Results: The application has been validated by two expert radiologists. This was carried out through a quantitative validation on 14 data sets in which the Euclidean distance between points marked by the radiologists and the corresponding points obtained by the application were measured. The results showed a mean error of 4.2 ± 1.9 mm for the MRI to CC registration, 4.8 ± 1.3 mm for the MRI to MLO registration, and 4.1 ± 1.3 mm for the CC and MLO to MRI registration. Conclusions: A complete software application that automatically registers XR and MR images of the breast has been implemented. The application permits radiologists to estimate the position of a lesion that is suspected of being a tumor in an imaging modality based on its position in another different modality with a clinically acceptable error. The results show that the application can accelerate the mammographic screening process for high risk populations or for dense breasts.

60 APPLIED LIFE SCIENCES↗

Dose Rate Analysis of the WCS Consolidated Interim Storage Facility

This report describes a confirmatory calculation for a proposed consolidated interim storage facility (CISF) to support the US Nuclear Regulatory Commission’s review of the license application submitted by the Interim Storage Partners. The scope of this report is limited to an assessment of the annual dose from direct radiation associated with Phase 1 of the Waste Control Specialists (WCS) CISF and the determination of the minimum controlled area boundary based on the 10 CFR 72.104 limit of 25 mrem annual dose to the whole body. The confirmatory dose rate calculations used the source term and shielding calculation capabilities of the SCALE 6.2.3 computer code system. The calculation method employed in the confirmatory calculations uses a detailed Monte Carlo radiation transport simulation from source to dose rate. This method differs from the two-step method used in the safety analysis report (SAR), which requires determination of the photon and neutron energy and angular distributions on the cask external surface and then the use of this surface source in a new radiation transport calculation to determine dose rate as a function of distance from the storage facility. A simulation was made using the complete site geometry (all casks present) but with only one cask containing source, and a total of 467 independent calculations were performed to obtain the dose rate maps produced by each storage cask. Basic input data and assumptions used in the SAR (e.g., cask design parameters and design basis assembly characteristics) were used in this confirmatory calculation for consistency with the SAR. Local atmospheric conditions for skyshine calculations were simulated using US standard atmosphere data, which describe average air temperature, pressure, and density as a function of altitude. The determined minimum required distances from the WCS Phase 1 CISF to a contour defined by the 25 mrem annual dose limit were approximately 623 m (0.387 miles) in the N-NE direction, 594 m (0.369 miles) in the SSW direction, and 533 m (0.331 miles) in the E-SE direction and the W-NW direction. The annual dose results of the confirmatory calculation are in relatively good agreement with the annual dose values presented in the WCS SAR. The annual dose produced by the Phase 1 CISF at the nearest site boundary, 0.75 miles away from the CISF, was estimated by this study to be approximately 6.02E-02 ± 2.0E-04 mrem, which is significantly smaller than the 25 mrem annual dose limit provided in 10 CFR 72.104. The SAR annual dose rate value for the nearest site boundary was 7.52E-2 mrem, which included 7.73E-03 mrem due to postulated leakage of the FO-, FC-, and FF-canisters. A comparison with the dose rate values reported in the SAR for 16 detector locations, D1 through D16, on the protected area boundary is also included in this report. This comparison shows that the SAR determined higher annual dose values for detectors D1 through D14 by 12% to 45% and lower annual dose values for detectors D15 and D16 by 18% and 8%, respectively, as compared with this ORNL confirmatory calculation. Sensitivity studies have been performed to support a comprehensive review of the application. In these studies, the effects on dose rate of the modeled air volume, density, and humidity were analyzed. The sensitivity calculations showed that the air volume used in the Monte Carlo simulations was adequate and that a larger air volume would not further increase the dose rate due to skyshine contributions. The sensitivity analysis showed that simulation of local atmospheric conditions is important for the accuracy of dose rate estimates at large distances from the facility. Air humidity was demonstrated to have negligible effects on dose rate values for the simulated local atmospheric conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗

Hydrodynamic limit of Wigner-Poisson kinetic theory: Revisited

In this paper, we revisit the hydrodynamic limit of the Langmuir wave dispersion relation based on the Wigner-Poisson model in connection with that obtained directly from the original Lindhard dielectric function based on the random-phase-approximation. It is observed that the (fourth-order) expansion of the exact Lindhard dielectric constant correctly reduces to the hydrodynamic dispersion relation with an additional term of fourth-order, beside that caused by the quantum diffraction effect. It is also revealed that the generalized Lindhard dielectric theory accounts for the recently discovered Shukla-Eliasson attractive potential (SEAP). However, the expansion of the exact Lindhard static dielectric function leads to a k{sup 4} term of different magnitude than that obtained from the linearized quantum hydrodynamics model. It is shown that a correction factor of 1/9 should be included in the term arising from the quantum Bohm potential of the momentum balance equation in fluid model in order for a correct plasma dielectric response treatment. Finally, it is observed that the long-range oscillatory screening potential (Friedel oscillations) of type cos(2k{sub F}r)/r{sup 3}, which is a consequence of the divergence of the dielectric function at point k = 2k{sub F} in a quantum plasma, arises due to the finiteness of the Fermi-wavenumber and is smeared out in the limit of very high electron number-densities, typical of white dwarfs and neutron stars. In the very low electron number-density regime, typical of semiconductors and metals, where the Friedel oscillation wavelength becomes much larger compared to the interparticle distances, the SEAP appears with a much deeper potential valley. It is remarked that the fourth-order approximate Lindhard dielectric constant approaches that of the linearized quantum hydrodynamic in the limit if very high electron number-density. By evaluation of the imaginary part of the Lindhard dielectric function, it is shown that the Landau-damping region in ω-k plane increases dramatically by increase of the electron number-density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of an Analytical Model for Therapeutic and Stray Dose Calculations of 6 MV Photon Beams

Cancer is one of the most significant health care problems today. It is estimated that 1 in 2 men and women will be diagnosed with some form of cancer in their lifetime. Nearly two thirds of cancer patients will receive radiation therapy at some point in their care and the vast majority of these radiotherapy treatments will be performed with megavoltage photon therapy. Photon treatment techniques typically devote considerable attention to limiting the exposure of healthy tissues outside of the intended treatment volume. Nonetheless, small levels of unwanted stray radiation are incident upon the patient's whole body, which can increase the risk for numerous late effects, including second malignancies. When one considers that second cancers account for approximately 15% of all cancers and that it is projected that by the year 2020, there will be almost 20 million cancer survivors in the United States and four million in the UK, the need for improved accuracy in stray dose calculations appears likely to increase dramatically. Despite the obvious benefits this capability would afford, contemporary commercial treatment planning systems (TPSs) lack the ability to accurately predict stray dose outside the treatment field, e.g., beyond the 5% isodose surface. Thus, there is a large, variable, and systematic underestimation of stray radiation exposures by commonly used clinical TPSs for photon therapy. To overcome this problem, in a companion paper we proposed a simple, physics-based analytical model that calculates the total therapeutic and stray dose from photon treatment fields. The purpose of this work was to evaluate the accuracy of the basic physics in the analytical absorbed dose model for several simple treatment fields. We assessed the performance of the basic physics in an analytical model that calculates the total absorbed dose deposited by primary, leakage, and scattered photons both inside and outside of square treatment fields of 6 MV photon conventional radiotherapy. The assessment was accomplished by comparing model predictions of total dose with measured dose profiles in a water-phantom. The major finding of this study is that the physics-based model performed well, predicting dose in water within 11.5% dose difference or 5 mm distance to agreement with measured values in the therapeutic region and within 1.1 mGy/Gy or 0.5 mm distance to agreement with measured values in the out-of-field region at 90% of the locations considered. Contemporary external-beam photon TPSs do not include the capability to accurately predict the absorbed dose beyond a few centimeters outside of the treatment field. Our work suggests that it is feasible to improve dosimetric accuracy in the out-of-field region using a simple and computationally fast, physics-based analytical algorithm. This is demonstrated by the measured absorbed dose as well as the corresponding values calculated by two contemporary commercial TPSs (Pinnacle, Philips Health care, Andover, MA and Eclipse, Varian Medical Systems, Inc., Palo Alto, CA) and the analytical model considered in this work. It is apparent that both TPSs underestimated the absorbed dose by greater than 3 mGy/Gy at locations outside of the 0.5% relative dose level. In contrast, the analytical model evaluated in this work calculated absorbed dose in the same region within 0.4 mGy/Gy of measured values. The physics-based analytical dose model could prove useful in a number of different applications. The accurate knowledge of absorbed dose deposited outside of the treatment field is especially important for predicting radiation late effects, such as second cancer, fertility complications, and cardiac toxicity. It may also facilitate predictions of the effects of in utero exposure. We conclude that it is possible to accurately calculate the therapeutic and stray absorbed dose from a conventional radiotherapy treatment field of 6 MV photons using a physics-based analytical model. This was demonstrated by the agreement between measured and calculated doses in water within 11.5% relative dose difference or 5 mm distance to agreement in the therapeutic region, and within 1.1 mGy/Gy absolute dose difference or 0.5 mm distance to agreement in the out-of-field region at 90% of positions for all of the fields considered in this work. This is also illustrated by Figure 4 which plots dose values from analytical model predictions and measurements at five representative irradiation conditions and shows that the analytical model agreement is good in all cases. Furthermore, the model validated in this work evidently can predict dose values with better accuracy at out-of-field locations than contemporary TPSs. This and previous works, taken together, suggest that analytical physics-based models will offer advantages over TPSs or Monte Carlo models for whole-body dose calculations for several radiotherapy applications.

61 RADIATION PROTECTION AND DOSIMETRY↗

Final Report. Analysis and Reduction of Complex Networks Under Uncertainty

The project was a collaborative effort among MIT, Sandia National Laboratories (local PI Dr. Habib Najm), the University of Southern California (local PI Prof. Roger Ghanem), and The Johns Hopkins University (local PI Prof. Omar Knio, now at Duke University). Our focus was the analysis and reduction of large-scale dynamical systems emerging from networks of interacting components. Such networks underlie myriad natural and engineered systems. Examples important to DOE include chemical models of energy conversion processes, and elements of national infrastructure—e.g., electric power grids. Time scales in chemical systems span orders of magnitude, while infrastructure networks feature both local and long-distance connectivity, with associated clusters of time scales. These systems also blend continuous and discrete behavior; examples include saturation phenomena in surface chemistry and catalysis, and switching in electrical networks. Reducing size and stiffness is essential to tractable and predictive simulation of these systems. Computational singular perturbation (CSP) has been effectively used to identify and decouple dynamics at disparate time scales in chemical systems, allowing reduction of model complexity and stiffness. In realistic settings, however, model reduction must contend with uncertainties, which are often greatest in large-scale systems most in need of reduction. Uncertainty is not limited to parameters; one must also address structural uncertainties—e.g., whether a link is present in a network—and the impact of random perturbations, e.g., fluctuating loads or sources. Research under this project developed new methods for the analysis and reduction of complex multiscale networks under uncertainty, by combining computational singular perturbation (CSP) with probabilistic uncertainty quantification. CSP yields asymptotic approximations of reduceddimensionality “slow manifolds” on which a multiscale dynamical system evolves. Introducing uncertainty in this context raised fundamentally new issues, e.g., how is the topology of slow manifolds transformed by parametric uncertainty? How to construct dynamical models on these uncertain manifolds? To address these questions, we used stochastic spectral polynomial chaos (PC) methods to reformulate uncertain network models and analyzed them using CSP in probabilistic terms. Finding uncertain manifolds involved the solution of stochastic eigenvalue problems, facilitated by projection onto PC bases. These problems motivated us to explore the spectral properties stochastic Galerkin systems. We also introduced novel methods for rank-reduction in stochastic eigensystems—transformations of a uncertain dynamical system that lead to lower storage and solution complexity. These technical accomplishments are detailed below. This report focuses on the MIT portion of the joint project.

97 MATHEMATICS AND COMPUTING↗

THE PROPER MOTION OF PALOMAR 5

Palomar 5 (Pal 5) is a faint halo globular cluster associated with narrow tidal tails. It is a useful system to understand the process of tidal dissolution, as well as to constrain the potential of the Milky Way. A well-determined orbit for Pal 5 would enable detailed study of these open questions. We present here the first CCD-based proper motion measurement of Pal 5 obtained using SDSS as a first epoch and new Large Binocular Telescope/Large Binocular Camera (LBC) images as a second, giving a baseline of 15 years. We perform relative astrometry, using SDSS as a distortion-free reference, and images of the cluster and also of the Pal 5 stream for the derivation of the distortion correction for LBC. The reference frame is made up of background galaxies. We correct for differential chromatic refraction using relations obtained from SDSS colors as well as from flux-calibrated spectra, finding that the correction relations for stars and for galaxies are different. We obtain μ{sub α} = −2.296 ± 0.186 mas yr{sup −1} and μ{sub δ} = −2.257 ± 0.181 mas yr{sup −1} for the proper motion of Pal 5. We use this motion, and the publicly available code galpy, to model the disruption of Pal 5 in different Milky Way models consisting of a bulge, a disk, and a spherical dark matter halo. Our fits to the observed stream properties (streak and radial velocity gradient) result in a preference for a relatively large Pal 5 distance of around 24 kpc. A slightly larger absolute proper motion than what we measure also results in better matches but the best solutions need a change in distance. We find that a spherical Milky Way model, with V{sub 0} = 220 km s{sup −1} and V{sub 20} {sub kpc}, i.e., approximately at the apocenter of Pal 5, of 218 km s{sup −1}, can match the data well, at least for our choice of disk and bulge parametrization.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-parameter distributed fiber optic sensing using double-Brillouin peak fiber in Brillouin optical time domain analysis

In this paper, we demonstrate a multi-parameter fiber sensing system based on stimulated Brillouin scattering in a double-Brillouin peak specialty fiber with enhanced Brillouin gain response. The amplitude level of the second Brillouin gain peak, which originated from the higher-order acoustic modes, has been improved with an approximately similar amplitude level to the first Brillouin gain peak from the fundamental acoustic mode. Compared to other multi-Brillouin peak fibers presented in the literature, the proposed fiber significantly reduces the measured Brillouin frequency shift error, thus improving strain and temperature accuracies. By utilizing the sensitivity values of the strain and temperature associated with each Brillouin gain spectrum (BGS) peak, a successful discriminative measurement of strain and temperature is performed with an accuracy of ±13 μɛ, and ±0.5 °C, respectively. The proposed double-Brillouin peak fiber appears to be a possible alternative to other multi-BGS peak fibers, for instance, large effective area fiber and dispersion compensating fibers, which are inherently accompanied by large measurement errors due to the weak Brillouin gain values originating from the higher-order acoustic modes. The demonstrated results show different strain and temperature coefficients of 47 kHz/µɛ, 1.15 MHz/°C for peak 1 and 51 kHz/µɛ, 1.37 MHz/°C for peak 2. Moreover, the enhanced BGS peak gains having nearly the same amplitude levels enable the discriminative measurement of strain and temperature. Such fibers in Brillouin interrogation eliminate the need for complex monitoring setups and reduce measurement errors. We recommend that for long-distance natural gas pipeline monitoring, where discriminative strain and temperature measurement is crucial, the proposed double-Brillouin peak fiber can be highly beneficial.

47 OTHER INSTRUMENTATION↗

Cybersecurity for Distance Relay Protection

This project is a DOE follow-up effort on the CREDC workshop held on September 13, 2018 in Cambridge, MA to discuss cybersecurity of distance relays, which considered the benefits, vulnerabilities and risk mitigations for the use of communication systems in power system protection. The objectives of this project are to define the taxonomy of relay protection and associated communications; define use cases describing approaches to reduce the cyber-attack surface on those protective relays; and evaluate the loss of operational functional capability from changes to communication coverage. Mitigating controls will also be evaluated to understand if there are other approaches to reduce attack surfaces while maintaining communications or partial communications. Distance relays are used to protect transmission lines of approximately 10 to 300 miles in length, by detecting short circuits (i.e., faults) on the lines and then tripping circuit breakers in the substation. Such protection systems are a subset of the power system and they incorporate sensing, logic and communication functions. Protection system exposure to cyberattack could be drastically limited by disconnecting relays from all vulnerable communication systems, but this may adversely impact overall power system performance in the absence of cyberattack. This project began with a use case analysis of protection systems with communications, as summarized in this report. It continued with modeling, testing and evaluation in a miniature power system (MPS), located in the Western Area Power Administration (WAPA) Electric Power Training Center (EPTC). The project also incorporated feedback from two industry meetings held in February and September 2019. The suggested next steps account for and complement the work already underway with DOE/CESER funding: 1. Study the performance of LCD and PC vs. PUTT, which is less reliant on communication system performance and GPS timing references. The PUTT scheme could prove to be more resilient to cyberattack or communications-related disruption. It could also be more tolerant of message re-routing with SDN/SDR communication systems. On the other hand, it will be more vulnerable to false tripping during dynamic events or to loss of the voltage signal. The optimum choice of scheme may depend on the specific power system and risk assessment. This study could provide a new template for evaluation based on business functions. 2. Research and develop new methods to detect and monitor distributed physical attacks, possibly using drones, video sensors, thermal sensors, machine learning and other advanced techniques. This will help mitigate the impact of cyberattack on the protection system, and will also help mitigate the impact of wild fires. 3. Implement a scalable PKI for use in electric utility protection systems. This will encourage widespread adoption of secure authentication methods that are already available, but not widely used at present. This will help secure engineering access to the relays. 4. Investigate the use of SDN in combination with SDR to achieve better cybersecurity and electromagnetic security of the network, incorporating path variability. This would help secure both engineering access and peer-to-peer GOOSE messaging. 5. Perform additional testing, with operator evaluation of “red button” scenarios, PUTT vs. LCD, relay mis-operations, and other cyberattacks in the EPTC. This is an important advantage of testing in the EPTC rather than by computer simulation or even hardware-in-the-loop simulation; the EPTC is already dedicated to managing the situational awareness, operator response times and other human impacts. One of the project objectives was to settle on a common nomenclature for this problem space. We have concluded that the OSI layer model, supplemented by ANSI device numbers and other IEEE standards, is already well-accepted by the industry. The IEEE PSRC knowledge base provides a great deal of public information

24 POWER TRANSMISSION AND DISTRIBUTION↗

Pilot study for compact microbeam radiation therapy using a carbon nanotube field emission micro-CT scanner

Purpose: Microbeam radiation therapy (MRT) is defined as the use of parallel, microplanar x-ray beams with an energy spectrum between 50 and 300 keV for cancer treatment and brain radiosurgery. Up until now, the possibilities of MRT have mainly been studied using synchrotron sources due to their high flux (100s Gy/s) and approximately parallel x-ray paths. The authors have proposed a compact x-ray based MRT system capable of delivering MRT dose distributions at a high dose rate. This system would employ carbon nanotube (CNT) field emission technology to create an x-ray source array that surrounds the target of irradiation. Using such a geometry, multiple collimators would shape the irradiation from this array into multiple microbeams that would then overlap or interlace in the target region. This pilot study demonstrates the feasibility of attaining a high dose rate and parallel microbeam beams using such a system. Methods: The microbeam dose distribution was generated by our CNT micro-CT scanner (100μm focal spot) and a custom-made microbeam collimator. An alignment assembly was fabricated and attached to the scanner in order to collimate and superimpose beams coming from different gantry positions. The MRT dose distribution was measured using two orthogonal radiochromic films embedded inside a cylindrical phantom. This target was irradiated with microbeams incident from 44 different gantry angles to simulate an array of x-ray sources as in the proposed compact CNT-based MRT system. Finally, phantom translation in a direction perpendicular to the microplanar beams was used to simulate the use of multiple parallel microbeams. Results: Microbeams delivered from 44 gantry angles were superimposed to form a single microbeam dose distribution in the phantom with a FWHM of 300μm (calculated value was 290 μm). Also, during the multiple beam simulation, a peak to valley dose ratio of ∼10 was found when the phantom translation distance was roughly 4x the beam width. The first prototype CNT-based x-ray tube dedicated to the development of compact MRT technology development was proposed and planned based on the preliminary experimental results presented here and the previous corresponding Monte Carlo simulations. Conclusions: The authors have demonstrated the feasibility of creating microbeam dose distributions at a high dose rate using a proposed compact MRT system. The flexibility of CNT field emission x-ray sources could possibly bring compact and low cost MRT devices to the larger research community and assist in the translational research of this promising new approach to radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗