Entropy Balance and Dispersive Oscillations in Lattice Boltzmann Models

We conduct an investigation into the dispersive post-shock oscillations in the entropic lattice-Boltzmann method (ELBM). To this end we use a root finding algorithm to implement the ELBM which displays fast cubic convergence and guaranties the proper sign of dissipation. The resulting simulation on the one-dimensional shock tube shows no benefit in terms of regularization from using the ELBM over the standard LBGK method. We also conduct an experiment investigating of the LBGK method using median filtering at a single point per time step. Here we observe that significant regularization can be achieved.Comment: 18 pages, 4 figures; 13/07/2009 Matlab code added to appendi

Spectral singularities and Bragg scattering in complex crystals

Spectral singularities that spoil the completeness of Bloch-Floquet states may occur in non-Hermitian Hamiltonians with complex periodic potentials. Here an equivalence is established between spectral singularities in complex crystals and secularities that arise in Bragg diffraction patterns. Signatures of spectral singularities in a scattering process with wave packets are elucidated for a PT-symmetric complex crystal.Comment: 6 pages, 5 figures, to be published in Phys. Rev.

A unified hyperbolic formulation for viscous fluids and elastoplastic solids

We discuss a unified flow theory which in a single system of hyperbolic partial differential equations (PDEs) can describe the two main branches of continuum mechanics, fluid dynamics, and solid dynamics. The fundamental difference from the classical continuum models, such as the Navier-Stokes for example, is that the finite length scale of the continuum particles is not ignored but kept in the model in order to semi-explicitly describe the essence of any flows, that is the process of continuum particles rearrangements. To allow the continuum particle rearrangements, we admit the deformability of particle which is described by the distortion field. The ability of media to flow is characterized by the strain dissipation time which is a characteristic time necessary for a continuum particle to rearrange with one of its neighboring particles. It is shown that the continuum particle length scale is intimately connected with the dissipation time. The governing equations are represented by a system of first order hyperbolic PDEs with source terms modeling the dissipation due to particle rearrangements. Numerical examples justifying the reliability of the proposed approach are demonstrated.Comment: 6 figure

A Multi-dimensional Code for Isothermal Magnetohydrodynamic Flows in Astrophysics

We present a multi-dimensional numerical code to solve isothermal magnetohydrodynamic (IMHD) equations for use in modeling astrophysical flows. First, we have built a one-dimensional code which is based on an explicit finite-difference method on an Eulerian grid, called the total variation diminishing (TVD) scheme. Recipes for building the one-dimensional IMHD code, including the normalized right and left eigenvectors of the IMHD Jacobian matrix, are presented. Then, we have extended the one-dimensional code to a multi-dimensional IMHD code through a Strang-type dimensional splitting. In the multi-dimensional code, an explicit cleaning step has been included to eliminate non-zero $\nabla\cdot B$ at every time step. To estimate the proformance of the code, one- and two-dimensional IMHD shock tube tests, and the decay test of a two-dimensional Alfv\'{e}n wave have been done. As an example of astrophysical applications, we have simulated the nonlinear evolution of the two-dimensional Parker instability under a uniform gravity.Comment: Accepted for publication in ApJ, using aaspp4.sty, 22 text pages with 10 figure

Quasinormal ringing of acoustic black holes in Laval nozzles: Numerical simulations

Quasinormal ringing of acoustic black holes in Laval nozzles is discussed. The equation for sounds in a transonic flow is written into a Schr\"{o}dinger-type equation with a potential barrier, and the quasinormal frequencies are calculated semianalytically. From the results of numerical simulations, it is shown that the quasinormal modes are actually excited when the transonic flow is formed or slightly perturbed, as well as in the real black hole case. In an actual experiment, however, the purely-outgoing boundary condition will not be satisfied at late times due to the wave reflection at the end of the apparatus, and a late-time ringing will be expressed as a superposition of "boxed" quasinormal modes. It is shown that the late-time ringing damps more slowly than the ordinary quasinormal ringing, while its central frequency is not greatly different from that of the ordinary one. Using this fact, an efficient way for experimentally detecting the quasinormal ringing of an acoustic black hole is discussed.Comment: 9 pages, 8 figures, accepted for publication in Physical Review

Relativistic MHD Simulations of Jets with Toroidal Magnetic Fields

This paper presents an application of the recent relativistic HLLC approximate Riemann solver by Mignone & Bodo to magnetized flows with vanishing normal component of the magnetic field. The numerical scheme is validated in two dimensions by investigating the propagation of axisymmetric jets with toroidal magnetic fields. The selected jet models show that the HLLC solver yields sharper resolution of contact and shear waves and better convergence properties over the traditional HLL approach.Comment: 12 pages, 5 figure

Classification of integrable two-component Hamiltonian systems of hydrodynamic type in 2+1 dimensions

Hamiltonian systems of hydrodynamic type occur in a wide range of applications including fluid dynamics, the Whitham averaging procedure and the theory of Frobenius manifolds. In 1+1 dimensions, the requirement of the integrability of such systems by the generalised hodograph transform implies that integrable Hamiltonians depend on a certain number of arbitrary functions of two variables. On the contrary, in 2+1 dimensions the requirement of the integrability by the method of hydrodynamic reductions, which is a natural analogue of the generalised hodograph transform in higher dimensions, leads to finite-dimensional moduli spaces of integrable Hamiltonians. In this paper we classify integrable two-component Hamiltonian systems of hydrodynamic type for all existing classes of differential-geometric Poisson brackets in 2D, establishing a parametrisation of integrable Hamiltonians via elliptic/hypergeometric functions. Our approach is based on the Godunov-type representation of Hamiltonian systems, and utilises a novel construction of Godunov's systems in terms of generalised hypergeometric functions.Comment: Latex, 34 page

Evolutions of Magnetized and Rotating Neutron Stars

We study the evolution of magnetized and rigidly rotating neutron stars within a fully general relativistic implementation of ideal magnetohydrodynamics with no assumed symmetries in three spatial dimensions. The stars are modeled as rotating, magnetized polytropic stars and we examine diverse scenarios to study their dynamics and stability properties. In particular we concentrate on the stability of the stars and possible critical behavior. In addition to their intrinsic physical significance, we use these evolutions as further tests of our implementation which incorporates new developments to handle magnetized systems.Comment: 12 pages, 8 figure

COSMOS: A Hybrid N-Body/Hydrodynamics Code for Cosmological Problems

We describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and individual galaxies. The code, named COSMOS, possesses several new features which distinguish it from other PM-PPM codes. In particular, to solve the Poisson equation we have written a new multigrid solver which can determine the gravitational potential of isolated matter distributions and which properly takes into account the finite-volume discretization required by PPM. All components of the code are constructed to work with a nonuniform mesh, preserving second-order spatial differences. The PPM code uses vacuum boundary conditions for isolated problems, preventing inflows when appropriate. The PM code uses a second-order variable-timestep time integration scheme. Radiative cooling and cosmological expansion terms are included. COSMOS has been implemented for parallel computers using the Parallel Virtual Machine (PVM) library, and it features a modular design which simplifies the addition of new physics and the configuration of the code for different types of problems. We discuss the equations solved by COSMOS and describe the algorithms used, with emphasis on these features. We also discuss the results of tests we have performed to establish that COSMOS works and to determine its range of validity.Comment: 43 pages, 14 figures, submitted to ApJS and revised according to referee's comment

Consistent thermodynamic derivative estimates for tabular equations of state

Numerical simulations of compressible fluid flows require an equation of state (EOS) to relate the thermodynamic variables of density, internal energy, temperature, and pressure. A valid EOS must satisfy the thermodynamic conditions of consistency (derivation from a free energy) and stability (positive sound speed squared). When phase transitions are significant, the EOS is complicated and can only be specified in a table. For tabular EOS's such as SESAME from Los Alamos National Laboratory, the consistency and stability conditions take the form of a differential equation relating the derivatives of pressure and energy as functions of temperature and density, along with positivity constraints. Typical software interfaces to such tables based on polynomial or rational interpolants compute derivatives of pressure and energy and may enforce the stability conditions, but do not enforce the consistency condition and its derivatives. We describe a new type of table interface based on a constrained local least squares regression technique. It is applied to several SESAME EOS's showing how the consistency condition can be satisfied to round-off while computing first and second derivatives with demonstrated second-order convergence. An improvement of 14 orders of magnitude over conventional derivatives is demonstrated, although the new method is apparently two orders of magnitude slower, due to the fact that every evaluation requires solving an 11-dimensional nonlinear system.Comment: 29 pages, 9 figures, 16 references, submitted to Phys Rev