Gaetano D'Avino

Rotation of a sphere in a viscoelastic liquid subjected to shear flow. Part I: Simulation results

In inertialess suspensions of rigid particles, the rotational motion of each particle is governed by the so-called freely rotating condition, whereby the total torque acting on the particle must be zero. In this work, we study the effect of viscoelasticity of the suspending liquid on the rotation period of a sphere by means of three-dimensional finite element simulations, for conditions corresponding to a macroscopic shear flow. The simulation results capture the slowing down of the rotation, relative to the Newtonian case, which was recently observed in experiments. It is shown that such a phenomenon depends on the specific constitutive equation adopted for the viscoelastic liquid. Analysis of transients shows a clear correlation between rotation rate and the development of first normal stress difference.

Rheology of viscoelastic suspensions of spheres under small and large amplitude oscillatory shear by numerical simulations

The dynamic response of a viscoelastic suspension of spheres under small and large amplitude oscillatory shear is investigated by three-dimensional direct numerical simulations. A sliding triperiodic domain is implemented whereby the computational domain is regarded as the bulk of an infinite suspension. A fictitious domain method is used to manage the particle motion. After the stress field is computed, the bulk properties are recovered by an averaging procedure. The numerical method is validated by comparing the computed linear viscoelastic response of Newtonian and non-Newtonian suspensions with previous theories and simulations. The numerical predictions are in very good quantitative agreement with experimental data for the Newtonian case, whereas deviations are found with respect to some sets of experiments for semidilute and concentrated viscoelastic suspensions. To investigate on such discrepancies, the effect of aggregates in the bulk of the suspension is examined. The simulations show that the pr...

Migration of a sphere in a viscoelastic fluid under planar shear flow: Experiments and numerical predictions

Solid–liquid disperse systems are complex fluids, extensively used for a variety of applications. The flow behavior of such materials is strongly dependent on the matrix rheological properties: when the suspending fluid shows viscoelastic behavior, isolated suspended particles in confined flow fields tend to migrate across the streamlines. In this work, we carry out a detailed experimental investigation of such migration phenomenon for a rigid, spherical particle in planar shear flow, and compare the experimental findings with 3D numerical simulations based on finite elements. A highly elastic non-Newtonian fluid is chosen as the suspending medium. The migration is studied by varying the confinement, the external shear rate and the fluid normal stresses. The particle is always found to move towards the closest wall, regardless of its initial position. For a fixed set of geometrical and flow parameters, the trajectories from different initial positions can be shifted in time to overlap on a single curve. Migration dynamics is governed by an exponential law as the particle is in the region around the midgap. Stronger confinements, higher shear rates and higher levels of normal stresses all enhance the migration speed. We found qualitative agreement between predictions from numerical simulations and experimental results.

Analysis of dynamic mechanical response in torsion

We investigate the dynamic response of industrial rubbers (styrene-butadiene random copolymers, SBR) in torsion and compare against common small amplitude oscillatory shear measurements by using a torsion rectangular fixture, a modified torsion cylindrical fixture, and a conventional parallel plate fixture, respectively, in two different rheometers (ARES 2kFRTN1 from TA Instruments, USA and MCR 702 from Anton Paar-Physica, Austria). The effects of specimen geometry (length-to-width aspect ratio) on storage modulus and level of clamping are investigated. For cylindrical specimens undergoing torsional deformation, we find that geometry and clamping barely affect the shear moduli, and the measurements essentially coincide with those using parallel plates. In contrast, a clear dependence of the storage modulus on the aspect ratio is detected for specimens having rectangular cross section. The empirical correction used routinely in this test is based on geometrical factors and can account for clamping effects, but works only for aspect ratios above a threshold value of 1.4. By employing a finite element analysis, we perform a parametric study of the effects of the aspect ratio in the cross-sectional stress distribution and the linear viscoelastic torsional response. We propose a new, improved empirical equation for obtaining accurate moduli values in torsion at different aspect ratios, whose general validity is demonstrated in both rheometers. These results should provide a guideline for measurements with different elastomers, for which comparison with dynamic oscillatory tests may not be possible due to wall slip issues.

Dynamics of prolate spheroidal elastic particles in confined shear flow

We investigate through numerical simulations the dynamics of a neo-Hookean elastic prolate spheroid suspended in a Newtonian fluid under shear flow. Both initial orientations of the particle within and outside the shear plane and both unbounded and confined flow geometries are considered. In unbounded flow, when the particle starts on the shear plane, two stable regimes of motion are found, i.e., trembling, where the particle shape periodically elongates and compresses in the shear plane and the angle between its major semiaxis and the flow direction oscillates around a positive mean value, and tumbling, where the particle shape periodically changes and its major axis performs complete revolutions around the vorticity axis. When the particle is initially oriented out of the shear plane, more complex dynamics arise. Geometric confinement of the particle between the moving walls also influences its deformation and regime of motion. In addition, when the particle is initially located in an asymmetric position with respect to the moving walls, particle lateral migration is detected. The effects on the particle dynamics of the geometric and physical parameters that rule the system are investigated.

On the choice of the optimal periodic operation for a continuous fermentation process

In this contribution we investigate the impact of the forcing waveform on the productivity of a continuous bioreactor governed by an unstructured, nonlinear kinetic model. The (periodic) forcing is applied on the substrate concentration in the feed. To this end, some alternative waveforms commonly encountered in practice are evaluated and their performance is compared. An analytical/numerical approach is used. The preliminary analytical step is based on the π‐criterion that gives useful information for small amplitudes. The extension to larger amplitudes, when significant improvements are expected, is then performed through a continuation‐optimization procedure. It is found that the choice of the specific waveform has an impact on the performance of the process and there is no unique best forcing for any process condition, but its choice depends on the operating parameters and the forcing amplitude and frequency values. Further, the influence of the waveform functions on the wash‐out conditions are extensively examined. The analysis shows that all the waveforms examined in this work may lead to significant enlargement of the nontrivial regime with respect to a steady state operation. In particular, square‐wave forcing leads in practice to the extinction of the wash‐out conditions for any feed substrate concentration and for a well defined choice of the forcing parameters.

Finite element formulation of fluctuating hydrodynamics for fluids filled with rigid particles using boundary fitted meshes

In this paper, we present a finite element implementation of fluctuating hydrodynamics with a moving boundary fitted mesh for treating the suspended particles. The thermal fluctuations are incorporated into the continuum equations using the Landau and Lifshitz approach 1. The proposed implementation fulfills the fluctuation-dissipation theorem exactly at the discrete level. Since we restrict the equations to the creeping flow case, this takes the form of a relation between the diffusion coefficient matrix and friction matrix both at the particle and nodal level of the finite elements. Brownian motion of arbitrarily shaped particles in complex confinements can be considered within the present formulation. A multi-step time integration scheme is developed to correctly capture the drift term required in the stochastic differential equation (SDE) describing the evolution of the positions of the particles.The proposed approach is validated by simulating the Brownian motion of a sphere between two parallel plates and the motion of a spherical particle in a cylindrical cavity. The time integration algorithm and the fluctuating hydrodynamics implementation are then applied to study the diffusion and the equilibrium probability distribution of a confined circle under an external harmonic potential.

Determination of the Optimal Periodic Waveform for a Continuous Fermentation Process

Abstract Forced bioreactors may lead to better performances with respect to stationary schemes as it is widely reported in the literature. In this work, we investigate the impact of the forcing waveform on the process productivity of a continuous bioreactor governed by an unstructured, non-linear kinetic model. The (periodic) forcing is applied on the substrate concentration in the feed. An analytical/numerical approach is used. The analytical step is based on the -criterion that gives useful information for small amplitudes. The extension to larger amplitudes, when significant improvements are expected, is performed through a continuation-optimization procedure. It is found that the choice of the specific waveform has an impact on the performance of the process and there is no unique best forcing for any process condition but its choice depends on the operating parameters. The influence of the waveform functions on the wash-out conditions are extensively examined. The analysis shows that all the waveforms examined in this work lead to a significant enlargement of the non-trivial regime with respect to a steady state operation. In particular, the maximum extension is predicted for a square-wave forcing.