arxiv:cond-mat/ v1 [cond-mat.soft] 5 May 2006
|
|
- Horace Dalton
- 5 years ago
- Views:
Transcription
1 Shear-induced crstallization of a dense rapid granular flow: hdrodnamics beond the melting point? Evgeni Khain and Baruch Meerson 2 Department of Phsics and Michigan Center for Theoretical Phsics, The Universit of Michigan, Ann Arbor, Michigan 489 and 2 Racah Institute of Phsics, Hebrew Universit of Jerusalem, Jerusalem 994, Israel arxiv:cond-mat/6534v cond-mat.soft] 5 Ma 26 We investigate shear-induced crstallization in a ver dense flow of mono-disperse inelastic hard spheres. We consider a stead plane Couette flow under constant pressure and neglect gravit. We assume that the granular densit is greater than the melting point of the equilibrium phase diagram of elastic hard spheres. We emplo a Navier-Stokes hdrodnamics with constitutive relations all of which except the shear viscosit) diverge at the crstal packing densit, while the shear viscosit diverges at a smaller densit. The phase diagram of the stead flow is described b three parameters: an effective Mach number, a scaled energ loss parameter, and an integer number m: the number of half-oscillations in a mechanical analog that appears in this problem. In a stead shear flow the viscous heating is balanced b energ dissipation via inelastic collisions. This balance can have different forms, producing either a uniform shear flow or a variet of more complicated, nonlinear densit, velocit and temperature profiles. In particular, the model predicts a variet of multi-laer two-phase stead shear flows with sharp interphase boundaries. Such a flow ma include a few zeroshear solid-like) laers, each of which moving as a whole, separated b fluid-like regions. As we are dealing with a hard sphere model, the granulate is fluidized within the solid laers: the granular temperature is non-zero there, and there is energ flow through the boundaries of the solid laers. A linear stabilit analsis of the uniform stead shear flow is performed, and a plausible bifurcation diagram of the sstem, for a fixed m, is suggested. The problem of selection of m remains open. PACS numbers: 45.7.Mg, 83.5.Ax I. INTRODUCTION In spite of extensive experimental and theoretical studies of dense granular flows, a theoretical description of these flows remains challenging ]. Multi-particle contacts and friction, intrinsic in slow dense flows, invalidate the kinetic theor 2, 3]. Furthermore, even rapid dense flowsthat is, flows dominated b binar collisions) present significant difficulties for analsis. Dilute and moderatel dense flows of mono-disperse inelastic hard sphere fluids are describable, for not too high inelasticit of collisions, b Navier-Stokes hdrodnamics 3], which can be derived in a sstematic wa from the kinetic theor: a one-particle kinetic equation properl generalized to account for inelastic collisions 2, 3]. Hdrodnamic equations, that is conservation laws for the mass and momentum of the granulate, and a balance equation for the energ, ma still be valid, for a hard sphere fluid, at higher densities, after the disorder-order transition occurs. However, the respective constitutive relations are not derivable from a kinetic equation anmore. In this work we attempt to address this difficult and suggest a possible hdrodnamic description of a sheared rapid granular flow that exhibits crstallization. Shearinduced ordering in a dense granular medium has attracted much recent attention. It was investigated experimentall usuall for slow flows) b man groups 4, 5, 6]. The crstallization dnamics in inelastic hard sphere fluids has been also extensivel studied in MD simulations 7, 8]. In this work we will be dealing with an idealized model of inelastic particle collisions characterized b a constant coefficient of normal restitution r, and focus on a plane shear flow. The present work is a next step in a series of recent attempts of extending granular hdrodnamics of inelastic hard sphere fluids to high densities 9,, ]. Grossman et al. 9] investigated a prototpical sstem of inelastic hard disks at zero gravit, placed in a two-dimensional rectangular box, one wall of which serving as a thermal wall. Grossman et al. suggested an equation of state, granular heat conductivit and inelastic energ loss rate which interpolated between the dilute limit and the close vicinit of the hexagonal close packing, where freevolume arguments are available. The densit profiles, obtained b solving the hdrostatic equations numericall, were in good agreement with the results of MD simulations 9]. Meerson et al. ] considered a similar twodimensional granular sstem, but with gravit. The sstemwasdrivenfrombelowba thermal base. Emploing the constitutive relations suggested b Grossman et al. 9], Meerson and coworkers found stead-state densit and temperature profiles and observed good agreement with MD simulations, including the region of the levitating cluster, where the densit is ver close to that of hexagonal close packing ]. Bocquet et al. ] emploed hdrodnamic equations to model a granular shear flow where the densit approached the random close packing densit. The also compared their theor with experiment in a circular Couette flow ]. As a shear flow was present, Bocquet and coworkers had to specif, in addition to the rest of the
2 2 constitutive relations, the coefficient of shear viscosit. In experiment the granular temperature had been found to decrease more slowl, with an increase of the distance from the shear surface, than the velocit. To account for this finding, Bocquet et al. assumed that the shear viscosit diverges more rapidl, at random close packing, than the rest of the transport coefficients. We emplo in this work Navier-Stokes granular hdrodnamics for a description of a stead crstallized shear flow of an assembl of mono-disperse inelastic hard spheres under constant pressure. Our approach is similar to that taken b Bocquet et al. ], with two important differences. First, we focus on the ordered crstalline phase which ends at the crstalhexagonal or facecentered cubic) close packing densit φ fcc, whereas Bocquet et al. considered the metastable disordered phase which ends at the random close packing densit φ rcp. Second, we assume that the shear viscosit diverges at a smaller densit than the rest of transport coefficients, see below. This assumption brings about the possibilit of a two-phase stead flow. In general, a granular shear flow reaches a stead state when the viscous heating makes up for the energ dissipation via inelastic collisions. As we show here, this balance can be achieved in different was, producing either a uniform shear flow, or a variet of flows with nonlinear densit, velocit and temperature profiles. Working in the range of densities beond the melting point, we determine the phase diagram of dense stead flows in terms of three parameters: an effective Mach number, a scaled energ loss parameter, and an integer number m: the number of half-oscillations in a mechanical analog that appears in this problem. There are regions on this phase diagram where two or more different stead flow solutions are possible for the same values of the parameters. To get an additional insight, we perform a linear stabilit analsis of the uniform stead shear flow. Based on these results, we suggest a plausible bifurcation diagram of the sstem which, in some region of the parameter space, describes bistabilit and hsteresis. The rest of the paper is organized as follows. In Sec. 2 we introduce our shear-induced crstallization scenario and the governing hdrodnamic equations and constitutive relations. Section 3 describes our model of a zerogravit constant-pressure shear flow and focuses on the analsis of a crstallized stead shear flow in different regimes. Section 4 presents a linear stabilit analsis of the uniform shear flow solution and suggests a plausible bifurcation diagram of the sstem. Details of the linear stabilit analsis are presented in the Appendix. Section 5 includes a brief discussion and summar of our results. II. PRESSURE-DENSITY DIAGRAM AND HYDRODYNAMIC EQUATIONS Figure depicts, in the coordinates volume fraction - pressure, the phase diagram of a homogeneous macroscopic sstem of elastic hard spheres 2]. The volume fraction φ = π/6)d 3 n where n is the particle number densit, and d is the particle diameter) changes from zero to φ fcc = 2π/6.74, the densit of crstal either hexagonal, or face-centered cubic) close packing, the densest possible packing of spheres in three dimensions. The phase diagram includes four branches. The disordered, or gas/liquid branch starts at φ = and continues until the freezing point that occurs at φ.494. Then this branch splits into two branches. The horizontal constant pressure) branch describes coexistence of the disordered and ordered phases. It starts at the freezing point and ends at the melting point at φ.545. At larger volume fractions the sstem is in the ordered crstalline phase that ends at the densit of the crstal close packing φ fcc. The last branch is a metastable extension of the gas/fluid branch. It starts at the freezing point and ends at random close packing at φ rcp.64. Clearl, each branch is described b a separate equation of state. The phase diagram, presented in Fig., appears in the context of a homogeneous sstem in either equilibrium, or a metastable state. Granular sstems are usuall inhomogeneous and, no less important, the are intrinsicall far from equilibrium, due to inelastic collisions between particles. We will deal, however, with a small inelasticit of particle collisions and, in the spirit of kinetic theor 2, 3], assume that the sstem is everwhere close to local thermodnamic equilibrium, so the phase diagramand the constitutive relations presented below) are valid locall. Let the fluid densit be sufficientl large, so that the volume fraction is everwhere above the melting point. There are two possible phases here. One of them is the disordered phase the the metastable branch in Fig. ), the other one is the crstallized phase. Each of the two phases, disordered and ordered, can have inhomogeneous temperature and densit profiles, without violating the constanc of the pressure. Furthermore, domains of disordered phase can in principle coexist with domains of the ordered phases again, without violating the constanc of the pressure). Shear-induced crstallization apparentl represents a non-equilibrium) phase transition, so that the metastable disordered branch gives wa, everwhere, to the stable ordered branch. This work addresses ver dense flows, where the volume fraction is close to that of the crstalline close packing, φ fcc φ φ fcc. We assume that the ordered crstalline branch of the phase diagram see Fig. ) has alread won in the competition with the disordered branch. To describe a flow in this phase, one needs hdrodnamic equations. These represent the mass conservation dn dt the momentum conservation n dv dt n v =, ) = P ng, 2)
3 3 Pressure P) Metastable branch Freezing Volume fraction φ) To RCP Melting.74 To FCC Close packed) FIG. : Phase diagram of elastic hard sphere fluid, see text. and the energ balance 3 2 n dt dt = QP : v Γ, 3) where Γ is the energ loss rate due to the inelasticit of binar collisions. Here nr, t) is the number densit of grains, Tr,t) is the granular temperature, vr,t) is the mean flow velocit, P is the stress tensor, g is the gravit acceleration, Q is the heat flux, and d/dt = / tv is the total derivative. In the following we will put the particlemasstounit. ThestresstensorPcanbewritten as P = pn,t)µn,t)trd)]i2ηn,t) ˆD, 4) where D = /2) v v) T] 5) is the rate of deformation tensor, ˆD = D trd)i 6) 3 is the deviatoric part of D, and I is the identit tensor. We assume that the heat flux Q is given b the Fourier law Q = κn,t) T. 7) Similarl to dilute 3] and moderatel dense 4] rapid granularflows, therecanbeanadditionaltermineq.7), proportional to the densit gradient. This term must vanish as r, and it can be neglected in the nearl elastic limit q that we are interested in throughout this paper. To make the formulation complete, one needs constitutive relations: the equation of state p = pn,t) and the dependences of the transport coefficients µ, η and κ, and of the energ loss rate Γ on n and T. For small and moderate densities, these relations can be derived from the Boltzmann or Enskog equation, properl generalized to account for inelastic collisions 2, 3]. For ver large densities, that we are interested in, one can use free volume arguments 5]. The resulting equation of state is p = p n 2 c T n c n, 8) where n c = 2/d 3 is the particle number densit at close packing, and p is a numerical factor of order unit. Here and in the followingwewill write n c insteadofn fcc facecentered cubic). The temperature dependence in Eq. 8) is exact. We assume that all transport coefficients except for the shear viscosit η, see below), diverge like n c n) near the close packing densit 9]. Therefore, one can write the bulk viscosit coefficient µn,t), the thermal conductivit κn, T), and the energ loss rate Γn,T) in the following form µ = µ n c T /2 n c n)d 2, κ = κ n c T /2 n c n)d 2, Γ = Γ r 2 ) n 2 c T 3/2 n c n)d. 9) Equations9)arevalidinthelimitn c n n c. Thetemperature dependences are exact. The numerical factors µ, κ, and Γ are of order unit and presentl unknown; the can be found in MD simulations. The same tpe of divergence of the transport coefficients again, except for the coefficient of shear viscosit η) was assumed in Ref. ], but in the vicinit of the random close packing densit. There is a recent evidence in the literature that the shear viscosit coefficient η of elastic hard disk fluid grows faster with the densit, at high densities, than the rest of the transport coefficients 6]. We believe that this behavior remains qualitativel correct also for three dimensional sstems. Bocquet et al. ] accounted for this fact in their description of the plane shear flow near the random close packing densit. The proposed, b analog with the behavior of supercooled liquids above the glass transition, that the shear viscosit coefficient diverges at random close packing densit, but with a larger exponent: η n rcp n) β, β >. In our model of a crstallized shear flow we suggest a different approach. We will accommodate a recent finding of Luding et al. 6] that the shear viscosit coefficient diverges like n c n) at a densit n c < n c : η = η n c T /2 n c n)d 2, ) where η is a numerical factor of order unit, which is presentl unknown. Divergence of the shear viscosit implies that the fluid is jammed on a macroscopic length scale. While a shear flow in a macroscopicall jammed sstem is impossible, the sstem of hard spheres ma still
4 4 have finite temperature, pressure, heat conduction, and collisional energ loss. There are arguments in the literature see, e.g., Santos et al. 7]) that stead sheared states of granular fluids are intrinsicall non-newtonian that is, require a description beond Navier-Stokes order). Wh are such effects not included in our description? The reason is that non-newtonian effects arise when one treats the inelasticit of particle collisions in a non-perturbative wa postulating that the Boltzmann equation, generalized to inelastic collisions, remains applicable for finite inelasticit). In this paper we work in the limit of nearl elastic particle collisions. In the first order in the inelasticit one does not need to take into account an inelasticit corrections to the transport coefficients for example, to the shear viscosit). The onl place where the inelasticit enters in this leading-order theor is the inelastic loss term Γ in the energ equation 3). III. STEADY SHEAR FLOW CLOSE TO CRYSTALLIZATION Now let us introduce the flowsetting we will be dealing with in the rest of the paper. We consider a plane Couette geometr. The model sstem is infinite in the horizontal x) direction and driven b the upper wall = H that moves horizontall with velocit u. The lower wall = is at rest of course, in the absence of gravit the upper and lower walls are interchangeable). The height H of the laer in this setting is not fixed. Instead, the upper wall is maintained at a constant pressure P, like in a recent experiment b Gollub s group 5]. As alread mentioned above, we assume a ver dense flow, so that the volume fraction is close everwhere to that of close packing: φ fcc φ φ fcc. As we consider a stead horizontal motion, the number densit n, the granular temperature T and the horizontal velocit u depend onl on the vertical coordinate. Then it follows from the -component of the momentum equation 2) that the pressure is constant throughout the sstem: p) = P. ) Now we write down the x-component of the momentum equation 2) and rewrite the energ balance equation 3): d η du ) =, d d d κ dt ) ) 2 du η Γn,T) =, 2) d d d where the constitutive relations are given b Eqs. 8)- ). Equations 8)-2) must be complemented b boundar conditions. We will assume rough walls and no-slip boundar conditions for the horizontal velocit: u = ) =, u = H) = u. The problem of evaluation of the granular heat flux at the rough walls was addressed b Chou 8]. He considered a model where inelastic spheres were driven b walls with attached halfspheres bump walls.) Extending an earlier treatment b Richman 9], Chou calculated the heat flux and the slip velocit at the boundaries. He showed that the heat flux at the boundaries can be positive or negative, depending on whether the slip work is larger or smaller than the energ loss due to inelastic particle collisions with the walls. For some values of parameters the total heat flux to the boundaries vanishes 8]. For simplicit, we will assume a vanishing heat flux and therefore prescribe dt/d = ) = dt/d = H) =. The total number of particles is conserved, which ields a normalization condition for the densit: H n)d = N, where N isthe numberofparticlesperunit areainthe xz plane. Let us rescale the vertical coordinate b the a priori unknown) sstem height H, the horizontal velocit u b the upper plate velocit u, the densit n b the closepacking densit n c, and the temperature T b the ratio P /n c of the constant applied pressure and close-packing densit. Rewriting the equations in the scaled form, we obtain η MT /2 κ λ n d T /2 d du d T p n =, ) du =, λ n d ) d T /2 dt d n d ) 2 Γ 2RT 3/2 =, 3) κ n where M = n c u 2 /P is the square of the effective Mach number of the flow, and λ = n c/n c is a small positive numerical factor. For simplicit, the three presentl unknown constants p = O), η /κ = O) and Γ /κ = O) are taken to be unit in the following. As the sstem height is unknown a priori, the scaled quantit R = r 2 )H 2 /2d 2 ) must be determined from the solution of the problem. From the first of Eqs. 3) we find a simple relation between the scaled densit and temperature: T) = n). Then, introducing a convenient auxiliar variable ǫ) = n)] /2, we rewrite the remaining Eqs. 3) as du 2c d = M ǫ 2 λ ǫ, 4) d 2 ǫ cλ c R)ǫ =, 5) d2 ǫ where c is an unknown constant to be found from the solution of the problem. The boundar and normalization conditions are u = ) =, u = ) =,
5 5 dǫ dǫ = ) = = ) =, d d ǫ) 2 ] d = f R/2, 6) where f = r 2 ) /2 Nd 2 /2. The stead flow equations 4)-6) include two known) scaled parameters: M and f. There are five unknown parameters in the problem: c, R and three arbitrar constants which determine the solutions of the ordinar differential equations 4) and 5) for u) and ǫ), respectivel. Correspondingl, there are five conditions 6) to determine the unknown parameters. In the following, when solving the equations numericall, we put λ =.5. Before we get into a detailed description of the stead flowsolutionsofeqs. 4)-6), hereisanoverviewofthe results. Strikingl, there is an infinite number of stead flow solutions. The all can be parameterized b three parameters: the scaled numbers M and f and an integer number m =,2,... For a fixed m, there are three possible tpes of solutions in different regions of the phase diagram M,f), see Fig. 2. The simplest solution is the uniform, or linear shear flow that exists for an values of M and f. Here the velocit gradient, densit and temperature are all constant. In region B between the solid and dashed lines of Fig. 2), there is an additional onephase solution: the one with nonlinear profiles of densit, temperature and velocit. Finall, there is a multi-laer two-phase solution which exists in regions A and B, that is below the dashed line. The flow there is organized in several distinct laers. In the solid-phase laers the densit is larger than the critical densit n c at which the shear viscosit diverges. Therefore, these laers move as a whole with a velocit u = const, while the temperature and densit profiles are non-trivial. In the fluid phase laers there is a mean flow with non-trivial profiles of the velocit, temperature and densit. The densit, temperature, heat flux and velocit are continuous at the interface between the laers. This means in particular, that there is no macroscopic flow in the bottom laer of the granulate, and the inelastic energ losses there are balanced b the conduction of heat from the flowing) top laer. Now let us consider these solutions in more detail. We start with the simplest one: the uniform shear flow. Here ǫ is independent of, and Eq. 5) ields ǫ = cλ c R. Substituting this value into Eq. 4) and integrating in, we obtain a linear velocit profile 2λ u) = R. 7) M c R) Using the normalization condition and the condition FIG. 2: Hdrodnamic phase diagram of the stead flow solutions for m = the upper panel) and m = 2 the lower panel) in terms of the scaled parameters f and M. In the regions A there are: a multi-laer flow and an unstable) uniform flow. In the regions B there are: a multi-laer flow, a non-uniform flow and an unstable) uniform flow. In the regions C there is onl a stable) uniform flow. In the forbidden regions the granulate is not dense enough for our theor to be valid. Plausible bifurcation diagrams, corresponding to the dotted arrows f = 2, are shown, for the cases of m = and m = 2 in Fig. 7. u = ) =, we obtain cλ c R = f R and c = R 2λR2 M. 8) Solving these two algebraic equations, we calculate R R = f 2 2 λ) 2 M λ) f 2 2M λ) ) ] /2 f 2 9) and obtain the scaled velocit profile and the constant value of ǫ: u ) =, ] /2 λ)k ǫ = ǫ = λ, 2) k 2k ) /2
6 6 where k = M λ)f 2. Notice that ǫ alwas satisfies the condition ǫ > λ /2 or, in the original variables, n < n c: the densit of the uniform flow is smaller than the critical densit at which the shear viscosit diverges. Note also that although the uniform shear flow solution Eq. 2)] formall exists everwhere, the assumption ǫ demands k. Therefore, the solution is valid for M f 2 / λ). This inequalit breaks down in the forbidden regions of the phase diagram, see Fig. 2. The nonlinear solution, which exists in region B, can be found numericall. We used the following numerical procedure, realized in Matlab. Let us denote ǫ ǫ = ). For fixed R and f, we first solve Eq. 5) b varing parameters ǫ and c and demanding dǫ/d = ) = and the normalization condition. This procedure ields ǫ). Then we solve Eq. 4) and find the velocit profile u) and the corresponding value of parameter M from the condition u = ) =. An important insight into the nature of this solution is provided b a mechanical analog following from Eq. 5). Let ǫ be a coordinate, while the vertical coordinate be time. Then Eq. 5) describes a Newtonian particle oscillating in the potential well Uǫ) = /2)c R)ǫ 2 cλlnǫ. The boundar conditions select a famil of solutions. A tpical solution includes an integer number m of halves of the oscillation period, and a completion of these oscillations takes the sstem a unit time, mτ/2 = see Eq. 2)], where τ is the oscillation period. Figure 3 shows an example of particle in the potential well, while the resulting spatial profiles of ǫ), u) and n) for a one-half of a full oscillation m = ), and for a full oscillation m = 2), are shown in Fig. 4. At a fixed m, this solution exists in region B of the phase plane of parameters f and M. This region is limited b two curves M = Mf) see Fig. 2). The first curve is obtained as one approaches the bottom of the potential well. Here the oscillation amplitude that is, the densit contrast in the sstem) vanishes. Expanding the potential Uǫ) near the bottom of the potential well ǫ = cλ/c R)] /2, we can calculate the oscillation period τ: τ = 2 ) /2 2 m = π. 2) c R The bottom of the well corresponds to a uniform shear flow, but Eq. 2) must be obeed arbitraril close to the bottom of the well. Using Eq. 2) and the second equation of Eqs. 8), we obtain R = mπ/2)m/λ) /2 and c = R m 2 π 2 /2. Substituting it into the first equation of Eqs. 8), we obtain the lower boundar of region B in terms of f and M: f = mπ ) /2 M 2 λ ) /4 ] λ λm)/2. 22) mπ For m = and m = 2 these curvesare shown b the solid lines in Fig. 2. The determine the boundar between the regions A and B for each m. U.9.8 Region B ε ε = λ / FIG. 3: The mechanical analog for Eq. 5). Shown is a trajector the dashed line) of a classical particle in the potential well Uǫ) the solid line, see text) which corresponds to regions B in Fig. 2. See Figure 4 for the resulting spatial profiles, which correspond to a one-half of the oscillation m = ) and a full oscillation m = 2, here the potential changes as the parameter M is different.) The oscillating solutions exist onl inside the potential well, and onl for ǫ = ) > λ /2, when the densit is everwhere below the critical densit at which the shear viscosit diverges. The parameters for this figure are M =.558 and f = 2. ε u n ε u n FIG. 4: The profiles of ǫ) = n)] /2 the upper row), the scaled velocit u) the middle row), and the scaled densit n) the bottom row), which correspond to region B in Figure 2. The profiles in the left column respectivel, in the right column) correspond to m =, a one-half of a full oscillation respectivel, to m = 2, a full oscillation) in the potential well, see Fig. 3. The parameters are f = 2 and M =.558 the left column), and f = 2 and M =.5 the right column). The second limiting curve plotted b the dashed lines in Fig. 2 for m = and m = 2) is obtained numericall b putting ǫ = λ /2 see Fig. 3). This curve determines the boundar between the regions B and C. Now let us consider the most interesting multi-laer two-phase famil of solutions, parameterized b the same
7 7 number m that now can be associated with the number of solid laers. Due to continuit of the velocit field, there is no macroscopic flow in the bottom laer of the granulate for an m, while the inelastic energ losses there are balanced b the heat conduction from the fluid-like top laer. A tpical flow here consists of m zero-shear solid-like) laers, each of which moving as a whole, separated b fluid-like regions. As we are dealing with a hard-sphere model, the granulate is fluidized within the solid laers: the granular temperature is non-zero, and there is energ flow through the boundaries of the laers. Let us first consider the m = solution. Here the onl solid-like laer is at the bottom, and it is at rest. Putting c = in Eqs. 5) and 4) we obtain ǫ) = ǫ coshr /2 ). The height of the bottom laer is determined from the condition ǫh) = λ /2. At the interface between the two laers we demand continuit of the densit and, therefore, of the temperature), of the heat flux, and of the velocit. Similarl to the procedure performed in region B, we solve the problem numericall b shooting in two parameters ǫ and c for the top laer) for fixed values of R and f, and then calculate the profiles and the respective value of the parameter M. Tpical profiles ǫ), u) and n) for m = are shown in Fig. 5, the left column. A multi-laer solution m > is obtained in a similar wa, b demanding that the velocit in an solid laer is constant, while the densit, temperature, heat flux and velocit are continuous at all interfaces between the laers. Tpical profiles of ǫ), u) and n) for m = 2 are shown in Fig. 5, right panel. ε u n ε u n FIG. 5: The profiles of ǫ) the upper row), the velocit u) the middle row) and the densit n) the bottom row) for the two-phase solution, for m = the left column) and m = 2 the right column). The solutions for the fluid solid) laers are shown b the solid dotted) lines. The parameters are f = 2 and M =.492 the left column), and f = 2 and M =.48 the right column.) At this point let us return, for a moment, to the case where the viscosit singularit is assumed to be at the same densit n c as the rest of transport coefficient. In this case λ =, and the onl possible solution is a uniform shear flow. Indeed, as ǫ must be positive, the onl acceptable solution of Eq. 5) with λ =, which obes the no-flux boundar conditions, is ǫ = const and c = R. Then the normalization condition the last of Eqs. 6)] ields ǫ 2 = f/r /2. One can obtain this solution directl b putting λ = in Eq. 2). Therefore, the assumption of a nonzero λ the specific form of viscosit divergence, Eq. )], is crucial for the existence of non-trivial solutions. One can see that, at fixed m, two different kinds of stead flow solutions exist in region A, for the same values of parameters M and f. Furthermore, three different kinds of stead flow solutions exist in region B, again for the same values of parameters M and f see Fig. 2). What is the selection rule for these solutions? We give a partial answerto this question in the next section b performing a linear stabilit analsis of the uniform dense shear flow. Then we suggest plausible bifurcation diagrams of the sstem for different m. IV. LINEAR STABILITY AND BIFURCATIONS Our linear stabilit analsis of the uniform flow deals with the full set oftime-dependent) hdrodnamic equations )-3) and constitutive relations 8)-). The details of linear stabilit analsis are shown in Appendix. Adding a small perturbations to the uniform shear flow and linearizing the equations, we finall arrive at a quadratic characteristic equation for the growth/damping rate Γ as a function of parameters and the wave number k Eq. A.4), see Appendix.] The two roots of this equation are real and correspond to two different collective modes of the sstem. One of them alwas decas. The other one is a purel growing mode for sufficientl small k, that is for long-wavelength perturbations. At small k we obtain: Γ = 2R 2 λk 2 >, 23) M Lǫ M f R /2 ) where L = 2M) /2 H/d. For sufficientl large k, Γ is negative: the heat conduction suppresses the instabilit. The critical wave number for the instabilit is k = 2Rλ /2 /M /2. Thedependenceofthescaledgrowth rate Γ on the scaled wave number k is shown in Fig. 6. Using Eq. 9), we obtain k = f 2 λ /2 2 λ) 2 M /2 M λ) f 2 2M λ) ) ] /2 f 2. 24) Notice that the time scale separation, emploed here for the reduction of the order of the dispersion equation see Appendix), breaks down when ǫ approaches λ /2. As one can see from Eq. 2), this happens, for a fixed
8 8..2. Γ k δ..2. a) M M * M FIG. 6: The scaled growth rate Γ, calculated from Eq. A.4), versus the scaled wave number k. The long-wavelength modes with k < k the value of k is indicated b the asterisk) are unstable. The parameters f = 2, M =.2, λ =.5, and L = correspond to region A in Fig. 2a and to region C in Fig. 2b. f, when M becomes sufficientl small. Importantl, the final results 23) and 24) remain valid in the general case, as we obtained from the full, unreduced fourthorder dispersion equation for Γk). The wave number k of the perturbation is quantized b the boundar conditions. Indeed, the velocit perturbations must vanish at the upper and lower plates: u =,t) = u =,t) = v =,t) = v =,t) =. These conditions ield a discrete spectrum of wave numbers: k = πm, where m =,2,3,... Therefore, the instabilit threshold is k = π. This equation determines a curve fm) on the f,m) plane. One can see that this curve coincides with the curve fm) determined b Eq. 22) for m = : the borderline between regions A and B in Fig. 2. Similarl, the threshold for m = 2 perturbation k = 2π coincides with the curve fm) determined b Eq. 22) for m = 2, and so on. This is not entirel surprising: as the instabilit is aperiodic, its threshold is provided b the marginal stabilit condition. It is convenient to characterize each stead flow solution b the maximum densit contrast it predicts. In terms of the auxiliar variable ǫ, we can define δ as the maximum change in ǫ throughout the sstem for m = this gives δ ǫ = ) ǫ = ).] We can calculate δ as a function of parameters f and M for each stead flow solution. Figure 7 shows all the resulting branches of the solution for m = and m = 2) for a fixed f, while M serves as a control parameter. Furthermore, there is an infinite number of bifurcation diagrams for multi-laer solutionsfor allintegersm > ), whichcan be computed in a similar wa. In order to determine the bifurcation diagram of the sstem at a fixed m, one needs to perform a linear stabilit analsis of each branch of the solution. Unfortunatel, such an analsis is quite cumbersome for the so- δ FIG. 7: A conjectured bifurcation diagrams of the stead flow for f = 2 and λ =.5. Shown is the densit contrast δ versus M for different stead flow solutions for m = the upper panel) and m = 2 the lower panel). The stable and unstable branches are indicated b the solid and dashed lines, respectivel. We proved that the uniform flow δ = ) is unstable at M < M. The dash-dotted line in the upper panel is the asmptote of the branch bifurcating from the point M = M : δ = M M ) /2. lutions with δ. It is natural to assume that bifurcation, which occurs in the vicinit of the point M = M, at which the uniform shear flow the one with δ = ) changes from unstable at M < M f)] to stable at M > M f)], is an inverse pitchfork bifurcation. Based on this assumption, we obtain a simple bifurcation diagram, presented in Fig. 7. The stable solutions are indicated b solid lines, the unstable solutions b dashed lines. Let us consider the case of m = and follow the uniform shear flow solution for which δ = ) as M increases, starting from a small value, at fixed f. We pass from region A to region B of the hdrodnamic phase diagram see the dotted arrow in Fig. 2) via an inverse pitchfork bifurcation. Here the uniform shear flow solution becomes stable, while the bifurcating second solution the one with nonlinear densit and velocit profiles) must be unstable. Exploiting the mechanical analog see Fig. 3), we can find the asmptote of the unstable branch in the close vicinit of the bifurcation point: δ = Af)M M ) /2. Now let us move along the unstable branch δ and increase M. At some critical value of M, which depends on f, we reach the border be- M b)
9 9 tween regions B and C. Here the unstable branch ends, as for larger M region C) the onl possible solution is the uniform shear flow. There is, however, another solution whichexistsat smallerm. Figure7showsthis stable branch which corresponds to the two-phase solution in regions A and B see Fig. 2). This simple scenario predicts bistabilit at sufficientl small values of the parameter M, and a hsteresis on the interval M f),m f)], see Fig. 7. Note that the parameter L does not affect the bifurcation diagram, it onl sets the time scales of transient motions in the sstem. V. SUMMARY AND DISCUSSION We have considered here some aspects of shear-induced crstallization in a dense but rapid mono-disperse granular shear flow. We focused on a stead crstallized flow under a constant pressure and zero gravit. Assuming ver high densities, n c n n c, we emploed a version of the Navier-Stokes hdrodnamics for inelastic hard sphere fluid with ad hoc constitutive relations based on the free volume argument. In contrast to earlier works on rapid granular shear flow, we assumed that the shear viscosit coefficient η diverges at a densit smaller than the close packing densit n c, while the rest of the constitutive relations diverge at n = n c. We have determined the phasediagramofthe steadflowin termsofthreeparameters: the effective Mach number, the scaled inelastic energ loss parameter, and an integer number m. In a stead flow the viscous heating of the granulate is balanced b energ dissipation through inelastic collisions. This balance is achieved, in different parts of the phase diagram, in different was, producing either a uniform shear flow with constant velocit gradient, densit and temperature), or a flow with nonlinear velocit, densit and temperature profiles. In some regions in the phase diagram two or three different stead flow solutions are possible for the same values of the parameters. We performed a linear stabilit analsis of the uniform flow, and suggested a plausible bifurcation diagram of the flow at a fixed m, which predicts bi-stabilit and hsteresis. We are unable as et to find a selection principle that would prefer certain stead state solutions out of a multitude of solutions at different m. This non-trivial selection problem should be addressed in the future work. One of the predictions of this work is the existence and, we conjecture, stabilit) of two-phase solutions. The simplest solution of this tpe consists of a zero-shear solid-like) laer at the bottom and a flowing top laer. Though there is no mean flow in the bottom laer, the particles there undergo thermal motion, and the granular temperature and pressure are non-zero. As a result, there is energ transfer through the bottom laer. There also exist two-phase multi-laer solutions, where solid-like laers, each of which moving as a whole, are separated b fluid-like regions with nonlinear velocit, densit and temperature profiles. The existence of these solutions is a direct consequence of our assumption that the coefficient of shear viscosit η in Eq. ) diverges at a densit which is smaller than the close packing densit n c. A comparison of our results with those of Alam and coworkers is in order now. Alam et al. 2] investigated a plane Couette shear flow of inelastic hard sphere fluid, in a sstem with a fixed height, in a wide interval of densities: from the dilute limit to the random close packing densit. The emploed constitutive relations, all of which including the viscosit) diverge at the random close packing densit n r. The found instabilit of the uniform shear flow when the inelasticit of the particle collisions becomes large enough or, alternativel, when the fixed) sstem height exceeds a critical value for a fixed inelasticit. The uniform shear flow instabilit considered in our work is quite different: it requires viscosit divergence at a smaller densit than the rest of constitutive relations. Indeed, Eqs. 22) and 24) show that this instabilit disappears when λ goes to zero. Although the parameter λ was identicall zero in Ref. 2], instabilit was observed. Where does the difference come from? To remind the reader, the constitutive relations that we used assume a ver dense sstem. Respectivel, our results are valid onl at leading order in the parameter n c n)/n c. We checked that, if one takes into account onl the leading order terms in n r n)/n r in the equations of Alam et al. 2], the instabilit disappears. Its worth mentioning that the stead flow equations 4)-6) would not change if the crstal close packing densit n c were replaced b the random close packing densit n r. Such a formulation would follow from the assumption that the shear viscosit diverges at a densit smallerthann r. Itwouldpredictavarietofsteadstate solutions on the metastable branch, in the vicinit of the random close packing. We hope that the basic assumptions of our modelincluding the specific form of viscosit divergence) and its non-trivial predictions will be tested in MD simulations and in experiment on dense but rapid granular flow. Our work focused on a dense but rapid granular flow, assuming that the granulate is full fluidized. In experiment it should be easier to achieve this regime when the shear is ver large, so that the effective Mach number is of order unit see Fig. 7.) Most experiments with dense granular flows are performed with slow flows, where the effective Mach number is small. For example, in the experiment of Gollub s group 5] the parameter M was about 5 5. In this regime the particles far from the moving boundar are in persistent contact with each other 5], inter-particle friction is important see also 2]), and the model of inelastic hard spheres and the Navier-Stokes hdrodnamics) is inapplicable. In general, the model of inelastic hard sphere fluid is considered as a good approximation for dilute and moderatel dense flows. Its validit range for dense flows is not well known 3]. Indeed, the particle collision rate increases with the densit, so the assumption of instan-
10 taneous collisions, intrinsic in the model of hard spheres, ma become restrictive at high densities. Being aware of this limitation, we still believe that the model of inelastic hard spheres can capture some of the phsics of dense flow. Like in man other problems, it is useful to push this model to an extreme and analze its predictions some of which being quite unexpected as we have shown) in detail. We notice in passing that the two-phase flow predicted in this work resembles experimentall observed shear bands: localized regions of ordered granular flow, coexisting with essentiall immobile solid-like regions 22]. The future experimental and theoretical work should elucidate the exact conditions under which shear-induced crstallization develops in granular flows. Though crstallization under shearing is well documented in MD simulations 7, 8], a recent experiment 23] showed that, under certain conditions, shearing can lead to disorder. Finall, wedid notattempt to describeinthis workthe shear-induced crstallization process. Such a description is beond the reach of theor as of present. A promising approach to this problem should deal, in addition to the hdrodnamic fields, with an order parameter field and its dnamics. For slow granular flows such a description is now emerging 24]. Acknowledgments We acknowledge useful discussions with J.P. Gollub, S. Luding, P.V. Sasorov, J.-C. Tsai, and A. Vilenkin. The work was supported b the the Israel Science Foundation grant No. 7/5) and b the German-Israel Foundation for Scientific Research and Development Grant I /23). APPENDIX: DETAILS OF LINEAR STABILITY ANALYSIS Let u and v be the velocit components in the x and directions, respectivel. We assume that the small perturbations to the uniform stead shear flow do not break the smmetr of the flow in the x and z directions: n,t) = n n,t), T,t) = T T,t), u,t) = u )u,t), v,t) = v,t), A.) where index denotes the uniform stead shear flow, see Eqs.2), whilet = n = ǫ 2. Linearizingthehdrodnamic equations with respect to small perturbations, we arrive at n t n v =, ) u Ln t v du = d { = T /2 λ n v Ln t = L p M T /2 n u 4T /2 T 2T 3 λ n ) 3 2 Ln T t L v = M T/2 du λ n d 2RT3/2 n 2 u 3T n 2T n ) n λ n v ] T /2 T n, ) ) ] } du, d ) )] T n du 2T λ n d ), A.2) where L = 2M) /2 H/d, the time is measured in the units of H/u, and the value of H is determined b the unperturbed flow, that is b the uniform stead flow solution. One can check that the new parameter L which is absent in the stead state problem, but enters the linear stabilit analsis), is full determined b the scaled parameters f and M and b the restitution coefficient r. A further simplification emplos time scale separation. Let us compare the characteristic time scales of the problem. The acoustic time scale is τ = Hǫ /P /n c ) /2, the heat conduction time scale is τ 2 = τ H/d), the energ loss time scale is τ 3 = τ H/dR)], and the viscous time is τ 4 = τ H/d)ǫ 2 λ)/ǫ 2. If the densit is not too close to the densit at which the shear viscosit diverges, that is if H/d)ǫ 2 λ)/ǫ 2, one can separate the different time scales and eliminate the acoustic-like) fast modes. This is equivalent to the assumption that the perturbations evolve in pressure equilibrium with the surroundings 25]. Using the condition p / t = instead of the full momentum equation the thirdequationin Eqs.A.2)], weobtaint = n. Then, differentiating the second equation of Eqs. A.2) with respecttoandsubstituting v / fromthefirstequation of Eqs. A.2), we arrive at L 2 ǫ ) ) u = 2 { ǫ 2 L n t M ǫ ǫ 2 λ t n ǫ 2 u ǫ 2 λ n 2ǫ 2 n )]}, ǫ 2 λ 2 ] 3 ǫ ) = 2 ǫ 2 ǫ 2 2 u n 2ǫ 2 n ǫ 2 λ ) n )] Rn ǫ =, where we substituted du /d = and T = n = ǫ 2. Consider a single Fourier mode of the perturbation: n,t) = ñexpγtik), u,t) = ũexpγtik),
11 where k is the wave number. Looking for nontrivial solutions, we obtain deta) =, where the elements of the matrix A are A = ΓL 4R2 ǫ k 2 M 2 Rk2, M ǫ A 2 = ΓLkf 2Rǫ k 3 R /2 M, 3f A 2 = ΓL A 22 = 4Rkǫ. R/2 2R/2 f ) k2 ǫ 4R2 ǫ M, A.3) where the coefficients are ) a 2 = Lk 2 3f2, 2R a = 4RLkǫ f ) R/2 M Lf k3 3Rǫ 2 ) R /2 ǫ M 2R2 2 ǫ M f 2, ) a = 2Rk3 k 2 4R2 λ. M M Equation A.3) ields a quadratic equation for Γ: a 2 Γ 2 a Γa =, A.4) ] Challenges in Granular Phsics, edited b T.C. Halse and A. Metha World Scientific, Singapore, 22). 2] J. T. Jenkins and M. W. Richman, Phs. Fluids 28, ); Arch. Ration. Mech. Anal. 87, ). 3] C. S. Campbell, Annu. Rev. Fluid Mech. 22, 57 99); I. Goldhirsch, Annu. Rev. Fluid Mech. 35, ). 4] M. Nicolas, P. Duru, and O. Pouliquen, Eur. Phs. J. E 3, 39 2). 5] J.-C. Tsai, G. A. Voth, and J. P. Gollub, Phs. Rev. Lett. 9, ); J.-C. Tsai and J. P. Gollub, Phs. Rev. E 7, ); ibid 72, ). 6] N. W. Mueggenburg, Phs. Rev. E 7, 33 25). 7] P. Zamankhan, T. Tnjala, W. Polashenski, P. Zamankhan, and P. Sarkomaa, Phs. Rev. E 6, ). 8] M. Alam and S. Luding, Phs. Fluids 5, ). 9] E.L. Grossman, T. Zhou, and E. Ben-Naim, Phs. Rev. E 55, ). ] L. Bocquet, W. Losert, D. Schalk, T. C. Lubensk, and J. P. Gollub, Phs. Rev. E 65, 37 22); L. Bocquet, J. Errami, and T. C. Lubensk, Phs. Rev. Lett. 89, ). ] B. Meerson, T. Pöschel, and Y. Bromberg, Phs. Rev. Lett. 9, ). 2] M.D.RintoulandS.Torquato, J. Chem.Phs.5, ); Phs. Rev. Lett. 77, ); Phs. Rev. E 58, ); R. K. Bowles and R. J. Speed, Phsica A 262, ). 3] J.J. Bre, F. Moreno, and J.W. Duft, Phs. Rev. E 54, ); N.SelaandI.Goldhirsch, J.FluidMech.36, 4 998); J.J. Bre, J.W. Duft, C.S. Kim, and A. Santos, Phs.Rev.E58, ); R.Soto, M.Mareschal, and D. Risso, Phs. Rev. Lett. 83, ). 4] J.F. Lutsko, Phs. Rev. E 72, ). 5] Z.W. Salsburg and W.W. Wood, J. Chem. Phs. 37, ), and references therein. 6] R. Garcia-Rojo, S. Luding, and J. J. Bre, cond-mat/567. 7] A. Santos, V. Garzo, and J.W. Duft, Phs. Rev. E 69, ). 8] C.S. Chou, Phsica A 287, 27 2). 9] M.W. Richman, Acta Mech. 75, ). 2] M. Alam and P. R. Nott, J. Fluid Mech 377, ); M. Alam, V. H. Arakeri, P. R. Nott, J. D. Goddard, and H. J. Herrmann, J. Fluid. Mech. 523, ). 2] S.J. Moon, J.B. Swift, and H.L. Swinne, Phs. Rev. E 69, 33 24). 22] D.M. Mueth, G.F. Debregeas, G.S. Karczmar, P.J. Eng, S.R. Nagel, and H.M. Jaeger, Nature 46, 385 2). 23] K.E. Daniels and R.P. Behringer, Phs. Rev. Lett. 94, 68 25). 24] D. Volfson, L.S. Tsimring, and I.S.Aranson, Phs. Rev. Lett. 9, ); Phs. Rev. E 68, 23 23); ibid 69, ). 25] B. Meerson, Rev. Mod. Phs. 68, ).
L O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S
L O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S Yoël FORTERRE, Olivier POULIQUEN Institut Universitaire des Sstèmes Thermiques et Industriels (UMR 6595 CNRS) 5, rue Enrico Fermi 13453
More informationModule 1 : The equation of continuity. Lecture 4: Fourier s Law of Heat Conduction
1 Module 1 : The equation of continuit Lecture 4: Fourier s Law of Heat Conduction NPTEL, IIT Kharagpur, Prof. Saikat Chakrabort, Department of Chemical Engineering Fourier s Law of Heat Conduction According
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 13 Nov 2003
Europhysics Letters PREPRINT arxiv:cond-mat/31133v1 cond-mat.soft] 13 Nov 23 Oscillatory instability in a driven granular gas Evgeniy Khain and Baruch Meerson Racah Institute of Physics, Hebrew University
More informationChapter 2 Basic Conservation Equations for Laminar Convection
Chapter Basic Conservation Equations for Laminar Convection Abstract In this chapter, the basic conservation equations related to laminar fluid flow conservation equations are introduced. On this basis,
More information7. TURBULENCE SPRING 2019
7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport
More informationUniversality of shear-banding instability and crystallization in sheared granular fluid
Under consideration for publication in J. Fluid Mech. 1 Universality of shear-banding instability and crystallization in sheared granular fluid MEEBOOB ALAM 1, PRIYANKA SUKLA 1 AND STEFAN LUDING 2 1 Engineering
More informationSubmarine sand ripples formation in a viscous fluid: 2D and 3D linear stability analysis
Marine Sandwave and River Dune Dnamics 1 & April 4 - Enschede, the Netherlands Submarine sand ripples formation in a viscous fluid: D and 3D linear stabilit analsis V. Langlois (1) and A. Valance (1) Groupe
More informationThe Attractors of the Rayleigh-Bénard Flowof ararefied Gas
The Attractors of the Raleigh-Bénard Flowof ararefied Gas S.Stefanov, V. Roussinov and C. Cercignani Inst. of Mech., Bulg. Acad. of Sciences, Sofia, Bulgaria Dipart. di Matematica, Polit. di Milano, Ital
More informationCONSERVATION OF ENERGY FOR ACONTINUUM
Chapter 6 CONSERVATION OF ENERGY FOR ACONTINUUM Figure 6.1: 6.1 Conservation of Energ In order to define conservation of energ, we will follow a derivation similar to those in previous chapters, using
More informationExtended kinetic theory applied to dense, granular, simple shear flows
Noname manuscript No. (will be inserted by the editor) Diego Berzi Extended kinetic theory applied to dense, granular, simple shear flows Received: date / Accepted: date Abstract We apply the Extended
More informationChapter 4 Transport of Pollutants
4- Introduction Phs. 645: Environmental Phsics Phsics Department Yarmouk Universit hapter 4 Transport of Pollutants - e cannot avoid the production of pollutants. hat can we do? - Transform pollutants
More informationComparison of Kinetic Theory and Hydrodynamics for Poiseuille Flow
Journal of Statistical Phsics, Vol. 109, Nos. 3/4, November 2002 ( 2002) Comparison of Kinetic Theor and Hdrodnamics for Poiseuille Flow Yihao Zheng, 1 Alejandro L. Garcia, 2, 3 and Berni J. Alder 1, 3
More informationConservation of Linear Momentum for a Differential Control Volume
Conservation of Linear Momentum for a Differential Control Volume When we applied the rate-form of the conservation of mass equation to a differential control volume (open sstem in Cartesian coordinates,
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 6 Oct 2002
Competitive Clustering in a Bi-disperse Granular Gas René Mikkelsen, Devaraj van der Meer, Ko van der Weele, and Detlef Lohse Department of Applied Physics and J.M. Burgers Centre for Fluid Dynamics, University
More informationSymmetry Breaking and Turbulence in Perturbed Plane Couette Flow
Theoret. Comput. Fluid Dnamics (2002) 16: 91 97 Digital Object Identifier (DOI) 10.1007/s00162-002-0064- Theoretical and Computational Fluid Dnamics Smmetr Breaking and Turbulence in Perturbed Plane Couette
More informationContinuum Model of Avalanches in Granular Media
Continuum Model of Avalanches in Granular Media David Chen May 13, 2010 Abstract A continuum description of avalanches in granular systems is presented. The model is based on hydrodynamic equations coupled
More informationSaturation of inertial instability in rotating planar shear flows
J. Fluid Mech. (27), vol. 583, pp. 413 422. c 27 Cambridge Universit Press doi:1.117/s2211276593 Printed in the United Kingdom 413 Saturation of inertial instabilit in rotating planar shear flows R. C.
More informationSpotlight on the Extended Method of Frobenius
113 Spotlight on the Extended Method of Frobenius See Sections 11.1 and 11.2 for the model of an aging spring. Reference: Section 11.4 and SPOTLIGHT ON BESSEL FUNCTIONS. Bessel functions of the first kind
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationScattering Properties of Gas Molecules on a Water Adsorbed Surface
Scattering Properties of Gas Molecules on a Water Adsorbed Surface Hideki Takeuchi a, Koji Yamamoto b, and Toru Hakutake c a Department of Mechanical Engineering, Kochi National College of Technolog, Nankoku
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More informationESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundary Layer
ESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundar Laer References: Descriptive Phsical Oceanograph, Pickard and Emer Introductor Dnamical Oceanograph, Pond and Pickard Principles of Ocean Phsics, Apel
More informationVISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL CHANNEL THROUGH A POROUS MEDIUM
Volume 2, No. 1, Januar 214 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at http://www.jgrma.info VISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL
More informationLECTURE NOTES - VIII. Prof. Dr. Atıl BULU
LECTURE NOTES - VIII «LUID MECHNICS» Istanbul Technical Universit College of Civil Engineering Civil Engineering Department Hdraulics Division CHPTER 8 DIMENSIONL NLYSIS 8. INTRODUCTION Dimensional analsis
More informationJoule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder
Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 7-3, 04. Available online at www.jafmonline.net, ISSN 735-357, EISSN 735-3645. Joule Heating Effects on MHD Natural Convection Flows in Presence of
More informationOn the Effect of the Boundary Conditions and the Collision Model on Rarefied Gas Flows
On the Effect of the Boundar Conditions and the Collision Model on Rarefied Gas Flows M. Vargas a, S. Stefanov a and D. Valougeorgis b a Institute of Mechanics, Bulgarian Academ of Sciences, Acad. G. Bonchev
More informationOscillatory heating in a microchannel at arbitrary oscillation frequency in the whole range of the Knudsen number
Journal of Phsics: Conference Series Oscillator heating in a microchannel at arbitrar oscillation freuenc in the whole range of the Knudsen number To cite this article: O Buchina and D Valougeorgis J.
More informationCONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 742-761, 21. c Association for Scientific Research CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL R. Naz 1,
More informationModelling Axial Flow Generation and Profile Evolution in the CSDX Linear Device R. J. Hajjar, P. H. Diamond, G. R. Tynan, R. Hong, S.
Modelling Axial Flow Generation and Profile Evolution in the CSDX Linear Device R. J. Hajjar, P. H. Diamond, G. R. Tnan, R. Hong, S. Thakur This material is based upon work supported b the U.S. Department
More informationInternational Journal of Mathematical Engineering and Science ISSN : Volume 1 Issue 2
ISSN : 77-698 Volume Issue The Mathematical Structure and Analsis of an MHD Flow in Porous Media () () Naji Qatanani and Mai Musmar () Department of Mathematics, An-Najah National Universit Nablus-Palestine
More informationCHAPTER 7 EARTHQUAKE RESPONSE OF INELASTIC SYSTEMS. Base shear force in a linearly elastic system due to ground excitation is Vb
CHAPTER 7 EARTHQUAKE RESPONSE OF INELASTIC SYSTEMS Base shear force in a linearl elastic sstem due to ground excitation is Vb = ( A/ g) w where A is the pseudo-acceleration corresponding to natural period
More informationInstabilities of coupled density fronts
and their nonlinear evolution in the two-laer rotating shallow water model. Influence of the lower laer and of the topograph. B. Ribstein and V. Zeitlin Laborator of Dnamical Meteorolog, Univ. P.et M.
More informationDynamics and Patterns in Sheared Granular Fluid : Order Parameter Description and Bifurcation Scenario
Dynamics and Patterns in Sheared Granular Fluid : Order Parameter Description and Bifurcation Scenario NDAMS Workshop @ YITP 1 st November 2011 Meheboob Alam and Priyanka Shukla Engineering Mechanics Unit
More informationNumerical Analysis of Heat Transfer in the Unsteady Flow of a non- Newtonian Fluid over a Rotating cylinder
Bulletin of Environment, Pharmacolog and Life Sciences Bull. Env.Pharmacol. Life Sci., Vol 4 [Spl issue 1] 215: 318-323 214 Academ for Environment and Life Sciences, India Online ISSN 2277-188 Journal
More informationSteady Flow and its Instability of Gravitational Granular Flow
Steady Flow and its Instability of Gravitational Granular Flow Namiko Mitarai Department of Chemistry and Physics of Condensed Matter, Graduate School of Science, Kyushu University, Japan. A thesis submitted
More informationUnsteady free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion *
CHAPTER 4 Unstead free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion * 4. Nomenclature A : Suction parameter C : Dimensionless species
More informationBOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES
BOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES Vicente Garzó Departamento de Física, Universidad de Extremadura Badajoz, SPAIN Collaborations Antonio Astillero, Universidad de Extremadura José
More informationFlow through a Variable Permeability Brinkman Porous Core
Journal of Applied Mathematics and Phsics, 06, 4, 766-778 Published Online April 06 in SciRes http://wwwscirporg/journal/jamp http://dxdoiorg/0436/jamp0644087 Flow through a Variable Permeabilit Brinman
More informationInternational Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India
International Journal of Applied Mathematics and Phsics, 3(), Jul-December 0, pp. 55-67 Global Research Publications, India Effects of Chemical Reaction with Heat and Mass Transfer on Peristaltic Flow
More informationUNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD
Proceedings of the 3rd ASME/JSME Joint Fluids Engineering Conference Jul 8-23, 999, San Francisco, California FEDSM99-8 UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX
More informationTwo-sided Mullins-Sekerka flow does not preserve convexity
Differential Equations and Computational Simulations III J. Graef, R. Shivaji, B. Soni, & J. Zhu (Editors) Electronic Journal of Differential Equations, Conference 01, 1997, pp. 171 179. ISSN: 1072-6691.
More informationSimple shear flow of collisional granular-fluid mixtures
Manuscript Click here to download Manuscript: Manuscript_r1.docx 1 Simple shear flow of collisional granular-fluid mixtures 2 3 4 5 D. Berzi 1 1 Department of Environmental, Hydraulic, Infrastructure,
More informationUNIT 4 HEAT TRANSFER BY CONVECTION
UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional
More informationInternational Journal of Scientific Research and Reviews
Research article Available online www.ijsrr.org ISSN: 2279 543 International Journal of Scientific Research and Reviews Soret effect on Magneto Hdro Dnamic convective immiscible Fluid flow in a Horizontal
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationMechanics Departmental Exam Last updated November 2013
Mechanics Departmental Eam Last updated November 213 1. Two satellites are moving about each other in circular orbits under the influence of their mutual gravitational attractions. The satellites have
More information3 First order nonlinear equations
Separable equations 3 First order nonlinear equations d f dx = where f(x,) is not (in general) a linear function of. ( x, ) The equation is autonomous if there is no explicit dependence on the independent
More informationSolitary Wave Solutions of KP equation, Cylindrical KP Equation and Spherical KP Equation
Commun. Theor. Phs. 67 (017) 07 11 Vol. 67 No. Februar 1 017 Solitar Wave Solutions of KP equation Clindrical KP Equation and Spherical KP Equation Xiang-Zheng Li ( 李向正 ) 1 Jin-Liang Zhang ( 张金良 ) 1 and
More informationSound Propagation in Ducts
Sound Propagation in Ducts Hongbin Ju Department of Mathematics Florida State Universit, Tallahassee, FL.3306 www.aeroacoustics.info Please send comments to: hju@math.fsu.edu In this section we will discuss
More informationA Practical Variance Reduced DSMC Method
A Practical Variance Reduced DSMC Method H.A. Al-Mohssen and N.G. Hadjiconstantinou Mechanical Engineering Dept., MIT, 77 Massachusetts Ave, Cambridge, MA 02139, USA Abstract. The goal of this work is
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationOn Heat Transfer in case of a Viscous Flow over a Plane Wall with Periodic Suction by Artificial Neural Network U. K. Tripathy #1, S. M.
On Heat Transfer in case of a Viscous Flow over a Plane Wall with Periodic Suction b Artificial Neural Network U. K. Tripath #, S. M. Patel * # td. Professor & Head Dept of Mathematics V.S.S.Universit,
More informationScale-Adaptive Simulation (SAS) Turbulence Modeling. F.R. Menter, ANSYS Germany GmbH
Scale-Adaptive Simulation (SAS) Turbulence Modeling F.. Menter, ANSYS German GmbH 1 Unstead ANS Based Models UANS (Unstead enolds averaged Navier Stoes) Methods UANS gives unphsical single mode unstead
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationNCHRP FY 2004 Rotational Limits for Elastomeric Bearings. Final Report APPENDIX E. John F. Stanton Charles W. Roeder Peter Mackenzie-Helnwein
NCHRP -68 FY 4 Rotational Limits for Elastomeric Bearings Final Report APPENDIX E John F. Stanton Charles W. Roeder Peter Mackenzie-Helnwein Department of Civil and Environmental Engineering Universit
More informationMagnetic Field and Chemical Reaction Effects on Convective Flow of
Communications in Applied Sciences ISSN 221-7372 Volume 1, Number 1, 213, 161-187 Magnetic Field and Chemical Reaction Effects on Convective Flow of Dust Viscous Fluid P. Mohan Krishna 1, Dr.V.Sugunamma
More informationDiscontinuous Galerkin method for a class of elliptic multi-scale problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 000; 00: 6 [Version: 00/09/8 v.0] Discontinuous Galerkin method for a class of elliptic multi-scale problems Ling Yuan
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationA MODEL FOR BUOYANT PUFF DISPERSION IN URBAN AREAS
A MODEL OR BUOYANT PU DISPERSION IN URBAN AREAS Steven Herring 1, David Hall 2 and Angela Spanton 2 1 Dstl Porton Down, Salisbur, Wiltshire, SP4 0JQ (England), 2 Envirobods Ltd, 13 Badminton Close, Stevenage,
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationA RESCALING VELOCITY METHOD FOR DISSIPATIVE KINETIC EQUATIONS APPLICATIONS TO GRANULAR MEDIA
A RESCALING VELOCITY METHOD FOR DISSIPATIVE KINETIC EQUATIONS APPLICATIONS TO GRANULAR MEDIA FRANCIS FILBET AND THOMAS REY Abstract. We present a new numerical algorithm based on a relative energ scaling
More informationFilm thickness Hydrodynamic pressure Liquid saturation pressure or dissolved gases saturation pressure. dy. Mass flow rates per unit length
NOTES DERITION OF THE CLSSICL REYNOLDS EQTION FOR THIN FIL FLOWS Te lecture presents te derivation of te Renolds equation of classical lubrication teor. Consider a liquid flowing troug a tin film region
More informationVisualizing Non-Equilibrium Flow Simulations using 3-D Velocity Distribution Functions
Visualizing Non-Equilibrium Flow Simulations using 3-D Velocit Distribution Functions A. Venkattraman and A.A. Alexeenko School of Aeronautics & Astronautics, Purdue Universit, West Lafaette IN 797 Abstract.
More informationStability of plane Poiseuille Couette flows of a piezo-viscous fluid
Stabilit of plane Poiseuille Couette flows of a piezo-viscous fluid Thien Duc Tran a and Serge A. Suslov b, a Computational Engineering and Science Research Centre, Universit of Southern Queensland, Toowoomba,
More informationPHYSICAL REVIEW E 78, Lattice Boltzmann method for simulation of compressible flows on standard lattices
PHYSICAL REVIEW E 7, 016704 00 Lattice Boltzmann method for simulation of compressible flows on standard lattices Nikolaos I. Prasianakis 1, * and Ilia V. Karlin 1,, 1 Aerothermochemistr and Combustion
More informationTurbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media
Copyright 011 Tech Science Press CMES, vol.71, no., pp.149-155, 011 Turbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media Zhi Yuan Cui 1, Jiu Hui Wu 1 and Di Chen Li 1 Abstract:
More informationMagnetohydrodynamic Fluid Flow Between Two Parallel Infinite Plates Subjected To An Inclined Magnetic Field Under Pressure Gradient
Journal of Multidisciplinar Engineering Science and Technolog (JMEST) Magnetohdrodnamic Fluid Flow Between Two Parallel Infinite Plates Subjected To An Inclined Magnetic Field Under Pressure Gradient Agnes
More informationSymmetry Arguments and the Role They Play in Using Gauss Law
Smmetr Arguments and the Role The la in Using Gauss Law K. M. Westerberg (9/2005) Smmetr plas a ver important role in science in general, and phsics in particular. Arguments based on smmetr can often simplif
More informationarxiv: v1 [cond-mat.stat-mech] 26 Jan 2008
Long Time Tail of the Velocity Autocorrelation Function in a Two-Dimensional Moderately Dense Hard Disk Fluid Masaharu Isobe Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555,
More information5.3.3 The general solution for plane waves incident on a layered halfspace. The general solution to the Helmholz equation in rectangular coordinates
5.3.3 The general solution for plane waves incident on a laered halfspace The general solution to the elmhol equation in rectangular coordinates The vector propagation constant Vector relationships between
More informationAvailable online at ScienceDirect. Procedia Engineering 103 (2015 ) 52 59
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 103 (2015 ) 52 59 The 13 th Hpervelocit Impact Smposium Guided Impact Mitigation in 2D and 3D Granular Crstals Haden A. Burgone
More information520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
More informationDynamics of quantum vortices in a toroidal trap
PHYSICAL REVIEW A 79, 4362 29 Dnamics of quantum vortices in a toroidal trap Peter Mason and Natalia G. Berloff Department of Applied Mathematics and Theoretical Phsics, Universit of Cambridge, Wilberforce
More informationSPECULAR REFLECTION BY CONTACTING CRACK FACES. J.D. Achenbach and A.N. Norris. The Technological Institute Northwestern University Evanston, IL 60201
SPECULAR REFLECTION BY CONTACTING CRACK FACES J.D. Achenbach and A.N. Norris The Technological Institute Northwestern Universit Evanston, IL 60201 ABSTRACT Non-linear relations are postulated between the
More informationRheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations
Rheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations José María Montanero Departamento de Electrónica e Ingeniería
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Dec 2002
arxiv:cond-mat/0212257v1 [cond-mat.stat-mech] 11 Dec 2002 International Journal of Modern Physics B c World Scientific Publishing Company VISCOUS FINGERING IN MISCIBLE, IMMISCIBLE AND REACTIVE FLUIDS PATRICK
More informationVelocity Distribution in a Viscous Granular Gas
Velocity Distribution in a Viscous Granular Gas Alexandre Rosas,, Daniel ben-avraham 3, and Katja Lindenberg,4 Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA
More informationENGI 9420 Engineering Analysis Solutions to Additional Exercises
ENGI 940 Engineering Analsis Solutions to Additional Exercises 0 Fall [Partial differential equations; Chapter 8] The function ux (, ) satisfies u u u + = 0, subject to the x x u x,0 = u x, =. Classif
More informationDivergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF)
Divergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF) arxiv:1608.07869, arxiv:1711.01657, arxiv:1709.06644 Frankfurt University 1.February.2018 Preview
More informationOcean Dynamics. The Equations of Motion 8/27/10. Physical Oceanography, MSCI 3001 Oceanographic Processes, MSCI dt = fv. dt = fu.
Phsical Oceanograph, MSCI 3001 Oceanographic Processes, MSCI 5004 Dr. Katrin Meissner k.meissner@unsw.e.au Ocean Dnamics The Equations of Motion d u dt = 1 ρ Σ F Horizontal Equations: Acceleration = Pressure
More informationPan Pearl River Delta Physics Olympiad 2005
1 Jan. 29, 25 Morning Session (9 am 12 pm) Q1 (5 Two identical worms of length L are ling on a smooth and horizontal surface. The mass of the worms is evenl distributed along their bod length. The starting
More informationGALLOPING OF SMALL ASPECT RATIO SQUARE CYLINDER
OL. NO. JANUARY 5 ISSN 89-668 6-5 Asian Research Publishing Network (ARPN). All rights reserved. GALLOPING OF SMALL ASPET RATIO SQUARE YLINDER A. N. Rabinin and. D. Lusin Facult of Mathematics and Mechanics
More informationControl Schemes to Reduce Risk of Extinction in the Lotka-Volterra Predator-Prey Model
Journal of Applied Mathematics and Phsics, 2014, 2, 644-652 Published Online June 2014 in SciRes. http://www.scirp.org/journal/jamp http://d.doi.org/10.4236/jamp.2014.27071 Control Schemes to Reduce Risk
More informationGranular Micro-Structure and Avalanche Precursors
Granular Micro-Structure and Avalanche Precursors L. Staron, F. Radjai & J.-P. Vilotte Department of Applied Mathematics and Theoretical Physics, Cambridge CB3 0WA, UK. Laboratoire de Mécanique et Génie
More informationarxiv:cond-mat/ v2 [cond-mat.soft] 4 Aug 1998
arxiv:cond-mat/98646v2 [cond-mat.soft] 4 Aug 1998 Transition to Centrifugal Particle Motion in Rotating Drums Gerald H. Ristow Fachbereich Physik, Philipps Universität Renthof 6, D 3532 Marburg, Germany
More informationIntroduction to Differential Equations. National Chiao Tung University Chun-Jen Tsai 9/14/2011
Introduction to Differential Equations National Chiao Tung Universit Chun-Jen Tsai 9/14/011 Differential Equations Definition: An equation containing the derivatives of one or more dependent variables,
More informationON THE CHARACTERISTICS OF SOFT SOIL INELASTIC DESIGN RESPONSE SPECTRAS FOR THE CITY OF GUAYAQUIL, ECUADOR
ON THE CHAACTEISTICS OF SOFT SOIL INELASTIC DESIGN ESPONSE SPECTAS FO THE CITY OF GUAYAQUIL, ECUADO 837 Otton LAA 1 And moncao MACELO SUMMAY Five accelerograph records obtained in the cit of Guaaquil,
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More informationEP225 Note No. 4 Wave Motion
EP5 Note No. 4 Wave Motion 4. Sinusoidal Waves, Wave Number Waves propagate in space in contrast to oscillations which are con ned in limited regions. In describing wave motion, spatial coordinates enter
More informationDry and wet granular flows. Diego Berzi
Dry and wet granular flows Diego Berzi Outline 2 What? Why? How? When? Who? Where? Then? What? Granular flows many solid moving particles 3 particle mass is large (at least 10 20 molecular masses). Hence,
More informationKinetic Theory for Binary Granular Mixtures at Low Density
10 Kinetic Theory for Binary Granular Mixtures at Low Density V. Garzó Departamento de Física, Universidad de Extremadura, E-06071 Badajoz, Spain vicenteg@unex.es Many features of granular media can be
More informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at.sciencedirect.com ScienceDirect Procedia Engineering 9 (14 383 388 1th International Conference on Mechanical Engineering, ICME 13 Effects of volumetric heat source and temperature
More informationMAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation
The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural
More informationSHORT PLANE BEARINGS LUBRICATION APPLIED ON SILENT CHAIN JOINTS
Bulletin of the Transilvania Universit of Braşov Vol. 9 (58) No. - Special Issue 016 Series I: Engineering Sciences SHORT PLANE BEARINGS LUBRICATION APPLIED ON SILENT CHAIN JOINTS L. JURJ 1 R. VELICU Abstract:
More informationSubband Coding and Wavelets. National Chiao Tung University Chun-Jen Tsai 12/04/2014
Subband Coding and Wavelets National Chiao Tung Universit Chun-Jen Tsai /4/4 Concept of Subband Coding In transform coding, we use N (or N N) samples as the data transform unit Transform coefficients are
More informationAvailable online at ScienceDirect. Procedia Engineering 103 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 103 (2015 ) 237 245 The 13 th Hypervelocity Impact Symposium Strength of Granular Materials in Transient and Steady State Rapid
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More information