Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical integration”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗
DOE OSTI2025

WEIS and FAST.Farm Advancements Beyond Wind Turbine Aeroelasticity

Presentation at the 7th Wind Energy Systems Engineering workshop from NREL principal engineer Jason Jonkman, Ph.D., on NREL's numerical tools Wind Energy with Integrated Servo-control (WEIS), which focuses on the integrated design of floating wind turbines, and FAST.Farm, a simulation tool for multi-turbine wind farms based on popular aero-servo-hydro-elastic solver OpenFAST. The presentation discusses recent advancements beyond wind turbine aero-elasticity.

17 WIND ENERGY↗
DOE OSTI2025

Coupled Hydrological-Thermal-Biogeochemical Modeling for Predicting Arctic Carbon Emissions (CH4PACE)

The Arctic, with nearly 1,700 billion metric tons of carbon stored in its soils, is experiencing significant warming trends that raise concerns about the potential release of greenhouse gases, particularly methane. This report investigates the conversion of carbon stored in Arctic soils and permafrost into methane and its subsequent atmospheric release over the coming decades. Recent data indicate a rise in global methane concentrations, highlighting the importance of understanding the mechanisms of methane generation under anaerobic conditions prevalent in thawing permafrost. Our study employs a multi-faceted approach that integrates fieldwork, benchtop experiments, and numerical simulations to explore the subsurface hydrological-thermal-biogeochemical systems influencing methane emissions. We focus on the relationship between microbial activity in decomposing organic matter under oxygen-limited conditions and the transport of methane out of these saturated wetland soils. Key hydrological parameters, including permeability, porosity, and thermal conductivity, were measured on soil samples from Fairbanks, Alaska, and used to parameterize PFLOTRAN, a coupled hydrological-thermal-reactive transport simulator. Our model enables prediction of the evolution of the active layer and its impact on methane production and release. Our findings underscore the critical need for a mechanistic understanding of methane dynamics in permafrost environments to improve predictions of methane emissions in Earth systems models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Two-loop master integrals for leading-color $$ pp\to t\overline{t}H $$ amplitudes with a light-quark loop

Abstract We compute the two-loop master integrals for leading-color QCD scattering amplitudes including a closed light-quark loop in$$ t\overline{t}H $$ t t ¯ H production at hadron colliders. Exploiting numerical evaluations in modular arithmetic, we construct a basis of master integrals satisfying a system of differential equations inϵ-factorized form. We present the analytic form of the differential equations in terms of a minimal set of differential one-forms. We explore properties of the function space of analytic solutions to the differential equations in terms of iterative integrals which can be exploited for studying the analytic form of related scattering amplitudes. Finally, we solve the differential equations using generalized series expansions to numerically evaluate the master integrals in physical phase space. As the first computation of a set of two-loop seven-scale master integrals, our results provide valuable input for analytic studies of scattering amplitudes in processes involving massive particles and a large number of kinematic scales.

Physics↗
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↗
DOE OSTI2026

Double-real corrections to color singlet decay in a parton-shower inspired scheme

We introduce a local infrared subtraction method for next-to-next-to-leading order QCD calculations in color singlet decays, with counterterms based on scalar radiators and pure splitting functions. Overlapping singularities in the multipole radiation pattern are disentangled by partial fractioning, and the kinematics mapping corresponds to iterated next-to-leading order kinematics. We verify that the double-real remainder to $e^+e^-\to\;q\bar{q}$ is rendered finite in the single and double unresolved limits and investigate the numerical convergence of the Monte-Carlo integral. We compute the phase-space integrals of the scalar counterterms in the back-to-back configuration, both analytically and with the help of numerical techniques based on sector decomposition.

Campbell, John M. [Fermilab]↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗
DOE OSTI2026

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Studying phonon coherence with a quantum sensor

Nanomechanical oscillators offer numerous advantages for quantum technologies. Their integration with superconducting qubits shows promise for hardware-efficient quantum error-correction protocols involving superpositions of mechanical coherent states. Limitations of this approach include mechanical decoherence processes, particularly two-level system (TLS) defects, which have been widely studied using classical fields and detectors. In this manuscript, we use a superconducting qubit as a quantum sensor to perform phonon number-resolved measurements on a piezoelectrically coupled phononic crystal cavity. This enables a high-resolution study of mechanical dissipation and dephasing in coherent states of variable size ($\overline{n}$ ≃ 1 – 10 phonons). We observe nonexponential relaxation and state size-dependent reduction of the dephasing rate, which we attribute to TLS. Using a numerical model, we reproduce the dissipation signatures (and to a lesser extent, the dephasing signatures) via emission into a small ensemble (N = 5) of rapidly dephasing TLS. Our findings comprise a detailed examination of TLS-induced phonon decoherence in the quantum regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗
DOE OSTI2026

ZEUS: An Efficient GPU Optimization Method Integrating PSO, BFGS, and Automatic Differentiation

We introduce a novel, efficient computational method, ZEUS, for numerical optimization, and provide an open-source implementation. It has four key ingredients: (1) particle swarm optimization (PSO), (2) the use of the Broyden-Fletcher-Goldfarb-Shanno (BFGS) method, (3) automatic differentiation (AD), and (4) GPUs. Our approach addresses the computational challenges inherent in high-dimensional, non-convex optimization problems. In the first phase of the algorithm, we get a potentially good set of starting points using PSO. Thereafter, we run BFGS independently in parallel from these starting points. BFGS is one of the best-performing algorithms for numerical optimization. However, it requires the gradient of the function being optimized. ZEUS integrates automatic differentiation into BFGS thus avoiding the need for the user to calculate derivatives explicitly. The use of GPUs allows ZEUS to speed up the calculations substantially. We carry out systematic studies to explore the trade-offs between the number of PSO iterations taken, starting points, and BFGS iteration depth. We show that a handful of iterations of PSO can improve global convergence when combined with BFGS. We also present performance studies using common test functions. The source code can be found at https://github.com/fnal-numerics/global-optimizer-gpu.

Soos, Dominik [Old Dominion U.]↗
DOE OSTI2024

From oxide epitaxy to freestanding membranes: Opportunities and challenges

Motivated by the growing demand to integrate functional oxides with dissimilar materials, numerous studies have been undertaken to detach a functional oxide film from its original substrate, effectively forming a membrane, which can then be affixed to the desired host material. This review article is centered on the synthesis of functional oxide membranes, encompassing various approaches to their synthesis, exfoliation, and transfer techniques. First, we explore the characteristics of thin-film growth techniques with emphasis on molecular beam epitaxy. We then examine the fundamental principles and pivotal factors underlying three key approaches of creating membranes: (i) chemical lift-off, (ii) the two-dimensional layer–assisted lift-off, and (iii) spalling. We review the methods of exfoliation and transfer for each approach. Last, we provide an outlook into the future of oxide membranes, highlighting their applications and emerging properties.

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

NORTH DAKOTA CARBONSAFE PHASE III: SITE CHARACTERIZATION AND PERMITTING OF GEOLOGIC STORAGE OF CARBON DIOXIDE

The Energy & Environmental Research Center (EERC), in partnership with Minnkota Power Cooperative Inc. (Minnkota), SLB, and Computer Modelling Group Ltd. (CMG), supported wide-scale deployment of carbon capture and storage (CCS) as part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative Phase III. This phase included the acquisition, analysis, and development of information to fully characterize two storage complexes to demonstrate viable storage resources for commercial volumes of CO2 (defined by DOE as a minimum of 50 million tonnes [MMt] of CO2 within a 30-year period) (National Energy Technology Laboratory, 2024). Phase III also involved the preparation, submission, and approval of North Dakota underground injection control (UIC) Class VI storage facility permits (SFPs)—required precursors to applications for Class VI injection well permits. The presumed viability of commercial-scale CCS, situated adjacent to Minnkota’s Milton R. Young Station (MRYS), is validated by Minnkota’s continued pursuit of Project Tundra—an initiative to build the world’s largest lignite-based CCS project in central North Dakota (www.projecttundrand.com). Project Tundra comprises two scopes of work, Tundra Capture (installation of postcombustion CO2 capture at MRYS) and Tundra SGS (secure geologic storage). The efforts of North Dakota CarbonSAFE Phase III, Site Characterization and Permitting, supported Tundra SGS. Extensive site-specific characterization activities included a successful multimeasurement geophysical approach and drilling a stratigraphic test well (J-ROC 1, subsequently renamed Liberty-1) adjacent to MRYS. Core collection and analyses, downhole testing and fluid sampling, and geophysical logging were performed on J-ROC 1 and on a nearby stratigraphic test well (known as J-LOC 1), which was drilled, cored, and tested under a complementary project funded by the North Dakota Lignite Research Program. The injection tests performed on J-LOC 1 positively impacted the CarbonSAFE project, resulting in fewer proposed injection wells and significant construction, operations, and monitoring cost savings. The characterization data collected and analyses performed were integrated into geologic models, and successive numerical simulations were run to determine CO2 plume extent and subsurface pressure buildup associated with the planned CO2 injection rate of nearly 4 MMt per year. The latter doubles the CarbonSAFE Initiative goal with an estimated 100 MMt of CO2 stored in 20 years. Application of the U.S. Environmental Protection Agency’s (EPA’s) method for estimating the Class VI Rule area of review (AOR) to the overpressurized Broom Creek Formation inspired an alternative method of calculation, called risk-based AOR delineation. This peer-reviewed method was applied for the first time during the storage facility-permitting process. The two SFP applications submitted in 2021 successfully resulted in North Dakota Industrial Commission (NDIC) orders in 2022 authorizing the creation of the storage facility areas and amalgamation of pore space as well as establishing financial responsibility requirements. After approval of the SFPs, Minnkota filed in 2022 applications for permits to reenter the J-ROC 1 well and to drill two new wells—all with the intended purpose to become Class VI injection wells. To establish eligibility under the Internal Revenue Code for Section 45Q tax incentives, a monitoring, reporting, and verification (MRV) plan was prepared and submitted by Minnkota to EPA in 2022, resulting in the first such plan approved in North Dakota. Also in 2022, under the National Environmental Policy Act (NEPA), Minnkota prepared and submitted an environmental information volume (EIV) describing the proposed CCS project and associated potential environmental impacts. Based on the EIV, DOE determined that the proposed construction project required an environmental assessment, and Minnkota submitted the first draft in 2023 and a revised draft in 2024. Both submissions were followed by a public comment period. Subsequently, DOE issued a finding of no significant impact (FONSI) on September 13, 2024. A successful outreach program, strongly based in the production, presentation, and dissemination of informational material, fostered an environment to aid stakeholders in making informed decisions regarding the planned project. Opportunities for public input were provided at various steps along the way, including at county planning and zoning meetings, before and during the SFP administrative hearing, and during environmental assessment public comment periods. In addition, land/pore space owners and mineral owners had various points of contact, including granting access rights, securing pore space leasing, and mineral owner notifications. Based upon the successful storage facility permitting issued by NDIC, approval of the MRV plan by EPA, and receipt of a FONSI under the NEPA, Minnkota is continuing its pursuit of Project Tundra. In December 2023, the Office of Clean Energy Demonstrations under its Carbon Capture Demonstrations Projects Program announced funding for the capture system (Office of Clean Energy Demonstrations, 2023) and a proposal for CarbonSAFE Phase IV: Construction funding was submitted in March 2024 for the storage project. A go/no-go decision to proceed with construction and operations in the Broom Creek Formation is anticipated in 2024. References National Energy Technology Laboratory, CarbonSafe Initiative, https://netl.doe.gov/carbon-management/carbon-storage/carbonsafe (accessed August 2024). Office of Clean Energy Demonstrations, 2023, OCED selects three projects in CA, ND, and TX to reduce harmful carbon pollution, create new economic opportunities, and advance carbon reducing technologies, December, www.energy.gov/oced/articles/oced-selects-three-projects-ca-nd-and-tx-reduce-harmful-carbon-pollution-create-new (accessed August 2024).

Peck, Wesley↗
DOE OSTI2024

Mitigation for Removal of Building 01-103 in the Area 1 Subdock, Nevada National Security Site, Nye County, Nevada

Building 01-103, the Drill Bit Repair Building was originally recorded as SHPO Resource No. B18848 in 2020 during an architectural survey of the Area 1 Subdock (Collins et al. 2022). It was recommended individually eligible for listing in the NRHP under Criteria A and C and contributing to the significance of the Area 1 Subdock Historic District. The SHPO concurred with these determinations in a letter dated April 29, 2022 (Reed 2022). The Drill Bit Repair Building retains all seven aspects of integrity. The location, materials, workmanship, feeling, and association have not changed since its original construction. The design was modified slightly in 1996 with the addition of a removable and non-destructive chain-link fence. Nearby buildings 01-101 and 01-102 were removed (Collins and Menocal 2020) but the other buildings and storage yards remain. Therefore, the integrity of setting has been changed, but not to the degree that the building no longer conveys its significance (see Figures 4 and 5 for comparison). Building 01-103 was constructed in 1985 to coincide with the relocation of the Area 3 Subdock to Area 1. The design for the new Subdock provided a centralized cluster of buildings surrounded by two large storage yards. The buildings included administrative and engineering office space, warehouses, a radiograph building, and Building 01-103 as the drill bit maintenance and repair building, among others. The Drill Bit Repair Building was critical to the success of underground nuclear testing and is eligible under Criterion A for its contribution to broad patterns of history by maintaining specially designed equipment for underground nuclear testing at the NNSS. It is also eligible under Criterion C as it is the only recorded example of a facility built specifically for the maintenance of big-hole drill bits on the NNSS. As it retains both significance and integrity, it remains individually eligible under Criteria A and C. In addition to its individual significance, the Drill Bit Repair Building also contributes to the significance of the Area 1 Subdock Historic District as a representative of a big-hole drilling support facility necessary for underground nuclear testing. The building and the equipment therein were designed specifically for testing at the NNSS, where much of the early technology and methods for big-hole drilling were developed. Its removal will result in not only the demolition of a contributing element but also a change in the physical features of the district’s setting that contribute to its historic significance. The change in setting will diminish the district’s integrity. However, the district still contains numerous examples of other buildings constructed for the purpose of supporting big-hole drilling and remains eligible under Criteria A and C.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

LLM integration into EPICS

The utilization of large language models (LLMs) such as ChatGPT has seen a remarkable increase in various fields over the past few years. These models have demonstrated their versatility and capability in understanding and generating human-like text, making them invaluable tools in numerous applications. In this project, we explore the integration of a LLM into the Experimental Physics and Industrial Control System (EPICS). The primary focus of this integration is to employ the LLM for advanced image processing and spatial analysis on images obtained from the beamlines. By leveraging the capabilities of the LLM, we aim to enhance the accuracy and efficiency of image interpretation, enabling more precise data analysis and decision-making within the EPICS framework. This integration not only showcases the potential of LLMs in scientific and industrial applications but also sets the stage for future advancements in automated control systems.

Adams, Ethan↗
DOE OSTI2025

Confined and deconfined chaos in classical spin systems

Weakly perturbed integrable many-body systems are typically chaotic, and thermal at late times. However, there are distinct relationships between the timescales for thermalization and chaos. The typical relationship is confined chaos: when trajectories are still confined to regions in phase space with constant conserved quantities (actions), the conjugate angle variables are already unstable. Chaotic instabilities thus far precede thermalization. In a different relationship, which we term deconfined chaos, chaotic instabilities and thermalization occur on the same timescale. We investigate these two qualitatively distinct scenarios through numerical and analytical studies of two perturbed integrable classical spin models: the Ishimori spin chain (confined chaos), and the central spin model with XX interactions (deconfined chaos). We analytically establish (super)-integrability in the latter model in a microcanonical shell. Deconfined chaos emerges through the separation of phase space into large quasi-integrable regions and a thin chaotic manifold. The latter leads to chaos and thermalization on the fastest possible timescale, which is proportional to the inverse perturbation strength. This behavior is reminiscent of the quantum SYK models and strange metals.

Chaotic Dynamics (nlin.CD)↗
DOE OSTI2025

Computation of generalised magnetic coordinates asymptotically close to the separatrix

Integrals to calculate generalised magnetic coordinates from an input magnetic flux function asymptotically close to the separatrix are presented, and implemented in the GPEC/DCON code suite. These integrals allow characterisation of the magnetic equilibrium of a diverted tokamak, in magnetic coordinates, arbitrarily close to the last closed flux surface, avoiding the numerical issues associated with calculating diverging field-line integrals near a magnetic x-point. Finally, these methods may assist ongoing efforts to develop robust asymptotic equilibrium behaviour for spectral 3D MHD codes at the separatrix.

equilibrium edge truncation↗