BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Zbigniew Koza

Zbigniew Koza appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dyncond-mat.stat-mechphysics.comp-phDistributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

16 published item(s)

preprint2022arXiv

Percolation in a triangle on a square lattice

Percolation on a plane is usually associated with clusters spanning two opposite sides of a rectangular system. Here we investigate three-leg clusters generated on a square lattice and spanning the three sides of equilateral triangles. If the position and orientation of the triangles relative to the lattice are uniformly randomized, one obtains an efficient method of determining the percolation threshold, on par with the most advanced Monte Carlo methods developed for the rectangular geometry. The universal crossing probability for three-leg clusters is geometry-independent, which opens a way for further improvements of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Jammed systems of oriented dimers always percolate on hypercubic lattices

Random sequential adsorption (RSA) is a standard method of modeling adsorption of large molecules at the liquid-solid interface. Here we consider jammed states of the RSA process of nonoverlapping dimers (objects occupying two nearest-neighbor lattice sites) in a hypercubic lattice of arbitrary space dimension $D \ge 2$. We show that each dimer in such a state belongs to a percolating cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From discrete to continuous percolation in dimensions 3 to 7

We propose a method of studying the continuous percolation of aligned objects as a limit of a corresponding discrete model. We show that the convergence of a discrete model to its continuous limit is controlled by a power-law dependency with a universal exponent $θ= 3/2$. This allows us to estimate the continuous percolation thresholds in a model of aligned hypercubes in dimensions $d = 3,\ldots,7$ with accuracy far better than that attained using any other method before. We also report improved values of the correlation length critical exponent $ν$ in dimensions $d = 4,5$ and the values of several universal wrapping probabilities for $d=4,\ldots,7$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Percolation of overlapping squares or cubes on a lattice

Porous media are often modelled as systems of overlapping obstacles, which leads to the problem of two percolation thresholds in such systems, one for the porous matrix and the other one for the void space. Here we investigate these percolation thresholds in the model of overlapping squares or cubes of linear size $k>1$ randomly distributed on a regular lattice. We find that the percolation threshold of obstacles is a nonmonotonic function of $k$, whereas the percolation threshold of the void space is well approximated by a function linear in $1/k$. We propose a generalization of the excluded volume approximation to discrete systems and use it to investigate the transition between continuous and discrete percolation, finding a remarkable agreement between the theory and numerical results. We argue that the continuous percolation threshold of aligned squares on a plane is the same for the solid and void phases and estimate the continuous percolation threshold of the void space around aligned cubes in a 3D space as 0.036(1). We also discuss the connection of the model to the standard site percolation with complex neighborhood.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The dual model for an Ising model with nearest and next-nearest neighbors

We construct and analyse a dual model to the Ising model with the nearest and next-nearest neighbors on the rectangular lattice (NNNI model). The Hamiltonian of the dual model turns out to contain two- and four-spin interactions. The free fermion approximation suggests that an increase in the critical temperature of the dual model caused by the four-spin interactions is limited to a finite range.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accelerating cellular automata simulations using AVX and CUDA

We investigated various methods of parallelization of the Frish-Hasslacher-Pomeau (FHP) cellular automata algorithm for modeling fluid flow. These methods include SSE, AVX, and POSIX Threads for central processing units (CPUs) and CUDA for graphics processing units (GPUs). We present implementation details of the FHP algorithm based on AVX/SSE and CUDA technologies. We found that (a) using AVX or SSE is necessary to fully utilize the potential of modern CPUs; (b) CPUs and GPUs are comparable in terms of computational and economic efficiency only if the CPU code uses AVX or SSE instructions; (c) AVX does not offer any substantial improvement relative to SSE.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The asymptotic behaviour of the initially separated A + B(static) -> 0 reaction-diffusion systems

We examine the long-time behaviour of A+B \to 0 reaction-diffusion systems with initially separated species A and B. All of our analysis is carried out for arbitrary (positive) values of the diffusion constant D_A of particles A and initial concentrations a_0 and b_0 of A's and B's. We derive general formulae for the location of the reaction zone centre, the total reaction rate, and the concentration profile of species A outside the reaction zone. The general properties of the reaction zone are studied with a help of the scaling ansatz. Using the mean-field approximation we find the functional forms of `tails' of the reaction rate R and the dependence of the width of the reaction zone on the external parameters of the system. We also study the change in the kinetics of the system with D_B > 0 in the limit D_B \to 0. Our results are supported by numerical solutions of the mean-field reaction-diffusion equation.

0 citations0 reviewsNo code linkedOpen access