No Arabic abstract
Fluid approximations to cosmic ray (CR) transport are often preferred to kinetic descriptions in studies of the dynamics of the interstellar medium (ISM) of galaxies, because they allow simpler analytical and numerical treatments. Magnetohydrodynamic (MHD) simulations of the ISM usually incorporate CR dynamics as an advection-diffusion equation for CR energy density, with anisotropic, magnetic field-aligned diffusion with the diffusive flux assumed to obey Ficks law. We compare test-particle and fluid simulations of CRs in a random magnetic field. We demonstrate that a non-Fickian prescription of CR diffusion, which corresponds to the telegraph equation for the CR energy density, can be easily calibrated to match the test particle simulations with great accuracy. In particular, we consider a random magnetic field in the fluid simulation that has a lower spatial resolution than that used in the particle simulation to demonstrate that an appropriate choice of the diffusion tensor can account effectively for the unresolved (subgrid) scales of the magnetic field. We show that the characteristic time which appears in the telegraph equation can be physically interpreted as the time required for the particles to reach a diffusive regime and we stress that the Fickian description of the CR fluid is unable to describe complex boundary or initial conditions for the CR energy flux.
In the absence of fractures, methane bubbles in deep-water sediments can be immovably trapped within a porous matrix by surface tension. The dominant mechanism of transfer of gas mass therefore becomes the diffusion of gas molecules through porewater. The accurate description of this process requires non-Fickian diffusion to be accounted for, including both thermodiffusion and gravitational action. We evaluate the diffusive flux of aqueous methane considering non-Fickian diffusion and predict the existence of extensive bubble mass accumulation zones within deep-water sediments. The limitation on the hydrate deposit capacity is revealed; too weak deposits cannot reach the base of the hydrate stability zone and form any bubbly horizon.
Many theoretical and experimental results show that solute transport in heterogeneous porous media exhibits multi-scaling behaviors. To describe such non-Fickian diffusions, this work provides a distributed order Hausdorff diffusion model to describe the tracer transport in porous media. This model is proved to be equivalent with the diffusion equation model with a nonlinear time dependent diffusion coefficient. In conjunction with the structural derivative, its mean squared displacement (MSD) of the tracer particles is explicitly derived as a dilogarithm function when the weight function of the order distribution is a linear function of the time derivative order. This model can capture both accelerating and decelerating anomalous and ultraslow diffusions by varying the weight parameter c. In this study, the tracer transport in water-filled pore spaces of two-dimensional Euclidean is demonstrated as a decelerating sub-diffusion, and can well be described by the distributed order Hausdorff diffusion model with c = 1.73. While the Hausdorff diffusion model can accurately fit the sub-diffusion experimental data of the tracer transport in the pore-solid prefractal porous media.
Cosmic ray propagation is diffusive because of pitch angle scattering by waves. We demonstrate that if the high-amplitude magnetohydrodynamic turbulence with $tilde B/langle Brangle sim 1$ is present on top of the mean field gradient, the diffusion becomes asymmetric. As an example, we consider the vertical transport of cosmic rays in our Galaxy propagating away from a point-like source. We solve this diffusion problem analytically using a one-dimensional Markov chain analysis. We obtained that the cosmic ray density markedly differs from the standard diffusion prediction and has a sizable effect on their distribution throughout the galaxy. The equation for the continuous limit is also derived, which shows limitations of the convection-diffusion equation.
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 particles 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.
Over the last years several experimental and theoretical studies of diffusion kinetics on the nanoscale have shown that the time evolution differs from the classical Fickian law (kc=0.5). However, all work was based on crystalline samples or models, so far. In this letter, we report on the diffusion kinetics of a thin amorphous-Si layer into amorphous-Ge to account for the rising importance of amorphous materials in nanodevices. Employing surface sensitive technics, the initial kc was found at 0.7+-0.1. Moreover, after some monolayers of Si dissolved into the Ge, kc changes to the generally expected classical Fickian law with kc=0.5.