Source author record

Chris H. Rycroft

Chris H. Rycroft 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

15works
14topics
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

15 published item(s)

preprint2022arXiv

A projection method for porous media flow

Flow through porous, elastically deforming media is present in a variety of natural contexts ranging from large-scale geophysics to cellular biology. In the case of incompressible constituents, the porefluid pressure acts as a Lagrange multiplier to satisfy the resulting constraint on fluid divergence. The resulting system of equations is a possibly non-linear saddle-point problem and difficult to solve numerically, requiring nonlinear implicit solvers or flux-splitting methods. Here, we present a method for the simulation of flow through porous media and its coupled elastic deformation. The pore pressure field is calculated at each time step by correcting trial velocities in a manner similar to Chorin projection methods. We demonstrate the method's second-order convergence in space and time and show its application to phase separating neo-Hookean gels.

preprint2022arXiv

Voronoi cell analysis: The shapes of particle systems

Many physical systems can be studied as collections of particles embedded in space, evolving through deterministic evolution equations. Natural questions arise concerning how to characterize these arrangements - are they ordered or disordered? If they are ordered, how are they ordered and what kinds of defects do they possess? Originally introduced to study problems in pure mathematics, Voronoi tessellations have become a powerful and versatile tool for analyzing countless problems in pure and applied physics. In this paper we explain the basics of Voronoi tessellations and the shapes they produce, and describe how they can be used to study many physical systems.

preprint2021arXiv

A model for the fragmentation kinetics of crumpled thin sheets

As a confined thin sheet crumples, it spontaneously segments into flat facets delimited by a network of ridges. Despite the apparent disorder of this process, statistical properties of crumpled sheets exhibit striking reproducibility. Experiments have shown that the total crease length accrues logarithmically when repeatedly compacting and unfolding a sheet of paper. Here, we offer insight to this unexpected result by exploring the correspondence between crumpling and fragmentation processes. We identify a physical model for the evolution of facet area and ridge length distributions of crumpled sheets, and propose a mechanism for re-fragmentation driven by geometric frustration. This mechanism establishes a feedback loop in which the facet size distribution informs the subsequent rate of fragmentation under repeated confinement, thereby producing a new size distribution. We then demonstrate the capacity of this model to reproduce the characteristic logarithmic scaling of total crease length, thereby supplying a missing physical basis for the observed phenomenon.

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

Computing the viscous effect in early-time drop impact dynamics

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. (2014b) show that the liquid-air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier-Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square root scaling relationship. We also predict the effect of interfacial tension at the liquid-gas interface on the drop impact, showing that it mediates the lift-off behavior.

preprint2021arXiv

Eulerian simulation of complex suspensions and biolocomotion in three dimensions

We present a numerical method specifically designed for simulating three-dimensional fluid--structure interaction (FSI) problems based on the reference map technique (RMT). The RMT is a fully Eulerian FSI numerical method that allows fluids and large-deformation elastic solids to be represented on a single fixed computational grid. This eliminates the need for meshing complex geometries typical in other FSI approaches, and greatly simplifies the coupling between fluid and solids. We develop the first three-dimensional implementation of the RMT, parallelized using the distributed memory paradigm, to simulate incompressible FSI with neo-Hookean solids. As part of our new method, we develop a new field extrapolation scheme that works efficiently in parallel. Through representative examples, we demonstrate the method's accuracy and convergence, as well as its suitability in investigating many-body and active systems. The examples include settling of a mixture of heavy and buoyant soft ellipsoids, lid-driven cavity flow containing a soft sphere, and swimmers actuated via active stress.

preprint2021arXiv

Manifold learning for coarse-graining atomistic simulations: Application to amorphous solids

We introduce a generalized machine learning framework to probabilistically parameterize upper-scale models in the form of nonlinear PDEs consistent with a continuum theory, based on coarse-grained atomistic simulation data of mechanical deformation and flow processes. The proposed framework utilizes a hypothesized coarse-graining methodology with manifold learning and surrogate-based optimization techniques. Coarse-grained high-dimensional data describing quantities of interest of the multiscale models are projected onto a nonlinear manifold whose geometric and topological structure is exploited for measuring behavioral discrepancies in the form of manifold distances. A surrogate model is constructed using Gaussian process regression to identify a mapping between stochastic parameters and distances. Derivative-free optimization is employed to adaptively identify a unique set of parameters of the upper-scale model capable of rapidly reproducing the system's behavior while maintaining consistency with coarse-grained atomic-level simulations. The proposed method is applied to learn the parameters of the shear transformation zone (STZ) theory of plasticity that describes plastic deformation in amorphous solids as well as coarse-graining parameters needed to translate between atomistic and continuum representations. We show that the methodology is able to successfully link coarse-grained microscale simulations to macroscale observables and achieve a high-level of parity between the models across scales.

preprint2017arXiv

Active elastohydrodynamics of vesicles in narrow, blind constrictions

Fluid-resistance limited transport of vesicles through narrow constrictions is a recurring theme in many biological and engineering applications. Inspired by the motor-driven movement of soft membrane-bound vesicles into closed neuronal dendritic spines, here we study this problem using a combination of passive three-dimensional simulations and a simplified semi-analytical theory for active transport of vesicles that are forced through such constrictions by molecular motors. We show that the motion of these objects is characterized by two dimensionless quantities related to the geometry and the strength of forcing relative to the vesicle elasticity. We use numerical simulations to characterize the transit time for a vesicle forced by fluid pressure through a constriction in a channel, and find that relative to an open channel, transport into a blind end leads to the formation of an effective lubrication layer that strongly impedes motion. When the fluid pressure forcing is complemented by forces due to molecular motors that are responsible for vesicle trafficking into dendritic spines, we find that the competition between motor forcing and fluid drag results in multistable dynamics reminiscent of the real system. Our study highlights the role of non-local hydrodynamic effects in determining the kinetics of vesicular transport in constricted geometries.

preprint2017arXiv

Density-equalizing maps for simply-connected open surfaces

In this paper, we are concerned with the problem of creating flattening maps of simply-connected open surfaces in $\mathbb{R}^3$. Using a natural principle of density diffusion in physics, we propose an effective algorithm for computing density-equalizing flattening maps with any prescribed density distribution. By varying the initial density distribution, a large variety of mappings with different properties can be achieved. For instance, area-preserving parameterizations of simply-connected open surfaces can be easily computed. Experimental results are presented to demonstrate the effectiveness of our proposed method. Applications to data visualization and surface remeshing are explored.

preprint2016arXiv

Notch fracture toughness of glasses: Rate, age and geometry dependence

Understanding the fracture toughness (resistance) of glasses is a fundamental problem of prime theoretical and practical importance. Here we theoretically study its dependence on the loading rate, the age (history) of the glass and the notch radius $ρ$. Reduced-dimensionality analysis suggests that the notch fracture toughness results from a competition between the initial, age- and history-dependent, plastic relaxation timescale $τ^{pl}_0$ and an effective loading timescale $τ^{ext}(\dot{K}_I,ρ)$, where $\dot{K}_I$ is the tensile stress-intensity-factor rate. The toughness is predicted to scale with $\sqrtρ$ independently of $ξ\!\equiv\!τ^{ext}\!/τ^{pl}_0$ for $ξ\!\ll\! 1$, to scale as $T\sqrtρ\,\log(ξ)$ for $ξ\!\gg\!1$ (related to thermal activation, where $T$ is the temperature) and to feature a non-monotonic behavior in the crossover region $ξ\!\sim\!{\cal O}(1)$ (related to plastic yielding dynamics). These predictions are verified using novel 2D computations, providing a unified picture of the notch fracture toughness of glasses. The theory highlights the importance of timescales competition and far from steady-state elasto-viscoplastic dynamics for understanding the toughness, and shows that the latter varies quite significantly with the glass age (history) and applied loading rate. Experimental support for bulk metallic glasses is presented.

preprint2016arXiv

Optimizing intermittent water supply in urban pipe distribution networks

In many urban areas of the developing world, piped water is supplied only intermittently, as valves direct water to different parts of the water distribution system at different times. The flow is transient, and may transition between free-surface and pressurized, resulting in complex dynamical features with important consequences for water suppliers and users. Here, we develop a computational model of transition, transient pipe flow in a network, accounting for a wide variety of realistic boundary conditions. We validate the model against several published data sets, and demonstrate its use on a real pipe network. The model is extended to consider several optimization problems motivated by realistic scenarios. We demonstrate how to infer water flow in a small pipe network from a single pressure sensor, and show how to control water inflow to minimize damaging pressure gradients.

preprint2015arXiv

An Eulerian projection method for quasi-static elastoplasticity

A well-established numerical approach to solve the Navier--Stokes equations for incompressible fluids is Chorin's projection method, whereby the fluid velocity is explicitly updated, and then an elliptic problem for the pressure is solved, which is used to orthogonally project the velocity field to maintain the incompressibility constraint. In this paper, we develop a mathematical correspondence between Newtonian fluids in the incompressible limit and hypo-elastoplastic solids in the slow, quasi-static limit. Using this correspondence, we formulate a new fixed-grid, Eulerian numerical method for simulating quasi-static hypo-elastoplastic solids, whereby the stress is explicitly updated, and then an elliptic problem for the velocity is solved, which is used to orthogonally project the stress to maintain the quasi-staticity constraint. We develop a finite-difference implementation of the method and apply it to an elasto-viscoplastic model of a bulk metallic glass based on the shear transformation zone theory. We show that in a two-dimensional plane strain simple shear simulation, the method is in quantitative agreement with an explicit method. Like the fluid projection method, it is efficient and numerically robust, making it practical for a wide variety of applications. We also demonstrate that the method can be extended to simulate objects with evolving boundaries. We highlight a number of correspondences between incompressible fluid mechanics and quasi-static elastoplasticity, creating possibilities for translating other numerical methods between the two classes of physical problems.

preprint2015arXiv

Asymmetric collapse by dissolution or melting in a uniform flow

An advection--diffusion-limited dissolution model of an object being eroded by a two-dimensional potential flow is presented. By taking advantage of the conformal invariance of the model, a numerical method is introduced that tracks the evolution of the object boundary in terms of a time-dependent Laurent series. Simulations of a variety of dissolving objects are shown, which shrink and then collapse to a single point in finite time. The simulations reveal a surprising exact relationship whereby the collapse point is the root of a non-analytic function given in terms of the flow velocity and the Laurent series coefficients describing the initial shape. This result is subsequently derived using residue calculus. The structure of the non-analytic function is examined for three different test cases, and a practical approach to determine the collapse point using a generalized Newton--Raphson root-finding algorithm is outlined. These examples also illustrate the possibility that the model breaks down in finite time prior to complete collapse, due to a topological singularity, as the dissolving boundary overlaps itself rather than breaking up into multiple domains (analogous to droplet pinch-off in fluid mechanics). In summary, the model raises fundamental mathematical questions about broken symmetries in finite-time singularities of both continuous and stochastic dynamical systems.

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.

preprint2012arXiv

Fracture Toughness of Metallic Glasses: Ductile-to-Brittle Transition?

Quantitative understanding of the fracture toughness of metallic glasses, including the associated ductile-to-brittle transitions, is not yet available. Here we use a simple model of plastic deformation in glasses, coupled to an advanced Eulerian level set formulation for solving complex free boundary problems, to calculate the fracture toughness of metallic glasses as a function of the degree of structural relaxation corresponding to different annealing times near the glass temperature. Our main result indicates the existence of an elasto-plastic crack tip instability for sufficiently relaxed glasses, resulting in a marked drop in the toughness, which we interpret as a ductile-to-brittle transition similar to experimental observations.