Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn angular momentum conserving Affine-Particle-In-Cell method
91
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a new technique for transferring momentum and velocity between particles and grid with Particle-In-Cell (PIC) calculations which we call Affine-Particle-In-Cell (APIC). APIC represents particle velocities as locally affine, rather than locally constant as in traditional PIC. We show that this representation allows APIC to conserve linear and angular momentum across transfers while also dramatically reducing numerical diffusion usually associated with PIC. Notably, conservation is achieved with lumped mass, as opposed to the more commonly used Fluid Implicit Particle (FLIP) transfers which require a full mass matrix for exact conservation. Furthermore, unlike FLIP, APIC retains a filtering property of the original PIC and thus does not accumulate velocity modes on particles as FLIP does. In particular, we demonstrate that APIC does not experience velocity instabilities that are characteristic of FLIP in a number of Material Point Method (MPM) hyperelasticity calculations. Lastly, we demonstrate that when combined with the midpoint rule for implicit update of grid momentum that linear and angular momentum are exactly conserved.
rate research
Read More
The recently developed energy conserving semi-implicit method (ECsim) for PIC simulation is applied to multiple scale problems where the electron-scale physics needs to be only partially retained and the interest is on the macroscopic or ion-scale processes. Unlike hybrid methods, the ECsim is capable of providing kinetic electron information, such as wave-electron interaction (Landau damping or cyclotron resonance) and non-Maxwellian electron velocity distributions. However, like hybrid, the ECsim does not need to resolve all electron scales, allowing time steps and grid spacing orders of magnitude larger than in explicit PIC schemes. The additional advantage of the ECsim is that the stability at large scale is obtained while conserving energy exactly. Three examples are presented: ion acoustic waves, electron acoustic instability and reconnection processes.
Based on the previously developed Energy Conserving Semi Implicit Method (ECsim) code, we present its cylindrical implementation, called ECsim-CYL, to be used for axially symmetric problems. The main motivation for the development of the cylindrical version is to greatly improve the computational speed by utilizing cylindrical symmetry. The ECsim-CYL discretizes the field equations in two-dimensional cylindrical coordinates using the finite volume method . For the particle mover, it uses a modification of ECsims mover for cylindrical coordinates by keeping track of all three components of velocity vectors, while only keeping radial and axial coordinates of particle positions. In this paper, we describe the details of the algorithm used in the ECsim-CYL and present a series of tests to validate the accuracy of the code including a wave spectrum in a homogeneous plasmas inside a cylindrical waveguide and free expansion of a spherical plasma ball in vacuum. The ECsim-CYL retains the stability properties of ECsim and conserves the energy within machine precision, while accurately describing the plasma behavior in the test cases.
The Levi-Malcev decomposition is applied to bosonic models of quantum mechanics based on unitary Lie algebras u(2), u(2)+u(2), u(3) and u(4) to clearly disentangle semisimple subalgebras. The theory of weighted Dynkin diagrams is then applied to identify conjugacy classes of relevant A_1 subalgebras allowing to introduce a complete classification of new angular momentum non conserving (AMNC) dynamical symmetries. The tensor analysis of the whole algebra based on the new angular momentum operators reveals unexpected spinors to occur in purely bosonic models. The new chains of subalgebra can be invoked to set up ANMC bases for diagonalization.
Particle-in-Cell (PIC) methods are widely used computational tools for fluid and kinetic plasma modeling. While both the fluid and kinetic PIC approaches have been successfully used to target either kinetic or fluid simulations, little was done to combine fluid and kinetic particles under the same PIC framework. This work addresses this issue by proposing a new PIC method, PolyPIC, that uses polymorphic computational particles. In this numerical scheme, particles can be either kinetic or fluid, and fluid particles can become kinetic when necessary, e.g. particles undergoing a strong acceleration. We design and implement the PolyPIC method, and test it against the Landau damping of Langmuir and ion acoustic waves, two stream instability and sheath formation. We unify the fluid and kinetic PIC methods under one common framework comprising both fluid and kinetic particles, providing a tool for adaptive fluid-kinetic coupling in plasma simulations.
We present a generic technique, automated by computer-algebra systems and available as open-source software cite{scuff-em}, for efficient numerical evaluation of a large family of singular and nonsingular 4-dimensional integrals over triangle-product domains, such as those arising in the boundary-element method (BEM) of computational electromagnetism. To date, practical implementation of BEM solvers has often required the aggregation of multiple disparate integral-evaluation schemes to treat all of the distinct types of integrals needed for a given BEM formulation; in contrast, our technique allows many different types of integrals to be handled by the emph{same} algorithm and the same code implementation. Our method is a significant generalization of the Taylor--Duffy approach cite{Taylor2003,Duffy1982}, which was originally presented for just a single type of integrand; in addition to generalizing this technique to a broad class of integrands, we also achieve a significant improvement in its efficiency by showing how the emph{dimension} of the final numerical integral may often be reduced by one. In particular, if $n$ is the number of common vertices between the two triangles, in many cases we can reduce the dimension of the integral from $4-n$ to $3-n$, obtaining a closed-form analytical result for $n=3$ (the common-triangle case).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
