BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Kacper Kornet appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
GRChombo is an open-source code for performing Numerical Relativity time evolutions, built on top of the publicly available Chombo software for the solution of PDEs. Whilst GRChombo uses standard techniques in NR, it focusses on applications in theoretical physics where adaptability, both in terms of grid structure, and in terms of code modification, are key drivers.
We present Direct Numerical Simulations of decaying Magnetohydrodynamic (MHD) turbulence at low magnetic Reynolds number. The domain considered is bounded by periodic boundary conditions in the two directions perpendicular to the magnetic field and by two plane Hartmann walls in the third direction. High magnetic fields (Hartmann number of up to 896) are considered thanks to a numerical method based on a spectral code using the eigenvectors of the dissipation operator. It is found that the decay proceeds through two phases: first, energy and integral lengthscales vary rapidly during a two-dimensionalisation phase extending over about one Hartmann friction time. During this phase, the evolution of the former appears significantly more impeded by the presence of walls than that of the latter. Once the large scales are close to quasi-two dimensional, the decay results from the competition of a two-dimensional dynamics driven by dissipation in the Hartmann boundary layers and the three-dimensional dynamics of smaller scales. In the later stages of the decay, three-dimensionality subsists under the form of barrel-shaped structures. A purely quasi-two dimensional decay dominated by friction in the Hartmann layers is not reached, because of residual dissipation in the bulk. However, this dissipation is not generated by the three-dimensionality that subsists, but by residual viscous friction due to horizontal velocity gradients. Also, the energy in the velocity component aligned with the magnetic field is found to be strongly suppressed, as is transport in this direction. This results reproduces the experimental findings of Kolesnikov & Tsinober (1974).
We put forward a new type of spectral method for the direct numerical simulation of flows where anisotropy or very fine boundary layers are present. The mean idea is to take advantage of the fact that such structures are dissipative and that their presence should reduce the number of degrees of freedom of the flow, when paradoxically, their fine resolution incurs extra computational cost in most current methods. The principle of this method is to use a functional basis with elements that already include these fine structure so as to avoid these extra costs. This leads us to develop an algorithm to implement a spectral method for arbitrary functional bases, and in particular, non-orthogonal ones. We construct a basic implementation of this algorithm to simulate Magnetohydrodynamic (MHD) channel flows with an externally imposed, transverse magnetic field, where very thin boundary layers are known to develop along the channel walls. In this case, the sought functional basis can be built out of the eigenfunctions of the dissipation operator, which incorporate these boundary layers, and it turns out to be non-orthogonal. We validate this new scheme against numerical simulations of freely decaying MHD turbulence based on a Finite Volume code and it is found to provide accurate results. Its ability to fully resolve wall-bounded turbulence with a number of modes close to that required by the dynamics is demonstrated on a simple example. This opens the way to full blown simulations of MHD turbulence under very high magnetic fields, which were until now too computationally expensive, as the computational cost of the proposed method, in contrast to traditional methods, does not depend on the intensity of the magnetic field.
We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin and Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we introduce the multifluid extension of Relaxing TVD scheme and present a test case of dust migration in a two-fluid disk composed of gas and dust. We demonstrate that due to the difference in azimuthal velocities of gas and dust and the drag force acting on both components dust drifts towards maxima of gas pressure distribution.