اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدGeneralized semi-analytical solution for coupled multispecies advection-dispersion equations in multilayer porous media
176
0
0.0
(
0
)
تأليف
Elliot J. Carr
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
Multispecies contaminant transport in the Earths subsurface is commonly modelled using advection-dispersion equations coupled via first-order reactions. Analytical and semi-analytical solutions for such problems are highly sought after but currently limited to either one species, homogeneous media, certain reaction networks, specific boundary conditions or a combination thereof. In this paper, we develop a semi-analytical solution for the case of a heterogeneous layered medium and a general first-order reaction network. Our approach combines a transformation method to decouple the multispecies equations with a recently developed semi-analytical solution for the single-species advection-dispersion-reaction equation in layered media. The generalized solution is valid for arbitrary numbers of species and layers, general Robin-type conditions at the inlet and outlet and accommodates both distinct retardation factors across layers or distinct retardation factors across species. Four test cases are presented to demonstrate the solution approach with the reported results in agreement with previously published results and numerical results obtained via finite volume discretisation. MATLAB code implementing the generalized semi-analytical solution is made available.
قيم البحث
اقرأ أيضاً
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.
Transport of viscous fluid through porous media is a direct consequence of the pore structure. Here we investigate transport through a specific class of two-dimensional porous geometries, namely those formed by fluid-mechanical erosion. We investigat
e the tortuosity and dispersion by analyzing the first two statistical moments of tracer trajectories. For most initial configurations, tortuosity decreases in time as a result of erosion increasing the porosity. However, we find that tortuosity can also increase transiently in certain cases. The porosity-tortuosity relationships that result from our simulations are compared with models available in the literature. Asymptotic dispersion rates are also strongly affected by the erosion process, as well as by the number and distribution of the eroding bodies. Finally, we analyze the pore size distribution of an eroding geometry. The simulations are performed by combining a high-fidelity boundary integral equation solver for the fluid equations, a second-order stable time stepping method to simulate erosion, and new numerical methods to stably and accurately resolve nearly-touching eroded bodies and particle trajectories near the eroding bodies.
In this paper, authors focus effort on improving the conventional discrete velocity method (DVM) into a multiscale scheme in finite volume framework for gas flow in all flow regimes. Unlike the typical multiscale kinetic methods unified gas-kinetic s
cheme (UGKS) and discrete unified gas-kinetic scheme (DUGKS), which concentrate on the evolution of the distribution function at the cell interface, in the present scheme the flux for macroscopic variables is split into the equilibrium part and the nonequilibrium part, and the nonequilibrium flux is calculated by integrating the discrete distribution function at the cell center, which overcomes the excess numerical dissipation of the conventional DVM in the continuum flow regime. Afterwards, the macroscopic variables are finally updated by simply integrating the discrete distribution function at the cell center, or by a blend of the increments based on the macroscopic and the microscopic systems, and the multiscale property is achieved. Several test cases, involving unsteady, steady, high speed, low speed gas flows in all flow regimes, have been performed, demonstrating the good performance of the multiscale DVM from free molecule to continuum Navier-Stokes solutions and the multiscale property of the scheme is proved.
The isentropic vortex problem is frequently solved to test the accuracy of numerical methods and verify corresponding code. Unfortunately, its existing solution was derived in the relativistic magnetohydrodynamics by numerically solving an ordinary d
ifferential equation. This note provides an analytical solution of the 2D isentropic vortex problem with explicit algebraic expressions in the special relativistic hydrodynamics and magnetohydrodynamics and extends it to the 3D case.
We are interested in the modeling of wave propagation in an infinite bilayered acoustic/poroelastic media. We consider the biphasic Biots model in the poroelastic layer. The first part is devoted to the calculation of analytical solution in two dimen
sions, thanks to Cagniard de Hoop method. In this second part we consider the 3D case.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
