Source author record

Indranil Saha Dalal

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

ResearcherUnclaimed source record

cond-mat.soft physics.comp-ph physics.flu-dyn

Catalog footprint

What is connected

3works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Formulation of an equivalent GNF model as an efficient approximation for flow of polymer solutions described by FENE-P

The molecular constitutive models, like FENE-P for polymer solutions, are known to have convergence issues at relatively larger flow rates. In this study, we investigate the possibility of a numerically efficient GNF-based approximation of FENE-P, which would closely approximate the flow field. Firstly, we compare the flow fields predicted by FENE-P and an equivalent GNF model. For these studies, we considered the flow around a sphere and selected the Carreau-Yasuda model as the representative GNF. This is made equivalent to the FENE-P by selecting parameters to equalize the viscosity-shear rate dependence. Our results show severe deficiencies of the GNF model, owing to its inability to account for chain stretching, particularly near the stagnation points. Next, the effect of extensional components on the local viscosity was added by formulating an equivalent GNF-X [Journal of Rheology 64, 493 (2020)] model. Even this failed to capture the asymmetry in the stress and flow profiles and predicted very large stresses at both stagnation points, relative to FENE-P. Hence, we proposed a novel modified formalism (denoted as GNF-XM) that was able to capture all trends successfully. The drag coefficients from GNF-XM agreed well with FENE-P predictions for all flow rates considered. Significantly, the computational times required to solve the flow field with GNF-XM is about an order of magnitude lower than that of FENE-P, especially at higher flow rates. Thus, we have successfully formulated a highly efficient GNF-based approximation to the FENE-P, whose formalism can be extended to other similarly complicated constitutive models.

preprint2022arXiv

Petascale Brownian dynamics simulations of highly resolved polymer chains with hydrodynamic interactions using modern GPUs

Brownian dynamics simulations of fairly long, highly detailed polymer chains, at the resolution of a single Kuhn step, remains computationally prohibitive even on the modern processors. This is especially true when the beads on the chain experience hydrodynamic interactions (HI), which requires the usage of methods like Cholesky decomposition of large matrices at every timestep. In this study, we perform Petascale BD simulations, with HI, of fairly long and highly resolved polymer chains on modern GPUs. Our results clearly highlight the inadequacies of the typical models that use beads connected by springs. In this manuscript, firstly, we present the details of a highly scalable, parallel hybrid code implemented on a GPU for BD simulations of chains resolved to a single Kuhn step. In this hybrid code using CUDA and MPI, we have incorporated HI using the Cholesky decomposition method. Next, we validate the GPU implementation extensively with theoretical expectations for polymer chains at equilibrium and in flow with results in the absence of HI. Further, our results in flow with HI show significantly different temporal variations of stretch, in both startup extensional and shear flows, relative to the conventional bead-spring models. In all cases investigated, the ensemble averaged chain stretch is much lower than bead-spring predictions. Also, quite remarkably, our GPU implementation shows a scaling of $\sim$$N^{1.2}$ and $\sim$$N^{2.2}$ of the computational times for shorter and longer chains in the most modern available GPU, respectively, which is significantly lower than the theoretically expected $\sim$$N^{3}$. We expect our methods and results to pave the way for further analysis of polymer physics in flow fields, with long and highly detailed chain models.

preprint2020arXiv

Seeding Method for Ice Nucleation under Shear

Hydrodynamic flow can have complex and far-reaching consequences on the rate of homogenous nucleation. We present a general formalism for calculating the nucleation rates of simply sheared systems. We have derived an extension to the conventional Classical Nucleation Theory, explicitly embodying the shear rate. Seeded Molecular Dynamics simulations form the backbone of our approach. The framework can be used for moderate supercoolings, at which temperatures brute-force methods are practically infeasible. The competing energetic and kinetic effects of shear arise naturally from the equations. We show how the theory can be used to identify shear regimes of ice nucleation behaviour for the mW water model, unifying disparate trends reported in the literature. At each temperature, we define a crossover shear rate in the limit of $1000-10,000 \ s^{-1}$, beyond which the nucleation rate increases steadily upto a maximum, at the optimal shear rate. For $235$, $240$, $255$ and $260 \ K$, the optimal shear rates are in the range of $\approx 10^6-10^7 \ s^{-1}$. For very high shear rates beyond $10^8 \ s^{-1}$, nucleation is strongly inhibited. Our results indicate that the shear-dependent nucleation rate curves have a non-monotonic dependence on temperature.