Search NASA⌕ Search

SEARCH · Search NASA

Results for “ANSTO”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Current Status of the Synroc Waste Treatment Facility - 20263

Construction of the Synroc Waste Treatment Facility is well underway at ANSTO. This plant has been designed to treat intermediate level liquid waste arising from the production of Mo-99 within the newly constructed ANSTO Nuclear Medicine (ANM) precinct [1]. The facility utilises ANSTO Synroc technology that has been tailored to the chemical, physical, and radiological properties of the waste. The result is a highly durable wasteform with a significant reduction in the final volume of the treated waste. This paper presents the construction status of the Synroc Waste Treatment Facility and the status of the associated technology maturation plan. Construction of the building is scheduled for completion in late 2020 followed by process installation and commissioning. A key component of the technology maturation plan has been the construction of an Inactive Engineering Demonstrator [2]. This has aided in the mitigation of risks with respect to technology selection, process integration, and process nuclearisation. Furthermore, the demonstration facility has provided an environment for the design and development of the instrumentation and control philosophy resulting in a seamlessly integrated plant control system. Successful demonstration of the technology by the inactive engineering demonstrator has significantly reduced the associated risks. Details of the Synroc Waste Treatment Facility project progress will also be presented. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Technology Maturation of Hot Isostatic Pressing for Nuclear Waste Treatment - 20259

ANSTO is constructing a nuclear waste treatment facility for the immobilization of intermediate-level liquid waste (ILLW) from its nuclear medicine production [1]. This facility shall deploy ANSTO Synroc process technology that converts a liquid waste stream into a durable solid wasteform. Hot Isostatic pressing (HIP) shall be employed to consolidate the dry powdered waste effectively reducing the waste volume and immobilizing the waste components. This paper presents the technology maturation for HIP technology for utilisation in a nuclear waste treatment facility. The HIP has been designed to operate in a remote shielded enclosure, fully integrated with the front-end powder production. Attention has been given in the HIP system design to ensure that it can operate and be maintained to achieve the required production output of the facility. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Thermal Process Technology for Nuclear Applications - 20379

ANSTO's Synroc technology has been developed to provide a safe, secure matrix for the immobilization and final disposal of radioactive waste. Synroc technology will be used to manage radioactive wastes from the production of the radioisotope Molybdenum-99 (Mo-99). This paper shall outline various stages of the process development with specific reference to the thermal treatment technology of calcination. Calcination is a key step in the Synroc process [1-2].The rotary thermal processing system includes: an advanced heating element design for increased robustness and ease of remote operation and maintenance, an enhanced modular design of components for ease of remote maintenance in a hot cell and in compliance with hot cell radioactive environment requirements for safety, reliability and maintainability. In addition to thermal treatment of waste from nuclear medicine production, this technology provides solutions for a variety of nuclear materials processing applications including sintering UO{sub 2} pellets for reactor fuel rods, oxidation of UO{sub 2} pellets, swarf, and powder to U{sub 3}O{sub 8}, de-nitration of Uranyl nitrate and hydrofluorination of UO{sub 2} pellets. The paper will also discuss thermal processing solutions for a range of nuclear applications. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Subventricular Accumulation of Cu in the Aging Mouse Brain Does Not Associate with Anticipated Increases in Markers of Oxidative Stress

Natural aging is associated with mild memory loss and cognitive decline, and age is the greatest risk factor for neurodegenerative diseases, such as Alzheimer’s disease. There is substantial evidence that oxidative stress is a major contributor to both natural aging and neurodegenerative disease, and coincidently, levels of redox active metals such as Fe and Cu are known to be elevated later in life. Recently, a pronounced age-related increase in Cu content has been reported to occur in mice and rats around a vital regulatory brain region, the subventricular zone of lateral ventricles. In our study herein, we have characterized lateral ventricle Cu content in a unique murine model of accelerated aging, senescence accelerated mouse-prone 8 (SAMP8) mice. Our results confirm an age-related increase in ventricle Cu content, consistent with the studies by others in wild-type mice and rats. Specifically, we observed Cu content to increase over the time frame 1 to 5 months and 5 to 9 months, but interestingly, no significant increase occurred between 9 and 12 months (although brain Cu content at 12 months was significantly elevated relative to 1 and 5 month-old animals). Despite the magnitude of Cu increase observed within the cells that comprise the subventricular zone of lateral ventricles (average 3 mM Cu, with isolated subcellular concentrations of 17 mM), we did not detect spectroscopic markers of thiol oxidation, protein aggregation, or lipid oxidation. The lack of evidence for oxidative stress in ex vivo animal tissue is in contrast to in vitro studies demonstrating that thiol, protein, and lipid oxidation is pronounced at these Cu concentrations. We suggest that our findings most likely indicate that the Cu ions in this brain region are sequestered in an unreactive form, possibly extended chains of Cu-thiolate complexes, which do not readily redox cycle in the aqueous cytosol. In conclusion, these results also appear to partially challenge the long-held view that age-related increases in brain metal content drive oxidative stress as we did not observe a concomitant association between age-related Cu increase and markers of oxidative stress, nor did we observe a net increase in Cu content between mice aged 9 and 12 months.

X-ray fluorescence↗
DOE OSTI2025

Multi-objective surrogate-assisted calibration of CPFEM models using macroscopic response and in situ EBSD measurements of grain reorientation trajectories

Crystal plasticity finite element method (CPFEM) models are widely used to simulate the deformation behaviour of polycrystalline materials, but their calibration is often limited by their high computational cost and the non-convexity of the optimisation landscape. Here, this study develops a multi-objective surrogate-assisted calibration workflow that couples a multi-objective genetic algorithm (MOGA) with an adaptively trained deep neural network (DNN) surrogate model to efficiently identify CPFEM parameters from experimental data. The workflow is demonstrated on three crystal plasticity (CP) formulations of increasing complexity — Voce hardening (VH), two-coefficient latent hardening (LH2), and six-coefficient latent hardening (LH6) — using in situ electron backscatter diffraction (EBSD) measurements of Alloy 617 under uniaxial tensile loading. The CPFEM models are calibrated against the experimentally observed stress–strain response and reorientation trajectories of eight grains, then validated against eight additional trajectories and overall texture evolution. Across the CP formulations, the macroscopic response was reproduced reliably, while differences emerged in the robustness and accuracy of the grain-scale predictions. Including grain reorientation trajectories in the multi-objective calibration improved texture evolution predictions and filtered out physically inconsistent parameter sets that can arise from calibrating against only the stress–strain data. The workflow also demonstrates good transferability of calibrated parameters from a low- to a high-fidelity microstructural model. These results provide practical guidance for integrating in situ microstructural data into CPFEM through efficient, repeatable, and physically meaningful multi-objective calibration.

Crystal plasticity finite element method↗
DOE OSTI2025

Benchmarking Monte Carlo codes for the modelling of low-energy neutron production target reactions

The increasing adoption of accelerator-based neutron sources (ABNS) for applications including neutron capture therapy (NCT) research has highlighted the need for accurate simulation tools. Precise modelling of the neutron production target is crucial to ensure that simulated predictions of neutron beam characteristics used for subsequent beam shaping assembly design are reliable. This work presents a comprehensive benchmarking of four widely-used Monte Carlo codes - Geant4, PHITS, FLUKA (CERN), and MCNP - for modelling low-energy neutron production target reactions. Using their recommended physics models and cross-section libraries, we evaluate each code’s performance in simulating four beam-target reactions: 7 Li(p,n) 7 Be, 9 Be(p,n) 9 B, 9 Be(d,n) 10 B, and C(d,n)N. Predictions of neutron yield, angular distributions, and energy spectra are compared against available thick target experimental data. Results show varying levels of agreement between the codes depending on the reaction type, energy range, and beam characteristics. Geant4, MCNP and PHITS are the overall best performing codes for the simulation of total neutron yield and yield in the forward direction across most reactions. Across energies where experimental benchmarks exist, inter-code discrepancies in total and forward-directed yield are typically 10 to 30%, with larger deviations at near-threshold incident ion energies. PHITS provides the best overall reproduction of experimental spectra, particularly for the 9 Be(p,n) 9 B reaction. Additionally, PHITS demonstrates superior computational performance for most reactions. These findings provide valuable guidance for ABNS design, highlighting the strengths and limitations of each code for the simulation of low-energy neutron production reactions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Thermal Stability of LiNi x Mn y Co z O 2 Cathode Materials

Here, the thermal evolution of LiNi x Mn y Co z O 2 (NMC) lithium-ion battery electrode materials is examined at various states of charge (SOC) or lithium concentrations for a variety of Ni:Mn:Co ratios or electrode compositions. Synchrotron X-ray diffraction (XRD) combined with Rietveld analysis shows the onset decomposition temperatures of phases, decomposition products, lattice parameters, and phase fractions as a function of composition and SOC. SOC impacts the lattice parameters of the NMC phase, where a collapse of the c-axis in the NMC phases is noted due to lithium extraction. Among the compositions examined, the low-Ni NMC111 uncycled sample (NMC111 0% SOC) exhibited the highest thermal stability, with a decomposition temperature approximately 250 °C higher than that of NMC532 0% and NMC811 0%. When the SOC exceeds 50% (i.e., more than 0.4 mol of Li ions extracted), the influence of Ni content on the decomposition temperature becomes negligible, with decomposition occurring around 250−300 °C for all compositions. Ni content also affects the decomposition pathways: NMC111 tends to first form a TM 3 O 4 -type phase, where TM represents transition metals, before transforming into a TMO-type phase, whereas most of the NMC811 samples directly decompose into the TMO phase. The presence of metallic phases was confirmed by both XRD and thermogravimetric-differential scanning calorimetry (TGA-DSC) analysis, as a result of heating under inert conditions. The TGA-DSC results suggest that metallic phase formation is favored at lower SOC in samples with a higher Ni content. This work provides comprehensive insight into the thermal degradation pathways of NMC materials as a function of composition, SOC, and temperature.

Peng, Jian [Univ. of New South Wales, Sydney, NSW ↗
DOE OSTI2025

Dynamic nanodomains dictate macroscopic properties in lead halide perovskites

Lead halide perovskites have emerged as promising materials for solar energy conversion and X-ray detection owing to their remarkable optoelectronic properties. However, the microscopic origins of their superior performance remain unclear. Here we show that low-symmetry dynamic nanodomains present in the high-symmetry average cubic phases, whose characteristics are dictated by the A-site cation, govern the macroscopic behaviour. We combine X-ray diffuse scattering, inelastic neutron spectroscopy, hyperspectral photoluminescence microscopy and machine-learning-assisted molecular dynamics simulations to directly correlate local nanoscale dynamics with macroscopic optoelectronic response. Our approach reveals that methylammonium-based perovskites form densely packed, anisotropic dynamic nanodomains with out-of-phase octahedral tilting, whereas formamidinium-based systems develop sparse, isotropic, spherical nanodomains with in-phase tilting, even when crystallography reveals cubic symmetry on average. We demonstrate that these sparsely distributed isotropic nanodomains present in formamidinium-based systems reduce electronic dynamic disorder, resulting in a beneficial optoelectronic response, thereby enhancing the performance of formamidinium-based lead halide perovskite devices. By elucidating the influence of the A-site cation on local dynamic nanodomains, and consequently, on the macroscopic properties, we propose leveraging this relationship to engineer the optoelectronic response of these materials, propelling further advancements in perovskite-based photovoltaics, optoelectronics and X-ray imaging.

Materials Science↗
DOE OSTI2025

Development, validation, and verification of multi-pass thermo-mechanical welding simulations using the open-source MOOSE framework: NeT TG4 benchmark weldment

This study develops and validates a sequentially coupled thermo-mechanical welding simulation for the three-pass 316L stainless steel NeT TG4 benchmark weldment using the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) and the Nuclear Engineering Material model Library (NEML). A diffused ellipsoidal heat source was calibrated against thermocouple data and weld macrographs to accurately model the fusion zone geometry and transient thermal fields. Material hardening is represented using the Lemaitre-Chaboche mixed isotropic-kinematic hardening model, while four annealing models - no annealing, single-stage at 1050 °C and 1300 °C, and two-stage at 800 °C/1300 °C - were implemented to assess the impact of annealing models on the accuracy of the predicted welding-induced plasticity, distortions, and residual stresses. The predictions were validated against experimental measurements and benchmarked against results from commercial software, demonstrating that thermo-mechanical MOOSE welding simulations achieve comparable accuracy with enhanced computational efficiency. This work highlights the potential of using open-source finite element frameworks like MOOSE for advanced manufacturing simulations.

Ji, Wendy [Australian Nuclear Science and Technolo↗
DOE OSTI2025

Probing the shear-induced microstructure of a smectite clay aqueous suspensions using rheo-USANS and rheo-SIPLI measurements

Hypothesis: The static microstructure of aqueous sodium-montmorillonite (Na-Mt) suspensions at low ionic strengths (where Particle Size/Debye Length ≈1) exhibits both the particle–particle ordering as well as aggregation with repulsive ordered domains having characteristic optical birefringence and attractive aggregated entities larger than 20 µm resulting in ever-increasing yield stresses also known as physical aging-rejuvenation behavior. We hypothesize that the attractive particle–particle aggregation is the underlying cause behind the physical aging-rejuvenation behavior observed in Na-Mt suspensions with no contribution from structural dynamics driven by repulsive particle–particle ordering or jamming. Experiments: We investigate the shear-induced microstructure of aqueous Na-Mt suspensions in the sol and gel state using rheo-ultra-small angle neutron scattering (rheo-USANS) experiments at shear rates of 1, 50, 500, and 2000 s −1 . We also perform rheo-shear-induced polarization light imaging (rheo-SIPLI) experiments to relate ordering with shearing and aging. Findings: Shearing the suspensions at low to moderate shear rates induces particle–particle aggregation and shearing at high shear rates induces the breakage of particle–particle aggregation in the sol and gel states, suggesting the microstructural aggregation in the sol and gel state is shear sensitive and a full rejuvenation or breakage of particle–particle aggregation is only achieved at a minimum critical shear rate. The rheo-SIPLI experiments reveal that the sol and gel state exhibited strong Maltese cross patterns at a shear rate of 1000 s −1 , indicating particle–particle ordering. Post shearing, the gel exhibited temporal evolution of storage modulus without any noticeable influence on the appearance of the Maltese cross pattern indicating physical aging and particle ordering are distinct length scale phenomena in Na-Mt suspensions and the physical aging-rejuvenation behaviour is a feature of particle–particle aggregation as opposed to ordering.

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

Room temperature sodium diffusion in sodium iron phosphate investigated with neutron scattering

Sodium iron phosphate, Na x FePO 4 , is currently one of the most promising cathode materials for sodium ion batteries. Here we employ a systematic approach using three neutron scattering techniques to investigate temperature dependent properties with a focus on the charge transfer processes. Two distinct Na + diffusive motions are observed with quasielastic neutron scattering, including room temperature motion. These processes correspond to rapid double jumps through extrinsic vacancy sites of ∼6 Å along with the more fundamental nearest neighbor jumps of 3.1 Å. Experimental observations of room temperature ionic diffusion are highly desirable, as the diffusion pathways are naturally well suited for efficient transport. Inelastic neutron scattering was used to probe the vibrational density of states that assist in the ionic and electronic (polaronic) transport processes, which are further supported with density functional theory. Neutron diffraction examines the disordering of the sodium sublattice upon transition to the solid-solution phase, which is believed to enhance the overall kinetics during sodium insertion/extraction. Finally, 57 Fe Mössbauer spectroscopy measurements probes polaron hopping and the local structural evolution that modifies the electron density surrounding the redox active iron sites. This study sheds light on the temperature induced atomic scale dynamics that occur in olivine Na x FePO 4 .

Novak, Eric [Swarthmore College, Swarthmore, PA (U↗
DOE OSTI2025

Assessment of irradiation embrittlement effect on fatigue life of a pressurized-water reactor pressure vessel using the fracture toughness master curve approach

The reactor pressure vessel (RPV) is a critical structural component in pressurized-water reactors, and it is designed to withstand extreme conditions, such as high pressures, elevated temperatures, and prolonged radiation exposure. Ensuring RPV integrity is essential for the safe and reliable long-term operation of nuclear power plants, especially as aging mechanisms such as fatigue and irradiation embrittlement pose increased risks. Fatigue, caused by cyclic thermal and mechanical loading, can lead to crack initiation in localized high-stress regions. Simultaneously, neutron irradiation, particularly in the beltline region, progressively reduces fracture toughness, increasing susceptibility to brittle fracture. These combined effects of fatigue and irradiation embrittlement potentially impact the RPV structural integrity, necessitating fitness-for-service assessments. This study applies the fracture toughness master curve approach to evaluate the impact of irradiation embrittlement on RPV fatigue life. A coupled thermo-mechanical stress analysis identifies critical stress locations under normal service transients, pinpointing regions most vulnerable to fatigue crack initiation and growth. Stress intensity factors for postulated flaws at these locations are calculated, enabling an assessment of fatigue life under irradiated and unirradiated conditions. The results indicate that neutron irradiation embrittlement accelerates the conditions in which a critical crack can form and lead to failure, particularly at lower temperatures. The failure occurs where reduced fracture toughness limits the material’s resistance to crack growth. Axial cracks at nozzle corners are the most life limiting without irradiation, while circumferential cracks demonstrate longer fatigue lives. The findings highlight the importance of incorporating irradiation effects into fatigue life predictions to ensure the long-term structural integrity of RPVs.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Reaction-Induced Fracturing of Porous Carbonate Rocks during Volume-Increasing Replacement by Witherite

Volume-increasing replacement reactions can lead either to fracturing of the parent phase or to the formation of cohesive layers that passivate further reaction, but the factors that drive one outcome or the other are not understood. In this experimental study, we investigated the volume-increasing replacement of carbonate rocks (consisting of CaCO 3 and CaMg(CO 3 ) 2 ) of different porosity by witherite (BaCO 3 ). Samples were characterized using scanning electron microscopy, Raman spectroscopy, small-angle neutron/X-ray scattering, X-ray tomography, and scanning transmission electron microscopy. We observed the formation of a witherite reaction layer and witherite formation within pores and along grain boundaries. Despite this being a volume-increasing replacement reaction, newly formed witherite was porous, potentially allowing further replacement. Filled fractures were observed in the low-porosity carbonates, whereas witherite formed within pores in high-porosity carbonates. We conclude that fracturing of the parent phase versus passivation is contingent on the initial microstructure of the parent with an optimal degree of porosity required for fracturing.

Calcite↗
DOE OSTI2026

Spin Texture Control and Magnetic Gap Engineering in a Ferromagnetic Insulator–Topological Insulator Sandwiched Heterostructure

Quantum materials that combine magnetism with topological order are emerging as key platforms for next-generation spintronics and low-energy electronics. They enable the realization of emergent quantum phenomena, such as the quantum anomalous Hall effect and axion insulator states. The ferromagnetic insulator (FMI)/topological insulator (TI)/FMI sandwich structure of a single-septuple layer (1SL) MnBi 2 Te 4 /four-quintuple layer (4QL) Bi 2 Te 3 /1SL MnBi 2 Te 4 holds great potential to achieve such desirable quantum phenomena at an elevated temperature, owing to its large Dirac point band gap and high Curie temperature. Here, for this work, spin- and angle-resolved photoemission spectroscopy (spin-ARPES) is employed to directly verify that the band gap arises from broken time-reversal symmetry via proximity-driven magnetization. This study demonstrates direct control of the spin state via external magnetic fields and unambiguously confirms the exchange interaction as the gap-opening mechanism. The robust magnetic gap and controllable spin texture make this heterostructure a suitable candidate for spintronic applications and magnetic topological quantum phases.

exchange gap↗
DOE OSTI2025

Bose–Einstein condensation of a two-magnon bound state in a spin-1 triangular lattice

In ordered magnets, the elementary excitations are spin waves (magnons), which obey Bose–Einstein statistics. Similarly to Cooper pairs in superconductors, magnons can be paired into bound states under attractive interactions. The Zeeman coupling to a magnetic field is able to tune the particle density through a quantum critical point, beyond which a ‘hidden order’ is predicted to exist. Here, in this work, we report direct observation of the Bose–Einstein condensation of the two-magnon bound state in Na 2 BaNi(PO 4 ) 2 . Comprehensive thermodynamic measurements confirmed the two-dimensional Bose–Einstein condensation quantum critical point at the saturation field. Inelastic neutron scattering experiments were performed to establish the microscopic model. An exact solution revealed stable two-magnon bound states that were further confirmed by electron spin resonance and nuclear magnetic resonance experiments, demonstrating that the quantum critical point is due to the pair condensation, and the phase below the saturation field is likely the long-sought-after spin nematic phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Accelerating ion transport by dynamic asymmetry of alternating polymer electrolytes

Polymer based electrolytes allow the absence of volatile components in batteries thus increasing their safety. Yet, they exhibit drawbacks based on their low conductivity. We have used an alternating polymer consisting of dimethyl siloxane (DMS) and ethylene glycol (EG) blocks to circumvent known disadvantages of the usually used polyethylene glycol (PEG). Incorporating dimethyl siloxane lowers the glass-transition temperature and thus reduces the segmental relaxation time, by dynamic asymmetry or internal plasticization of the constituting polymer blocks. The alternating structure ensures miscibility of the different components and hinders crystallization. Furthermore, the pure polymer, P(DMS 3 -alt-EG 4 ), shows a segmental relaxation time well in the range needed for polymer electrolytes. Mixtures of LiClO 4 and P(DMS 3 -alt-EG 4 ) show a drastically reduced temperature dependence of their DC conductivity in comparison to PEG based systems, resulting in an increase by two orders of magnitude at T = 5 °C and even three to four orders of magnitude at T = 0 °C. Addition of coordinating (acetonitrile) or non-coordinating (toluene) solvent increases conductivity either via additional plasticization or by weakening the Li-binding yet looking at the dynamics at low concentrations of additional solvent the mobility of the polymer is reduced. In conclusion, the solvent addition leads only at higher solvent concentration to a reduction in relaxation time.

Jakobi, Bruno [Louisiana State Univ., Baton Rouge,↗
DOE OSTI2025

Accelerated and heterogeneous corrosion of Cr-doped uranium nitride fuel pellets

The steam oxidation of Cr-doped UN fuel pellets is analyzed during sequential isothermal holds up to 720 °C. In situ neutron diffraction results show how Cr is accommodated in a secondary U 2 CrN 3 phase, leading to the formation of a duplex UN/U 2 CrN 3 microstructure. Under corrosion, the oxidation of the two phases begins at 400 °C for UN and 430 °C for U 2 CrN 3 , respectively. Because the UN phase is preferentially oxidized in the presence of U 2 CrN 3 , addition of Cr in UN based nuclear fuel is found to accelerate the corrosion rate. At 430 °C the oxidation of UN in the UN/U 2 CrN 3 microstructure is ~ 5 times faster than pure UN, increasing to ~19 times faster at 460 °C. The oxidation of U 2 CrN 3 produces UO 2 via the formation of two transient intermediate phases. In situ neutron diffraction enables oxidation processes of UN and U 2 CrN 3 components to be followed separately within the two-phase system.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗