14,896 research outputs found
A new numerical mesoscopic scale one-domain approach solver for free fluid/porous medium interaction
A new numerical continuum one-domain approach (ODA) solver is presented for the simulation of the transfer processes between a free fluid and a porous medium. The solver is developed in the \textit{mesoscopic} scale framework, where a continuous variation of the physical parameters of the porous medium (e.g., porosity and permeability) is assumed. The Navier--Stokes--Brinkman equations are solved along with the continuity equation, under the hypothesis of incompressible fluid. The porous medium is assumed to be fully saturated and can potentially be anisotropic. The domain is discretized with unstructured meshes allowing local refinements. A fractional time step procedure is applied, where one predictor and two corrector steps are solved within each time iteration. The predictor step is solved in the framework of a marching in space and time procedure, with some important numerical advantages. The two corrector steps require the solution of large linear systems, whose matrices are sparse, symmetric and positive definite, with -matrix property over Delaunay-meshes. A fast and efficient solution is obtained using a preconditioned conjugate gradient method. The discretization adopted for the two corrector steps can be regarded as a Two-Point-Flux-Approximation (TPFA) scheme, which, unlike the standard TPFA schemes, does not require the grid mesh to be K-orthogonal, (with {K the anisotropy tensor). As demonstrated with the provided test cases, the proposed scheme correctly retains the anisotropy effects within the porous medium. Furthermore, it overcomes the restrictions of existing mesoscopic scale one-domain approaches proposed in the literature
Biot-JKD model: simulation of 1D transient poroelastic waves with fractional derivatives
A time-domain numerical modeling of Biot poroelastic waves is presented. The
viscous dissipation occurring in the pores is described using the dynamic
permeability model developed by Johnson-Koplik-Dashen (JKD). Some of the
coefficients in the Biot-JKD model are proportional to the square root of the
frequency: in the time-domain, these coefficients introduce order 1/2 shifted
fractional derivatives involving a convolution product. Based on a diffusive
representation, the convolution kernel is replaced by a finite number of memory
variables that satisfy local-in-time ordinary differential equations. Thanks to
the dispersion relation, the coefficients in the diffusive representation are
obtained by performing an optimization procedure in the frequency range of
interest. A splitting strategy is then applied numerically: the propagative
part of Biot-JKD equations is discretized using a fourth-order ADER scheme on a
Cartesian grid, whereas the diffusive part is solved exactly. Comparisons with
analytical solutions show the efficiency and the accuracy of this approach.Comment: arXiv admin note: substantial text overlap with arXiv:1210.036
Finite element modeling of free surface flow in variable porosity media
The aim of the present work is to present an overview of some numerical procedures for the simulation of free surface flows within a porous structure. A particular algorithm developed by the authors for solving this type of problems is presented. A modified form of the classical Navier–Stokes equations is proposed, with the principal aim of simulating in a unified way the seepage flow inside rockfill-like porous material and the free surface flow in the clear fluid region. The problem is solved using a semi-explicit stabilized fractional step algorithm where velocity is calculated using a 4th order Runge–Kutta scheme. The numerical formulation is developed in an Eulerian framework using a level set technique to track the evolution of the free surface. An edge-based data structure is employed to allow an easy OpenMP parallelization of the resulting finite element code. The numerical model is validated against laboratory experiments on small scale rockfill dams and is compared with other existing methods for solving similar problems.Peer ReviewedPostprint (author’s final draft
Numerical Methods for the Fractional Laplacian: a Finite Difference-quadrature Approach
The fractional Laplacian is a non-local operator which
depends on the parameter and recovers the usual Laplacian as . A numerical method for the fractional Laplacian is proposed, based on
the singular integral representation for the operator. The method combines
finite difference with numerical quadrature, to obtain a discrete convolution
operator with positive weights. The accuracy of the method is shown to be
. Convergence of the method is proven. The treatment of far
field boundary conditions using an asymptotic approximation to the integral is
used to obtain an accurate method. Numerical experiments on known exact
solutions validate the predicted convergence rates. Computational examples
include exponentially and algebraically decaying solution with varying
regularity. The generalization to nonlinear equations involving the operator is
discussed: the obstacle problem for the fractional Laplacian is computed.Comment: 29 pages, 9 figure
- …