Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

PETSc Users Manual (Rev. 3.5)

This manual describes the use of PETSc for the numerical solution of partial differential equations and related problems on high-performance computers. The Portable, Extensible Toolkit for Scientific Computation (PETSc) is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all message-passing communication. PETSc includes an expanding suite of parallel linear, nonlinear equation solvers and time integrators that may be used in application codes written in Fortran, C, C++, Python, and MATLAB (sequential). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++ or Fortran and experience using a debugger such as gdb or dbx, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. ;For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. PETSc also has a MATLAB interface, so portions of your code can be written in MATLAB to “try out” the PETSc solvers. The resulting code will not be scalable however because currently MATLAB is inherently not scalable; and PETSc should not be used to attempt to provide a “parallel linear solver” in an otherwise sequential code. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect any kind of reasonable performance. Do not expect to generate your matrix sequentially and then “use PETSc” to solve the linear system in parallel. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc is supported; see the web site http://www.mcs.anl.gov/petsc for information on contacting support. A http://www.mcs.anl.gov/petsc/publications may be found a list of publications and web sites that feature work involving PETSc. We welcome any reports of corrections for this document.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A split-step method to include electron–electron collisions via Monte Carlo in multiple rate equation simulations

A split-step numerical method for calculating ultrafast free-electron dynamics in dielectrics is introduced. The two split steps, independently programmed in C++11 and FORTRAN 2003, are interfaced via the presented open source wrapper. The first step solves a deterministic extended multi-rate equation for the ionization, electron–phonon collisions, and single photon absorption by free-carriers. The second step is stochastic and models electron–electron collisions using Monte-Carlo techniques. This combination of deterministic and stochastic approaches is a unique and efficient method of calculating the nonlinear dynamics of 3D materials exposed to high intensity ultrashort pulses. Results from simulations solving the proposed model demonstrate how electron–electron scattering relaxes the non-equilibrium electron distribution on the femtosecond time scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

An Optimization-Based Planning Tool for On-Demand Mobility Service Operations

Regions worldwide are adopting and exploring low-speed automated electric shuttle (AES) service as an on-demand shared mobility service in dense geofenced urban areas. Building on this concept, the National Renewable Energy Laboratory (NREL) recently developed the Automated Mobility District (AMD) toolkit. The AMD toolkit—comprising of a travel micro-simulation model and an energy estimation model—estimates the mobility and energy impacts of a given shuttle configuration within an AMD. Early-stage AMD deployments need to find optimal operational configurations that include: (a) passenger capacity of an AES, (b) time-dependent routes, and (c) fleet size (AES units) to satisfy the demand for the region. This research extends the AMD toolkit functionality by developing an optimization-based planning module that will assist in the operations of AES units. We developed a constrained mixed-integer program accounting for passenger waiting time, battery range, and passenger capacity of AES units. For scalability, we demonstrated the Tabu search-based solution technique for a real-world network—a proposed AMD deployment in Greenville, South Carolina, USA. Compared to rule-based operations, our developed solution yields higher travel time and energy savings for the network at different demand levels. The sensitivity analyses for waiting time thresholds indicate nonlinearity in the system performance, underscoring the need to meet shared-use mobility user-level expectations. The developed optimization framework can be adapted and extended to accommodate different categories of shared-use on-demand mobility services.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗
DOE OSTI2019

Additive Manufacturing of Biomechanically Tailored Meshes for Compliant Wearable and Implantable Devices

Abstract Additive manufacturing (AM) of medical devices such as orthopedic implants and hearing aids is highly attractive because of the potential of AM to match the complex form and mechanics of individual human bodies. Externally worn and implantable tissue‐support devices, such as ankle or knee braces, and hernia repair mesh, offer a new opportunity for AM to mimic tissue‐like mechanics and improve both patient outcomes and comfort. Here, it is demonstrated how explicit programming of the toolpath in an extrusion AM process can enable new, flexible mesh materials having digitally tailored mechanical properties and geometry. Meshes are fabricated by extrusion of thermoplastics, optionally with continuous fiber reinforcement, using a continuous toolpath that tailors the elasticity of unit cells of the mesh via incorporation of slack and modulation of filament–filament bonding. It is shown how the tensile mesh mechanics can be engineered to match the nonlinear response of muscle. An ankle brace with directionally specific inversion stiffness arising from embedded mesh is validated, and further concepts for 3D mesh devices are prototyped.

Pattinson, Sebastian W.↗
DOE OSTI2015

Method for conducting nonlinear electrochemical impedance spectroscopy

A method for conducting nonlinear electrochemical impedance spectroscopy. The method includes quantifying the nonlinear response of an electrochemical system by measuring higher-order current or voltage harmonics generated by moderate-amplitude sinusoidal current or voltage perturbations. The method involves acquisition of the response signal followed by time apodization and fast Fourier transformation of the data into the frequency domain, where the magnitude and phase of each harmonic signal can be readily quantified. The method can be implemented on a computer as a software program.

25 ENERGY STORAGE↗
DOE OSTI2019

Designing Metabolic Division of Labor in Microbial Communities

ABSTRACT Microbes face a trade-off between being metabolically independent and relying on neighboring organisms for the supply of some essential metabolites. This balance of conflicting strategies affects microbial community structure and dynamics, with important implications for microbiome research and synthetic ecology. A “gedanken” (thought) experiment to investigate this trade-off would involve monitoring the rise of mutual dependence as the number of metabolic reactions allowed in an organism is increasingly constrained. The expectation is that below a certain number of reactions, no individual organism would be able to grow in isolation and cross-feeding partnerships and division of labor would emerge. We implemented this idealized experiment using in silico genome-scale models. In particular, we used mixed-integer linear programming to identify trade-off solutions in communities of Escherichia coli strains. The strategies that we found revealed a large space of opportunities in nuanced and nonintuitive metabolic division of labor, including, for example, splitting the tricarboxylic acid (TCA) cycle into two separate halves. The systematic computation of possible solutions in division of labor for 1-, 2-, and 3-strain consortia resulted in a rich and complex landscape. This landscape displayed a nonlinear boundary, indicating that the loss of an intracellular reaction was not necessarily compensated for by a single imported metabolite. Different regions in this landscape were associated with specific solutions and patterns of exchanged metabolites. Our approach also predicts the existence of regions in this landscape where independent bacteria are viable but are outcompeted by cross-feeding pairs, providing a possible incentive for the rise of division of labor. IMPORTANCE Understanding how microbes assemble into communities is a fundamental open issue in biology, relevant to human health, metabolic engineering, and environmental sustainability. A possible mechanism for interactions of microbes is through cross-feeding, i.e., the exchange of small molecules. These metabolic exchanges may allow different microbes to specialize in distinct tasks and evolve division of labor. To systematically explore the space of possible strategies for division of labor, we applied advanced optimization algorithms to computational models of cellular metabolism. Specifically, we searched for communities able to survive under constraints (such as a limited number of reactions) that would not be sustainable by individual species. We found that predicted consortia partition metabolic pathways in ways that would be difficult to identify manually, possibly providing a competitive advantage over individual organisms. In addition to helping understand diversity in natural microbial communities, our approach could assist in the design of synthetic consortia.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗
DOE OSTI2014

Understanding the many-body expansion for large systems. I. Precision considerations

Electronic structure methods based on low-order “n-body” expansions are an increasingly popular means to defeat the highly nonlinear scaling of ab initio quantum chemistry calculations, taking advantage of the inherently distributable nature of the numerous subsystem calculations. Here, we examine how the finite precision of these subsystem calculations manifests in applications to large systems, in this case, a sequence of water clusters ranging in size up to (H{sub 2}O){sub 47}. Using two different computer implementations of the n-body expansion, one fully integrated into a quantum chemistry program and the other written as a separate driver routine for the same program, we examine the reproducibility of total binding energies as a function of cluster size. The combinatorial nature of the n-body expansion amplifies subtle differences between the two implementations, especially for n ⩾ 4, leading to total energies that differ by as much as several kcal/mol between two implementations of what is ostensibly the same method. This behavior can be understood based on a propagation-of-errors analysis applied to a closed-form expression for the n-body expansion, which is derived here for the first time. Discrepancies between the two implementations arise primarily from the Coulomb self-energy correction that is required when electrostatic embedding charges are implemented by means of an external driver program. For reliable results in large systems, our analysis suggests that script- or driver-based implementations should read binary output files from an electronic structure program, in full double precision, or better yet be fully integrated in a way that avoids the need to compute the aforementioned self-energy. Moreover, four-body and higher-order expansions may be too sensitive to numerical thresholds to be of practical use in large systems.

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

Analysis of SNL/MSU/DOE Fatigue Database Trends for Wind Turbine Blade Materials 2010-2015.

Wind turbine blades are designed to several major structural conditions, including tip deflection, strength and b uckling during severe loading, as well as very high numbers of fatigue cycles and various service environments. The MSU Database Program has, since 1989, addressed the broad range of properties needed for current and potential blade materials through stati c and fatigue testing and test development in cooperation with Sandia National Laboratories and wind industry and supplier partners. This report is the latest in a series, giving test results and analysis for the period 2010 - 2015. Program data are compiled in a public database [1] and other reports and publications given in the cited references. The report begins with an executive summary and introductory material including background discussion of previous related studies. Section 3 describes experimental methods including processing, test methods, instrumentation and test development. Section 4 provides static tension, compression and shear stress - strain properties in three directions using coupons sectioned from a thick infused unidirectional glass/epoxy laminate. The nonlinear, shear dominated static properties were characterized with loading - u nloading - reloading (LUR) tests in tension and compression to increasing load levels, for +-45O laminates. Section 5 explores the origins of tensile fatigue sensitivity in glass fiber dominated laminates with variations in fabric architecture including speci ally prepared fabrics and aligned strand laminates. Several types of resins are considered, with variations in resin toughness and bonding to fibers, as well as cure cycle variations for an epoxy. Conclusions are drawn as to the limits of tensile fatigue r esistance and the effects of resin type and fabric architecture, including the behavior of a commercial aligned glass strand product. Interactions between cyclic fatigue response and creep are addressed for off - axis (+-45O) glass/epoxy laminates in Sectio n 6. The nonlinear fatigue and creep stress - strain and cumulative strain response are characterized in tension and compression as a function of stress level, cycles and cumulative time, using square and sinewave loading over a broad range of frequency. The results are analyzed in terms of the cycles and cumulative time under load. A cumulative strain failure criterion is established, and used to construct shear and tension constant life diagrams (CLD's) with data for nine R - values. The effects of a more duc tile urethne resin are also explored. A previous study of thick adhesives testing is extended to mixed mode fracture mechanics testing in Section 7. Mechanisms of static and fatigue crack extension near the laminate adherend interface are reported in deta il. Data are presented for mixed mode adhesive fracture, compared to mixed mode fracture in ply delamination. Fatigue crack growth exponents are also developed for a mixed mode cracked lap shear coupon. The data for fatigue trends and relative failure stra ins and exponents are compared for various blade component materials in Section 8. The effects of temperature and seawater saturation are considered for selected materials of interest for wind and hydrokinetic turbine blades in Section 9. Section 10 gives detailed conclusions for each section. A cknowledgements The research presented in this report was carried out under Sandia National Laboratories purchase orders 1325028 an d 1543945 between 2010 and 2015, with support from the DOE Wind and Water Technologies Office . In addition to the authors listed, significant contributions were made by Patrick Flaherty, Pancastya Agastra, Michael Schuster, and Michael Voth. Industry m aterials suppliers include Vectorply, Saertex, OCV, AGY, Bayer, Ashland, 3M and Nextel. Industry suppliers with significant contributions to the study were Hexion, PPG, Reichhold, Gurit and NEPTCO. Intentionally Left Blank

17 WIND ENERGY↗
DOE OSTI2012

Potential for luminosity improvement for low-energy RHIC operation

At the Brookhaven National Laboratory, a physics program, motivated by the search of the QCD phase transition critical point, requires operation of the Relativistic Heavy Ion Collider (RHIC) with heavy ions at very low beam energies corresponding to 2.5-20 GeV/n. Several physics runs were already successfully performed at these low energies. However, the luminosity is very low at lowest energies of interest (< 10 GeV/n) limited by the intra-beam scattering and space-charge, as well as by machine nonlinearities. At these low energies, electron cooling could be very effective in counteracting luminosity degradation due to the IBS, while it is less effective against other limitations. Overall potential luminosity improvement for low-energy RHIC operation from cooling is summarized for various energies, taking into account all these limitations as well as beam lifetime measured during the low-energy RHIC runs. We also explore a possibility of further luminosity improvement under the space-charge limitation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Advances in Multi-Physics Modeling for Parallel Finite Element Code Suite ACE3P

SLAC developed ACE3P code suite is a comprehensive set of parallel finite-element codes for multi-physics modeling of particle accelerators. Running on massively parallel computer platforms for high fidelity and high accuracy simulation, ACE3P enables rapid virtual prototyping of accelerator and RF component design, optimization and analysis. Furthermore, recent advances of ACE3P have been focused on the implementation of advanced numerical algorithms, enhancement of multi-physics capabilities, and improvement of code performance on state-of-the-art high performance computing (HPC) platforms for large-scale RF modeling in accelerator applications and beyond. In this paper, we will present a nonlinear electromagnetic (EM) eigensolver for damping calculation of resonant modes for cavities coupled with waveguides, the thermal solver for accurate evaluation of heat loads in superconducting RF (SRF) cavity power coupler, the elastic solver for investigation of SRF cavity cryomodule (CM) bowing due to thermal effect, the frequency domain mechanical solver for studying SRF cavity mechanical oscillation modes, the particle tracking code for dark current simulation in SRF CM using a hybrid MPI+OpenMP parallel programing model, the time domain solver for high fidelity modeling of EM wave propagation for medical device in human body.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗
DOE OSTI2021

Solving the Dynamics-Aware Economic Dispatch Problem with the Koopman Operator

The dynamics-aware economic dispatch (DED) problem embeds low-level generator dynamics and operational constraints to enable near real-time scheduling of generation units in a power network. DED produces a more dynamic supervisory control policy than traditional economic dispatch (T-ED) that reduces overall generation costs. However, in contrast to T-ED, DED is a nonlinear, non-convex optimization problem that is computationally prohibitive to solve. We introduce a machine learning-based operator-theoretic approach for solving the DED problem efficiently. Specifically, we develop a novel discrete-time Koopman Operator (KO) formulation that embeds domain information into the structure of the KO to learn high-fidelity approximations of the generator dynamics. Using the KO approximation, the DED problem can be reformulated as a computationally tractable linear program (abbreviated DED-KO). We demonstrate the high solution quality and computational-time savings of the DED-KO model over the original DED formulation on a 9-bus test system.

King, Ethan↗
DOE OSTI2019

Economic Optimal Control for Minimizing Fuel Consumption of Heavy-Duty Trucks in a Highway Environment

This paper provides a comparative assessment of three economic optimal control strategies, aimed at minimizing the fuel consumption of heavy-duty trucks in a highway environment, under a representative lead vehicle model informed by traffic data. These strategies fuse a global, offline dynamic programming (DP) optimization with online model predictive control (MPC). We then show how two of the three strategies can be adapted to accommodate the presence of traffic and optimally navigate signalized intersections using infrastructure-to-vehicular (I2V) communication. The MPC optimization, which is local in nature, makes refinements to a coarsely (but globally, subject to grid resolution) optimized target velocity profile from the DP optimization. The three candidate economic MPC formulations that are evaluated include: a nonlinear time-based formulation that directly penalizes predicted fuel consumption, a nonlinear time-based formulation that penalizes braking effort as a surrogate for fuel consumption, and a linear distance-based convex formulation that maintains a tradeoff between energy expenditure and tracking of the coarsely optimized velocity profile obtained from DP. Using a medium-fidelity Simulink model, based on a Volvo truck's longitudinal and engine dynamics, we analyze the optimization's performance on four highway routes under various traffic scenarios. Results demonstrate 3.7-8.3% fuel economy improvement on highway routes without traffic and 6.5-10% on the same routes with traffic included. Furthermore, we present a detailed analysis of energy usage by "type" (aerodynamic losses, braking losses, and comparison of brake-specific fuel consumption), under each candidate control strategy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

One- and Two-Qubit Gates: Rabi Technique and Single Unipolar Pulses

Methods for controlling states of interacting superconducting flux qubits using power-efficient devices of fast single-quantum logic (Josephson nonlinearity cavities) are studied. One- and two-qubit quantum logical operations performed within the conventional control technique using Rabi pulses and using picosecond single unipolar magnetic field pulses are comparatively analyzed. It is shown that all main operations can be implemented with an accuracy of better than 97% due to optimization of the shape and parameters of unipolar control pulses (associated with, e.g., propagation of fluxons in transmission lines). The efficiency of the developed technique for programming a two-qubit quantum processor implementing the simplest Deutsch–Jozsa algorithm is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Seismic Retrofit for Electric Power Systems

Our paper develops a two-stage stochastic program and solution procedure to optimize the selection of seismic retrofit strategies to increase the resilience of electric power systems against earthquake hazards. The model explicitly considers the range of earthquake events that are possible and, for each, an approximation of the distribution of damage experienced. Furthermore, this is important because electric power systems are spatially distributed and so their performance is driven by the distribution of component damage. We also test this solution procedure against the nonlinear integer solver in LINGO 13 and apply the formulation and solution strategy to the Eastern Interconnection, where seismic hazard stems from the New Madrid seismic zone.

58 GEOSCIENCES↗
DOE OSTI2011

CESR Conversion Damping Ring Studies of Electron Cloud Instabilities (CESR-TA)

In the International Linear Collider, two linear accelerators will accelerate bunches of positrons and electrons to over a hundred billion electron volts and collide them in a central detector. In order to obtain useful collision rates, the bunches, each containing twenty billion particles, must be compressed to a cross section of a few nanometers by a few hundred nanometers. In order to prepare these ultra high density bunches, damping rings (DRs) are employed before the linear accelerators. The DRs take the high emittance bunches that are provided by the electron and positron sources and, through the process of radiation damping, squeeze them into ultra low emittance beams that are ready for the main linear accelerators. In the damping rings, a number of effects can prevent the successful preparation of the beams. In the electron ring, an effect known as the fast ion instability can lead to beam growth and, in the positron ring, the build-up of an electron cloud (EC), which interacts with the circulating bunches, can produce the same effect. EC build-up and the subsequent interaction of the cloud with the positron beam in the DR have been identified as major risks for the successful construction of a linear collider. The CESRTA research program at the Cornell Electron Storage Ring (CESR) was developed in order to study the build-up of the EC, the details of its impact on ultra low emittance beams, as well as methods to mitigate the impact of the cloud. In the DR, the EC forms when synchrotron photons radiated from the circulating beam strike the walls of the vacuum chamber, resulting in the emission of photoelectrons. These low energy electrons can be accelerated across the vacuum chamber by the electric field of the beam, and strike the walls, causing the emission of secondary electrons. The secondary electrons are subsequently accelerated into the walls yet again via the same mechanism. The result is that the EC can rapidly begin to fill the vacuum chamber. In an electron DR, the EC build-up is limited by the Coulomb repulsion. But in a positron ring, the electrons are pulled into the potential well of the beam. The resulting interaction of the circulating bunches of positrons with the EC that ultimately limits DR performance. Typically we store long trains of closely space bunches in the damping ring. The interaction of stored beam and the EC that is generated by the long train is manifested by different mechanisms: (1) The cloud focuses the beam, which causes a tune shift that increases along the bunch train as the cloud density increases. (2) The cloud electrons couple the motion of bunches along a train. Transverse motion of a leading bunch is transferred to the cloud, and subsequently to a trailing bunch which can result in a multi-bunch instability. (3) The cloud couples the positrons in the head of the bunch to those in the tail of the same bunch, which can excite a 'head-tail' instability. (4) The nonlinear fields of the EC can lead to emittance growth before the onset of instabilities. The CESRTA collaboration, which includes researchers from Cornell University's Laboratory for Elementary-Particle Physics as well as more than 50 senior staff members from over a dozen accelerator laboratories and universities around the world, has operated CESR as a damping ring for the past three years to study these EC effects. A range of specialized instrumentation has been deployed to study the local build-up of the cloud in the vacuum chambers as well as the complicated dynamics exhibited when the beam and the EC interact. The program has significantly advanced our understanding of these issues and has helped identify the most promising methods to mitigate the impact of the EC on the DR beams. It has pointed the way towards a DR design that can meet the stringent specifications of the ILC.

43 PARTICLE ACCELERATORS↗