Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

THE TRANSITION REGION RESPONSE TO A CORONAL NANOFLARE: FORWARD MODELING AND OBSERVATIONS IN SDO/AIA

The corona and transition region (TR) are fundamentally coupled through the processes of thermal conduction and mass exchange. It is not possible to understand one without the other. Yet the temperature-dependent emissions from the two locations behave quite differently in the aftermath of an impulsive heating event such as a coronal nanoflare. Whereas the corona cools sequentially, emitting first at higher temperatures and then at lower temperatures, the TR is multithermal and the emission at all temperatures responds in unison. We have previously applied the automated time lag technique of Viall and Klimchuk to disk observations of an active region (AR) made by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory. Lines of sight passing through coronal plasma show clear evidence for post-nanoflare cooling, while lines of sight intersecting the TR footpoints of coronal strands show zero time lag. In this paper, we use the EBTEL hydrodynamics code to demonstrate that this is precisely the expected behavior when the corona is heated by nanoflares. We also apply the time lag technique for the first time to off-limb observations of an AR. Since TR emission is not present above the limb, the occurrence of zero time lags is greatly diminished, supporting the conclusion that zero time lags measured on the disk are due to TR plasma. Lastly, we show that the ''coronal'' channels in AIA can be dominated by bright TR emission. When defined in a physically meaningful way, the TR reaches a temperature of roughly 60% the peak temperature in a flux tube. The TR resulting from impulsive heating can extend to 3 MK and higher, well within the range of the ''coronal'' AIA channels.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Saturation mechanisms of backward stimulated Raman scattering in a one-dimensional geometry

In this paper, we investigate the saturation mechanisms of backward stimulated Raman scattering (BSRS) induced by nonlinear kinetic effects. In particular, we stress the importance of accounting for both the nonlinear frequency shift of the electron plasma wave and the growth of sidebands, in order to understand what stops the coherent growth of Raman scattering. Using a Bernstein-Greene-Kruskal approach, we provide an estimate for the maximum amplitude reached by a BSRS-driven plasma wave after the phase of monotonic growth. This estimate is in very good agreement with the results from kinetic simulations of stimulated Raman scattering using both a Vlasov and a Particle in Cell code. Our analysis, which may be generalized to a multidimensional geometry, should provide a means to estimate the limits of backward Raman amplification or the effectiveness of strategies that aim at strongly reducing Raman reflectivity in a fusion plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Numerical simulations of impulsively generated Alfvén waves in solar magnetic arcades

We perform numerical simulations of impulsively generated Alfvén waves in an isolated solar arcade, which is gravitationally stratified and magnetically confined. We study numerically the propagation of Alfvén waves along the magnetic structure that extends from the lower chromosphere, where the waves are generated, to the solar corona, and analyze the influence of the arcade size and the width of the initial pulses on the wave propagation and reflection. Our model of the solar atmosphere is constructed by adopting the temperature distribution based on the semi-empirical VAL-C model and specifying the curved magnetic field lines that constitute the asymmetric magnetic arcade. The propagation and reflection of Alfvén waves in this arcade is described by 2.5-dimensional magnetohydrodynamic equations that are numerically solved by the FLASH code. Our numerical simulations reveal that the Alfvén wave amplitude decreases as a result of a partial reflection of Alfvén waves in the solar transition region, and that the waves that are not reflected leak through the transition region and reach the solar corona. We also find the decrement of the attenuation time of Alfvén waves for wider initial pulses. Moreover, our results show that the propagation of Alfvén waves in the arcade is affected by the spatial dependence of the Alfvén speed, which leads to phase mixing that is stronger for more curved and larger magnetic arcades. We discuss the processes that affect the Alfvén wave propagation in an asymmetric solar arcade and conclude that besides phase mixing in the magnetic field configuration, the plasma properties of the arcade, the size of the initial pulse, and the structure of the solar transition region all play a vital role in the Alfvén wave propagation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Understanding detachment of the W7-X island divertor

The fundamental behavior of the W7-X island divertor under detached conditions, which has been theoretically predicted with the EMC3-Eirene code, is re-examined here under the experimental conditions achieved so far and compared with the first experimental results. Both simulations and experiments cover a range of divertor configurations and plasma parameters, and show the following common trends: (1) with rising impurity radiation, the target heat load decreases 'uniformly' over the entire target surface in the sense that both the peak and average heat loads can drop by an order of magnitude. Impurity radiation (mainly from intrinsic carbon) occurs primarily at the plasma edge and the resulting negative impact on the stored energy is less than 10%. (2) When the total radiation exceeds a critical level, the target particle flux (the recycling flux Γrecy) begins to fall and can drop by a factor of 3–5 at high radiation levels without an obvious indication of significant volume recombination. (3) While Γ recy decreases, the divertor neutral pressure continues to build up and reaches a maximum, at which point Γ recy has declined significantly. (4) During detachment, the electron temperature at the last closed flux surface falls in a way that is not quantitatively understandable from parallel classical heat conduction processes. This paper presents a physical explanation of the numerical/experimental results described above. Furthermore, using the EMC3-Eirene code as a diagnostic tool, we are able, apparently for the first time, to provide a full quantitative analysis of each transport channel in the island divertor, aiming to clarify how the island divertor plasma self-regulates to maintain particle, energy, and momentum balance under detached conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Infrastructure for Multiphysics Software Integration in High Performance Computing-Aided Science and Engineering

The project described in this report constructed and exercised an innovative multiphysics coupling toolkit called the Illinois Rocstar MultiPhysics Application Coupling Toolkit (IMPACT). IMPACT is an open source, flexible, natively parallel infrastructure for coupling multiple uniphysics simulation codes into multiphysics computational systems. IMPACT works with codes written in several high-performance-computing (HPC) programming languages, and is designed from the beginning for HPC multiphysics code development. It is designed to be minimally invasive to the individual physics codes being integrated, and has few requirements on those physics codes for integration. The goal of IMPACT is to provide the support needed to enable coupling existing tools together in unique and innovative ways to produce powerful new multiphysics technologies without extensive modification and rewrite of the physics packages being integrated. There are three major outcomes from this project: 1) construction, testing, application, and open-source release of the IMPACT infrastructure, 2) production of example open-source multiphysics tools using IMPACT, and 3) identification and engagement of interested organizations in the tools and applications resulting from the project. This last outcome represents the incipient development of a user community and application echosystem being built using IMPACT. Multiphysics coupling standardization can only come from organizations working together to define needs and processes that span the space of necessary multiphysics outcomes, which Illinois Rocstar plans to continue driving toward. The IMPACT system, including source code, documentation, and test problems are all now available through the public gitHUB.org system to anyone interested in multiphysics code coupling. Many of the basic documents explaining use and architecture of IMPACT are also attached as appendices to this document. Online HTML documentation is available through the gitHUB site. There are over 100 unit tests provided that run through the Illinois Rocstar Application Development (IRAD) lightweight testing infrastructure that is also supplied along with IMPACT. The package as a whole provides an excellent base for developing high-quality multiphysics applications using modern software development practices. To facilitate understanding how to utilize IMPACT effectively, two multiphysics systems have been developed and are available open-source through gitHUB. The simpler of the two systems, named ElmerFoamFSI in the repository, is a multiphysics, fluid-structure-interaction (FSI) coupling of the solid mechanics package Elmer with a fluid dynamics module from OpenFOAM. This coupling illustrates how to combine software packages that are unrelated by either author or architecture and combine them into a robust, parallel multiphysics system. A more complex multiphysics tool is the Illinois Rocstar Rocstar Multiphysics code that was rebuilt during the project around IMPACT. Rocstar Multiphysics was already an HPC multiphysics tool, but now that it has been rearchitected around IMPACT, it can be readily expanded to capture new and different physics in the future. In fact, during this project, the Elmer and OpenFOAM tools were also coupled into Rocstar Multiphysics and demonstrated. The full Rocstar Multiphysics codebase is also available on gitHUB, and licensed for any organization to use as they wish. Finally, the new IMPACT product is already being used in several multiphysics code coupling projects for the Air Force, NASA and the Missile Defense Agency, and initial work on expansion of the IMPACT-enabled Rocstar Multiphysics has begun in support of a commercial company. These initiatives promise to expand the interest and reach of IMPACT and Rocstar Multiphysics, ultimately leading to the envisioned standardization and consortium of users that was one of the goals of this project.

42 ENGINEERING↗
DOE OSTI2013

Electron acceleration in cavitated laser produced ion channels

This paper is concerned with the channeling of a relativistic laser pulse in an underdense plasma and with the subsequent generation of fast electrons in the cavitated ion channel. The laser pulse has a duration of several hundreds of femtoseconds and its power P{sub L} exceeds the critical power for laser channeling P{sub ch}, with P{sub ch}≈1.1P{sub c}, P{sub c} denoting the critical power for relativistic self-focusing. The laser pulse is focused in a plasma of electron density n{sub 0} such that the ratio n{sub 0}/n{sub c} lies in the interval [10{sup −3},10{sup −1}], n{sub c} denoting the critical density. The laser-plasma interaction under such conditions is investigated by means of three dimensional Particle-In-Cell (PIC) simulations. It is observed that the steep laser front gives rise to the excitation of a surface wave which propagates along the sharp radial boundaries of the electron free channel created by the laser pulse. The mechanism responsible for the generation of relativistic electrons observed in the PIC simulations is also analyzed by means of a test particles code. The fast electrons are found to be generated by the combination of a surface wave and of the betatron resonance. The maximum electron energy observed in the simulations is scaled as a function of P{sub L}/P{sub c}; it reaches 350–600 MeV for P{sub L}/P{sub c} = 70–140.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Towards a novel stellar opacity measurement scheme using stability properties of double ablation front structures

A novel method for bringing sample elements to hydrodynamic conditions relevant to the base of the solar convection zone is investigated. The method is designed in the framework of opacity measurements and exploits the temporal and spatial stability of hydrodynamic parameters in counter-propagating Double Ablation Front (DAF) structures. The physics of symmetric DAF structures is first studied in 1D geometries to assess the influence of tracer layers in the target. These results are used to propose an experimental design compatible with the OMEGA [Boehly et al., Opt. Commun. 133(1-6), 495–506 (1997)] laser. Radiative-hydrodynamic simulations conducted using the Chic code [Breil et al., Comput. Fluids 46, 161–167 (2011).] in 2D-axisymmetric geometries suggest that a Fe sample can be brought to an electron temperature of ∼160 eV and electron number density of ∼1.35 × 1023 cm−3. These parameters are reached during a 500 ps window with temporal variations of the order of 10 eV and 1022 cm−3, respectively. This allows for potential time-integrated spectral measurements. During that time, the sample is almost at local thermal equilibrium and 2D spatial gradients in the sample are less than 5% in a 360 μm diameter cylindrical volume, including the potential effects of Hot Electrons (HE) and typical uncertainties related to target fabrication and laser performances. The effects of HEs are assessed using an inline model in Chic. The HEs are found to deposit most of their energy in the cold and dense ablator between the two fronts, leading to a small efficiency loss on the DAF parameters. The calculations also suggest that negligible amounts of unabsorbed HEs are present in the probed volume, thus not affecting the atomic properties of the sample. Potential extensions of the current design to higher sample temperatures within the OMEGA capabilities are briefly discussed.

Colaïtis, A.↗
DOE OSTI2025

EMPOWERED Distributed Energy Resources Permit Accelerator Pilot

This report summarizes the EMPOWERED Distributed Energy Resources (DER) Permit Accelerator Pilot project. Its primary objective was to streamline the design, installation, permitting, and inspections of targeted Distributed Energy Resource (DER) solutions. Streamlining these processes is critical to strengthening and accelerating the adoption of DER solutions nationwide. Moreover, streamlined processes help expedite the clean energy transition and expand energy resilience by improving understanding and capability among key stakeholders like property owners, code officials, and the general workforce. The project team accomplished this by developing and deploying a set of design and permitting guides for simple DER solutions, which they rolled out via a permitting pilot program. The project team developed nine permitting and inspection guides for simple DER solutions, as follows: 1. Single-Family and Duplex: Electric Vehicle Service Equipment (EVSE), and Storage, Solar + Storage (each in two code cycles) 2. Multi-family and Office: EVSE, Storage, Solar + Storage (each in two code cycles) 3. A permitting process guide to support the implementation of the permitting guides. The guides cover key code requirements via plan review and field inspection checklists for code officials and building owners to follow to streamline the permitting and inspections process for DER solutions. Following the development of the initial set of guides, the project team formed two cohorts, each made up of four jurisdictions, to participate in the Permit Accelerator Pilot. Each cohort received training and technical assistance to support them in incorporating the guides into use. The project successfully engaged jurisdictions in two disparate US regions: Chelsea, Somerville, Natick, and Norwood in metro Boston; and Pima County, Town of Gilbert, Flagstaff, and Sedona in Arizona. Feedback gathered from the project was valuable in refining the guides to better address jurisdictional needs. After the pilot period, the project team developed final versions of the guides based on extensive feedback from the cohorts and from external technical peer review. Feedback included needs and barriers related to implementation, as well as detailed technical improvements from external peer review. The final guides were presented in a series of webinars with a national reach, with 192 attendees from 35 states in the live session of the final webinar, and the final guides are posted online and available for use nationwide. Key findings: • The most receptive jurisdictions were those adjacent to other jurisdictions that had already adopted similar initiatives. The least receptive jurisdictions were those in the midst of adopting new code cycles. • The most common barrier to jurisdictional adoption of new guidelines were internal administrative and process delays, more so than technical barriers. • The main technical barrier to adoption proved to be the concern of fire hazard risks posed by energy storage systems. There is some apprehension in the code enforcement community about the safety of the batteries. Education coordinated with fire and safety services will help alleviate concerns and increase confidence in the acceptance of innovative technologies. • Opportunities to further leverage and advance the guides and related resources include outreach, education, and jurisdictional support; updates to the resources to align with code cycle changes and local requirements; and advancements supporting emerging technologies.

14 SOLAR ENERGY↗
DOE OSTI2022

Nonlinear MHD modeling of n = 1 RMP -induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

We report fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

mode coupling↗
DOE OSTI2014

MAGNETIC FIELDS IN COSMOLOGICAL SIMULATIONS OF DISK GALAXIES

Observationally, magnetic fields reach equipartition with thermal energy and cosmic rays in the interstellar medium of disk galaxies such as the Milky Way. However, thus far cosmological simulations of the formation and evolution of galaxies have usually neglected magnetic fields. We employ the moving-mesh code AREPO to follow for the first time the formation and evolution of a Milky Way-like disk galaxy in its full cosmological context while taking into account magnetic fields. We find that a prescribed tiny magnetic seed field grows exponentially by a small-scale dynamo until it saturates around z = 4 with a magnetic energy of about 10% of the kinetic energy in the center of the galaxy's main progenitor halo. By z = 2, a well-defined gaseous disk forms in which the magnetic field is further amplified by differential rotation, until it saturates at an average field strength of ∼6 μG in the disk plane. In this phase, the magnetic field is transformed from a chaotic small-scale field to an ordered large-scale field coherent on scales comparable to the disk radius. The final magnetic field strength, its radial profile, and the stellar structure of the disk compare well with observational data. A minor merger temporarily increases the magnetic field strength by about a factor of two, before it quickly decays back to its saturation value. Our results are highly insensitive to the initial seed field strength and suggest that the large-scale magnetic field in spiral galaxies can be explained as a result of the cosmic structure formation process.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Stochastic Gain Degradation in III-V Heterojunction Bipolar Transistors due to Single Particle Displacement Damage

As device dimensions decrease single displacement effects are becoming more important. We measured the gain degradation in III-V Heterojunction Bipolar Transistors due to single particles using a heavy ion microbeam. Two devices with different sizes were irradiated with various ion species ranging from oxygen to gold to study the effect of the irradiation ion mass on the gain change. From the single steps in the inverse gain (which is proportional to the number of defects) we calculated Cumulative Distribution Functions to help determine design margins. The displacement process was modeled using the Marlowe Binary Collision Approximation (BCA) code. The entire structure of the device was modeled and the defects in the base-emitter junction were counted to be compared to the experimental results. While we found good agreement for the large device, we had to modify our model to reach reasonable agreement for the small device.

42 ENGINEERING↗
DOE OSTI2017

Improved understanding of the hot cathode current modes and mode transitions [Mechanism of the hot cathode current mode transitions]

Hot cathodes are crucial components in a variety of plasma sources and applications, but they induce mode transitions and oscillations that are not fully understood. It is often assumed that negatively biased hot cathodes have a space-charge limited (SCL) sheath whenever the current is limited. Here, we show on theoretical grounds that a SCL sheath cannot persist. First, charge-exchange ions born within the virtual cathode (VC) region get trapped and build up. After the ion density reaches the electron density at a point in the VC, a new neutral region is formed and begins growing in space. In planar geometry, this 'new plasma' containing cold trapped ions and cold thermoelectrons grows towards the anode and fills the gap, leaving behind an inverse cathode sheath. This explains how transitions from temperature-limited mode to anode glow mode occur in thermionic discharge experiments with magnetic fields. If the hot cathode is a small filament in an unmagnetized plasma, the trapped ion region is predicted to grow radially in both directions, get expelled if it reaches the cathode, and reform periodically. Filament-induced current oscillations consistent with this prediction have been reported in experiments. Here, we set up planar geometry simulations of thermionic discharges and demonstrate several mode transition phenomena for the first time. Lastly, our continuum kinetic code lacks the noise of particle simulations, enabling a closer study of the temporal dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Early experiences on the OLCF Frontier system with AthenaPK and Parthenon–Hydro

The Oak Ridge Leadership Computing Facility (OLCF) has been preparing the nation's first exascale system, Frontier, for production and end users. Frontier is based on HPE Cray's new EX architecture and Slingshot interconnect and features 74 cabinets of optimized 3rd Gen AMD EPYC CPUs for HPC and AI and AMD Instinct 250X accelerators. As a part of this preparation, “real-world” user codes have been selected to help assess the functionality, performance, and usability of the system. This article describes early experiences using the system in collaboration with the Hamburg Observatory for two selected codes, which have since been adopted in the OLCF test harness. Experiences discussed include efforts to resolve performance variability and per-cycle slowdowns. Results are shown for a performance portable astrophysical magnetohydronamics code, AthenaPK, and a mini-application stressing the core functionality of a performance portable block-structured adaptive mesh refinement framework, Parthenon-Hydro. Here, these results show good scaling characteristics to the full system. At the largest scale, the Parthenon-Hydro miniapp reaches a total of $1.7$ $\times$ $10^{13}$ zone-cycles/s on 9216 nodes (73,728 logical GPUs) at ≈92% weak scaling parallel efficiency (starting from a single node using a second-order, finite-volume method).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Investigation of a K-band large coaxial relativistic backward wave oscillator

A K-band large coaxial relativistic backward wave oscillator has been investigated by the 2.5-D particle-in-cell code. This device can generate high-power microwave at a constant frequency with a constant efficiency by increasing the radius of the electron beam and the average radius of the slow-wave structure. The simulation results show that the power conversion efficiency can reach 38.8% at the frequency of 25.48 GHz with the output power of 1.65 GW, while the electron beam has the energy of 196 kV and carries the current of 21.6 kA guided by the magnetic field of 2.5 T.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hydrologic Model Data for the East Fork Poplar Creek Watershed Simulated with the Advanced Terrestrial Simulator (ATS): Streamflow and Network Expansion–Contraction Dynamics

This dataset supports hydrologic modeling and stream network expansion–contraction analysis for the East Fork Poplar Creek (EFPC) Watershed in Tennessee. It includes a Jupyter notebook for model setup, model configuration files, simulation outputs, and derived products used to evaluate model performance and investigate stream dynamics under varying hydrologic conditions. The dataset was generated using the Watershed Workflow Python package and the Advanced Terrestrial Simulator (ATS), enabling integrated surface–subsurface hydrologic simulations using a stream-aligned mesh. Outputs include high-resolution time series of streamflow, active network length, water table depth, and related hydrologic variables. Also included are spatially explicit stream persistency indices and classifications of reaches as perennial or non-perennial. These data facilitate reproducibility and support further research on stream intermittency and variability in network extent.The model data archive is organized in following directories:1) model_setup_inputsContains the Watershed Workflow Jupyter notebooks (accessed through any open source code editor), selected input datasets, and resulting ATS input files, including XML files (access through any open source code editor), computational mesh (.exo files can be viewed using Paraview), and meteorological forcing files (.h5 files can be accessed through h5py python package and HDFView open source software). 2) model_outputsIncludes ATS simulation outputs relevant to this study. Time series of spatially integrated or averaged variables (e.g., streamflow, water table depth) are provided as CSV files. Select spatial fields (e.g., ponded depth and water table depth) are saved as pickled Python objects to reduce file size, and can be accessed through pickle package in Python. Key geometry objects from Watershed Workflow—such as the surface mesh and river tree—are also included to support analysis of streamflow persistency and expansion–contraction dynamics. These files can also be accessed through Watershed Workflow Python package.3) model_evaluationProvides observed streamflow time series and field survey-based flow regime classifications used to evaluate model performance. Jupyter notebooks for processing ATS outputs and comparing model predictions with observations to build confidence in the model prior to scientific analysis are also included.4) Q_L_relationshipsContains workflows for generating time series of discharge, active network length, and related hydrologic variables used in the stream network expansion–contraction analysis. Includes routines for delineating baseflow-dominated periods. For each catchment, notebooks and processed data (as pickled DataFrames accessed through Pandas Python package) are provided. 5) figure_scriptsProvides the Jupyter notebooks used to generate the figures presented in the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Description and evaluation of the Community Ice Sheet Model (CISM) v2.1

Abstract. We describe and evaluate version 2.1 of the Community Ice Sheet Model (CISM). CISM is a parallel, 3-D thermomechanical model, written mainly in Fortran, that solves equations for the momentum balance and the thickness and temperature evolution of ice sheets. CISM's velocity solver incorporates a hierarchy of Stokes flow approximations, including shallow-shelf, depth-integrated higher order, and 3-D higher order. CISM also includes a suite of test cases, links to third-party solver libraries, and parameterizations of physical processes such as basal sliding, iceberg calving, and sub-ice-shelf melting. The model has been verified for standard test problems, including the Ice Sheet Model Intercomparison Project for Higher-Order Models (ISMIP-HOM) experiments, and has participated in the initMIP-Greenland initialization experiment. In multimillennial simulations with modern climate forcing on a 4 km grid, CISM reaches a steady state that is broadly consistent with observed flow patterns of the Greenland ice sheet. CISM has been integrated into version 2.0 of the Community Earth System Model, where it is being used for Greenland simulations under past, present, and future climates. The code is open-source with extensive documentation and remains under active development.

58 GEOSCIENCES↗
DOE OSTI2014

Improving the hot-spot pressure and demonstrating ignition hydrodynamic equivalence in cryogenic deuterium–tritium implosions on OMEGA

Reaching ignition in direct-drive (DD) inertial confinement fusion implosions requires achieving central pressures in excess of 100 Gbar. The OMEGA laser system [T. R. Boehly et al., Opt. Commun. 133, 495 (1997)] is used to study the physics of implosions that are hydrodynamically equivalent to the ignition designs on the National Ignition Facility (NIF) [J. A. Paisner et al., Laser Focus World 30, 75 (1994)]. It is shown that the highest hot-spot pressures (up to 40 Gbar) are achieved in target designs with a fuel adiabat of α ≃ 4, an implosion velocity of 3.8 × 107 cm/s, and a laser intensity of ∼1015 W/cm2. These moderate-adiabat implosions are well understood using two-dimensional hydrocode simulations. The performance of lower-adiabat implosions is significantly degraded relative to code predictions, a common feature between DD implosions on OMEGA and indirect-drive cryogenic implosions on the NIF. Simplified theoretical models are developed to gain physical understanding of the implosion dynamics that dictate the target performance. These models indicate that degradations in the shell density and integrity (caused by hydrodynamic instabilities during the target acceleration) coupled with hydrodynamics at stagnation are the main failure mechanisms in low-adiabat designs. To demonstrate ignition hydrodynamic equivalence in cryogenic implosions on OMEGA, the target-design robustness to hydrodynamic instability growth must be improved by reducing laser-coupling losses caused by cross beam energy transfer.

Goncharov, V. N.↗
DOE OSTI2017

Warthog: Coupling Status Update

The Warthog code was developed to couple codes that are developed in both the Multi-Physics Object-Oriented Simulation Environment (MOOSE) from Idaho National Laboratory (INL) and SHARP from Argonne National Laboratory (ANL). The initial phase of this work, focused on coupling the neutronics code PROTEUS with the fuel performance code BISON. The main technical challenge involves mapping the power density solution determined by PROTEUS to the fuel in BISON. This presents a challenge since PROTEUS uses the MOAB mesh format, but BISON, like all other MOOSE codes, uses the libMesh format. When coupling the different codes, one must consider that Warthog is a light-weight MOOSE-based program that uses the Data Transfer Kit (DTK) to transfer data between the various mesh types. Users set up inputs for the codes they want to run, and then Warthog transfers the data between them. Currently Warthog supports XSProc from SCALE or the Sub-Group Application Programming Interface (SGAPI) in PROTEUS for generating cross sections. It supports arbitrary geometries using PROTEUS and BISON. DTK will transfer power densities and temperatures between the codes where the domains overlap. In the past fiscal year (FY), much work has gone into demonstrating two-way coupling for simple pin cells of various materials. XSProc was used to calculate the cross sections, which were then passed to PROTEUS in an external file. PROTEUS calculates the fission/power density, and Warthog uses DTK to pass this information to BISON, where it is used as the heat source. BISON then calculates the temperature profile of the pin cell and sends it back to XSProc to obtain the temperature corrected cross sections. This process is repeated until the convergence criteria (tolerance on BISON solve, or number of time steps) is reached. Models have been constructed and run for both uranium oxide and uranium silicide fuels. These models demonstrate a clear difference in power shape that is not accounted for in a stand-alone BISON run. Future work involves improving the user interface (UI), likely through integration with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Workbench. Furthermore, automating the input creation would ease the user experience. The next priority is to continue coupling the work with other codes in the SHARP package. Efforts on other projects include work to couple the Nek5000 thermo-hydraulics code to MOOSE, but this is in the preliminary stages.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗