Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed computation”

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 541 records · Page 30

On Eigenfunctions of the Fourier Transform

A nontrivial example of an eigenfunction in the sense of the theory of distributions for the planar Fourier transform was described by the authors in their previous work. In this paper, a method for obtaining other eigenfunctions is proposed. Positive homogeneous distributions in ℝ{sup n} of order −n/2 are considered, and it is shown that F(ω)|x|{sup −n/2}, |ω| = 1, is an eigenfunction in the sense of the theory of distributions of the Fourier transform if and only if F(ω) is an eigenfunction of a certain singular integral operator on the unit sphere of ℝ{sup n}. Since Y{sub m,n}{sup (k)}(ω)|x|{sup −n/2}, where Y{sub m,n}{sup (k)} denote the spherical functions of order m in ℝ{sup n}, are eigenfunctions of the Fourier transform, it follows that Y{sub m,n}{sup (k)} are eigenfunctions of the above-mentioned singular integral operator. In the planar case, all eigenfunctions of the Fourier transform of the form F(ω)|x|{sup −1} are described by means of the Fourier coefficients of F(ω).

97 MATHEMATICS AND COMPUTING↗

Clustering at Massive Scale

ClaMS provides hierarchical clustering technology for use on massive, high-dimensional datasets that require distributed memory for processing. The algorithm employed is inspired by the popular HDBSCAN algorithm but makes use of computational kernels better suited for distributed computing. ClaMS is built on scalable nearest neighbor graph construction, metric forest completion, and approximate minimum spanning tree techniques.

Stanley, ThomasA [Lawrence Livermore National Labo↗

Using Computational Fluid Dynamics of Wind Simulations Coupled With Weather Data to Calculate Dynamic Line Ratings

Dynamic line rating is a technology that allows for the rating of electrical conductors to be calculated based on local weather conditions rather than using ”worst-case” assumptions of weather conditions. The static ratings are typically based on weather conditions that are conservative in nature, and thus dynamic line rating may provide room to increase the current capacity. The local weather conditions may estimate conditions more accurately near the transmission lines accounting for terrain such as nearby hills. One way to estimate the local weather conditions is through the use of computational fluid dynamics. The results of a wind field simulation can be coupled with sparsely located weather stations to provide an accurate assessment of the wind speeds along the path of the conductor. Here, the wind field simulations are calculated using a steady state Reynolds-averaged Navier-Stokes model. These results along with local solar radiance measurements can be used to provide estimates of the steady state ampacity with the standard IEEE equations. An advanced test case for this process involves performing these simulations in a complex terrain, Hells Canyon, the deepest river gorge in North America. The results show that by using weather-based sensors without considering localized wind conditions, the available ampacity may be over-predicted significantly.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Modular chip-integrated photonic control of artificial atoms in diamond waveguides

A central goal in creating long-distance quantum networks and distributed quantum computing is the development of interconnected and individually controlled qubit nodes. Atom-like emitters in diamond have emerged as a leading system for optically networked quantum memories, motivating the development of visible-spectrum, multi-channel photonic integrated circuit (PIC) systems for scalable atom control. However, it has remained an open challenge to realize optical programmability with a qubit layer that can achieve high optical detection probability over many optical channels. Here, we address this problem by introducing a modular architecture of piezoelectrically actuated atom-control PICs (APICs) and artificial atoms embedded in diamond nanostructures designed for high-efficiency free-space collection. The high-speed four-channel APIC is based on a splitting tree mesh with triple-phase shifter Mach–Zehnder interferometers. This design simultaneously achieves optically broadband operation at visible wavelengths, high-fidelity switching (>40dB) at low voltages, submicrosecond modulation timescales (>30MHz), and minimal channel-to-channel crosstalk for repeatable optical pulse carving. Via a reconfigurable free-space interconnect, we use the APIC to address single silicon vacancy color centers in individual diamond waveguides with inverse tapered couplers, achieving efficient single photon detection probabilities (∼15%) and second-order autocorrelation measurements g (2) (0)<0.14 for all channels. The modularity of this distributed APIC–quantum memory system simplifies the quantum control problem, potentially enabling further scaling to thousands of channels.

47 OTHER INSTRUMENTATION↗

Analysis and Visualization of Discrete Fracture Networks Using a Flow Topology Graph

We present an analysis and visualization prototype using the concept of a flow topology graph (FTG) for characterization of flow in constrained networks, with a focus on discrete fracture networks (DFN), developed collaboratively by geoscientists and visualization scientists. Our method allows users to understand and evaluate flow and transport in DFN simulations by computing statistical distributions, segment paths of interest, and cluster particles based on their paths. The new approach enables domain scientists to evaluate the accuracy of the simulations, visualize features of interest, and compare multiple realizations over a specific domain of interest. Geoscientists can simulate complex transport phenomena modeling large sites for networks consisting of several thousand fractures without compromising the geometry of the network. However, few tools exist for performing higher-level analysis and visualization of simulated DFN data. The prototype system we present addresses this need. Here, we demonstrate its effectiveness for increasingly complex examples of DFNs, covering two distinct use cases - hydrocarbon extraction from unconventional resources and transport of dissolved contaminant from a spent nuclear fuel repository.

97 MATHEMATICS AND COMPUTING↗

Rare behavior of growth processes via umbrella sampling of trajectories

We compute probability distributions of trajectory observables for reversible and irreversible growth processes. These results reveal a correspondence between reversible and irreversible processes, at particular points in parameter space, in terms of their typical and atypical trajectories. Thus key features of growth processes can be insensitive to the precise form of the rate constants used to generate them, recalling the insensitivity to microscopic details of certain equilibrium behavior. We obtained these results using a sampling method, inspired by the "s-ensemble" large-deviation formalism, that amounts to umbrella sampling in trajectory space. The method is a simple variant of existing approaches, and applies to ensembles of trajectories controlled by the total number of events. It can be used to determine large-deviation rate functions for trajectory observables in or out of equilibrium.

97 MATHEMATICS AND COMPUTING↗

Non-Diffusive Volume Advection with A High Order Interface Reconstruction Method

We show that non-diffusive volume advection in two-dimensions is achieved with several benchmark problems using a newly developed high-order volume of fluids (VOF) interface reconstruction method. (1) A new VOF interface reconstruction method using circular/corner facets (linear facets are a degenerate case of arcs). We create a circular interface facet in each mixed zone by matching neighbor volume with a hybrid Newton’s-bisection method and the local solution is final. In the general case, the new VOF interface reconstruction has 3rd order accuracy and can be easily made seamless. The new method addresses intrinsic issues with Young’s method such as gaps between interface facets in the case of a curved interface, and inability to define curvature nor identify corners. (2) A non-diffusive volume advection scheme. In an ALE advection step, a well-defined interface can be carried over through a Lagrange step and used to compute volume distribution into a relaxed mesh. Then, an interface reconstruction step is performed to redefine the interface in the relaxed mesh. We must point out that the interface carried over is also a solution of interface reconstruction because all the volume fractions in the relaxed mesh are naturally matched. We provide an interface tracking method compatible with our reconstruction scheme, where it is granted to use the prior info as an initial guess to capture sub-mesh resolution features. As a result, we are able to treat multiple facets inside a single mixed cell and obtain highly accurate, non-diffusive solution for advection problems with rather coarse meshes. We show our solutions for two-dimensional incompressible flows with two materials with a) the X + O diagonal translation; b) the Zalesak rotational test; and c) the single vortex spiral test.

97 MATHEMATICS AND COMPUTING↗

SU-E-I-07: An Improved Technique for Scatter Correction in PET

Purpose: In positron emission tomography (PET), the single scatter simulation (SSS) algorithm is widely used for scatter estimation in clinical scans. However, bias usually occurs at the essential steps of scaling the computed SSS distribution to real scatter amounts by employing the scatter-only projection tail. The bias can be amplified when the scatter-only projection tail is too small, resulting in incorrect scatter correction. To this end, we propose a novel scatter calibration technique to accurately estimate the amount of scatter using pre-determined scatter fraction (SF) function instead of the employment of scatter-only tail information. Methods: As the SF depends on the radioactivity distribution and the attenuating material of the patient, an accurate theoretical relation cannot be devised. Instead, we constructed an empirical transformation function between SFs and average attenuation coefficients based on a serious of phantom studies with different sizes and materials. From the average attenuation coefficient, the predicted SFs were calculated using empirical transformation function. Hence, real scatter amount can be obtained by scaling the SSS distribution with the predicted SFs. The simulation was conducted using the SimSET. The Siemens Biograph™ 6 PET scanner was modeled in this study. The Software for Tomographic Image Reconstruction (STIR) was employed to estimate the scatter and reconstruct images. The EEC phantom was adopted to evaluate the performance of our proposed technique. Results: The scatter-corrected image of our method demonstrated improved image contrast over that of SSS. For our technique and SSS of the reconstructed images, the normalized standard deviation were 0.053 and 0.182, respectively; the root mean squared errors were 11.852 and 13.767, respectively. Conclusion: We have proposed an alternative method to calibrate SSS (C-SSS) to the absolute scatter amounts using SF. This method can avoid the bias caused by the insufficient tail information and therefore improve the accuracy of scatter estimation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

External beam planning module of Eclipse for external beam radiation therapy

Eclipse is a 3-dimensional (3D) treatment planning system for radiation therapy offered by Varian Medical Systems, Inc. The system has the network connectivity for the electronic transfer of image datasets and digital data communication among different equipment. The scope of this project for this special issue of Medical Dosimetry on 3D treatment planning systems is the assessment of planning tools in the external beam planning module of Eclipse to generate optimized treatment plans for patients undergoing external beam radiation therapy. This treatment planning system is relatively mature to be able to generate (1) simple treatment plans, (2) conformal radiation therapy plans, (3) static intensity-modulated radiation therapy (IMRT) plans, (4) volumetric-modulated arc therapy (VMAT) plans, and (5) treatment plans for electron beam therapy. The treatment planning tools are relatively plentiful to assist in the radiation therapy treatment planning. Some new features have been incorporated in the latest version and are helpful for making high-quality treatment plans. However, the location of the tools is not intuitive, and hence, familiarity with the user interface is essential to the efficient use of the treatment planning system. In addition, there are a number of dose algorithms available for the computation of dose distributions. The understanding of each dose computation algorithm is essential for the optimal use of this treatment planning system.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Coupling computational fluid dynamics with the high resolution rapid refresh model for forecasting dynamic line ratings

This study looks at forecasted dynamic line ratings in southern Idaho by using data from the High Resolution Rapid Refresh (HRRR) model for forecasted weather conditions. The HRRR model can provide accurate 18- hour forecasts with a 15-minute temporal resolution. Typical static ratings used for over head transmission lines use overly conservative assumptions for local weather conditions, such as using the maximum solar irradiance or temperature measured during the summer combined with low wind speed for an entire season. The HRRR forecast model used here has a high spatial resolution to provide local forecast conditions along the entire length of a transmission line. The area that is of interest in this study is in southern Idaho spanning a total of 15,000 square kilometers. The forecasted weather data is coupled with a computational fluid dynamics model of the wind in the region to improve the resolution to an even finer scale for resolving local convective cooling rates on each individual transmission line midpoint span. This high fidelity approach can be used to find the minimum ampacity limit across all given midpoints of a transmission line to determine the absolute minimum to be used for the line rating. The ability to increase the line rating above the conservative static approach with the forecasted weather data provides a large potential to alleviate congestion and provide data for the utility market transactions. Furthermore, this study shows that for the region of interest, using forecasted weather data coupled with CFD modeling will calculate DLR ampacity rating above static about 90% of the time, with only a small relative error in the forecasted ampacity over time.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Assessing DER Network Cybersecurity Defences in a Power-Communication Co-Simulation Environment

Increasing penetrations of interoperable distributed energy resources (DER) in the electric power system are expanding the power system attack surface. Maloperation or malicious control of DER equipment can now cause substantial disturbances to grid operations. Fortunately, many options exist to defend and limit adversary impact on these newly-created DER communication networks, which typically traverse the public internet. However, implementing these security features will increase communication latency, thereby adversely impacting real-time DER grid support service effectiveness. In this work, a collection of software tools called SCEPTRE were used to create a co-simulation environment where SunSpec-compliant PV inverters were deployed as virtual machines and interconnected to simulated communication network equipment. Network segmentation, encryption, and moving target defence security features were deployed on the control network to evaluate their influence on cybersecurity metrics and power system performance. The results indicated that adding these security features did not impact DER-based grid control systems but improved the cybersecurity posture of the network when implemented appropriately.

97 MATHEMATICS AND COMPUTING↗

Fluctuating magnetic field induced resonant activation

In this paper, we have studied the properties of a Brownian particle at stationary state in the presence of a fluctuating magnetic field. Time dependence of the field makes the system thermodynamically open. As a signature of that the steady state distribution function becomes function of damping strength, intensity of fluctuations and constant parts of the applied magnetic field. It also depends on the correlation time of the fluctuating magnetic field. Our another observation is that the random magnetic field can induce the resonant activation phenomenon. Here correlation time is increased under the fixed variance of the fluctuating field. But if the correlation time (τ) increases under the fixed field strength then the mean first passage time rapidly grows at low τ and it almost converges at other limit. This is sharp contrast to the usual colored noise driven open system case where the mean first passage time diverges exponentially. We have also observed that a giant enhancement of barrier crossing rate occurs particularly at large strength of constant parts of the applied magnetic field even for very weak fluctuating magnetic field. Finally, break down of the Arrhenius result and disappearance of the Kramers’ turn over phenomenon may occur in the presence of a fluctuating magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Forecasting Dynamic Line Rating with Spatial Variation Considerations

Dynamic line rating (DLR) is a technology that allows the ampacity of an electrical conductor to be calculated using real-time or forecasted weather conditions. Historically, the ampacity of a conductor has been determined using a static line rating method which assumes conservative weather assumptions. Therefore, not only can DLR give a more accurate measurement of the true ampacity of a conductor, but it can also increase its ampacity during weather conditions with greater thermal mitigations. The two primary cooling factors in the ampacity calculations are wind speed and direction. In complex terrain, wind speed and direction can have large variations over short distances. Therefore, accurately identifying the limiting span of a transmission line requires high spatial resolution of the wind along its path. One solution is to install dense weather stations along their path, though this can become costly over long distances. Therefore, researchers have investigated the use of Computation Fluid Dynamic (CFD) simulations to accurately compute the wind field along the path of a transmission line and use these results to identify the limiting section of the conductor. This work presents a case study that evaluates the coupling of CFD simulations and forecasted weather simulations using the High-Resolution Rapid Refresh (HRRR) model points over a 2-year span within a region in south eastern Idaho. The primary goal of the work is the evaluation of the number of HRRR model points used, i.e., weather stations, along the path of the line and the accuracy of the resulting ampacity. This was done using 4, 10, 17, 26, and 35 HRRR model points along two transmission line paths. The results indicate that as the number of model points are increased, the DLR ampacity of the lines decrease, yet converge as more points are added and demonstrate little change with additional HRRR points. It is expected that these results can help transmission line operators identify the number of weather stations that must be installed when coupled with CFD simulations and DLR ampacity to ensure accurate ratings and safe operations.

17 WIND ENERGY↗

Using Computational Fluid Dynamics to Assess Dynamic Line Rating in Southern Idaho

The overall goal of the study is to combine computational fluid dynamics (CFD) simulations with weather data that is collected over a 1-year long period across southern Idaho to calculate the dynamic line ratings (DLR) of several transmission lines. These ratings can be compared to the static rating assumptions to show the potential for additional current carrying capacity along transmission lines through account for additional effects of convective cooling from the wind. The region of interest is quite large, sized at 144 km by 98 km; in order to cover this large region, the CFD simulations are split into four separate regions of over 80 million computational cells each. The weather data that was collected suggest initial assumptions used to make the static ratings were vastly over predicted with regards to summer temperatures and solar irradiance. The weather data is ran through the General Line Ampacity State Solver (GLASS) tool that was developed by Idaho National Laboratory (INL) to parse large quantities of weather data across thousands of transmission midpoint spans in the region. For shorter lines, the dynamic line ratings often show large improvement over static ratings. The case study here shows over 95% of the time DLR values are above static values on the short transmission lines. However, for long lines, a conundrum occurs, where due to the large number of weather stations associated with the line, there is often a single weather station which may read local wind speeds of zero. This often defaults the line rating to be limited by natural convective cooling and does not show as much improvement over static as might be expected from the weather data if DLR calculations are used with a minimum ampacity across all midpoints. Other methods with transient state calculations should be utilized for long lines.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SU-F-19A-01: APBI Brachytherapy Treatment Planning: The Impact of Heterogeneous Dose Calculations

Purpose: To quantify the difference in dose to PTV-EVAL and OARs (skin and rib) as calculated by (TG43) and heterogeneous calculations (CCC). Methods: 25 patient plans (5 Contura and 20 SAVI) were selected for analysis. Clinical dose distributions were computed with a commercially available treatment planning algorithm (TG43-D-(w,w)) and then recomputed with a pre-clinical collapsed cone convolution algorithm (CCCD-( m,m)). PTV-EVAL coverage (V90%, V95%), and rib and skin maximum dose were compared via percent difference. Differences in dose to normal tissue (V150cc, V200cc of PTV-EVAL) were also compared. Changes in coverage and maximum dose to organs at risk are reported in percent change, (100*(TG43 − CCC) / TG43)), and changes in maximum dose to normal tissue are absolute change in cc (TG43 − CCC). Results: Mean differences in V90, V95, V150, and V200 for the SAVI cases were −0.2%, −0.4%, −0.03cc, and −0.14cc, respectively, with maximum differences of −0.78%, −1.7%, 1.28cc, and 1.01cc, respectively. Mean differences in the 0.1cc dose to the rib and skin were −1.4% and −0.22%, respectively, with maximum differences of −4.5% and 16%, respectively. Mean differences in V90, V95, V150, and V200 for the Contura cases were −1.2%, −2.1%, −1.8cc, and −0.59cc, respectively, with maximum differences of −2.0%, −3.16%, −2.9cc, and −0.76cc, respectively. Mean differences in the 0.1cc dose to the rib and skin were −2.6% and −3.9%, respectively, with maximum differences of −3.2% and −5.7%, respectively. Conclusion: The effects of translating clinical knowledge based on D-(w,w) to plans reported in D-(m,m) are minimal (2% or less) on average, but vary based on the type and placement of the device, source, and heterogeneity information.

60 APPLIED LIFE SCIENCES↗

The cross-correlation between 21 cm intensity mapping maps and the Lyα forest in the post-reionization era

We investigate the cross-correlation signal between 21cm intensity mapping maps and the Lyα forest in the fully non-linear regime using state-of-the-art hydrodynamic simulations. The cross-correlation signal between the Lyα forest and 21cm maps can provide a coherent and comprehensive picture of the neutral hydrogen (HI) content of our Universe in the post-reionization era, probing both its mass content and volume distribution. We compute the auto-power spectra of both fields together with their cross-power spectrum at z = 2.4 and find that on large scales the fields are completely anti-correlated. This anti-correlation arises because regions with high (low) 21cm emission, such as those with a large (low) concentration of damped Lyα systems, will show up as regions with low (high) transmitted flux. We find that on scales smaller than k ≅ 0.2 h Mpc{sup −1} the cross-correlation coefficient departs from −1, at a scale where non-linearities show up. We use the anisotropy of the power spectra in redshift-space to determine the values of the bias and of the redshift-space distortion parameters of both fields. We find that the errors on the value of the cosmological and astrophysical parameters could decrease by 30% when adding data from the cross-power spectrum, in a conservative analysis. Our results point out that linear theory is capable of reproducing the shape and amplitude of the cross-power up to rather non-linear scales. Finally, we find that the 21cm-Lyα cross-power spectrum can be detected by combining data from a BOSS-like survey together with 21cm intensity mapping observations by SKA1-MID with a S/N ratio higher than 3 in k element of [0.06,1] h Mpc{sup −1}. We emphasize that while the shape and amplitude of the 21cm auto-power spectrum can be severely affected by residual foreground contamination, cross-power spectra will be less sensitive to that and therefore can be used to identify systematics in the 21cm maps.

79 ASTRONOMY AND ASTROPHYSICS↗

Prostate-Specific Membrane Antigen Positron Emission Tomography–Computed Tomography for Prostate Cancer: Distribution of Disease and Implications for Radiation Therapy Planning

Purpose: To explore the prostate-specific membrane antigen (PSMA)–avid distribution of prostate cancer (PC) on positron emission tomography (PET), both at the time of initial diagnosis and at the time of relapse after definitive local treatment. Methods and Materials: A total of 179 PSMA PET scans in patients with nil or ≤3 lesions on conventional imaging were retrospectively categorized into 3 subgroups: group A, high-risk PC with no prior definitive therapy (n=34); group B, prior prostatectomy (n=75); and group C, prior radiation therapy (n=70). The numbers and locations of the PSMA-avid lesions were mapped. The PSMA-positive lesions were identified subjectively by a nuclear medicine physician on the basis of clinical experience and taking into account the recent literature and artefacts. Results: A total of 893 PSMA-avid lesions were identified; at least 1 lesion was detected in 80% of all scans. A high detection rate was present even at very low serum PSA levels (eg, at PSA ≤0.20 ng/mL in group B, the detection rate was 46%). Thirty-eight percent of studies revealed extrapelvic disease (41%, 31%, and 46% in groups A, B, and C, respectively). Almost one-third of all studies showed only oligometastases (24%, 36%, and 31% in groups A, B, and C, respectively). A large proportion of these (40%) were a solitary lesion. Conclusions: Prostate-specific membrane antigen PET demonstrated a large number of otherwise unknown metastatic lesions. Therefore we recommend PSMA PET for more accurate assessment of disease burden in initial staging of high-risk PC, as well as for restaging in patients with prostate-specific antigen relapse after primary therapies. Furthermore, a high proportion of oligometastases on PSMA PET provides a prime opportunity to investigate the role of targeted local therapies for oligometastatic PCs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-T-618: Plan Robustness Study of Volumetric-Modulated Arc Therapy Vs. Intensity-Modulated Radiation Therapy for Head and Neck Cancer

Purpose: Lack of plan robustness may contribute to local failure in volumetric-modulated arc therapy (VMAT) to treat head and neck (H&N) cancer. Thus we compared plan robustness of VMAT with intensity-modulated radiation therapy (IMRT). Methods: VMAT and IMRT plans were created for 9 H&N cancer patients. For each plan, six new perturbed dose distributions were computed — one each for ± 3mm setup deviations along the S-I, A-P and L-R directions. We used three robustness quantification tools: (1) worst-case analysis (WCA); (2) dose-volume histograms (DVHs) band (DVHB); and (3) root-mean-square-dose deviation (RMSD) volume histogram (DDVH). DDVH represents the relative volume (y) on the vertical axis and the RMSD (x) on the horizontal axis. Similar to DVH, this means that y% of the volume of the indicated structure has the RMSD at least x Gy[RBE].The width from the first two methods at different target DVH indices (such as D95 and D5) and the area under the DDVH curves (AUC) for the target were used to indicate plan robustness. In these robustness quantification tools, the smaller the value, the more robust the plan is. Plan robustness evaluation metrics were compared using Wilcoxon test. Results: DVHB showed the width at D95 from IMRT to be larger than from VMAT (unit Gy) [1.59 vs 1.18 (p=0.49)], while the width at D5 from IMRT was found to be slightly larger than from VMAT [0.59 vs 0.54 (p=0.84)]. WCA showed similar results [D95: 3.28 vs 3.00 (p=0.56); D5: 1.68 vs 1.95 (p=0.23)]. DDVH showed the AUC from IMRT to be slightly smaller than from VMAT [1.13 vs 1.15 (p=0.43)]. Conclusion: VMAT plan robustness is comparable to IMRT plan robustness. The plan robustness conclusions from WCA and DVHB are DVH parameter dependent. On the other hand DDVH captures the overall effect of uncertainties on the dose to a volume of interest. NIH/NCI K25CA168984; Eagles Cancer Research Career Development; The Lawrence W. and Marilyn W. Matteson Fund for Cancer Research Mayo ASU Seed Grant; The Kemper Marley Foundation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗