In magnetized plasmas gravitational and electromagnetic waves may interact coherently and exchange energy between themselves and with plasma flows. We derive the wave interaction equations for these processes in the case of waves propagating perpendicular or parallel to the plasma background magnetic field. In the latter case, the electromagnetic waves are taken to be circularly polarized waves of arbitrary amplitude. We allow for a background drift flow of the plasma components which increases the number of possible evolution scenarios. The interaction equations are solved analytically and the characteristic time scales for conversion between gravitational and electromagnetic waves are found. In particular, it is shown that in the presence of a drift flow there are explosive instabilities resulting in the generation of gravitational and electromagnetic waves. Conversely, we show that energetic waves can interact to accelerate particles and thereby produce a drift flow. The relevance of these results for astrophysical and cosmological plasmas is discussed.

SUMMARYThe smooth and nonsmooth approaches to the discrete element method (DEM) are examined from a computational perspective. The main difference can be understood as using explicit versus implicit time integration. A formula is obtained for estimating the computational effort depending on error tolerance, system geometric shape and size, and on the dynamic state. For the nonsmooth DEM (NDEM), a regularized version mapping to the Hertz contact law is presented. This method has the conventional nonsmooth and smooth DEM as special cases depending on size of time step and value of regularization. The use of the projected Gauss-Seidel solver for NDEM simulation is studied on a range of test systems. The following characteristics are found. First, the smooth DEM is computationally more efficient for soft materials, wide and tall systems, and with increasing flow rate. Secondly, the NDEM is more beneficial for stiff materials, shallow systems, static or slow flow, and with increasing error tolerance. Furthermore, it is found that just as pressure saturates with depth in a granular column, due to force arching, also the required number of iterations saturates and become independent of system size. This effect make the projected Gauss-Seidel solver scale much better than previously thought.

The exact 1+3 covariant dynamical fluid equations for a multi-component plasma, together with Maxwell's equations are presented in such a way as to make them suitable for a gauge-invariant analysis of linear density and velocity perturbations of the Friedmann-Robertson-Walker model. In the case where the matter is described by a two component plasma where thermal effects are neglected, a mode representing high-frequency plasma oscillations is found in addition to the standard growing and decaying gravitational instability picture. Further applications of these equations are also discussed.

We consider propagation of gravitational radiation in a magnetized multicomponent plasma. It is shown that large density perturbations can be generated, even for small deviations from flat space, provided the cyclotron frequency is much larger than the plasma frequency. Furthermore, the induced density gradients can generate frequency conversion of electromagnetic radiation, which may give rise to indirect observational effect of the gravitational waves.

Abstract-We present a fluid simulation method based on Smoothed Particle Hydrodynamics (SPH) in which incompressibility and boundary conditions are enforced using holonomic kinematic constraints on the density. This formulation enables systematic multiphysics integration in which interactions are modeled via similar constraints between the fluid pseudo-particles and impenetrable surfaces of other bodies. These conditions embody Archimede's principle for solids and thus buoyancy results as a direct consequence. We use a variational time stepping scheme suitable for general constrained multibody systems we call SPOOK. Each step requires the solution of only one Mixed Linear Complementarity Problem (MLCP) with very few inequalities, corresponding to solid boundary conditions. We solve this MLCP with a fast iterative method. Overall stability is vastly improved in comparison to the unconstrained version of SPH, and this allows much larger time steps, and an increase in overall performance by two orders of magnitude. Proof of concept is given for computer graphics applications and interactive simulations.