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

Piotr Szymczak

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

ResearcherUnclaimed source record
physics.geo-phphysics.flu-dyncond-mat.softnlin.PScond-mat.stat-mechastro-ph.EPBiological Physicscond-mat.mtrl-scimath-phmath.GTmath.MGmath.MPPopulations and Evolution

Catalog footprint

What is connected

12works
13topics
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

12 published item(s)

preprint2026arXiv

Self-assembled versus biological pattern formation in geology

Both abiotic self-organization and biological mechanisms have been put forward as the origin of a number of geological patterns. It is important to comprehend the formation mechanisms of such structures both to understand geological self-organization and in order to differentiate them from biological patterns -- fossils and bio-influenced structures -- seen in geological systems. Being able to distinguish the traces of biological activity from geological self-organization is fundamental both for understanding the origin of life on Earth and for the search for life beyond Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chirality for crooked curves

Chiral objects rotate when placed in a collimated flow or wind. We exploit this hydrodynamic intuition to construct a tensorial chirality measure for rigid filaments and curves. This tensor is trace-free, so if a curve has a right-handed twist about some axis, there is a perpendicular axis about which the twist is left-handed. Our measure places minimal requirements on the smoothness of the curve, hence it can be readily used to quantify chirality for biomolecules and polymers, polygonal and rectifiable curves, and other discrete geometrical structures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stabilizing effect of tip splitting on the interface motion

Pattern-forming processes, such as electrodeposition, dielectric breakdown or viscous fingering are often driven by instabilities. Accordingly, the resulting growth patterns are usually highly branched, fractal structures. However, in some of the unstable growth processes the envelope of the structure grows in a highly regular manner, with the perturbations smoothed out over the course of time. In this paper, we show that the regularity of the envelope growth can be connected to small-scale instabilities leading to the tip splitting of the fingers at the advancing front of the structure. Whenever the growth velocity becomes too large, the finger splits into two branches. In this way it can absorb an increased flux and thus damp the instability. Hence, somewhat counterintuitively, the instability at a small scale results in a stability at a larger scale. The quantitative analysis of these effects is provided by means of the Loewner equation, which one can use to reduce the problem of the interface motion to that of the evolution of the conformal mapping onto the complex plane. This allows an effective analysis of the multifingered growth in a variety of different geometries. We show how the geometry impacts the shape of the envelope of the growing pattern and compare the results with those observed in natural systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Steadily translating parabolic dissolution fingers

Dissolution fingers (or wormholes) are formed during the dissolution of a porous rock as a result of nonlinear feedbacks between the flow, transport and chemical reactions at pore surfaces. We analyze the shapes and growth velocities of such fingers within the thin-front approximation, in which the reaction is assumed to take place instantaneously with the reactants fully consumed at the dissolution front. We concentrate on the case when the main flow is driven by the constant pressure gradient far from the finger, and the permeability contrast between the inside and the outside of the finger is finite. Using Ivantsov ansatz and conformal transformations we find the family of steadily translating fingers characterized by a parabolic shape. We derive the reactant concentration field and the pressure field inside and outside of the fingers and show that the flow within them is uniform. The advancement velocity of the finger is shown to be inversely proportional to its radius of curvature in the small Péclet number limit and constant for large Péclet numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Where do rivers grow? Path selection and growth in a harmonic field

River networks exhibit a complex ramified structure that has inspired decades of studies. Yet, an understanding of the propagation of a single stream remains elusive. Here we invoke a criterion for path selection from fracture mechanics and apply it to the growth of streams in a diffusion field. We show that a stream will follow local symmetry in order to maximize the water flux and that its trajectory is defined by the local field in its vicinity. We also study the growth of a real network. We use this principle to construct the history of a network and to find a growth law associated with it. The results show that the deterministic growth of a single channel based on its local environment can be used to characterize the structure of river networks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generalization of the Rotne-Prager-Yamakawa mobility and shear disturbance tensors

Rotne-Prager-Yamakawa approximation is one of the most commonly used methods of including hydrodynamic interactions in modelling of colloidal suspensions and polymer solutions. The two main merits of this approximation is that it includes all long-range terms (i.e. decaying as R^-3 or slower in interparticle distances) and that the diffusion matrix is positive definite, which is essential for Brownian dynamics modelling. Here, we extend the Rotne-Prager-Yamakawa approach to include both translational and rotational degrees of freedom, and derive the regularizing corrections to account for overlapping particles. Additionally, we show how the Rotne-Prager-Yamakawa approximation can be generalized for other geometries and boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interacting length scales in the reactive-infiltration instability

The reactive-infiltration instability, which develops when a porous matrix is dissolved by a flowing fluid, contains two important length scales. Here we outline a linear stability analysis that simultaneously incorporates both scales. We show that the commonly used "thin-front" model is a limiting case of a more general theory, which also includes convection-dominated dissolution as another special case. The wavelength of the instability is bounded from below, and lies in the range 1mm to 1km for physically reasonable flow rates and reaction rates. We obtain a closed form for the growth rate when the change in porosity is small.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Network models of dissolution of porous media

We investigate the chemical dissolution of porous media using a network model in which the system is represented as a series of interconnected pipes with the diameter of each segment increasing in proportion to the local reactant consumption. Moreover, the topology of the network is allowed to change dynamically during the simulation: as the diameters of the eroding pores become comparable with the interpore distances, the pores are joined together thus changing the interconnections within the network. With this model, we investigate different growth regimes in an evolving porous medium, identifying the mechanisms responsible for the emergence of specific patterns. We consider both the random and regular network and study the effect of the network geometry on the patterns. Finally, we consider practically important problem of finding an optimum flow rate that gives a maximum increase in permeability for a given amount of reactant.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Instabilities in the dissolution of a porous matrix

A reactive fluid dissolving the surrounding rock matrix can trigger an instability in the dissolution front, leading to spontaneous formation of pronounced channels or wormholes. Theoretical investigations of this instability have typically focused on a steadily propagating dissolution front that separates regions of high and low porosity. In this paper we show that this is not the only possible dissolutional instability in porous rocks; there is another instability that operates instantaneously on any initial porosity field, including an entirely uniform one. The relative importance of the two mechanisms depends on the ratio of the porosity increase to the initial porosity. We show that the "inlet" instability is likely to be important in limestone formations where the initial porosity is small and there is the possibility of a large increase in permeability. In quartz-rich sandstones, where the proportion of easily soluble material (e.g. carbonate cements) is small, the instability in the steady-state equations is dominant.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Structural transitions in hypersphere fluids: predictions of Kirkwood's approximation

We use an analytic criterion for vanishing of exponential damping of correlations developed previously (Piasecki et al, J. Chem. Phys., 133, 164507, 2010) to determine the threshold volume fractions for structural transitions in hard sphere systems in dimensions D=3,4,5 and 6, proceeding from the YBG hierarchy and using the Kirkwood superposition approximation. We conclude that the theory does predict phase transitions in qualitative agreement with numerical studies. We also derive, within the superposition approximation, the asymptotic form of the analytic condition for occurence of a structural transition in the D->Infinity limit .

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Prediction of a Structural Transition in the Hard Disk Fluid

Starting from the second equilibrium equation in the BBGKY hierarchy under the Kirkwood superposition closure, we implement a new method for studying the asymptotic decay of correlations in the hard disk fluid in the high density regime. From our analysis and complementary numerical studies, we find that exponentially damped oscillations can occur only up to a packing fraction η*~0.718, a value which is in substantial agreement with the packing fraction, η~0.723, believed to characterize the transition from the ordered solid phase to a dense fluid phase, as inferred from Mak's Monte Carlo simulations [Phys. Rev. E 73, 065104 (2006)]. We next show that the same method of analysis predicts that exponential damping of oscillations in the hard sphere fluid becomes impossible when λ= 4nπσ^3 [1 + H(1)]>/- 34.81, where H(1) is the contact value of the correlation function, n is the number density and σ is the sphere diameter, in exact agreement with the condition, λ>/- 34.8, first reported in a numerical study of the Kirkwood equation by Kirkwood et al. [J. Chem. Phys. 18, 1040 (1950)]. Finally, we show that our method confirms the absence of any structural transition in hard rods for the entire range of densities below close packing.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The initial stages of cave formation: Beyond the one-dimensional paradigm

The solutional origin of limestone caves was recognized over a century ago, but the short penetration length of an undersaturated solution made it seem impossible for long conduits to develop. This is contradicted by field observations, where extended conduits, sometimes several kilometers long, are found in karst environments. However, a sharp drop in the dissolution rate of CaCO_3 near saturation provides a mechanism for much deeper penetration of reactant. The notion of a "kinetic trigger" - a sudden change in rate constant over a narrow concentration range - has become a widely accepted paradigm in speleogenesis modeling. However, it is based on one-dimensional models for the fluid and solute transport inside the fracture, assuming that the dissolution front is planar in the direction perpendicular to the flow. Here we show that this assumption is incorrect; a planar dissolution front in an entirely uniform fracture is unstable to infinitesimal perturbations and inevitably breaks up into highly localized regions of dissolution. This provides an alternative mechanism for cave formation, even in the absence of a kinetic trigger. Our results suggest that there is an inherent wavelength to the erosion pattern in dissolving fractures, which depends on the reaction rate and flow rate, but is independent of the initial roughness. In contrast to one-dimensional models, two-dimensional simulations indicate that there is only a weak dependence of the breakthrough time on kinetic order; localization of the flow tends to keep the undersaturation in the dissolution front above the threshold for non-linear kinetics.

0 citations0 reviewsNo code linkedOpen access