Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical solver”

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 955 records · Page 53

Computational Evaluation of Thermal Response of Open-Cell Foam With Circular Pore (Computional Evaluation of Thermal Response of Open-Cell Foam with Circular Pore)

The evaluation of effective material properties in heterogeneous materials (e.g., composites or multicomponent structures) critically relevant to a wide spectrum of applications, including nuclear power, electronic packaging, flame retardants, hypersonics, and gas turbine power. The work described in this paper is centered around the numerical assessment of the thermal behavior of porous materials obtained from finite element thermal modeling and simulation. Here, two-dimensional, steady state analyses were performed on unit cells with centered, circular pores using a second order accurate Galerkin finite element method (FEM). The effective thermal conductivities of the porous systems were examined, encompassing a range of porosities from 4.9% to 60.1%. The geometries of the models were generated based on ordered circular pores for each modeled porosity level. The system response quantity (SRQ) under investigation was the dimensionless effective thermal conductivity across the unit cell. The dimensionless effective thermal conductivity was compared across all simulated cases, producing a trend between porosity and effective thermal conductivity. In the presented investigation, the method of manufactured solutions (MMS) was used to perform code verification, and the grid convergence index (GCI) was employed to estimate discretization uncertainty as solution verification. Code verification concluded an approximately second order accurate Galerkin FEM solver. It was found that the introduction of porosity to the unit cell material structure reduces effective thermal conductivity, as anticipated. Numerical results obtained in this study are compared to an analytical solution and to a sample of empirical data.

36 MATERIALS SCIENCE↗

A Butterfly-Accelerated Volume Integral Equation Solver for Broad Permittivity and Large-Scale Electromagnetic Analysis

In this work, a butterfly-accelerated volume integral equation (VIE) solver is proposed for fast and accurate electromagnetic (EM) analysis of scattering from heterogeneous objects. The proposed solver leverages the hierarchical off-diagonal butterfly (HOD-BF) scheme to construct the system matrix and obtain its approximate inverse, used as a preconditioner. Complexity analysis and numerical experiments validate the O(N log 2 N) construction cost of the HOD-BF-compressed system matrix and O(N log 1.5 N) inversion cost for the preconditioner, where N is the number of unknowns in the high-frequency EM scattering problem. For many practical scenarios, the proposed VIE solver requires less memory and computational time to construct the system matrix and obtain its approximate inverse compared to a H matrix-accelerated VIE solver. The accuracy and efficiency of the proposed solver have been demonstrated via its application to the EM analysis of large-scale canonical and real-world structures comprising of broad permittivity values and involving millions of unknowns.

42 ENGINEERING↗

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗

A Fast Solver for Implicit Integration of the Vlasov--Poisson System in the Eulerian Framework

In this paper, we present a domain decomposition algorithm to accelerate the solution of Eulerian-type discretizations of the linear, steady-state Vlasov equation. The steady-state solver then forms a key component in the implementation of fully implicit or nearly fully implicit temporal integrators for the nonlinear Vlasov--Poisson system. The solver relies on a particular decomposition of phase space that enables the use of sweeping techniques commonly used in radiation transport applications. The original linear system for the phase space unknowns is then replaced by a smaller linear system involving only unknowns on the boundary between subdomains, which can then be solved efficiently with Krylov methods such as GMRES. Steady-state solves are combined to form an implicit Runge--Kutta time integrator, and the Vlasov equation is coupled self-consistently to the Poisson equation via a linearized procedure or a nonlinear fixed-point method for the electric field. Finally, numerical results for standard test problems demonstrate the efficiency of the domain decomposition approach when compared to the direct application of an iterative solver to the original linear system.

97 MATHEMATICS AND COMPUTING↗

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗

Characteristics of temperature rise in variable inductor employing magnetorheological fluid driven by a high-frequency pulsed voltage source

A variable inductor with magnetorheological (MR) fluid has been successfully applied to power electronics applications; however, its thermal characteristics have not been investigated. To evaluate the performance of the variable inductor with respect to temperature, we measured the characteristics of temperature rise and developed a numerical analysis technique. The characteristics of temperature rise were determined experimentally and verified numerically by adopting a multiphysics analysis technique. In order to accurately estimate the temperature distribution in a variable inductor with an MR fluid-gap, the thermal solver should import the heat source from the electromagnetic solver to solve the eddy current problem. To improve accuracy, the B–H curves of the MR fluid under operating temperature were obtained using the magnetic property measurement system. In addition, the Steinmetz equation was applied to evaluate the core loss in a ferrite core. The predicted temperature rise for a variable inductor showed good agreement with the experimental data and the developed numerical technique can be employed to design a variable inductor with a high-frequency pulsed voltage source.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Improved numerical methods for turbulent viscous recirculating flows

The performance of discrete methods for the prediction of fluid flows can be enhanced by improving the convergence rate of solvers and by increasing the accuracy of the discrete representation of the equations of motion. This report evaluates the gains in solver performance that are available when various acceleration methods are applied. Various discretizations are also examined and two are recommended because of their accuracy and robustness. Insertion of the improved discretization and solver accelerator into a TEACH mode, that has been widely applied to combustor flows, illustrates the substantial gains to be achieved.

Turan, A.↗

A three-dimensional application with the numerical grid generation code: EAGLE (utilizing an externally generated surface)

Program EAGLE (Eglin Arbitrary Geometry Implicit Euler) is a multiblock grid generation and steady-state flow solver system. This system combines a boundary conforming surface generation, a composite block structure grid generation scheme, and a multiblock implicit Euler flow solver algorithm. The three codes are intended to be used sequentially from the definition of the configuration under study to the flow solution about the configuration. EAGLE was specifically designed to aid in the analysis of both freestream and interference flow field configurations. These configurations can be comprised of single or multiple bodies ranging from simple axisymmetric airframes to complex aircraft shapes with external weapons. Each body can be arbitrarily shaped with or without multiple lifting surfaces. Program EAGLE is written to compile and execute efficiently on any CRAY machine with or without Solid State Disk (SSD) devices. Also, the code uses namelist inputs which are supported by all CRAY machines using the FORTRAN Compiler CF177. The use of namelist inputs makes it easier for the user to understand the inputs and to operate Program EAGLE. Recently, the Code was modified to operate on other computers, especially the Sun Spare4 Workstation. Several two-dimensional grid configurations were completely and successfully developed using EAGLE. Currently, EAGLE is being used for three-dimension grid applications.

Houston, Johnny L.↗

Validation of a Multiphase Computational Fluid Dynamics Model for Vapor Pull-Through in Normal and Low Gravity

On-orbit fluid transfer such as refueling of propellant tanks and life-support systems can enable long-duration space missions. For safe and efficient liquid transfer operations, prior knowledge of liquid positioning and liquid-vapor interface behavior while draining in a low-gravity environment is required. Numerical models capable of predicting vapor ingestion (or vapor pull-through) can be used to design liquid transfer operations while reducing liquid residuals, mission risk and settling thrust required to prevent vapor ingestion. An experimental program conducted in the 2.2 Second Drop Tower facility at NASA Glenn Research Center investigated the vapor ingestion phenomenon for a range of outflow rates and tank sizes providing a database for validation. This study presents a Computational Fluid Dynamics model capable of accurately predicting the vapor ingestion using the Volume-of-Fluid multiphase solver in commercial code STAR-CCM+. A description of the experimental setup and general trends from similar studies are presented. Comparisons of the numerical prediction and test data in normal and low gravity show good agreement and give confidence in pursuing design of full-scale propellant transfer systems.

Computational Fluid Dynamics↗

Diffusion Synthetic Acceleration for High Order S{sub N} Transport on Meshes with Curved Surfaces

Thermal radiation transport is integral to multiphysics applications such as inertial confinement fusion and astrophysics. Materials in these problems can be exceptionally optically thick and diffusive. A transport solver must be capable of accurate and efficient solutions to such problems. Previously, Woods et al. established a proof-of-concept for high order (HO) discontinuous Galerkin finite element (DGFEM) transport on meshes with curved surfaces. They numerically demonstrated the spatial convergence is order p + 1. Remarkably, these convergence rates are relatively invariant to the amount of curvature of the spatial mesh. We present here an extension of this work that demonstrates a rapidly-convergent iterative solver for the FEM transport equations on meshes with curved sides. We use the modified interior penalty (MIP) diffusion synthetic acceleration (DSA) established by Wang and Ragusa. The partially consistent MIP equations stem from the symmetric interior penalty (IP) derivation of the diffusion equation. The IP method is not stable for optically thick cells so it was combined with the diffusion conforming form (DCF), a discretized diffusion equation derived from the discretized S N equations that is not stable for intermediate optical thicknesses. The MIP equations adapt a 'switch' to the IP method in optically thinner regions and the DCF method in optically thicker regions. The literature does not reveal any Fourier analyses performed for iterative transport on meshes with curved surfaces. Wang and Ragusa performed a Fourier analysis for orthogonal meshes and our results appear to be as expected. Previous literature has reported on some of the sensitivities that affect the spectral radius of DSA schemes: polynomial order, and scattering ratios. Presently, we focus on the sensitivity of the spectral radius to polynomial order, scattering ratio, the constant C, and mesh curvature for varying cell thicknesses. We also numerically study the behavior of the spectral radius to verify our DSA implementation. A sequence of problems is designed to exercise the acceleration. As the problems become more optically thick and diffusive, the accelerated solution continues to converge toward the analytic solution. We employ the Modular Finite Element Methods (MFEM) package, developed at Lawrence Livermore National Laboratory, to generate the linear algebraic system of equations for each angle in the discrete ordinates quadrature set, and use UMFPack to solve this system of equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Adjoint sensitivity analysis of chaotic dynamical systems with non-intrusive least squares shadowing

This paper presents a discrete adjoint version of the recently developed non-intrusive least squares shadowing (NILSS) algorithm, which circumvents the instability that conventional adjoint methods encounter for chaotic systems. The NILSS approach involves solving a smaller minimization problem than other shadowing approaches and can be implemented with only minor modifications to preexisting tangent and adjoint solvers. Adjoint NILSS is demonstrated on a small chaotic ODE, a one-dimensional scalar PDE, and a direct numerical simulation (DNS) of the minimal flow unit, a turbulent channel flow on a small spatial domain. This is the first application of an adjoint shadowing-based algorithm to a three-dimensional turbulent flow. - Highlights: • A discrete adjoint non-intrusive least squares shadowing (NILSS) is presented. • The NILSS approach is closely related to multiple shooting shadowing (MSS). • Adjoint NILSS prevents exponential growth in time of the adjoint field. • Adjoint NILSS is demonstrated on a simulation of wall-bounded turbulent flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Grain2mesh: A Python and cubit mesh generator from unprocessed mesoscale images

Predicting bulk behavior from microscale features constitutes a key objective in multiscale modeling research, often involving numerical models composed of finite elements that capture the diversity of constituent phases, shapes, and orientations within the material. The Grain2mesh toolbox allows the user to input unprocessed mesoscopic images for automatic segmentation, pre-processing, quality control, and numerical mesh generation. The numerical mesh generation incorporates Cubit routines to generate robust multi-phase mesh structure for use in computational mechanics solvers. The python classes developed contain detailed documentation and examples to support standard usage and case-specific alternative options.

58 GEOSCIENCES↗

Aerodynamic and aeroacoustic for wind turbine

This paper describes a hybrid approach forpredicting noise radiated from the rotating Wind Turbine (HAWT) blades, where the sources are extracted from an unsteady Reynolds-Averaged-Navier Stocks (URANS) simulation, ANSYS CFX 11.0, was used to calculate The near-field flow parameters around the blade surface that are necessary for FW-H codes. Comparisons with NREL Phase II experimental results are presented with respect to the pressure distributions for validating a capacity of the solver to calculate the near-field flow on and around the wind turbine blades, The results show that numerical data have a good agreement with experimental. The acoustic pressure, presented as a sum of thickness and loading noise components, is analyzed by means of a discrete fast Fourier transformation for the presentation of the time acoustic time histories in the frequency domain. The results convincingly show that dipole source noise is the dominant noise source for this wind turbine.

17 WIND ENERGY↗

Acoustic Codes in 2D Spherical Coordinate

Finite-difference methods are widely used to simulate infrasound propagation in the atmosphere. Flexibility of finite-difference scheme allows implementation of highly heterogeneous media for sound propagation as well as complex source models for sound generation. While full 3-D finite-difference methods have been utilized for local infrasound propagation with pronounced topography, 2-D modeling approach has been preferred for regional and global propagation as full 3-D methods generally require enormous computational resources. Infrasound propagation is often simulated with a second-order finite difference scheme. This lowest-order finite-difference scheme is robust and straightforward to implement complex boundary conditions, but the solution includes large error with numerical dispersion and dissipation. This large numerical error may make the second-order finite-difference not optimal for long range infrasound propagation modeling as the numerical dispersion degrades the accuracy of the solution unacceptably. Here, we developed a high-order finite-difference solver for long-range infrasound simulation. The high-order scheme is particularly popular for linear wave modeling in aeroacoustics owing to its low-dispersive and low-dissipative behavior. We develop and evaluate a high-order finite difference scheme in 2-D axisymmetric coordinates. The axisymmetry allows to approximate 3-D spherical sound propagation and amplitude attenuation by a 2-D method. AC2Dr is developed to simulate infrasound propagation in realistic atmosphere, but can be used for linear acoustic waves in general materials with background flow. AC2Dr in an axisymmetric coordinates allows for spherical radiation of acoustic waves from compact sources.

Sjogreen, Bjorn↗

Investigation of Microbunching Instabilities in Modern Recirculating Accelerators

Particle accelerators are machines to accelerate and store charged particles, such as electrons or protons, to the energy levels for various scientific applications. A collection of charged particles usually forms a particle beam. There are three basic types of particle accelerators: linear accelerators (linac), storage-ring (or circular) accelerators, and recirculating accelerators. In a linac, particles are accelerated and pass through once along a linear or straight beamline. Storage-ring accelerators propel particles around a circular track and repetitively append the energy to the stored beam. The third type, also the most recent one in chronology, the recirculating accelerator, is designed to accelerate the particle beam in a short section of linac, circulate the beam, and then either continue to accelerate for energy boost or decelerate it for energy recovery. The beam properties of a linac machine are set at best by the initial particle sources. For storage rings, the beam equilibria are instead determined by the overall machine design. The modern recirculating machines share with linacs the advantages to both accelerate and preserve the beam with high beam quality, as well as efficiently reuse the accelerating components. The beamline design in such a machine configuration can however be much more complicated than that of linacs. As modern accelerators push toward the high-brightness or high-intensity frontier by demanding particles in a highly charged bunch (about nano-Coulomb per bunch) to concentrate in an ever-decreasing beam phase space (transverse normalized emittance about 1 μm and relative energy spread of the order of 10^-5 in GeV beam energy), the interaction amongst particles via their self-generated electromagnetic fields can potentially lead to coherent instabilities of the beam and thus pose significant challenges to the machine design and operation. In the past decade and a half, microbunching instability (MBI) has been one of the most challenging issues for such high-brightness or high-intensity beam transport, as it would degrade lasing performance in the fourth-generation light sources, reduce cooling efficiency in electron cooling facilities, and compromise the luminosity of colliding beams in lepton or lepton-hadron colliders. The dissertation work will focus on the MBI in modern recirculating electron accelerators. It has been known that the collective interactions, the coherent synchrotron radiation (CSR) and the longitudinal space charge (LSC) forces, can drive MBI. The CSR effect is a collective phenomenon in which the electrons in a curved motion, e.g. a bending dipole, emit radiation at a scale comparable to the micro-bunched structure of the bunch distribution. The LSC effect stems from non-uniformity of the charge distribution, acts as plasma oscillation, and can eventually accumulate an amount of energy modulation when the beam traverses a long section of a beamline. MBI can be seeded by non-uniformity or shot noise of the beam, which originates from granularity of the elementary charge. Through the aforementioned collective effects, the modulation of the bunch sub-structure can be amplified and, once the beam-wave interaction formed a positive feedback, can result in MBI. The problem of MBI has been intensively studied for linac-based facilities and for storage-ring accelerators. However, systematic studies for recirculation machines are still very limited and form a knowledge gap. Because of the much more complicated machine configuration of the recirculating accelerators than that of linacs, the existing MBI analysis needs to be extended to accommodate the high-brightness particle beam transport in modern recirculating accelerators. This dissertation is focused on theoretical investigation of MBI in such machine configuration in the following seven themes: (1) Development and generalization of MBI theory The theoretical formulation has been extended so as to be applicable to a general linear beamline lattice including horizontal and vertical transport bending elements, and beam acceleration or deceleration. These featured generalizations are required for MBI analysis in recirculation accelerators. (2) Construction of CSR impedance models In addition to the steady-state CSR interaction, it has been found that the exit transient effect (or CSR drift) can even result in more serious MBI in high-brightness recirculation arcs. The onedimensional free-space CSR impedances, especially the exit transients, are derived. The steady-state CSR impedance is also extended to non-ultrarelativistic beam energy for MBI analysis of low-energy merger sections in recirculating accelerators. (3) Numerical implementation of the derived semi-analytical formulation This includes the development of a semi-analytical Vlasov solver for MBI analysis, and also benchmarking of the solver against massive particle tracking simulations. (4) Exploration of multistage amplification behavior of CSR microbunching development The CSR-induced MBI acts as an amplifier, which amplifies the sub-bunch modulation of a beam. The amplification is commonly quantified by the amplification gain. A beam transport system can be considered as a cascaded amplifier. Unlike the two-stage amplification of four-dipole bunch compressor chicanes employed in linacs, the recirculation arcs, which are usually constituted by several tens of bending magnets, show a distinguishing feature of up to six-stage microbunching amplification for our example arc lattices. That is, the maximal CSR amplification gain can be proportional to the peak bunch current up to sixth power. A method to compare lattice performance has been developed in terms of gain coefficients, which nearly depend on the lattice properties only. This method has also proven to be an effective way to quantify the current dependence of the maximal (5) Control of CSR MBI in multibend transport or recirculation arcs The existing mitigation schemes of MBI mostly aim to linac-based accelerators and may not be practical to the recirculating accelerator facilities. Thus a set of conditions for suppression of CSR MBI was proposed and examined for example lattices from low (~100 MeV) to high (~1 GeV) energies. (6) Study of more aspects of microbunched structures in beam phase spaces For a cascaded amplifier in circuit electronics, the total amplification gain can be estimated as the product of individual gains. In a beam transport line of an accelerator, the (scalar) gain multiplication was examined and found to under-estimate the overall microbunching amplification. The concept of gain matrix was developed, which includes the density, energy and transverse-longitudinal modulations in a beam phase space, and used to analyze MBI for a proposed recirculating machine. Throughout the gain matrix approach, it reasonably gives the upper limit of spectral MBI gain curves. This extended analysis can be employed to study multi-pass recirculation. (7) Study of MBI for magnetized beams Driven by a recent energy-recovery-linac based cooler design for electron cooling at Jefferson Lab Electron-Ion Collider Project, the generalized theoretical formulation for MBI to a transversely coupled beam has been developed and applied to this study. A magnetized beam in general features non-zero canonical angular momentum, thus considered to be a transversely coupled beam. A novel idea of utilizing magnetized beam transport was proposed for improvement of cooling efficiency and possible mitigation of collective effects. A concern of MBI regarding this design was studied and excluded. The large transverse beam size associated with the beam magnetization is found to help suppress MBI via the transverse-longitudinal correlation.

43 PARTICLE ACCELERATORS↗

Finite dynamic element formulation for a plane triangular element

Higher order correction terms for the stiffness and inertia matrices associated with a triangular plane stress-strain finite dynamic element are developed in detail. Numerical results are presented which indicate that the adoption of these matrices along with a suitable quadratic matrix eigenproblem solver effects a significant economy in the free vibration solution of structure when compared with the analysis based on the usual finite element procedure. Finally, a FORTRAN IV computer program listing of the various relevant matrices is given.

Gupta, K. K.↗