Source author record

Renato Feres

Renato Feres 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

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

11 published item(s)

preprint2020arXiv

Entropy Production in Random Billiards

We introduce a class of random mechanical systems called random billiards to study the problem of quantifying the irreversibility of nonequilibrium macroscopic systems. In a random billiard model, a point particle evolves by free motion through the interior of a spatial domain, and reflects according to a reflection operator, specified in the model by a Markov transition kernel, upon collision with the boundary of the domain. We derive a formula for entropy production rate that applies to a general class of random billiard systems. This formula establishes a relation between the purely mathematical concept of entropy production rate and textbook thermodynamic entropy, recovering in particular Clausius' formulation of the second law of thermodynamics. We also study an explicit class of examples whose reflection operator, referred to as the Maxwell-Smolukowski thermostat, models systems with boundary thermostats kept at possibly different temperatures. We prove that, under certain mild regularity conditions, the class of models are uniformly ergodic Markov chains and derive formulas for the stationary distribution and entropy production rate in terms of geometric and thermodynamic parameters.

preprint2018arXiv

Rolling and no-slip bouncing in cylinders

The purpose of this paper is to compare a classical non-holonomic system---a sphere rolling against the inner surface of a vertical cylinder under gravity---and a class of discrete dynamical systems known as no-slip billiards in similar configurations. A well-known notable feature of the non-holonomic system is that the rolling sphere does not fall; its height function is bounded and oscillates harmonically up and down. The central issue of the present work is whether similar bounded behavior can be observed in the no-slip billiard counterpart. Our main results are as follows: for circular cylinders in dimension $3$, the no-slip billiard has the bounded orbits property, and very closely approximates rolling motion, for a class of initial conditions which we call transversal rolling impact. When this condition does not hold, trajectories undergo vertical oscillations superimposed to an overall downward acceleration. Considering cylinders with different cross-section shapes, we show that no-slip billiards between two parallel hyperplanes in Euclidean space of arbitrary dimension are always bounded even under a constant force parallel to the plates; for general cylinders, when the orbit of the transverse system (a concept that depends on a factorization of the motion into transversal and longitudinal components) has period two---a very common occurrence in planar no-slip billiards---the motion in the longitudinal direction, under no forces, is generically not bounded. This is shown using a formula for a longitudinal linear drift that we prove in arbitrary dimensions. While the systems for which we can prove the existence of bounded orbits have relatively simple transverse dynamics, we also briefly explore numerically a no-slip billiard system, namely the stadium cylinder billiard, that can exhibit chaotic transversal dynamics.

preprint2016arXiv

No-slip billiards in dimension two

We investigate the dynamics of no-slip billiards, a model in which small rotating disks may exchange linear and angular momentum at collisions with the boundary. We give new results on periodicity and boundedness of orbits which suggest that a class of billiards (including all polygons) is not ergodic. Computer generated phase portraits demonstrate non-ergodic features, suggesting chaotic no-slip billiards cannot readily be constructed using the common techniques for generating chaos in standard billiards.

preprint2016arXiv

Stability of periodic orbits in no-slip billiards

Rigid bodies collision maps in dimension two, under a natural set of physical requirements, can be classified into two types: the standard specular reflection map and a second which we call, after Broomhead and Gutkin, no-slip. This leads to the study of no-slip billiards--planar billiard systems in which the moving particle is a disc (with rotationally symmetric mass distribution) whose translational and rotational velocities can both change at each collision with the boundary of the billiard domain. In this paper we greatly extend previous results on boundedness of orbits (Broomhead and Gutkin) and linear stability of periodic orbits for a Sinai-type billiard (Wojtkowski) for no-slip billiards. We show among other facts that: (i) for billiard domains in the plane having piecewise smooth boundary and at least one corner of inner angle less than $π$, no-slip billiard dynamics will always contain elliptic period-$2$ orbits; (ii) polygonal no-slip billiards always admit small invariant open sets and thus cannot be ergodic with respect to the canonical invariant billiard measure; (iii) the no-slip version of a Sinai billiard must contain linearly stable periodic orbits of period $2$ and, more generally, we provide a curvature threshold at which a commonly occurring period-$2$ orbit shifts from being hyperbolic to being elliptic; (iv) finally, we make a number of observations concerning periodic orbits in a class of polygonal billiards.

preprint2015arXiv

Differential Geometry of Rigid Bodies Collisions and Non-standard Billiards

The configuration manifold $M$ of a mechanical system consisting of two unconstrained rigid bodies in $\mathbb{R}^n$, $n\geq 1$, is a manifold with boundary (typically with singularities.) A complete description of the system requires boundary conditions that specify how orbits should be continued after collisions. A boundary condition is the assignment of a collision map at each tangent space on the boundary of $M$ that gives the post-collision state of the system as a function of the pre-collision state. Our main result is a complete description of the space of linear collision maps satisfying energy and (linear and angular) momentum conservation, time reversibility, and the natural requirement that impulse forces only act at the point of contact of the colliding bodies. These assumptions can be stated in geometric language by making explicit a family of vector subbundles of the tangent bundle to the boundary of $M$: the diagonal, non-slipping, and impulse subbundles. Collision maps at a boundary configuration are shown to be the isometric involutions that restrict to the identity on the non-slipping subspace. The space of such maps is naturally identified with the union of Grassmannians of $k$-dimensional subspaces of $\mathbb{R}^{n-1}$, $0\leq k\leq n-1$, each subspace specifying the directions of contact roughness. We then consider non-standard billiard systems, defined by fixing the position of one of the bodies and allowing boundary conditions different from specular reflection. We also make a few observations of a dynamical nature for simple examples of non-standard billiards and provide a sufficient condition for the billiard map on the space of boundary states to preserve the canonical (Liouville) measure on constant energy hypersurfaces.

preprint2015arXiv

Explicit formulas for reaction probability in reaction-diffusion experiments

A computational procedure is developed for determining the conversion probability for reaction-diffusion systems in which a first-order catalytic reaction is performed over active particles. We apply this general method to systems on metric graphs, which may be viewed as 1-dimensional approximations of 3-dimensional systems, and obtain explicit formulas for conversion. We then study numerically a class of 3-dimensional systems and test how accurately they are described by model formulas obtained for metric graphs. The optimal arrangement of active particles in a 1-dimensional multiparticle system is found, which is shown to depend on the level of catalytic activity: conversion is maximized for low catalytic activity when all particles are bunched together close to the point of gas injection, and for high catalytic activity when the particles are evenly spaced.

preprint2015arXiv

Reaction-diffusion on metric graphs and conversion probability

Motivated by a problem in heterogeneous catalysis, we study a model for irreversible first-order reactions in which gas transport occurs only by diffusion, and reaction occurs only at a small number of well-localized sites. The main problem is to determine fractional conversion in terms of the diffusion coefficient and geometric properties of the reactor. We formulate an appropriate stochastic model for this problem, and then show that, when the domain is composed of a network of thin tubes, reasonably explicit formulas can be obtained.

preprint2012arXiv

Multiple scattering in random mechanical systems and diffusion approximation

This paper is concerned with stochastic processes that model multiple (or iterated) scattering in classical mechanical systems of billiard type, defined below. From a given (deterministic) system of billiard type, a random process with transition probabilities operator P is introduced by assuming that some of the dynamical variables are random with prescribed probability distributions. Of particular interest are systems with weak scattering, which are associated to parametric families of operators P_h, depending on a geometric or mechanical parameter h, that approaches the identity as h goes to 0. It is shown that (P_h -I)/h converges for small h to a second order elliptic differential operator L on compactly supported functions and that the Markov chain process associated to P_h converges to a diffusion with infinitesimal generator L. Both P_h and L are selfadjoint (densely) defined on the space L2(H,η) of square-integrable functions over the (lower) half-space H in R^m, where η is a stationary measure. This measure's density is either (post-collision) Maxwell-Boltzmann distribution or Knudsen cosine law, and the random processes with infinitesimal generator L respectively correspond to what we call MB diffusion and (generalized) Legendre diffusion. Concrete examples of simple mechanical systems are given and illustrated by numerically simulating the random processes.

preprint2012arXiv

Random billiards with wall temperature and associated Markov chains

By a random billiard we mean a billiard system in which the standard specular reflection rule is replaced with a Markov transition probabilities operator P that, at each collision of the billiard particle with the boundary of the billiard domain, gives the probability distribution of the post-collision velocity for a given pre-collision velocity. A random billiard with microstructure (RBM) is a random billiard for which P is derived from a choice of geometric/mechanical structure on the boundary of the billiard domain. RBMs provide simple and explicit mechanical models of particle-surface interaction that can incorporate thermal effects and permit a detailed study of thermostatic action from the perspective of the standard theory of Markov chains on general state spaces. We focus on the operator P itself and how it relates to the mechanical/geometric features of the microstructure, such as mass ratios, curvatures, and potentials. The main results are as follows: (1) we characterize the stationary probabilities (equilibrium states) of P and show how standard equilibrium distributions studied in classical statistical mechanics, such as the Maxwell-Boltzmann distribution and the Knudsen cosine law, arise naturally as generalized invariant billiard measures; (2) we obtain some basic functional theoretic properties of P. Under very general conditions, we show that P is a self-adjoint operator of norm 1 on an appropriate Hilbert space. In a simple but illustrative example, we show that P is a compact (Hilbert-Schmidt) operator. This leads to the issue of relating the spectrum of eigenvalues of P to the features of the microstructure;(3) we explore the latter issue both analytically and numerically in a few representative examples;(4) we present a general algorithm for simulating these Markov chains based on a geometric description of the invariant volumes of classical statistical mechanics.

preprint2009arXiv

Higher order approximation of isochrons

Phase reduction is a commonly used techinque for analyzing stable oscillators, particularly in studies concerning synchronization and phase lock of a network of oscillators. In a widely used numerical approach for obtaining phase reduction of a single oscillator, one needs to obtain the gradient of the phase function, which essentially provides a linear approximation of isochrons. In this paper, we extend the method for obtaining partial derivatives of the phase function to arbitrary order, providing higher order approximations of isochrons. In particular, our method in order 2 can be applied to the study of dynamics of a stable oscillator subjected to stochastic perturbations, a topic that will be discussed in a future paper. We use the Stuart-Landau oscillator to illustrate the method in order 2.