Excessive Memory Usage of the ELLPACK Sparse Matrix Storage Scheme throughout the Finite Element Computations

Sparse matrices are occasionally encountered during solution of various problems by means of numerical methods, particularly the finite element method. ELLPACK sparse matrix storage scheme, one of the most widely used methods due to its implementation ease, is investigated in this study. The scheme uses excessive memory due to its definition. For the conventional finite element method, where the node elements are used, the excessive memory caused by redundant entries in the ELLPACK sparse matrix storage scheme becomes negligible for large scale problems. On the other hand, our analyses show that the redundancy is still considerable for the occasions where facet or edge elements have to be used

Spin-Mediated Consciousness: Theory, Experimental Studies, Further Development & Related Topics

We postulate that consciousness is intrinsically connected to quantum spin since the latter is the origin of quantum effects in both Bohm and Hestenes quantum formulisms and a fundamental quantum process associated with the structure of space-time. Applying these ideas to the particular structures and dynamics of the brain, we have developed a detailed model of quantum consciousness. We have also carried out experiments from the perspective of our theory to test the possibility of quantum-entangling the quantum entities inside the brain with those of an external chemical substance. We found that applying magnetic pulses to the brain when an anaesthetic was placed in between caused the brain to feel the effect of said anaesthetic as if the test subject had actually inhaled the same. We further found that drinking water exposed to magnetic pulses, laser light or microwave when an anaesthetic was placed in between also causes brain effects in various degrees. Additional experiments indicate that the said brain effect is indeed the consequence of quantum entanglement. Recently we have studied non-local effects in simple physics systems. We have found that the pH value, temperature and gravity of a liquid in the detecting reservoirs can be non-locally affected through manipulating another liquid in a remote reservoir quantum-entangled with the former. In particular, the pH value changes in the same direction as that being manipulated; the temperature can change against that of local environment; and the gravity can change against local gravity. We suggest that they are mediated by quantum entanglement between nuclear and/or electron spins in treated liquid and discuss the profound implications of these results. This paper now also includes materials on further development of the theory and related topics.Comment: 92 pages; expanded content; minor corrections; for additional information, please visit http://quantumbrain.or

A Critical Behaviour of Anomalous Currents, Electric-Magnetic Universality and CFT_4

We discuss several aspects of superconformal field theories in four dimensions (CFT_4), in the context of electric-magnetic duality. We analyse the behaviour of anomalous currents under RG flow to a conformal fixed point in N=1, D=4 supersymmetric gauge theories. We prove that the anomalous dimension of the Konishi current is related to the slope of the beta function at the critical point. We extend the duality map to the (nonchiral) Konishi current. As a byproduct we compute the slope of the beta function in the strong coupling regime. We note that the OPE of $T_{\mu\nu}$ with itself does not close, but mixes with a special additional operator $\Sigma$ which in general is the Konishi current. We discuss the implications of this fact in generic interacting conformal theories. In particular, a SCFT_4 seems to be naturally equipped with a privileged off-critical deformation $\Sigma$ and this allows us to argue that electric-magnetic duality can be extended to a neighborhood of the critical point. We also stress that in SCFT_4 there are two central charges, c and c', associated with the stress tensor and $\Sigma$, respectively; c and c' allow us to count both the vector multiplet and the matter multiplet effective degrees of freedom of the theory.Comment: harvmac tex, 28 pages, 3 figures. Version to be published in Nucl. Phys.

A unified gas kinetic scheme for transport and collision effects in plasma

In this study, the Vlasov-Poisson equation with or without collision term for plasma is solved by the unified gas kinetic scheme (UGKS). The Vlasov equation is a differential equation describing time evolution of the distribution function of plasma consisting of charged particles with long-range interaction. The distribution function is discretized in discrete particle velocity space. After the Vlasov equation is integrated in finite volumes of physical space, the numerical flux across a cell interface and source term for particle acceleration are computed to update the distribution function at next time step. The flux is decided by Riemann problem and variation of distribution function in discrete particle velocity space is evaluated with central difference method. A electron-ion collision model is introduced in the Vlasov equation. This finite volume method for the UGKS couples the free transport and long-range interaction between particles. The electric field induced by charged particles is controlled by the Poisson's equation. In this paper, the Poisson's equation is solved using the Green's function for two dimensional plasma system subjected to the symmetry or periodic boundary conditions. Two numerical tests of the linear Landau damping and the Gaussian beam are carried out to validate the proposed method. The linear electron plasma wave damping is simulated based on electron-ion collision operator. Compared with previous methods, it is shown that the current method is able to obtain accurate results of the Vlasov-Poisson equation with a time step much larger than the particle collision time. Highly non-equilibrium and rarefied plasma flows, such as electron flows driven by electromagnetic field, can be simulated easily.Comment: 33 pages, 13 figure

Current Flow and Pair Creation at Low Altitude in Rotation Powered Pulsars' Force-Free Magnetospheres: Space-Charge Limited Flow

(shortened) We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of Rotation Powered Pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We observe novel behavior. a) When the current density is less than the Goldreich-Julian (GJ) value (0<j/j_{GJ}<1), space charge limited acceleration of the current carrying beam is mild, with the full GJ charge density being comprised of the charge density of the beam, co-existing with a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are on the order of mc^2/e, and pair creation is absent. b) When the current density exceeds the GJ value (j/j_{GJ}>1), the system develops high voltage drops, causing emission of gamma rays and intense bursts of pair creation. The bursts exhibit limit cycle behavior, with characteristic time scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). c) In return current regions, where j/j_{GJ}<0, the system develops similar bursts of pair creation. In cases b) and c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady, finding that such steady flows can occupy only a small fraction of the current density parameter space of the force-free magnetospheric model. The generic polar flow dynamics and pair creation is strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage as a function of the applied current.Comment: 35 pages, 29 figures. Accepted for publication in MNRAS. Added new appendix, several minor changes in the tex

Distributed-memory parallelization of an explicit time-domain volume integral equation solver on Blue Gene/P

Two distributed-memory schemes for efficiently parallelizing the explicit marching-on in-time based solution of the time domain volume integral equation on the IBM Blue Gene/P platform are presented. In the first scheme, each processor stores the time history of all source fields and only the computationally dominant step of the tested field computations is distributed among processors. This scheme requires all-to-all global communications to update the time history of the source fields from the tested fields. In the second scheme, the source fields as well as all steps of the tested field computations are distributed among processors. This scheme requires sequential global communications to update the time history of the distributed source fields from the tested fields. Numerical results demonstrate that both schemes scale well on the IBM Blue Gene/P platform and the memory efficient second scheme allows for the characterization of transient wave interactions on composite structures discretized using three million spatial elements without an acceleration algorithm

Solving for Micro- and Macro- Scale Electrostatic Configurations Using the Robin Hood Algorithm

We present a novel technique by which highly-segmented electrostatic configurations can be solved. The Robin Hood method is a matrix-inversion algorithm optimized for solving high density boundary element method (BEM) problems. We illustrate the capabilities of this solver by studying two distinct geometry scales: (a) the electrostatic potential of a large volume beta-detector and (b) the field enhancement present at surface of electrode nano-structures. Geometries with elements numbering in the O(10^5) are easily modeled and solved without loss of accuracy. The technique has recently been expanded so as to include dielectrics and magnetic materials.Comment: 40 pages, 20 figure

Two-dimensional turbulence in magnetised plasmas

In an inhomogeneous magnetised plasma the transport of energy and particles perpendicular to the magnetic field is in general mainly caused by quasi two-dimensional turbulent fluid mixing. The physics of turbulence and structure formation is of ubiquitous importance to every magnetically confined laboratory plasma for experimental or industrial application. Specifically, high temperature plasmas for fusion energy research are also dominated by the properties of this turbulent transport. Self-organisation of turbulent vortices to mesoscopic structures like zonal flows is related to the formation of transport barriers that can significantly enhance the confinement of a fusion plasma. This subject of great importance in research is rarely touched on in introductory plasma physics or continuum dynamics courses. Here a brief tutorial on 2D fluid and plasma turbulence is presented as an introduction to the field, appropriate for inclusion in undergraduate and graduate courses.Comment: This is an author-created, un-copyedited version of an article published in European Journal of Physics. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The definitive publisher authenticated version is available online at doi: 10.1088/0143-0807/29/5/00