Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic Programming”

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 1,567 records · Page 87

Power requirements and environmental impact of a pedelec. A case study based on real-life applications

This paper describes the methodologies to appraise the power requests and environmental analysis of an electrically assisted bicycle under real driving conditions, also containing regulations and technical-science-related aspects. For this purpose, in this study, the on-road test program of an electrically assisted bicycle was executed in the urban area of Naples on different test tracks, so a general assessment about its driving behavior under several driving conditions was performed. The power requirements in different typical riding situations were estimated by a procedure based on the experimental kinematic parameters that characterize the driving dynamics collected during the real-life applications. An environmental analysis was also performed, with a methodology that takes into account the environmental assessment of a moped by measuring the experimental moped exhaust emissions of the regulated pollutants. Starting from the results acquired during the different test samples, besides, an assessment of the electric traction offered by this pedelec on the driving comfort was evaluated for different riding situations. - Highlights: • The power requirements of a pedelec in typical riding conditions were identified. • The estimated electricity consumption for battery recharging was defined. • An environmental valuation of the tested pedelec and of a moped was performed. • Emissions that could be saved utilizing a pedelec instead of a moped were derived.

54 ENVIRONMENTAL SCIENCES↗

Empower Wall: Active insulation system leveraging additive manufacturing and model predictive control

Buildings are one of the largest energy consumers worldwide, using large amounts of energy during their construction and for climate control during operation. Active insulation systems (AIS) have been shown to reduce the energy needed for climate control in buildings by dynamically regulating the heat transferred between a building’s interior and exterior. Infrastructure-scale additive manufacturing (AM) has the potential to reduce the resources needed for building construction. Combining these two technologies into a single building envelope would create a path towards more sustainable buildings. A test was conducted for the Federal Energy Management Program (FEMP) Energy Exchange training and trade show, in August 2021, to investigate a new building envelope design, termed the Empower Wall, that utilized an AIS and was constructed using AM. Model predictive control was implemented to manage operation of the Empower Wall in concert with the existing HVAC system. Finally, the prototype system demonstrated that the Empower Wall lowered total energy consumption and reduced the cost of energy used.

36 MATERIALS SCIENCE↗

Singular linear-quadratic control problem for systems with linear delay

A singular linear-quadratic optimization problem on the trajectories of non-autonomous linear differential equations with linear delay is considered. The peculiarity of this problem is the fact that this problem has no solution in the class of integrable controls. To ensure the existence of solutions is required to expand the class of controls including controls with impulse components. Dynamical systems with linear delay are used to describe the motion of pantograph from the current collector with electric traction, biology, etc. It should be noted that for practical problems fact singularity criterion of quality is quite commonly occurring, and therefore the study of these problems is surely important. For the problem under discussion optimal programming control contained impulse components at the initial and final moments of time is constructed under certain assumptions on the functional and the right side of the control system.

97 MATHEMATICS AND COMPUTING↗

ESnet/JLab FPGA Accelerated Transport

To increase the science rate for high data rates/volumes, Thomas Jefferson National Accelerator Facility (JLab) has partnered with Energy Sciences Network (ESnet) to define an edge to data center traffic shaping / steering transport capability featuring data event aware network shaping and forwarding. The keystone of this ESnet+JLab FPGA Accelerated Transport (EJFAT) is the joint development of an AI/ML directed dynamic compute work Load Balancer (LB) of UDP streamed data. The LB is a suite consisting of a Field Programmable Gate Array (FPGA) executing the dynamically configurable, low fixed latency LB data plane featuring real-time packet redirection and high throughput, and a control plane running on the FPGA host computer that monitors network and compute farm telemetry in order to make dynamic AI/ML guided decisions for destination compute host redirection/load balancing and destination resource provisioning. The LB provides for three-tier horizontal scaling across LB suites, core compute hosts, and CPUs within a host. The LB effectively provides seamless integration of edge/core computing to support direct experimental data processing for immediate use by JLab science programs and others such as the EIC as well as data centers of the future requiring high throughput and low latency for both hot and cooled data for both running experiment data acquisition systems and data center use cases.

97 MATHEMATICS AND COMPUTING↗

Modeling transmission line constraints in two-stage robust unit commitment problem

Integration of renewable energy sources and demand response poses new challenges to system operators as they increase the uncertainty of the power supply and demand. Recently, robust optimization techniques are applied to the unit commitment problem with uncertainty as an alternative to the stochastic programming approaches. However, it remains challenging to solve the robust unit commitment model with full transmission line constraints. In this paper, we propose novel acceleration techniques for solving two-stage robust unit commitment problem with consideration of full transmission line constraints. We use 1) the cutting-plane algorithm for the master problem, which dynamically includes critical transmission line constraints, and 2) column-generation methods, including the branch-and-price-and-cut algorithm and heuristic approaches, for the subproblems, which add only necessary transmission line dual variables on the fly. Furthermore, computational results for the modified IEEE 118-bus system show that the combination of the cutting-plane algorithm and the heuristic column-generation approach greatly reduces the total solution time of the two-stage robust unit commitment problem.

Heuristic algorithms↗

CICE, The Los Alamos Sea Ice Model

The Los Alamos sea ice model (CICE) is the result of an effort to develop a computationally efficient sea ice component for a fully coupled atmosphere–land–ocean–ice global climate model. It was originally designed to be compatible with the Parallel Ocean Program (POP), an ocean circulation model developed at Los Alamos National Laboratory for use on massively parallel computers. CICE has several interacting components: a vertical thermodynamic model that computes local growth rates of snow and ice due to vertical conductive, radiative and turbulent fluxes, along with snowfall; an elastic-viscous-plastic model of ice dynamics, which predicts the velocity field of the ice pack based on a model of the material strength of the ice; an incremental remapping transport model that describes horizontal advection of the areal concentration, ice and snow volume and other state variables; and a ridging parameterization that transfers ice among thickness categories based on energetic balances and rates of strain. It also includes a biogeochemical model that describes evolution of the ice ecosystem. The CICE sea ice model is used for climate research as one component of complex global earth system models that include atmosphere, land, ocean and biogeochemistry components. It is also used for operational sea ice forecasting in the polar regions and in numerical weather prediction models.

Hunke, Elizabeth↗

Multiscale Modeling of Deformation, Fracture and Failure of Fusion Materials and Structures Final Report

The main objective of the research performed at UCLA under grant number DE-FG02-03ER54708 was to develop a physical understanding of the underlying mechanisms that control high-temperature deformation and fracture of fusion structural materials. The key aspect is the synergistic interaction between thermo-mechanical damage phenomena (e.g. creep, fatigue, and fracture), and irradiation damage phenomena (helium and hydrogen bubbles, and irradiation-induced defects). Another thrust area of the research was on the development of constitutive equations for the mechanical behavior of irradiated structures in fusion components (mainly the First Wall and Plasma-Facing components), and in utilizing them in descriptions in large-scale Finite Element Modeling of fusion structures. A parallel experimental effort on thermo-mechanical and plasma-induced damage is designed to validate and confront the developed theoretical and computational predictions. The models have been experimentally validated through experimental tests in simulated plasma transients on specially-designed W samples. This grant started at UCLA on November 14, 2003 and has had successive continuations till November 14, 2017. Research performed under the grant was to develop a range of hierarchical models for the post-elastic deformation, fracture and failure of fusion structural materials. The developed multiscale modeling approach has been based on rigorous mathematical, physical, and computational methods at the forefront of computational materials science. At the fundamental and smallest length scale (nm-μm), we developed advanced rate theory and Monte Carlo approaches to model microstructure evolution, non-equilibrium phase transformations, and dislocation-defect interactions. Microscopic and mesoscopic models of radiation hardening, ductile-to-brittle-transition, post-elastic deformation, plastic instabilities, and fracture processes are based on Dislocation Dynamics (DD) and Grain Growth Dynamics (GGD). At the continuum level, we focused on modeling fracture and failure mechanisms for Virtual Integrated Testing (VISTA), and for understanding the limits of both materials-by-design versus structural component design through advancing the new approach of multiscale-multiphysics. The overall research has been strongly coupled with the U.S. national experimental program to ascertain and verify the range of investigated phenomena. In the following, we give a summary of research progress supported by this grant since its inception till November 14, 2017. Detailed research results are available in our publications, listed in the publications section of this report.

36 MATERIALS SCIENCE↗

Relativistically Induced Transparency in Plasma (Final Research Report)

The major goal of this project was to study the process of relativistic transparency in plasmas. This process occurs when a very high intensity laser pulse interacts with a plasma. The electric field of the laser oscillates the electron speeds approaching the speed of light, increasing the effective mass of the electrons by the time-averaged Lorentz factor <γ>. The increase in the effective mass alters the critical density, the density at which the plasma becomes opaque to the electromagnetic wave of the laser. This project used experiments and particle-in-cell simulations to study this process in ultra-thin foil interactions. Experiments were performed through the LaserNetUS program on the Ohio State University laser system Scarlet. The transmitted and reflected light was studied as a function of the target thickness, the near-field profiles, the total energy and pulse duration were measured. Particle-in-cell simulations, using the OSIRIS code were used to model the interaction and better understand the dynamics. The results of both show current relativistic transparency theory is insufficient to predict the results, this is likely because the theories do not account for the laser energy absorption into the target and global plasma fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrating Machine-learning-assisted Computer Vision with RICH System

Developments in artificial intelligence have vastly expanded the capabilities of robots. Currently, the Spallation Neutron Source (SNS) beamlines at Oak Ridge National Lab (ORNL) have robotic sample loaders to increase the efficiency of running experiments. However, they require retraining if anything about the situation changes, e.g., where the samples are, and cannot notice if errors occur. So, the viability of using computer vision and machine learning to enhance these sample loaders’ functionality was investigated. In this project, the RICH system with a Dobot CR3 6-axis robot present at the VULCAN beamline assisted by an Intel Realsense D435i camera, a unique camera that enables convenient translation of 2D pixel coordinates to 3D world points, was programmed to load ceramic crucibles into a thermogravimetric analyzer (TGA) furnace. An algorithm was constructed in Python with three major phases planned: (1) obtaining a sample, (2) moving it to the target location, and then (3) bringing the sample back to its original location once the experiment finished. In the first phase, the algorithm would dynamically detect sample locations using ArUco markers to recognize the samples’ general location and a custom-trained yolov5 object detection model to locate the crucibles’ centers. Afterward, the robot would be directed to pick up samples based on the crucibles’ calculated positions. In the second phase, the robot would move the sample to a secondary point, reorient its grip, and place the sample at the target location. In the final phase, the robot would determine whether the sample was intact and would bring it back to its original place if it was or raise an alarm. Using this algorithm, the robot was able to pick up different types of crucibles at varying positions. These results indicate that integrating machine-learning-assisted computer vision with robotic sample loaders can result in effective autonomous detection of samples.

97 MATHEMATICS AND COMPUTING↗

Fingerprints of a riming event on cloud radar Doppler spectra: observations and modeling

In this paper, Radar Doppler spectra measurements are exploited to study a riming event when precipitating ice from a seeder cloud sediment through a supercooled liquid water (SLW) layer. The focus is on the "golden sample" case study for this type of analysis based on observations collected during the deployment of the Atmospheric Radiation Measurement Program's (ARM) mobile facility AMF2 at Hyytiälä, Finland, during the Biogenic Aerosols – Effects on Clouds and Climate (BAECC) field campaign. The presented analysis of the height evolution of the radar Doppler spectra is a state-of-the-art retrieval with profiling cloud radars in SLW layers beyond the traditional use of spectral moments. Dynamical effects are considered by following the particle population evolution along slanted tracks that are caused by horizontal advection of the cloud under wind shear conditions. In the SLW layer, the identified liquid peak is used as an air motion tracer to correct the Doppler spectra for vertical air motion and the ice peak is used to study the radar profiles of rimed particles. A 1-D steady-state bin microphysical model is constrained using the SLW and air motion profiles and cloud top radar observations. The observed radar moment profiles of the rimed snow can be simulated reasonably well by the model, but not without making several assumptions about the ice particle concentration and the relative role of deposition and aggregation. In conclusion, this suggests that in situ observations of key ice properties are needed to complement the profiling radar observations before process-oriented studies can effectively evaluate ice microphysical parameterizations.

54 ENVIRONMENTAL SCIENCES↗

CICE, The Los Alamos Sea Ice Model

The Los Alamos sea ice model (CICE) is the result of an effort to develop a computationally efficient sea ice component for a fully coupled atmosphere–land–ocean–ice global climate model. It was originally designed to be compatible with the Parallel Ocean Program (POP), an ocean circulation model developed at Los Alamos National Laboratory for use on massively parallel computers. CICE has several interacting components: a vertical thermodynamic model that computes local growth rates of snow and ice due to vertical conductive, radiative and turbulent fluxes, along with snowfall; an elastic-viscous-plastic model of ice dynamics, which predicts the velocity field of the ice pack based on a model of the material strength of the ice; an incremental remapping transport model that describes horizontal advection of the areal concentration, ice and snow volume and other state variables; and a ridging parameterization that transfers ice among thickness categories based on energetic balances and rates of strain. It also includes a biogeochemical model that describes evolution of the ice ecosystem. The CICE sea ice model is used for climate research as one component of complex global earth system models that include atmosphere, land, ocean and biogeochemistry components. It is also used for operational sea ice forecasting in the polar regions and in numerical weather prediction models.

Hunke, Elizabeth↗

The Alliance of Advanced Process Control and Accountability – A Future Safeguards-By-Design Tool

For any chemical separation process producing a valuable product, a material balance is an important process control measurement. That is particularly true for the separation of actinides from irradiated nuclear fuel, not only for their intrinsic value but also because an incomplete material balance may indicate diversion for unauthorized use. The DOE Office of Nuclear Energy is currently carrying out at the Pacific Northwest National Laboratory an experimental measurement of how well and with what precision current technologies can implement near real-time actinide material balances. This measurement effort is called the CoDCon project. It involves the separation of a product with a 70/30 uranium/plutonium mass ratio. Initial tests will use dissolved fuel simulants prepared with pure uranium and plutonium nitrates at the same input ratios as irradiated fuel. Subsequent testing with actual irradiated fuel would be done to verify the results obtained with simulants. The experiments will use advanced on-line instrumentation supported by dynamic process models. Since accountability uncertainties could mask diversions, the aim of the project is not only to measure present-day capabilities but also, through sensitivity analyses, to identify those measurements with the greatest potential for overall material-balance improvements. The latter results will help identify priorities for future fuel cycle R&D programs. Advanced separations process control and material accountability technologies thus have a common goal: to provide the best tools available for safeguards-by-design [defined by the International Atomic Energy Agency (IAEA) as the integration of the design of a new nuclear facility through planning, construction, operation and decommissioning]. Since the potential domestic use of CoDCon results may be later than their possible foreign applications, arrangements may be feasible for possible bilateral or multinational cooperation in the CoDCon project.

Lumetta, Gregg J.↗

Volume III. DUNE far detector technical coordination

The preponderance of matter over antimatter in the early universe, the dynamics of the supernovae that produced the heavy elements necessary for life, and whether protons eventually decay -- these mysteries at the forefront of particle physics and astrophysics are key to understanding the early evolution of our universe, its current state, and its eventual fate. The Deep Underground Neutrino Experiment (DUNE) is an international world-class experiment dedicated to addressing these questions as it searches for leptonic charge-parity symmetry violation, stands ready to capture supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. The DUNE far detector technical design report (TDR) describes the DUNE physics program and the technical designs of the single- and dual-phase DUNE liquid argon TPC far detector modules. Volume III of this TDR describes how the activities required to design, construct, fabricate, install, and commission the DUNE far detector modules are organized and managed. This volume details the organizational structures that will carry out and/or oversee the planned far detector activities safely, successfully, on time, and on budget. It presents overviews of the facilities, supporting infrastructure, and detectors for context, and it outlines the project-related functions and methodologies used by the DUNE technical coordination organization, focusing on the areas of integration engineering, technical reviews, quality assurance and control, and safety oversight. Because of its more advanced stage of development, functional examples presented in this volume focus primarily on the single-phase (SP) detector module.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Detection and Perception of Sound by Eagles and Surrogate Raptors

One overarching objective of this program of study was the accumulation of objective, scientifically valid information relating to auditory performance of bald and golden eagles that may be used to guide the development of acoustic alerting/deterrence technologies intended to discourage encroachment into wind energy air spaces. To that end, analyses aimed at the characterization of sensitivity to sound in bald and golden eagles, along with findings in the supra-threshold, dynamic frequency spaces related to response latencies and amplitudes, leads us to conclude that bald, and golden eagles navigate the same basic working auditory space, as in other known and thus far characterized members of the diurnal raptor family. Specifically, bald and golden eagles, along with other raptor species within the group, operate in an auditory space characterized by a frequency band at least four octaves wide and centered on 2 kHz, with an upper frequency limit between 6 and 10 kHz at 80 dB SPL and a lower frequency limit that almost certainly extends below 0.2 kHz. Consequently, we recommend that signal designers use these data as a guideline in efforts to design effective and efficient acoustic alerting/deterrent systems. It is important to note that signal energy broadcast outside of this frequency band at moderate levels will not contribute to the efficacy of a deterrent but will add an unnecessary fraction to the overall acoustic pollution budget. The importance of this consideration is heightened by contemporaneous concerns related to the transmission of noise broadcast by wind energy farms. In addition, based on analyses of data acquired from red-tailed hawks using the same experimental paradigm and data acquisition system, we conclude that auditory function in the red-tailed hawk is sufficiently like that observed in bald and golden eagles to permit its use as a surrogate species. Response waveforms, threshold-frequency curves, and input-output characteristics match those of eagles closely. It should be noted however, that differences in sensitivity and slightly extended high-frequency limits of hearing should be taken into account when extrapolating findings from one species to the others. Although the inclusion of behavioral tests of red-tailed hawks to acoustic stimuli was beyond the scope of this investigation, future efforts to assess response parameters like signal-type preference and habituation rate will further elucidate their suitability to serve as eagle surrogates in behavioral studies; nonetheless, the species in question are well matched with respect to basic auditory performance. A second essential objective of this program of study was the acoustic characterization of a subset of calls comprising the vocal repertoires of bald and golden eagles that may be used to supplement auditory performance findings in the effort to guide the development of acoustic alerting signals. With regard to that objective, the vocal repertoires of both bald and golden eagle species are rich and varied. While similar in spectrographic structure, distinctive differences are also clear. Generally, golden eagles produce some calls with shorter durations, and similar “sounding” calls exhibit distinctively different spectrographic patterns than those of bald eagles. Both species produce calls that contain a wide variety of nonlinear elements that operate to enhance the rich and varied nature of commonly observed vocal products. Comparison of the average power spectra of commonly observed bald and golden eagle calls with threshold-frequency curves leads to the conclusion that call energies fall within the frequency bounds of hearing. Further, the acoustic energy of calls considered in this report tend to fall into overlapping, but different frequency ranges of the acoustic sensitivity curve. This condition may encourage signal designers to vary the frequency content of acoustic deterrence signals in the field. Finally, preliminary observations relating to the tendencies and proclivities of bald eagles to attend to the acoustic landscape lead to the conclusion that eagles monitor their immediate sound environment assiduously. Individuals respond to a variety of natural and synthetic sound signals reliably and, perhaps most relevant in the context of the engineering of acoustic alerting/deterrence technologies, habituation to most sounds considered in this effort was minimal. These preliminary results, while calling for extended behavioral testing, are promising and set the stage for the exportation of behavioral studies into real world scenarios.

17 WIND ENERGY↗

Defect Production and Microstructural Feature Impact for Radiation Damage in Additively Manufactured 316 Stainless Steel

This milestone presents multi-scale modeling research results for additively manufactured 316 stainless steel. A combination of phase field, cluster dynamics, molecular dynamics, and density functional theory with machine learning is used, allowing for predictions of radiation-driven microstructural evolution in additively manufactured 316 stainless steel over a range of temperatures, damage rates, neutron spectra, and microstructures, and supporting the development of combined ion and neutron irradiation for material qualification. Informed by ion irradiation and neutron irradiation results across the Advanced Materials and Manufacturing Technologies program, we investigate the unique aspects of radiation-driven microstructure evolution in additively manufactured 316 stainless steel. In particular, we focus on understanding the impact of carbon concentration (varying, for example, between the 316L and 316H standards) on void formation; radiation-induced segregation at dislocation cells and grain boundaries; and the evolution of dislocation loops and network dislocation populations. We find that the unique characteristics of the additively manufactured microstructures must be accounted for in understanding the evolution of dislocation populations under thermal and irradiation conditions, such as the variation in sink strengths arising due to the variation in dislocation density. We also find that the radiation-induced segregation of Cr and Ni to grain boundaries and cell walls differs due to the differences in their defect sink biases. We also find that increasing the Ni content can slow vacancy diffusion, which may provide a mechanism for the observed reduction in transient swelling rate for austenitic Fe-Cr-Ni alloys with increasing Ni content. In addition, ion irradiations have shown that increasing carbon content in 316 SS results in a larger population of smaller voids, suggesting reduced vacancy diffusion. Our results show that the carbon content of additively manufactured 316 SS has a significant impact on the migration rate of defect clusters. The presence of carbon atoms results in carbon-vacancy trapping, significantly reducing the diffusion rate of vacancies. Carbon atoms may also be trapped near the surface of a void, which may reduce void growth by trapping vacancies that diffuse toward the void.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SU-F-T-563: Delivered Dose Reconstruction of Moving Targets for Gated Volumetric Modulated Arc Therapy (VMAT)

Purpose: Actual delivered dose of moving tumors treated with gated volumetric arc therapy (VMAT) may significantly differ from the planned dose assuming static target. In this study, we developed a method which reconstructs actual delivered dose distribution of moving target by taking into account both tumor motion and dynamic beam delivery of gated VMAT, and applied to abdominal tumors. Methods: Fifteen dual-arc VMAT plans (Eclipse, Varian Medical Systems) for 5 lung, 5 pancreatic, and 5 liver cancer patients treated with gated VMAT stereotactic body radiotherapy (SBRT) were studied. For reconstruction of the delivered dose distribution, we divided each original arc beam into control-point-wise sub-beams, and applied beam isocenter shifting to each sub-beam to reflect the tumor motion. The tumor positions as a function of beam delivery were estimated by synchronizing the beam delivery with the respiratory signal which acquired during treatment. For this purpose, an in-house program (MATLAB, Mathworks) was developed to convert the original DICOM plan data into motion-involved treatment plan. The motion-involved DICOM plan was imported into Eclipse for dose calculation. The reconstructed delivered dose was compared to the plan dose using the dose coverage of gross tumor volume (GTV) and dose distribution of organs at risk (OAR). Results: The mean GTV dose coverage difference between the reconstructed delivered dose and the plan dose was 0.2 % in lung and pancreas cases, and no difference in liver cases. Mean D1000cc of ipsilateral lungs was reduced (0.8 ± 1.4cGy). Conclusion: We successfully developed a method of delivered dose reconstruction taking into account both respiratory tumor motion and dynamic beam delivery, and applied it to abdominal tumors treated with gated VAMT. No significant deterioration of delivered dose distribution indicates that interplay effect would be minimal even in the case of gated SBRT. This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (No. 2015038710)

60 APPLIED LIFE SCIENCES↗

Hydrocode Development on the Ncube and the Connection Machine Hypercubes [Book Chapter]

Hydrocode simulations constitute an important tool at Sandia National Laboratories and elsewhere for analyzing complex two- and three-dimensional systems. However, current vector supercomputers do not provide a growth path to enable fast, routine, and cost-effective simulations of large problems. Future, massively-parallel computers will provide a solution. Sandia has already developed simplified versions of the production hydrocode CTH for the Connection Machine and the nCUBE massively-parallel supercomputers. The parallel versions solve problems in two-dimensional, multi-fluid, shock-wave physics. Code development strategy, coding methodology, visualization techniques and performance results for this work are described.

640410 -- Fluid Physics-- General Fluid Dynamics↗

JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

Efficient accelerator modeling and particle tracking are key for the design and configuration of modern particle accelerators. In this work, we present JuTrack, a nested accelerator modeling package developed in the Julia programming language and enhanced with compiler-level automatic differentiation (AD). With the aid of AD, JuTrack enables rapid derivative calculations in accelerator modeling, facilitating sensitivity analyses and optimization tasks. Here we demonstrate the effectiveness of AD-derived derivatives through several practical applications, including sensitivity analysis of space-charge-induced emittance growth, nonlinear beam dynamics analysis for a synchrotron light source, and lattice parameter tuning of the future Electron-Ion Collider (EIC). Through the incorporation of automatic differentiation, this package opens up new possibilities for accelerator physicists in beam physics studies and accelerator design optimization.

43 PARTICLE ACCELERATORS↗