Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 1. DNS and fluctuating force model

The fluctuating force model is developed and applied to the turbulent flow of a gas-particle suspension in a channel in the limit of high Stokes number, where the particle relaxation time is large compared to the fluid correlation time, and low particle Reynolds number where the Stokes drag law can be used to describe the interaction between the particles and fluid. In contrast to the Couette flow, the fluid velocity variances in the different directions in the channel are highly non-homogeneous, and they exhibit significant variation across the channel. First, we analyse the fluctuating particle velocity and acceleration distributions at different locations across the channel. The distributions are found to be non-Gaussian near the centre of the channel, and they exhibit significant skewness and flatness. However, acceleration distributions are closer to Gaussian at locations away from the channel centre, especially in regions where the variances of the fluid velocity fluctuations are at a maximum. The time correlations for the fluid velocity fluctuations and particle acceleration fluctuations are evaluated, and it is found that the time correlation of the particle acceleration fluctuations is close to the time correlations of the fluid velocity in a `moving Eulerian' reference, moving with the mean fluid velocity. The variances of the fluctuating force distributions in the Langevin simulations are determined from the time correlations of the fluid velocity fluctuations and the results are compared with direct numerical simulations. Quantitative agreement between the two simulations are obtained provided the particle viscous relaxation time is at least five times larger than the fluid integral time.

Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 2. Comparison of fluctuating force simulations and experiments

The particle and fluid velocity fluctuations in a turbulent gas-particle suspension are studied experimentally using two-dimensional particle image velocimetry with the objective of comparing the experiments with the predictions of fluctuating force simulations. Since the fluctuating force simulations employ force distributions which do not incorporate the modification of fluid turbulence due to the particles, it is of importance to quantify the turbulence modification in the experiments. For experiments carried out at a low volume fraction of 9.15 x 10(-5) (mass loading is 0.19), where the viscous relaxation time is small compared with the time between collisions, it is found that the gas-phase turbulence is not significantly modified by the presence of particles. Owing to this, quantitative agreement is obtained between the results of experiments and fluctuating force simulations for the mean velocity and the root mean square of the fluctuating velocity, provided that the polydispersity in the particle size is incorporated in the simulations. This is because the polydispersity results in a variation in the terminal velocity of the particles which could induce collisions and generate fluctuations; this mechanism is absent if all of the particles are of equal size. It is found that there is some variation in the particle mean velocity very close to the wall depending on the wall-collision model used in the simulations, and agreement with experiments is obtained only when the tangential wall-particle coefficient of restitution is 0.7. The mean particle velocity is in quantitative agreement for locations more than 10 wall units from the wall of the channel. However, there are systematic differences between the simulations and theory for the particle concentrations, possibly due to inadequate control over the particle feeding at the entrance. The particle velocity distributions are compared both at the centre of the channel and near the wall, and the shape of the distribution function near the wall obtained in experiments is accurately predicted by the simulations. At the centre, there is some discrepancy between simulations and experiment for the distribution of the fluctuating velocity in the flow direction, where the simulations predict a bi-modal distribution whereas only a single maximum is observed in the experiments, although both distributions are skewed towards negative fluctuating velocities. At a much higher particle mass loading of 1.7, where the time between collisions is smaller than the viscous relaxation time, there is a significant increase in the turbulent velocity fluctuations by similar to 1-2 orders of magnitude. Therefore, it becomes necessary to incorporate the modified fluid-phase intensity in the fluctuating force simulation; with this modification, the mean and mean-square fluctuating velocities are within 20-30% of the experimental values.

Group Actions and Elementary Number Theory

Abstract | In this article the shuffling of cards is studied by using the concept of a group action. We use some fundamental results in Elementary Number Theory to obtain formulas for the orders of some special shufflings, namely the Faro and Monge shufflings and give necessary and sufficient conditions for the Monge shuffling to be a cycle. In the final section we extend the considerations to the shuffling of multisets.

Synthesis, Reactivity Studies, Structural Aspects, and Solution Behavior of Half Sandwich Ruthenium(II) N,N `,N `'-Triarylguanidinate Complexes

[(eta(6)-C(10)H(14))RuCl(mu-Cl)](2) (eta(6)-C(10)H(14) = eta(6)-p-cymene) was subjected to a bridge-splitting reaction with N,N',N `'-triarylguanidines, (ArNH)(2)C=NAr, in toluene at ambient temperature to afford [(eta(6)-C(10)H(14))RuCl{kappa(2)(N,N')((ArN)(2)C-N(H)Ar)}] (Ar = C(6)H(4)Me-4 (1), C(6)H(4)(OMe)-2 (2), C(6)H(4)Me-2 (3), and C(6)H(3)Me(2)-2,4 (4)) in high yield with a view aimed at understanding the influence of substituent(s) on the aryl rings of the guanidine upon the solid-state structure, solution behavior, and reactivity pattern of the products. Complexes 1-3 upon reaction with NaN(3) in ethanol at ambient temperature afforded [(eta(6)-C(10)H(14))RuN(3){kappa(2)(N,N')((ArN)(2)C-N(H)Ar)}] (Ar = C(6)H(4)Me-4 (5), C(6)H(4)(OMe)-2 (6), and C(6)H(4)Me-2 (7)) in high yield. [3 + 2] cycloaddition reaction of 5-7 with RO(O)C-C C-C(O)OR (R = Et (DEAD) and Me (DMAD)) (diethylacetylenedicarboxylate, DEAD; dimethylacetylenedicarboxylate, DMAD) in CH(2)Cl(2) at ambient temperature afforded [(eta(6)-C(10)H(14))Ru{N(3)C(2)(C(O)OR)(2)}{kappa(2)(N,N')((ArN)(2) C-N(H)Ar)}center dot xH(2)O (x = 1, R = Et, Ar = C(6)H(4)Me-4 (8 center dot H(2)O); x = 0, R = Me, Ar = C(6)H(4)(OMe)-2 (9), and C(6)H(4)Me-2 (10)) in moderate yield. The molecular structures of 1-6, 8 center dot H(2)O, and 10 were determined by single crystal X-ray diffraction data. The ruthenium atom in the aforementioned complexes revealed pseudo octahedral ``three legged piano stool'' geometry. The guanidinate ligand in 2, 3, and 6 revealed syn-syn conformation and that in 4, and 10 revealed syn-anti conformation, and the conformational difference was rationalized on the basis of subtle differences in the stereochemistry of the coordinated nitrogen atoms caused by the aryl moiety in 3 and 4 or steric overload caused by the substituents around the ruthenium atom in 10. The bonding pattern of the CN(3) unit of the guanidinate ligand in the new complexes was explained by invoking n-pi conjugation involving the interaction of the NHAr/N(coord)Ar lone pair with C=N pi* orbital of the imine unit. Complexes 1, 2, 5, 6, 8 center dot H(2)O, and 9 were shown to exist as a single isomer in solution as revealed by NMR data, and this was ascribed to a fast C-N(H)Ar bond rotation caused by a less bulky aryl moiety in these complexes. In contrast, 3 and 10 were shown to exist as a mixture of three and five isomers in about 1:1:1 and 1.0:1.2:2:7:3.5:6.9 ratios, respectively in solution as revealed by a VT (1)H NMR, (1)H-(1)H COSY in conjunction with DEPT-90 (13)C NMR data measured at 233 K in the case of 3. The multiple number of isomers in solution was ascribed to the restricted C-N(H)(o-tolyl) bond rotation caused by the bulky o-tolyl substituent in 3 or the aforementioned restricted C-NH(o-tolyl) bond rotation as well as the restricted ruthenium-arene(centroid) bond rotation caused by the substituents around the ruthenium atom in 10.

Mathematical Model of Viral Kinetics In Vitro Estimates the Number of E2-CD81 Complexes Necessary for Hepatitis C Virus Entry

Interaction between the hepatitis C virus (HCV) envelope protein E2 and the host receptor CD81 is essential for HCV entry into target cells. The number of E2-CD81 complexes necessary for HCV entry has remained difficult to estimate experimentally. Using the recently developed cell culture systems that allow persistent HCV infection in vitro, the dependence of HCV entry and kinetics on CD81 expression has been measured. We reasoned that analysis of the latter experiments using a mathematical model of viral kinetics may yield estimates of the number of E2-CD81 complexes necessary for HCV entry. Here, we constructed a mathematical model of HCV viral kinetics in vitro, in which we accounted explicitly for the dependence of HCV entry on CD81 expression. Model predictions of viral kinetics are in quantitative agreement with experimental observations. Specifically, our model predicts triphasic viral kinetics in vitro, where the first phase is characterized by cell proliferation, the second by the infection of susceptible cells and the third by the growth of cells refractory to infection. By fitting model predictions to the above data, we were able to estimate the threshold number of E2-CD81 complexes necessary for HCV entry into human hepatoma-derived cells. We found that depending on the E2-CD81 binding affinity, between 1 and 13 E2-CD81 complexes are necessary for HCV entry. With this estimate, our model captured data from independent experiments that employed different HCV clones and cells with distinct CD81 expression levels, indicating that the estimate is robust. Our study thus quantifies the molecular requirements of HCV entry and suggests guidelines for intervention strategies that target the E2-CD81 interaction. Further, our model presents a framework for quantitative analyses of cell culture studies now extensively employed to investigate HCV infection.

Rotation field in wedge indentation of metals

A study is made of the rotation field in wedge indentation of metals using copper as the model material system. Wedges with apical angles of 60 and 120 are used to indent annealed copper, and the deformation is mapped using image correlation. The indentation of annealed and strain-hardened copper is simulated using finite element analysis. The rotation field, derived from the deformation measurements, provides a clear way of distinguishing between cutting and compressive modes of deformation. Largely unidirectional rotation on one side of the symmetry line with small spatial rotation gradients is characteristic of compression. Bidirectional rotation with neighboring regions of opposing rotations and locally high rotation gradients characterizes cutting. In addition, the rotation demarcates such characteristic regions as the pile-up zone in indentation of a strain-hardened metal. The residual rotation field obtained after unloading is essentially the same as that at full load, indicating that it is a scalar proxy for plastic deformation as a whole.

Rainbow Connection Number and Connectivity

The rainbow connection number, rc(G), of a connected graph G is the minimum number of colors needed to color its edges, so that every pair of vertices is connected by at least one path in which no two edges are colored the same. Our main result is that rc(G) <= inverted right perpendicularn/2inverted left perpendicular for any 2-connected graph with at least three vertices. We conjecture that rc(G) <= n/kappa + C for a kappa-connected graph G of order n, where C is a constant, and prove the conjecture for certain classes of graphs. We also prove that rc(G) < (2 + epsilon)n/kappa + 23/epsilon(2) for any epsilon > 0.

Solar dynamo models with realistic internal rotation

Solar dynamo models based on differential rotation inferred from helioseismology tend to produce rather strong magnetic activity at high solar latitudes, in contrast to the observed fact that sunspots appear at low latitudes. We show that a meridional circulation penetrating below the tachocline can solve this problem.

Real-time monitoring of critical nodes with minimal number of Phasor Measurement Units

This paper presents a method for placement of Phasor Measurement Units, ensuring the monitoring of vulnerable buses which are obtained based on transient stability analysis of the overall system. Real-time monitoring of phase angles across different nodes, which indicates the proximity to instability, the very purpose will be well defined if the PMUs are placed at buses which are more vulnerable. The issue is to identify the key buses where the PMUs should be placed when the transient stability prediction is taken into account considering various disturbances. Integer Linear Programming technique with equality and inequality constraints is used to find out the optimal placement set with key buses identified from transient stability analysis. Results on IEEE-14 bus system are presented to illustrate the proposed approach.