Researcher profile

Maciej Matyka

Maciej Matyka contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2020arXiv

Predicting Porosity, Permeability, and Tortuosity of Porous Media from Images by Deep Learning

Convolutional neural networks (CNN) are utilized to encode the relation between initial configurations of obstacles and three fundamental quantities in porous media: porosity ($φ$), permeability $k$, and tortuosity ($T$). The two-dimensional systems with obstacles are considered. The fluid flow through a porous medium is simulated with the lattice Boltzmann method. It is demonstrated that the CNNs are able to predict the porosity, permeability, and tortuosity with good accuracy. With the usage of the CNN models, the relation between $T$ and $φ$ has been reproduced and compared with the empirical estimate. The analysis has been performed for the systems with $φ\in (0.37,0.99)$ which covers five orders of magnitude span for permeability $k \in (0.78, 2.1\times 10^5)$ and tortuosity $T \in (1.03,2.74)$.

preprint2020arXiv

Pushing a soft body droplet through porous medium

We investigate the problem of transport of a single fluid droplet through a non-wettable superhydrophobic porous medium. A mechanical soft body model is developed and used to simulate the process of pushing fluid droplets through pore space of a random porous medium. Path lines of the center of mass of each droplet are used to calculate tortuosity of the media. Our results show that droplet based tortuosity increases with decreasing porosity of the porous samples. Although qualitatively this agrees with the behaviour observed for tortuosity derived from the fluid flow, the form of this relation is different.

preprint2016arXiv

Multiparticle collision dynamics in porous media

Multiparticle collision dynamics (MPCD) is a relatively new algorithm of fluid flow simulations that has been applied mostly to flows around simple objects. One might ask how it behaves in more complex flows. Therefore, we extend MPCD to simulate transport in porous media. For this, a particle-level drag force is introduced into the original algorithm. The force hinders the flow, which results in global resistance to flow and decrease of permeability. The extended algorithm is validated in the flow through a porous channel and compared with an analytical solution. Some basic properties of the solver are investigated.

preprint2015arXiv

Power exponential velocity distributions in disordered porous media

Velocity distribution functions link the micro- and macro-level theories of fluid flow through porous media. Here we study them for the fluid absolute velocity and its longitudinal and lateral components relative to the macroscopic flow direction in a model of a random porous medium. We claim that all distributions follow the power exponential law controlled by an exponent $γ$ and a shift parameter $u_0$ and examine how these parameters depend on the porosity. We find that $γ$ has a universal value $1/2$ at the percolation threshold and grows with the porosity, but never exceeds 2.

preprint2014arXiv

Compressed Multi-Row Storage Format for Sparse Matrices on Graphics Processing Units

A new format for storing sparse matrices is proposed for efficient sparse matrix-vector (SpMV) product calculation on modern graphics processing units (GPUs). This format extends the standard compressed row storage (CRS) format and can be quickly converted to and from it. Computational performance of two SpMV kernels for the new format is determined for over 130 sparse matrices on Fermi-class and Kepler-class GPUs and compared with that of five existing generic algorithms and industrial implementations, including Nvidia cuSparse CSR and HYB kernels. We found the speedup of up to $\approx 60%$ over the best of the five alternative kernels.

preprint2013arXiv

Anisotropy of flow in stochastically generated porous media

Models of porous media are often applied to relatively small systems, which leads not only to system-size-dependent results, but also to phenomena that would be absent in larger systems. Here we investigate one such finite-size effect: anisotropy of the permeability tensor. We show that a non-zero angle between the external body force and macroscopic flux vector exists in three-dimensional periodic models of sizes commonly used in computer simulations and propose a criterion, based on the system size to the grain size ratio, for this phenomenon to be relevant or negligible. The finite-size anisotropy of the porous matrix induces a pressure gradient perpendicular to the axis of a porous duct and we analyze how this effect scales with the system and grain sizes.

preprint2012arXiv

Complete PISO and SIMPLE solvers on Graphics Processing Units

We implemented the pressure-implicit with splitting of operators (PISO) and semi-implicit method for pressure-linked equations (SIMPLE) solvers of the Navier-Stokes equations on Fermi-class graphics processing units (GPUs) using the CUDA technology. We also introduced a new format of sparse matrices optimized for performing elementary CFD operations, like gradient or divergence discretization, on GPUs. We verified the validity of the implementation on several standard, steady and unsteady problems. Computational effciency of the GPU implementation was examined by comparing its double precision run times with those of essentially the same algorithms implemented in OpenFOAM. The results show that a GPU (Tesla C2070) can outperform a server-class 6-core, 12-thread CPU (Intel Xeon X5670) by a factor of 4.2.

preprint2012arXiv

How to Calculate Tortuosity Easily?

Tortuosity is one of the key parameters describing the geometry and transport properties of porous media. It is defined either as an average elongation of fluid paths or as a retardation factor that measures the resistance of a porous medium to the flow. However, in contrast to a retardation factor, an average fluid path elongation is difficult to compute numerically and, in general, is not measurable directly in experiments. We review some recent achievements in bridging the gap between the two formulations of tortuosity and discuss possible method of numerical and an experimental measurements of the tortuosity directly from the fluid velocity field.

preprint2012arXiv

Wall Orientation and Shear Stress in the Lattice Boltzmann Model

The wall shear stress is a quantity of profound importance for clinical diagnosis of artery diseases. The lattice Boltzmann is an easily parallelizable numerical method of solving the flow problems, but it suffers from errors of the velocity field near the boundaries which leads to errors in the wall shear stress and normal vectors computed from the velocity. In this work we present a simple formula to calculate the wall shear stress in the lattice Boltzmann model and propose to compute wall normals, which are necessary to compute the wall shear stress, by taking the weighted mean over boundary facets lying in a vicinity of a wall element. We carry out several tests and observe an increase of accuracy of computed normal vectors over other methods in two and three dimensions. Using the scheme we compute the wall shear stress in an inclined and bent channel fluid flow and show a minor influence of the normal on the numerical error, implying that that the main error arises due to a corrupted velocity field near the staircase boundary. Finally, we calculate the wall shear stress in the human abdominal aorta in steady conditions using our method and compare the results with a standard finite volume solver and experimental data available in the literature. Applications of our ideas in a simplified protocol for data preprocessing in medical applications are discussed.

preprint2011arXiv

Hydraulic tortuosity in arbitrary porous media flow

Tortuosity ($T$) is a parameter describing an average elongation of fluid streamlines in a porous medium as compared to free flow. In this paper several methods of calculating this quantity from lengths of individual streamlines are compared and their weak and strong features are discussed. An alternative method is proposed, which enables one to calculate $T$ directly from the fluid velocity field, without the need of determining streamlines, which greatly simplifies determination of tortuosity in complex geometries, including those found in experiments or 3D computer models. Numerical results obtained with this method suggest that (a) the hydraulic tortuosity of an isotropic fibrous medium takes on the form $T = 1 + p\sqrt{1-ϕ}$, where $ϕ$ is the porosity and $p$ is a constant and (b) the exponent controlling the divergence of $T$ with the system size at percolation threshold is related to an exponent describing the scaling of the most probable traveling length at bond percolation.

preprint2009arXiv

Finite-size anisotropy in statistically uniform porous media

Anisotropy of the permeability tensor in statistically uniform porous media of sizes used in typical computer simulations is studied. Although such systems are assumed to be isotropic by default, we show that de facto their anisotropic permeability can give rise to significant changes of transport parameters such as permeability and tortuosity. The main parameter controlling the anisotropy is $a/L$, being the ratio of the obstacle to system size. Distribution of the angle $α$ between the external force and the volumetric fluid stream is found to be approximately normal, and the standard deviation of $α$ is found to decay with the system size as $(a/L)^{d/2}$, where $d$ is the space dimensionality. These properties can be used to estimate both anisotropy-related statistical errors in large-scale simulations and the size of the representative elementary volume.

preprint2008arXiv

Tortuosity-porosity relation in the porous media flow

We study numerically the tortuosity--porosity relation in a microscopic model of a porous medium arranged as a collectin of freely overlapping squares. It is demonstrated that the finite-size effects and the discretization errors, which were ignored in previous studies, may cause significant underestimation of tortuosity. The simple tortuosity calculation method proposed here eliminates the need for using complicated, weighted averages. The numerical results presented here are in good agreement with an empirical relation between tortuosity ($T$) and porosity ($ϕ$) given by $T-1\propto \lnϕ$, that was found by others experimentally in granule packings and sediments. This relation can be also written as $T-1\propto R S/ϕ$ with $R$ and $S$ denoting the hydraulic radius of granules and the specific surface area, respectively.