Source author record

Lou Kondic

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

21 published item(s)

preprint2026arXiv

Topological Data Analysis combined with Machine Learning for Predicting Permeability of Porous Media

Flow in porous media is difficult to address using standard analytical or numerical methods due to its complexity. However, since synthetic representations of porous media are easy to produce and data from physical experiments are becoming more widely available, the problem is well-suited to studies that include machine learning (ML) techniques. We discuss a number of features that can be extracted from such data, and their utility as input variables into a standard ML algorithm. These features include structural measures describing the geometry of the porous media, topological measures describing the connectivity, and network measures obtained by modeling the porous media as simplified pore networks. These features enable the prediction of the permeability of the considered (synthetic) porous materials using ML techniques that also leverage the separately computed exact permeability (ground truth). Comparing results obtained using different input variables helps develop a better understanding of the utility of various measures for predicting permeability based on the porous media structure. We show, in particular, that topological data analysis (TDA) provides a useful set of features that can be easily combined with ML to yield meaningful results.

preprint2022arXiv

A Graphical Representation of Membrane Filtration

We analyze the performance of membrane filters represented by pore networks using two criteria: 1) total volumetric throughput of filtrate over the filter lifetime and 2) accumulated foulant concentration in the filtrate. We first formulate the governing equations of fluid flow on a general network, and we model transport and adsorption of particles (foulants) within the network by imposing an advection equation with a sink term on each pore (edge) as well as conservation of fluid and foulant volumetric flow rates at each pore junction (network vertex). Such a setup yields a system of partial differential equations on the network. We study the influence of three geometric network parameters on filter performance: 1) average number of neighbors of each vertex; 2) initial total void volume of the pore network; and 3) tortuosity of the network. We find that total volumetric throughput depends more strongly on the initial void volume than on average number of neighbors. Tortuosity, however, turns out to be a universal parameter, leading to almost perfect collapse of all results for a variety of different network architectures. In particular, the accumulated foulant concentration in the filtrate shows an exponential decay as tortuosity increases.

preprint2022arXiv

Network-based membrane filters: Influence of network and pore size variability on filtration performance

We model porous membrane filters as networks of connected cylindrical pores via a random network generation protocol, and their initial pore radii via a uniform distribution of widths that vary about some mean value. We investigate the influence of network and pore size (radius) variations on the performance of membrane filters that undergo adsorptive fouling. We find that membrane porosity variations, independently of whether induced by variations of the pore radii or of the random pore network, are an important factor determining membrane filter performance. Network and pore size variations still play a role, in particular if pore radii variations are significant. To quantify the influence of these variations, we compare the performance metrics of networks built from pores of variable radii to their (equal porosity) counterparts built from pores of uniform radius. We show that the effect of pore radii variations is to increase throughput, but also to reduce foulant control.

preprint2022arXiv

On intermittency in sheared granular systems

We consider a system of granular particles, modeled by two dimensional frictional elastic disks, that is exposed to externally applied time-dependent shear stress in a planar Couette geometry. We concentrate on the external forcing that produces intermittent dynamics of stick-slip type. In this regime, the top wall remains almost at rest until the applied stress becomes sufficiently large, and then it slips. We focus on the evolution of the system as it approaches a slip event. Our main finding is that there are two distinct groups of measures describing system behavior before a slip event. The first group consists of global measures defined as system-wide averages at a fixed time. Typical examples of measures in this group are averages of the normal or tangent forces acting between the particles, system size and number of contacts between the particles. These measures do not seem to be sensitive to an approaching slip event. On average, they tend to increase linearly with the force pulling the spring. The second group consists of the time-dependent measures that quantify the evolution of the system on a micro (particle) or mesoscale. Measures in this group first quantify the temporal differences between two states and only then aggregate them to a single number. For example, Wasserstein distance quantitatively measures the changes of the force network as it evolves in time while the number of broken contacts quantifies the evolution of the contact network. The behavior of the measures in the second group changes dramatically before a slip event starts. They increase rapidly as a slip event approaches, indicating a significant increase in fluctuations of the system before a slip event is triggered.

preprint2022arXiv

Thin viscoelastic dewetting films of Jeffreys type subjected to gravity and substrate interactions

This work presents a study of the interfacial dynamics of thin viscoelastic films subjected to the gravitational force and substrate interactions induced by the disjoining pressure, in two spatial dimensions. The governing equation is derived as a long-wave approximation of the Navier-Stokes equations for incompressible viscoelastic liquids under the effect of gravity, with the Jeffreys model for viscoelastic stresses. For the particular cases of horizontal or inverted planes, the linear stability analysis is performed to investigate the influence of the physical parameters involved on the growth rate and length scales of instabilities. Numerical simulations of the nonlinear regime of the dewetting process are presented for the particular case of an inverted plane. Both gravity and the disjoining pressure are found to affect not only the length scale of instabilities, but also the final configuration of dewetting, by favoring the formation of satellite droplets, that are suppressed by the slippage with the solid substrate.

preprint2020arXiv

Failure of confined granular media due to pullout of an intruder: From force networks to a system wide response

We investigate computationally the pullout of a spherical intruder initially buried at the bottom of a granular column. The intruder starts to move out of the granular bed once the pulling force reaches a critical value, leading to material failure. The failure point is found to depend on the diameter of the granular column, pointing out the importance of particle-wall interaction in determining the material response. Discrete element simulations show that prior to failure, the contact network is essentially static, with only minor rearrangements of the particles. However, the force network, which includes not only the contact information, but also the information about the interaction strength, undergoes a nontrivial evolution. An initial insight is reached by considering the relative magnitudes of normal and tangential forces between the particles, and in particular the proportion of contacts that reach Coulomb threshold. More detailed understanding of the processes leading to failure is reached by the analysis of both spatial and temporal properties of the force network using the tools of persistent homology. We find that the forces between the particles undergo intermittent temporal variations ahead of the failure. In addition to this temporal intermittency, the response of the force network is found to be spatially dependent and influenced by proximity to the intruder. Furthermore, the response is modified significantly by the interaction strength, with the relevant measures describing the response showing differing behavior for the contacts characterized by large interaction forces.

preprint2019arXiv

Interaction network analysis in shear thickening suspensions

Dense, stabilized, frictional particulate suspensions in a viscous liquid undergo increasingly strong continuous shear thickening (CST) as the solid packing fraction, $ϕ$, increases above a critical volume fraction, and discontinuous shear thickening (DST) is observed for even higher packing fractions. Recent studies have related shear thickening to a transition from mostly lubricated to predominantly frictional contacts with the increase in stress. The rheology and networks of frictional forces from two and three-dimensional simulations of shear-thickening suspensions are studied. These are analyzed using measures of the topology of the network, including tools of persistent homology. We observe that at low stress the frictional interaction networks are predominantly quasi-linear along the compression axis. With an increase in stress, the force networks become more isotropic, forming loops in addition to chain-like structures. The topological measures of Betti numbers and total persistence provide a compact means of describing the mean properties of the frictional force networks and provide a key link between macroscopic rheology and the microscopic interactions. A total persistence measure describing the significance of loops in the force network structure, as a function of stress and packing fraction, shows behavior similar to that of relative viscosity and displays a scaling law near the jamming fraction for both dimensionalities simulated.

preprint2019arXiv

Intruder in a two-dimensional granular system: Effects of dynamic and static basal friction on stick-slip and clogging dynamics

We discuss the results of simulations of an intruder pulled through a two-dimensional granular system by a spring, using a model designed to lend insight into the experimental findings described by Kozlowski et al. [Phys. Rev. E, 100, 032905 (2019)]. In that previous study the presence of basal friction between the grains and the base was observed to change the intruder dynamics from clogging to stick-slip. Here we first show that our simulation results are in excellent agreement with the experimental data for a variety of experimentally accessible friction coefficients governing interactions of particles with each other and with boundaries. Then, we use simulations to explore a broader range of parameter space, focusing on the friction between the particles and the base. We consider a range of both static and dynamic basal friction coefficients, which are difficult to vary smoothly in experiments. The simulations show that dynamic friction strongly affects the stick-slip behaviour when the coefficient is decreased below 0.1, while static friction plays only a marginal role in the intruder dynamics.

preprint2019arXiv

Modeling and design optimization for pleated membrane filters

Pleated membrane filters, which offer larger surface area to volume ratios than unpleated membrane filters, are used in a wide variety of applications. However, the performance of the pleated filter, as characterized by a flux-throughput plot, indicates that the equivalent unpleated filter provides better performance under the same pressure drop. Earlier work (Sanaei & Cummings 2016) used a highly-simplified membrane model to investigate how the pleating effect and membrane geometry affect this performance differential. In this work, we extend this line of investigation and use asymptotic methods to couple an outer problem for the flow within the pleated structure to an inner problem that accounts for the pore structure within the membrane. We use our new model to formulate and address questions of optimal membrane design for a given filtration application.

preprint2016arXiv

Granular Impact: A Grain-scale Approach

This work summarizes a series of studies on two-dimensional granular impact, where an intruding object strikes a granular material at high speed. Many previous studies on granular impact have used a macroscopic force law, which is dominated by an inertial drag term proportional to the intruder velocity squared. The primary focus here is on the microscopic force response of the granular material, and how the grain-scale effects give rise to this inertial drag term. We show that the inertial drag arises from intermittent collisions with force-chain-like structures. We construct a simple collisional model to explain the inertial drag, as well as off-axis instability and rotations. Finally, we show how the granular response changes when the intruder speed approaches $d/t_c$, leading to a failure of the inertial drag description in this regime. Here, $d$ is the mean particle diameter and $t_c$ the characteristic momentum-transfer time between two grains.

preprint2016arXiv

Instability of nanometric fluid films on a thermally conductive substrate

We consider thin fluid films placed on thermally conductive substrates and exposed to time-dependent spatially uniform heat source. The evolution of the films is considered within the long-wave framework in the regime such that both fluid/substrate interaction, modeled via disjoining pressure, and Marangoni forces, are relevant. We analyze the problem by the means of linear stability analysis as well as by time-dependent nonlinear simulations. The main finding is that when self-consistent computation of the temperature field is performed, a complex interplay of different instability mechanisms results. This includes either monotonous or oscillatory dynamics of the free surface. In particular, we find that the oscillatory behavior is absent if the film temperature is assumed to be slaved to the current value of the film thickness. The results are discussed within the context of liquid metal films, but are of relevance to dynamics of any thin film involving variable temperature of the free surface, such that the temperature and the film interface itself evolve on comparable time scales.

preprint2016arXiv

Interfacial Dynamics of Thin Viscoelastic Films and Drops

We present a computational investigation of thin viscoelastic films and drops on a solid substrate subject to the van der Waals interaction force, in two spatial dimensions. The governing equations are obtained within a long-wave approximation of the Navier-Stokes equations with Jeffreys model for viscoelastic stresses. We investigate the effects of viscoelasticity, Newtonian viscosity, and the substrate slippage on the dynamics of thin viscoelastic films. We also study the effects of viscoelasticity on drops that spread or recede on a prewetted substrate. For dewetting films, the numerical results show the presence of multiple secondary droplets for higher values of elasticity, consistently with experimental findings. For drops, we find that elastic effects lead to deviations from the Cox-Voinov law for partially wetting fluids. In general, elastic effects enhance spreading, and suppress retraction, compared to Newtonian ones.

preprint2016arXiv

Steady Flow Dynamics during Granular Impact

We study experimentally and computationally the dynamics of granular flow during impacts, where intruders strike a collection of disks from above. In the regime where granular force dynamics are much more rapid than the intruder motion, we find that the particle flow near the intruder is proportional to the instantaneous intruder speed; it is essentially constant when normalized by that speed. The granular flow is nearly divergence-free and remains in balance with the intruder, despite the latter's rapid deceleration. Simulations indicate that this observation is insensitive to grain properties, which can be explained by the separation of time scales between intergrain force dynamics and intruder dynamics. Assuming there is a comparable separation of time scales, we expect our results are applicable to a broad class of dynamic or transient granular flows. Our results suggest that descriptions of static-in-time granular flows might be extended or modified to describe these dynamic flows. Additionally, we find that accurate grain-grain interactions are not necessary to correctly capture the granular flow in this regime.

preprint2015arXiv

A volume of fluid method for simulating fluid/fluid interfaces in contact with solid boundaries

In this paper, we present a novel approach to model the fluid/solid interaction forces in a direct solver of the Navier-Stokes equations based on the volume of fluid interface tracking method. The key ingredient of the model is the explicit inclusion of the fluid/solid interaction forces into the governing equations. We show that the interaction forces lead to a partial wetting condition and in particular to a natural definition of the equilibrium contact angle. We present two numerical methods to discretize the interaction forces that enter the model; these two approaches differ in complexity and convergence. To validate the computational framework, we consider the application of these models to simulate two-dimensional drops at equilibrium, as well as drop spreading. We demonstrate that the model, by including the underlying physics, captures contact line dynamics for arbitrary contact angles. More generally, the approach permits novel means to study contact lines, as well as a diverse range of phenomena that previously could not be addressed in direct simulations.

preprint2015arXiv

Nonlinear Force Propagation during Granular Impact

We experimentally study nonlinear force propagation into granular material during impact from an intruder, and we explain our observations in terms of the nonlinear grain-scale force relation. Using high-speed video and photoelastic particles, we determine the speed and spatial structure of the force response just after impact. We show that these quantities depend on a dimensionless parameter, $M'=t_c v_0/d$, where $v_0$ is the intruder speed at impact, $d$ is the particle diameter, and $t_c$ is the collision time for a pair of grains impacting at relative speed $v_0$. The experiments access a large range of $M'$ by using particles of three different materials. When $M' \ll 1$, force propagation is chain-like with a speed, $v_f$, satisfying $v_f \propto d/t_c$. For larger $M'$, the force response becomes spatially dense and the force propagation speed departs from $v_f\propto d/t_c$, corresponding to collective stiffening of a strongly compressed packing of grains.

preprint2015arXiv

Scaling Properties of Force Networks for Compressed Particulate Systems

We consider, computationally and experimentally, the scaling properties of force networks in the systems of circular particles exposed to compression in two spatial dimensions. The simulations consider polydisperse and monodisperse particles, both frictional and frictionless, and in experiments we use monodisperse and bidisperse particles. While for some of the considered systems we observe consistent scaling exponents describing the behavior of the force networks, we find that this behavior is it not universal. In particular, monodisperse frictionless systems that partially crystallize under compression, show scaling properties that are significantly different compared to the other considered systems. The findings of non-universality are confirmed by explicitly computing fractal dimension for the considered systems. The results of physical experiments are consistent with the results obtained in simulations of frictional particles.

preprint2014arXiv

Evolution of Force Networks in Dense Particulate Media

We introduce novel sets of measures with the goal of describing dynamical properties of force networks in dense particulate systems. The presented approach is based on persistent homology and allows for extracting precise, quantitative measures that describe the evolution of geometric features of the interparticle forces, without necessarily considering the details related to individual contacts between particles. The networks considered emerge from discrete element simulations of two dimensional particulate systems consisting of compressible frictional circular disks. We quantify the evolution of the networks for slowly compressed systems undergoing jamming transition. The main findings include uncovering significant but localized changes of force networks for unjammed systems, global (system-wide) changes as the systems evolve through jamming, to be followed by significantly less dramatic evolution for the jammed states. We consider both connected components, related in loose sense to force chains, and loops, and find that both measures provide a significant insight into the evolution of force networks. In addition to normal, we consider also tangential forces between the particles and find that they evolve in the consistent manner. Consideration of both frictional and frictionless systems leads us to the conclusion that friction plays a significant role in determining the dynamical properties of the considered networks. We find that the proposed approach describes the considered networks in a precise yet tractable manner, allowing to identify novel features which could be difficult or impossible to describe using other approaches.

preprint2013arXiv

Modelling spreading dynamics of liquid crystals in three spatial dimensions

We study spreading dynamics of nematic liquid crystal droplets within the framework of the long-wave approximation. A fourth order nonlinear parabolic partial differential equation governing the free surface evolution is derived. The influence of elastic distortion energy and of imposed anchoring variations at the substrate are explored through linear stability analysis and scaling arguments, which yield useful insight and predictions for the behaviour of spreading droplets. This behaviour is captured by fully nonlinear time-dependent simulations of three dimensional droplets spreading in the presence of anchoring variations that model simple defects in the nematic orientation at the substrate.

preprint2013arXiv

Note on the hydrodynamic description of thin nematic films: strong anchoring model

We discuss the long-wave hydrodynamic model for a thin film of nematic liquid crystal in the limit of strong anchoring at the free surface and at the substrate. We rigorously clarify how the elastic energy enters the evolution equation for the film thickness in order to provide a solid basis for further investigation: several conflicting models exist in the literature that predict qualitatively different behaviour. We consolidate the various approaches and show that the long-wave model derived through an asymptotic expansion of the full nemato-hydrodynamic equations with consistent boundary conditions agrees with the model one obtains by employing a thermodynamically motivated gradient dynamics formulation based on an underlying free energy functional. As a result, we find that in the case of strong anchoring the elastic distortion energy is always stabilising. To support the discussion in the main part of the paper, an appendix gives the full derivation of the evolution equation for the film thickness via asymptotic expansion.

preprint2012arXiv

Particle Scale Dynamics in Granular Impact

We perform an experimental study of granular impact, where intruders strike 2D beds of photoelastic disks from above. High-speed video captures the intruder dynamics and the local granular force response, allowing investigation of grain-scale mechanisms in this process. We observe rich acoustic behavior at the leading edge of the intruder, strongly fluctuating in space and time, and we show that this acoustic activity controls the intruder deceleration, including large force fluctuations at short time scales. The average intruder dynamics match previous studies using empirical force laws, suggesting a new microscopic picture, where acoustic energy is carried away and dissipated.

preprint2010arXiv

Thin films flowing down inverted substrates: two dimensional flow

We consider free surface instabilities of films flowing on inverted substrates within the framework of lubrication approximation. We allow for the presence of fronts and related contact lines, and explore the role which they play in instability development. It is found that a contact line, modeled by a commonly used precursor film model, leads to free surface instabilities without any additional natural or excited perturbations. A single parameter D=(3Ca)^{1/3}\cotα$, where $Ca$ is the capillary number and $α$ is the inclination angle, is identified as a governing parameter in the problem. This parameter may be interpreted to reflect the combined effect of inclination angle, film thickness, Reynolds number and the fluid flux. Variation of $D$ leads to change of the wave-like properties of the instabilities, allowing us to observe traveling wave behavior, mixed waves, and the waves resembling solitary ones.