Source author record

Sarthok Sircar

Sarthok Sircar 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

8works
9topics
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

8 published item(s)

preprint2023arXiv

Implicit-explicit time integration method for fractional advection-reaction-diffusion equations

We propose a novel family of asymptotically stable, implicit-explicit, adaptive, time integration method (denoted with the $θ$-method) for the solution of the fractional advection-diffusion-reaction (FADR) equations. This family of time integration method generalized the computationally explicit $L_1$-method adopted by Brunner (J. Comput. Phys. {\bf 229} 6613-6622 (2010)) as well as the fully implicit method proposed by Jannelli (Comm. Nonlin. Sci. Num. Sim., {\bf 105}, 106073 (2022)). The spectral analysis of the method (involving the group velocity and the phase speed) indicates a region of favorable dispersion for a limited range of Peclet number. The numerical inversion of the coefficient matrix is avoided by exploiting the sparse structure of the matrix in the iterative solver for the Poisson equation. The accuracy and the efficacy of the method is benchmarked using (a) the two-dimensional (2D) fractional diffusion equation, originally proposed by Brunner, and (b) the incompressible, subdiffusive dynamics of a planar viscoelastic channel flow of the Rouse chain melts (FADR equation with fractional time-derivative of order $α=\nicefrac{1}{2}$) and the Zimm chain solution ($α=\nicefrac{2}{3}$). Numerical simulations of the viscoelastic channel flow effectively capture the non-homogeneous regions of high viscosity at low fluid inertia (or the so-called `spatiotemporal macrostructures'), experimentally observed in the flow-instability transition of subdiffusive flows.

preprint2023arXiv

Spatiotemporal linear stability of viscoelastic subdiffusive channel flows: a fractional calculus framework

The temporal and spatiotemporal linear stability analyses of viscoelastic, subdiffusive, plane Poiseuille and Couette flows obeying the Fractional Upper Convected Maxwell (FUCM) equation in the limit of low to moderate Reynolds number ($Re$) and Weissenberg number ($We$), is reported to identify the regions of topological transition of the advancing flow interface. In particular, we demonstrate how the exponent in the subdiffusive power{\color{black}-}law scaling ($t^α$, with $0 {\color{black} < } α\le 1$) of the mean square displacement of the tracer particle, in the microscale [Mason and Weitz, Phys. Rev. Lett. {\bf 74}, 1250-1253 (1995)] is related to the fractional order of the derivative, $α$, of the corresponding non-linear stress constitutive equation in the continuum. The stability studies are limited to two exponents: monomer diffusion in Rouse chain melts, $α=\nicefrac{1}{2}$, and in Zimm chain solutions, $α=\nicefrac{2}{3}$. The temporal stability analysis indicates that with decreasing order of the fractional derivative: (a) the most unstable mode decreases, (b) the peak of the most unstable mode shifts to lower values of $Re$, and (c) the peak of the most unstable mode, for the Rouse model precipitates towards the limit $Re \rightarrow 0$. The Briggs idea of analytic continuation is deployed to classify regions of temporal stability, absolute and convective instabilities and evanescent modes. The spatiotemporal phase diagram indicates an abnormal region of temporal stability at high fluid inertia, revealing the presence of a non-homogeneous environment with hindered flow, thus highlighting the potential of the model to effectively capture certain experimentally observed, flow-instability transition in subdiffusive flows.

preprint2016arXiv

Dynamics of suspended rigid aggregating particles in flowing medium: theory, analysis and scientific computing

We develop and present a unified multi-scale model (involving three scales of spatial organisation) to study the dynamics of rigid aggregating particles suspended in a viscous fluid medium and subject to a steady poiseuille flow. At micro-level, the theory of adhesion describing the attachment / detachment kinetics of two rigid spheres coated with binding ligands, is utilized to describe the collision frequency function. The meso-scale dynamics is outlined through a continuous general dynamic equation governing the time rate of change of the particle size distribution function. The micro-meso coupling is achieved via the balancing of the mesoscale drag forces and couples with the micro-scale forces associated with the binder kinetics. Inside the macro domain (i.e., a long pipe), the model is equation free and divided into equal sized patches. The macroscale solution within each patch is obtained via appropriate (extrapolatory) coupling and amplitude conditions.

preprint2015arXiv

Ion mediated cross-shield driven mucous swelling kinetics

We present an experimentally guided, bi-phasic, multi-species ionic gel model to compare and make quantitative predictions on the viscoelastic properties of healthy mucus versus Cystic-Fibrosis infected mucus. The mixture theory consists of the mucus (polymer phase) and water (solvent phase) as well as several different ions: H$^+$, Na$^+$, Cl$^-$ and Ca$^{2+}$. The model is calibrated with the experimental data of mucus diffusivities. Numerical studies reveal that the Donnan forces are the dominating mechanism driving the mucus swelling/deswelling transition. However, these transitions could be continuous or discontinuous, depending on the complex interplay between ionic bath concentrations, pH, polymer mass in the solvent and the average charge per monomer.

preprint2015arXiv

Ligand mediated adhesive mechanics of two deformed spheres

A self-consistent model is developed to investigate attachment / detachment kinetics of two soft, deformable microspheres with irregular surface and coated with flexible binding ligands. The model highlights how the microscale binding kinetics of these ligands as well as the attractive/repulsive potential of the charged surface affects the static deformed configuration of the spheres. It is shown that in the limit of smooth, neutral charged surface (i.e., the Debye length, $κ\rightarrow \infty $), interacting via elastic binders (i.e., the stiffness coefficient, $λ\rightarrow 0$) the adhesion mechanics approaches the regime of application of the JKR theory, and in this particular limit, the contact radius scales with the particle radius, according to the scaling law, $R_c\propto R^{\frac{2}{3}}$. We show that adhesion dominates in larger particles with highly charged surface and with resilient binders. Normal stress distribution within the contact area fluctuates with the binder stiffness coefficient, from a maximum at the center to a maximum at the periphery of the region. Surface heterogeneities result in a diminished adhesion with a distinct reduction in the pull off force, larger separation gap, weaker normal stress and limited area of adhesion. These results are in agreement with the published experimental findings.

preprint2015arXiv

Surface deformation and shear flow in ligand mediated cell adhesion

We present a single, unified, multi-scale model to study the attachment\detachment dynamics of two deforming, near spherical cells, coated with binding ligands and subject to a slow, homogeneous shear flow in a viscous fluid medium. The binding ligands on the surface of the cells experience attractive and repulsive forces in an ionic medium and exhibit finite resistance to rotation via bond tilting. The macroscale drag forces and couples describing the fluid flow inside the small separation gap between the cells, are calculated using a combination of methods in lubrication theory and previously published numerical results. For a select range of material and fluid parameters, a hysteretic transition of the sticking probability curves between the adhesion and fragmentation domain is attributed to a nonlinear relation between the total microscale binding forces and the separation gap between the cells. We show that adhesion is favored in highly ionic fluids, increased deformability of the cells, elastic binders and a higher fluid shear rate (until a critical value). Continuation of the limit points predict a bistable region, indicating an abrupt switching between the adhesion and fragmentation regimes at critical shear rates, and suggesting that adhesion of two deformable surfaces in shearing fluids may play a significant dynamical role in some cell adhesion applications.

preprint2013arXiv

Determining equilibrium osmolarity in Poly(ethylene glycol) / Chondrotin sulfate gels mimicking articular cartilage

We present an experimentally guided, multi-phase, multi-species polyelectrolyte gel model to make quantitative predictions on the electro-chemical properties of articular cartilage. The mixture theory consists of two different types of polymers: Poly(ethylene gylcol) (PEG), Chondrotin sulfate (ChS), water (acting as solvent) and several different ions: H$^+$, Na$^+$, Cl$^-$. The polymer chains have covalent cross-links modeled using Doi rubber elasticity theory. Numerical studies on polymer volume fraction and net osmolarity (difference in the solute concentration across the gel) show the interplay between ionic bath concentrations, pH, polymer mass in the solvent and the average charge per monomer; governing the equilibrium swelled / de-swelled state of the gel. We conclude that swelling is aided due to a higher average charge per monomer (or a higher percentage of charged ChS component of the polymer), low solute concentration in the bath, a high pH or a low cross-link fraction. However, the swelling-deswelling transitions could be continuous or discontinuous depending upon the relative influence of the various competing forces.

preprint2013arXiv

Sticky Surfaces: Sphere-Sphere Adhesion Dynamics

We present a multi-scale model to study the attachment of spherical particles with a rigid core, coated with binding ligands and in equilibrium with the surrounding, quiescent fluid medium. This class of fluid-immersed adhesion is widespread in many natural and engineering settings. Our theory highlights how the micro-scale binding kinetics of these ligands, as well as the attractive / repulsive surface potential in an ionic medium effects the eventual macro-scale size distribution of the particle aggregates (flocs). The results suggest that the presence of elastic ligands on the particle surface allow large floc aggregates by inducing efficient inter-floc collisions (i.e., a large, non-zero collision factor). Strong electrolytic composition of the surrounding fluid favors large floc formation as well.