Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FALLON”

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 55 records · Page 3

SU-G-TeP1-15: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation

Purpose: To develop a Graphic Processor Unit (GPU) accelerated deterministic solution to the Linear Boltzmann Transport Equation (LBTE) for accurate dose calculations in radiotherapy (RT). A deterministic solution yields the potential for major speed improvements due to the sparse matrix-vector and vector-vector multiplications and would thus be of benefit to RT. Methods: In order to leverage the massively parallel architecture of GPUs, the first order LBTE was reformulated as a second order self-adjoint equation using the Least Squares Finite Element Method (LSFEM). This produces a symmetric positive-definite matrix which is efficiently solved using a parallelized conjugate gradient (CG) solver. The LSFEM formalism is applied in space, discrete ordinates is applied in angle, and the Multigroup method is applied in energy. The final linear system of equations produced is tightly coupled in space and angle. Our code written in CUDA-C was benchmarked on an Nvidia GeForce TITAN-X GPU against an Intel i7-6700K CPU. A spatial mesh of 30,950 tetrahedral elements was used with an S4 angular approximation. Results: To avoid repeating a full computationally intensive finite element matrix assembly at each Multigroup energy, a novel mapping algorithm was developed which minimized the operations required at each energy. Additionally, a parallelized memory mapping for the kronecker product between the sparse spatial and angular matrices, including Dirichlet boundary conditions, was created. Atomicity is preserved by graph-coloring overlapping nodes into separate kernel launches. The one-time mapping calculations for matrix assembly, kronecker product, and boundary condition application took 452±1ms on GPU. Matrix assembly for 16 energy groups took 556±3s on CPU, and 358±2ms on GPU using the mappings developed. The CG solver took 93±1s on CPU, and 468±2ms on GPU. Conclusion: Three computationally intensive subroutines in deterministically solving the LBTE have been formulated on GPU, resulting in two orders of magnitude speedup. Funding support from Natural Sciences and Engineering Research Council and Alberta Innovates Health Solutions. Dr. Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗

SU-G-JeP2-09: Minimal Skin Dose Increase in Longitudinal Rotating Biplanar Linac-MR Systems: Examination of Radiation Energy and Flattening Filter Design

Purpose: To quantify increase in entrance skin-dose due to magnetic fields of the Alberta longitudinal linac-MR by examining the effect of radiation energy and flattening filter, using Monte Carlo calculations and accurate 3-D models of the magnetic field. Methods: The 3-D magnetic fields generated by the bi-planar Linac-MR are calculated with FEM using Opera-3D. BEAMnrc simulates the particle phase-space in the presence of the rapidly decaying fringe field of 0.5T MRI assembled with a Varian 600C linac with an isocentre distance of 130 cm for 6 MV and 10 MV beams. Skin doses are calculated at an average depth of 70 µm using DOSXYZnrc with varying SSDs and field sizes. Furthermore, flattening filters are reshaped to compensate for the significant drop in dose rate due to increased SAD of 130 cm and skin-doses are evaluated. Results: The confinement effect of the MRI fringe field on the contaminant electrons is minimal. For SSDs of 100 – 120 cm the increase in skin dose is ∼6% – 19% and ∼1% – 9% for the 6 and 10 MV beams, respectively. For 6MV, skin dose increases from ∼10.5% to 1.5%. for field-size increases of 5×5 cm2 to 20×20 cm2. For 10 MV, skin dose increases by ∼6% for a 5×5 cm2 field, and decreases by ∼1.5% for a 20×20 cm2 field. The reshaped flattening filter increases the dose rate from 355 MU/min to 529 MU/min (6 MV) or 604 MU/min (10 MV), while the skin-dose increases by only an additional ∼2.6% (all percent increases in skin dose are relative to Dmax). Conclusion: There is minimal increase in the entrance skin dose and minimal/no decrease in the dose rate of the Alberta longitudinal linac-MR system. There is even lower skin-dose increase at 10 MV. Funding: Alberta Innovates - Health Solutions (AIHS) Conflict of Interest: Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization)

60 APPLIED LIFE SCIENCES↗

SU-G-JeP2-15: Proton Beam Behavior in the Presence of Realistic Magnet Fields

Purpose: To investigate the effects of magnetic fields on proton therapy beams for integration with MRI. Methods: 3D magnetic fields from an open-bore superconducting MRI model (previously developed by our group) and 3D magnetic fields from an in-house gradient coil design were applied to various mono energetic proton pencil beam (80MeV to 250MeV) simulations. In all simulations, the z-axis of the simulation geometry coincided with the direction of the B0 field and magnet isocentre. In each simulation, the initial beam trajectory was varied. The first set of simulations performed was based on analytic magnetic force equations (analytic simulations), which could be rapidly calculated yet were limited to propagating proton beams in vacuum. The second set is full Monte Carlo (MC) simulations, which used GEANT4 MC toolkit. Metrics such as the beam position and dose profiles were extracted. Comparisons between the cases with and without magnetic fields present were made. Results: The analytic simulations served as verification checks for the MC simulations when the same simulation geometries were used. The results of the analytic simulations agreed with the MC simulations performed in vacuum. The presence of the MRI’s static magnetic field causes proton pencil beams to follow a slight helical trajectory when there were some initial off-axis components. The 80MeV, 150MeV, and 250MeV proton beams rotated by 4.9o, 3.6o, and 2.8o, respectively, when they reached z=0cm. The deflections caused by gradient coils’ magnetic fields show spatially invariant patterns with a maximum range of 0.5mm at z=0cm. Conclusion: This investigation reveals that both the MRI’s B0 and gradient magnetic fields can cause small but observable deflections of proton beams at energies studied. The MRI’s static field caused a rotation of the beam while the gradient coils’ fields effects were spatially invariant. Dr. B Gino Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization)

60 APPLIED LIFE SCIENCES↗

TU-H-BRA-04: A Novel Superconducting Magnet Design for Optimized Patient Access and Minimal SSD for Use in a Linac-MR Hybrid

Purpose: A prototype rotating hybrid MR imaging system and linac has been developed to allow for simultaneous imaging and radiation delivery parallel to B{sub 0}. However, the design of a compact magnet capable of rotation in a small vault with sufficient patient access and a typical clinical source-to-surface distance (SSD) is challenging. This work presents a novel superconducting magnet design that allows for a reduced SSD and ample patient access by moving the superconducting coils to the side of the yoke. The yoke and pole-plate structures are shaped to direct the magnetic flux appropriately. Methods: The surface of the pole plate for the magnet assembly is optimized. The magnetic field calculations required in this work are performed with the 3D finite element method software package Opera-3D. Each tentative design strategy is virtually modeled in this software package and externally controlled by MATLAB, with its key geometries defined as variables. The particle swarm optimization algorithm is used to optimize the variables subject to the minimization of a cost function. At each iteration, Opera-3D will solve the magnetic field solution over a field-of-view suitable for MR imaging and the degree of field uniformity will be assessed to calculate the value of the cost function associated with that iteration. Results: An optimized magnet assembly that generates a homogenous 0.2T magnetic field over an ellipsoid with large axis of 30 cm and small axes of 20 cm is obtained. Conclusion: The distinct features of this model are the minimal distance between the yoke’s top and the isocentre and the improved patient access. On the other hand, having homogeneity over an ellipsoid give us a larger field-of-view, essential for geometric accuracy of the MRI system. The increase of B{sub 0} from 0.2T in the present model to 0.5T is the subject of future work. Funding Sources: Alberta Innovates - Health Solutions (AIHS)| Disclosure and Conflict of Interest: B. Gino Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta biplanar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗

TU-H-BRA-09: Relationship Between B0 and the Contrast-To-Noise Ratio (CNR) of Tumour to Background for MRI/Radiotherapy Hybrids

Purpose: To investigate the relationship in MRI between B{sub 0} and the contrast-to-noise ratio (CNR) of various tumour/normal tissue pairs. This study is motivated by the current interest in MRI/radiotherapy hybrids, for which multiple magnetic field strengths have been proposed. CNR is the single most important parameter governing the ability of a system to identify a tumour in real time for treatment guidance. The MRI community has long since recognized that the SNR of a well-designed MR system is roughly proportional to B{sub 0}, the polarizing magnetic field. However, the CNR between two tissues is much more complicated - dependent not only on this signal behavior, but also on the different relaxation properties of the tissues. Methods: Experimentally-based models of B{sub 0}-dependant relaxation for various tumour and normal tissues from the literature were used in conjunction with signal equations for MR sequences suitable for rapid realtime imaging to develop field-dependent predictions for CNR. These CNR models were developed for liver, lung, breast, glioma, and kidney tumours for spoiled-gradient echo (SGE) and balanced steady-state free precession (bSSFP) sequences. Results: In all cases there was an improved CNR at lower fields compared to linear dependency. Further, in some tumour sites, the CNR at lower fields was found to be comparable to, or sometimes higher than those at higher fields (i.e. bSSFP CNR for glioma, kidney and liver tumours). Conclusion: Due to the variation of tissue relaxation parameters with field, lower B{sub 0} fields have been shown to perform as well or better (in terms of CNR) than higher fields for some tumour sites. In other sites this effect was less pronounced. It is the complex relationship between CNR and B{sub 0} that leads to greater CNR at 0.5 T for certain tumour types studied here for fast imaging. B. Gino Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization)

60 APPLIED LIFE SCIENCES↗

TH-AB-BRA-09: Stability Analysis of a Novel Dose Calculation Algorithm for MRI Guided Radiotherapy

Purpose: To determine the iterative deterministic solution stability of the Linear Boltzmann Transport Equation (LBTE) in the presence of magnetic fields. Methods: The LBTE with magnetic fields under investigation is derived using a discrete ordinates approach. The stability analysis is performed using analytical and numerical methods. Analytically, the spectral Fourier analysis is used to obtain the convergence rate of the source iteration procedures based on finding the largest eigenvalue of the iterative operator. This eigenvalue is a function of relevant physical parameters, such as magnetic field strength and material properties, and provides essential information about the domain of applicability required for clinically optimal parameter selection and maximum speed of convergence. The analytical results are reinforced by numerical simulations performed using the same discrete ordinates method in angle, and a discontinuous finite element spatial approach. Results: The spectral radius for the source iteration technique of the time independent transport equation with isotropic and anisotropic scattering centers inside infinite 3D medium is equal to the ratio of differential and total cross sections. The result is confirmed numerically by solving LBTE and is in full agreement with previously published results. The addition of magnetic field reveals that the convergence becomes dependent on the strength of magnetic field, the energy group discretization, and the order of anisotropic expansion. Conclusion: The source iteration technique for solving the LBTE with magnetic fields with the discrete ordinates method leads to divergent solutions in the limiting cases of small energy discretizations and high magnetic field strengths. Future investigations into non-stationary Krylov subspace techniques as an iterative solver will be performed as this has been shown to produce greater stability than source iteration. Furthermore, a stability analysis of a discontinuous finite element space-angle approach (which has been shown to provide the greatest stability) will also be investigated. Dr. B Gino Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization)

60 APPLIED LIFE SCIENCES↗

Training dataset and results for geothermal exploration artificial intelligence, applied to Brady Hot Springs and Desert Peak

The submission includes the labeled datasets, as ESRI Grid files (.gri, .grd) used for training and classification results for our machine leaning model: - brady_som_output.gri, brady_som_output.grd, brady_som_output.* - desert_som_output.gri, desert_som_output.grd, desert_som_output.* The data corresponds to two sites: Brady Hot Springs and Desert Peak, both located near Fallon, NV. Input layers include: - Geothermal: Labeled data (0: Non-geothermal; 1: Geothermal) - Minerals: Hydrothermal mineral alterations, as a result of spectral analysis using Chalcedony, Kaolinite, Gypsum, Hematite and Epsomite - Temperature: Land surface temperature (% of times a pixel was classified as "Hot" by K-Means) - Faults: Fault density with a 300mradius - Subsidence: PSInSAR results showing subsidence displacement of more than 5mm - Uplift: PSInSAR results showing subsidence displacement of more than 5mm Also, the results of the classification using Brady and Desert Peak to build 2 Convolutional Neural Networks. These were applied to the training site as well as the other site, the results are in GeoTiff format. - brady_classification: Results of classification of the Brady-trained model - desert_classification: Results of classification of the Desert Peak-trained model - b2d_classification: Results of classification of Desert Peak using the Brady-trained model - d2b_classification: Results of classification of Brady using the Desert Peak-trained model

15 GEOTHERMAL ENERGY↗

A Physical Basis for M s -Yield Scaling in Hard Rock and Implications for Late-Time Damage of the Source Medium

Surface wave magnitude M s for a compilation of 72 nuclear tests detonated in hard rock media for which yields and burial depths have been reported in the literature is shown to scale with yield W as a + b × log[W], where a = 2.50 ± 0.08 and b = 0.80 ± 0.05. While the exponent b is consistent with an M s scaling model for fully coupled, normal containment-depth explosions, the intercept a is offset 0.45 magnitude units lower than the model. The cause of offset is important to understand in terms of the explosion source. Hard rock explosions conducted in extensional and compressional stress regimes show similar offsets, an indication that the tectonic setting in which an explosion occurs plays no role causing the offset. The scaling model accounts for the effects of source medium material properties on the generation of 20-s period Rayleigh wave amplitudes. Aided by thorough characterizations of the explosion and tectonic release sources, an extensive analysis of the 1963 October 26 Shoal nuclear test detonated in granite 27 miles southeast of Fallon NV shows that the offset is consistent with the predictions of a material damage source model related to non-linear stress wave interactions with the free surface. This source emits Rayleigh waves with polarity opposite to waves emitted by the explosion. The Shoal results were extended to analyse surface waves from the 1962 February 15 Hardhat nuclear test, the 1988 September 14 Soviet Joint Verification Experiment, and the anomalous 1979 August 18 northeast Balapan explosion which exhibits opposite polarity, azimuth-independent source component U1 compared to an explosion. Modelling these tests shows that Rayleigh wave amplitudes generated by the damage source are nearly as large as or larger than amplitudes from the explosion. As such, destructive interference can be drastic, introducing metastable conditions due to the sensitivity of reduced amplitudes to Rayleigh wave initial phase angles of the explosion and damage sources. This meta-stability is a likely source of scatter in M s -yield scaling observations. The agreement of observed scaling exponent b with the model suggests that the damage source strength does not vary much with yield, in contrast to explosions conducted in weak media where Ms scaling rates are greater than the model predicts, and the yield dependence of the damage source strength is significant. This difference in scaling behaviour is a consequence of source medium material properties.

58 GEOSCIENCES↗

Closure Report for the Project Shoal Area: Subsurface Corrective Action Unit 447 Shoal,Nevada, Site

The U.S. Department of Energy (DOE) Office of Legacy Management (LM) has prepared this Closure Report for subsurface Corrective Action Unit (CAU) 447 at the Shoal, Nevada, Site, which was formerly called the Project Shoal Area. The Shoal site was the location of an underground nuclear test in 1963 that resulted in residual radionuclide contamination at the detonation depth of 1211 feet (ft) (DOE 2015a). The environmental restoration process and Corrective Action strategy are conducted in accordance with the Federal Facility Agreement and Consent Order (FFACO) (FFACO 1996, as amended) and all applicable Nevada Division of Environmental Protection (NDEP) statues and regulations. Responsibility for the environmental site restoration was transferred from the DOE, National Nuclear Security Administration, Nevada Field Office, to LM on October 1, 2006. The Shoal site is south of U.S. Highway 50, approximately 30 miles southeast of Fallon, in Churchill County, Nevada (Figure 1). The site was selected for underground nuclear testing in 1961, as part of the Vela Uniform Program. Vela was a research and development program directed toward locating, detecting, and identifying underground and high-altitude nuclear detonations. The Vela Uniform program was specific to underground detonations (AEC 1964). The objective of Project Shoal was to detonate a nuclear device underground in an active seismic area to improve the United States’ ability to detect, identify, and locate underground nuclear detonations. The underground nuclear test at Shoal was a joint effort between the U.S. Department of Defense (DOD) and the U.S. Atomic Energy Commission (AEC), which was the predecessor agency to DOE. Figure 1 is a map showing the site location.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Enrichments for phototrophic bacteria and characterization by morphology and pigment analysis

The purpose of this investigation was to examine several sulfide containing environments for the presence of phototrophic bacteria and to obtain enriched cultures of some of the bacteria present. The field sites were Alum Rock State Park, the Palo Alto salt marsh, the bay area salt ponds, and Big Soda Lake (near Fallon, Nevada). Bacteria from these sites were characterized by microscopic examination, measurement of in vitro absorption spectra, and analysis of carotenoid pigments. Field observations at one of the bay area salt ponds, in which the salt concentration was saturating (about 30 percent NaCl) and the sediments along the shore of the pond covered with a gypsum crust, revealed a layer of purple photosynthetic bacteria under a green layer in the gypsum crust. Samples of this gypsum crust were taken to the laboratory to measure light transmission through the crust and to try to identify the purple photosynthetic bacteria present in this extremely saline environment.

Brune, D.↗

Terminating states in the positive-parity structures of As67

The energy levels and γ-ray decay scheme of the positive-parity states in the Tz=12 nucleus 67As have been studied by using the 40Ca(36Ar,2αp)67As reaction at a beam energy of 145 MeV. Two new band structures have been identified which can be connected to the previously known levels. The results for these bands are compared with configuration-dependent cranked Nilsson–Strutinsky calculations. The good level of agreement between theory and experiment suggests that these structures can be interpreted in terms of configurations that involve three g92 particles and that both possess noncollective terminating states.

59 ≤ A ≤ 89↗

Characterization of a gamma-ray tracking array: A comparison of GRETINA and Gammasphere using a 60 Co source

In this paper, we provide a formalism for the characterization of the tracking arrays with emphasis on the proper corrections required to extract their photopeak efficiencies and peak-to-total ratios. The methods are first applied to Gammasphere, a well characterized 4pi array based on the principle of Compton suppression, and subsequently to GRETINA. The tracking efficiencies are then discussed and some guidelines as to what clustering angle to use in the tracking algorithm are presented. It was possible, using GEANT4 simulations, to scale the measured efficiencies up to the expected values for the full 4pi implementation of GRETA.

efficiency measurements↗

Configuration mixing in 28 Mg and the 26 Mg (t,p) 28 Mg reaction

We have studied the 26 Mg(t,p) 28 Mg reaction. The nucleus 28 Mg lies between 24 Mg and 32 Mg, the latter of which is central to the “island of inversion” characterized by low-lying neutron fp-shell configurations. The present results show that the second-excited 0 + state in 28 Mg contains large fp-shell occupation for the neutrons. The experiment was performed in inverse kinematics using the HELIcal Orbit Spectrometer (HELIOS) at Argonne National Laboratory. Shell-model calculations using the SDPF-MU interaction yielded level energies and two-neutron transfer amplitudes that were used in one-step distorted-wave Furthermore, born approximation calculations to provide theoretical predictions of the proton angular distributions. In many, but not all cases, the data were in good agreement with theoretical predictions.

20 ≤ A ≤ 38↗

Rotational band structure in Mg 32

There is significant evidence supporting the existence of deformed ground states within the neutron-rich N ≈ 20 neon, sodium, and magnesium isotopes that make up what is commonly called the “island of inversion.” However, the rotational band structures, which are a characteristic fingerprint of a rigid nonspherical shape, have yet to be observed. Here, in this work, we report on a measurement and analysis of the yrast (lowest lying) rotational band in 32 Mg up to spin I = 6 + produced in a two-step projectile fragmentation reaction and observed using the state-of-the-art γ-ray tracking detector array, GRETINA (γ-ray energy tracking in-beam nuclear array). Large-scale shell-model calculations using the SDPF-U-MIX effective interaction show excellent agreement with the new data. Finally, a theoretical analysis of the spectrum of rotational states as a function of the pairing gap, together with cranked-shell-model calculations, provides intriguing evidence for a reduction in pairing correlations with increased angular momentum, also in line with the shell-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Strongly coupled rotational band in Mg 33

The “island of inversion” at N≈20 for the neon, sodium, and magnesium isotopes has long been an area of interest both experimentally and theoretically due to the subtle competition between 0p-0h and np-nh configurations leading to deformed shapes. However, the presence of rotational band structures, which are fingerprints of deformed shapes, have only recently been observed in this region. In this work, we report on a measurement of the low-lying level structure of Mg33 populated by a two-stage projectile fragmentation reaction and studied with the Gamma Ray Energy Tracking In-Beam Nuclear Array (GRETINA). The experimental level energies, ground-state magnetic moment, intrinsic quadrupole moment, and γ-ray intensities show good agreement with the strong-coupling limit of a rotational model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗