Source author record

G. R. Sarson

G. R. Sarson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.CO astro-ph.HE astro-ph.IM physics.flu-dyn Populations and Evolution

Catalog footprint

What is connected

5works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2021arXiv

Mean fields and fluctuations in compressible magnetohydrodynamic flows

We apply Gaussian smoothing to obtain mean magnetic field, density, velocity, and magnetic and kinetic energy densities from our numerical model of the interstellar medium, based on three-dimensional magnetohydrodynamic equations in a shearing box $1 \times 1 \times 2 \, {\rm kpc}$ in size. The interstellar medium is highly compressible, as the turbulence is transonic or supersonic; it is thus an excellent context in which to explore the use of smoothing to represent physical variables in a compressible medium in terms of their mean and fluctuating parts. Unlike alternative averaging procedures, such as horizontal averaging, Gaussian smoothing retains the three-dimensional structure of the mean fields. Although Gaussian smoothing does not obey the Reynolds rules of averaging, physically meaningful and mathematically sound central statistical moments are defined as suggested by Germano (1992). We discuss methods to identify an optimal smoothing scale $\ell$ and the effects of this choice on the results. From spectral analysis of the magnetic, density and velocity fields, we find a suitable smoothing length for all three fields, of $\ell \approx 75 \, {\rm pc}$. Such a smoothing scale is likely to be sensitive to the choice of simulation parameters; this may be considered in future work, but here we just explore the methodology. We discuss the properties of third-order statistical moments in fluctuations of kinetic energy density in compressible flows, and suggest their physical interpretation. The mean magnetic field, amplified by a mean-field dynamo, significantly alters the distribution of kinetic energy in space and between scales, reducing the magnitude of kinetic energy at intermediate scales. This intermediate-scale kinetic energy is a useful diagnostic of the importance of SN-driven outflows.

preprint2020arXiv

The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained

The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid.

preprint2016arXiv

Global diffusion of cosmic rays

The propagation of charged particles, including cosmic rays, in a partially ordered magnetic field is characterized by a diffusion tensor whose components depend on the particle's Larmor radius $R_L$ and the degree of order in the magnetic field. Most studies of the particle diffusion presuppose a scale separation between the mean and random magnetic fields (e.g., there being a pronounced minimum in the magnetic power spectrum at intermediate scales). Scale separation is often a good approximation in laboratory plasmas, but not in most astrophysical environments such as the interstellar medium (ISM). Modern simulations of the ISM have numerical resolution of order 1 pc, so the Larmor radius of the cosmic rays that dominate in energy density is at least $10^{6}$ times smaller than the resolved scales. Large-scale simulations of cosmic ray propagation in the ISM thus rely on oversimplified forms of the diffusion tensor. We take the first steps towards a more realistic description of cosmic ray diffusion for such simulations, obtaining direct estimates of the diffusion tensor from test particle simulations in random magnetic fields (with the Larmor radius scale being fully resolved), for a range of particle energies corresponding to $10^{-2}\lesssim R_L/l_c \lesssim 10^{3}$, where $l_c$ is the magnetic correlation length. We obtain explicit expressions for the cosmic ray diffusion tensor for $R_L/l_c \ll 1$, that might be used in a sub-grid model of cosmic ray diffusion. The diffusion coefficients obtained are closely connected with existing transport theories that include the random walk of magnetic lines.

preprint2015arXiv

Productivity of pre-modern agriculture in the Cucuteni-Trypillia area

(Abridged) We present palaeoeconomy reconstructions for pre-modern agriculture, with the Cucuteni-Trypillia Cultural unity (5,400-2,700 BC, modern Ukraine, Moldova and Romania) as example. The starting point of our analysis is the palaeodiet structure suggested by archaeological data, stable isotope analyses of human remains, and palynology. We allow for the archeologically attested contributions of domesticated and wild animal products to the diet, develop plausible estimates of the yield of ancient cereal varieties cultivated with ancient techniques, and quantify the yield dependence on the time after initial planting and on rainfall (as a climate proxy). Our conclusions involve analysis of the labour costs of the agricultural cycle of both an individual and a farmer's family. Finally, we put our results into the context of the exploitation territory and catchment analysis. The simplest economic complex based on cereals, domestic and wild animal products, with fallow cropping, appears to be capable of supporting an isolated, relatively small farming community of 50-300 people (2-10 ha). Our results strongly suggest that dairy products played a significant role in the dietary and labour balance. The smaller settlements are typical of the earliest Trypillia A but remain predominant at the later stages. A larger settlement of several hundred people could function in isolation only with technological innovations, such as manure fertiliser and ard tillage. Very large settlements of a few hundred hectares could function only if supported by satellite farming villages and stable exchange networks. We also discuss, quantify and assess some strategies to mitigate the risks of arable agriculture associated with strong temporal fluctuations in the cereal yield, such as manure fertilisation, increased fraction of cereals in the diet combined with producing grain surplus for emergency storage.

preprint2012arXiv

The supernova-regulated ISM. II. The mean magnetic field

The origin and structure of the magnetic fields in the interstellar medium of spiral galaxies is investigated with 3D, non-ideal, compressible MHD simulations, including stratification in the galactic gravity field, differential rotation and radiative cooling. A rectangular domain, 1x1x2 kpc^{3} in size, spans both sides of the galactic mid-plane. Supernova explosions drive transonic turbulence. A seed magnetic field grows exponentially to reach a statistically steady state within 1.6 Gyr. Following Germano (1992) we use volume averaging with a Gaussian kernel to separate magnetic field into a mean field and fluctuations. Such averaging does not satisfy all Reynolds rules, yet allows a formulation of mean-field theory. The mean field thus obtained varies in both space and time. Growth rates differ for the mean-field and fluctuating field and there is clear scale separation between the two elements, whose integral scales are about 0.7 kpc and 0.3 kpc, respectively.