Rheo-NMR in food science—Recent opportunities

For over 25 years, nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) techniques have been used to study materials under mechanical deformation. Collectively, these methods are referred to as Rheo-NMR. In many cases, it provides spatially and temporally resolved maps of NMR spectra, intrinsic NMR parameters (such as relaxation times), or motion (such as diffusion or flow). Therefore, Rheo-NMR is complementary to conventional rheological measurements. This review will briefly summarize current capabilities and limitations of Rheo-NMR in the context of material science and food science in particular. It will report on recent advances such as the incorporation of torque sensors or the implementation of large amplitude oscillatory shear and point out future opportunities for Rheo-NMR in food science.     

Reconciling different formulations of viscous water waves and their mass conservation

The viscosity of water induces a vorticity near the free surface boundary. The resulting rotational component of the fluid velocity vector greatly complicates the water wave system. Several approaches to close this system have been proposed. Our analysis compares three common sets of model equations. The first set has a rotational kinematic boundary condition at the surface. In the second set, a gauge choice for the velocity vector is made that cancels the rotational contribution in the kinematic boundary condition, at the cost of rotational velocity in the bulk and a rotational pressure. The third set circumvents the problem by introducing two domains: the irrotational bulk and the vortical boundary layer. This comparison puts forward the link between rotational pressure on the surface and vorticity in the boundary layer, addresses the existence of nonlinear vorticity terms, and shows where approximations have been used in the models. Furthermore, we examine the conservation of mass for the three systems, and how this can be compared to the irrotational case.     

Steady shear banding in complex fluids

Shear banding (SB) is manifested by the abrupt “demixing” of the flow into regions of high and low shear rate. In this paper, we first relate analytically the rheological parameters of the fluid with the range of shear rates and stresses of SB occurrence. For this, we accept that the origin of shear banding is constitutive, and adopt a non-linear viscoelastic expression able to accommodate the double-valuedness of the stress with flow intensity, under certain conditions. We then implement the model for the case of pressure driven flow through a cylindrical pipe; we derive approximate expressions for the velocity profile in the two-banded regions (core and outer annular), the overall throughput in the presence or absence of “spurt”, and the radial location limits of the shear rate discontinuity.     

Espaces d'écoulements dits « universels », suite 2Spaces of universal flows,part 2

At small velocities, every universal motion can be performed in a fluid whenever the constitutive equation approaches a most general linear model. Both the celebrate theorems of Appell, Cauchy, Helmholtz, Kelvin hold true when dealing with the universal motions. Here the theorem of Bernoulli admits four generalisations. At last, the Helmholtz–Rayleigh theorem about dissipation may be regarded as defining the universal motions. To cite this article: M. Bouthier, C. R. Mecanique 332 (2004).     

Dynamics of vorticity stretching and breakup of isolated viscoelastic droplets in an immiscible viscoelastic matrix

The effect of droplet elasticity on transient deformation of isolated droplets in immiscible polymer blends of equal viscosity was investigated. In terms of the deformation parameter, Def*=a*–c/a*+c where a* and c are apparent drop principal axes, it undergoes two cycles of positive oscillations before reaching a negative value, followed by one cycle of oscillation before attaining a steady state negative value. This behavior was observed when Capillary number, Ca, was varied between 3 and 9 at a fixed Weissenber number, Wi, of 0.31, and when Ca number was held fixed at 8 and Wi number was varied between 0.21 and 0.40. In another blend of relatively lower Wi number of 0.21, one cycle of oscillation in Def* was observed before reaching steady state negative values when Ca number was varied between 3 and 14. The steady state Def* varies inversely with Ca number, with a stronger dependence for the blend with higher Wi number. The magnitude of oscillation increases with increasing Ca and Wi numbers. The critical Ca was found to be 12 and 14 for the two blends studied; these values are about 30 times greater than that of Newtonian blends.     

Estimation of shock induced vorticity on irregular gaseous interfaces: a wavelet-based approach

We study the interaction of a shock with a density-stratified gaseous interface (Richtmyer–Meshkov instability) with localized jagged and irregular perturbations, with the aim of developing an analytical model of the vorticity deposition on the interface immediately after the passage of the shock. The jagged perturbations, meant to simulate machining errors on the surface of a laser fusion target, are characterized using Haar wavelets. Numerical solutions of the Euler equations show that the vortex sheet deposited on the jagged interface rolls into multiple mushroom-shaped dipolar structures which begin to merge before the interface evolves into a bubble-spike structure. The peaks in the distribution of x-integrated vorticity (vorticity integrated in the direction of the shock motion) decay in time as their bases widen, corresponding to the growth and merger of the mushrooms. However, these peaks were not seen to move significantly along the interface at early times i.e. t < 10 τ, where τ is the interface traversal time of the shock. We tested our analytical model against inviscid simulations for two test cases – a Mach 1.5 shock interacting with an interface with a density ratio of 3 and a Mach 10 shock interacting with a density ratio of 10. We find that this model captures the early time (t/τ ∼ 1) vorticity deposition (as characterized by the first and second moments of vorticity distributions) to within 5% of the numerical results. PACS 47.40.Nm; 47.20.Ma     

Some uses of Green's theorem in solving the Navier–Stokes equations

This paper gives a review of methods where Green's theorem may be employed in solving numerically the Navier–Stokes equations for incompressible fluid motion. They are based on the concept of using the theorem to transform local boundary conditions given on the boundary of a closed region in the solution domain into global, or integral, conditions taken over it. Two formulations of the Navier–Stokes equations are considered: that in terms of the streamfunction and vorticity for two-dimensional motion and that in terms of the primitive variables of the velocity components and the pressure. In the first formulation overspecification of conditions for the streamfunction is utilized to obtain conditions of integral type for the vorticity and in the second formulation integral conditions for the pressure are found. Some illustrations of the principle of the method are given in one space dimension, including some derived from two-dimensional flows using the series truncation method. In particular, an illustration is given of the calculation of surface vorticity for two-dimensional flow normal to a flat plate. An account is also given of the implementation of these methods for general two-dimensional flows in both of the mentioned formulations and a numerical illustration is given.     

A note on spin,vorticity and the deformation-rate tensor

The validity of statements in the literature concerning the angular velocity of the principal axes of the stretching tensor is examined and an alternative expression is provided for it.     

A remark on vorticity maximal function

In this note the following inequality is proved. For any nonnegative measure μH⁻¹(R²), xR² and 0<r<1, there holds

(1)

where C is a positive constant. Using (1) an estimate for the vorticity maximal function similar to the estimate in Majda [A. Majda, Remarks on weak solutions for vortex sheets with a distinguished sign, Indiana Univ. Math. J. 42 (1993) 921–939] is established without assuming the initial vorticity having compact support. From this a more simple proof of the Delort's celebrated theorem [J.M. Delort, Existence de mappes de fourbillon en dimension deux, J. Amer. Math. Soc. 4 (1991) 553–586] is presented.     

A spectral solver for the Navier–Stokes equations in the velocity–vorticity formulation

A numerical algorithm intended for the study of flows in a cylindrical container under laminar flow conditions is proposed. High resolution of the flow field, governed by the Navier–Stokes equations in velocity–vorticity formulation relative to a cylindrical frame of reference, is achieved through spatial discretisation by means of the spectral method. This method is based on a Fourier expansion in the azimuthal direction and an expansion in Chebyshev polynomials in the (nonperiodic) radial and axial directions. Several regularity constraints are used to take care of the coordinate singularity. These constraints are implemented, together with the boundary conditions at the top, bottom and mantle of the cylinder, via the tau method. The a priori unknown boundary values of the vorticity are evaluated by means of the influence-matrix technique. The compatibility between the mathematical and numerical formulation of the Navier–Stokes equations is established through a tau-correction procedure. The resolved flow field exhibits high-precision satisfaction of the incompressibility constraints for velocity and vorticity and the definition of the vorticity. The performance of the solver is illustrated by resolution of several configurations representative of generic three-dimensional laminar flows.