Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Compact Nuclear Simulator”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Accelerated radiation transport modeling techniques for pencil beam computed tomography using gamma rays

Monte Carlo radiation transport modeling studies were performed for a compact, and high-resolution gamma-ray computed tomography system designed for imaging irradiated nuclear fuel. The system comprises a 60 Co source – chosen for its highly penetrating 1173 keV and 1332 keV gamma rays – a pair of high-aspect-ratio pencil beam collimators, and an inorganic scintillator detector. Here two acceleration methods are proposed to rapidly model a transmission type gamma-ray tomography system. The first, a variance reduction technique, is based on performing Monte Carlo simulations with a monodirectionally-biased source, sampled from a characteristic sub-volume of the full source volume. The second acceleration method is based on the deterministic calculations using the Beer–Lambert law and detector response characteristics. Comparison of simulations using acceleration approaches with analog simulations of the fully isotropic, full-volume equivalent, show that the Monte Carlo variance reduction technique gives quantitatively accurate predictions for large collimator aspect ratios while the deterministic calculations are semi-quantitative but converge close to the correct result as the collimator aspect ratio increases. As such, these techniques can be used to reduce the computational cost in generating simulated radiographs and tomographs by several orders of magnitude. Experimental validation efforts are currently underway and will be demonstrated in future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Neural Thermal Scattering (NeTS) Modules for Graphite & Beryllium [Slides]

This presentation discusses the motivations behind the project, which include facilitating compact formulation for TSL data, extending AI knowledge, and providing advanced reactor simulation framework. It also discusses how the project accounts for beryllium and graphite complexity, including atomistic, dynamical, and neural. Additionally, new material-informed neural thermal scattering (NeTS) Modules are examined. The implications, conclusions, and plan future work for the project are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Physics and Components syntax to enable a systems-based approach to multiphysics

Simulations in MOOSE have traditionally used kernel and boundary condition classes to describe the equations. Downstream applications leveraged a system called Actions to define a pre-packaged discretization of the equations they solve. Unfortunately, the Action base class was very limited, and most applications implemented the same concepts in their Actions. This led to an increased maintenance burden and a reduction in coupling opportunities, save for the use of MultiApps which renders each input mostly independent. With the introduction of multi-system capabilities in MOOSE, there is growing interest in defining entire simulations of complex multiphysics systems in a single input file. By introducing a new Physics system, with its dedicated syntax and a new base class providing wide-ranging capabilities, we are now able to define multiple equations in a single input file in a compact and user-friendly way. With new interactions between Physics and the Component system, these equations can be defined on each component of a complex system. In this talk, we will present the capabilities of these new systems, their interactions, and how to define complex systems multiphysics simulations with Physics and Components.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Employing MACS/ViBRANT as a Surrogate MARVEL Reactor for Startup Reactivity Tuning and Supervisory Control Processes

Advanced nuclear reactors are a key part of the future of nuclear energy both in the United States and globally. They offer unique benefits for various energy-demanding applications, including use in remote locations, compact size, modular manufacturing, remote monitoring, low and/or variable power rating operation, and reliance on novel technologies to enhance operational safety. To achieve economic feasibility, advanced reactors must significantly reduce their workforces in comparison with the current fleet. Achieving this reduction will occur through reducing staff workloads using technology to achieve autonomous or semi-autonomous operations, demonstrated by comprehensive testing and validation activities. These operations will require both software and hardware platforms during the design and testing phases. While simulations are useful during the design phase, their performance can significantly deviate during actual deployment on hardware. This report presents the outcomes of a collaborative technical initiative between the U.S. Department of Energy (DOE) Microreactor Program (MRP) and Advanced Sensors and Instrumentation (ASI) Program. The collaboration utilized the Microreactor Automated Control System (MACS) hardware platform to bridge the gap between theoretical reactor design and actual startup and control operations. Two key use cases were investigated: facilitating the startup testing period and demonstrating supervisory control. The first use case details the key Microreactor Applications Research Validation and Evaluation (MARVEL) reactor startup physics testing activities conducted using the MACS platform. These activities included drum worth measurements, shutdown margin assessment, temperature feedback analysis, and scram time evaluation, as well as unique testing that would apply to the MARVEL reactor to demonstrate the testing methodologies in a low-risk environment. The MACS platform, serving as a surrogate representation of the MARVEL reactor, proved instrumental in performing these tests. The exercise revealed aspects that led to optimized processes, refined hardware design, and enhanced base software capabilities. By maturing methods and technologies in this manner, the initiative promises to reduce wasted time in the actual on-site reactor deployment effort, thereby saving significant time and resources. The second use case focuses on the development and implementation of supervisory control methods aimed at managing core tilt, which can result from asymmetrical operations or manufacturing imperfections in fuel rods or reactivity control devices. A key objective was to assess and compare the use of artificial intelligence (AI) for supervisory control. The effort aimed to define the role of supervisory control to enhance performance without risking control instability. This effort explored three distinct approaches: rules-based (RB) methods, optimization techniques, and reinforcement learning (RL) algorithms. Each approach was evaluated for its ease of implementation, its usability, and its effectiveness in responding to asymmetries in neutron flux. Comparative analysis of these approaches provided valuable insights into their applicability and effectiveness, offering a robust framework for advanced reactor operations. Together, these two use cases highlight the potential of hardware test beds to help streamline the design, operation, and control of advanced nuclear reactors. This collaborative effort underscores the importance of continued innovation and experimentation in achieving the next generation of safe, reliable, and economically viable nuclear energy solutions.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

FGMS-poster

Idaho National Laboratory (INL) performs post irradiation examination (PIE) of tri-structural isotropic (TRISO)-coated particle fuel to help qualify it for high temperature gas cooled reactors. TRISO fuel compacts are re-irradiated in the Neutron Radiography Reactor (NRAD) to generate the short lived fission products needed for fission product release testing. The Fuel Accident Condition Simulator (FACS) furnace and newly added Screen Neutron Irradiated Fuel for Failure (SNIFF) furnace heat the compacts in helium to temperatures of up to 2,000°C, prompting fission product release—predominantly gaseous xenon and krypton isotopes and condensable products such as cesium—from failed particles. These released isotopes are transported to a fission gas monitoring system (FGMS 1 or FGMS 3), where they accumulate in cryogenic cold traps and are quantified using high-purity germanium (HPGe) detectors. The addition of SNIFF and FGMS 3 increases throughput by enabling simultaneous testing of multiple compacts. Furthermore, automated INL developed software provides continuous, near-real time monitoring of fission product inventories and manages the liquid nitrogen cooling of the traps. These system enhancements improve the efficiency, data quality, and testing capacity of TRISO fuel performance evaluations.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2022

Application of recoil-imaging time projection chambers to directional neutron background measurements in the SuperKEKB accelerator tunnel

Gaseous time projection chambers (TPCs) with high readout segmentation are capable of reconstructing detailed 3D ionization distributions of nuclear recoils resulting from elastic neutron scattering. Here, using a system of six compact TPCs with pixel ASIC readout, filled with a 70:30 mixture of He:CO 2 gas, we analyze the first directional measurements of beam-induced neutron backgrounds in the tunnel regions surrounding the Belle II detector at the SuperKEKB e + e – collider. With the use of 3D recoil tracking, we show that these TPCs are capable of maintaining nearly 100% nuclear recoil purity to reconstructed ionization energies (E reco ) as low as 5 keV ee . Using a large sample of Monte-Carlo (MC)-simulated 4 He, 12 C, and 16 O recoil tracks, we find consistency between predicted and measured recoil energy spectra in five of the six TPCs, providing useful validation of the neutron production mechanisms modeled in simulation. Restricting this sample to 4 He recoil tracks with E reco > 40 keV ee , we further demonstrate axial angular resolutions within 8° and we introduce a procedure that under suitable conditions, correctly assigns the vector direction to 91% of these simulated 4 He recoils. Applying this procedure to assign vector directions to measured 4 He recoil tracks, we observe consistency between the angular distributions of observed and simulated recoils, providing first experimental evidence of localized neutron “hotspots” in the accelerator tunnel. Observed rates of nuclear recoils in these TPCs suggest that simulation overestimates the neutron flux from these hotspots. Despite this, we estimate these hotspots to produce the majority of neutron backgrounds in the accelerator tunnel at SuperKEKB’s target luminosity of 6.3 x 10 35 cm –2 s –1 , making them important regions to continue to monitor.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Correlated Signal Analysis for Nuclear Emergency Response

In the context of nuclear emergency response (NER) scenarios, the critical task is to quickly identify a "black box" as a potential threat, as failing to do so could have catastrophic consequences. Techniques for passive assay of a “black box” typically include gamma-ray spectroscopy and neutron coincidence/multiplicity counting. However, there are significant challenges associated with these type of measurements. First, the presence of intervening materials can obstruct the detection of relevant signatures. Second, the presence of strong non-fission neutron sources, like (α, n) emitters, can add uncertainties to the neutron multiplicity analysis. In our LDRD-MFR Phase II work, a portable neutron spectrometer, called the Compact Fast Neutron Spectrometer (CFNS), was developed for NER applications. The CFNS system leverages information-rich neutron energy spectra to derive actionable information. This work investigates the use of correlated signals in the CFNS from special nuclear material (SNM) to characterize physical properties, such as intervening shielding material and fission to non-fission neutron contributions. In this report, the Phase II results from bulk SNM measurements at the National Criticality Experiments Research Center (NCERC) are briefly discussed along with the motivation for this work. Afterwards, simulations using the MCNPX-PoliMi transport code are discussed, which were used to expand our correlated signal study. Signal triggered analysis for neutron multiplicity extraction will also be discussed. Lastly, the use neutron-photon correlations for intervening material identification are shown, along with the use of correlated neutron energy spectra for α-ratio extraction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Updating and Refining BISON TRISO Assessment

Several advanced reactor concepts consider tristructural isotropic (TRISO) fuel in their designs. To support the advanced reactor industry, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program within the U.S. Department of Energy (DOE) is developing modeling and simulation capabilities within the BISON fuel performance code. Assessment of these capabilities against available experi- mental data, such as the Advanced Gas Reactor (AGR) program, is essential in instilling confidence in the code predictions. Associated uncertainties in model inadequacy, experimental measurements, and manufacturing are also of importance. This report discusses the updating and standardization of the existing assessment suite in BISON for TRISO fuel performance analysis. The cases considered the AGR-1, AGR-2, and AGR-3/4 experimental programs. The report also includes a preliminary investigation into the sources of model inadequacy associated with silver (Ag) release fractions using the AGR-2 Compact 2-1-2 as a case study. The results indicate that uncertainties in the operating temperatures supplied to the model and the form of those temperature conditions (spatially averaged versus spatially resolved) both affect Ag release predictions. This suggests that previously observed model inadequacy can be attributed to both physical models and the operating conditions supplied to them. These findings will be leveraged to guide development of assessment cases for additional experiments and refinement of physical models for behaviors that potentially influence Ag release, such as palladium penetration.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Occurrence of Gravitational Collapse in the Accreting Neutron Stars of Binary-driven Hypernovae

The binary-driven hypernova (BdHN) model proposes long gamma-ray bursts (GRBs) originate in binaries composed of a carbon–oxygen (CO) star and a neutron star (NS) companion. The CO core collapse generates a newborn NS and a supernova that triggers the GRB by accreting onto the NSs, rapidly transferring mass and angular momentum to them. This article aims to determine the conditions under which a black hole (BH) forms from NS collapse induced by the accretion and the impact on the GRB's observational properties and taxonomy. We perform three-dimensional, smoothed particle hydrodynamics simulations of BdHNe using up-to-date NS nuclear equations of state, with and without hyperons, and calculate the structure evolution in full general relativity. We assess the binary parameters leading either NS in the binary to the critical mass for gravitational collapse into a BH and its occurrence time, t col . We include a nonzero angular momentum of the NSs and find that t col ranges from a few tens of seconds to hours for decreasing NS initial angular momentum values. BdHNe I are the most compact (about 5 minute orbital period), promptly form a BH, and release ≳10 52 erg of energy. They form NS–BH binaries with tens of kiloyears merger timescales by gravitational-wave emission. BdHNe II and III do not form BHs, and release ~10 50 –10 52 erg and ≲10 50 erg of energy, respectively. They form NS–NS binaries with a range of merger timescales larger than for NS–BH binaries. In some compact BdHNe II, either NS can become supramassive, i.e., above the critical mass of a nonrotating NS. Magnetic braking by a 10 13 G field can delay BH formation, leading to BH–BH or NS–BH with tens of kiloyears merger timescales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Assessment and Improvement of Fission Product Transport Predictions of Particle Fuel in BISON

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This study was motivated by the need to incorporate more physics-based models in BISON in order to foster tri-structural isotropic (TRISO) applications. This document details the integration of new modeling capabilities in BISON, including (1) development of pyrolytic carbon (PyC) and silicon carbide (SiC) layer anisotropic thermal and mass transport capabilities, (2) verification of the mass diffusion solution in TRISO modeling, (3) calibration of fission product diffusivity using Advanced Gas Reactor (AGR) experiments, (4) improved fission product release modeling by developing compact diffusion modeling capabilities, and (5) documentation of accelerated failure analysis on the BISON website. Improvements made to the diffusion models and parameters were documented and validated against AGR-1 and -2 experiment data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Validation of the Porous Medium Approximation for Hydrodynamics Analysis in Compact Heat Exchangers

Compact heat exchangers (HXs) have gained attention in recent years in various fields such as solar and nuclear power generation, oil and gas, and refrigeration due to their low cost, high power density, and robustness in high-pressure and/or high-temperature environments. However, the large difference between a compact HX's overall dimensions (~m) and the much smaller scale of its channels (~mm) makes it challenging to model the entire HX at once, due to computational limitations. In this work, we treat the channeled region of a compact HX as a porous medium (PM) to circumvent the need to model/mesh each individual channel. This allows us to simulate the entire HX, including both the header and channeled regions while maintaining the computational cost at a practical level. Although the porous medium approach has been used to model heat exchangers, its validity is still questionable because (1) the resistance coefficients are heavily data-based and thus difficult to be applied to new heat exchangers and (2) the validation has been focused on matching the overall pressure drop in the channel region, which does not address whether such model can predict detailed pressure and velocity field. For the first time, this work addresses under what circumstances and with what uncertainty does the PM approach work for hydrodynamics modeling in compact HXs. By answering these questions, we introduce the PM approach as a powerful tool for HX hydrodynamics modeling that can predict not only the overall pressure drop but also the detailed pressure and velocity distributions.

Engineering↗
DOE OSTI2022

Cell Dancoff-based embedded self-shielding capability for doubly heterogeneous particulate fuels in SCALE/Polaris

Polaris is a two-dimensional (2D) transport lattice physics capability in the SCALE code system for light-water reactor (LWR) analysis. SCALE/Polaris is being extended to treat doubly heterogeneous particulate fuels in LWRs with accident-tolerant fuels and prismatic high-temperature gas-cooled reactor (HTGR) fuel. Recently, a Dancoff-based Wigner-Seitz approximation capability based on the embedded self-shielding method (DWA-ESSM) was implemented into Polaris for efficient computation. A new double heterogeneity treatment capability based on DWA-ESSM in conjunction with the Hebert's collision probability method for double heterogeneity was developed and implemented for Polaris. The new capability was verified by performing benchmark calculations for the prismatic HTGR fuel compact problems with various design parameters through a code-to-code comparison between the Polaris and continuous-energy Monte Carlo results. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

High-Fidelity CFD Assessments of Flow Resistance in a 61-Pin Wire-Wrapped Assembly with Partially Blocked Channels

The examination of thermal-hydraulic behaviors in wire-wrapped rod bundles continues to be an active area of research. The sodium fast reactor, a prominent candidate in next-generation nuclear designs, utilizes a hexagonal configuration of wire-wrapped fuel pins. Here, the potential for channel blockage within this compact arrangement poses a significant safety challenge, spurring a number of recent experimental and computational investigations to evaluate its impact on flow and heat transfer. The present work aims to benchmark the high-fidelity NekRS computational fluid dynamics (CFD) solver in predicting the pressure drops associated with substantial blockages, using available experimental data as a reference. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at a range of low to moderate Reynolds numbers (487 ≤ Re ≤ 14 600). The NekRS solver demonstrates an exponential reduction of spatial discretization error with increasing polynomial order. The high level of agreement between the numerical results and measured data confirms the accuracy and consistency of the present numerical approach. This benchmark study establishes the capability of NekRS to perform reliable hydrodynamic simulations for sodium fast reactor applications and supports its use in design, licensing, and safety analyses.

CFD Benchmarking↗
NASA NTRS2005

Heat Pipe Reactor Dynamic Response Tests: SAFE-100 Reactor Core Prototype

The SAFE-I00a test article at the NASA Marshall Space Flight Center was used to simulate a variety of potential reactor transients; the SAFEl00a is a resistively heated, stainless-steel heat-pipe (HP)-reactor core segment, coupled to a gas-flow heat exchanger (HX). For these transients the core power was controlled by a point kinetics model with reactivity feedback based on core average temperature; the neutron generation time and the temperature feedback coefficient are provided as model inputs. This type of non-nuclear test is expected to provide reasonable approximation of reactor transient behavior because reactivity feedback is very simple in a compact fast reactor (simple, negative, and relatively monotonic temperature feedback, caused mostly by thermal expansion) and calculations show there are no significant reactivity effects associated with fluid in the HP (the worth of the entire inventory of Na in the core is .<el, so fluid movement and temperature changes will cause very minor effects). In previous SAFE-100 tests, the point kinetics model was based on core thermal expansion via deflection measurements. It was found that core deflection was a strung function of how the SAFE-100 modules were fabricated and assembled (in terms of straightness, gaps, and other tolerances). To remove the added variable of how this particular core expands as compared to a different concept, it was decided to use a temperature based feedback model (based on several thermocouples placed throughout the core).

Bragg-Sitton, Shannon M.↗
DOE OSTI2022

Design of a neutral thermal scattering (NeTS) module for hydrogen in light water

The accurate representation of thermal scattering law (TSL) data is integral to the design and characterization of many modern nuclear systems, particularly those using light water as a moderator/coolant. As a material-dependent distribution over energy-momentum phase space, the TSL may exhibit a variety of unique and relevant conditional (e.g., temperature, pressure) and compositional (e.g., porosity, stoichiometry, radiation damage) dependencies. Currently, there are various approaches to incorporate temperature dependence, which is especially important in coupled neutronic-thermal hydraulic simulations. In each approach, there is an inherent tradeoff between memory consumption and accuracy. Some techniques require tens to hundreds of MBs or more, while others fail to reproduce the underlying data to within 10% error over the considered input domain, despite having a reduced storage burden. This work aims to address both sides of the tradeoff simultaneously by implementing a novel deep learning (DL) approach to TSL representation. The neural thermal scattering (NeTS) concept, which is amenable to an arbitrary number of dependencies, is demonstrated via the inclusion of temperature dependence into a highly compact, highly accurate functional form of the multi- variate TSL (i.e., S(α, β, T)) for hydrogen in light water. Resulting storage requirements are on the order of 100 kB, and median and maximum percent deviations are on the order of 0.1% and 1%, respectively. These measures represent a step improvement over previous techniques. Notably, the developed neural network and feature methodology build on those employed in prior work on beryllium oxide. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Square Pulse Thermoreflectance Technique for the Measurement of Thermal Properties

We report on a laser-based square pulse thermoreflectance (SPTR) technique for the measurement of thermal properties for a wide range of materials. SPTR adopts the pump-probe thermoreflectance principle to monitor the evolution of local temperature after square pulse excitation. The technique features a compact setup, high spatial resolution, and fast data collection. By comparing the acquired SPTR signals with a continuum heat transfer model, material thermal properties can be obtained. Taking advantage of various spot sizes and modulation frequencies, SPTR can measure both the thermal diffusivity and thermal conductivity of poorly to moderately conductive materials and the thermal conductivity of conductive materials with satisfactory accuracy, with potential to be applied to more conductive materials. The technique was validated on three materials: fused silica, single crystal CaF2 and single crystal nickel (with conductivities ranging from 1 W·m -1 ·K -1 to 100 W·m -1 ·K -1 ) with typical measurement errors of 5 % to 20 %. The leading sources of error have been identified by Monte Carlo simulations, and the primary limitations of SPTR are discussed. The compact, fiberized platform we describe here will allow instruments based on this methodology to be deployed in complex, multi-analytical environments for the type of high-throughput correlative analyses that are key to materials design and discovery.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Modeling and analysis of the tritium fuel cycle for ARC- and STEP-class D-T fusion power plants

Abstract The limited tritium resources available for the first fusion power plants (FPPs) make fuel self-sufficiency and tritium inventory minimization leading issues in FPP design. This work builds on the model proposed by Abdou et al (2020 Nucl. Fusion 61 013001), which analyzed the fuel cycle (FC) of a DEMOnstration nuclear FPP-class FPP with a time-dependent system-level model. Here, we use a modified version of their model to analyze the FC of an Affordable, Robust, Compact (ARC)-class tokamak and two versions of a Spherical Tokamak for Energy Production (STEP)-class tokamak. The ARC-class tokamak breeds tritium in a 2LiF + BeF 2 liquid immersion blanket, while the STEP-class tokamak breeds tritium utilizing either a liquid-lithium blanket design or an encapsulated breeding blanket. A time-dependent system-level model is developed in Matlab Simulink ® to simulate the evolution of tritium flows and tritium inventories in the FC. The main goals of this work are to assess tritium self-sufficiency of the ARC- and STEP-class designs and to determine quantitative design requirements that can be used to analyze the adequacy of a proposed FC system. These design requirements are aimed at achieving a low tritium inventory doubling time ( t d ) and a low start-up inventory ( I s t a r t u p ) while keeping the required tritium breeding ratio (TBR r ) as low as possible. We also consider how improvements in FC technology and POs affect TBR r and I s t a r t u p . The model results show that TBR r for ARC- and STEP-class FPPs should be achievable if the tritium burn efficiency (TBE) reaches 0.5%–1% (TBR r < 1.2). This assumes significant, but attainable, improvements over current abilities. However, the model results indicate that an FPP must achieve ambitious performance targets, including FPP availability > 70%, tritium processing time < 4 h, and the implementation of direct internal recycling (DIR). If future research yields major improvements to achievable TBE, it may be possible to achieve tritium self-sufficiency while operating at lower availability and without implementing DIR.

Physics↗
DOE OSTI2022

Long-term 3D MHD simulations of black hole accretion discs formed in neutron star mergers

ABSTRACT We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in flash 4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disc simulations that vary the initial magnetic field geometry and disc compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modelling the BH’s gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively inefficient state in the disc. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disc compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively inefficient state at late times. The electron fraction distribution of the disc outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding ${\sim} 50{{\ \rm per\ cent}}$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

Fahlman, Steven (ORCID:0000000348759940)↗