Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Trajectory Generation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Viral but not bacterial community successional patterns reflect extreme turnover shortly after rewetting dry soils

As central members of soil trophic networks, viruses have the potential to drive substantial microbial mortality and nutrient turnover. Pinpointing viral contributions to terrestrial ecosystem processes remains a challenge, as temporal dynamics are difficult to unravel in the spatially and physicochemically heterogeneous soil environment. In Mediterranean grasslands, the first rainfall after seasonal drought provides an ecosystem reset, triggering microbial activity during a tractable window for capturing short-term dynamics. In this report we simulated precipitation in microcosms from four distinct dry grassland soils and generated 144 viromes, 84 metagenomes and 84 16S ribosomal RNA gene amplicon datasets to characterize viral, prokaryotic and relic DNA dynamics over 10 days. Vastly different viral communities in each soil followed remarkably similar successional trajectories. Wet-up triggered a significant increase in viral richness, followed by extensive compositional turnover. Temporal succession in prokaryotic communities was much less pronounced, perhaps suggesting differences in the scales of activity captured by viromes (representing recently produced, ephemeral viral particles) and total DNA. Still, differences in the relative abundances of Actinobacteria (enriched in dry soils) and Proteobacteria (enriched in wetted soils) matched those of their predicted phages, indicating viral predation of dominant bacterial taxa. Rewetting also rapidly depleted relic DNA, which subsequently reaccumulated, indicating substantial new microbial mortality in the days after wet-up, particularly of the taxa putatively under phage predation. Production of abundant, diverse viral particles via microbial host cell lysis appears to be a conserved feature of the early response to soil rewetting, and results suggest the potential for ‘Cull-the-Winner’ dynamics, whereby viruses infect and cull but do not decimate dominant host populations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Probabilistic error estimation for non-intrusive reduced models learned from data of systems governed by linear parabolic partial differential equations

This work derives a residual-based a posteriori error estimator for reduced models learned with non-intrusive model reduction from data of high-dimensional systems governed by linear parabolic partial differential equations with control inputs. It is shown that quantities that are necessary for the error estimator can be either obtained exactly as the solutions of least-squares problems in a non-intrusive way from data such as initial conditions, control inputs, and high-dimensional solution trajectories or bounded in a probabilistic sense. Here, the computational procedure follows an offline/online decomposition. In the offline (training) phase, the high-dimensional system is judiciously solved in a black-box fashion to generate data and to set up the error estimator. In the online phase, the estimator is used to bound the error of the reduced-model predictions for new initial conditions and new control inputs without recourse to the high-dimensional system. Numerical results demonstrate the workflow of the proposed approach from data to reduced models to certified predictions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Data-driven Alignment of the Tracker for the Mu2e Experiment

An alignment of the Mu2e Tracker is investigated using the Mu2e Offline Software and Millepede-II. Simulated misalignments of Tracker planes were inserted between Monte-Carlo event generation and before reconstruction using randomly generated $(x,y,z)$-translations and rotations of the rigid-body planes in the detector model. Cosmic ray muons, with no magnetic field present in the Detector Solenoid, will follow a straight-line trajectory. A maximum-likelihood fit is used to reconstruct a track and determine its parameters. Tracks traversing at least four planes, with ten registered straw hits, and a maximum time residual of $20$ ns are selected to be used in an alignment sample. The cuts were chosen to discard tracks with low information content, such as those with high levels of multiple scattering or those with low numbers of measurements contributing to a single track fit. Millepede-II produces the resulting alignment constants according to a simultaneous least-squares fit of the alignment parameters and the track parameters for all input tracks. The improved alignment constants are used to repeat the process and recover an increased number of suitable tracks, further improving the estimate of the alignment constants in subsequent iterations. A performance study of the process was conducted on a misaligned Tracker in the case of translation-only, and translation plus rotational misalignments. In both cases, the correct alignment constants were recovered to within statistical error after three global fit iterations. Additionally, a continuous integration system was developed for the \texttt{Mu2e/Offline} GitHub repository, introducing a convenient automated testing process for all proposed changes (``pull requests"). The tests developed included a full build of the software, followed by tests of the validation routines, geometry overlap checking, and multi-threading mode. Additionally, code quality checks and low-statistics physics comparison checks a re also available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Validation and Demonstration of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an integrated energy system (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. The new RAVEN plugin for grid and capacity optimization (HERON) only addresses the limits that affect the production variables and the corresponding rates of variation (explicit constraints). However, other variables are subject to constraints, and the associated limits should be accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response affects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. To obtain a more accurate approximation of the nonlinear dynamic behavior, a parametric version of the dynamic mode decomposition with control (DMDc) algorithm was developed to derive the state-space representation matrices of the IES components at different scheduling parameter. Thanks to this approach, a more accurate approximation of the system response can be obtained, the limits imposed by thermal mechanical implicit constraints can be translated into power dispatch limits, and the feedbacks to HERON power dispatcher can be provided. To assess the developed methodology, a power dispatching test case composed of three power generating and storage units (Balance of Plant, Secondary Energy Source, Thermal Energy Storage) was developed. The power output of each one of the three units was optimized to meet the imposed time-dependent load demand trajectory and to maximize the IES profitability by meeting both the explicit and implicit constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Toward a Generalizable Framework of Disturbance Ecology Through Crowdsourced Science

Disturbances fundamentally alter ecosystem functions, yet predicting their impacts remains a key scientific challenge. While the study of disturbances is ubiquitous across many ecological disciplines, there is no agreed-upon, cross-disciplinary foundation for discussing or quantifying the complexity of disturbances, and no consistent terminology or methodologies exist. This inconsistency presents an increasingly urgent challenge due to accelerating global change and the threat of interacting disturbances that can destabilize ecosystem responses. By harvesting the expertise of an interdisciplinary cohort of contributors spanning 42 institutions across 15 countries, we identified an essential limitation in disturbance ecology: the word ‘disturbance’ is used interchangeably to refer to both the events that cause, and the consequences of, ecological change, despite fundamental distinctions between the two meanings. In response, we developed a generalizable framework of ecosystem disturbances, providing a well-defined lexicon for understanding disturbances across perspectives and scales. The framework results from ideas that resonate across multiple scientific disciplines and provides a baseline standard to compare disturbances across fields. This framework can be supplemented by discipline-specific variables to provide maximum benefit to both inter- and intra-disciplinary research. To support future syntheses and meta-analyses of disturbance research, we also encourage researchers to be explicit in how they define disturbance drivers and impacts, and we recommend minimum reporting standards that are applicable regardless of scale. Finally, we discuss the primary factors we considered when developing a baseline framework and propose four future directions to advance our interdisciplinary understanding of disturbances and their social-ecological impacts: integrating across ecological scales, understanding disturbance interactions, establishing baselines and trajectories, and developing process-based models and ecological forecasting initiatives. Our experience through this process motivates us to encourage the wider scientific community to continue to explore new approaches for leveraging Open Science principles in generating creative and multidisciplinary ideas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

GANISP: A GAN-Assisted Importance Splitting Probability Estimator

Designing manufacturing processes with high yield and strong reliability relies on effective methods for rare event estimation. Genealogical importance splitting reduces the variance of rare event probability estimators by iteratively selecting and replicating realizations that are headed towards a rare event. The replication step is difficult when applied to deterministic systems where the initial conditions of the offspring realizations need to be modified. Typically, a random perturbation is applied to the offspring to differentiate their trajectory from the parent realization. However, this random perturbation strategy may be effective for some systems while failing for others, preventing variance reduction in the probability estimate. This work seeks to address this limitation using a generative model such as a Generative Adversarial Network (GAN) to generate perturbations that are consistent with the attractor of the dynamical system. The proposed GAN-assisted Importance SPlitting method (GANISP) improves the variance reduction for the system targeted. An implementation of the method is available in a companion repository (https://github.com/NREL/GANISP).

generative adversarial network↗
DOE OSTI2024

The OSIRIS-REx Sample Return Capsule re-entry: A coordinated seismo-acoustic observational campaign for the study of meteor phenomena

Objects entering Earth’s atmosphere at hypersonic velocities may generate powerful acoustic waves. Depending on atmospheric conditions and acoustic propagation paths, they may be captured using ground or balloon-borne microbarometers and microphones. Impacts into the Earth’s atmosphere by asteroids in a meter-size range are sporadic and unannounced, and thus, well-documented scientific observations of asteroids are rare and generally happen by chance. Arriving from interplanetary space at hypervelocity, spacecraft are considered analogues for low velocity meteoroids and asteroids impacting the Earth’s atmosphere, and as such provide unprecedented and unique opportunity to carry out planned observational campaigns and perform detailed studies of meteor phenomena. However, since the end of the Apollo era, only four instances of a hypersonic re-entry of an artificial body from interplanetary space with an incident speed of 11-12 km/s have been observed and studied. These were the Sample Return Capsules (SRCs) that brought physical samples of extraterrestrial material back to Earth (Genesis, Stardust, Hayabusa 1, Hayabusa 2). These re-entries were also detected by infrasound and/or seismic sensors. The next opportunity to observe an artificial meteor will be on 24 September 2023 with the re-entry of NASA’s OSIRIS-REx SRC that will bring samples of the carbonaceous near-Earth asteroid Bennu. The re-entry consists of several flight phases, including hypersonic and supersonic, providing a unique opportunity to observe these by infrasound and seismic sensors. We will discuss observational campaign efforts aimed to capture geophysical signals generated by the OSIRIS-REx SRC re-entry. This campaign utilizes an unparalleled number of instruments, both ground and airborne, strategically positioned in the immediate and extended region around the projected re-entry trajectory to maximize the scientific output.

47 - OTHER INSTRUMENTATION↗
DOE OSTI2020

Experimental Demonstration of Hadron Beam Cooling Using Radio-Frequency Accelerated Electron Bunches

Cooling of beams of gold ions using electron bunches accelerated with radio-frequency systems was recently experimentally demonstrated in the Relativistic Heavy Ion Collider at Brookhaven National Laboratory. Such an approach is new and opens the possibility of using this technique at higher energies than possible with electrostatic acceleration of electron beams. The challenges of this approach include generation of electron beams suitable for cooling, delivery of electron bunches of the required quality to the cooling sections without degradation of beam angular divergence and energy spread, achieving the required small angles between electron and ion trajectories in the cooling sections, precise velocity matching between the two beams, high-current operation of the electron accelerator, as well as several physics effects related to bunched-beam cooling. Here we report on the first demonstration of cooling hadron beams using this new approach.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Investigation of plasmon relaxation mechanisms using nonadiabatic molecular dynamics

Hot carriers generated from the decay of plasmon excitation can be harvested to drive a wide range of physical or chemical processes. However, their generation efficiency is limited by the concomitant phonon-induced relaxation processes by which the energy in excited carriers is transformed into heat. However, simulations of dynamics of nanoscale clusters are challenging due to the computational complexity involved. Here, in this paper, we adopt our newly developed Trajectory Surface Hopping (TSH) nonadiabatic molecular dynamics algorithm to simulate plasmon relaxation in Au 20 clusters, taking the atomistic details into account. The electronic properties are treated within the Linear Response Time-Dependent Tight-binding Density Functional Theory (LR-TDDFTB) framework. The relaxation of plasmon due to coupling to phonon modes in Au 20 beyond the Born–Oppenheimer approximation is described by the TSH algorithm. The numerically efficient LR-TDDFTB method allows us to address a dense manifold of excited states to ensure the inclusion of plasmon excitation. Starting from the photoexcited plasmon states in Au 20 cluster, we find that the time constant for relaxation from plasmon excited states to the lowest excited states is about 2.7 ps, mainly resulting from a stepwise decay process caused by low-frequency phonons of the Au 20 cluster. Furthermore, our simulations show that the lifetime of the phonon-induced plasmon dephasing process is ~10.4 fs and that such a swift process can be attributed to the strong nonadiabatic effect in small clusters. Our simulations demonstrate a detailed description of the dynamic processes in nanoclusters, including plasmon excitation, hot carrier generation from plasmon excitation dephasing, and the subsequent phonon-induced relaxation process.

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

Wassersplines for Neural Vector Field-Controlled Animation

Much of computer-generated animation is created by manipulating meshes with rigs. While this approach works well for animating articulated objects like animals, it has limited flexibility for animating less structured free-form objects. Here we introduce Wassersplines, a novel trajectory inference method for animating unstructured densities based on recent advances in continuous normalizing flows and optimal transport. The key idea is to train a neurally-parameterized velocity field that represents the motion between keyframes. Trajectories are then computed by advecting keyframes through the velocity field. We solve an additional Wasserstein barycenter interpolation problem to guarantee strict adherence to keyframes. Our tool can stylize trajectories through a variety of PDE-based regularizers to create different visual effects. We demonstrate our tool on various keyframe interpolation problems to produce temporally-coherent animations without meshing or rigging.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Method for Spatiotemporal Solar Power Profile Estimation for a Proposed U.S.–Caribbean–South America Super Grid under Hurricanes

Solar photovoltaic (PV) generation technology stands out as a scalable and cost-effective solution to enable the transition toward decarbonization. However, PV solar output, beyond the daily solar irradiance variability and unavailability during nights, is very sensitive to weather events like hurricanes. Hurricanes nucleate massive amounts of clouds around their centers, shading hundreds of kilometers in their path, reducing PV power output. This research proposes a spatiotemporal method, implemented in MATLAB R2023b coding, to estimate the shading effect of hurricanes over a wide distribution of PV solar plants connected to a high-voltage power infrastructure called the U.S.–Caribbean–South America super grid. The complete interconnection of the U.S., the Caribbean, and South America results in the lowest power valley levels, i.e., an overall percentual reduction in PV power output caused by hurricane shading. The simulations assess the impact of hurricanes in 10 synthetic trajectories spanning from Texas to Florida. The Caribbean would also experience lower power valleys with expanded interconnectivity schemes. The U.S.–Caribbean–South America super grid reduces Caribbean variability from 37.8% to 8.9% in the case study. The proposed spatiotemporal method for PV power profile estimation is a valuable tool for future solar power generation expansion, transmission planning, and system design considering the impact of hurricanes.

14 SOLAR ENERGY↗
DOE OSTI2023

Iridescence from Total Internal Reflection at 3D Microscale Interfaces: Mechanistic Insights and Spectral Analysis

An experimental investigation and the optical modeling of the structural coloration produced from total internal reflection interference within 3D microstructures are described. Ray-tracing simulations coupled with color visualization and spectral analysis techniques are used to model, examine, and rationalize the iridescence generated for a range of microgeometries, including hemicylinders and truncated hemispheres, under varying illumination conditions. An approach to deconstruct the observed iridescence and complex far-field spectral features into its elementary components and systematically link them to ray trajectories that emanate from the illuminated microstructures is demonstrated. The results are compared with experiments, wherein microstructures are fabricated with methods such as chemical etching, multiphoton lithography, and grayscale lithography. Microstructure arrays patterned on surfaces with varying orientation and size lead to unique color-traveling optical effects and highlight opportunities for how total internal reflection interference can be used to create customizable reflective iridescence. The findings herein provide a robust conceptual framework for rationalizing this multibounce interference mechanism and establish approaches for characterizing and tailoring the optical and iridescent properties of microstructured surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

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

Non-Reciprocity, Metastability, and Dynamic Reconfiguration in Co-Assembly of Active and Passive Particles

Living organisms often exhibit non-reciprocal interactions where the forces acting on the objects are not equal in magnitude or opposite in direction. The combination of reciprocal and non-reciprocal interactions between synthetic building blocks remains largely unexplored. Here, out-of-equilibrium assemblies of non-motile isotropic passive and metal-patched motile active particles are formed by overlapping bulk interactions with directed self-propulsion. An external alternating current (AC) electric field generates concurrent dipolar and induced-charge electrophoretic forces between the particles which are evaluated using microscopy. The interaction force measurements allow to determine the degree of reciprocity in interactions, which is tunable by designing the active particle and its trajectory. While linearly-propelled active particles evade assembly with passive particles, helically propelled active particles form active-passive clusters with dynamic reconfiguration and long-lived metastability. Large clusters display programmable fluctuations and reconfigurability by controlling the fraction of active particles. The study establishes principles of integrating reciprocal and non-reciprocal interactions in guided colloidal assembly of reconfigurable metastable structures.

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

Augmented signal processing in Liquid Argon Time Projection Chambers with a deep neural network

The Liquid Argon Time Projection Chamber (LArTPC) is an advanced neutrino detector technology widely used in recent and upcoming accelerator neutrino experiments. It features a low energy threshold and high spatial resolution that allow for comprehensive reconstruction of event topologies. In current-generation LArTPCs, the recorded data consist of digitized waveforms on wires produced by induced signal on wires of drifting ionization electrons, which can also be viewed as two-dimensional (2D) (time versus wire) projection images of charged-particle trajectories. For such an imaging detector, one critical step is the signal processing that reconstructs the original charge projections from the recorded 2D images. For the first time, we introduce a deep neural network in LArTPC signal processing to improve the signal region of interest detection. By combining domain knowledge (e.g., matching information from multiple wire planes) and deep learning, this method shows significant improvements over traditional methods. This work details the method, software tools, and performance evaluated with realistic detector simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Life cycle assessment of utilizing freshly cut urban wood: A case study

Large amounts of urban wood resources are generated through tree removals in urban areas. Therefore, there is a growing interest to improve the environmental performance of the urban wood supply chain by quantifying the environmental impacts of creating high-value products from urban tree removals to enhance its utilization. By surveying existing urban wood utilization operations for primary data in two major cities: Baltimore, Maryland, and Milwaukee, Wisconsin, a new life cycle inventory (LCI) dataset per m 3 of dry urban hardwood lumber and live edge slabs was created. Incorporating the new LCI data into life cycle assessment method, the total global warming (GW) impact of converting urban trees to kiln-dried lumber was estimated to be 122 kg CO 2 eq and 336 kg CO 2 eq per m 3 of lumber produced for Baltimore and Milwaukee, respectively. In both cases, the total environmental impact of the product was dominated by the conversion of urban trees to kiln dried lumber at the sawmill processing stage. Using scenario analysis for the source of heat generation used in kiln (force) drying wood, substituting wood fuel from coproducts for natural gas in a boiler for Milwaukee resulted in a substantial reduction in fossil-derived GW impact. Urban and traditional hardwood lumber production follow similar GW impact trajectories over their life-cycle stages with the lumber processing stage having the highest environmental impacts for both. Finally, increasing the overall lumber yield and decreasing kiln (force) drying would substantially improve the environmental performance of urban hardwood lumber and help make it comparable to traditional lumber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluation of radiation hardness of semiconductor materials against alpha particles for an API detector

An Associated Particle Imaging (API) system consists of a high vacuum chamber where deuterium ions are accelerated on a tritium-loaded target, resulting in alpha particles and neutrons emitted back-to-back, the latter escaping the chamber and used, for example, for tomography of a high-Z sample. An accurate determination of the alpha position is crucial to determine the trajectory of the neutron. Existing API systems have several limitations which a semiconductor-based API detector placed inside the vacuum chamber should not present. The semiconductor material can be either silicon or diamond. In particular, large and fast signals are generated by the alpha interactions in these materials. A pixelated API detector will measure the time and hit position of the alphas produced in the reaction. However, during the lifetime of the generator, the semiconductor detector will be exposed to an intense flux of alpha particles that will degrade its performance over time. To assess the radiation hardness of silicon and diamond against the alpha particles, we exposed single diodes made of either of the two materials to a 5 MeV alpha flux generated by an 241Am radioactive source. During irradiation, the diodes were biased and mounted on current-sensitive preamplifier boards. We have therefore been able to measure insitu the evolution of the waveforms as the damage was accumulating in their substrates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Precise stellarator quasi-symmetry can be achieved with electromagnetic coils

Magnetic fields with quasi-symmetry are known to provide good confinement of charged particles and plasmas, but the extent to which quasi-symmetry can be achieved in practice has remained an open question. Recent work [M. Landreman and E. Paul, Phys. Rev. Lett. 128, 035001, 2022] reports the discovery of toroidal magnetic fields that are quasi-symmetric to orders-of-magnitude higher precision than previously known fields. We show that these fields can be accurately produced using electromagnetic coils of only moderate engineering complexity, that is, coils that have low curvature and that are sufficiently separated from each other. Our results demonstrate that these new quasi-symmetric fields are relevant for applications requiring the confinement of energetic charged particles for long time scales, such as nuclear fusion. The coils’ length plays an important role for how well the quasi-symmetric fields can be approximated. For the longest coil set considered and a mean field strength of 1 T, the departure from quasi-symmetry is of the order of Earth’s magnetic field. Additionally, we find that magnetic surfaces extend far outside the plasma boundary used by Landreman and Paul, providing confinement far from the core. Simulations confirm that the magnetic fields generated by the new coils confine particles with high kinetic energy substantially longer than previously known coil configurations. In particular, when scaled to a reactor, the best found configuration loses only 0.04% of energetic particles born at midradius when following guiding center trajectories for 200 ms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗