Elongational behaviour of electrostatically stabilised and concentrated CNF and CNC hydrogels: Experiments and modelling.

TEMPO-oxidized cellulose nanofibril (CNF) hydrogels or cellulose nanocrystal (CNC) hydrogels can now be obtained at high concentrations (>10 wt%) and used to fabricate biobased materials and structures. Thus, it is required to control and model their rheology in process-induced multiaxial flow conditions using 3D tensorial models. For that purpose, it is necessary to investigate their elongational rheology. Thus, concentrated TEMPO-oxidized CNF and CNC hydrogels were subjected to monotonic and cyclic lubricated compression tests. These tests revealed for the first time that the complex compression rheology of these two electrostatically stabilised hydrogels combines viscoelasticity and viscoplasticity. The effect of their nanofibre content and aspect ratio on their compression response was clearly emphasised and discussed. The ability of a non-linear elasto-viscoplastic model to reproduce the experiments was assessed. Even if some discrepancies were observed at low or high strain rates, the model was consistent with the experiments.

Rheology of concentrated and highly concentrated enzymatic cellulose nanofibril hydrogels during lubricated compression.

Processing cellulose nanofibril (CNF) hydrogels with a high concentration is a solution to reduce logistics costs and drying energy and to produce CNF-based materials with good dimensional stability. However, the rheology of concentrated and highly concentrated CNF hydrogels is poorly understood due to the difficulties to characterise them using standard shear rheometers. In this study, enzymatic CNF hydrogels in the concentrated and highly concentrated regimes (3-13.6 wt%) were subjected to lubricated compression at various strain rates. At low strains, compression curves exhibited a linear regime. At higher strains and low strain rates, a heterogeneous and marked hardening of stress levels was observed and accompanied with a two-phase flow with significant fluid segregation and network consolidation. At high strain rates, a homogeneous and incompressible one-phase plateau-like regime progressively established. In this regime, a yield stress was measured and compared with literature data, showing a good agreement with them.

A fibrous cellulose paste formulation to manufacture structural parts using 3D printing by extrusion.

An optimized paste based on short natural cellulose fibers combined with carboxymethyl cellulose at a high dry content (42 wt.%) was implemented as a bio-based material for 3D printing by extrusion. The homogeneous paste exhibited a pronounced thinning behavior and yield stress; it was extruded using a screw extrusion-based direct ink writing system and could easily flow through a small nozzle. The optimized formulation enabled accurate additive manufacturing of parts using a natural air-drying process with or without an ethanol bath. We characterized the anisotropic shrinkage that occurred during the drying of 3D printed parts and proposed a compensation method to account for it. The obtained results emphasized that cellulose had a strong potential to be used as a raw material for 3D printing of cheap, lightweight, robust, and recyclable parts.

Micro-mechanics of electrostatically stabilized suspensions of cellulose nanofibrils under steady state shear flow.

In this study, we characterized and modeled the rheology of TEMPO-oxidized cellulose nanofibril (NFC) aqueous suspensions with electrostatically stabilized and unflocculated nanofibrous structures. These colloidal suspensions of slender and wavy nanofibers exhibited a yield stress and a shear thinning behavior at low and high shear rates, respectively. Both the shear yield stress and the consistency of these suspensions were power-law functions of the NFC volume fraction. We developed an original multiscale model for the prediction of the rheology of these suspensions. At the nanoscale, the suspensions were described as concentrated systems where NFCs interacted with the Newtonian suspending fluid through Brownian motion and long range fluid-NFC hydrodynamic interactions, as well as with each other through short range hydrodynamic and repulsive colloidal interaction forces. These forces were estimated using both the experimental results and 3D networks of NFCs that were numerically generated to mimic the nanostructures of NFC suspensions under shear flow. They were in good agreement with theoretical and measured forces for model colloidal systems. The model showed the primary role played by short range hydrodynamic and colloidal interactions on the rheology of NFC suspensions. At low shear rates, the origin of the yield stress of NFC suspensions was attributed to the combined contribution of repulsive colloidal interactions and the topology of the entangled NFC networks in the suspensions. At high shear rates, both concurrent colloidal and short (in some cases long) range hydrodynamic interactions could be at the origin of the shear thinning behavior of NFC suspensions.

Heterogeneous flow kinematics of cellulose nanofibril suspensions under shear.

The rheology of NFC suspensions that exhibited different microstructures and colloidal stability, namely TEMPO and enzymatic NFC suspensions, was investigated at the macro and mesoscales using a transparent Couette rheometer combined with optical observations and ultrasonic speckle velocimetry (USV). Both NFC suspensions showed a complex rheology, which was typical of yield stress, non-linear and thixotropic fluids. Hysteresis loops and erratic evolutions of the macroscale shear stress were also observed, thereby suggesting important mesostructural changes and/or inhomogeneous flow conditions. The in situ optical observations revealed drastic mesostructural changes for the enzymatic NFC suspensions, whereas the TEMPO NFC suspensions did not exhibit mesoscale heterogeneities. However, for both suspensions, USV measurements showed that the flow was heterogeneous and exhibited complex situations with the coexistence of multiple flow bands, wall slippage and possibly multidimensional effects. Using USV measurements, we also showed that the fluidization of these suspensions could presumably be attributed to a progressive and spatially heterogeneous transition from a solid-like to a liquid-like behavior. As the shear rate was increased, the multiple coexisting shear bands progressively enlarged and nearly completely spanned over the rheometer gap, whereas the plug-like flow bands were eroded.

In-plane mechanics of soft architectured fibre-reinforced silicone rubber membranes.

Silicone rubber membranes reinforced with architectured fibre networks were processed with a dedicated apparatus, allowing a control of the fibre content and orientation. The membranes were subjected to tensile loadings combined with continuous and discrete kinematical field measurements (DIC and particle tracking). These tests show that the mechanical behaviour of the membranes is hyperelastic at the first order. They highlight the influence of the fibre content and orientation on both the membrane in-plane deformation and stress levels. They also prove that for the considered fibrous architectures and mechanical loadings, the motion and deformation of fibres is an affine function of the macroscale transformation. These trends are fairly well described by the micromechanical model proposed recently in Bailly et al. (JMBBM, 2012). This result proves that these materials are very good candidates for new biomimetic membranes, e.g. to improve aortic analogues used for in vitro experiments, or existing textiles used for vascular (endo)prostheses.

Towards a biomimetism of abdominal healthy and aneurysmal arterial tissues.

The aim of this work is to develop a new hyperelastic and anisotropic material mimicking histological and mechanical features of healthy and aneurysmal arterial tissues. The material is constituted by rhombic periodic lattices of hyperelastic fibres embedded into a soft elastomer membrane. To fit bi-axial experimental data obtained from the literature, with normal or pathologic human abdominal aortic tissues, the microstructure of the periodic lattices (fibre length, angle between fibres) together with the mechanical behaviour of the fibres (fibre tension-elongation curve) were optimised by using theoretical results arising from a multi-scale homogenisation process. It is shown that (i) a material constituted by only one periodic lattice of fibres is clearly not sufficient to describe all the experimental data set, (ii) a quantitative agreement between measurements and theoretical predictions is obtained by using a material with two fibre lattices, (iii) the optimised microstructures and mechanical properties of the fibrous lattices are strongly different for the abdominal healthy and aneurysmal arterial tissues, (iv) the anisotropic mechanical behaviour of the optimised material is described by only five parameters and (v) the optimal angles between fibres in the case of the healthy aorta are consistent with histological data. Several technical solutions of fibres can be considered as relevant candidates: this is illustrated in the particular cases of straight and wavy fibres.

Upscaling the diffusion equations in particulate media made of highly conductive particles. I. Theoretical aspects.

Many analytical and numerical works have been devoted to the prediction of macroscopic effective transport properties in particulate media. Usually, structure and properties of macroscopic balance and constitutive equations are stated a priori. In this paper, the upscaling of the transient diffusion equations in concentrated particulate media with possible particle-particle interfacial barriers, highly conductive particles, poorly conductive matrix, and temperature-dependent physical properties is revisited using the homogenization method based on multiple scale asymptotic expansions. This method uses no a priori assumptions on the physics at the macroscale. For the considered physics and microstructures and depending on the order of magnitude of dimensionless Biot and Fourier numbers, it is shown that some situations cannot be homogenized. For other situations, three different macroscopic models are identified, depending on the quality of particle-particle contacts. They are one-phase media, following the standard heat equation and Fourier's law. Calculations of the effective conductivity tensor and heat capacity are proved to be uncoupled. Linear and steady state continuous localization problems must be solved on representative elementary volumes to compute the effective conductivity tensors for the two first models. For the third model, i.e., for highly resistive contacts, the localization problem becomes simpler and discrete whatever the shape of particles. In paper II [Vassal, Phys. Rev. E 77, 011303 (2008)], diffusion through networks of slender, wavy, entangled, and oriented fibers is considered. Discrete localization problems can then be obtained for all models, as well as semianalytical or fully analytical expressions of the corresponding effective conductivity tensors.

Upscaling the diffusion equations in particulate media made of highly conductive particles. II. Application to fibrous materials.

In paper I [Vassal, Phys. Rev. E77, 011302 (2008)] of this contribution, the effective diffusion properties of particulate media with highly conductive particles and particle-particle interfacial barriers have been investigated with the homogenization method with multiple scale asymptotic expansions. Three different macroscopic models have been proposed depending on the quality of contacts between particles. However, depending on the nature and the geometry of particles contained in representative elementary volumes of the considered media, localization problems to be solved to compute the effective conductivity of the two first models can rapidly become cumbersome, time and memory consuming. In this second paper, the above problem is simplified and applied to networks made of slender, wavy and entangled fibers. For these types of media, discrete formulations of localization problems for all macroscopic models can be obtained leading to very efficient numerical calculations. Semianalytical expressions of the effective conductivity tensors are also proposed under simplifying assumptions. The case of straight monodisperse and homogeneously distributed slender fibers with a circular cross section is further explored. Compact semianalytical and analytical estimations are obtained when fiber-fiber contacts are perfect or very poor. Moreover, two discrete element codes have been developed and used to solve localization problems on representative elementary volumes for the same types of contacts. Numerical results underline the significant roles of the fiber content, the orientation of fibers as well as the relative position and orientation of contacting fibers on the effective conductivity tensors. Semianalytical and analytical predictions are discussed and compared with numerical results.