Search NASA⌕ Search

SEARCH · Search NASA

Results for “test functions”

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

NREL's Advanced Distribution Management System (ADMS) Test Bed

Transform utility electric distribution management systems to enable the integration and management of all assets and functions across the utility enterprise regardless of vendor or technology. Four program areas: Platform: Develop an open-source platform; evaluate advanced applications. Test bed: Build a vendor-neutral test bed to evaluate existing and future advanced distribution management system (ADMS) functionalities in a realistic setting. Applications: Develop an initial suite of ADMS applications. Advanced control: Develop new integrated optimization and control solutions.

ADMS test bed↗
DOE OSTI2023

A publicly available PyTorch-$\mathrm{ABAQUS}$ $\mathrm{UMAT}$ deep-learning framework for level-set plasticity

Here this paper introduces a publicly available PyTorch-ABAQUS deep-learning framework of a family of plasticity models where the yield surface is implicitly represented by a scalar-valued function. In particular, our focus is to introduce a practical framework that can be deployed for engineering analysis that employs a user-defined material subroutine (UMAT/VUMAT) for ABAQUS, which is written in FORTRAN. To accomplish this task while leveraging the back-propagation learning algorithm to speed up the neural-network training, we introduce an interface code where the weights and biases of the trained neural networks obtained via the PyTorch library can be automatically converted into a generic FORTRAN code that can be a part of the UMAT/VUMAT algorithm. To enable third-party validation, we purposely make all the data sets, source code used to train the neural-network-based constitutive models, and the trained models available in a public repository. Furthermore, the practicality of the workflow is then further tested on a dataset for anisotropic yield function to showcase the extensibility of the proposed framework. A number of representative numerical experiments are used to examine the accuracy, robustness and reproducibility of the results generated by the neural network models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

In-Sodium Testing of a Prototype Thermoacoustic Power Sensor for Sodium-Cooled Fast Reactors

The ultimate goal of this project is to develop and demonstrate a thermoacoustic power sensor (TAPS) for Sodium-Cooled Fast Reactors (SFRs), with potential application also envisioned to other nuclear technologies such as Lead-Cooled Fast Reactors (LFRs), Molten Salt Reactors (MSRs), in addition to Light Water Reactors (LWRs). The project was led by Westinghouse Electric Company, LLC (Westinghouse) and carried out in collaboration with Argonne National Laboratory (ANL) and the University of Pittsburgh. A TAPS is a passive (self-powered), non-invasive (wireless) sensor envisioned for measuring key parameters, such as local temperature and neutron flux, in a nuclear reactor core. The sensor generates pressure waves (i.e., sound waves) with a frequency and amplitude dependent upon nuclear operating conditions (coolant temperature or power changes). The acoustic waves are able to travel through the core and associated structures, and reach to the sensor network placed outside and/or inside of the reactor vessel. These sensors require a very small amount of power which, during loss of power events, can be provided, for example, by harvesting gamma radiation energy, thus resulting in a monitoring system that can function both during normal operation and during loss of power events. Westinghouse and the University of Pittsburgh designed and fabricated TAPS prototypes for Argonne National Laboratory (ANL) to carry out in-sodium testing to evaluate the effects of sodium on the TAPS and the performance of the TAPS technique in sodium. Argonne received a TAPS prototype from Westinghouse, and the prototype was modified such that it can be installed into a test vessel and function in sodium at elevated temperature without potentially leaking. A water mockup test apparatus was constructed to validate proper working of the prototype. An instrumentation and control (I&C) system, running on the National Instruments (NI) LabView platform, was developed to: 1) operate both the water mockup test and the in-sodium test facility; and 2) process and analyze the received acoustic signals from an array of accelerometers and the Argonne sodium-submersible high-temperature acoustic sensor. The prototype was successfully tested in a water bath at different temperatures. Water mockup tests demonstrated that the TAPS prototype is working properly and its resonance frequency changes linearly with the coolant (water) temperature. A TAPS test apparatus was constructed and integrated with the upgraded Under-Sodium Viewing (USV) sodium test facility. The integrated USV-TAPS sodium test facility has been operational. The TAPS prototype and a high-temperature sodium-submersible acoustic sensor (SSAS) developed by Argonne were both installed inside the TAPS test vessel. Being operated within argon cover gas under ambient conditions, the TAPS prototype demonstrated that it was functioning properly with a resonance frequency at 1407.2 Hz, which was successfully detected by the accelerometers mounted on the external wall of the vessel and the high-temperature SSAS installed inside the vessel. After successfully transferring sodium into the vessel, in-sodium tests of the prototype were conducted. Tests of the TAPS prototype demonstrated that the resonance frequency of the TAPS changes linearly with respect to the temperature difference between the interior of the TAPS and bulk sodium. The early tests showed that the TAPS prototype could not establish a continuous and consistent resonance in sodium. The resonance diminished before the prototype reached its operating temperature. A signal postprocessor was added to the DAQ module latterly to isolate interferences, enhance signal conditioning, improve peak detection, and generate resonance frequency versus temperature plots. After testing in molten sodium and immersion at higher temperature for several weeks, the TAPS prototype was able to establish a continuous resonance. Performance evaluation of the TAPS prototype was then conducted in sodium. The tests included the investigation of 1) the effects of the temperature difference between the TAPS and bulk sodium, 2) the effects of sodium flowrate; and 3) the performance of the different sensor-receiver systems positioned inside or outside the vessel. Results of a test demonstrated that, with limited sodium circulation, a continuous and consistent resonance of the TAPS prototype was established occasionally. The tests also demonstrated that, because of the nature of detection principles and mounting methods, the high temperature SSAS is more affected by acoustic noise, while accelerometers are more affected by vibrations, in the test environment. It is unknown why the TAPS prototype only occasionally established a continuous and consistent resonance when the TAPS temperature reached its operating temperature in molten sodium, and why it ultimately failed to resonate at all. A failure modes assessment was conducted and a few potential causes of failure were identified. Different post in-sodium tests were conducted to obtain additional information potentially relevant to the cause of the failure. Nondestructive evaluation techniques are suggested to examine the internal integrity as well as the gas mixture of the prototype. If they prove inconclusive, the prototype should be cut open to conduct a thorough inspection of its internal integrity and determine the state of the gas mixture.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗
DOE OSTI2025

Mortality correlates with tree functional traits across a wood density gradient in the Central Amazon

Introduction: Understanding the mechanisms of tree mortality in tropical ecosystems remains challenging, in part due to the high diversity of tree species and the inherently stochastic nature of mortality. Plant functional traits offer a mechanistic link between plant physiology and performance, yet their ability to predict growth and mortality remains poorly understood. Given recent increases in tree mortality rates in the Amazon forest following extreme drought and wind events, we tested if lower wood density and acquisitive plant functional traits were associated with increased growth and mortality for common co-occurring trees in the Central Amazon. Methods: Seventeen trees of different species with similar sizes but a range in wood density (WD) and wood traits were felled, then assessed for 27 different individual functional parameters, including whole tree architecture, stem xylem anatomical and hydraulic traits and leaf traits. Traits of the individual trees were related to stand-level growth and mortality rates collected periodically over 30 years from nearby permanent inventory plots. Results: Higher wood density was associated with smaller leaf size, lower foliar base cations, lower stem water content and sapwood fraction, in agreement with the fast-slow plant economics spectrum. Lower wood density was associated with more acquisitive characteristics with greater hydraulic capacity and foliar nutrient concentrations, correlating with greater growth and mortality rates. Discussion: Our results show that lower wood density is part of a coordinated suite of traits linked to high resource acquisition, fast growth, and increased mortality risk, providing a functional framework for predicting species performance and forest vulnerability under future climate stress.

demographics↗
DOE OSTI2026

Port scanner and Testing Suite

This project addresses the challenge of identifying and managing open network ports across physical and virtual hosts. The current form of verifying ports in use required manually searching individual ports - a process that was both time- consuming and a potential bottleneck for deployment timelines. To resolve this, an automated port scanning tool was developed in Python. The tool supports simultaneous multiple port scans. To ensure functionality and long-term maintainability, a comprehensive testing suite was implemented using Python’s unittest framework. Edge cases, including valid port numbers, reversed ranges, and closed ports, were explicitly tested to ensure robust handling of real-world scenarios. The resulting tool reduces the time required to verify port security across a network, supporting both targeted and host checks and broader Classless Inter-Domain Routing (CIDR) -based network scans. This work demonstrates the value of automation and test-driven development in strengthening network security practices, and provides a foundation for future enhancements.

Rivera, Linda [Fermilab]↗
DOE OSTI2020

Development and Demonstration of a Class 6 Range-Extended Electric Vehicle for Commercial Pickup and Delivery Operation: Preprint

Range-extended hybrids are an attractive option for medium- and heavy-duty commercial vehicle fleets because they offer the efficiency of an electrified powertrain with the driving range of a conventional diesel powertrain. The vehicle essentially operates as if it was purely electric for most trips, while ensuring that all commercial routes can be completed in any weather conditions or geographic terrain. Fuel use and point-source emissions can be significantly reduced, and in some cases eliminated, as many shorter routes can be fully electrified with this architecture. Under a U.S. Department of Energy (DOE)-funded project for Medium- and Heavy-Duty Vehicle Powertrain Electrification, Cummins has developed a plug-in hybrid electric Class 6 truck with a range-extending engine designed for pickup and delivery application. The National Renewable Energy Laboratory (NREL) assisted by developing representative workday drive cycles for Class 6 operation and an adapted cycle to enable vehicle track testing. A novel, automated driving system was utilized by Southwest Research Institute (SwRI) to improve the repeatability of the track testing conducted to quantify vehicle energy consumption. Cummins used a drivetrain model to design the hybrid control system for increased fuel savings. The control system functionality and fuel savings objective of 50% or more were confirmed by track testing described in this paper.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Pathways to IEC Standards for Heliostat Design Qualification and Site Acceptance in Central Receiver CSP Applications

This paper surveys the existing landscape of standards relevant to heliostats, identifies their gaps, and proposes a path forward to a comprehensive set of heliostat guidelines, technical specifications, and standards under the framework of International Electrotechnical Commission (IEC) TC 117. Gaps in existing guidelines and standards are surveyed using a three-tiered taxonomy: component-level, heliostat-level, and field-level. At each level, the gap analysis is followed by a proposal for a coordinated path forward on the development of standards. At the component level, advances in the understanding of wind loading should inform a technical specification for drives and structures. Reflectors require consolidation of measurement guidelines into existing standards documents. Communications & controls require technical standards to inform their selection and secure implementation. At the heliostat level, IEC 62817 (solar trackers) adequately characterizes drive systems, structures, and electronics, but requires adaptation to heliostats' use patterns, operating modes, and expected life cycles. IEC 62817 does not address heliostat beam quality and pointing accuracy, but the process for determining both is elaborated in the SolarPACES Guideline for Heliostat Performance Testing. This SolarPACES document requires two main modifications: adaptation to IEC language and inclusion of testing after heliostats which have undergone accelerated weathering and mechanical cycling (to understand performance degradation). At the field level, IEC 62862-4-2 addresses the function and control of heliostat fields but does not cover the statistically rigorous testing of heliostat groups, or field performance factors like security and soiling. The addition of documents under IEC-62862-4 is proposed to address this gap.

beam↗
DOE OSTI2021

Measurement of internal dark current in a 17 GHz accelerator structure with an elliptical sidewall

A 17 GHz single cell, standing wave, copper accelerator structure with an axisymmetric elliptical central cell sidewall was tested for internal and downstream dark current as a function of gradient up to 93 MV/m. The elliptical sidewall was predicted to suppress the internal dark current and the lower order multipactor modes as compared with a previously tested structure having a straight sidewall. During the conditioning phase of the elliptical sidewall structure, strong internal dark current generated by an N = 1 multipactor mode was observed at gradients in the 10 to 20 MV/m range. After conditioning with 2.2x10^5 pulses to 93 MV/m, the N = 1 mode was completely suppressed and no multipactor resonances were observed. The internal dark current was reduced to a comparatively low level, much smaller than in the previously studied, straight sidewall structure, in good agreement with simulations. The energy spectrum of the electrons colliding with the sidewall was measured using an isolated side dark current monitor and a bias voltage. As the conditioning progressed, the electron energy spectrum showed an increase in the concentration of lower energy electrons, also in good agreement with simulations. Studies of internal dark current may help to understand the rf conditioning and ultimate performance of high gradient accelerator structures.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Plutonium-238 Loadout Calorimeter Return To Service

The calorimeter in the transfer area in Building 7920 at the US Department of Energy’s Oak Ridge National Laboratory is used to measure the amount of 238 Pu in the Plutonium Heat Source Oxide. The calorimetry measurement is taken after the oxide is packaged in an ORISC (Oak Ridge Inner Shipping Capsule)/FSO (Fuel Storage Outer) assembly before loading into the Containment Vessel and shipping container. The equipment is used in three specific configurations: zero measurement (nothing in the sample cup), FSO measurement (actual measurement of 238 Pu), and reference calibration measurement (using the reference FSO as a known heat source). System checks were performed in July 2023 in preparation for scheduled material loadout. During these checks, personnel discovered that the reference FSO had experienced an electrical short in the heater cable and was no longer functional. The calorimeter was in the reference calibration measurement configuration at this time. Facility management was notified when the damage was discovered. The data from the Yokogawa recorder were pulled and analyzed for a span of time between September 2022 and July 2023. An unusual spike was found in the sample cup temperature data that started on October 19, 2022, and ended on October 20, 2022. The sample cup reached a temperature of 191.1°C. The maximum allowable temperature for the calorimeter thermopile to reach is 220°C based on the manufacturer statement provided in Appendix B of this report [1]. The calorimeter thermopile and the sample cup thermocouple are physically close to each other; therefore, the data suggest that the calorimeter itself never saw a temperature in excess of 191.1°C. Physical damage present on the cable, tape, and insulation suggests that the temperature reached a higher value in some locations within the sample cup than the data shows. Because of this discrepancy, testing of the calorimeter was required to confirm that the calorimeter is functioning the same as it was before the electrical short. The intent of this report is to assess the functionality of the calorimeter and determine if it is acceptable to use for 238 Pu loadout material measurement. Testing showed that the calorimeter is functionally equivalent to the initial calibration that was performed in the summer of 2021, and thus, the calorimeter was cleared to be used for its intended purposes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Test–retest reliability for a social discounting of personal information task

Increasing cybercrime rates means identifying potential victims is critically important. Social discounting tasks show that individuals share less personally identifying information as social distance increases. However, the test–retest reliability and uniqueness of this measure is unclear. The current study assessed social discounting for personally identifying information (SDPII), delay discounting, risk taking, and personality at two measurement waves 30 days apart for 64 undergraduate students. Test–retest reliability was statistically significant for the SDPII and all other measures, replicating previous studies. SDPII rates were not significantly correlated with other measures during both measurement waves, showing discriminant validity. SDPII rates were lower than those reported in a previous study but were still well described by a hyperbolic discounting function, suggesting replicability across studies. Furthermore, the high test–retest reliability, uniqueness, and replicability of the SDPII suggests that it may quantitatively identify cybercrime victimization. Future research should test which measure or combination of measures can accurately predict scam and cybercrime victimization to inform data-based interventions.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Feasibility study for test rig assessments of fish passage conditions in a Kaplan turbine

The assessment of fish passage conditions in hydroelectric turbines consists of identifying and quantifying physical magnitudes leading to increased risks of injury of fish passing through turbines in operation. Such assessments are usually carried out either with the use of computer-based methods during design or with field testing of live fish and sensors passing through prototypes. A method in between consists of test rig experimentation, which is critical for testing fish-focused design concepts and offers the opportunity for implementing the most effective design measures for improved fish survivability. However, fish-related assessments in test rigs are not sufficiently documented for industrial applications. This work presents the main findings of an experimental campaign to quantify fish-related hydraulic magnitudes in a physical model of a Kaplan turbine in a commercial test rig. Two operating conditions were tested by releasing miniaturized autonomous sensor devices (termed Sensor Fish Mini) at the turbine intake flow, passing them through the runner in motion and recovering them at the draft tube exit. During passage, time series of acceleration, absolute pressure and rotational velocity were recorded. The recordings were then interpreted to determine the magnitude and likely location of hydraulic stressors hazardous to fish. The statistical tests on the reported measurements indicated that low pressure, collisions on the runner and rotations in the draft tube were not different between the two tested operating points. On the other hand, pressure drop and collision rates on the distributor differed considerably as a function of net head. The outcomes of this investigation showed that test rig evaluations of fish-related properties with Sensor Fish Mini can contribute to an evidence-based development of turbine geometries designed for providing safer passage conditions. Future work will investigate the scaling of test rig measurements to hydraulically equivalent magnitudes in the prototype and their biological consequences.

13 HYDRO ENERGY↗
DOE OSTI2026

Secure Time Synchronization in Power Grids and Network HIL Synchrophasor Testing

Reliable and secure time synchronization underpins the monitoring and control functions of modern power grids. As GPS-based timing infrastructures grow more susceptible to spoofing and jamming, their vulnerabilities pose escalating risks to grid stability. This work investigates a secure, resilient timing framework that can serve as an alternative or redundant source for power grids, with a particular focus on synchrophasorbased applications. A candidate timing system architecture is evaluated to guarantee trustworthy time dissemination, even in degraded conditions. A network hardware-in-the-loop testing of two synchrophasors validates the concepts, demonstrating enhanced timing integrity, improved detection of timing anomalies, and sustained observability during adverse timing events.

Wu, Ori [ORNL] (ORCID:0000000326723410)↗
DOE OSTI2021

Extending a 10‐node composite tetrahedral finite element for solid mechanics

Abstract We propose to extend the composite tetrahedral finite element first introduced by Thoutireddy et al. (2002) and recently reformulated by Ostien et al. (2016). We generalize the gradient operator and mass matrix to curved domains through analytical expressions weighted by subtetrahedra Jacobians. Optimal integration weights are constructed to increase the accuracy of Gaussian quadrature. We preserve the variational structure of the formulation through a new five‐field functional with additional, independent fields for the Jacobian and the pressure. This approach not only obviates volumetric locking but also yields symmetry. A deleterious soft mode, common to both the quadratic and composite tetrahedral element with constant pressure formulations is effectively stabilized through a novel convex energy penalty function. Numerous numerical examples spanning a patch test to the impact of a Taylor bar demonstrate the accuracy, robustness, and convergence of the extended composite tetrahedral element for application to structural metals.

Foulk III, James W.↗
DOE OSTI2026

Constraints on long-range forces in de Sitter space

The representation theory of de Sitter space admits partially massless (PM) particles, but whether such particles can participate in consistent interacting theories remains unclear. We investigate the consistency of theories containing PM fields, particularly when these fields are coupled to gravity. Our strategy exploits the fact that PM fields correspond to partially conserved currents on the spacetime boundary, which generate symmetries. These symmetries place stringent constraints on correlation functions of charged operators, allowing us to test the consistency of a proposed bulk spectrum. When the assumed operator content violates these constraints, the corresponding bulk theory is ruled out. Applying this framework, we show that, in four-dimensional de Sitter space, PM fields of spin 2 or 3 (at depth 0) cannot couple consistently to gravity: such couplings necessitate additional massive fields, which are inevitably non-unitary. In higher dimensions, however, the constraints can be satisfied without violating unitarity if further PM fields are included. The resulting structure leads to additional charge conservation laws, which suggests that consistency may ultimately require an infinite tower of higher-spin PM fields, akin to the situation for ordinary higher-spin symmetries. The methods developed here provide powerful constraints on possible long-range interactions in de Sitter space and delineate the landscape of consistent quantum field theories in cosmological spacetimes.

AdS-CFT Correspondence↗
DOE OSTI2023

Corrosion property of Alloy 625 in Molten FLiNaK salt according to the Tellurium Concentrations

For this work, the corrosion properties of Alloy 625 have been studied in molten FLiNaK salt for 48 h at 700 °C as a function of tellurium (Te) concentrations using immersion tests and electrochemical methods. Under specific Te concentrations, we found that Te-induced corrosion can be suppressed and help protect the alloy against corrosion. The sample after the immersion test containing 0.1 wt.% Te (0.1 wt.% Te sample) had the lowest corrosion depth and mass loss, coupled with the highest charge transfer resistance, obtained using electrochemical impedance spectroscopy (EIS). On the other hand, the alloy with Te content above 0.1 wt.% exhibited severe corrosion penetration and lower charge transfer resistance in FLiNaK. Based on thermodynamical and electrochemical analyses, a corrosion mechanism of Alloy 625 in the molten salt containing Te was theorized involving the formation of stable corrosion products such as Ni/Fe rich layer depleted in Cr, nickel telluride, and chromium oxide.

Hong, Minsung↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗