Parallel Computation of Intergrated Electronmagnetic, Thermal and Structural Effects for Accelerator Cavities

The successful operation of accelerator cavities has to satisfy both rf and mechanical requirements. It is highly desirable that electromagnetic, thermal and structural effects such as cavity wall heating and Lorentz force detuning in superconducting rf cavities can be addressed in an integrated analysis. Based on the SLAC parallel finite-element code infrastructure for electromagnetic modeling, a novel multi-physics analysis tool has been developed to include additional thermal and mechanical effects. The parallel computation enables virtual prototyping of accelerator cavities on computers, which would substantially reduce the cost and time of a design cycle. The multi-physics tool is applied to the LCLS rf gun for electromagnetic, thermal and structural analyses

Thermal Analysis of SRF Cavity Couplers Using Parallel Multi physics Tool TEM3P

SLAC has developed a multi-physics simulation code TEM3P for simulating integrated effects of electromagnetic, thermal and structural loads. TEM3P shares the same software infrastructure with SLAC's parallel finite element electromagnetic codes, thus enabling all physics simulations within a single framework. The finite-element approach allows high-fidelity, high-accuracy simulations and the parallel implementation facilitates large-scale computation with fast turnaround times. In this paper, TEM3P is used to analyze thermal loading at coupler end of the JLAB SRF cavity

Inversion of Airborne Contaminants in a Regional Model

We are interested in a DDDAS problem of localization of airborne contaminant releases in regional atmospheric transport models from sparse observations. Given measurements of the contaminant over an observation window at a small number of points in space, and a velocity field as predicted for example by a mesoscopic weather model, we seek an estimate of the state of the contaminant at the beginning of the observation interval that minimizes the least squares misfit between measured and predicted contaminant field, subject to the convection-diffusion equation for the contaminant. Once the ''initial'' conditions are estimated by solution of the inverse problem, we issue predictions of the evolution of the contaminant, the observation window is advanced in time, and the process repeated to issue a new prediction, in the style of 4D-Var. We design an appropriate numerical strategy that exploits the spectral structure of the inverse operator, and leads to efficient and accurate resolution of the inverse problem. Numerical experiments verify that high resolution inversion can be carried out rapidly for a well-resolved terrain model of the greater Los Angeles area

Computational Science Research in Support of Petascale Electromagnetic Modeling

Computational science research components were vital parts of the SciDAC-1 accelerator project and are continuing to play a critical role in newly-funded SciDAC-2 accelerator project, the Community Petascale Project for Accelerator Science and Simulation (ComPASS). Recent advances and achievements in the area of computational science research in support of petascale electromagnetic modeling for accelerator design analysis are presented, which include shape determination of superconducting RF cavities, mesh-based multilevel preconditioner in solving highly-indefinite linear systems, moving window using h- or p- refinement for time-domain short-range wakefield calculations, and improved scalable application I/O

Analysis of the Cause of High External Q Modes in the JLab High Gradient Prototype Cryomodule Renascence

The Renascence cryomodule [1] installed in CEBAF in 2007 consists of 8 cavities as shown in Figure 1. The first three cavities (No.1-No.3) in the upstream end are of the Low Loss (LL) shape design, and the remaining 5 cavities (No.4-No.8) on the beam downstream end are the High Gradient (HG) shape design. The fundamental power couplers (FPCs) are the rectangular waveguides, and the little cylindrical structures are the HOM couplers. The locations of the FPC in the last four cavities are mirrored about the beam z axis. Cavities No.4 and No.5 form a back-to-back cavity pair. Among the HG cavities installed in the Renascence cryomodule, the only identifiable difference from their fabrication documentation is that cavity No.5 received an extra EBW pass on one equator weld, specifically cell 5. The non-uniform mechanical tuning required to compensate the fundamental mode tune and flatness for the extra shrinkage of this cell is believed to contribute the most significant differences from the other HG cavities. Beam based instability studies on this cryomodule in CEBAF have shown a significant beam breakup (BBU) threshold current reduction, well below design value. Frequency spectrum peaked by the off-sided beam power indicated the cause is due to abnormal high Q modes in the cavity No.5. Measured beam off-axis position at the cavity No.5 does not correspond to the shunt impedances calculated for an ideal cavity. Low power RF measurements have identified that the problematic modes are in the second dipole band (TM110 like). Three of the modes have external Qs two orders magnitude higher than the others, while the rest of modes in the first two dipole bands are normal in terms of the design values. The cause of this abnormality and the future impact on the BBU was not able to be resolved due to the limitations of information that can be obtained from the measurements. It is important to understand the cause of this abnormality so that effective QA/QC measures can be implemented to avoid such problem in the final upgrade design and manufacture. The goal of this work is to utilize advanced simulation tools to understand the high external Q (Q{sub ext}) problem observed in the Renascence cryomodule. In the past years, SLAC has built a set of state-of-the-art advanced simulation tools based on finite-element unstructured meshes and parallel computation implementations on supercomputers [2, 3]. The codes are capable of simulating large complex RF systems with unprecedented resolution and turnaround time. They have been successfully applied to many existing and future accelerator R&D projects to improve the machine performance and to optimize the designs. These tools are essential to perform accurate full system analyses such as the JLab's SRF cavities. We will use the simulation results and the data from the RF measurements to gain a better understanding of the cavity performance and tolerance issues and provide a solid foundation to do the BBU simulation and prediction for the 12GeV Upgrade project by using JLab's BBU codes. In this report, we will focus on the following two main tasks: (1) Ideal cavity simulation--to evaluate the effectiveness of the damping by the higher-order-mode (HOM) couplers, and search for possible trapped modes in a back-to-back cavity pair (e.g. cavity No.4 & No.5). (2) Abnormal cavity study--to understand the cause of the high Q{sub ext} modes in cavity No.5 using an advanced Shape Determination Tool