Heat transfer in plane Couette flow of a rarefied gas using Bhatnagar-Gross-Krook model

Using the linearized BGK model and the method of moments of half-range distribution functions the temperature jumps at two plates are determined, and it is found that the results are in fair agreement with those of Gross and Ziering, and Ziering.

Heat transfer in rarefied couette flow on the basis of discrete ordinate method

The method of discrete ordinates, in conjunction with the modified "half-range" quadrature, is applied to the study of heat transfer in rarefied gas flows. Analytic expressions for the reduced distribution function, the macroscopic temperature profile and the heat flux are obtained in the general n-th approximation. The results for temperature profile and heat flux are in sufficiently good accord both with the results of the previous investigators and with the experimental data.

Non-linear couette flow with heat transfer using the B-G-K model

In recent years a large number of investigators have devoted their efforts to the study of flow and heat transfer in rarefied gases, using the BGK [1] model or the Boltzmann kinetic equation. The velocity moment method which is based on an expansion of the distribution function as a series of orthogonal polynomials in velocity space, has been applied to the linearized problem of shear flow and heat transfer by Mott-Smith [2] and Wang Chang and Uhlenbeck [3]. Gross, Jackson and Ziering [4] have improved greatly upon this technique by expressing the distribution function in terms of half-range functions and it is this feature which leads to the rapid convergence of the method. The full-range moments method [4] has been modified by Bhatnagar [5] and then applied to plane Couette flow using the B-G-K model. Bhatnagar and Srivastava [6] have also studied the heat transfer in plane Couette flow using the linearized B-G-K equation. On the other hand, the half-range moments method has been applied by Gross and Ziering [7] to heat transfer between parallel plates using Boltzmann equation for hard sphere molecules and by Ziering [83 to shear and heat flow using Maxwell molecular model. Along different lines, a moment method has been applied by Lees and Liu [9] to heat transfer in Couette flow using Maxwell's transfer equation rather than the Boltzmann equation for distribution function. An iteration method has been developed by Willis [10] to apply it to non-linear heat transfer problems using the B-G-K model, with the zeroth iteration being taken as the solution of the collisionless kinetic equation. Krook [11] has also used the moment method to formulate the equivalent continuum equations and has pointed out that if the effects of molecular collisions are described by the B-G-K model, exact numerical solutions of many rarefied gas-dynamic problems can be obtained. Recently, these numerical solutions have been obtained by Anderson [12] for the non-linear heat transfer in Couette flow,

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear ow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear flow of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Linear stability, transient growth and the role of viscosity stratification in compressible Couette flow

Linear stability and the nonmodal transient energy growth in compressible plane Couette flow are investigated for two prototype mean flows: (a) the uniform shear flow with constant viscosity, and (b) the nonuniform shear flow with stratified viscosity. Both mean flows are linearly unstable for a range of supersonic Mach numbers (M). For a given M, the critical Reynolds number (Re) is significantly smaller for the uniform shear flow than its nonuniform shear counterpart; for a given Re, the dominant instability (over all streamwise wave numbers, α) of each mean flow belongs to different modes for a range of supersonic M. An analysis of perturbation energy reveals that the instability is primarily caused by an excess transfer of energy from mean flow to perturbations. It is shown that the energy transfer from mean flow occurs close to the moving top wall for “mode I” instability, whereas it occurs in the bulk of the flow domain for “mode II.” For the nonmodal transient growth analysis, it is shown that the maximum temporal amplification of perturbation energy, Gmax, and the corresponding time scale are significantly larger for the uniform shear case compared to those for its nonuniform counterpart. For α=0, the linear stability operator can be partitioned into L∼L̅ +Re2 Lp, and the Re-dependent operator Lp is shown to have a negligibly small contribution to perturbation energy which is responsible for the validity of the well-known quadratic-scaling law in uniform shear flow: G(t∕Re)∼Re2. In contrast, the dominance of Lp is responsible for the invalidity of this scaling law in nonuniform shear flow. An inviscid reduced model, based on Ellingsen-Palm-type solution, has been shown to capture all salient features of transient energy growth of full viscous problem. For both modal and nonmodal instability, it is shown that the viscosity stratification of the underlying mean flow would lead to a delayed transition in compressible Couette flow.

Rheometry of granular materials in cylindrical Couette cells: Anomalous stress caused by gravity and shear

The cylindrical Couette device is commonly employed to study the rheology of fluids, but seldom used for dense granular materials. Plasticity theories used for granular flows predict a stress field that is independent of the shear rate, but otherwise similar to that in fluids. In this paper we report detailed measurements of the stress as a function of depth, and show that the stress profile differs fundamentally from that of fluids, from the predictions of plasticity theories, and from intuitive expectation. In the static state, a part of the weight of the material is transferred to the walls by a downward vertical shear stress, bringing about the well-known Janssen saturation of the stress in vertical columns. When the material is sheared, the vertical shear stress changes sign, and the magnitudes of all components of the stress rise rapidly with depth. These qualitative features are preserved over a range of the Couette gap and shear rate, for smooth and rough walls and two model granular materials. To explain the anomalous rheological response, we consider some hypotheses that seem plausibleapriori, but showthat none survive after careful analysis of the experimental observations. We argue that the anomalous stress is due to an anisotropic fabric caused by the combined actions of gravity, shear, and frictional walls, for which we present indirect evidence from our experiments. A general theoretical framework for anisotropic plasticity is then presented. The detailed mechanics of how an anisotropic fabric is brought about by the above-mentioned factors is not clear, and promises to be a challenging problem for future investigations. (C) 2013 AIP Publishing LLC.

Nonlinear viscoelasticity of entangled wormlike micellar fluid under large-amplitude oscillatory shear: Role of elastic Taylor-Couette instability

The role of elastic Taylor-Couette flow instabilities in the dynamic nonlinear viscoelastic response of an entangled wormlike micellar fluid is studied by large-amplitude oscillatory shear (LAOS) rheology and in situ polarized light scattering over a wide range of strain and angular frequency values, both above and below the linear crossover point. Well inside the nonlinear regime, higher harmonic decomposition of the resulting stress signal reveals that the normalized third harmonic I-3/I-1 shows a power-law behavior with strain amplitude. In addition, I-3/I-1 and the elastic component of stress amplitude sigma(E)(0) show a very prominent maximum at the strain value where the number density (n(v)) of the Taylor vortices is maximum. A subsequent increase in applied strain (gamma) results in the distortions of the vortices and a concomitant decrease in n(v), accompanied by a sharp drop in I-3 and sigma(E)(0). The peak position of the spatial correlation function of the scattered intensity along the vorticity direction also captures the crossover. Lissajous plots indicate an intracycle strain hardening for the values of gamma corresponding to the peak of I-3, similar to that observed for hard-sphere glasses.

Classical and Cosserat plasticity and viscoplasticity models for slow granular flow

We consider models for the rheology of dense, slowly deforming granular materials based of classical and Cosserat plasticity, and their viscoplastic extensions that account for small but finite particle inertia. We determine the scale for the viscosity by expanding the stress in a dimensionless parameter that is a measure of the particle inertia. We write the constitutive relations for classical and Cosserat plasticity in stress-explicit form. The viscoplastic extensions are made by adding a rate-dependent viscous stress to the plasticity stress. We apply the models to plane Couette flow, and show that the classical plasticity and viscoplasticity models have features that depart from experimental observations; the prediction of the Cosserat viscoplasticity model is qualitatively similar to that of Cosserat plasticity, but the viscosities modulate the thickness of the shear layer.

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Stability and transition in the flow of polymer solutions

The linear stability analysis of a plane Couette flow of viscoelastic fluid have been studied with the emphasis on two dimensional disturbances with wave number k similar to Re-1/2, where Re is Reynolds number based on maximum velocity and channel width. We employ three models to represent the dilute polymer solution: the classical Oldroyd-B model, the Oldroyd-B model with artificial diffusivity and the non-homogeneous polymer model. The result of the linear stability analysis is found to be sensitive to the polymer model used. While the plane Couette flow is found to be stable to infinitesimal disturbances for the first two models, the last one exhibits a linear instability.

The stress in a slowly sheared granular column

We report measurements of the wall stress in a granular material sheared in a cylindrical Couette cell, as a function of the distance from the free surface. Our results shows that when the material is static, all components of the stress saturate to constant values within a short distance from the free surface, in conformity with earlier experiments and theoretical predictions. When the material is sheared by rotating the inner cylinder at a constant rate, the stresses are remarkably altered. The radial normal stress does not saturate, and increases even more rapidly with depth than the linear hydrostatic pressure profile. The axial shear stress changes sign on shearing, and its magnitude increases with depth. These results are discussed in the context of the predictions of the classical and Cosserat plasticity theories.

Semi-similar solutions of the unsteady compressible second-order boundary layer flow at the stagnation point

Semi-similar solutions of the unsteady compressible laminar boundary layer flow over two-dimensional and axisymmetric bodies at the stagnation point with mass transfer are studied for all the second-order boundary layer effects when the free stream velocity varies arbitrarily with time. The set of partial differential equations governing the unsteady compressible second-order boundary layers representing all the effects are derived for the first time. These partial differential equations are solved numerically using an implicit finite-difference scheme. The results are obtained for two particular unsteady free stream velocity distributions: (a) an accelerating stream and (b) a fluctuating stream. It is observed that the total skin friction and heat transfer are strongly affected by the surface mass transfer and wall temperature. However, their variation with time is significant only for large times. The second-order boundary layer effects are found to be more pronounced in the case of no mass transfer or injection as compared to that for suction. Résumé Des solutions semi-similaires d'écoulement variable compressible de couche limite sur des corps bi-dimensionnels thermique, sont étudiées pour tous les effets de couche limite du second ordre, lorsque la vitesse de l'écoulement libre varie arbitrairement avec le temps. Le systéme d'équations aux dérivées partielles représentant tous les effets est écrit pour la premiére fois. On le résout numériquement á l'aide d'un schéma implicite aux différences finies. Les résultats sont obtenus pour deux cas de vitesse variable d'écoulement libre: (a) un écoulement accéléré et (b) un écoulement fluctuant. On observe que le frottement pariétal total et le transfert de chaleur sont fortement affectés par le transfert de masse et la température pariétaux. Néanmoins, leur variation avec le temps est sensible seulement pour des grandes durées. Les effets sont trouvés plus prononcés dans le cas de l'absence du transfert de masse ou de l'injection par rapport au cas de l'aspiration.

Vapor phase oxidation of p-xylene

The kinetics of the vapor phase oxidation of p-xylene over ferric molybdate catalyst were studied in an isothermal, differential, tubular flow reactor in the temperature range of 360 to 420° C. The major product obtained was p-tolualdehyde with small amounts of maleic anhydride and p-toluic acid. No terephthalic acid or CO2 were observed. The reaction rate data collected fit the redox model given by Equation 1. The values of activation energies Ex, Eo and frequency factors Ax, Ao obtained are 72, 63 kJ/mol and 0.64, 2.89 m3/kg catalyst s respectively. The reaction mechanism was established by studying the oxidation of p-tolualdehyde, toluic and terephthalic acids. It is concluded that the reaction follows a parallel-consecutive scheme. On a étudié la cinétique de l'oxydation, en phase gazeuse, du para-xylène sur un catalyseur consistant en molybdate ferrique; cette oxydation s'est faite dans un réacteur à écoulement tubulaire, isothermique et différentiel, dans une échelle de températures comprises entre 360°C et 420°C. Le produit principal obtenu a été le para-tolualdéhyde; on a aussi trouvé de faibles quantités d'anhydride maléique et d'acide para-toluique, mais on n'a pas noté la présence d'acide téréphtalique ni d'anhydride carbonique (CO2). Les résultats obtenus en ce qui a trait à la vitesse de réaction concordent bien avec les données du modèle redox indiquées par l'équation 1. Les valeurs des énergies d'activation Ex et Eo ainsi que des facteurs de fréquence Ax et Ao obtenus sont respectivement 72 et 63 kilojoules/mol. et 0.64 × 103 et 2.89 m3/kg de catalyseur. On a établi le mécanisme de la réaction en étudiant l'oxydation du para-tolualdéhyde et des acides toluique et téréphtalique. On conclut que la réaction se fait d'une manière parallèle et consécutive.