Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible 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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2010

A Theoretical Mechanism of Szilard Engine Function in Nucleic Acids and the Implications for Quantum Coherence in Biological Systems

Nucleic acids theoretically possess a Szilard engine function that can convert the energy associated with the Shannon entropy of molecules for which they have coded recognition, into the useful work of geometric reconfiguration of the nucleic acid molecule. This function is logically reversible because its mechanism is literally and physically constructed out of the information necessary to reduce the Shannon entropy of such molecules, which means that this information exists on both sides of the theoretical engine, and because information is retained in the geometric degrees of freedom of the nucleic acid molecule, a quantum gate is formed through which multi-state nucleic acid qubits can interact. Entangled biophotons emitted as a consequence of symmetry breaking nucleic acid Szilard engine (NASE) function can be used to coordinate relative positioning of different nucleic acid locations, both within and between cells, thus providing the potential for quantum coherence of an entire biological system. Theoretical implications of understanding biological systems as such 'quantum adaptive systems' include the potential for multi-agent based quantum computing, and a better understanding of systemic pathologies such as cancer, as being related to a loss of systemic quantum coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Inverse scattering theory: Inverse scattering series method for one dimensional non-compact support potential

The reversion of the Born-Neumann series of the Lippmann-Schwinger equation is one of the standard ways to solve the inverse acoustic scattering problem. One limitation of the current inversion methods based on the reversion of the Born-Neumann series is that the velocity potential should have compact support. However, this assumption cannot be satisfied in certain cases, especially in seismic inversion. Based on the idea of distorted wave scattering, we explore an inverse scattering method for velocity potentials without compact support. The strategy is to decompose the actual medium as a known single interface reference medium, which has the same asymptotic form as the actual medium and a perturbative scattering potential with compact support. After introducing the method to calculate the Green’s function for the known reference potential, the inverse scattering series and Volterra inverse scattering series are derived for the perturbative potential. Analytical and numerical examples demonstrate the feasibility and effectiveness of this method. Besides, to ensure stability of the numerical computation, the Lanczos averaging method is employed as a filter to reduce the Gibbs oscillations for the truncated discrete inverse Fourier transform of each order. Our method provides a rigorous mathematical framework for inverse acoustic scattering with a non-compact support velocity potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Nonlinear saturation of reverse shear Alfvén eigenmodes induced by kinetic electrons

A strong nonlinear saturation mechanism of the reverse shear Alfvén eigenmode (RSAE) induced by thermal electrons is observed in gyrokinetic $\delta$$f$ particle-in-cell simulation. This nonlinear effect occurs at moderate mode amplitude, $\delta$$B$ ⟂ /$B$ < 10 -3 , and is associated with the electron parallel streaming along the perturbed magnetic field lines, the term v $\|$ $\frac{δΒ_{⟂}}{B}$ ∙ ∇$\delta$$f$, in the electron drift-kinetic equation. In the case of an n = 4 RSAE, where n is the toroidal mode number, this magnetic fluttering nonlinearity leads to strong mode coupling and broadening of the saturation spectrum. The nonlinear components, which are highly damped through the interaction between thermal electrons and the parallel electric field, modifies the RSAE mode structure and suppresses energetic particle drive significantly. The split-weight scheme (Chen and Parker 2007 J. Comput. Phys. 220 839–55) is used to simulate kinetic electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fundamentals of hydrogen storage in nanoporous materials

Abstract Physisorption of hydrogen in nanoporous materials offers an efficient and competitive alternative for hydrogen storage. At low temperatures (e.g. 77 K) and moderate pressures (below 100 bar) molecular H 2 adsorbs reversibly, with very fast kinetics, at high density on the inner surfaces of materials such as zeolites, activated carbons and metal–organic frameworks (MOFs). This review, by experts of Task 40 ‘Energy Storage and Conversion based on Hydrogen’ of the Hydrogen Technology Collaboration Programme of the International Energy Agency, covers the fundamentals of H 2 adsorption in nanoporous materials and assessment of their storage performance. The discussion includes recent work on H 2 adsorption at both low temperature and high pressure, new findings on the assessment of the hydrogen storage performance of materials, the correlation of volumetric and gravimetric H 2 storage capacities, usable capacity, and optimum operating temperature. The application of neutron scattering as an ideal tool for characterising H 2 adsorption is summarised and state-of-the-art computational methods, such as machine learning, are considered for the discovery of new MOFs for H 2 storage applications, as well as the modelling of flexible porous networks for optimised H 2 delivery. The discussion focuses moreover on additional important issues, such as sustainable materials synthesis and improved reproducibility of experimental H 2 adsorption isotherm data by interlaboratory exercises and reference materials.

08 HYDROGEN↗
NASA NTRSDate not supplied

Variance Decomposition of MEDLI2 Reconstructed Heating Using Neural Networks

The Mars Entry, Descent, and Landing Instrumentation (MEDLI2) sensor suite collected data during entry of the Mars 2020 Perseverance rover into Mars’ atmosphere. This suite included a network of MEDLI2 Instrumented Sensor Plugs (MISPs). Each MISP was comprised of a cylinder made of Thermal Protection System (TPS) material with 1-3 embedded thermocouples (TCs), and it was flush mounted into the heatshield or backshell. Data from these in-depth TCs were used to reconstruct the aeroheating environment of the vehicle throughout entry. Surface heating was posed as an inverse problem, with the goal of estimating the surface heating by minimizing an objective function of the difference between MISP temperature measurements during flight and the temperature predictions derived from the Fully Implicit Ablation and Thermal response (FIAT) program. Given an aerothermal environment, FIAT calculates the material response and provides in-depth temperatures throughout the TPS material. To achieve the reverse, an internal tool called FIAT_Opt runs through multiple different environments until the output temperature at the TC depth closely matches the flight data. 95% confidence intervals on the reconstructed surface heating were obtained using Monte Carlo analysis, in which uncertainties in the thermocouple depth and the TPS material properties (e.g., density, thermal conductivity, heat capacity, emissivity) based on flight-lot material testing were included. A variance decomposition method using Sobol indices was employed to assess the sensitivity of the reconstructed peak heating to the TC placement and material property uncertainties. Variance decomposition was found to require tens of thousands of FIAT_Opt runs in order for the Sobol indices to converge. With a single FIAT_Opt run taking on the order of 40 minutes, the required number of computations would take months to complete, even if using multiple CPUs. To mitigate this problem, three machine learning models (ridge regression with cross-validation, random forest regression, and a deep neural network) were trained and tested using the 2000 Monte Carlo runs that were already completed. A subset of 1600 runs were used to train the model (i.e., training set), while the remaining 400 runs were used as the test set. The predictions from the deep neural network (DNN) on the test set showed nearly perfect agreement to the actual values computed with FIAT_Opt (R2 > 0.99). Using the DNN as a surrogate model, the variance decomposition using 50,000 runs was completed within minutes. The resulting Sobol indices showed that the reconstructed peak surface heating was most sensitive to the uncertainties in the thermal conductivity (ST = 0.37) and heat capacity (ST = 0.26). This method can be leveraged to provide requirements for material property measurements needed to improve the accuracy of surface heating prediction and ultimately lead to the reduction of design margins in the future. This presentation will include background on the MEDLI2 suite; the method used for inverse heating estimation; the way that material property uncertainties were accounted for using Monte Carlo analysis; a brief background on variance decomposition; the motivation for using machine learning in this context; how a neural network was trained on the data to enable variance decomposition in a fraction of the time; and the variance decomposition results for one of the MISPs.

Hannah Alpert↗
DOE OSTI2020

Amphiphile Organization in Organic Solutions: An Alternative Explanation for Small-Angle X-ray Scattering Features in Malonamide/Alkane Mixtures

In this work, the role of different intermolecular interactions in the aggregation of amphiphiles in an organic solvent is studied for systems of relevance to liquid-liquid extraction (LLE), a chemical process used to selectively recover metals from complex mixtures. Of specific interest is the role, or lack thereof, of hydrogen bonding, which is often assumed to be a main driver of the organic phase structural organization that has been linked to separation efficacy. Toward that end, a series of malonamide extractants in n-dodecane have been studied in the absence of any extracted aqueous solutes, including water. The series of extractants includes N,N'-dimethyl-N,N'-dibutyltetradecylmalonamide (DMDBTDMA), two of its homologs, and N,N'-dimethyl-N,N'-dioctylhexylethoxymalonamide (DMDOHEMA). This simplified model LLE system enables systematic investigation of the role of dipole-dipole and alkyl tail steric interactions in amphiphile aggregation. Small-angle X-ray scattering (SAXS) profiles computed from molecular dynamics trajectories are in good agreement with experimental SAXS data. Molecular dynamics simulations show that malonamide aggregation results from dipole-driven self-association and lacks characteristic aggregate sizes. Mid-q correlation peaks in the SAXS profiles emerge at high concentration for each malonamide. In those densely packed solutions, the correlation peaks are observed to result from alkyl tail-induced spacing between electron-rich polar head groups, with peak positions determined by the different alkyl tail lengths present in the malonamide molecule. This explanation of the SAXS correlation peaks contrasts with the prevailing literature, which attributes mesoscale features observed in small-angle scattering to the formation of microemulsions. Instead, this work finds that these features are present in the absence of water or any reverse micellar organization of the malonamides. As such, molecular-scale malonamide self-association and packing, rather than microemulsion-based colloidal-scale descriptions, is a more appropriate framework for these LLE systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Visualizing group II intron dynamics between the first and second steps of splicing

Group II introns are ubiquitous self-splicing ribozymes and retrotransposable elements evolutionarily and chemically related to the eukaryotic spliceosome, with potential applications as gene-editing tools. Recent biochemical and structural data have captured the intron in multiple conformations at different stages of catalysis. Here, we employ enzymatic assays, X-ray crystallography, and molecular simulations to resolve the spatiotemporal location and function of conformational changes occurring between the first and the second step of splicing. We show that the first residue of the highly-conserved catalytic triad is protonated upon 5’-splice-site scission, promoting a reversible structural rearrangement of the active site (toggling). Protonation and active site dynamics induced by the first step of splicing facilitate the progression to the second step. Our insights into the mechanism of group II intron splicing parallels functional data on the spliceosome, thus reinforcing the notion that these evolutionarily-related molecular machines share the same enzymatic strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Simulation of injector dynamics during steady inductive helicity injection current drive in the HIT-SI experiment

We present simulations of inductive helicity injection in the Helicity Injected Torus with Steady Inductive helicity injection (HIT-SI) device that treats the entire plasma volume in a single dynamic MHD model. A new fully 3D numerical tool, the PSI-center TETrahedral mesh code, was developed that provides the geometric flexibility required for this investigation. Implementation of a zero-β Hall MHD model using PSI-TET will be presented including formulation of a new self-consistent magnetic boundary condition for the wall of the HIT-SI device. Results from simulations of HIT-SI are presented focusing on injector dynamics that are investigated numerically for the first time. Asymmetries in the plasma loading between the two helicity injectors and progression of field reversal in each injector are observed. Analysis indicates cross-coupling between injectors through confinement volume structures. Injector impedance is found to scale with toroidal current at fixed density, consistent with experimental observation. Comparison to experimental data with an injector drive frequency of 14.5 kHz shows good agreement with magnetic diagnostics. Global mode structures from Bi-Orthogonal decomposition agree well with experimental data for the first four modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

First-principles theory of the Dirac semimetal Cd 3 As 2 under Zeeman magnetic field

Time-reversal broken Weyl semimetals have attracted much attention recently, but certain aspects of their behavior, including the evolution of their Fermi surface topology and anomalous Hall conductivity with Fermi-level position, have remained underexplored. A promising route to obtain such materials may be to start with a nonmagnetic Dirac semimetal and break time-reversal symmetry via magnetic doping or magnetic proximity. Here we explore this scenario in the case of the Dirac semimetal Cd 3 As 2 , based on first-principles density-functional calculations and subsequent low-energy modeling of Cd 3 As 2 in the presence of a Zeeman field applied along the symmetry axis. We clarify how each four-fold degenerate Dirac node splits into four Weyl nodes, two with chirality ±1 and two higher-order nodes with chirality ±2. Using a minimal k · p model Hamiltonian whose parameters are fit to the first-principles calculations, we detail the evolution of the Fermi surfaces and their Chern numbers as the Fermi energy is scanned across the region of the Weyl nodes at fixed Zeeman field. We also compute the intrinsic anomalous Hall conductivity as a function of Fermi-level position, finding a characteristic inverted-dome structure. Cd 3 As 2 is especially well suited to such a study because of its high mobility, but the qualitative behavior revealed here should be applicable to other Dirac semimetals as well.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Low-Cost Parallel Algorithms for 2:1 Octree Balance

The logical structure of a forest of octrees can be used to create scalable algorithms for parallel adaptive mesh refinement (AMR), which has recently been demonstrated for several petascale applications. Among various frequently used octree-based mesh operations, including refinement, coarsening, partitioning, and enumerating nodes, ensuring a 2:1 size balance between neighboring elements has historically been the most expensive in terms of CPU time and communication volume. The 2:1 balance operation is thus a primary target to optimize. One important component of a parallel balance algorithm is the ability to determine whether any two given octants have a consistent distance/size relation. Based on new logical concepts we propose fast algorithms for making this decision for all types of 2:1 balance conditions in 2D and 3D. Since we are able to achieve this without constructing any parent nodes in the tree that would otherwise need to be sorted and communicated, we can significantly reduce the required memory and communication volume. In addition, we propose a lightweight collective algorithm for reversing the asymmetric communication pattern induced by non-local octant interactions. We have implemented our improvements as part of the opensource “p4est” software. Benchmarking this code with both synthetic and simulation-driven adapted meshes we are able to demonstrate much reduced runtime and excellent weak and strong scalability. On our largest benchmark problem with 5.13 × 10 11 octants the new 2:1 balance algorithm executes in less than 8 seconds on 112,128 CPU cores of the Jaguar Cray XT5 supercomputer.

Computer Science↗
DOE OSTI2021

Embodied carbon in the U.S. manufacturing and trade

On average, one quarter of the global carbon footprint is embodied in traded goods. These emissions are a growing issue for global efforts to decarbonize the world economy. Embodied emissions in trade are not accounted for by existing greenhouse gas accounting systems. For example, countries only report their domestic carbon emissions (also known as production-based or territorial accounting) to the Intergovernmental Panel on Climate Change (IPCC). If the embodied carbon in trade were accounted for and reported, the promising climate trends depicted by many countries would be negated or reversed. For example, many achievements of reducing emissions by developed countries under the Kyoto Protocol would actually appear as emissions outsourced to developing countries. Like many other developed countries, the United States is a net importer of embodied greenhouse gas (GHG) emissions in trade. The goal of this study was to calculate the carbon footprint of various sectors of the U.S. economy using the latest available data and based on that, estimate the embodied carbon in certain manufacturing sector products imported and exported by the U.S. First, we conducted an input-output analysis to calculate the carbon footprint of 401 aggregate level sectors of the U.S. economy. After that, for a selected number of manufacturing subsectors /products(hereafter referred to as products),we collected their trade(import/export) data and combined that with their computed carbon footprint to estimate the embodied carbon in import and export of these products for the U.S. While we conducted the carbon footprint analysis for sectors covering the entire U.S. economy (401 sectors), we only conducted the embodied carbon in trade analysis for a few products within specific aggregate level manufacturing sectors,such as computer and electronic product manufacturing (NAICS 334), transport equipment manufacturing (NAICS 336), machinery manufacturing (NAICS 333), and chemical manufacturing (NAICS 325). It is possible to do the detailed embodied carbon analysis for all other products within the various aggregate level manufacturing sectors of the U.S. economy, but given the scope of this study and space constraint of this report, we only present the results for a few products. The results of our analysis for the selected products show the embodied carbon in trade for each product by country and also show the rankings of countries in terms of top importers and exporters of embodied carbon in trade with the U.S. Our results show that in majority of cases, the U.S. is a net importer of embodied carbon in trade. Unless consumption-based accounting is used, the U.S. and other net embodied carbon importing countries may continue to outsource their emissions to meet their climate change mitigation targets under the Paris Agreement, as observed previously with the Kyoto Protocol. Some countries reported reductions that exceeded their Kyoto targets, however, the changes in emissions embodied in imports were comparable to or larger than changes in domestic emissions. Traded emissions undermined emissions reductions in the Kyoto Protocol, and threaten to continue to do so for the Paris Agreement. New climate policies such as California’s Buy Clean act, EU’s proposed border carbon tax adjustment, as well as efforts in the private sector can help address the embodied carbon in traded products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Defining the S-Glutathionylation Proteome by Biochemical and Mass Spectrometric Approaches

Protein S-glutathionylation (SSG) is a reversible post-translational modification (PTM) featuring the conjugation of glutathione to a protein cysteine thiol. SSG can alter protein structure, activity, subcellular localization, and interaction with small molecules and other proteins. Thus, it plays a critical role in redox signaling and regulation in various physiological activities and pathological events. In this review, we summarize current biochemical and analytical approaches for characterizing SSG at both the proteome level and at individual protein levels. To illustrate the mechanism underlying SSG-mediated redox regulation, we highlight recent examples of functional and structural consequences of SSG modifications. Finally, we discuss the analytical challenges in characterizing SSG and the thiol PTM landscape, future directions for understanding of the role of SSG in redox signaling and regulation and its interplay with other PTMs, and the potential role of computational approaches to accelerate functional discovery.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

Adaptive multigrid domain decomposition solutions for viscous interacting flows

Several viscous incompressible flows with strong pressure interaction and/or axial flow reversal are considered with an adaptive multigrid domain decomposition procedure. Specific examples include the triple deck structure surrounding the trailing edge of a flat plate, the flow recirculation in a trough geometry, and the flow in a rearward facing step channel. For the latter case, there are multiple recirculation zones, of different character, for laminar and turbulent flow conditions. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable, non-staggered grid computation.

Rubin, Stanley G.↗
NASA NTRS1997

Computer-aided design of a proton pump

The use of transmembrane proton gradients in energy transduction is an almost universal feature of life on earth. These proton gradients are established and maintained by specialized assemblies of proteins which actively pump protons across membranes. One broad class of proton pumps uses captured light energy to drive the proton pumping. Our goal is to elucidate the minimum structural requirements of a light-driven proton-pump. There are two basic components to a simple light-driven proton pump: a source of photo-generated protons and a "gate-keeper" which prevents these protons from reattaching themselves to their source. A wide variety of molecules in the membrane, even as simple as polycyclic aromatic hydrocarbons, are capable of releasing protons when illuminated. Our work is therefore focused on the design of the "gate-keeper." Our initial model involves a pair of proton acceptors, coupled to each other by a transient water bridge, and supported in the membrane by a small bundle of peptide helices. Upon illumination, the proton source transfers its proton to the:- first acceptor of the gate-keeper. While the reverse reaction is highly probable, all that is needed to ensure irreversibility is a nonvanishing probability that the proton will be transferred to the second acceptor across a transient water bridge. Back transfer of the proton to the first acceptor, and thence to the proton source, is impeded by the free energy required to move the proton uphill towards the. proton source and by the disruption of the transient water bridge. As a prototypical water-bridged proton transfer system, we are studying the transfer of a proton across a water bridge from a formic acid to a formate anion. With a pK(sub alpha), of 3.7. formic acid is a good model for the acidic amino acids glutamate and aspartate which are good candidates for gate-keeper proton acceptors. Simulations of proton transfer reactions in a membrane are complicated by the quantum mechanical nature of the process breaking and forming chemical bonds.

New, Michael H.↗
DOE OSTI2011

Electron Bernstein Wave Studies in MST

The overdense plasma in an RFP prevents electromagnetic waves from propagating past the edge, however use of the electron Bernstein wave (EBW) has the potential to heat and drive current in the plasma. MHD simulations have demonstrated that resistive tearing mode stability is very sensitive to gradients in the edge current density profile allowing EBW to potentially be a stabilizing influence. A new MW level experiment is being commissioned on MST to evaluate the potential use of the EBW for current profile control on the RFP. The development of new equipment includes a 5.5 GHz klystron driven by a novel switchmode power supply. A quartz window has been constructed and coupling with a cylindrical molybdenum wave guide antenna has been studied. Due to the steep edge density gradient in the RFP, it is possible to efficiently couple to the EBW with O or X mode launch. The EBW is strongly damped at the electron cyclotron resonance where it couples to the electron gyromotion and alters the electron distribution. Either Fisch-Boozer or Ohkawa current drive mechanisms can be activated to drive off axis current in the plasma. Preliminary experiments have been performed to verify high power coupling and understand heating via observed x-ray emission when compared to Fokker-Plank modeling in CQL3D.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Role of resistivity and viscosity in the excitation of stable m = 0 modes during the RFP sawtooth crash

Visco-resistive magnetohydrodynamic simulations are used herein to investigate the role of resistivity and viscosity in the reversed field pinch sawtooth crash. Specifically, we examine the temporal behavior of the linearly stable ($m, n$) = (0, 1) mode. Both volume-averaged dissipation and dissipation levels in the region near the $m$ = 0 resonance are found to influence characteristic time scales. Increasing resistivity causes the mode rise time, fall time, and total crash duration to decrease, while increasing viscosity has the opposite effect. Examination of energy flow during the sawtooth crash reveals that despite the linear stability of the (0, 1) mode, it receives most of its energy from the mean current profile during the crash rather than from nonlinear interactions with other modes. Resistivity and viscosity do not impact mode behavior directly through dissipative energy loss but rather through modification of the large scale current profile evolution and radial structure of the mode. Computational results are compared to experimental data from the Madison Symmetric Torus and found to largely agree when resistivity and viscosity are similar.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

MO-A-BRD-07: Feasibility of X-Ray Acoustic Computed Tomography as a Tool for Calibration and In Vivo Dosimetry of Radiotherapy Electron and Photon Beams

Purpose: This work simulates radiation-induced acoustic waves to assess the feasibility of x-ray acoustic computed tomography (XACT) as a dosimeter. XACT exploits the phenomenon that acoustic waves with amplitude proportional to the dose deposited are induced following a radiation pulse. After detecting these acoustic waves with an ultrasound transducer, an image of the dose distribution can be reconstructed in realtime. Methods: Monte Carlo was used to simulate the dose distribution for monoenergetic 6 MeV photon and 9 MeV electron beams incident on a water tank. The dose distribution for a prostate patient planned with a photon 4-field box technique was calculated using clinical treatment planning software. All three dose distributions were converted into initial pressure distributions, and transportation of the induced acoustic waves was simulated using an open-source toolkit. Ideal transducers were placed around the circumference of the target to detect the acoustic waves, and a time reversal reconstruction algorithm was used to obtain an XACT image of the dose for each radiation pulse. Results: For the photon water tank relative dosimetry case, it was found that the normalized acoustic signal amplitude agreed with the normalized dose at depths from 0 cm to 10 cm, with an average percent difference of 0.5%. For the reconstructed in-plane dose distribution of an electron water tank irradiation, all pixels passed a 3%–3 mm 2D gamma test. The reconstructed prostate dose distribution closely resembled the plan, with 89% of pixels passing a 3%–3 mm 2D gamma test. For all situations, the amplitude of the induced acoustic waves ranged from 0.01 Pa to 1 Pa. Conclusion: Based on the amplitude of the radiation-induced acoustic waves and accuracy of the reconstructed dose distributions, XACT is a feasible technique for dosimetry in both calibration and in vivo environments for photon and electron beams and merits further investigation. Funding from NSERC, CIHR and McGill University. SH acknowledges partial support by the CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council (Grant number: 432290)

60 APPLIED LIFE SCIENCES↗