Source author record

Ken Kamrin

Ken Kamrin 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

17works
6topics
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

17 published item(s)

preprint2021arXiv

An incompressible Eulerian method for fluid-structure interaction with mixed soft and rigid solids

We present a general simulation approach for incompressible fluid--structure interactions in a fully Eulerian framework using the reference map technique (RMT). The approach is suitable for modeling one or more rigid or finitely-deformable objects or soft objects with rigid components interacting with the fluid and with each other. It is also extended to control the kinematics of structures in fluids. The model is based on our previous Eulerian fluid--soft solver, and generalized to rigid structures by constraining the deformation-rate tensor in a projection framework. Several numerical examples are presented to illustrate the capability of the method.

preprint2021arXiv

Modeling silo clogging with nonlocal granular rheology

Granular flow in a silo demonstrates multiple nonlocal rheological phenomena due to the finite size of grains. We solve the Nonlocal Granular Fluidity (NGF) continuum model in quasi-2D silo geometries and evaluate its ability to predict these nonlocal effects, including flow spreading and, importantly, clogging (arrest) when the opening is small enough. The model is augmented to include a free-separation criterion and is implemented numerically with an extension of the trans-phase granular flow solver described in arXiv:1411.5447, to produce full-field solutions. The implementation is validated against analytical results of the model in the inclined chute geometry, such as the solution for the $H_{\mathrm{stop}}$ curve for size-dependent flow arrest, and the velocity profile as a function of layer height. We then implement the model in the silo geometry and vary the apparent grain size. The model predicts a jamming criterion when the opening competes with the scale of the mean grain size, which agrees with previous experimental studies, marking the first time to our knowledge that silo jamming has been achieved with a continuum model. For larger openings, the flow within the silo obtains a diffusive characteristic whose spread depends on the model's nonlocal amplitude and the mean grain size. The numerical tests are controlled for grid effects and a comparison study of coarse vs refined numerical simulations shows agreement in the pressure field, the shape of the arch in a jammed silo configuration, and the velocity field in a flowing configuration.

preprint2021arXiv

Surprising simplicity in the modeling of dynamic granular intrusion

Granular intrusions, such as dynamic impact or wheel locomotion, are complex multiphase phenomena where the grains exhibit solid-like and fluid-like characteristics together with an ejected gas-like phase. Despite decades of modeling efforts, a unified description of the physics in such intrusions is as yet unknown. Here we show that a continuum model based on the simple notions of frictional flow and tension-free separation describes complex granular intrusions near free surfaces. This model captures dynamics in a variety of experiments including wheel locomotion, plate intrusions, and running legged robots. The model reveals that three effects (a static contribution and two dynamic ones) primarily give rise to intrusion forces in such scenarios. Identification of these effects enables the development of a further reduced-order technique (Dynamic Resistive Force Theory) for rapid modeling of granular locomotion of arbitrarily shaped intruders. The continuum-motivated strategy we propose for identifying physical mechanisms and corresponding reduced-order relations has potential use for a variety of other materials.

preprint2020arXiv

Power-law scaling in granular rheology across flow geometries

Based on discrete element method simulations, we propose a new form of the constitution equation for granular flows independent of packing fraction. Rescaling the stress ratio $μ$ by a power of dimensionless temperature $Θ$ makes the data from a wide set of flow geometries collapse to a master curve depending only on the inertial number $I$. The basic power-law structure appears robust to varying particle properties (e.g. surface friction) in both 2D and 3D systems. We show how this rheology fits and extends frameworks such as kinetic theory and the Nonlocal Granular Fluidity model.

preprint2019arXiv

A general constitutive model for dense, fine particle suspensions validated in many geometries

Fine particle suspensions (such as cornstarch mixed with water) exhibit dramatic changes in viscosity when sheared, producing fascinating behaviors that captivate children and rheologists alike. Recent examination of these mixtures in simple flow geometries suggests inter-granular repulsion is central to this effect --- for mixtures at rest or shearing slowly, repulsion prevents frictional contacts from forming between particles, whereas, when sheared more forcefully, granular stresses overcome the repulsion allowing particles to interact frictionally and form microscopic structures that resist flow. Previous constitutive studies of these mixtures have focused on particular cases, typically limited to two-dimensional, steady, simple shearing flows. In this work, we introduce a predictive and general, three-dimensional continuum model for this material, using mixture theory to couple the fluid and particle phases. Playing a central role in the model, we introduce a micro-structural state variable, whose evolution is deduced from small-scale physical arguments and checked with existing data. Our space- and time-dependent model is implemented numerically in a variety of unsteady, non-uniform flow configurations where it is shown to accurately capture a variety of key behaviors: (i) the continuous shear thickening (CST) and discontinuous shear thickening (DST) behavior observed in steady flows, (ii) the time-dependent propagation of `shear jamming fronts', (iii) the time-dependent propagation of `impact activated jamming fronts', and (iv) the non-Newtonian, `running on oobleck' effect wherein fast locomotors stay afloat while slow ones sink.

preprint2016arXiv

Coupled dynamics of flow, microstructure, and conductivity in sheared suspensions

We propose a model for the evolution of the conductivity tensor for a flowing suspension of electrically conductive particles. We use discrete particle numerical simulations together with a continuum physical framework to construct an evolution law for the suspension microsutructure during flow. This model is then coupled with a relationship between the microstructure and the electrical conductivity tensor. The parameters of the joint model are fit experimentally using rheo- electrical conductivity measurements of carbon black suspensions under flow over a range of shear rates. The model is applied to the case of steady shearing as well as time-varying conductivity of unsteady flow experiments. We find that the model prediction agrees closely with the measured experimental data in all cases.

preprint2015arXiv

Continuum modelling and simulation of granular flows through their many phases

We propose and numerically implement a constitutive framework for granular media that allows the material to traverse through its many common phases during the flow process. When dense, the material is treated as a pressure sensitive elasto-viscoplastic solid obeying a yield criterion and a plastic flow rule given by the $μ(I)$ inertial rheology of granular materials. When the free volume exceeds a critical level, the material is deemed to separate and is treated as disconnected, stress-free media. A Material Point Method (MPM) procedure is written for the simulation of this model and many demonstrations are provided in different geometries. By using the MPM framework, extremely large strains and nonlinear deformations, which are common in granular flows, are representable. The method is verified numerically and its physical predictions are validated against known results.

preprint2015arXiv

Eulerian method for multiphase interactions of soft solid bodies in fluids

We introduce an Eulerian approach for problems involving one or more soft solids immersed in a fluid, which permits mechanical interactions between all phases. The reference map variable is exploited to simulate finite-deformation constitutive relations in the solid(s) on the same fixed grid as the fluid phase, which greatly simplifies the coupling between phases. Our coupling procedure, a key contribution in the current work, is shown to be computationally faster and more stable than an earlier approach, and admits the ability to simulate both fluid--solid and solid--solid interaction between submerged bodies. The interface treatment is demonstrated with multiple examples involving a weakly compressible Navier--Stokes fluid interacting with a neo-Hookean solid, and we verify the method's convergence. The solid contact method, which exploits distance-measures already existing on the grid, is demonstrated with two examples. A new, general routine for cross-interface extrapolation is introduced and used as part of the new interfacial treatment.

preprint2015arXiv

Intrusion in heterogeneous materials: Simple global rules from complex micro-mechanics

The interaction of intruding objects with deformable materials is a common phenomenon, arising in impact and penetration problems, animal and vehicle locomotion, and various geo-space applications. The dynamics of arbitrary intruders can be simplified using Resistive Force Theory (RFT), an empirical framework originally used for fluids but works surprisingly well, better in fact, in granular materials. That such a simple model describes behavior in dry grains, a complex nonlinear material, has invigorated a search to determine the underlying mechanism of RFT. We have discovered that a straightforward friction-based continuum model generates RFT, establishing a link between RFT and local material behavior. Our theory reproduces experimental RFT data without any parameter fitting and generates RFT's key simplifying assumption: a geometry-independent local force formula. Analysis of the system explains why RFT works better in grains than in viscous fluids, and leads to an analytical criterion to predict RFT's in other materials.

preprint2015arXiv

Modeling Growth Paths of Interacting Crack Pairs in Elastic Media

The problem of predicting the growth of a system of cracks, each crack influencing the growth of the others, arises in multiple fields. We develop an analytical framework toward this aim, which we apply to the `En-Passant' family of crack growth problems, in which a pair of initially parallel, offset cracks propagate nontrivially toward each other under far-field opening stress. We utilize boundary integral and perturbation methods of linear elasticity, Linear Elastic Fracture Mechanics, and common crack opening criteria to calculate the first analytical model for curved En-Passant crack paths. The integral system is reduced under a hierarchy of approximations, producing three methods of increasing simplicity for computing crack paths. The last such method is a major highlight of this work, using an asymptotic matching argument to predict crack paths based on superposition of simple, single-crack fields. Within the corresponding limits of the three methods, all three are shown to agree with each other. We provide comparisons to exact results and existing experimental data to verify certain approximation steps.

preprint2015arXiv

Modeling tensorial conductivity of particle suspension networks

Significant microstructural anisotropy is known to develop during shearing flow of attractive particle suspensions. These suspensions, and their capacity to form conductive networks, play a key role in flow-battery technology, among other applications. Herein, we present and test an analytical model for the tensorial conductivity of attractive particle suspensions. The model utilizes the mean fabric of the network to characterize the structure, and the relationship to the conductivity is inspired by a lattice argument. We test the accuracy of our model against a large number of computer-generated suspension networks, based on multiple in-house generation protocols, giving rise to particle networks that emulate the physical system. The model is shown to adequately capture the tensorial conductivity, both in terms of its invariants and its mean directionality.

preprint2014arXiv

Continuum modeling of mechanically-induced creep in dense granular materials

Recently, a new nonlocal granular rheology was successfully used to predict steady granular flows, including grain-size-dependent shear features, in a wide variety of flow configurations, including all variations of the split-bottom cell. A related problem in granular flow is that of mechanically-induced creep, in which shear deformation in one region of a granular medium fluidizes its entirety, including regions far from the sheared zone, effectively erasing the yield condition everywhere. This enables creep deformation when a force is applied in the nominally quiescent region through an intruder such as a cylindrical or spherical probe. We show that the nonlocal fluidity model is capable of capturing this phenomenology. Specifically, we explore creep of a circular intruder in a two-dimensional annular Couette cell and show that the model captures all salient features observed in experiments, including both the rate-independent nature of creep for sufficiently slow driving rates and the faster-than-linear increase in the creep speed with the force applied to the intruder.

preprint2014arXiv

Modeling the nonlocal behavior of granular flows down inclines

Flows of granular media down a rough inclined plane demonstrate a number of nonlocal phenomena. We apply the recently proposed nonlocal granular fluidity model to this geometry and find that the model captures many of these effects. Utilizing the model's dynamical form, we obtain a formula for the critical stopping height of a layer of grains on an inclined surface. Using an existing parameter calibration for glass beads, the theoretical result compares quantitatively to existing experimental data for glass beads. This provides a stringent test of the model, whose previous validations focused on driven steady-flow problems. For layers thicker than the stopping height, the theoretical flow profiles display a thickness-dependent shape whose features are in agreement with previous discrete particle simulations. We also address the issue of the Froude number of the flows, which has been shown experimentally to collapse as a function of the ratio of layer thickness to stopping height. While the collapse is not obvious, two explanations emerge leading to a revisiting of the history of inertial rheology, which the nonlocal model references for its homogeneous flow response.

preprint2014arXiv

Symmetry Relations in Viscoplastic Drag Laws

The following note shows that the symmetry of various resistance formulae, often based on Lorentz reciprocity for linearly viscous fluids, applies to a wide class of non-linear viscoplastic fluids. This follows from Edelen's non-linear generalization of the Onsager relation for the special case of \emph{strongly dissipative} rheology, where constitutive equations are derivable from his dissipation potential. For flow domains with strong dissipation in the interior and on a portion of the boundary this implies strong dissipation on the remaining portion of the boundary, with strongly dissipative traction-velocity response given by a dissipation potential. This leads to a non-linear generalization of Stokes resistance formulae for a wide class of viscoplastic fluid problems. We consider the application to non-linear Darcy flow and to the effective slip for viscoplastic flow over textured surfaces.

preprint2014arXiv

Two-temperature continuum thermomechanics of deforming amorphous solids

There is an ever-growing need for predictive models for the elasto-viscoplastic deformation of solids. Our goal in this paper is to incorporate recently developed out-of-equilibrium statistical concepts into a thermodynamically consistent, finite-deformation, continuum framework for deforming amorphous solids. The basic premise is that the configurational degrees of freedom of the material --- the part of the internal energy/entropy that corresponds to mechanically stable microscopic configurations --- are characterized by a configurational temperature that might differ from that of the vibrational degrees of freedom, which equilibrate rapidly with an external heat bath. This results in an approximate internal energy decomposition into weakly interacting configurational and vibrational subsystems, which exchange energy following a Fourier-like law, leading to a thermomechanical framework permitting two well-defined temperatures. In this framework, internal variables that carry information about the state of the material equilibrate with the configurational subsystem, are explicitly associated with energy and entropy of their own, and couple to a viscoplastic flow rule. The coefficients that determine the rate of flow of entropy and heat between different internal systems are proposed to explicitly depend on the rate of irreversible deformation. As an application of this framework, we discuss two constitutive models for the response of glassy materials, a simple phenomenological model and a model related to the concept of Shear-Transformation-Zones as the basis for internal variables. The models account for several salient features of glassy deformation phenomenology. Directions for future investigation are briefly discussed.

preprint2010arXiv

Effective slip boundary conditions for arbitrary periodic surfaces: The surface mobility tensor

In a variety of applications, most notably microfluidic design, slip-based boundary conditions have been sought to characterize fluid flow over patterned surfaces. We focus on laminar shear flows over surfaces with periodic height fluctuations and/or fluctuating Navier scalar slip properties. We derive a general formula for the "effective slip", which describes equivalent fluid motion at the mean surface as depicted by the linear velocity profile that arises far from it. We show that the slip and the applied stress are related linearly through a tensorial mobility matrix, and the method of domain perturbation is then used to derive an approximate formula for the mobility law directly in terms of surface properties. The specific accuracy of the approximation is detailed, and the mobility relation is then utilized to address several questions, such as the determination of optimal surface shapes and the effect of random surface fluctuations on fluid slip.

preprint2006arXiv

A stochastic flow rule for granular materials

There have been many attempts to derive continuum models for dense granular flow, but a general theory is still lacking. Here, we start with Mohr-Coulomb plasticity for quasi-2D granular materials to calculate (average) stresses and slip planes, but we propose a "stochastic flow rule" (SFR) to replace the principle of coaxiality in classical plasticity. The SFR takes into account two crucial features of granular materials - discreteness and randomness - via diffusing "spots" of local fluidization, which act as carriers of plasticity. We postulate that spots perform random walks biased along slip-lines with a drift direction determined by the stress imbalance upon a local switch from static to dynamic friction. In the continuum limit (based on a Fokker-Planck equation for the spot concentration), this simple model is able to predict a variety of granular flow profiles in flat-bottom silos, annular Couette cells, flowing heaps, and plate-dragging experiments -- with essentially no fitting parameters -- although it is only expected to function where material is at incipient failure and slip-lines are inadmissible. For special cases of admissible slip-lines, such as plate dragging under a heavy load or flow down an inclined plane, we postulate a transition to rate-dependent Bagnold rheology, where flow occurs by sliding shear planes. With different yield criteria, the SFR provides a general framework for multiscale modeling of plasticity in amorphous materials, cycling between continuum limit-state stress calculations, meso-scale spot random walks, and microscopic particle relaxation.