1
|
Moritz C, Sega M, Innerbichler M, Geissler PL, Dellago C. Weak scaling of the contact distance between two fluctuating interfaces with system size. Phys Rev E 2020; 102:062801. [PMID: 33465946 DOI: 10.1103/physreve.102.062801] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2020] [Accepted: 10/30/2020] [Indexed: 06/12/2023]
Abstract
A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ^{*} depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ^{*} scales weakly with system size, i.e., the dependence of ℓ^{*} on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim_{L→∞}(ℓ^{*}/L)=0. In particular, for (2+1)-dimensional interfaces ℓ^{*} is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. This weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.
Collapse
Affiliation(s)
- Clemens Moritz
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Marcello Sega
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Fürther Straße 248, 90429 Nürnberg, Germany
| | - Max Innerbichler
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Phillip L Geissler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
- Erwin Schrödinger Institute for Mathematics and Physics, Boltzmanngasse 9, 1090, Vienna, Austria
| |
Collapse
|
2
|
Abstract
Rough surfaces are in contact locally by the peaks of roughness. At this local scale, the pressure of contact can be sharply superior to the macroscopic pressure. If the roughness is assumed to be a random morphology, a well-established observation in many practical cases, mechanical indicators built from the contact zone are then also random variables. Consequently, the probability density function (PDF) of any mechanical random variable obviously depends upon the morphological structure of the surface. The contact pressure PDF, or the probability of damage of this surface can be determined for example when plastic deformation occurs. In this study, the contact pressure PDF is modeled using a particular probability density function, the generalized Lambda distributions (GLD). The GLD are generic and polymorphic. They approach a large number of known distributions (Weibull, Normal, and Lognormal). The later were successfully used to model damage in materials. A semi-analytical model of elastic contact which takes into account the morphology of real surfaces is used to compute the contact pressure. In a first step, surfaces are simulated by Weierstrass functions which have been previously used to model a wide range of surfaces met in tribology. The Lambda distributions adequacy is qualified to model contact pressure. Using these functions, a statistical analysis allows us to extract the probability density of the maximal pressure. It turns out that this density can be described by a GLD. It is then possible to determine the probability that the contact pressure generates plastic deformation.
Collapse
|
3
|
Sahli R, Pallares G, Papangelo A, Ciavarella M, Ducottet C, Ponthus N, Scheibert J. Shear-Induced Anisotropy in Rough Elastomer Contact. PHYSICAL REVIEW LETTERS 2019; 122:214301. [PMID: 31283347 DOI: 10.1103/physrevlett.122.214301] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/06/2019] [Indexed: 05/22/2023]
Abstract
True contact between randomly rough solids consists of myriad individual microjunctions. While their total area controls the adhesive friction force of the interface, other macroscopic features, including viscoelastic friction, wear, stiffness, and electric resistance, also strongly depend on the size and shape of individual microjunctions. We show that, in rough elastomer contacts, the shape of microjunctions significantly varies as a function of the shear force applied to the interface. This process leads to a growth of anisotropy of the overall contact interface, which saturates in the macroscopic sliding regime. We show that smooth sphere-plane contacts have the same shear-induced anisotropic behavior as individual microjunctions, with a common scaling law over 4 orders of magnitude in the initial area. We discuss the physical origin of the observations in light of a fracture-based adhesive contact mechanics model, described in the companion article, which captures the smooth sphere-plane measurements. Our results shed light on a generic, overlooked source of anisotropy in rough elastic contacts, not taken into account in current rough contact mechanics models.
Collapse
Affiliation(s)
- R Sahli
- Univ Lyon, Ecole Centrale de Lyon, ENISE, ENTPE, CNRS, Laboratoire de Tribologie et Dynamique des Systèmes LTDS, UMR 5513, F-69134 Ecully, France
| | - G Pallares
- Univ Lyon, Ecole Centrale de Lyon, ENISE, ENTPE, CNRS, Laboratoire de Tribologie et Dynamique des Systèmes LTDS, UMR 5513, F-69134 Ecully, France
- CESI, LINEACT, Zone Aéroportuaire Méditerranée, 34130 Mauguio, France
| | - A Papangelo
- Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari, Viale Japigia 182, 70126 Bari, Italy
- Hamburg University of Technology, Department of Mechanical Engineering, Am Schwarzenberg-Campus 1, 21073 Hamburg, Germany
| | - M Ciavarella
- Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari, Viale Japigia 182, 70126 Bari, Italy
- Hamburg University of Technology, Department of Mechanical Engineering, Am Schwarzenberg-Campus 1, 21073 Hamburg, Germany
| | - C Ducottet
- Univ Lyon, UJM-Saint-Etienne, CNRS, IOGS, Laboratoire Hubert Curien UMR5516, F-42023 Saint-Etienne, France
| | - N Ponthus
- Univ Lyon, Ecole Centrale de Lyon, ENISE, ENTPE, CNRS, Laboratoire de Tribologie et Dynamique des Systèmes LTDS, UMR 5513, F-69134 Ecully, France
| | - J Scheibert
- Univ Lyon, Ecole Centrale de Lyon, ENISE, ENTPE, CNRS, Laboratoire de Tribologie et Dynamique des Systèmes LTDS, UMR 5513, F-69134 Ecully, France
| |
Collapse
|