BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ivan C. Christov

Ivan C. Christov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dyncond-mat.softphysics.comp-phComputational Engineering, Finance, and Sciencemath.APmath.NAnlin.PS

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

10 published item(s)

preprint2025arXiv

Oscillatory flows in three-dimensional deformable microchannels

Deformable microchannels emulate a key characteristic of soft biological systems and flexible engineering devices: the flow-induced deformation of the conduit due to slow viscous flow within. Elucidating the two-way coupling between oscillatory flow and deformation of a three-dimensional (3D) rectangular channel is crucial for designing lab- and organ-on-a-chip microsystems and eventually understanding flow-structure instabilities that can enhance mixing and transport. To this end, we determine the axial variations of the primary flow, pressure, and deformation for Newtonian fluids in the canonical geometry of a slender (long) and shallow (wide) 3D rectangular channel with a deformable top wall under the assumption of weak compliance and without restriction on the oscillation frequency (\textit{i.e.}, on the Womersley number). Unlike rigid conduits, the pressure distribution is not linear with the axial coordinate. To validate this prediction, we design a PDMS-based experimental platform with a speaker-based flow-generation apparatus and a pressure acquisition system with multiple ports along the axial length of the channel. The experimental measurements show good agreement with the predicted pressure profiles across a wide range of the key dimensionless quantities: the Womersley number, the compliance number, and the elastoviscous number. Finally, we explore how the nonlinear flow-deformation coupling leads to self-induced streaming (rectification of the oscillatory flow). Following Zhang and Rallabandi (\textit{J.\ Fluid Mech.}, vol.~996, 2024, A16), we develop a theory for the cycle-averaged pressure based on the primary problem's solution, and we validate the predictions for the axial distribution of the streaming pressure against the experimental measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hydrodynamic Bulge Testing: Materials Characterization without Measuring Deformation

Characterizing the elastic properties of soft materials through bulge testing relies on accurate measurement of deformation, which is experimentally challenging. To avoid measuring deformation, we propose a hydrodynamic bulge test for characterizing the material properties of thick, pre-stressed elastic sheets via their fluid--structure interaction with a steady viscous fluid flow. Specifically, the hydrodynamic bulge test relies on a pressure drop measurement across a rectangular microchannel with a deformable top wall. We develop a mathematical model using first-order shear-deformation theory of plates with stretching, and the lubrication approximation for Newtonian fluid flow. Specifically, a relationship is derived between the imposed flow rate and the total pressure drop. Then, this relationship is inverted numerically to yield estimates of the Young's modulus (given the Poisson ratio), if the pressure drop is measured (given the steady flow rate). Direct numerical simulations of two-way-coupled fluid--structure interaction are carried out in ANSYS to determine the cross-sectional membrane deformation and the hydrodynamic pressure distribution. Taking the simulations as ``ground truth,'' a hydrodynamic bulge test is performed using the simulation data to ascertain the accuracy and validity of the proposed methodology for estimating material properties. An error propagation analysis is performed via Monte Carlo simulation to characterize the susceptibility of the hydrodynamic bulge test estimates to noise. We find that, while a hydrodynamic bulge test is less accurate in characterizing material properties, it is less susceptible to noise, in the input (measured) variable, than a hydrostatic bulge test.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Long-time asymptotics of non-degenerate non-linear diffusion equations

We study the long-time asymptotics of prototypical non-linear diffusion equations. Specifically, we consider the case of a non-degenerate diffusivity function that is a (non-negative) polynomial of the dependent variable of the problem. We motivate these types of equations using Einstein's random walk paradigm, leading to a partial differential equation in non-divergence form. On the other hand, using conservation principles leads to a partial differential equation in divergence form. A transformation is derived to handle both cases. Then, a maximum principle (on both an unbounded and a bounded domain) is proved, in order to obtain bounds above and below for the time-evolution of the solutions to the non-linear diffusion problem. Specifically, these bounds are based on the fundamental solution of the linear problem (the so-called Aronson's Green function). Having thus sandwiched the long-time asymptotics of solutions to the non-linear problems between two fundamental solutions of the linear problem, we prove that, unlike the case of degenerate diffusion, a non-degenerate diffusion equation's solution converges onto the linear diffusion solution at long times. Select numerical examples support the mathematical theorems and illustrate the convergence process. Our results have implications on how to interpret asymptotic scalings of potentially anomalous diffusion processes (such as in the flow of particulate materials) that have been discussed in the applied physics literature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revisiting steady viscous flow of a generalized Newtonian fluid through a slender elastic tube using shell theory

A flow vessel with an elastic wall can deform significantly due to viscous fluid flow within it, even at vanishing Reynolds number (no fluid inertia). Deformation leads to an enhancement of throughput due to the change in cross-sectional area. The latter gives rise to a non-constant pressure gradient in the flow-wise direction and, hence, to a nonlinear flow rate--pressure drop relation (unlike the Hagen--Poiseuille law for a rigid tube). Many biofluids are non-Newtonian, and are well approximated by generalized Newtonian (say, power-law) rheological models. Consequently, we analyze the problem of steady low Reynolds number flow of a generalized Newtonian fluid through a slender elastic tube by coupling fluid lubrication theory to a structural problem posed in terms of Donnell shell theory. A perturbative approach (in the slenderness parameter) yields analytical solutions for both the flow and the deformation. Using matched asymptotics, we obtain a uniformly valid solution for the tube's radial displacement, which features both a boundary layer and a corner layer caused by localized bending near the clamped ends. In doing so, we obtain a ``generalized Hagen--Poiseuille law'' for soft microtubes. We benchmark the mathematical predictions against three-dimensional two-way coupled direct numerical simulations (DNS) of flow and deformation performed using the commercial computational engineering platform by ANSYS. The simulations show good agreement and establish the range of validity of the theory. Finally, we discuss the implications of the theory on the problem of the flow-induced deformation of a blood vessel, which is featured in some textbooks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Static response of deformable microchannels: A comparative modelling study

We present a comparative modelling study of fluid-structure interactions in microchannels. Through a mathematical analysis based on plate theory and the lubrication approximation for low-Reynolds-number flow, we derive models for the flow rate-pressure drop relation for long shallow microchannels with both thin and thick deformable top walls. These relations are tested against full three-dimensional two-way-coupled fluid-structure interaction simulations. Three types of microchannels, representing different elasticity regimes and having been experimentally characterized previously, are chosen as benchmarks for our theory and simulations. Good agreement is found in most cases for the predicted, simulated and measured flow rate-pressure drop relationships. The numerical simulations performed allow us to also carefully examine the deformation profile of the top wall of the microchannel in any cross section, showing good agreement with the theory. Specifically, the prediction that span-wise displacement in a long shallow microchannel decouples from the flow-wise deformation is confirmed, and the predicted scaling of the maximum displacement with the hydrodynamic pressure and the various material and geometric parameters is validated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time-averaged transport in oscillatory squeeze flow of a viscoelastic fluid

Periodically-driven flows are known to generate non-zero, time-averaged fluxes of heat or solute species, due to the interactions of out-of-phase velocity and temperature/concentration fields, respectively. Herein, we investigate such transport (a form of the well-known Taylor--Aris dispersion) in the gap between two parallel plates, one of which oscillates vertically, generating a time-periodic squeeze flow of either a newtonian or Maxwellian fluid. Using the method of multiple time-scale homogenization, the mass/heat balance equation describing transport in this flow is reduced to a one-dimensional advection--diffusion--reaction equation. This result indicates three effective mechanisms in the mass/heat transfer in the system: an effective diffusion that spreads mass/heat along the concentration/temperature gradient, an effective advective flux, and an effective reaction that releases or absorbs mass/heat - in the time-averaged frame. Our results demonstrate that there exist resonant modes under which the velocity peaks when the dimensionless plate oscillation frequency (embodied by the Womersley number, the ratio of the transient inertia to viscous forces) approaches specific values. As a result, transport in this flow is significantly influenced by the dimensionless frequency. On the one hand, the effective, time-averaged dispersion coefficient is always larger than the molecular diffusivity, and is sharply enhanced near resonance. The interaction between fluid elasticity and the oscillatory forcing enhances the efficiency of transport in the system. On the other hand, the identified effective advection and reaction mechanisms may transport mass/heat from regions of high concentration/temperature to those of low concentration/temperature, or vice versa, depending on the value of dimensionless frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tuning a magnetic field to generate spinning ferrofluid droplets with controllable speed via nonlinear periodic interfacial waves

Two dimensional free surface flows in Hele-Shaw configurations are a fertile ground for exploring nonlinear physics. Since Saffman and Taylor's work on linear instability of fluid--fluid interfaces, significant effort has been expended to determining the physics and forcing that set the linear growth rate. However, linear stability does not always imply nonlinear stability. We demonstrate how the combination of a radial and an azimuthal external magnetic field can manipulate the interfacial shape of a linearly unstable ferrofluid droplet in a Hele-Shaw configuration. We show that weakly nonlinear theory can be used to tune the initial unstable growth. Then, nonlinearity arrests the instability, and leads to a permanent deformed droplet shape. Specifically, we show that the deformed droplet can be set into motion with a predictable rotation speed, demonstrating nonlinear traveling waves on the fluid-fluid interface. The most linearly unstable wavenumber and the combined strength of the applied external magnetic fields determine the traveling wave shape, which can be asymmetric.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unsteady fluid--structure interactions in a soft-walled microchannel: A one-dimensional lubrication model for finite Reynolds number

We develop a one-dimensional model for the unsteady fluid--structure interaction (FSI) between a soft-walled microchannel and viscous fluid flow within it. A beam equation, which accounts for both transverse bending rigidity and nonlinear axial tension, is coupled to a one-dimensional fluid model obtained from depth-averaging the two-dimensional incompressible Navier--Stokes equations across the channel height. Specifically, the Navier--Stokes equations are scaled in the viscous lubrication limit relevant to microfluidics. The resulting set of coupled nonlinear partial differential equations is solved numerically through a segregated approach employing fully-implicit time stepping. We explore both the static and dynamic FSI behavior of this example microchannel system by varying a reduced Reynolds number $Re$, which necessarily changes the Strouhal number $St$, while we keep the geometry and a modified dimensionless Young's modulus $Σ$ fixed. At steady state, an order-of-magnitude analysis (balancing argument) shows that the axially-averaged pressure in the flow, $\langle P\rangle$, exhibits two different scaling regimes, while the maximum deformation of the top wall of the channel, $H_{\mathrm{max}}$, can fall into four different regimes, depending on the magnitudes of $Re$ and $Σ$. These regimes are physically explained as resulting from the competition between the inertial and viscous forces in the fluid flow as well as the bending resistance and tension in the elastic wall. Finally, the linear stability of the steady inflated microchannel shape is assessed via a modal analysis, showing the existence of many highly oscillatory but stable modes, which further highlights the computational challenge of simulating unsteady FSIs.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Computational Analysis of Interfacial Dynamics in Angled Hele-Shaw Cells: Instability Regimes

We present a theoretical and numerical study on the (in)stability of the interface between two immiscible liquids, i.e., viscous fingering, in angled Hele-Shaw cells across a range of capillary numbers ($Ca$). We consider two types of angled Hele-Shaw cells: diverging cells with a positive depth gradient and converging cells with a negative depth gradient, and compare those against parallel cells without a depth gradient. A modified linear stability analysis is employed to derive an expression for the growth rate of perturbations on the interface and for the critical capillary number ($Ca_c$) for such tapered Hele-Shaw cells with small gap gradients. Based on this new expression for $Ca_c$, a three-regime theory is formulated to describe the interface (in)stability: (i) in Regime I, the growth rate is always negative, thus the interface is stable; (ii) in Regime II, the growth rate remains zero (parallel cells), changes from negative to positive (converging cells), or from positive to negative (diverging cells), thus the interface (in)stability possibly changes type at some location in the cell; (iii) in Regime III, the growth rate is always positive, thus the interface is unstable. We conduct three-dimensional direct numerical simulations of the full Navier--Stokes equations, using a phase field method to enforce surface tension at the interface, to verify the theory and explore the effect of depth gradient on the interface (in)stability. We demonstrate that the depth gradient has only a slight influence in Regime I, and its effect is most pronounced in Regime III. Finally, we provide a critical discussion of the stability diagram derived from theoretical considerations versus the one obtained from direct numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Solving Nonlinear Parabolic Equations by a Strongly Implicit Finite-Difference Scheme

We discuss the numerical solution of nonlinear parabolic partial differential equations, exhibiting finite speed of propagation, via a strongly implicit finite-difference scheme with formal truncation error $\mathcal{O}\left[(Δx)^2 + (Δt)^2 \right]$. Our application of interest is the spreading of viscous gravity currents in the study of which these type of differential equations arise. Viscous gravity currents are low Reynolds number (viscous forces dominate inertial forces) flow phenomena in which a dense, viscous fluid displaces a lighter (usually immiscible) fluid. The fluids may be confined by the sidewalls of a channel or propagate in an unconfined two-dimensional (or axisymmetric three-dimensional) geometry. Under the lubrication approximation, the mathematical description of the spreading of these fluids reduces to solving the so-called thin-film equation for the current's shape $h(x,t)$. To solve such nonlinear parabolic equations we propose a finite-difference scheme based on the Crank--Nicolson idea. We implement the scheme for problems involving a single spatial coordinate (i.e., two-dimensional, axisymmetric or spherically-symmetric three-dimensional currents) on an equispaced but staggered grid. We benchmark the scheme against analytical solutions and highlight its strong numerical stability by specifically considering the spreading of non-Newtonian power-law fluids in a variable-width confined channel-like geometry (a "Hele-Shaw cell") subject to a given mass conservation/balance constraint. We show that this constraint can be implemented by re-expressing it as nonlinear flux boundary conditions on the domain's endpoints. Then, we show numerically that the scheme achieves its full second-order accuracy in space and time. We also highlight through numerical simulations how the proposed scheme accurately respects the mass conservation/balance constraint.

0 citations0 reviewsNo code linkedOpen access