Search NASA⌕ Search

SEARCH · Search NASA

Results for “reach code”

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 721 records · Page 40

CASSIS: THE CORNELL ATLAS OF SPITZER/INFRARED SPECTROGRAPH SOURCES

We present the spectral atlas of sources observed in low resolution with the Infrared Spectrograph on board the Spitzer Space Telescope. More than 11,000 distinct sources were extracted using a dedicated algorithm based on the SMART software with an optimal extraction (AdOpt package). These correspond to all 13,000 low-resolution observations of fixed objects (both single source and cluster observations). The pipeline includes image cleaning, individual exposure combination, and background subtraction. Particular attention is given to bad pixel and outlier rejection at the image and spectra levels. Most sources are spatially unresolved so that optimal extraction reaches the highest possible signal-to-noise ratio. For all sources, an alternative extraction is also provided that accounts for all of the source flux within the aperture. CASSIS provides publishable quality spectra through an online database together with several important diagnostics, such as the source spatial extent and a quantitative measure of detection level. Ancillary data such as available spectroscopic redshifts are also provided. The database interface will eventually provide various ways to interact with the spectra, such as on-the-fly measurements of spectral features or comparisons among spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Determination of transit dose profile for a {sup 192}Ir HDR source

Purpose: Several studies have reported methodologies to calculate and correct the transit dose component of the moving radiation source for high dose rate (HDR) brachytherapy planning systems. However, most of these works employ the average source speed, which varies significantly with the measurement technique used, and does not represent a realistic speed profile, therefore, providing an inaccurate dose determination. In this work, the authors quantified the transit dose component of a HDR unit based on the measurement of the instantaneous source speed to produce more accurate dose values. Methods: The Nucletron microSelectron-HDR Ir-192 source was characterized considering the Task Group 43 (TG-43U1) specifications. The transit dose component was considered through the calculation of the dose distribution using a Monte Carlo particle transport code, MCNP5, for each source position and correcting it by the source speed. The instantaneous source speed measurements were performed in a previous work using two optical fibers connected to a photomultiplier and an oscilloscope. Calculated doses were validated by comparing relative dose profiles with those obtained experimentally using radiochromic films. Results: TG-43U1 source parameters were calculated to validate the Monte Carlo simulations. These agreed with the literature, with differences below 1% for the majority of the points. Calculated dose profiles without transit dose were also validated by comparison with ONCENTRA{sup Registered-Sign} Brachy v. 3.3 dose values, yielding differences within 1.5%. Dose profiles obtained with MCNP5 corrected using the instantaneous source speed profile showed differences near dwell positions of up to 800% in comparison to values corrected using the average source speed, but they are in good agreement with the experimental data, showing a maximum discrepancy of approximately 3% of the maximum dose. Near a dwell position the transit dose is about 22% of the dwell dose delivered by the source dwelling 1 s and reached 104.0 cGy per irradiation in a hypothetical clinical case studied in this work. Conclusions: The present work demonstrated that the transit dose correction based on average source speed fails to accurately correct the dose, indicating that the correct speed profile should be considered. The impact on total dose due to the transit dose correction near the dwell positions is significant and should be considered more carefully in treatments with high dose rate, several catheters, multiple dwell positions, small dwell times, and several fractions.

60 APPLIED LIFE SCIENCES↗

Modeling of small tungsten dust grains in EAST tokamak with NDS-BOUT ++

In order to investigate the transport of small dusts as well as their evolution property along their trajectories, the NDS module is developed under the BOUT++ framework, a highly desirable C++ code package to perform parallel plasma fluid simulations with an arbitrary number of equations in three-dimensional curvilinear coordinates. Due to the severe dust ablation in fusion plasmas, the dust size would decrease from micrometer to nanometer, resulting in impurities. Small dusts in the simulations here are specified as tungsten spheres with the radii on or below the order of submicrometer. The Rayleigh limit is included in the charging process when the dust is ablated to the droplet phase. The simulation results from the NDS module show that a 200 nm radius spherical tungsten dust originated from upper divertor region of EAST Tokamak is ablated completely due to the intense heating from the incoming plasma inside the core region, well consistent with the CCD footage of EAST shot # 81459. Furthermore it is found that the magnetic field dominates the dust transport when the dust radius is below 100 nm during the ablation along the trajectory. Our simulations predict that a 10 nm radius spherical tungsten dust injected from the inner midplane is well constrained by the magnetic field, and it reaches the inner divertor target with a velocity on the order of km/s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory comparison and numerical benchmarking on neoclassical toroidal viscosity torque

Systematic comparison and numerical benchmarking have been successfully carried out among three different approaches of neoclassical toroidal viscosity (NTV) theory and the corresponding codes: IPEC-PENT is developed based on the combined NTV theory but without geometric simplifications [Park et al., Phys. Rev. Lett. 102, 065002 (2009)]; MARS-Q includes smoothly connected NTV formula [Shaing et al., Nucl. Fusion 50, 025022 (2010)] based on Shaing's analytic formulation in various collisionality regimes; MARS-K, originally computing the drift kinetic energy, is upgraded to compute the NTV torque based on the equivalence between drift kinetic energy and NTV torque [J.-K. Park, Phys. Plasma 18, 110702 (2011)]. The derivation and numerical results both indicate that the imaginary part of drift kinetic energy computed by MARS-K is equivalent to the NTV torque in IPEC-PENT. In the benchmark of precession resonance between MARS-Q and MARS-K/IPEC-PENT, the agreement and correlation between the connected NTV formula and the combined NTV theory in different collisionality regimes are shown for the first time. Additionally, both IPEC-PENT and MARS-K indicate the importance of the bounce harmonic resonance which can greatly enhance the NTV torque when E×B drift frequency reaches the bounce resonance condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Community-driven development for computational biology at Sprints, Hackathons and Codefests

Background: Computational biology comprises a wide range of technologies and approaches. Multiple technologies can be combined to create more powerful workflows if the individuals contributing the data or providing tools for its interpretation can find mutual understanding and consensus. Much conversation and joint investigation are required in order to identify and implement the best approaches. Traditionally, scientific conferences feature talks presenting novel technologies or insights, followed up by informal discussions during coffee breaks. In multi-institution collaborations, in order to reach agreement on implementation details or to transfer deeper insights in a technology and practical skills, a representative of one group typically visits the other. However, this does not scale well when the number of technologies or research groups is large. Conferences have responded to this issue by introducing Birds-of-a-Feather (BoF) sessions, which offer an opportunity for individuals with common interests to intensify their interaction. However, parallel BoF sessions often make it hard for participants to join multiple BoFs and find common ground between the different technologies, and BoFs are generally too short to allow time for participants to program together. Results: This report summarises our experience with computational biology Codefests, Hackathons and Sprints, which are interactive developer meetings. They are structured to reduce the limitations of traditional scientific meetings described above by strengthening the interaction among peers and letting the participants determine the schedule and topics. These meetings are commonly run as loosely scheduled “unconferences” (self-organized identification of participants and topics for meetings) over at least two days, with early introductory talks to welcome and organize contributors, followed by intensive collaborative coding sessions. We summarise some prominent achievements of those meetings and describe differences in how these are organised, how their audience is addressed, and their outreach to their respective communities. Conclusions: Hackathons, Codefests and Sprints share a stimulating atmosphere that encourages participants to jointly brainstorm and tackle problems of shared interest in a self-driven proactive environment, as well as providing an opportunity for new participants to get involved in collaborative projects.

59 BASIC BIOLOGICAL SCIENCES↗

OECD MCCI Project: Category 4 Integral Test to Validate Severe Accident Codes (Rev. 1, Final Report)

Ex-vessel debris coolability is an important light water reactor (LWR) technical challenge. For existing plants, resolution of this challenge will confirm the technical basis for severe accident management guidelines (SAMGs). For new reactors, understanding this challenge will help confirm the effectiveness of the design and implementation of new accident mitigation features and severe accident management design alternatives (SAMDAs). The first OECD-MCCI program conducted reactor material experiments focused on achieving the following technical objectives: i) provide confirmatory evidence and data for various cooling mechanisms through separate effect tests for severe accident model development, and ii) provide long-term 2-D coreconcrete interaction data for code assessment and improvement. In response to the PRG mandate, a summary test plan for this sixth Core-Concrete Interaction (CCI-6) experiment was prepared reviewed by the PRG, and subsequently approved for test conduct. Detailed specifications that were developed as part of the test planning process are provided. The experiment approach was to incorporate an array of water nozzles into the core-concrete interaction test apparatus. The test section was enlarged to incorporate a 70 cm x 70 cm basemat. The initial melt depth was increased to 28 cm,1 which corresponds to a thermite charge mass of 900 kg. The melt depth is noted to be slightly less than the 30 cm level requested by the PRG. The thermite charge was redesigned to include a reduced amount of concrete decomposition products (i.e., 6 wt %) relative to the previously performed CCI tests. The concrete sidewalls allow for 24 cm of lateral erosion, and the sidewall design was modified to include refractory MgO with a UO 2 pellet and U 3 O 8 powder liner to stop the ablation if the 24 cm ablation depth was reached. The basemat and sidewalls also included an array of 45 concrete nozzles that were arranged on a 14 cm lateral pitch and were located at a uniform depth of 27.5 cm below the initial concrete basemat surface elevation. The nozzles were originally sealed at the top, and were fed with water from a common water supply tank that provided a static water head of nominally 10 kPa measured at the upper surface of the nozzles. The CCI-6 experiment was carried out on 18 March 2010. Following a description of related studies in the next section, the balance of this final report provides a summary description of the test facility as well a presentation and discussion of results from the experiment.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Rapid Sample Insertion/Extraction Systems for Gamma Irradiation

The Rapid Sample Insertion/Extraction System for Gamma Irradiation, otherwise known as the "rabbit" system, was a four-week long project which included many different aspects; from coding an Arduino to building PVC piping to 3-D printing the "rabbit" capsules. The "rabbit" system is a system of PVC piping that allows a quick and efficient transfer of materials into/out of one of the irradiation chambers in the Gamma Irradiation Facility (GIF) with the use of a 3-D printed "rabbit." This "rabbit" encapsulates material to be irradiated and carries it from a position outside of the irradiation chamber to the basket inside of the chamber. The main purpose of this system is to save time and provide more exact data without any delays that normally occur when a person has to enter the chamber, retrieve data, and then analyze the data. This system should take measurements and retrieve the data instantaneously. The way in which the "rabbit" is sent through the PVC piping is with an advanced bi-directional, high-throughput pneumatic system, or a shop vacuum cleaner. When the vacuum is set to blow or suck then the "rabbit" will be pulled through the PVC piping to its intended destination and will hit sensors along the sides of the tubing when it reaches the end of the piping. These sensors tell the Arduino that the "rabbit" is finished moving throughout the tubing and stops a timer. Another timer is used to see how long the "rabbit" is being irradiated so when the "rabbit" reaches the sensors in the basket in the irradiation chamber another timer is started and it ends when the sensors no longer detect the "rabbit," which means that it has begun its journey back to the starting point.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DIRECT NUMERICAL SIMULATION OF RADIATION PRESSURE-DRIVEN TURBULENCE AND WINDS IN STAR CLUSTERS AND GALACTIC DISKS

The pressure exerted by the radiation of young stars may be an important feedback mechanism that drives turbulence and winds in forming star clusters and the disks of starburst galaxies. However, there is great uncertainty in how efficiently radiation couples to matter in these high optical depth environments. In particular, it is unclear what levels of turbulence the radiation can produce, and whether the infrared radiation trapped by the dust opacity can give rise to heavily mass-loaded winds. In this paper, we report a series of two-dimensional flux-limited diffusion radiation-hydrodynamics calculations performed with the code ORION in which we drive strong radiation fluxes through columns of dusty matter confined by gravity in order to answer these questions. We consider both systems where the radiation flux is sub-Eddington throughout the gas column, and those where it is super-Eddington at the midplane but sub-Eddington in the atmosphere. In the latter, we find that the radiation-matter interaction gives rise to radiation-driven Rayleigh-Taylor instability, which drives supersonic turbulence at a level sufficient to fully explain the turbulence seen in Galactic protocluster gas clouds, and to make a non-trivial contribution to the turbulence observed in starburst galaxy disks. However, the instability also produces a channel structure in which the radiation-matter interaction is reduced compared to time-steady analytic models because the radiation field is not fully trapped. For astrophysical parameters relevant to forming star clusters and starburst galaxies, we find that this effect reduces the net momentum deposition rate in the dusty gas by a factor of {approx}2-6 compared to simple analytic estimates, and that in steady state the Eddington ratio reaches unity and there are no strong winds. We provide an approximation formula, appropriate for implementation in analytic models and non-radiative simulations, for the force exerted by the infrared radiation field in this regime.

79 ASTRONOMY AND ASTROPHYSICS↗

Rapid Sample Insertion/Extraction Systems for Gamma Irradiation

The Rapid Sample Insertion/Extraction System for Gamma Irradiation, otherwise known as the "rabbit" system, was a four-week long project which included many different aspects such as coding an Arduino, building PVC piping, and 3-D printing the "rabbit" capsules. The "rabbit" system is a system of PVC piping that allows a quick and efficient transfer of materials into and out of one of the irradiation chambers in the Gamma Irradiation Facility (GIF) with the use of a 3-D printed "rabbit." This "rabbit" encapsulates material to be irradiated and carries it from a position outside of the irradiation chamber to the basket inside of the chamber. The main purpose of this system is to save time and provide more exact data without any delays that normally occur when a person has to enter the chamber, retrieve data, and then analyze the data. This system should take measurements and retrieve the data instantaneously. The way in which the "rabbit" is sent through the PVC piping is with an advanced bi-directional, high-throughput pneumatic system, or a shop vacuum cleaner. When the vacuum is set to blow or suck then the "rabbit" will be pulled or pushed through the PVC piping to its intended destination and will hit sensors along the sides of the tubing when it reaches the end of the piping. These sensors tell the Arduino that the "rabbit" is finished moving throughout the tubing and stops a timer. Another timer is used to see how long the "rabbit" is being irradiated so when the "rabbit" reaches the sensors in the basket in the irradiation chamber another timer is started and it ends when the sensors no longer detect the "rabbit," which means that it has begun its journey back to the starting point. These times, as well as the temperatures, are the data necessary for the project and must be extremely accurate.

61 RADIATION PROTECTION AND DOSIMETRY↗

Model Validation and Uncertainty Quantification on the KRUSTY Microreactor Design Using GRIFFIN Neutron Transport Code [Poster]

Argonne National Laboratory (ANL) and INL have developed a GRIFFIN steady state neutronics model for the multiphysics simulations of the Kilopower Reactor Using Sterling TechnologY (KRUSTY) microreactor in the Multiphysics Object Oriented Simulation Environment (MOOSE). The reliability of such deterministic neutronics models can be validated by comparing with computations from Monte Carlo codes (e.g. MCNP, SERPENT, OpenMC, Shift, etc). Furthermore, potential modeling/design improvements can be identified by incorporating uncertainty quantification (UQ), which can be performed by MOOSE’s Stochastic Tools Module (STM). KRUSTY is a prototype for a 5-kW thermal nuclear-powered space reactor. Its primary components consist of nuclear fuel, heat pipes, a control rod, a reflector, and the shielding. The fuel consists of 3 stacked U-7.65Mo cylinders with a hole in the center for the control rod. 8 liquid sodium heat pipes transfer fission energy from the solid fuel block to the Sterling power conversion system where the energy is extracted, and the cooled sodium flows back to the core via capillary action . The movable Boron Carbide control rod regulates the neutron population during startup or when a reactor temperature boost is needed . The beryllium oxide reflector is in 3 places in the reactor; it surrounds the core axially, it lies beneath the core on a platen, and it is present in the shim. The axial and lower reflectors rest on an adjustable stainless-steel platen that moves upward to cover the fuel and help the reactor reach criticality. Lastly, radial stainless steel surrounds the core offering protection from radiation exposure .

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Data for: Time to anoxia: Observations and predictions of oxygen drawdown following coastal flood events

This dataset contains data and code associated with the paper Patel et al. "Time to anoxia: Observations and predictions of oxygen drawdown following coastal flood events".The coastal terrestrial-aquatic interface (TAI) is a highly dynamic system characterized by strong physical, chemical, and biological gradients. In particular, shifting soil redox conditions and consumption of terminal electron acceptors, due in part to dynamic hydrologic conditions, is a strong driver of carbon availability and transformations across TAIs. However, while redox dynamics are well described, our ability to quantitatively forecast rates of oxic to anoxic shifts in soils with different characteristics and inundation regimes is limited. We integrated field measurements, laboratory incubations, and model simulations to improve mechanistic understanding of oxygen consumption dynamics in coastal soils. Continuous in situ monitoring unexpectedly revealed that flooding caused temporary spikes in subsurface dissolved oxygen followed by rapid consumption in the wetlands. To further investigate these mechanisms in a controlled setting, we performed laboratory incubations using surface and subsurface soils from a TAI gradient (defined here as upland forest to transitional forest to wetland) in Western Lake Erie to measure oxygen consumption rates in TAI soils during flood events. In our experiments, wetland soils reached anoxia the fastest, in ∼ 9 h on average, whereas upland soils turned anoxic in ∼ 18 h. Subsurface upland soils did not turn anoxic even after two weeks of saturation in the lab, and their oxygen consumption patterns suggested carbon and/or nutrient limitation. These results are consistent with in-situ groundwater redox and oxygen measurements in the field, where wetland soils exhibited the highest rates of oxygen consumption along the TAI. Model simulations of oxygen consumption suggested that oxygen consumption had stronger abiotic controls in wetland soils but stronger biotic controls in upland soils, providing a useful framework for future incubation experiments. Microbial activity is a strong driver of oxygen consumption in TAI soils, although it is constrained by the availability of dissolved carbon in subsurface soils.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/COMPASS-DOE/cmps-time_to_anoxia (DOI: 10.5281/zenodo.10815374).

54 ENVIRONMENTAL SCIENCES↗

Radiative properties of astrophysical matter : a quest to reproduce astrophysical conditions on earth.

Experiments in terrestrial laboratories can be used to evaluate the physical models that interpret astronomical observations. The properties of matter in astrophysical objects are essential components of these models, but terrestrial laboratories struggle to reproduce the extreme conditions that often exist. Megajoule-class DOE/NNSA facilities such as the National Ignition Facility and Z can create unprecedented amounts of matter at extreme conditions, providing new capabilities to test astrophysical models with high accuracy. Experiments at these large facilities are challenging, and access is very competitive. However, the cylindrically-symmetric Z source emits radiation in all directions, enabling multiple physics experiments to be driven with a single Z discharge. This helps ameliorate access limitations. This article describes research efforts under way at Sandia National Laboratories Z facility investigating radiation transport through stellar interior matter, population kinetics of atoms exposed to the intense radiation emitted by accretion powered objects, and spectral line formation in white dwarf (WD) photospheres. Opacity quantifies the absorption of radiation by matter and strongly influences stellar structure and evolution, since radiation dominates energy transport deep inside stars. Opacity models have become highly sophisticated, but laboratory tests at the conditions existing inside stars have not been possible - until now. Z research is presently focused on measuring iron absorption at conditions relevant to the base of the solar convection zone, where the electron temperature and density are 190 eV and 9 x 10{sup 22} e/cc, respectively. Creating these conditions in a sample that is sufficiently large, long-lived, and uniform is extraordinarily challenging. A source of radiation that streams through the relatively-large samples can produce volumetric heating and thus, uniform conditions, but to achieve high temperatures a strong source is required. Z dynamic hohlraums provide such a megajoule-class source. Initial Z experiments measured transmission through iron samples ionized to the same charge states that exist at the solar convection zone base. The resulting data made it possible to test challenging aspects of the opacity calculations such as the ionization balance and the completeness and accuracy of the atomic energy level description. However, the density was too low to provide a definitive test of the physics at the solar convection zone base. Recent experiments have reached higher densities, and opacity model tests for stellar interiors now appear within reach. Accretion powered objects, including active galactic nuclei, x-ray binaries, and black hole accretion disks, are the most luminous objects in the universe. Astrophysical models for these objects rely largely on comparing spectroscopic predictions with observations. A dilemma arises because the spectra originate from plasmas that are bathed in the enormous photon flux from the accretion disk and photoionization dominates the atomic ionization and energy level populations. Thus, constraining astrophysical models depends on accurate atomic models for photoionized plasmas. Unfortunately, to date the ionization in almost all laboratory experiments is collision-dominated and very few tests of photoionized plasma atomic kinetics exist. Megajoule class high-energy-density facilities can help because they generate higher x-ray fluence over larger spatial scales and longer times. Expanded iron foils and pre-filled neon gas cells have been used in experiments at Z to study photoionized atomic kinetics in two elements commonly observed in astrophysical objects. In these experiments, low density samples are exposed to a measured intense x-ray spectrum, emergent emission or absorption spectra are recorded, and the results are compared to predictions made with spectral synthesis codes used by astrophysicists. Initial experiments focused on testing models used to interpret spectra from the 'warm absorber,' a plasma observed in the vicinity of some active galactic nuclei. In future experiments, spectra from plasmas exposed to higher radiation intensities will be measured, possibly leading to improved understanding of the plasma in the immediate vicinity of the accretion disk. WDs are the oldest stars and can serve as cosmic clocks, since the universe must be at least as old as the objects within it. Astrophysicists determine WD ages using stellar models combined with effective temperature (T{sub eff}) and mass inferred from spectral observations of WD photospheres. Many line profiles in the observed spectra are dominated by Stark broadening, a process sensitive to the photosphere density and related to the total mass through the stellar model. Accurate Stark broadening theory is, therefore, critical to the precise determination of the WD properties and the inferred ages.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A KINE-CHEMICAL INVESTIGATION OF THE AB DOR MOVING GROUP 'STREAM'

The AB Dor Moving Group consists of a 'nucleus' of {approx}10 stars at d {approx_equal} 20 pc, along with dozens of purported 'stream' members distributed across the sky. We perform a chemical and kinematic analysis of a subsample of AB Dor stream stars to test whether they constitute a physical stellar group. We use the NEMO Galactic kinematic code to investigate the orbits of the stream members, and perform a chemical abundance analysis using high resolution spectra taken with the Magellan Clay 6.5 m telescope. Using a {chi}{sup 2} test with the measured abundances for 10 different elements, we find that only half of the purported AB Dor stream members could possibly constitute a statistically chemically homogeneous sample. Some stream members with three-dimensional velocities were hundreds of parsecs from the AB Dor nucleus {approx}10{sup 8} yr ago, and hence were unlikely to share a common origin. We conclude that the published lists of AB Dor moving group stream members are unlikely to represent the dispersed remnant of a single star formation episode. A subsample of the stream stars appears to be both statistically chemically homogeneous and in the vicinity of the AB Dor nucleus at birth. Their mean metallicity is [Fe/H] = 0.02 {+-} 0.02 dex, which we consider representative for the AB Dor group. Finally, we report a strong lower limit on the age of the AB Dor nucleus of >110 Myr based on the pre-main sequence contraction times for K-type members which have reached the main sequence.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct-drive double-shell implosion: A platform for burning-plasma physics studies

Double-shell ignition designs have been investigated with the indirect-drive inertial confinement fusion (ICF) scheme in both simulations and experiments in which the inner-shell kinetic energy was limited to ~10 to 15 kJ, even driven by megajoule-class lasers such as the National Ignition Facility. Since direct-drive ICF can couple more energy to the imploding shells, we have performed a detailed study on direct-drive double-shell (D 3 S) implosions with state-of-the-art physics models implemented in radiation-hydrodynamic codes (LILAC and DRACO), including nonlocal thermal transport, cross-beam energy transfer (CBET), and first-principles-based material properties. To mitigate classical unstable interfaces, we have proposed the use of a tungsten/beryllium–mixed inner shell with gradient-density layers that can be made by magnetron sputtering. In our D 3 S designs, a 70-μm-thick beryllium outer shell is driven symmetrically by a high-adiabat (α ≥ 10), 1.9-MJ laser pulse to a peak velocity of ~240 km/s. Upon spherical impact, the outer shell transfers ~30 to 40 kJ of kinetic energy to the inner shell filled with deuterium–tritium gas or liquid, giving neutron-yield energies of ~6 MJ in 1-D simulations. Two-dimensional high-mode DRACO simulations indicated that such high-adiabat D 3 S implosions are not susceptible to laser imprint, but the long-wavelength perturbations from the laser port configuration along with CBET can be detrimental to the target performance. Yet, neutron yields of ~0.3- to 1.0-MJ energies can still be obtained from our high-mode DRACO simulations. The robust α-particle bootstrap is readily reached, which could provide a viable platform for burning-plasma physics studies. Once CBET mitigation and/or more laser energy becomes available, we believe that breakeven or moderate energy gain might be feasible with the proposed D 3 S scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aspect ratio effects on neoclassical tearing modes from comparison between DIII-D and National Spherical Torus Experiment

Neoclassical tearing mode islands are sustained by helically perturbed bootstrap currents arising at finite beta from toroidal effects that trap a fraction of the particles in non-circulating orbits. DIII-D and NSTX are here operated with similar shape and cross-sectional area but almost a factor of two difference in inverse aspect ratio a/R. In these experiments, destabilized n=1 tearing modes were self-stabilized (reached the 'marginal point') by reducing neutral-beam power and thus beta. The measure of the marginal island gives information on the small-island stabilizing physics that in part (with seeding) governs onset. The marginal island width on NSTX is found to be about three times the ion banana width and agrees with that measured in DIII-D, except for DIII-D modes closer to the magnetic axis, which are about two times the ion banana width. There is a balance of the helically perturbed bootstrap term with small island effects with the sum of the classical and curvature terms in the modified Rutherford equation for tearing-mode stability at the experimental marginal point. Empirical evaluation of this sum indicates that while the stabilizing effect of the curvature term is negligible in DIII-D, it is important in NSTX. The mode temporal behavior from the start of neutral-beam injection reduction also suggests that NSTX operates closer to marginal classical tearing stability; this explains why there is little hysteresis in beta between mode onset, saturation, and self-stabilization (while DIII-D has large hysteresis in beta). NIMROD code module component calculations based on DIII-D and NSTX reconstructed experimental equilibria are used to diagnose and confirm the relative importance of the stabilizing curvature effect, an advantage for low aspect ratio; the relatively greater curvature effect makes for less susceptibility to NTM onset even if the classical tearing stability index is near marginal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Neutrino Distributions in Phase Space for the Rotating Core-collapse Supernova Simulated with a Boltzmann-neutrino-radiation-hydrodynamics Code

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, in this work we simulate the collapse of a rotating core of the progenitor with a zero-age- main-sequence mass of 11.2 M ⊙ and a shellular rotation of $1\,\mathrm{rad}\,{{\rm{s}}}^{-1}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches ~20% in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.

79 ASTRONOMY AND ASTROPHYSICS↗

X-RAY SPECTRA FROM MAGNETOHYDRODYNAMIC SIMULATIONS OF ACCRETING BLACK HOLES

We present the results of a new global radiation transport code coupled to a general relativistic magnetohydrodynamic simulation of an accreting, non-rotating black hole. For the first time, we are able to explain from first principles in a self-consistent way all the components seen in the X-ray spectra of stellar-mass black holes, including a thermal peak and all the features associated with strong hard X-ray emission: a power law extending to high energies, a Compton reflection hump, and a broad iron line. Varying only the mass accretion rate, we are able to reproduce a wide range of X-ray states seen in most galactic black hole sources. The temperature in the corona is T{sub e} {approx} 10 keV in a boundary layer near the disk and rises smoothly to T{sub e} {approx}> 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to Almost-Equal-To 6M as the accretion rate decreases, we find that the shape of the Fe K{alpha} line is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

79 ASTRONOMY AND ASTROPHYSICS↗

Partnership for Edge Physics Simulation

A major goal of our participation in the Edge Physics Simulation project has been to contribute to the understanding of the self-organization of tokamak turbulence fluctuations resulting in the formation of a staircase structure in the ion temperature. A second important goal is to demonstrate how small scale turbulence in plasmas self-organizes with dynamically driven quasi-stationary flow shear. These goals have been accomplished through the analyses of the statistical properties of XGC1 flux driven Gyrokinetic electrostatic ion temperature gradient (ITG) turbulence simulation data in which neutrals are included. The ITG turbulence data, and in particular fluctuation data, were obtained from a massively parallel flux-driven gyrokinetic full-f particle-in-cell simulation of a DIII-D like equilibrium. Below some the findings are summarized. It was observed that the emergence of staircase structure is related to the variations in the normalized temperature gradient length (R/LT) and the poloidal flow shear. Average turbulence intensity is found to be large in the vicinity of minima in R/LTi, where ITG growth is expected to be lower. The distributions of the occurrences of potential fluctuation are found to be Gaussian away from the staircase-step locations, but they are found to be non-Gaussian in the vicinity of staircase-step locations. The results of analytically derived expressions for the distribution of the occurrences of turbulence intensity and intensity flux were compared with the corresponding quantities computed using XGC1 simulation data and good agreement is found. The derived expressions predicts inward and outward propagation of turbulence intensity flux in an intermittent fashion. The outward propagation of turbulence intensity flux occurs at staircase-step locations and is related to the change in poloidal flow velocity shear and to the change in the ion temperature gradient. The standard deviation, skewness and kurtosis for turbulence quantities were computed and found to be large in the vicinity of the staircase-step structures. Large values of skewness and kurtosis can be explained by a temporary opening and closing of the structure which allows turbulence intensity events to propagate. The staircase patterns may reduce the ion heat transport and a manipulation of these patterns may be used to optimize heat transport in tokamaks. An additional objective of the research in support of the Edge Physics Simulation initiative has been to improve the understanding of scrape-off layer thermal transport. In planning experiments and designing future tokamaks, it is important to understand the physical effects that contribute to divertor heat-load fluxes. The research accomplished will contribute to developing new models for the scrape-off layer region. The XGC0 code was used to compute the heat fluxes and the heat-load width in the outer divertor plates of C-Mod and DIII-D tokamaks. It was observed that the width of the XGC0 neoclassical heat-load was approximately inversely proportional to the total plasma current. Anomalous transport in the H-mode pedestal region of five Alcator C-Mod discharges, representing a collisionality scan, was analyzed. The understanding of anomalous transport in the pedestal region is important for the development of a comprehensive model for the H-mode pedestal slope. It was found that the electron thermal anomalous diffusivities at the pedestal top increase with the electron collisionality. This dependence can point to the DRIBM as the modes that drive the anomalous transport in the plasma edge of highly collisional discharges. The effects of plasma shaping on the H-mode pedestal structure was also investigated. The differences in the predicted H-mode pedestal width and height for the DIII-D discharges with different elongation and triangularities were discussed. For the discharges with higher elongation, it was found that the gradients of the plasma profiles in the H-mode pedestal reach semi-steady states. In these simulations, the pedestal slowly continued to evolve to higher pedestal pressures and bootstrap currents until the peeling ballooning stability conditions were satisfied. The discharges with lower elongation do not reach the semi-steady state, and ELM crashes were triggered at earlier times. The plasma elongation was found to have a stronger stabilizing effect than the plasma triangularity. For the discharges with lower elongation and lower triangularity, the ELM frequency was large, and the H-mode pedestal evolves rapidly. It was found that the temperature of neutrals in the scrape-off-layer region can affect the dynamics of the H-mode pedestal buildup. However, the final pedestal profiles were nearly independent of the neutral temperature. The elongation and triangularity affected the pedestal widths of plasma density and electron temperature profiles differently. This study illustrated a new mechanism for controlling the pedestal bootstrap current and the pedestal stability.

97 MATHEMATICS AND COMPUTING↗