The maximum spreading of drops impacting on smooth and rough surfaces is measured from low to high impact velocity for liquids with different surface tensions and viscosities. We demonstrate that dynamic wetting plays an important role in the spreading at low velocity, characterized by the dynamic contact angle at maximum spreading. In the energy balance, we account for the dynamic wettability by introducing the capillary energy at zero impact velocity, which relates to the spreading ratio at zero impact velocity. Correcting the measured spreading ratio by the spreading ratio at zero velocity, we find a correct scaling behaviour for low and high impact velocity and, by interpolation between the two, we find a universal scaling curve. The influence of the liquid as well as the nature and roughness of the surface are taken into account properly by rescaling with the spreading ratio at zero velocity, which, as demonstrated, is equivalent to accounting for the dynamic contact angle.

Droplet impact has been imaged on different rigid, smooth, and rough substrates for three liquids with different viscosity and surface tension, with special attention to the lower impact velocity range. Of all studied parameters, only surface tension and viscosity, thus the liquid properties, clearly play a role in terms of the attained maximum spreading ratio of the impacting droplet. Surface roughness and type of surface (steel, aluminum, and parafilm) slightly affect the dynamic wettability and maximum spreading at low impact velocity. The dynamic contact angle at maximum spreading has been identified to properly characterize this dynamic spreading process, especially at low impact velocity where dynamic wetting plays an important role. The dynamic contact angle is found to be generally higher than the equilibrium contact angle, showing that statically wetting surfaces can become less wetting or even nonwetting under dynamic droplet impact. An improved energy balance model for maximum spreading ratio is proposed based on a correct analytical modeling of the time at maximum spreading, which determines the viscous dissipation. Experiments show that the time at maximum spreading decreases with impact velocity depending on the surface tension of the liquid, and a scaling with maximum spreading diameter and surface tension is proposed. A second improvement is based on the use of the dynamic contact angle at maximum spreading, instead of quasi-static contact angles, to describe the dynamic wetting process at low impact velocity. This improved model showed good agreement compared to experiments for the maximum spreading ratio versus impact velocity for different liquids, and a better prediction compared to other models in literature. In particular, scaling according to We(1/2) is found invalid for low velocities, since the curves bend over to higher maximum spreading ratios due to the dynamic wetting process.

Hysteresis is observed in sorption-induced swelling in various soft nanoporous polymers. The associated coupling mechanism responsible for the observed sorption-induced swelling and associated hysteresis needs to be unraveled. Here we report a microscopic scenario for the molecular mechanism responsible for hysteresis in sorption-induced swelling in natural polymers such as cellulose using atom-scale simulation; moisture content and swelling exhibit hysteresis upon ad- and desorption but not swelling versus moisture content. Different hydrogen bond networks are examined; cellulose swells to form water–cellulose bonds upon adsorption but these bonds do not break upon desorption at the same chemical potential. These findings, which are supported by mechanical testing and cellulose textural assessment upon sorption, shed light on experimental observations for wood and other related materials.

No abstract

The quintessential form of cellulose in wood consists of microfibrils that have high aspect ratio crystalline domains embedded within an amorphous cellulose domain. In this study, we apply unitedatom molecular dynamics simulations to quantify changes in different morphologies of cellulose. We compare the structure of crystalline cellulose with paracrystalline and amorphous phases that are both obtained by high temperature equilibration followed by quenching at room temperature. Our study reveals that the paracrystalline phase may be an intermediate, kinetically arrested phase formed upon amorphisation of crystalline cellulose. The quenched structures yield isotropic amorphous polymer domains consistent with experimental results, thereby validating a new computational protocol for achieving amorphous cellulose structure. The non-crystalline cellulose compared to crystalline structure is characterized by a dramatic decrease in elastic modulus, thermal expansion coefficient, bond energies, and number of hydrogen bonds. Analysis of the lattice parameters shows that Ib cellulose undergoes a phase transition into hightemperature phase in the range of 450-550 K. The mechanisms of the phase transition elucidated here present an atomistic view of the temperature dependent dynamic structure and mechanical properties of cellulose. The paracrystalline state of cellulose exhibits intermediate mechanical properties, between crystalline and amorphous phases, that can be assigned to the physical properties of the interphase regions between crystalline and amorphous cellulose in wood microfibrils. Our results suggest an atomistic structural view of amorphous cellulose which is consistent with Electronic supplementary material The online version of this article (

The interaction of water with many biopolymers is known to rearrange their internal structure, make them moisture sensitive, and influence their physical properties. We study amorphous cellulose and hemicellulose, two hydrophilic biopolymers, using molecular dynamics simulations, and we analyze their structural and physical properties over the full range of moisture content. We find a quasi-linear dependence of volumetric strain on moisture content, and a linear scaling between volumetric strain and porosity, showing that swelling is directly related to the space created by adsorbed water molecules. The interaction of water with the polymer structure results in a weakening of the mechanical properties, leading to rubberlike behavior at high moisture content. Weakening is caused by a decrease in the number of hydrogen bonds that follow exponential scaling. Breaking of the hydrogen bonds system is found to control not only the mechanical response but also the evolution of porosity and the volumetric strain.

No abstract

A two-phase model of a wood microfibril consisting of crystalline cellulose and amorphous hemicellulose is investigated with molecular dynamics in full range of sorption to understand the molecular origin of swelling and weakening of wood. Water is adsorbed in hemicellulose, and an excess of sorption is found at the interface, while no sorption occurs within cellulose. Water molecules adsorbed on the interface push away polymer chains, forcing the two phases to separate and causing breaking of h-bonds, particularly pronounced on the interface. Existence of two different regions in moisture response is demonstrated. At low moisture content, water is uniformly adsorbed within hemicellulose, breaking a small amount of hydrogen bonds. Microfibril does not swell, and the porosity does not change. As moisture content increases, water is adsorbed preferentially at the interface, which leads to additional swelling and porosity increase at the interface. Young's and shear moduli decrease importantly due to breaking of h-bonds and screening of the long-range interactions.