Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamics testing”

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 919 records · Page 51

Exploitation of Defects in High Entropy Ceramic Barrier Materials

A critical mission need exists to develop new materials that can withstand extreme environments and multiple sequential threats. High entropy materials, those containing 5 or more metals, exhibit many exciting properties which would potentially be useful in such situations. However, a particularly hard challenge in developing new high entropy materials is determining a priori which compositions will form the desired single phase material. The project outlined here combined several modeling and experimental techniques to explore several structure-property-relationships of high entropy ceramics in an effort to better understand the connection between their compositional components, their observed properties, and stability. We have developed novel machine learning algorithms which rapidly predict stable high entropy ceramic compositions, identified the stability interplay between configurational entropy and cation defects, and tested the mechanical stability of high entropy oxides using the unique capabilities at the Dynamic Compression Sector facility and the Saturn accelerator.

36 MATERIALS SCIENCE↗

Halo Nuclei from Ab Initio Nuclear Theory

A realistic description of halo nuclei, characterized by low-lying breakup thresholds, requires a proper treatment of continuum effects. We have developed an ab initio approach, the No-Core Shell Model with Continuum (NCSMC), capable of describing both bound and unbound states in light nuclei in a unified way. With chiral two- and three-nucleon interactions as the only input, we can predict the structure and dynamics of halo and other light nuclei and, by comparing to available experimental data, test the quality of chiral nuclear forces. We review NCSMC calculations of weakly bound states and resonances of the exotic halo nuclei 6He, 8B, 11Be, and 15C. For the latter, we discuss its production in the capture reaction 14C(n,𝛾 )15C. We highlight the challenges of a description of 6He as a Borromean n-n-4He system. Finally, we present our calculations of excited states in 10Be exhibiting a one-neutron halo structure and a large scale No-Core Shell Model investigation of 11Li as a precursor of a full n-n-9Li NCSMC study.

Navrátil, Petr↗

Model, Design, and Control for Power Conversion in Wave Energy Converter System

Wave energy, if all been harvested along the U.S. coastline, can power around 65% of the energy consumption in U.S.. Comparing to other renewable energy sources like solar or wind, ocean wave can provide up to 90% of steady uptime. With the high energy density (2-3 kW/m 2 ), it can produce more energy with the same amount of installation area comparing to the energy density of wind turbine (0.6 kW/m 2 ) and solar panel(0.2 kW/m 2 ). The predictability of wave provides advantages like planning installation, power dispatching, and maintenance activities. Although with all these advantages, wave energy converter system is still in the research stage due to its high system cost per electricity production. One of the challenges that need to be solved is the irregularity from the wave motion that leads to high instantaneous peak power into the wave energy converter, which usually reaches up to 10 - 20 times of the average power. The high peak power will not only bring high mechanical/electrical stress but also result in an overrating design of the components in the system. Another obstacle that prevents the wave energy system from moving forward is the high testing cost from the validations in wave-energy-test sites or tank-test sites. A high-fidelity multi-disciplinary system model, including hydrodynamics, mechanical dynamics, electromagnetics, and power electronics, is needed to predict the behavior of the system and reduce the cost of design validation.

16 TIDAL AND WAVE POWER↗

Modeling of H2 Dispersion at ARIES

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower and help decarbonize electricity grids, industry, and transportation. Using the Hydrogen Research Facility under Advanced Research on Integrated Energy Systems (ARIES) at the National Renewable Energy Laboratory's (NREL) Flatirons campus as a test bench, the study examines the feasibility, useability, and value of using computational fluid dynamics (CFD) techniques to model hydrogen dispersion. The ARIES facility was chosen because controlled hydrogen releases can be performed at a rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind speed and temperature were used as inputs to the model. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. Wind conditions are found to significantly impact the release behavior, including the hydrogen cloud's direction and concentrations. At low wind speeds (below 1 mph), hydrogen forms a cloud and at higher wind speeds (> 2-4 mph) hydrogen plume stretches in the direction of wind momentum. From >100 simulations for ARIES site-specific conditions, statistical quantities combined with a clustering algorithm were used to propose sensor location at various elevations from ground.

dispersion↗

Frequency Control and Disturbance Containment Using Grid-Forming Embedded Storage Networks

The paper presents a distributed approach for operating a network of inverter-based energy storage resources embedded in a bulk power system. Departing from their traditional role of steady-state reserves, the storage assets in the network are utilized as frequency-responsive resources shaping system dynamics. The power electronics converter systems interfacing the storage resources are equipped with local controllers designed to respond under transient disturbances. To this end, a safety-constrained distributed control strategy is explored. The paper compares the performance of converter-interfaced grid-forming and grid-following storage networks for fast frequency control and disturbance containment/localization. Sensitivity studies are performed to study the impact of storage size, steady-state dispatch, and controller design on dynamic performance. The findings are presented through case studies with results from the IEEE test systems.

Chatterjee, Kaustav (ORCID:0000000153273860)↗

Coupled Long-Term Experiment and Model Investigation of the Differential Response of Plants and Soil Microbes in a Changing Permafrost Tundra Ecosystem

The major research goal of this project was to understand and quantify the fate of carbon stored in permafrost ecosystems using a combination of field and laboratory experiments to measure and model isotope ratios and carbon fluxes in a tundra ecosystem exposed to experimental warming. Field measurements centered on a two-factor experimental warming to increase air and soil temperatures alone, and in combination, at a tundra field site in Alaska. A second manipulation of water table was embedded within the warming treatment such that both major environmental factors (temperature, moisture) controlling ecosystem carbon dynamics were experimentally altered. The experiment was interfaced with modeling activities using both data assimilation and forward modeling approaches. Models were used to make forecasts of ecosystem carbon dynamics beyond the time frame and environmental space of the experiment itself. This project was a follow-on to previous work that made use of the same field manipulation; here we continued the experimental manipulation in order to test new hypotheses regarding long-term effects of warming and permafrost thaw and subsequent changes in moisture availability on ecosystem carbon dynamics. These results were linked with ongoing synthesis and model intercomparison activities through the Permafrost Carbon Network, which led to additional synthesis publications as a result of collaboration and data sharing from this project. As key outcomes, this project has produced new papers published in the peer-reviewed literature, archived datasets used for model-data intercomparisons, trained graduate students and postdoctoral researchers, and has raised awareness about the vulnerability of permafrost carbon to a wider audience.

54 ENVIRONMENTAL SCIENCES↗

DEFECT DETECTION USING DYNAMIC ANALYSIS FOR ADDITIVE MANUFACTURED METALS

Additive manufacturing (AM) has the ability to produce parts with complex geometries and internal features, however, for demanding applications such as the automotive and aerospace industries, it is crucial that the parts can meet the demanding functional and geometric requirements. Quality control for AM parts focuses on nondestructive methods of testing, but many of the current methods are expensive and time-consuming. The research presented in this report explores various methods of nondestructive evaluation (NDE) using dynamic analysis on stainless steel parts produced with selective laser melting (SLM). Methods include, but are not limited to, frequency response functions (FRF), impedance-based measurements, and scanning laser doppler vibrometry. Additionally, mode shape analysis was performed in MATLAB and FEA simulations were used for comparison with experimental results. The results indicate that dynamic analysis has the potential to be a feasible method of defect detection and NDE in AM parts and future work should focus on refining these methods, such as optimizing test parameters to improve sensitivity to defects.

Deonarain, Gita↗

Multi-Timescale Integrated Dynamics and Scheduling for Solar (MIDAS-Solar) (Final Technical Report)

Solar photovoltaic (PV) installations have experienced unprecedented growth in the United States. PV will become not only an energy producer but also a necessary provider of ancillary services at multiple timescales. Conventional methods to simulate power systems operations - such as long-term production simulation (which typically considers schedules from hours to minutes by using an optimization framework) and short-term transient studies (which simulate dynamics from seconds to sub-seconds using state variables and differential equations) - are not sufficient for studying the multiple-timescale variation of solar generation and its impact on system reliability. Long-term system economics and short-term system dynamics are highly coupled, particularly when the penetration level of renewable generation is extremely high, because the uncertainty and variability of solar generation will impact both power system steady-state and dynamic performance. This project helps meet and exceed the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Solar Energy Technologies Office goal of systems integration by directly addressing this stability and reliability challenge for power grid planning and operation. We have developed a temporally comprehensive, closed-loop simulation model, named Multi-timescale Integrated Dynamics and Scheduling (MIDAS), that seamlessly simulates power system operations from economic scheduling (day-ahead to hours) to dynamic response analysis (seconds to sub-seconds). For schedules with very high levels of inverter-based resources (IBRs), up to and including 100%, the stability of grid controls has been evaluated through electromagnetic transient (EMT) simulations and power-hardware-in-the-loop (PHIL) simulations of key transient events at key schedule points. Specifically, MIDAS provides: 1) a closed-loop simulation framework for simulating timescales from economic scheduling to dynamic stability analysis; 2) machine learning-based stability assessment; 3) EMT modeling and analysis for large-scale power systems; 4) MIDAS PHIL test bed. We worked with Hawaii Electric Companies to apply the MIDAS study framework to a Maui grid study. The entire island's transmission system was modeled in detail - from a yearly scheduling model, to a second-level frequency dynamic model, down to a sub-second-scale EMT model to address critical stability issues. The project demonstrated how MIDAS can help system planners and operators assess system reliability and stability while the power grid is marching toward a high-renewable, high-IBR future. In this Maui grid study, we found that 100% instantaneous IBR operation is achievable in EMT simulation and PHIL testing, and grid planners and operators might need new analysis/simulation tools to assess grid reliability and stability in the scheduling stage. MIDAS will bring Maui and other systems closer to 100% clean and stable energy futures. (In this study, we examined transient stability. Other topics necessary for 100% IBR operation, such as protection and resource adequacy, were not examined.)

100% Renewables↗

In situ multimodal imaging for nanoscale visualization of tribofilm formation

A third of the energy from fuel combustion in passenger car gasoline engines is lost due to friction. Carefully designed engine lubricants can recover some of these losses by reducing friction and wear by forming a nanometer-scale chemico-physico tribofilm between surfaces. Accordingly, attention has focused on developing oil formulations that form low-friction tribofilms. However, analyses of resultant tribofilms are typically conducted after tribo-tests with conventional characterization tools and do not offer insights into tribofilm formation and evolution, precluding information critical to tuning tribofilm properties. In this work, we developed a unique multimodal methodology based on Atomic Force Microscopy (AFM) with local probe heating for in situ tribological studies that activates friction modifiers and simultaneously captures the evolution of friction and surface roughness, with nanometer resolution. As a platform to demonstrate the ability of this methodology to visualize dynamics of tribofilm formation in situ, we apply it to molybdenum-based friction modifiers to distinguish key factors in their functionality and correlate nanoscale AFM and Friction Force Microscopy data to bench tribo-tests used in the industry. To decode the formation mechanisms observed in situ and underlying chemistry of tribofilms, we performed ab initio Molecular Dynamics (AIMD) simulations at comparable conditions. AIMD simulations confirmed both nanoscale and bench tribo-test results and showed deviations in molecular organization in tribofilms that are formed from different molybdenum-based friction modifiers that corroborates with surface functionality. With this innovative methodology, we demonstrate proof-of-principle in situ formation of molybdenum-based tribofilms directly on steel surfaces that could be applied generally to studying tribofilm formation.

36 MATERIALS SCIENCE↗

Stability and pulsation of the first dark stars

ABSTRACT The first bright objects to form in the Universe might not have been ‘ordinary’ fusion-powered stars, but ‘dark stars’ (DSs) powered by the annihilation of dark matter (DM) in the form of weakly interacting massive particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of the order of M⊙ and may grow to a few million solar masses; in this work we investigate the properties of early DSs with masses up to $\sim \! 1000 \, \mathrm{ M}_\odot$, fueled by WIMPS weighing 100 GeV. We improve the previous implementation of the DM energy source into the stellar evolution code mesa. We show that the growth of DSs is not limited by astrophysical effects: DSs up to $\sim \!1000 \, \mathrm{ M}_{\odot }$ exhibit no dynamical instabilities; DSs are not subject to mass-loss driven by super-Eddington winds. We test the assumption of previous work that the injected energy per WIMP annihilation is constant throughout the star; relaxing this assumption does not change the properties of the DSs. Furthermore, we study DS pulsations, for the first time investigating non-adiabatic pulsation modes, using the linear pulsation code gyre. We find that acoustic modes in DSs of masses smaller than $\sim \! 200 \, \mathrm{ M}_\odot$ are excited by the κ − γ and γ mechanism in layers where hydrogen or helium is (partially) ionized. Moreover, we show that the mass-loss rates potentially induced by pulsations are negligible compared to the accretion rates.

79 ASTRONOMY AND ASTROPHYSICS↗

Photovoltaic Analysis and Response Support (PARS) Platform for Solar Situational Awareness and Resiliency Services

The project's primary objective is to develop a digital-twin based Photovoltaic (PV) Analysis and Response Support (PARS) platform, which aims to provide real-time situational awareness and optimal response plans. This platform is designed to enhance the performance of hybrid PV systems, making them competitive with or even superior to conventional generation resources. The PARS platform enabled the project team to develop and evaluate an extensive suite of grid support functionalities for the hybrid PV systems to enhance grid performance, across key areas including visibility, dispatchability, security, resilience, and reliability. Given the global push toward achieving 100% clean energy by 2035, there is a significant increase in the integration of inverter-based resources (IBRs) throughout the energy grid. Effectively managing the inherent variability and uncertainty associated with IBRs is crucial for ensuring cost-effectiveness, reliability, and security in both the main grid and islanded microgrids. Constrained to a limited array of IEEE test systems or standard feeder models, traditional IBR modeling struggles to assimilate new field data, accurately reflect system dynamics, and adapt to the evolving energy landscape. In our project, we embraced a Digital Twin (DT) strategy for crafting the PARS platform. A digital twin acts as a precise virtual counterpart of a physical system, built on historical data and continuously honed with real-time insights. This enables the high-fidelity DT to accurately mirror current system operations and forecast future scenarios. Consequently, the PARS platform becomes an ideal environment for testing and refining monitoring, control, power, and energy management algorithms designed to boost hybrid PV system performance. The defining feature of the PARS platform, distinguishing it from other advanced simulation tools, is its exceptional adaptability. This is achieved by employing actual network topologies and utilizing real-time field data for fine-tuning and calibration, ensuring a close emulation of real-world conditions. The project deliverables include: 1) High-fidelity IBR models and tools for real-time parameterization, utilizing real-time field measurements to refine IBR models for enhanced accuracy and performance; 2) Grid-forming and Grid-following capabilities to deliver resilience services, including blackstart, voltage and frequency support, cold-load pick-up, power reserves, and three-phase load balancing across grid-connected and microgrid settings; 3) Machine learning-based forecasting tools and methods for generating synthetic data and topologies, creating diverse and realistic simulation environments for evaluating varied operational scenarios; 4) Advanced microgrid power and energy management algorithms for optimizing the integration and operation of PV, storage, and demand response resources within both feeder and community scales. The power grid data sets are provided by four utility companies in North Carolina and the New York Power Administration. Acting as industry advisors, our industry partners communicated stakeholder needs and regulatory standards to the research teams, aiding technology transfer by incorporating the developed methodologies into their daily operations. This collaboration ensures that the PARS platform, functioning as a power system digital twin, enhances our understanding of IBR dynamic behaviors and enables the development and evaluation of IBR control functions that match or exceed the capabilities of conventional synchronous generators.

14 SOLAR ENERGY↗

CICE Consortium Progress Report, June 2018 - December 2019

The CICE Consortium accomplished a number of major upgrades to our sea ice model code, scripts and testing infrastructure, while strengthening our visibility in the international modeling community. We released CICE versions 6.0.0, 6.0.1, 6.0.2 and 6.1.0, all of which include updated releases of Icepack. Major new capabilities include a floe size distribution, new coupling infrastructure, improvements to all of CICE’s dynamics options, a new atmospheric data set for historical simulations, overhauled documentation, and more robust and flexible scripts for running and testing the code. Several key testing processes were automated this year. In addition, the team is actively engaged in planning our first user workshop and tutorial in February 2020. In addition to overseeing modifications to the Consortium’s resources and code releases, management tasks included meeting with and updating the Consortium’s Sponsors group in October 2018, prioritizing tasks by mapping them onto potential release targets, and rescuing early CICE releases (version 1 through 3.14) from a disappearing LANL archive.

54 ENVIRONMENTAL SCIENCES↗

Model reduction for hypersonic aerodynamics via conservative LSPG projection and hyper-reduction

High-speed aerospace engineering applications rely heavily on computational fluid dynamics (CFD) models for design and analysis due to the expense and difficulty of flight tests and experiments. This reliance on CFD models necessitates performing accurate and reliable uncertainty quantification (UQ) of the CFD models. However, it is very computationally expensive to run CFD for hypersonic flows due to the fine grid resolution required to capture the strong shocks and large gradients that are typically present. Furthermore, UQ approaches are “many-query” problems requiring many runs with a wide range of input parameters. One way to enable computationally expensive models to be used in such many-query problems is to employ projection-based reduced-order models (ROMs) in lieu of the (high-fidelity) full-order model. In particular, the least-squares Petrov–Galerkin (LSPG) ROM (equipped with hyper-reduction) has demonstrated the ability to significantly reduce simulation costs while retaining high levels of accuracy on a range of problems including subsonic CFD applications. This allows computationally inexpensive LSPG ROM simulations to replace the full-order model simulations in UQ studies, which makes this many-query task tractable, even for large-scale CFD models. This work presents the first application of LSPG to a hypersonic CFD application. In particular, we present results for LSPG ROMs of the HIFiRE-1 in a three-dimensional, turbulent Mach 7.1 flow, showcasing the ability of the ROM to significantly reduce computational costs while maintaining high levels of accuracy in computed quantities of interest.

Fluid Dynamics↗

Weak-Form Latent Space Dynamics Identification

This software showcases the enhanced capabilities of the Latent Space Dynamics Identification (LaSDI) algorithm through the application of the weak form, resulting in WLaSDI. WLaSDI first compresses the data, then projects it onto test functions, and subsequently learns the local latent space models. Notably, WLaSDI demonstrates significantly improved robustness to noise. Using weak-form equation learning techniques, WLaSDI achieves local latent space modeling. Compared to the standard sparse identification of nonlinear dynamics (SINDy) used in LaSDI, the variance reduction of the weak form ensures robust and precise latent space recovery, enabling fast, robust, and accurate simulations. We demonstrate the efficacy of WLaSDI against LaSDI using several common benchmark examples, including viscid and inviscid Burgers', radial advection, and heat conduction. For instance, in 1D inviscid Burgers' simulations with up to 100% Gaussian white noise, WLaSDI maintains relative errors consistently below 6%, whereas LaSDI errors can exceed 10,000%. Similarly, in radial advection simulations, WLaSDI keeps relative errors below 16%, compared to potential errors of up to 10,000% with LaSDI. Additionally, WLaSDI achieves significant speedups, such as a 140X speedup in 1D Burgers' simulations compared to the corresponding full order model.

Choi, Youngsoo↗

Scalable modular dynamic molten salt reactor system model with decay heat

Nodal dynamic models allow rapid reduced-order simulations of complex systems. This paper presents a publicly available system dynamic model of a molten salt reactor system with a topology of the Molten Salt Reactor Experiment (MSRE), in MATLAB/Simulink, building upon a model previously published and validated with MSRE data. It adds scaling of nominal power, dynamic representation of decay heat generation, a generic decay heat removal system, improved documentation, modularity, and test cases. One hour of system transient runs in less than 10 min on a laptop. This model is intended for first order engineering calculations, developing insights into system transient behavior, equipment sizing, parameter sensitivity studies. The model contains an easy to use toolkit which can be customized and expanded to describe any molten salt reactor system. Finally, scenarios involving off-normal transients are presented. Relation between shutdown reactivity insertions, decay heat removal, core recriticality, and safe shutdown conditions are discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Pivotal Roles of Triple Screening-Topological, Electrostatic, and Hydrodynamic-On Dynamics in Semidilute Polyelectrolyte Solutions

For semidilute polyelectrolyte solutions, it is generally assumed that topological, electrostatic, and hydrodynamic interactions are screened (called triple screening). Despite a large body of research focused on polyelectrolyte solutions, the concept of triple screening has never been rigorously verified. In this work, we test the concept by probing concentration fluctuations in aqueous solutions containing a well-studied polyelectrolyte, sodium poly(styrenesulfonate) (NaPSS) with neutron scattering, theory, and molecular dynamics simulations. Neutron spin–echo (NSE) and small-angle neutron scattering (SANS) data from semidilute solutions of NaPSS are presented at different polymer and salt (NaCl) concentrations. A combined theory for structure (J. Chem. Phys. 105, 5183 (1996)) and dynamics (J. Chem. Phys. 107, 2619 (1997)), which captures effects of hydrodynamic, topological, and electrostatic screening, is used to interpret the experimental results. Here, the theory quantitatively predicts the decay rate obtained from the NSE measurements while capturing the shape and concentration dependencies of the polyelectrolyte peak observed in the SANS spectra. Detailed comparisons of the theory and the experiments reveal that the wavevector-dependent decay rate of concentration fluctuations in semidilute solutions of polyelectrolytes is dictated by the screening of hydrodynamic, topological, and electrostatic interactions. This conclusion is corroborated by coarse-grained molecular dynamics simulations, executed without any hydrodynamic interactions, which fail to capture the correct wavevector dependence of the decay rate. These results highlight that the theories based on the concept of triple screening provide a quantitative framework for predicting a relation between the structure and dynamics of polyelectrolyte solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dual-Purpose Canister Filling Demonstration Project Progress Report at ORNL, 2023

The US DOE Office of Nuclear Energy is investigating the feasibility of direct disposal of dual-purpose canisters (DPCs) in a hypothetical geological repository to offset the potential requirement to repackage spent nuclear fuel (SNF) from existing DPCs into smaller, disposal-ready canisters. Oak Ridge National Laboratory (ORNL) is currently evaluating the feasibility of filling void space in loaded DPCs with an engineered material to prevent a criticality event caused by groundwater/moderator intrusion. Metal alloys are being investigated as a filler material because of their relatively low viscosities when molten, which may facilitate their injection via an existing drainpipe that runs almost the full length of the DPC. ORNL’s strategy for evaluating filler viability includes simulations and physical demonstrations of filling and casting behavior, as well as evaluations of materials for compatibility in the repository environment. Filling of DPCs in this manner is expected to mitigate the risk associated with a post-closure criticality event during the repository performance assessment time frame (10,000 years or greater). Efforts in this fiscal year focused on (1) destructive analysis of experimental filler castings made in FY 2022, (2) a report outlining a conceptual design of a DPC filling facility (Fortner et al., 2023, M3SF 23OR010305044/ ORNL/SPR-2023/2921, May 31, 2023), (3) a report on affected features, events, and processes (FEPs) due to DPC filling (Price et al., 2023, M3SF-23SN010305093, issuance pending), (4) developing and testing a more practical alloy filler based upon a Sn-Al eutectic, and (5) modeling the heating and cooling dynamics of a DPC subjected to molten metal filling. The preliminary results from each of these tasks support the feasibility of filling DPCs with metal as a strategy against the possibility of criticality in the repository. The FY 2022 casting was found to penetrate even very small orifices in the mold and internal structures. Preliminary testing of the Sn-Al eutectic indicate little interaction of the melt with Zircaloy cladding. Thermal modelling shows that a filled DPC will cool to manageable temperatures within 2-3 days, which is likely manageable in a facility.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Symplectic Particle Tracking in a Thick Nonlinear McMillan Lens for the Fermilab Integrable Optics Test Accelerator (IOTA)

The McMillan system is a novel method to increase the tune spread of a beam without decreasing its dynamic aperture due to the systems integrability. While the ideal system is based on an infinitely thin kick, the physical design requires a thick electron lens, including a solenoid. Particle transport through the lens is difficult to simulate due to the nature of the force on the circulating beam. This paper demonstrates accurate simulation of a thick McMillan lens in a solenoid using symplectic integrators derived from Yoshidas method.

43 PARTICLE ACCELERATORS↗