Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel performance simulation”

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 325 records · Page 18

Realistic fuel spray modeling for gasoline direct injection engine applications

Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.

33 ADVANCED PROPULSION SYSTEMS↗

Modeling a Thermal-Spectrum LEU-fueled Molten Salt Reactor Co-fueled with Thorium in SCALE 6.3.1 and Serpent-2

Interest in the development of advanced nuclear fission reactors for commercial electricity production has risen in recent years, with Generation IV reactor designs offering numerous advantages in safety, efficiency, fuel cycle sustainability, and waste management. In particular, thorium-fueled molten salt reactors (MSRs) are considerably promising for enhancing fuel cycle sustainability in their ability to breed fissile 233U fuel from thorium, a presently untapped and widely abundant resource. A novel MSR fuel cycle concept, the ”Sourdough” refueling and waste management strategy, has previously been demonstrated with a traditional uranium-based fuel cycle in a thermal-spectrum MSR operating with low enriched uranium (LEU) fuel with favorable neutronic performance. However, the ability to use this unique fuel cycle approach with thorium-based molten salt fuels has not yet been studied. In this work, the Sourdough fuel cycle was implemented in a small, thermal-spectrum MSR fueled with high assay low enriched uranium (HALEU) and fertile 232Th for breeding 233U fuel. Relevant neutronic data, including fuel and isothermal temperature feedback behavior, was studied using the SCALE 6.3.1 and Serpent-2 code systems, and keff data was measured during simulated depletion at 400 MWth. The Sourdough fuel cycle concept is shown to perform favorably with a thorium-fueled MSR model, thus warranting further study into its use in other MSR designs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Determination of the spent fuel decay heat with the VERA core simulator

VERA core simulator capabilities to predict spent fuel decay heat have been explored. The comparison with the Serpent2 and NEWT/TRITON sequence from the SCALE package, performed on the NPP Krško fuel test cases, has enabled independent verification of the successful implementation of the VERA depletion/decay calculation sequence applied on the extended nuclide set. Since the VERA is a 3-D core simulator, the focus in the rest of the paper was on the investigation of the 3-D effects and non-linearity of the soluble boron concentration averaging process, which are difficult to assess with most of available computational tools. A difference in the decay heat prediction between the 3-D and 2-D model approaches was analyzed. In addition, the impact of the Inconel grids on the decay heat was determined. Averaging the concentration of soluble boron over the fuel irradiation period has been found to be a viable approach, since it results in slight decay heat overprediction, which is considered conservative for most applications. Finally, demonstration on the Watts Bar Unit 1 cycle 1 depletion has shown formidable VERA capabilities to accurately predict spent fuel decay heat.

Kromar, Marjan↗

GOTHIC Aerosol Source Depletion Studies

Pacific Northwest National Laboratory has continued work to develop an aerosol-laden flow modeling capability with the Generation of Thermal Hydraulic Information in Containment (GOTHIC TM ) computer code to perform simulations for thermal hydraulic conditions and aerosol transport and deposition in spent fuel casks. This report describes our recent work to expand our model that was originally developed in 2019 to allow for thermal characterization, carrier gas flow determination, and tracking of particulate behavior throughout the entire canister volume. This was achieved through conversion of decay heat source in all fuel tubes within the model from heaters to thermal conductors in GOTHIC as well as remeshing the internal volume of the canister. These model improvements allowed for simulations of base case scenarios for comparison of similar efforts by our collaborators with alternative tools at Sandia and Oak Ridge. Several characteristics of the GOTHIC code were elucidated by this effort and were documented for consideration in how the code is used for this effort moving forward. For instance, treatment of particle size distributions, initial spatial distributions, mesh gradients, and minimum volume concentrations were identified as important issues for consideration. Ultimately, the current version of the code is capable of tracking temperatures, flow rates, and particle behavior throughout the canister internal volume. This report presents predicted depletion times and preferred deposition patterns within the canister. Additionally, this report documents code performance and insights for further model development and use for comparisons to other codes and for predicting system behavior in experiments being performed and planned by our collaborators.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A simulation-based parametric study of CLOU chemical looping reactor performance

Chemical looping with oxygen uncoupling (CLOU) is a variant on chemical looping combustion in which the oxygen carrier releases gaseous O 2 in the fuel reactor, making it well-suited for solid fuels, since the released gaseous O 2 readily reacts with solid char. This work presents several computational fluid dynamic (CFD) simulations of copper-based CLOU in a dual fluidized bed system, each with different operating conditions. The modeling predicted that coal particle sizes as large as 1000 μm did not significantly affect performance. Increased oxygen carrier copper loading resulted in an excess of gaseous oxygen in the product gas stream. Decreasing the oxygen carrier bed mass as well as reducing the air reactor fluidizing velocity did not supply enough oxygen to the fuel reactor to complete combustion of the coal. This generated a failure state in which the temperature continued to decrease in the fuel reactor from the lack of combustion, which in turn reduced the O 2 equilibrium partial pressure, further lowering the amount of combustion possible. Sufficient O 2 can be maintained in the fuel reactor by ensuring a high enough air reactor velocity and a large enough supply of oxygen carrier inventory to handle the chosen coal feed rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of Reduced Order Model for HT-9 Creep and Modifications to Current HT-9 Creep Model in BISON

In order to leverage existing reduced order models (ROM) for the modeling of LANL-developed HT-9 cladding constitutive behavior, mechanistic-based models have been implemented into BISON. It is posited that the increased fidelity of the ROM will show improvement over engineering-scale models in terms of comparison against experimental measurements. INL will assist LANL in implementing an HT9 ROM in BISON. Once the HT-9 mechanistic constitutive model can be leveraged in nuclear performance simulations, the numerical results of fast reactor models will be analyzed and compared against the Fuels Irradiation & Physics Database (FIPD) and separate effects tests to attempt a validation process. INL and LANL will work together to perform any necessary improvement that is identified during the implementation and use processes. Successful completion of the milestone will enable the general use of validated mechanistic ROMs for HT-9 cladding in BISON. This report provides a brief introduction to the Los Alamos Reduced Order Model Applied to Nonlinear Constitutive Equations (LAROMANCE) code, mechanical testing used to illustrate its creep and plastic deformation predictions for HT-9 cladding, and a comparison with existing models implemented into BISON. Current models have been modified to more appropriately account for primary thermal creep and results for corresponding testing are included in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulated Service and Stress Corrosion Cracking Testing for Friction Stir Welded Spun Form Domes

Damage tolerance testing development was required to help qualify a new spin forming dome fabrication process for the Ares 1 program at Marshall Space Flight Center (MSFC). One challenge of the testing was due to the compound curvature of the dome. The testing was developed on a sub-scale dome with a diameter of approximately 40 inches. The simulated service testing performed was based on the EQTP1102 Rev L 2195 Aluminum Lot Acceptance Simulated Service Test and Analysis Procedure generated by Lockheed Martin for the Space Shuttle External Fuel Tank. This testing is performed on a specimen with an induced flaw of elliptical shape generated by Electrical Discharge Machining (EDM) and subsequent fatigue cycling for crack propagation to a predetermined length and depth. The specimen is then loaded in tension at a constant rate of displacement at room temperature until fracture occurs while recording load and strain. An identical specimen with a similar flaw is then proof tested at room temperature to imminent failure based on the critical offset strain achieved by the previous fracture test. If the specimen survives the proof, it is then subjected to cryogenic cycling with loads that are a percentage of the proof load performed at room temperature. If all cryogenic cycles are successful, the specimen is loaded in tension to failure at the end of the test. This standard was generated for flat plate, so a method of translating this to a specimen of compound curvature was required. This was accomplished by fabricating a fixture that maintained the curvature of the specimen rigidly with the exception of approximately one-half inch in the center of the specimen containing the induced flaw. This in conjunction with placing the center of the specimen in the center of the load train allowed for successful testing with a minimal amount of bending introduced into the system. Stress corrosion cracking (SCC) tests were performed using the typical double beam assembly and with 4-point loaded specimens under alternate immersion conditions in a 3.5% NaCl environment for 90 days. In addition, experiments were conducted to determine the threshold stress intensity factor for SCC (K1SCC) of Al-Li 2195 which to our knowledge has not been determined previously. The successful simulated service and stress corrosion testing helped to provide confidence to continue to Ares 1 scale dome fabrication.

Stewart, Thomas J.↗

Hydride Reprecipitation in Prototypical 17x17 PWR Fuel Rods after Simulated Vacuum Drying Conditions

Light water reactor (LWR) fuel rod cladding picks up hydrogen as it oxidizes in reactor. Once the hydrogen concentration exceeds the solubility limit, zirconium hydride platelets precipitate, generally in the circumferential direction. In preparation for dry storage, elevated temperatures can dissolve some of the circumferential hydrides, which during cooling can then reprecipitate in a direction perpendicular to the hoop stress (radially). From 2019 to 2024, ORNL published a series of reports that included metallographic observations of cladding after a thermal transient to 400°C followed by slow cooling. The simulated dry storage process was performed on three full-length 17x17 PWR fuel rods at their as-discharged rod internal pressure. There were two important observations from the metallographic examinations of these heat-treated rods: (1) the propensity for hydride reorientation is strongly influenced by nearby as-discharged precipitated hydrides, (2) the radial hydrides observed near the cladding inner diameter were often associated with cracks in the pellet. Subsequent modeling using the BISON fuel performance code that considered the effects of pellet cracks and pellet-cladding bonding observed in high-burnup fuel support the conclusion that pellet cracks enhance the local hoop stress in the cladding creating preferential sites for hydride reorientation, and that the influence of the pellet cannot be ignored in accurate fuel rod performance predictions.

Montgomery, Rose [ORNL] (ORCID:0000000286038936)↗

Energy Consumption and Cost Reduction of Future Light-Duty Vehicles through Advanced Vehicle Technologies: A Modeling Simulation Study Through 2050

The U.S. Department of Energy’s (DOE’s) Vehicle Technologies Office (VTO) and Hydrogen and Fuel Cell Technologies Office (HFTO) aim to develop sustainable, affordable and efficient technologies for transportation of goods and people. Translating investments in advanced transportation component technologies and powertrains to estimate vehicle-level fuel savings potential is critical for understanding DOE’s impact. In this work, we simulated technologies funded by VTO and HFTO for light duty vehicles. The simulations were performed across: Multiple powertrain configurations (i.e., conventional, power-split, extended-range electric vehicle, battery electric drive, and fuel-cell vehicles), Vehicle classes (i.e., compact car, midsize car, small sport utility vehicle [SUV], midsize SUV, and pickup trucks); and Fuels (i.e., gasoline, diesel, hydrogen, and battery electricity). These various technologies are assessed for six different timeframes: laboratory years 2015, 2020, 2025, 2030, and 2045. A delay of 5 years is assumed between laboratory year and model year (year technology is introduced into production). Finally, uncertainties are included for both technology performance and cost aspects by considering two cases: Low case, aligned with DOE technology manager estimates of expected original equipment manufacturer (OEM) improvements based on regulations, business as usual; and High case, aligned with aggressive technology advancements based on R&D targets developed through support by VTO & HFTO. These scenarios are not intended as predictions of future performances. The energy and cost impact of different technologies were estimated using Autonomie (www.autonomie.net), Argonne vehicle system simulation tool. Autonomie is a state-of-the-art vehicle system simulation tool used to assess the energy consumption, performance and cost of multiple advanced vehicle technologies across classes (from light to heavy duty), powertrains (from conventional to HEVs, FCEVs, PHEVs and BEVs), components and control strategies. Autonomie is packaged with a complete set of vehicle models for a wide range of vehicle classes, powertrain configurations and component technologies, including vehicle level and component level controls. These controls were developed and calibrated using dynamometer test data. Autonomie has been used to support a wide range of studies including analyzing various component technologies, sizing powertrains components for different vehicle requirements, comparing the benefits of powertrain configurations, optimizing both heuristic and route based vehicle energy control and predicting transportation energy use when paired with a traffic modeling tool such as POLARIS. This report documents the assumptions and estimates the vehicle-level energy consumption benefits and associated technology costs for the various types of light duty vehicles. All details of vehicle assumptions and simulation results are available in the spreadsheets accompanying this report.

33 ADVANCED PROPULSION SYSTEMS↗

In-House Developed Multiphysics Simulation for the Performance of Solid Oxide Cells (SOCs)

Reversible solid oxide cells(rSOCs) are an enabling energy storage/production technology for a dynamic grid environment. Reversible operation requires strong working knowledge of fuel cell (SOFC) and electrolyzer (SOEC). Performance degradation of SOECs has been observed; however, the details of physical processes related to the performance degradation remain unknown. Multiphysics simulations were performed to investigate the performance degradation of solid oxide electrolysis cells under various working conditions.

Yang, Tao↗

Advanced Turboprop Transport Aircraft Modeling for the Electrified Powertrain Flight Demonstration Project

The objective of this study is to establish models for state of the art, representative regional non-electrified turboprop configurations in the 19 to 50 passenger capability range and assess performance impacts from advanced aircraft technologies. Two of these models will serve as reference baseline configurations for comparison to electrified aircraft propulsion (EAP) demonstrator concepts for the Electrified Powertrain Flight Demonstration (EPFD) project. The other two non-electrified configurations will primarily assess the impacts of advanced technologies on fuel consumption and CO₂ emissions. This paper details an aircraft synthesis and performance analysis approach using the General Aviation Synthesis Program (GASP) which is used to establish baseline performance capabilities, calibrate airplane and engine models using published manufacturer data and conduct technology sensitivity studies. Advanced technologies are selected from the categories of propulsion, aerodynamics, and flight systems. A GASP-Monte Carlo simulation framework for uncertainty propagation is used to obtain performance distributions of the expected fuel burn, weight reduction and CO₂ emissions reductions from advanced aircraft technologies. Significant fuel burn savings from application of aerodynamics and propulsion technologies was observed.

Dahlia D.V. Pham↗

Effects of Internal and External Heat Sources on Cladding Microstructure and Rupture Performance

In the event of a Loss of Coolant Accident (LOCA), the primary supply of cooling water for a nuclear reactor is lost, leading to a significant pressure differential across the cladding wall. Without adequate cooling, the fuel rods continue to heat as a result of fission reactions. Research at Oak Ridge National Laboratory’s (ORNL’s) Severe Accident Test Station (SATS) is currently focused on evaluating fuel cladding performance using an external infrared lamp as a heat source, whereas legacy testing primarily utilized an internal heating approach. While external heating may better simulate the effect of neighboring fuel rods heating a central rod, it may not accurately represent the internal heat absorption from the fuel during accident transients. The heating dynamics depend greatly on the fuel rod's position within the bundle and the reactor. To fully understand the implications of a LOCA event and assess the influence of internal heating on cladding performance, combined internal heating and pressurization capability was developed at ORNL. Tests were conducted to compare cladding segments heated internally (representing heat from the fuel within the rod) and externally (representing heat from adjacent fuel rods). The findings indicate that both internal and external heating result in comparable rupture temperatures during 5°C/s laboratory LOCA tests and also agree with legacy test data. Axial temperature gradients and internal heat source dispersal were found to significantly impact cladding deformation and rupture geometry. Clear modifications were outlined to further improve the capability with heating rates above 5°C/s.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effects of Internal and External Heat Sources on Cladding Microstructure and Rupture Performance

In the event of a Loss of Coolant Accident (LOCA), the primary supply of cooling water for a nuclear reactor is lost, leading to a significant pressure differential across the cladding wall. Without adequate cooling, the fuel rods continue to heat as a result of fission reactions. Research at Oak Ridge National Laboratory’s (ORNL’s) Severe Accident Test Station (SATS) is currently focused on evaluating fuel cladding performance using an external infrared lamp as a heat source, whereas legacy testing primarily utilized an internal heating approach. While external heating may better simulate the effect of neighboring fuel rods heating a central rod, it may not accurately represent the internal heat absorption from the fuel during accident transients. The heating dynamics depend greatly on the fuel rod's position within the bundle and the reactor. To fully understand the implications of a LOCA event and assess the influence of internal heating on cladding performance, combined internal heating and pressurization capability was developed at ORNL. Tests were conducted to compare cladding segments heated internally (representing heat from the fuel within the rod) and externally (representing heat from adjacent fuel rods). The findings indicate that both internal and external heating result in comparable rupture temperatures during 5°C/s laboratory LOCA tests and also agree with legacy test data. Axial temperature gradients and internal heat source dispersal were found to significantly impact cladding deformation and rupture geometry. Clear modifications were outlined to further improve the capability with heating rates above 5°C/s.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Achieving record hot spot energies with large HDC implosions on NIF in HYBRID-E

HYBRID-E is an inertial confinement fusion implosion design that increases energy coupled to the hot spot by increasing the capsule scale in cylindrical hohlraums while operating within the current experimental limits of the National Ignition Facility. HYBRID-E reduces the hohlraum scale at a fixed capsule size compared to previous HYBRID designs, thereby increasing the hohlraum efficiency and energy coupled to the capsule, and uses the cross-beam energy transfer (CBET) to control the implosion symmetry by operating the inner (23° and 30°) and outer (44° and 50°) laser beams at different wavelengths (Δλ> 0). Small case to capsule ratio designs can suffer from insufficient drive at the waist of the hohlraum. We show that only a small amount of wavelength separation between the inner and outer beams (Δλ1-2 Å) is required to control the symmetry in low-gas-filled hohlraums (0.3 mg/cm3 He) with enough drive at the waist of the hohlraum to symmetrically drive capsules 1180 μm in outer radius. This campaign is the first to use the CBET to control the symmetry in 0.3 mg/cm3 He-filled hohlraums, the lowest gas fill density yet fielded with Δλ> 0. We find a stronger sensitivity of hot spot P2 in μm per Angstrom (40–50 μm/Å wavelength separation) than observed in high-gas-filled hohlraums and previous longer pulse designs that used a hohlraum gas fill density of 0.6 mg/cm3. There is currently no indication of transfer roll-off with increasing Δλ, indicating that even longer pulses or larger capsules could be driven using the CBET in cylindrical hohlraums. We show that the radiation flux symmetry is well controlled during the foot of the pulse, and that the entire implosion can be tuned symmetrically in the presence of the CBET in this system, with low levels of laser backscatter out of the hohlraum and low levels of hot electron production from intense laser–plasma interactions. Radiation hydrodynamic simulations can accurately represent the early shock symmetry and be used as a design tool, but cannot predict the late-time radiation flux symmetry during the peak of the pulse, and semi-empirical models are used to design the experiments. Deuterium–tritium (DT)-layered tests of 1100 μm inner radius implosions showed performance close to expectations from simulations at velocities up to ~360 km/s, and record yields at this velocity, when increasing the DT fuel layer thickness to mitigate hydrodynamic mixing of the ablator into the hot spot as a result of defects in the ablator. However, when the implosion velocity was increased, mixing due to these defects impacted performance. The ratio of measured to simulated yield for these experiments was directly correlated with the level of observed mixing. These simulations suggest that reducing the mixing, e.g., by improving the capsule defects, could result in higher performance. In addition, future experiments are planned to reduce the coast time at this scale, delay between the peak compression and the end of the laser, to increase the hot spot convergence and pressure. To reduce the coast time by several hundred ps compared to the 1100 μm inner radius implosions, HYBRID-E has also fielded 1050 μm inner radius capsules, which resulted in higher hot spot pressure and a fusion energy yield of ~170 kJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detailed Simulation Datasets Quantifying U.S. DOE VTO/HFTO R&D Benefits Across Light- to Heavy-Duty Vehicles

For more than 20 years, Argonne National Laboratory’s Vehicle & Mobility Systems Department has assessed how R&D investments by the U.S. Department of Energy’s Transportation Technologies Office and Alternative Fuels and Feedstocks Office affect vehicle energy use and cost. The analyses are performed using Autonomie, Argonne’s full-vehicle simulation tool for energy consumption, performance, and cost. The study covers five time frames ranging from present day through 2050, with more than 30 vehicle classes and applications (10 light duty and >20 medium and heavy duty), as well as six powertrain configurations (conventional, start-stop, hybrid electric vehicle, plug-in hybrid electric vehicle, battery-electric vehicle, and fuel cell electric vehicle) and five fuels (gasoline, diesel, natural gas, hydrogen, and electricity). Low and high technology uncertainty scenarios have been considered to capture a realistic range of outcomes. The resulting datasets include the assumptions used (i.e., efficiency, $/kWh), vehicle-level data (power, energy, weight, and cost), and outputs such as energy consumption, manufacturer’s suggested retail price, and total cost of ownership. These data are critical to stakeholders working in transportation, technology assessment, and long-term R&D planning.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

BISON validation to in situ cladding burst test and high-burnup LOCA experiments

The process for developing and qualifying nuclear fuels for commercial nuclear application requires fundamental material development, characterization, and design; out-of-pile testing on unirradiated materials; integral fuel rod irradiations, testing, and postirradiation examinations; and transient analyses. The historical approach depends on the generation of large empirical datasets and series of integral fuel rod irradiations, and this approach ultimately takes ~20 years—or sometimes longer—to acquire data through extensive sequential testing. Thus, the qualification and eventual deployment of new fuel systems constitute a long process. However, recent technological advancements have provided researchers the opportunity to perform out-of-cell, in situ measurements to assess material performance for the duration of the experiment. One such example of this capability is the use of digital image coordination and thermal imaging to assess Zircaloy cladding performance under a simulated loss-of-coolant accident (LOCA) transient condition. In situ measurements generally provide high-fidelity strain, strain rates, and temperature surface maps. This is critical for the US nuclear industry, which is actively developing a technical basis to support extending the peak rod average burnup from 62 to ~75 GWd/tU and the deployment of accident-tolerant fuel. However, the US Nuclear Regulatory Commission (NRC) outlined in its research information letter several technical issues that the industry must address before extending burnup. One topic of specific interest is understanding the cladding balloon and rupture geometry during the LOCA heat-up phase. By leveraging these advanced in situ capabilities, this work used in situ data generated from a simulated LOCA to better understand high-temperature creep and its effect on Zircaloy balloon and rupture performance. Here, this work used the BISON fuel performance code to assess the high-temperature creep model predictions with in situ data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Fully Explicit Integrator for Modeling Astrophysical Reactive Flows

Simulating complex astrophysical reacting flows is computationally expensive—reactions are stiff and typically require implicit integration methods. The reaction update is often the most expensive part of a simulation, which motivates the exploration of more economical methods. In this research note, we investigate how the explicit Runge–Kutta–Chebyshev (RKC) method performs compared to an implicit method when applied to astrophysical reactive flows. These integrators are applied to simulations of X-ray bursts arising from unstable thermonuclear burning of accreted fuel on the surface of neutron stars. We show that the RKC method performs with similar accuracy to our traditional implicit integrator, but is more computationally efficient when run on CPUs.

GPU computing↗