Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDistributed order Hausdorff derivative diffusion model to characterize non-Fickian diffusion in porous media
339
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
This paper presents a time-space Hausdorff derivative model for depicting solute transport in aquifers or water flow in heterogeneous porous media. In this model, the time and space Hausdorff derivatives are defined on non-Euclidean fractal metrics with power law scaling transform which, respectively, connect the temporal and spatial complexity during transport. The Hausdorff derivative model can be transformed to an advection-dispersion equation with time- and space-dependent dispersion and convection coefficients. This model is a fractal partial differential equation (PDE) defined on a fractal space and differs from the fractional PDE which is derived for non-local transport of particles on a non-fractal Euclidean space. As an example of applications of this model, an explicit solution with a constant diffusion coefficient and flow velocity subject to an instantaneous source is derived and fitted to the breakthrough curves of tritium as a tracer in porous media. These results are compared with those of a scale-dependent dispersion model and a time-scale dependent dispersion model. Overall, it is found that the fractal PDE based on the Hausdorff derivatives better captures the early arrival and heavy tail in the scaled breakthrough curves for variable transport distances. The estimated parameters in the fractal Hausrdorff model represent clear physical mechanisms such as linear relationships between the orders of Hausdorff derivatives and the transport distance. The mathematical formulation is applicable to both solute transport and water flow in porous media.
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.
It is of great concern to produce numerically efficient methods for moisture diffusion through porous media, capable of accurately calculate moisture distribution with a reduced computational effort. In this way, model reduction methods are promising approaches to bring a solution to this issue since they do not degrade the physical model and provide a significant reduction of computational cost. Therefore, this article explores in details the capabilities of two model-reduction techniques - the Spectral Reduced-Order Model (Spectral-ROM) and the Proper Generalised Decomposition (PGD) - to numerically solve moisture diffusive transfer through porous materials. Both approaches are applied to three different problems to provide clear examples of the construction and use of these reduced-order models. The methodology of both approaches is explained extensively so that the article can be used as a numerical benchmark by anyone interested in building a reduced-order model for diffusion problems in porous materials. Linear and non-linear unsteady behaviors of unidimensional moisture diffusion are investigated. The last case focuses on solving a parametric problem in which the solution depends on space, time and the diffusivity properties. Results have highlighted that both methods provide accurate solutions and enable to reduce significantly the order of the model around ten times lower than the large original model. It also allows an efficient computation of the physical phenomena with an error lower than 10^{-2} when compared to a reference solution.
In this paper we investigate the solution of generalized distributed order diffusion equations with composite time fractional derivative by using the Fourier-Laplace transform method. We represent solutions in terms of infinite series in Fox $H$-functions. The fractional and second moments are derived by using Mittag-Leffler functions. We observe decelerating anomalous subdiffusion in case of two composite time fractional derivatives. Generalized uniformly distributed order diffusion equation, as a model for strong anomalous behavior, is analyzed by using Tauberian theorem. Some previously obtained results are special cases of those presented in this paper.
We investigate the elastoviscoplastic flow through porous media by numerical simulations. We solve the Navier-Stokes equations combined with the elastoviscoplastic model proposed by Saramito for the stress tensor evolution. In this model, the material behaves as a viscoelastic solid when unyielded, and as a viscoelastic Oldroyd-B fluid for stresses higher than the yield stress. The porous media is made of a symmetric array of cylinders, and we solve the flow in one periodic cell. We find that the solution is time-dependent even at low Reynolds numbers as we observe oscillations in time of the unyielded region especially at high Bingham numbers. The volume of the unyielded region slightly decreases with the Reynolds number and strongly increases with the Bingham number; up to 70% of the total volume is unyielded for the highest Bingham numbers considered here. The flow is mainly shear dominated in the yielded region, while shear and elongational flow are equally distributed in the unyielded region. We compute the relation between the pressure drop and the flow rate in the porous medium and present an empirical closure as function of the Bingham and Reynolds numbers. The apparent permeability, normalized with the case of Newtonian fluids, is shown to be greater than 1 at low Bingham numbers, corresponding to lower pressure drops due to the flow elasticity, and smaller than 1 for high Bingham numbers, indicating larger dissipation in the flow owing to the presence of the yielded regions. Finally we investigate the effect of the Weissenberg number on the distribution of the unyielded regions and on the pressure gradient.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
