Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userImproved Numerical Method for Solution of the Haissinski Equation
110
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The longitudinal charge density of an electron beam in its equilibrium state is given by the solution of the Haissinski equation, which provides a stationary solution of the Vlasov-Fokker-Planck equation. The physical input is the longitudinal wake potential. We formulate the Haissinski equation as a nonlinear integral equation with the normalization integral stated as a functional of the solution. This equation can be solved in a simple way by the matrix version of Newtonss iteration, beginning with the Gaussian as a first guess. We illustrate for several quasi-realistic wake potentials. Convergence is extremely robust, even at currents much higher than nominal for the storage rings considered. The method overcomes limitations of earlier procedures, and provides the convenience of automatic normalization of the solution.
rate research
Read More
An effective algorithm is presented for solving the Beltrami equation fzbar = mu fz in a planar disk. The algorithm involves no evaluation of singular integrals. The strategy, working in concentric rings, is to construct a piecewise linear mu-conformal mapping and then correct the image using a known algorithm for conformal mappings. Numerical examples are provided and the computational complexity is analyzed.
We numerically solve the Boltzmann equation for trapped fermions in the normal phase using the test-particle method. After discussing a couple of tests in order to estimate the reliability of the method, we apply it to the description of collective modes in a spherical harmonic trap. The numerical results are compared with those obtained previously by taking moments of the Boltzmann equation. We find that the general shape of the response function is very similar in both methods, but the relaxation time obtained from the simulation is significantly longer than that predicted by the method of moments. It is shown that the result of the method of moments can be corrected by including fourth-order moments in addition to the usual second-order ones and that this method agrees very well with our numerical simulations.
Two numerical methods are used to evaluate the relativistic spectrum of the two-centre Coulomb problem (for the $H_{2}^{+}$ and $Th_{2}^{179+}$ diatomic molecules) in the fixed nuclei approximation by solving the single particle time-independent Dirac equation. The first one is based on a min-max principle and uses a two-spinor formulation as a starting point. The second one is the Rayleigh-Ritz variational method combined with kinematically balanced basis functions. Both methods use a B-spline basis function expansion. We show that accurate results can be obtained with both methods and that no spurious states appear in the discretization process.
We describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for an interacting Bose gas in a harmonic oscillator potential, with the inclusion of a long-range dipolar interaction. The central difficulty in solving this equation is the requirement that the field is restricted to a small set of prescribed modes that constitute the low energy c-field region of the system. We present a scheme, using a Hermite-polynomial based spectral representation, that precisely implements this mode restriction and allows an efficient and accurate solution of the dipolar PGPE. We introduce a set of auxiliary oscillator states to perform a Fourier transform necessary to evaluate the dipolar interaction in reciprocal space. We extensively characterize the accuracy of our approach, and derive Ehrenfest equations for the evolution of the angular momentum.
The validation and parallel implementation of a numerical method for the solution of the time-dependent Dirac equation is presented. This numerical method is based on a split operator scheme where the space-time dependence is computed in coordinate space using the method of characteristics. Thus, most of the steps in the splitting are calculated exactly, making for a very efficient and unconditionally stable method. We show that it is free from spurious solutions related to the fermion-doubling problem and that it can be parallelized very efficiently. We consider a few simple physical systems such as the time evolution of Gaussian wave packets and the Klein paradox. The numerical results obtained are compared to analytical formulas for the validation of the method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
