Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code verification”

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 595 records · Page 33

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Summary of FY2024 Experimental Results to Support Development of New Inelastic Material Models and Validation of Section III, Division 5, Class A Design Rules

To address the data gap crucial for updating existing viscoplastic constitutive material models for the Class A materials in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code, Section III, Division 5, High Temperature Reactors (ASME, 2023) for inelastic design analysis, Oak Ridge National Laboratory (ORNL) conducted experimental studies on Alloy 800H and Alloy 617 across a range of temperatures up to their maximum temperature limits in Section III, Division 5. The studies aim to characterize the materials' deformation behavior under both strain-controlled mechanical cyclic loading and thermal cycling conditions. This report summarizes ORNL's FY 2024 experimental findings on Alloy 800H and Alloy 617, focusing on pure fatigue tests, cyclic stress-strain curves, thermomechanical fatigue experiments, and verification of the high temperature cyclic damage summation design rules.

36 MATERIALS SCIENCE↗

Multiphysics Analysis of the MSRE Experiment Using Griffin-SAM Coupled Code System

Molten salt reactors (MSRs) with flowing fuel have the unique feature of utilizing the fuel salt for heat generation and extraction at the same time since the fuel salt is circulating through the whole primary loop of the reactor. This movement of the fuel salt results in a partial decay of the delayed neutron precursors (DNPs) outside the core and corresponding redistribution in the active core region. To capture this phenomenon accurately, neutronics and thermal hydraulics computational tools need to be modified to handle the movement of the DNPs and their decay. Idaho and Argonne National Laboratories are actively working on developing neutronics and thermal hydraulics tools to model MSRs with flowing fuel. The neutronics code Griffin and the thermal hydraulics code SAM were extended to handle flowing fuel with the drift of the DNPs under the MOOSE framework. In this work, a simplified Multiphysics model of the Molten Salt Reactor Experiment (MSRE) is developed and utilized to perform steady-state and transient analyses. The MSRE experiments were designed to show the impact of the DNPs losses on core reactivity. The following section provides a description of the MSRE experiment with the core main parameters used to develop the model. Then, the developed Multiphysics model of the MSRE is presented. Finally, steady-state verification tests and unprotected loss of flow transient tests are discussed followed by a summary and conclusions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Serpent - Bison - THM Preliminary Multiphysics Modeling of a Nuclear Thermal Propulsion System Fuel Assembly

This work demonstrates the Monte Carlo neutronic and Thermo-Hydraulic coupling scheme using the Serpent code and MOOSE application Bison and Thermo Hydraulic Module. The coupling scheme is then applied to the reference BWX Technologies Nuclear Ther- mal Propulsion system at he fuel assembly level where it’s used to perform an analysis of the isothermal material coefficients and potential material reactivity worth. A method is developed to isolate which feedback effects should be considered for proceeding with reduced order deterministic neutronic modeling where branch off analysis must be con- ducted. The convergence behavior of the coupling scheme is demonstrated where it fol- lows the standard Picard iteration approach. Verification studies for the method of deduc- ing relevant feedback effects is also demonstrated.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Particulate Fuel Modeling of MC 2 -3 using Iterative Local Spatial Self-shielding Method

We report a new spatial self-shielding method for particulate fuels has been developed based on disadvantage factors and implemented in the MC 2 -3 code. This method named the iterative local spatial self-shielding (ILSS) method considers the shadowing effect of randomly distributed particles on spatial self-shielding in particles through a homogenized composition region added outside the particle of interest at the center. The self-shielded cross sections of the central particle are determined iteratively since they are used in determining the cross sections of the homogenized composition region. The ILSS method was verified for infinite stochastic medium problems of single and multiple types of particles, VHTR unit cell problems, and HTTR assembly problems. The verification test results show that the ILSS method accurately predicts the stochastic particle shadowing effect and reaction rates in particles, whereas the regular array model and the stochastic collision probability method underpredict the particle shadowing effect and overestimate reaction rates in particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗

Free-boundary MRxMHD equilibrium calculations using the stepped-pressure equilibrium code

Here the stepped-pressure equilibrium code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 , 112 502) is extended to enable free-boundary multi-region relaxed magnetohydrodynamic (MRxMHD) equilibrium calculations. The vacuum field surrounding the plasma inside an arbitrary 'computational boundary', , is computed, and the virtual-casing principle is used iteratively to compute the normal field on so that the equilibrium is consistent with an externally produced magnetic field. Recent modifications to SPEC are described, such as the use of Chebyshev polynomials to describe the radial dependence of the magnetic vector potential, and a variety of free-boundary verification calculations are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efforts to enhance reproducibility in a human performance research project

Background: Ensuring the validity of results from funded programs is a critical concern for agencies that sponsor biological research. In recent years, the open science movement has sought to promote reproducibility by encouraging sharing not only of finished manuscripts but also of data and code supporting their findings. While these innovations have lent support to third-party efforts to replicate calculations underlying key results in the scientific literature, fields of inquiry where privacy considerations or other sensitivities preclude the broad distribution of raw data or analysis may require a more targeted approach to promote the quality of research output. Methods: We describe efforts oriented toward this goal that were implemented in one human performance research program, Measuring Biological Aptitude, organized by the Defense Advanced Research Project Agency's Biological Technologies Office. Our team implemented a four-pronged independent verification and validation (IV&V) strategy including 1) a centralized data storage and exchange platform, 2) quality assurance and quality control (QA/QC) of data collection, 3) test and evaluation of performer models, and 4) an archival software and data repository. Results: Our IV&V plan was carried out with assistance from both the funding agency and participating teams of researchers. QA/QC of data acquisition aided in process improvement and the flagging of experimental errors. Holdout validation set tests provided an independent gauge of model performance. Conclusions: In circumstances that do not support a fully open approach to scientific criticism, standing up independent teams to cross-check and validate the results generated by primary investigators can be an important tool to promote reproducibility of results.

59 BASIC BIOLOGICAL SCIENCES↗

Passive method to measure strength of turbulence

Disclosed is a method to passively measure and calculate the strength of turbulence via the index of refraction structure constant Cn2 from video imagery gathered by an imaging device, such as a video camera. Processing may occur with any type computing device utilizing a processor executing machine executable code stored on memory. This method significantly simplifies instrumentation requirements, reduces cost, and provides rapid data output. This method combines an angle of arrival methodology, which provides scale factors, with a new spatial/temporal frequency domain method. As part of the development process, video imagery from high speed cameras was collected and analyzed. The data was decimated to video rates such that statistics could be computed and used to confirm that this passive method accurately characterizes the atmospheric turbulence. Cn2 accuracy from this method compared well with scintillometer data through two full orders of magnitude and more capability is expected beyond this verification.

O'Neill, Mary Morabito↗

FAUST - Benchmark and Validation Framework [Slides]

This presentation discusses Python packages for nuclear data applications and for benchmarking. It also discusses the goals and objectives for FAUST which include providing input and output processing for different calculation codes, allowing for exchanging results between different applications and codes, running benchmarks and processing the results, automating and simplify plot and report generation, and providing a basis for developing applications useful for nuclear data evaluators. The presentation also discusses how the packages can be used both inside and outside LANL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation of an Advanced High-Pressure PEM Electrolyzer and Composite Hydrogen Storage

One of the challenges associated with transitioning our nation’s transportation economy to any alternative fuel is bringing fuel production, storage, transportation, and end-customer delivery infrastructure online. There are challenges associated with planning and implementing hydrogen fueling stations at locations relevant to the parallel effort of bringing the related transportation vehicles into the market. There are additional challenges associated with the operation of hydrogen fueling stations and reducing operating costs through energy reduction. Through calculation and design efforts it was concluded that elevating the pressure of a PEM cell stack and further developing hydrogen processing equipment that support elevated pressures while maintaining similar characteristics to their current production counterparts will yield in energy savings and by extension cost reductions. In this report, the design of a 57-bar (827-psig) electrolysis system and electrolysis cell stack is described. Furthermore, the fabrication, component verification, and integrated testing are summarized. Prior to this development effort, Proton had demonstrated thousand of hours of durability on 30-bar (435-psig) differential pressure PEM electrolysis hardware. Proton has sold and installed many of these 30-bar devices at customer sites. Through this project, Proton has developed and tested 57-bar PEM electrolysis, thereby extending the proven capability of PEM electrolysis as a hydrogen generation and electrochemical compression device. A design and validation effort to arrange hydrogen processioning equipment for refueling OEM vehicles was also undertaken within this project. By arranging equipment within two 6-meter ISO shipping containers and further developing NFPA 2 Hydrogen Technologies Code to support containerized hydrogen equipment, efforts to site, install, and commission, and approve hydrogen fueling stations were reduced and now provide a groundwork for future installations.

08 HYDROGEN↗

Analysis with the ORIGEN Module of PNAR Spent Fuel Measurements for Nuclear Safeguards Applications

This report documents the analysis of the Passive Neutron Albedo Reactivity (PNAR) measurements for 23 boiling water reactor spent fuel assemblies that were performed in Finland under Action Sheet 65, which is an international collaboration on spent fuel safeguards verification methods in the context of the Finnish spent fuel encapsulation/repository system. PNAR measures the passive neutron and gamma emission rates from each of the spent fuel assemblies like a Fork detector, and it also measures the PNAR ratio, which is expected to correlate with the net neutron multiplication of the measured fuel assembly. The analysis was performed with the Oak Ridge Isotope Generation and Depletion (ORIGEN) Data Analysis Module, which was originally developed for predicting Fork detector spent fuel measurement signals in real time and has been integrated into the Integrated Review and Analysis Package developed by Euratom and the International Atomic Energy Agency. The module includes the ORIGEN burnup analysis code and integrates detector response functions pregenerated using the MCNP Monte Carlo transport code to predict the detector signals in several seconds per assembly. In this study, new response functions specific to the analyzed PNAR measurements were generated for the ORIGEN Module. The study also analyzes the impact of using detailed fuel design and operation information vs. standard safeguards information on the results calculated with the ORIGEN Module. Using detailed information reduced the standard deviation of the relative differences between calculated and measured neutron count rates among the 23 assemblies from ~10% to ~4%. The results obtained using standard safeguards information for these PNAR measurements were similar to those obtained for the Fork detector. A clear trend was found between the calculated net neutron multiplications and the measured PNAR ratios of the 23 assemblies. An uncertainty assessment of the calculations was also performed to estimate the potential impact of the accuracy and the completeness of declaration information on the calculated results for the count rates and net neutron multiplications. The ORIGEN Module predictions for PNAR signals and net neutron multiplication are expected to directly support the safeguards inspectors’ efforts to verify operator declarations of a spent fuel assembly in real time in safeguards practices.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Improvements of SAM Heat Transfer Models for Molten Salt-Cooled Pebble Bed Reactors

To support the development and utilization of the SAM code for molten salt-cooled pebble bed reactor safety analyses and licensing, selected improvements were made to the one-dimensional system-level heat transfer modeling in SAM with the goal of better capturing heat transfer phenomena specific to these designs. In consideration of molten salt as a coolant, new models were introduced to capture the enhanced radiative heat transfer that must be considered at the higher working temperatures. Radiative heat transfer is handled not only between two structural surfaces, but also in combination between the surfaces and the molten salt fluid, which has a higher opacity and absorption of radiative heat as compared to water. Several examples of analytical verification of the new radiative heat transfer methods are included in this report. Additionally, convective heat transfer correlations were selected from literature and included in SAM to account for the enhanced heat transfer expected in pebble bed geometry flows. These correlations are verified and demonstrated in the present work with a reference FHR example.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Draft Rules for Alloy 617 Creep-Fatigue Design Using an EPP+SMT Approach

The report describes an improved new creep-fatigue design method for structural components in high temperature nuclear service. The new method uses an integrated elastic-perfectly plastic (EPP) analysis and Simplified Model Test (SMT) approach that greatly simplifies the design process by avoiding the separate evaluation of creep and fatigue damages and eliminating the requirement of stress classification. The report includes draft design rules, presented in a format compatible with an ASME nuclear Code Case, as well as a commentary on the rules and the validation of the rules using pressurized SMT (p-SMT) test data. Areas of potential improvement have been identified for further evaluation. The report also develops EPP+SMT design charts for Alloy 617 at temperatures between 800°C and 950°C by extrapolating high strain range, short hold test data to low strain range and long hold times, typical for structural components in high temperature nuclear service. The design charts are currently limited to 10,000 cycles to avoid excessive extrapolation outside the test database. More robust alternate extrapolation procedures will be considered in future work. The report includes several verification problems comparing the new EPP+SMT design method with current ASME creep-fatigue design methods. The comparison demonstrates that the EPP+SMT method significantly reduces the over-conservatism in current methods. Finally, the report includes a fully-documented sample problem with multiple load cases detailing the application of proposed EPP+SMT design rules.

36 MATERIALS SCIENCE↗

Analysis of PNAR Spent Fuel Safeguards Measurements using the ORIGEN Data Analysis Module

This paper summarizes the analysis of the Passive Neutron Albedo Reactivity (PNAR) measurements using the ORIGEN data analysis module for 23 boiling water reactor spent fuel assemblies that were performed in Finland under an international collaboration on spent fuel safeguards verification methods. PNAR is part of the proposed integrated nondestructive system to be used for safeguards verifications at the planned Finnish encapsulation facility. Besides measuring passive neutron and gamma emission rates from an assembly like a Fork detector, PNAR also measures the PNAR ratio, which is expected to correlate with the fissile content in the assembly. The emission rates and PNAR ratio can be used to verify the operator declarations and the fissile content of an assembly, respectively. The analysis was performed using the ORIGEN Data Analysis Module, which was originally developed for predicting Fork detector neutron and gamma signals for spent fuel measurements and has been integrated into the Integrated Review and Analysis Package developed by Euratom and the IAEA. The Module includes the ORIGEN burnup analysis code and integrates detector response functions pre-generated using MCNP to predict detector signals in several seconds per assembly. In this study, new response functions specific to PNAR measurements were generated for ORIGEN Module. The study also analyzes impacts of using detailed fuel design and operation information vs. standard safeguards information on results calculated by ORIGEN Module. Using detailed information reduced the standard deviation of relative differences between calculated and measured neutron count rates among the 23 assemblies from ~10% to ~4%. The results obtained using standard safeguards information for these PNAR measurements were similar to those obtained for the Fork detector. A clear trend was found between the calculated net neutron multiplications and the measured PNAR ratios of the 23 assemblies. This paper describes how ORIGEN Module calculates the expected PNAR neutron and gamma signals and PNAR ratio and how they compare with corresponding measured values.

Ilas, Germina↗

NREL Hydrogen Sensor Testing Laboratory

Relevance: Detection is recognized as a critical element for hydrogen facility safety design and supports risk mitigation. Detection methodologies will support validation of H2 behavior research. Hydrogen point sensors play a critical role for safety and process monitoring, but other methodologies can be developed. Approach: NREL Sensor Laboratory tests and verifies sensor performance for manufacturers, developers, end-users, regulatory agencies and SDOs/CDOs NREL deployment activity supports regulatory requirement verification, hydrogen behavior models, and method development for use by stakeholders. Accomplishments and Progress: NREL's R&D accomplishments have supported developers, industry, and SDOs by providing sensor performance and deployment expertise not otherwise available. Development of alternative detection strategies for hydrogen applications have been initiated. HyWAM and advanced detection methodologies deployments at H2@Scale Facilities are being implemented. Collaborations: Collaboration with government laboratories, universities, private organizations and regulatory agencies has leveraged the NREL Sensor Laboratory's success in advancing hydrogen safety sensors and process control. Proposed Future Work: NREL will support hydrogen deployment by the proper implementation of hydrogen sensors and advanced detection strategies. NREL will continue to support science-based codes and standards. This effort will be guided by the needs of the hydrogen community.

codes and standards↗

Develop and verify soil/structure interaction for pile/foundation interaction

Phase II of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project was used to verify the implementation of a new soil-structure interaction (SSI) model for use within offshore wind turbine modeling software. The REDWIN Macro-element model implemented and verified in this study enables a computationally efficient way to model the linear and nonlinear SSI problem, including hysteretic damping, of a monopile structure. The modeling approach was integrated into several modeling tools and a series of increasingly complex simulations was conducted using the IEA 10MW reference turbine mounted on a monopile support structure to verify the coupling between the tools and the REDWIN Macro-element SSI model. This campaign includes only numerical verification between various software and modeling approaches so no experimental measurements are available. The load cases (LC) considered include: LC1 – static response of the tower and substructure LC2 – frequency and mode-shape analysis of the tower and substructure LC3 – response of the tower and substructure due to wind-only loading LC4 – response of the tower and substructure due to wave-only loading LC5 – response of the tower and substructure due to wind and wave loading. Detailed properties of the modeled system are found in the following reference, “Bergua, Roger, Amy Robertson, Jason Jonkman, and Andy Platt. 2021. "Specification Document for OC6 Phase II: Verification of an Advanced Soil-Structure Interaction Model for Offshore Wind Turbines.” Golden, CO: National Renewable Energy Laboratory. NREL/TP-5000-79938. https://www.nlr.gov/docs/fy21osti/79938.pdf. Details on the results from the OC6 Phase II project can be found in the following reference, “Bergua R, Robertson A, Jonkman J, et al. OC6 Phase II: Integration and verification of a new soil–structure interaction model for offshore wind design.” Wind Energy. 2022;25(5):793-810. doi:10.1002/we.2698

17 WIND ENERGY↗

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗