Time average shear rate, equation

This additional Eq. (18) was discretized at the same resolution as the flow equations, typical grids comprising 1203 and 1803 nodes. At every time step, the local particle concentration is transported within the resolved flow field. Furthermore, the local flow conditions yield an effective 3-D shear rate that can be used for estimating the local agglomeration rate constant /10. Fig. 10 (from Hollander et al., 2003) presents both instantaneous and time-averaged spatial distributions of /i0 in vessels agitated by two different impellers color versions of these plots can be found in Hollander (2002) and in Hollander et al. (2003). 

Clearly flow aligning behavior of the director is present and do increases linearly with the tilt angle, do. Above a threshold in the Spain rate, y 0.011, undulations in vorticity direction set in. In Fig. 14 the results of simulations for y 0.015 are shown. In Fig. 15 we have plotted the undulation amplitude obtained as a function of the shear rate. The dashed line indicates a square root behavior corresponding to a forward bifurcation near the onset of undulations. This is, indeed, what is expected, when a weakly nonlinear analysis based on the underlying macroscopic equations is performed [54], In Fig. 16 we have plotted an example for the dynamic behavior obtained from molecular dynamics simulations. It shows the time evolution after a step-type start for two shear rates below the onset of undulations. The two solid lines correspond to a fit to the data using the solutions of the averaged linearized form of (27). The shear approaches its stationary value for small tilt angle (implied by the use of the linearized equation) with a characteristic time scale t = fi/Bi. 

The flow of viscoplastic fluids through beds of particles has not been studied as extensively as that of power-law fluids. However, since the expressions for the average shear stress and the nominal shear rate at the wall, equations (5.41) and (5.42), are independent of fluid model, they may be used in conjimction with any time-independent behaviour fluid model, as illuslrated here for the streamline flow of Bingham plastic fluids. The mean velocity for a Bingham plastic fluid in a circular tube is given by equation (3.13) ... 

A closer look at Fig.3 shows that the relations (6) and (8) are violated in certain ranges of shear rates. This is expected in real physical systems and for data inferred from NEMD simulations as presented in [1, 6], as well as for the model considered here since the use of a constant friction coefficient in the equations of motion discussed above is an approximation. Physical reality is closer to a time dependent effective friction coefficient as used here. Assuming, as above, that the long time average of the time change of the angular momentum vanishes for a stationary state and that there is no torque associated with the force P, (3) has to be replaced by... 

Experiments are not only made for imposed shear rates but also for imposed shear stress. Calculations intended to give long-time averages for the latter case are performed by replacing the constant shear rate 7 in (21) by the dynamic shear rate g = g t) which obeys the equation... 

The derivation of equation (4.46) is based on the shear force on an element cube of liquid, causing torsional work which is evaluated as work/time = power. Consequently, the value of G assumes that all elemental cubes of liquid in the volume V are being sheared at the same rate (on average) and that G is also a time-averaged value. ... 

The program uses the input value of concentration and relative flow times to calculate a molecular weight according to Equations 10.17 and 10.9. This value is used to calculate a zero shear viscosity, and average values for the other parameters of the HN model discussed in the text are adopted to calculate the frequency dependence (i.e., to calculate the viscoelastic moduli as a function of shear rate). All the values adopted in the model are disclosed on the worksheet page. 

Empirical rate equations have been reported for the mastication degradation of many polymers, for example, polyethylene [29], poly(methyl methacrylate) [30], styrene rubbers [31, 32], polychloroprene [33], and EPDM [34]. Reaction rate generally depends on [35]. Experiments on degradation rate have been performed also on the molten state [15, 16, 35]. Pohl and co-workers [15, 16] reported their data as the fraction of bonds broken (Eq. 2.17) as a function of capillary residence time. They found that the average rate of bond breaking b depends on shear rate y and absolute temperature T by the following equation for polyisobutylene ... 

As in equilibrium molecular dynamics, the equations of motion have to be solved for a system with periodic boundaries. For shear, the boundaries are modified to become the Lees-Edwards sliding brick conditions (Lees Edwards 1972), in which periodic images of the simulation cell above and below the unit cell are moved in opposite directions at a velocity determined by the imposed shear rate (see Fig. 9.9). The properties of the system follow firom the appropriate time averages, <. . >, usually (but not necessarily) after the system has reached the steady state. Given, for example, a system at a number density, n = N/V, under an applied shear rate, the kinetic temperature is constrained with an appropriate thermostat Different properties can then be evaluated, for example, the internal energy. 

The discussion on the distorted fluid microstructure demonstrates the utility of NEMD as a practical tool and as a means to get an insight into the physics of viscous phenomena. The shear-rate dependent non-Newtonian properties of a liquid can be calculated from the appropriate steady-state averages, but the mechanical transport properties can be obtained indirectly from the nonequilibrium pair distribution function using equation (9.49), and a more physical picture of the behavior of a nonequilibrium fluid is obtained. As an example, the phenomenological theory of Hanley et al. (1987,1988) is outlined where the functionals are evaluated from a relaxation time approximation. The... 

To close this Section we comment on two papers that do not fit under any neat heading. The first of these is by Xiao et al,261 who study the final stages of the collapse of an unstable bubble or cavity using MD simulations of an equilibrated Lennard-Jones fluid from which a sphere of molecules has been removed. They find that the temperature inside this bubble can reach up to an equivalent of 6000 K for water. It is at these temperatures that sonolumines-cence is observed experimentally. The mechanism of bubble collapse is found to be oscillatory in time, in agreement with classical hydrodynamics predictions and experimental observation. The second paper, by Lue,262 studies the collision statistics of hard hypersphere fluids by MD in 3, 4 and 5 dimensions. Equations of state, self-diffusion coefficients, shear viscosities and thermal conductivities are determined as functions of density. Exact expressions for the mean-free path in terms of the average collision time and the compressibility factor in terms of collision rate are also derived. Work such as this, abstract as it may appear, may be valuable in the development of microscopic theories of fluid transport as well as provide insight into transport processes in general. 

---