We present a new method for the numerical solution of the radiative-transfer equation (RTE) in multidimensional scenarios commonly encountered in computational astrophysics. The method is based on the direct solution of the Boltzmann equation via an extension of the Lattice Boltzmann (LB) methods and allows to model the evolution of the radiation field as it interacts with a background fluid, via absorption, emission, and scattering. As a first application of this method, we restrict our attention to a frequency independent (``grey'') formulation within a special-relativistic framework, which can be employed also for classical computational astrophysics. For a number of standard tests that consider the performance of the method in optically thin, optically thick and intermediate regimes with a static fluid, we show the ability of the LB method to produce accurate and convergent results matching the analytic solutions. We also contrast the LB method with commonly employed moment-based schemes for the solution of the RTE, such as the M1 scheme. In this way, we are able to highlight that the LB method provides the correct solution for both non-trivial free-streaming scenarios and the intermediate optical-depth regime, for which the M1 method either fails or provides inaccurate solutions. When coupling to a dynamical fluid, on the other hand, we present the first self-consistent solution of the RTE with LB methods within a relativistic-hydrodynamic scenario. Finally, we show that besides providing more accurate results in all regimes, the LB method features smaller or comparable computational costs compared to the M1 scheme. We conclude that LB methods represent a competitive and promising avenue to the solution of radiative transport, one of the most common and yet important problems in computational astrophysics.

We formulate a smoothed-particle hydrodynamics numerical method, traditionally used for the Euler equations for fluid dynamics in the context of astrophysical simulations, to solve the non-linear Schrodinger equation in the Madelung formulation. The probability density of the wavefunction is discretized into moving particles, whose properties are smoothed by a kernel function. The traditional fluid pressure is replaced by a quantum pressure tensor, for which a novel, robust discretization is found. We demonstrate our numerical method on a variety of numerical test problems involving the simple harmonic oscillator, Bose-Einstein condensates, collapsing singularities, and dark matter halos governed by the Gross-Pitaevskii-Poisson equation. Our method is conservative, applicable to unbounded domains, and is automatically adaptive in its resolution, making it well suited to study problems with collapsing solutions.

We revisit the integer lattice (IL) method to numerically solve the Vlasov-Poisson equations, and show that a slight variant of the method is a very easy, viable, and efficient numerical approach to study the dynamics of self-gravitating, collisionless systems. The distribution function lives in a discretized lattice phase-space, and each time-step in the simulation corresponds to a simple permutation of the lattice sites. Hence, the method is Lagrangian, conservative, and fully time-reversible. IL complements other existing methods, such as N-body/particle mesh (computationally efficient, but affected by Monte-Carlo sampling noise and two-body relaxation) and finite volume (FV) direct integration schemes (expensive, accurate but diffusive). We also present improvements to the FV scheme, using a moving mesh approach inspired by IL, to reduce numerical diffusion and the time-step criterion. Being a direct integration scheme like FV, IL is memory limited (memory requirement for a full 3D problem scales as N^6, where N is the resolution per linear phase-space dimension). However, we describe a new technique for achieving N^4 scaling. The method offers promise for investigating the full 6D phase-space of collisionless systems of stars and dark matter.

Fluid dynamics in intrinsically curved geometries is encountered in many physical systems in nature, ranging from microscopic bio-membranes all the way up to general relativity at cosmological scales. Despite the diversity of applications, all of these systems share a common feature: the free motion of particles is affected by inertial forces originating from the curvature of the embedding space. Here we reveal a fundamental process underlying fluid dynamics in curved space: the free motion of fluids, in the complete absence of solid walls or obstacles, exhibits loss of energy due exclusively to the intrinsic curvature of space. We find that local sources of curvature generate viscous stresses as a result of the inertial forces. The curvature-induced viscous forces are shown to cause hitherto unnoticed and yet appreciable energy dissipation, which might play a significant role for a variety of physical systems involving fluid dynamics in curved spaces.

We investigate Poiseuille channel flow through intrinsically curved media, equipped with localized metric perturbations. To this end, we study the flux of a fluid driven through the curved channel in dependence of the spatial deformation, characterized by the parameters of the metric perturbations (amplitude, range and density). We find that the flux depends only on a specific combination of parameters, which we identify as the average metric perturbation, and derive a universal flux law for the Poiseuille flow. For the purpose of this study, we have improved and validated our recently developed lattice Boltzmann model in curved space by considerably reducing discrete lattice effects.

In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfv\'en waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to $\sigma^{-\frac{1}{2}}$, $\sigma$ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.

We have found that the relation between the flow through campylotic (generically curved) media, consisting of randomly located curvature perturbations, and the average Ricci scalar of the system exhibits two distinct functional expressions (hysteresis), depending on whether the typical spatial extent of the curvature perturbation lies above or below the critical value maximizing the overall Ricci curvature. Furthermore, the flow through such systems as a function of the number of curvature perturbations presents a sublinear behavior for large concentrations due to the interference between curvature perturbations that, consequently, produces a less curved space. For the purpose of this study, we have developed and validated a lattice kinetic model capable of describing fluid flow in arbitrarily curved manifolds, which allows to deal with highly complex spaces in a very compact and efficient way.

We provide numerical evidence that the Richtmyer-Meshkov (RM) instability contributes to the cooling of a relativistic fluid. Due to the presence of jet particles traveling throughout the medium, shock waves are generated in the form of Mach cones. The interaction of multiple shock waves can trigger the RM instability, and we have found that this process leads to a down-cooling of the relativistic fluid. To confirm the cooling effect of the instability, shock tube Richtmyer-Meshkov instability simulations are performed. Additionally, in order to provide an experimental observable of the RM instability resulting from the Mach cone interaction, we measure the two particle correlation function and highlight the effects of the interaction. The simulations have been performed with an improved version of the relativistic lattice Boltzmann model, including general equations of state and external forces.

We present a new approach to find accurate solutions to the Poisson equation, as obtained from the steady-state limit of a diffusion equation with strong source terms. For this purpose, we start from Boltzmann's kinetic theory and investigate the influence of higher order terms on the resulting macroscopic equations. By performing an appropriate expansion of the equilibrium distribution, we provide a method to remove the unnecessary terms up to a desired order and show that it is possible to find, with high level of accuracy, the steady-state solution of the diffusion equation for sizeable Knudsen numbers. In order to test our kinetic approach, we discretise the Boltzmann equation and solve the Poisson equation, spending up to six order of magnitude less computational time for a given precision than standard lattice Boltzmann methods.

Theoretical and numerical analysis of the relativistic effects on the Richtmyer-Meshkov (RM) instability reveals new and potentially very useful effects. We find that, in contrast with the non- relativistic case, the growth rate of the RM instability depends strongly on the equation of state of the fluid, opening up the possibility to infer equations of state from experimental observations of the RM instability. As opposed to the non-relativistic case, we also discover that, above a critical value of the fluid velocity, the growth rate of the instability counter-intuitively decreases due to the Lorentz's factor, and vanishes in the ultrarelativistic limit, as the speed of the particles approaches the speed of light. Both effects might prove very useful for leading-edge applications, such as the study of the equation of state of quark-gluon matter, and the design of fast ignition inertial confinement fusion (ICF) schemes. We perform a linear stability analysis to characterize the instability, for an arbitrary equation of state, and implement numerical simulations to study the instability in the non-linear regime, using the equation of state of an ideal gas. Furthermore, based on the numerical results, we propose a general expression that characterizes the long term evolution of the instability.

We develop a relativistic lattice Boltzmann model capable of describing relativistic fluid dynamics at ultra-high velocities, with Lorentz factors up to $\gamma \sim 10$. To this purpose, we first build a new lattice kinetic scheme by expanding the Maxwell-J\"uttner distribution function in an orthogonal basis of polynomials and applying an appropriate quadrature, providing the discrete versions of the relativistic Boltzmann equation and the equilibrium distribution. To achieve ultra-high velocities, we include a flux limiter scheme, and introduce the bulk viscosity by a suitable extension of the discrete relativistic Boltzmann equation. The model is validated by performing simulations of shock waves in viscous quark-gluon plasmas and comparing with existing models, finding very good agreement. To the best of our knowledge, we for the first time successfully simulate viscous shock waves in the highly relativistic regime. Moreover, we show that our model can also be used for near-inviscid flows even at very high velocities. Finally, as an astrophysical application, we simulate a relativistic shock wave, generated by, say, a supernova explosion, colliding with a massive interstellar cloud, e.g. molecular gas.

We present a new Lattice Boltzmann (LB) formulation to solve the Maxwell equations for electromagnetic (EM) waves propagating in a heterogeneous medium. By using a pseudo-vector discrete Boltzmann distribution, the scheme is shown to reproduce the continuum Maxwell equations. The technique compares well with a pseudo-spectral method at solving for two-dimensional wave propagation in a heterogeneous medium, which by design contains substantial contrasts in the refractive index. The extension to three dimensions follows naturally and, owing to the recognized efficiency of LB schemes for parallel computation in irregular geometries, it gives a powerful method to numerically simulate a wide range of problems involving EM wave propagation in complex media.

A detailed derivation of the Lattice Boltzmann (LB) scheme for relativistic fluids recently proposed in Ref. [1], is presented. The method is numerically validated and applied to the case of two quite different relativistic fluid dynamic problems, namely shock-wave propagation in quark-gluon plasmas and the impact of a supernova blast-wave on massive interstellar clouds. Close to second order convergence with the grid resolution, as well as linear dependence of computational time on the number of grid points and time-steps, are reported.