SEARCH · Search NASA

Results for “collocation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Use of TEMPO as a Proxy for Hyperspectral Geostationary Ocean Color Measurements from the GeoXO OCX Instrument: Harnessing Machine Learning and Principal Component Techniques for Atmospheric and Glint Correction

Retrievals of ocean color from space are important for better understanding the ocean ecosystem. The launch of atmospheric geostationary hyperspectral sensors such as TEMPO, provides a unique opportunity to examine the diurnal variability in ocean ecology. While TEMPO does not have as high spatial resolution or full spectral coverage as planned coastal ocean sensors such as the Geosynchronous Littoral Imaging and Monitoring Radiometer (GLIMR) or GeoXO Ocean Color instrument (OCX), its hourly measurements provide coverage of regions such as Lake Erie and the Gulf of Mexico at spatial scales of approximately 5 km. These data can be useful for testing new algorithms. We will apply our newly developed machine learning based atmospheric correction approach for ocean color retrievals to TEMPO data. Our approach begins by decomposing measured radiances from hyperspectral sensors into spectral features that describe the scattering and absorption of the atmosphere as well as the underlying surface reflectance. The coefficients of the principal components are then used to train a neural network to predict ocean color properties derived from collocated MODIS/VIIRS physically-based retrievals. This machine learning approach does not rely on radiative transfer modeling, and the use of MODIS/VIIRS data for training accounts for possible calibration b in hyperspectral data. Previously, we applied our approach using blue and UV wavelengths with the Ozone Monitoring Instrument (OMI) and TROPOspheric Monitoring Instrument (TROPOMI) to show that it can estimate ocean color properties in less-than-ideal conditions such as lightly to moderately clouded conditions as well as sun glint and thus improve the spatial coverage of ocean color measurements. TEMPO provides an opportunity to improve on this approach since it will provide collocated measurements at green and red wavelengths that were not available from OMI and TROPOMI and are important particularly for coastal waters. Additionally, our technique can be applied early in the mission and has potential to demonstrate the value of near real time ocean color products that are important for monitoring of harmful algae blooms and other oceanic phenomena.

Zachary Fasnacht↗

NASA NTRSDate not supplied

Creating Regional and Seasonal Climatologies of Marine and Dusty Marine Aerosol Lidar Ratios using MODIS AOD Constrained Retrievals and GOCART Model Simulations

The CALIPSO aerosol algorithms currently assign one lidar ratio (LR) value globally for each of the seven tropospheric aerosol types. In this study, MODIS total column aerosol optical depths (AODs) are used to constrain collocated CALIOP backscatter profiles in a Fernald inversion that infers aerosol LRs for CALIOP-classified marine and dusty marine aerosols. The GOCART aerosol model is leveraged to estimate the sea salt volume fractions (SSVFs) that are collocated with the CALIOP+MODIS LR retrievals. An inverse empirical relationship is found between the SSVFs and LRs (i.e., smaller SSVFs and larger LRs near coastlines, but the opposite in the remote oceans). This SSVF/LR relationship is applied to create regional and seasonal hybrid (i.e., retrieval & model-assisted) climatological LR maps so as to develop more robust LR selections for marine and dusty marine aerosols in the CALIPSO algorithms. These analyses also provide critical LR information for the next generation of spaceborne elastic backscatter lidars.

Travis D. Toth↗

NASA NTRSDate not supplied

Creating Regional and Seasonal Climatologies of Marine and Dusty Marine Aerosol Lidar Ratios using MODIS AOD Constrained Retrievals and GOCART Model Simulations

Travis D. Toth↗

NASA NTRS2024

Retrievals of Aerosol Optical Depth Over the Western North Atlantic Ocean During ACTIVATE

Aerosol optical depth was retrieved from two airborne remote sensing instruments, the Research Scanning Polarimeter (RSP) and Second Generation High Spectral Resolution Lidar (HSRL-2), during the National Aeronautics and Space Administration (NASA) Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE). The field campaign offers a unique opportunity to evaluate an extensive 3-year dataset under a wide range of meteorological conditions from two instruments on the same platform. However, a long-standing issue in atmospheric field studies is that there is a lack of reference datasets for properly validating field measurements and estimating their uncertainties. Here we address this issue by using the triple collocation method, in which a third collocated satellite dataset from the Moderate Resolution Imaging Spectroradiometer (MODIS) is introduced for comparison. HSRL-2 is found to provide a more accurate retrieval than RSP over the study region. The error standard deviation of HSRL-2 with respect to the ground truth is 0.027. Moreover, this approach enables us to develop a simple, yet efficient, quality control criterion for RSP data. The physical reasons for the differences in two retrievals are determined to be cloud contamination, aerosols near the surface, multiple aerosol layers, absorbing aerosols, non-spherical aerosols, and simplified retrieval assumptions. These results demonstrate the pathway for optimal aerosol retrievals by combining information from both lidars and polarimeters for future airborne and satellite missions.

Aerosol optical depth↗

NASA NTRSDate not supplied

Geodetic Reference Instrument Transponder for Small Satellites

The local ties between Satellite Laser Ranging stations and other collocated geodetic stations are an important component of the realization of the International Terrestrial Reference Frame but the survey accuracies are limited due to the inaccessibility of the instrument reference points. The NASA Geodetic Reference Instrument Transponder for Small Satellites (GRITSS) technology demonstration mission is designed to overcome this limitation by using a small satellite as a space-based reference point for tying together the measurement points of collocated geodetic stations. The GRITSS instrument upconverts the Global Positioning System (GPS) signals received at the satellite in real-time to S and X band and transmits them to a VLBI Global Observing System (VGOS) antenna ground station. The measurement observables do not require the satellite to be in view of more than one VLBI station at a time. The GRITSS demonstration mission utilizes a 12-UL CubeSat in a Low Earth, Sun Synchronous orbit. Tracking by the global ILRS network throughout the life of the mission is critical to establishing the space-based reference along with the precision orbit determination required by the measurement concept. Due to limitations in spacecraft power and other constraints, the initial phase of the mission will only broadcast the GRITSS VLBI signals twice a day to each of the NASA VGOS stations in Maryland, Hawaii, and Texas, but other international VGOS stations may be added during the second half of the mission. An overview of the GRITSS mission and measurement concept will be presented.

Stephen M Merkowitz↗

NASA NTRSDate not supplied

Variable Lidar Ratios for Marine and Dusty Marine Using MODIS AOD Constrained Retrievals and GOCART Model Simulations and their Impact on Level 2 CALIPSO Data Products in Version 5

The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard CALIPSO provided global measurements of attenuated backscatter profiles of various tropospheric aerosol species from June 2006 through June 2023. Extinction profiles are retrieved from these backscatter measurements by assuming a value for the extinction-to-backscatter or lidar ratio (LR). Up until version V4.51 of the data products, only a constant global value for each of the species has been used, spatially and temporally. However, it has been well known from various ground-based measurements that LRs of each aerosol species can potentially vary spatially as well seasonally. For the forthcoming Version 5 of the CALIPSO products, the CALIPSO team will implement variable lidar ratios at 532 nm for the species classified as “marine” and “dusty marine”. Here we present the essential elements of this scheme, which involves estimating lidar ratios by constraining Fernald solutions for the CALIOP particulate backscatter profiles using aerosol optical depths (AOD) from collocated MODIS retrievals. To account for the lack of measurements in various regions and seasons (e.g., Arctic winters), we leverage collocated model sea-salt volume fractions (SSVF) simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model and use an empirical relationship between the SSVF and LR to build hybrid climatological maps for the entire globe, i.e., using the modelled values when reliable values of the constrained LRs are not available. Initial evaluation of the resulting extinction profiles and corresponding AODs, both globally and in specific regions are presented using a limited amount of data. The primary improvement in AOD occurs in the coastal areas, like the Arabian Sea and Bay of Bengal. The estimated marine LR values in these regions are significantly higher than the constant value of 23 sr used in previous data releases, likely due to mixing with offshore pollution which results in higher AODs. Initial validation results using the CALIPSO Ocean Derived Column Optical Depth (ODCOD) retrievals and ground-based AERONET AOD measurements are presented.

CALIPSO↗

DOE OSTI2024

Estimation of General Area Dose Rates During Closure Weld Operations on the DOE Standard Canister Pacific Basin Nuclear Conference 2024

The Department of Energy (DOE) Spent Nuclear Fuel (SNF) Road Ready Demonstration Project is an effort to package DOE-managed SNF for road-ready dry storage. A key operation in this demonstration is performing the closure weld on the DOE Standard Canister after it has been loaded with the SNF. Planning the closure weld operations is ongoing and must consider what parts of the operations (if any) can be done manually and what parts must be done remotely. A key parameter to inform the evolving planning of these operations is the dose rate in the vicinity of the canister and the dose hazard posed to the closure weld operators and collocated workers. This study developed a Monte Carlo N-Particle (MCNP) model of the closure weld operation, which was used to estimate general area dose rates as well as dose rates to the welding location, weld operators, and collocated workers. These estimates provide some insight into the potential radiation exposure hazard for workers during the operations and are a factor to be considered in the ongoing development of designs, technology, processes, and administrative controls for packaging DOE-managed SNF for road-ready dry storage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Accurate data-driven surrogates of dynamical systems for forward propagation of uncertainty

Stochastic collocation (SC) is a well-known non-intrusive method of constructing surrogate models for uncertainty quantification. In dynamical systems, SC is especially suited for full-field uncertainty propagation that characterizes the distributions of the high-dimensional solution fields of a model with stochastic input parameters. However, due to the highly nonlinear nature of the parameter-to-solution map in even the simplest dynamical systems, the constructed SC surrogates are often inaccurate. Here, this work presents an alternative approach, where we apply the SC approximation over the dynamics of the model, rather than the solution. By combining the data-driven sparse identification of nonlinear dynamics framework with SC, we construct dynamics surrogates and integrate them through time to construct the surrogate solutions. We demonstrate that the SC-over-dynamics framework leads to smaller errors, both in terms of the approximated system trajectories as well as the model state distributions, when compared against full-field SC applied to the solutions directly. We present numerical evidence of this improvement using three test problems: a chaotic ordinary differential equation, and two partial differential equations from solid mechanics.

42 ENGINEERING↗

DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗

DOE OSTI2025

Quantifying UAS Observation Error Variance Used in Data Assimilation Systems and Its Impact on Predictive Skill

Observation error determines the weights of the observations and background state used in data assimilation to generate analyses. Quantifying observation error is critical for the optimal assimilation of observational data sets. Uncrewed Aircraft System (UAS) observations have shown potential benefits in filling observational gaps in the lower atmosphere; however, characterization of their error characteristics has been limited. To optimize the use of UAS observations in numerical weather prediction, UAS observation error is estimated based on the 3‐cornered hat diagnostic approach which uses three independent estimates of the atmospheric state. This approach is applied to data from the 2018 Lower Atmospheric Profiling Studies at Elevation‐a Remotely‐piloted Aircraft Team Experiment field campaign using collocated UAS and rawinsonde observations along with output from a set of convection‐permitting model simulations. The estimated observation error values for UAS temperature, wind, and relative humidity measurements were found to be only weakly dependent on height AGL with mean values equal to 0.5°C, 0.8 m s −1 , and 3%, respectively. Only the newly estimated observation error for temperature differed from that previously used to assimilate commercial aircraft observations into global models (1.0°C). However, using this reduced temperature observation error produced more accurate mesoscale analyses and forecasts of both terrain‐driven flows and convection initiation generated by colliding outflow boundaries within the San Luis Valley of Colorado.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

A Model Predictive Control to Improve Grid Resilience

The following article details a model predictive control (MPC) to improve grid resilience when faced with variable generation resources. This topic is of significant interest to utility power systems where distributed intermittent energy sources will increase significantly and be relied on for electric grid ancillary services. Previous work on MPCs has focused on narrowly targeted control applications such as improving electric vehicle (EV) charging infrastructure or reducing the cost of integrating Energy Storage Systems (ESSs) into the grid. In contrast, this article develops a comprehensive treatment of the construction of an MPC tailored to electric grids and then applies it integration of intermittent energy resources. To accomplish this, the following article includes a description of a reduced order model (ROM) of an electric power grid based on a circuit model, an optimization formulation that describes the MPC, a collocation method for solving linear time-dependent differential algebraic equations (DAEs) that result from the ROM, and an overall strategy for iteratively refining the behavior of the MPC. Next, the algorithm is validated using two separate numerical experiments. First, the algorithm is compared to an existing MPC code and the results are verified by a numerically precise simulation. It is shown that this algorithm produces a control comparable to existing algorithms and the behavior of the control carefully respects the bounds specified. Second, the MPC is applied to a small nine bus system that contains a mix of turbine-spinning-machine-based and intermittent generation in order to demonstrate the algorithm’s utility for resource planning and control of intermittent resources. This study demonstrates how the MPC can be tuned to change the behavior of the control, which can then assist with the integration of intermittent resources into the grid. The emphasis throughout the paper is to provide systematic treatment of the topic and produce a novel nonlinear control compatible design framework applicable to electric grids and the control of variable resources. This differs from the more targeted application-based focus in most presentations.

microgrid↗

NASA NTRS2021

Vegetation Monitoring Optimization with Normalized Difference Vegetation Index and Evapotranspiration Using Remote Sensing Measurements and Land Surface Models Over East Africa

The majority of people in East Africa rely on the agro-pastoral system for their livelihood, which is highly vulnerable to droughts and flooding. Agro-pastoral droughts are endemic to the region and are considered the main natural hazard that contributes to food insecurity. Drought begins with rainfall deficit, gradually leading to soil moisture deficit, higher land surface temperature, and finally impacts to vegetation growth. Therefore, monitoring vegetation conditions is essential in understanding the progression of drought, potential effects on food security, and providing early warning information needed for drought mitigation decisions. Because vegetation processes couple the land and atmosphere, monitoring of vegetation conditions requires consideration of both water provision and demand. While there is consensus in using either the Normalized Difference Vegetation Index (NDVI) or evapotranspiration(ET) for vegetation monitoring, a comprehensive assessment optimizing the use of both has not yet been done. Moreover, the evaluation methods for understanding the relationships between NDVI and ET for vegetation monitoring are also limited. Taking these gaps into account we have developed a framework to optimize vegetation monitoring using both NDVI and ET by identifying where they perform the best by using triple collocation and cross-correlation methods. We estimated the random error structure in Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI; ET from the Operational Simplified Surface Energy Balance (SSEBop) model; and ET from land surface models (LSMs). LSM ET and SSEBop ET have been found to be better indicators for vegetation monitoring during extreme drought events, while NDVI could provide better information on vegetation condition during wetter than normal conditions. The random error structures of these variables suggest that LSM ET is most likely to provide important information for vegetation monitoring over low and high ends of the vegetation fraction areas. Over moderate vegetative areas, any of these variables could provide important vegetation information for drought characterization and food security assessments. While this study provides a framework for optimizing vegetation monitoring for drought and food security assessments over East Africa, the framework can be adopted to optimize vegetation monitoring over any other drought and food insecure region of the world.

triple collocation↗

NASA NTRS2024

Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods

To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.

Diagonally-Implicit Runge-Kutta↗

NASA NTRS2024

Deducing Land-Atmosphere Coupling Regimes from SMAP Soil Moisture

In recent years, there has been a growing recognition of the significance of Land-Atmosphere (L-A) interactions and feedback mechanisms in understanding and predicting Earth’s water and energy cycles. Soil moisture plays a critical role in mediating the strength of L-A interactions and is important for understanding the complex and governing processes across this interface. This study aims to identify the significance of soil moisture in identifying L-A coupling strength within the Convective Triggering Potential (CTP) and Humidity Index (HI) framework. To address this, a consistent and reliable dataset of atmospheric profiles is created by merging CTP and HI using Triple Collocation (TC) with three reanalysis datasets. The merged CTP and HI product demonstrates enhanced performance globally as compared to the individual datasets when validated with radiosonde and satellite observations. This merged product of CTP and HI is then used to compare the L-A coupling strength based on Soil Moisture Active Passive Level 3 (SMAPL3) and SMAP Level 4 (SMAPL4) over two decades (2003-2022) where L-A coupling strength is defined as the persistence probability within the dry and wet coupling regimes. Results indicate that the persistency-based coupling strength is related to the ability of soil moisture to predict future atmospheric humidity and dry vs. wet coupling state. The coupling strength in SMAPL4 is consistently stronger than in SMAPL3 and is likely due to its reliance on a land surface model and reduced susceptibility to random noise. The difference in coupling strength based on the same CTP-HI underscores the importance of soil moisture data in estimating coupling strength within the CTP-HI framework. These findings lay the groundwork for understanding the role of L-A interactions and drought evolution due to soil moisture variations, by providing insight into the quantification of coupling strength and its role in drought monitoring and forecast efforts.

Land-atmosphere coupling↗

DOE OSTI2025

An end-to-end deep learning method for solving nonlocal Allen–Cahn and Cahn–Hilliard phase-field models

Here, we propose an efficient end-to-end deep learning method for solving nonlocal Allen–Cahn (AC) and Cahn–Hilliard (CH) phase-field models. One motivation for this effort emanates from the fact that discretized partial differential equation-based AC or CH phase-field models result in diffuse interfaces between phases, with the only recourse for remediation is to severely refine the spatial grids in the vicinity of the true moving sharp interface whose width is determined by a grid-independent parameter that is substantially larger than the local grid size. In this work, we introduce non-mass conserving nonlocal AC or CH phase-field models with regular, logarithmic, or obstacle double-well potentials. Because of non-locality, some of these models feature totally sharp interfaces separating phases. The discretization of such models can lead to a transition between phases whose width is only a single grid cell wide. Another motivation is to use deep learning approaches to ameliorate the otherwise high cost of solving discretized nonlocal phase-field models. To this end, loss functions of the customized neural networks are defined using the residual of the fully discrete approximations of the AC or CH models, which results from applying a Fourier collocation method and a temporal semi-implicit approximation. To address the long-range interactions in the models, we tailor the architecture of the neural network by incorporating a nonlocal kernel as an input channel to the neural network model. We then provide the results of extensive computational experiments to illustrate the accuracy, predictive capabilities, and cost reductions of the proposed method.

42 ENGINEERING↗

DOE OSTI2025

Physics-informed latent neural operator for real-time predictions of time-dependent parametric PDEs

Deep operator network (DeepONet) has shown significant promise as surrogate models for systems governed by partial differential equations (PDEs), enabling accurate mappings between infinite-dimensional function spaces. However, when applied to systems with high-dimensional input-output mappings arising from large numbers of spatial and temporal collocation points, these models often require heavily overparameterized networks, leading to long training times. Latent DeepONet addresses some of these challenges by introducing a two-step approach: first learning a reduced latent space using a separate model, followed by operator learning within this latent space. While efficient, this method is inherently data-driven and lacks mechanisms for incorporating physical laws, limiting its robustness and generalizability in data-scarce settings. Here, in this work, we propose PI-Latent-NO, a physics-informed latent neural operator framework that integrates governing physics directly into the learning process. Our architecture features two coupled DeepONets trained end-to-end: a Latent-DeepONet that learns a low-dimensional representation of the solution, and a Reconstruction-DeepONet that maps this latent representation back to the physical space. By embedding PDE constraints into the training via automatic differentiation, our method eliminates the need for labeled training data and ensures physics-consistent predictions. The proposed framework is both memory and compute-efficient, exhibiting near-constant scaling with problem size and demonstrating significant speedups over traditional physics-informed operator models. We validate our approach on a range of parametric PDEs, showcasing its accuracy, scalability, and suitability for real-time prediction in complex physical systems.

Latent representations↗

DOE OSTI2025

Geometry-aware framework for deep energy method: An application to structural mechanics with hyperelastic materials

Here, in this work, we introduce a novel physics-informed framework named the Geometry-Aware Deep Energy Method (GADEM) for solving structural mechanics problems on different geometries. As the weak form of the physical system equation (or the energy-based approach) has demonstrated clear advantages compared to the strong form for solving solid mechanics problems, GADEM employs the weak form and aims to infer the solution on multiple shapes of geometries. Integrating a geometry-aware framework into an energy-based method results in an effective physics-informed deep learning model in terms of accuracy and computational cost. Different ways to represent the geometric information and to encode the geometric latent vectors are investigated in this work. We introduce a loss function of GADEM which is minimized based on the potential energy of all considered geometries. An adaptive learning method is also employed for the sampling of collocation points to enhance the performance of GADEM. We present some applications of GADEM to solve solid mechanics problems, including a loading simulation of a toy tire involving contact mechanics and large deformation hyperelasticity. The numerical results of this work demonstrate the remarkable capability of GADEM to infer the solution on various and new shapes of geometries using only one trained model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗