BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Richard M. Lueptow

Richard M. Lueptow appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softphysics.flu-dynmath.DSphysics.chem-phphysics.comp-phnlin.CD

Catalog footprint

What is connected

17works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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)

preprint2022arXiv

Designing minimally-segregating granular mixtures for gravity-driven surface flows

In dense flowing bidisperse particle mixtures varying in size or density alone, smaller particles sink (driven by percolation) and lighter particles rise (driven by buoyancy). But when the particle species differ from each other in both size and density, percolation and buoyancy can either enhance (large/light and small/heavy) or oppose (large/heavy and small/light) each other. In the latter case, a local equilibrium condition can exist in which the two segregation mechanisms balance and particles remain mixed: this allows the design of minimally-segregating mixtures by specifying particle size ratio, density ratio, and mixture concentration. Using experimentally validated DEM simulations, we show that mixtures specified by the methodology remain relatively well-mixed in the thin rapid surface flows characteristic of heaps and tumblers commonly used in industry. Furthermore, minimally-segregating particle mixtures prepared in a fully segregated state in a tumbler mix over time and eventually reach a state of nearly uniform species concentration.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modeling Stratified Segregation in Periodically Driven Granular Heap Flow

We present a continuum approach to model segregation of size-bidisperse granular materials in unsteady bounded heap flow as a prototype for modeling segregation in other time varying flows. In experiments, a periodically modulated feed rate produces stratified segregation like that which occurs due to intermittent avalanching, except with greater layer-uniformity and higher average feed rates. Using an advection-diffusion-segregation equation and characterizing transient changes in deposition and erosion after a feed rate change, we model stratification for varying feed rates and periods. Feed rate modulation in heap flows can create well-segregated layers, which effectively mix the deposited material normal to the free surface at lengths greater than the combined layer-thickness. This mitigates the strong streamwise segregation that would otherwise occur at larger particle-size ratios and equivalent steady feed rates and can significantly reduce concentration variation during hopper discharge. Coupling segregation, deposition and erosion is challenging but has many potential applications.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Segregation forces in dense granular flows: Closing the gap from single intruders to mixtures

Using simulations and a virtual-spring-based approach, we measure the segregation force, Fseg, over a range of size-bidisperse mixture concentrations, particle size ratios, and shear rates to develop a model for Fseg that extends its applicability from the well-studied non-interacting intruders regime to finite-concentration mixtures where cooperative phenomena occur. The model predicts the concentration below which the single intruder assumption applies and provides an accurate description of the pressure partitioning between species.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring segregation characteristics of industrially relevant granular mixtures: Part I -- A continuum model approach

We present a method to estimate the segregation parameter, $S,$ a key input in a continuum transport model of particulate flows. $S$ is determined by minimizing the difference between measured and model-predicted concentration profiles. To validate the approach, we conduct discrete element method simulations of size-bidisperse mixtures in quasi-2D bounded heap flow; the resulting data show that $S$ calculated from concentration profiles is consistent with the directly measured value. The method's accuracy depends critically on the velocity profile during filling, but only weakly on the diffusion coefficient. When the velocity profile is nominally spanwise invariant, the error between estimated and measured $S$ is $10\%$. This method is intended for practical application (described in Part II), so we restrict characterization of the velocity profile to that which can be readily determined experimentally, and explore the sensitivity of concentration profiles to variation of the gap between the sidewalls of the heap.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Identifying Invariant Ergodic Subsets and Barriers to Mixing by Cutting and Shuffling: Study in a Bi-rotated Hemisphere

Mixing by cutting-and-shuffling can be mathematically described by the dynamics of piecewise isometries (PWIs), higher dimensional analogs of one-dimensional interval exchange transformations. In a two-dimensional domain under a PWI, the exceptional set, $\bar{E}$, which is created by the accumulation of cutting lines (the union of all iterates of cutting lines and all points that pass arbitrarily close to a cutting line), defines where mixing is possible but not guaranteed. There is structure within $\bar{E}$ that directly influences the mixing potential of the PWI. Here we provide new computational and analytical formalisms for examining this structure by way of measuring the density and connectivity of $\varepsilon$-fattened cutting lines that form an approximation of $\bar{E}$. For the example of a PWI on a hemispherical shell studied here, this approach reveals the subtle mixing behaviors and barriers to mixing formed by invariant ergodic subsets (confined orbits) within the fractal structure of the exceptional set. Some PWIs on the shell have provably non-ergodic exceptional sets, which prevent mixing, while others have potentially ergodic exceptional sets where mixing is possible since ergodic exceptional sets have uniform cutting line density. For these latter exceptional sets, we show the connectivity of orbits in the PWI map through direct examination of orbit position and shape and through a two-dimensional return plot to explain the necessity of orbit connectivity for mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Modified Archimedes' principle predicts rising and sinking of intruders in sheared granular flows

We computationally determine the force on single spherical intruder particles in sheared granular flows as a function of particle size, particle density, shear rate, overburden pressure, and gravitational acceleration. The force scales similarly to, but deviates from, the buoyancy force predicted by Archimedes' principle. The deviation depends only on the intruder to bed particle size ratio, but not the density ratio or flow conditions. We propose a simple force model that successfully predicts whether intruders rise or sink, knowing only the size and density ratios, for a variety of flow configurations in physical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of water and solute transport in polymeric reverse osmosis membranes via molecular dynamics simulations

The Angstrom-scale transport characteristics of water and six different solutes, methanol, ethanol, 2-propanol, urea, Na+, and Cl-, were studied for a polyamide reverse osmosis (RO) membrane, FT-30, using non-equilibrium molecular dynamics (NEMD) simulations. Results indicate that water transport increases with an increasing fraction of connected percolated free volume, or water-accessible open space, in the membrane polymer structure. This free volume is enhanced by the dynamic structure of the membrane at the molecular level as it swells when hydrated and vibrates due to molecular collisions allowing a continuous path connecting the opposite membrane surfaces. The tortuous paths available for transport of solutes result in Brownian motion of solute molecules and hopping from pore to pore as they pass through the polymer network structure of the membrane. The transport of alcohol solutes decreases for solutes with larger Van der Waals volume, which corresponds to less available percolated free volume, or solute-accessible space, within the membrane polymer structure. However, the Van der Waals size of the dehydrated solutes is generally not a good measure to predict solute transport or rejection. Urea has reduced transport compared to ethanol, most likely due to more complex chemistry, even though urea has a smaller Van der Waals volume than ethanol. Na+ and Cl- experience the lowest transport, likely due to strong ion-water and ion-ion electrostatic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rejection mechanisms for contaminants in polymeric reverse osmosis membranes

Despite the success of reverse osmosis (RO) for water purification, the molecular-level physico-chemical processes of contaminant rejection are not well understood. Here we carry out NEMD simulations on a model polyamide RO membrane to understand the mechanisms of transport and rejection of both ionic and neutral contaminants in water. We observe that the rejection changes non-monotonously with ion sizes. In particular, the rejection of urea, 2.4 A radius, is higher than ethanol, 2.6 A radius, and the rejections for organic solutes, 2.2-2.8 A radius, are lower than Na+, 1.4 A radius, or Cl-, 2.3 A radius. We show that this can be explained in terms of the solute accessible intermolecular volume in the membrane and the solute-water pair interaction energy. If the smallest open spaces in the membrane's molecular structure are all larger than the hydrated solute, then the solute-water pair interaction energy does not matter. However, when the open spaces in the polymeric structure are such that solutes have to shed at least one water molecule to pass through a portion of the membrane molecular structure, the pair interaction energy governs solute rejection. The high pair interaction energy for water molecules in the solvation shell for ions makes the water molecules difficult to shed, thus enhancing the rejection of ions. On the other hand, the organic solute-water interaction energies are governed by the water molecules that are hydrogen bonded to the solute. Urea molecules have more hydrogen-bonding sites than alcohol molecules, leading to a higher rejection of urea than occurs for ethanol, a molecule of similar size but with fewer hydrogen bonding sites. These findings underline the importance of the solute's solvation shell and solute-water-membrane chemistry in the context of reverse osmosis, thus providing new insights into solute transport and rejection in RO membranes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Parametric Study of Mixing in a Granular Flow a Bi-Axial Spherical Tumbler

We report on a computational parameter space study of mixing protocols for a half-full bi-axial spherical granular tumbler. The quality of mixing is quantified via the intensity of segregation (concentration variance) and computed as a function of three system parameters: angles of rotation about each tumbler axis and the flowing layer depth. Only the symmetric case is considered in which the flowing layer depth is the same for each rotation. We also consider the dependence on $\bar{R}$, which parametrizes the concentric spheroids ("shells") that comprise the volume of the tumbler. The intensity of segregation is computed over 100 periods of the mixing protocol for each choice of parameters. Each curve is classified via a time constant, $τ$, and an asymptotic mixing value, $bias$. We find that most choices of angles and most shells throughout the tumbler volume mix well, with mixing near the center of the tumbler being consistently faster (small $τ$) and more complete (small $bias$). We conclude with examples and discussion of the pathological mixing behaviors of the outliers in the so-called $τ$-$bias$ scatterplots.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Influence of Rough and Smooth Walls on Macroscale Flows in Tumblers

Walls in discrete element method simulations of granular flows are sometimes modeled as a closely packed monolayer of fixed particles, resulting in a rough wall rather than a geometrically smooth wall. An implicit assumption is that the resulting rough wall differs from a smooth wall only locally at the particle scale. Here we test this assumption by considering the impact of the wall roughness at the periphery of the flowing layer on the flow of monodisperse particles in a rotating spherical tumbler. We find that varying the wall roughness significantly alters average particle trajectories even far from the wall. Rough walls induce greater poleward axial drift of particles near the flowing layer surface, but decrease the curvature of the trajectories. Increasing the volume fill level in the tumbler has little effect on the axial drift for rough walls, but increases the drift while reducing curvature of the particle trajectories for smooth walls. The mechanism for these effects is related to the degree of local slip at the bounding wall, which alters the flowing layer thickness near the walls, affecting the particle trajectories even far from the walls near the equator of the tumbler. Thus, the proper choice of wall conditions is important in the accurate simulation of granular flows, even far from the bounding wall.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Influence of Rough and Smooth Walls on Macroscale Granular Segregation Patterns

Size bidisperse granular materials in a spherical tumbler segregate into two different patterns of three bands with either small particles at the equator and large particles at the poles or vice versa, depending upon the fill level in the tumbler. Here we use discrete element method (DEM) simulations with supporting qualitative experiments to explore the effect of the tumbler wall roughness on the segregation pattern, modeling the tumbler walls as either a closely packed monolayer of fixed particles resulting in a rough wall, or as a geometrically smooth wall. Even though the tumbler wall is in contact with the flowing layer only at its periphery, the impact of wall roughness is profound. Smooth walls tend toward a small-large-small (SLS) band pattern at the pole-equator-pole at all but the highest fill fractions; rough walls tend toward a large-small-large (LSL) band pattern at all but the lowest fill fractions. This comes about because smooth walls induce poleward axial drift of small particles and an equator-directed drift for large particles, resulting in an SLS band pattern. On the other hand, rough walls result in both sizes of particles moving poleward at the surface of the flow, but due to radial segregation, small particles percolate lower in the flowing layer where there is a return drift toward the equator while large particles remain at the surface near the pole, resulting in an LSL band pattern. The tendency toward either of the two band patterns depends on the fill level in the tumbler and the roughness of the tumbler's bounding wall.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Study in Three-Dimensional Chaotic Dynamics: Granular Flow and Transport in a Bi-Axial Spherical Tumbler

We study 3D chaotic dynamics through an analysis of transport in a granular flow in a half-full spherical tumbler rotated sequentially about two orthogonal axes (a bi-axial "blinking" tumbler). The flow is essentially quasi-2D in any vertical slice of the sphere during rotation about a single axis, and we provide an explicit exact solution to the model in this case. Hence, the cross-sectional flow can be represented by a twist map, allowing us to express the 3D flow as a linked twist map (LTM). We prove that if the rates of rotation about each axis are equal, then (in the absence of stochasticity) particle trajectories are restricted to 2D surfaces consisting of a portion of a hemispherical shell closed by a "cap"; if the rotation rates are unequal, then particles can leave the surface they start on and traverse a volume of the tumbler. The period-one structures of the governing LTM are examined in detail: analytical expressions are provided for the location of period-one curves, their extent into the bulk of the granular material, and their dependence on the protocol parameters (rates and durations of rotations). Exploiting the restriction of trajectories to 2D surfaces in the case of equal rotation rates about the axes, a method is proposed for identifying and constructing 3D Kolmogorov-Arnold-Moser (KAM) tubes around the normally elliptic period-one curves. The invariant manifold structure arising from the normally hyperbolic period-one curves is also examined. When the motion is restricted to 2D surfaces, the structure of manifolds of the hyperbolic points in the bulk differs from that corresponding to hyperbolic points in the flowing layer. Each is reminiscent of a template provided by a non-integrable perturbation to a Hamiltonian system, though the governing LTM is not. This highlights the novel 3D chaotic behaviors observed in this model dynamical system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling size segregation of granular materials: the roles of segregation, advection and diffusion

Predicting segregation of granular materials composed of different-sized particles is a challenging problem. In this paper, we develop and implement a theoretical model that captures the interplay between advection, segregation, and diffusion in size bidisperse granular materials. The fluxes associated with these three driving factors depend on the underlying kinematics, whose characteristics play key roles in determining particle segregation configurations. Unlike previous models for segregation, our model uses parameters based on kinematic measures from discrete element method simulations instead of arbitrarily adjustable fitting parameters, and it achieves excellent quantitative agreement with both experimental and simulation results when applied to quasi-two-dimensional bounded heaps. The model yields two dimensionless control parameters, both of which are only functions of physically control parameters (feed rate, particle sizes, and system size) and kinematic parameters (diffusion coefficient, flowing layer depth, and percolation velocity). The Péclet number, $Pe$, captures the interplay of advection and diffusion, and the second dimensionless parameter, $Λ$, describes the interplay between segregation and advection. A parametric study of $Λ$ and $Pe$ demonstrates how the particle segregation configuration depends on the interplay of advection, segregation, and diffusion. The model can be readily adapted to other flow geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Slow axial drift in three-dimensional granular tumbler flow

Models of monodisperse particle flow in partially filled three-dimensional tumblers often assume that flow along the axis of rotation is negligible. We test this assumption, for spherical and double cone tumblers, using experiments and discrete element method simulations. Cross sections through the particle bed of a spherical tumbler show that, after a few rotations, a colored band of particles initially perpendicular to the axis of rotation deforms: particles near the surface drift toward the pole, while particles deeper in the flowing layer drift toward the equator. Tracking of mm-sized surface particles in tumblers with diameters of 8-14 cm shows particle axial displacements of one to two particle diameters, corresponding to axial drift that is 1-3% of the tumbler diameter, per pass through the flowing layer. The surface axial drift in both double cone and spherical tumblers is zero at the equator, increases moving away from the equator, and then decreases near the poles. Comparing results for the two tumbler geometries shows that wall slope causes axial drift, while drift speed increases with equatorial diameter. The dependence of axial drift on axial position for each tumbler geometry is similar when both are normalized by their respective maximum values.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cutting and Shuffling a Line Segment: Mixing by Interval Exchange Transformations

We present a computational study of finite-time mixing of a line segment by cutting and shuffling. A family of one-dimensional interval exchange transformations is constructed as a model system in which to study these types of mixing processes. Illustrative examples of the mixing behaviors, including pathological cases that violate the assumptions of the known governing theorems and lead to poor mixing, are shown. Since the mathematical theory applies as the number of iterations of the map goes to infinity, we introduce practical measures of mixing (the percent unmixed and the number of intermaterial interfaces) that can be computed over given (finite) numbers of iterations. We find that good mixing can be achieved after a finite number of iterations of a one-dimensional cutting and shuffling map, even though such a map cannot be considered chaotic in the usual sense and/or it may not fulfill the conditions of the ergodic theorems for interval exchange transformations. Specifically, good shuffling can occur with only six or seven intervals of roughly the same length, as long as the rearrangement order is an irreducible permutation. This study has implications for a number of mixing processes in which discontinuities arise either by construction or due to the underlying physics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stratification, segregation and mixing of granular materials in quasi-2D bounded heaps

Segregation and mixing of granular mixtures during heap formation have important consequences in industry and agriculture. This research investigates three different final particle configurations of bi-disperse granular mixtures -- stratified, segregated and mixed -- during filling of quasi-two dimensional silos. We consider a larger number and relatively wider range of control parameters than previous studies, including particle size ratio, flow rate, system size and heap rise velocity. The boundary between stratified and unstratified states is primarily controlled by the two-dimensional flow rate, with the critical flow rate for the transition depending weakly on particle size ratio and flowing layer length. In contrast, the transition from segregated to mixed states is controlled by the rise velocity of the heap, a control parameter not previously considered. The critical rise velocity for the transition depends strongly on the particle size ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stretching and folding versus cutting and shuffling: An illustrated perspective on mixing and deformations of continua

We compare and contrast two types of deformations inspired by mixing applications -- one from the mixing of fluids (stretching and folding), the other from the mixing of granular matter (cutting and shuffling). The connection between mechanics and dynamical systems is discussed in the context of the kinematics of deformation, emphasizing the equivalence between stretches and Lyapunov exponents. The stretching and folding motion exemplified by the baker's map is shown to give rise to a dynamical system with a positive Lyapunov exponent, the hallmark of chaotic mixing. On the other hand, cutting and shuffling does not stretch. When an interval exchange transformation is used as the basis for cutting and shuffling, we establish that all of the map's Lyapunov exponents are zero. Mixing, as quantified by the interfacial area per unit volume, is shown to be exponentially fast when there is stretching and folding, but linear when there is only cutting and shuffling. We also discuss how a simple computational approach can discern stretching in discrete data.

0 citations0 reviewsNo code linkedOpen access