Korteweg-de Vries equation for the propagation of fast ion-acoustic waves in multi-dimensions

Abstract is not available.

Heat transfer to power-law fluids in thermal entrance region with viscous dissipation for constant heat flux conditions

An analytical solution of the heat transfer problem with viscous dissipation for non-Newtonian fluids with power-law model in the thermal entrance region of a circular pipe and two parallel plates under constant heat flux conditions is obtained using eigenvalue approach by suitably replacing one of the boundary conditions by total energy balance equation. Analytical expressions for the wall and the bulk temperatures and the local Nusselt number are presented. The results are in close agreement with those obtained by implicit finite-difference scheme. It is found that the role of viscous dissipation on heat transfer is completely different for heating and cooling conditions at the wall. The results for the case of cooling at the wall are of interest in the design of the oil pipe line.

Simulation of laminar flow in a three-dimensional lid-driven cavity by lattice Boltzmann method

Purpose - The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid-driven flow in a three-dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two-dimensional (2D) cavity. Design/methodology/approach - The second-order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25-4) and Reynolds numbers (0.01-1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier-Stokes (NS) equations and with available experimental data. Findings - The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side-walls in transverse direction (in 3D), due to which the vorticity in the core-region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry-plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid-motion are observed; these are weak at low. Reynolds numbers, but become quite strong at high Reynolds numbers. Originality/value - The findings will be useful in the study of variety of enclosed fluid flows.

The alpha method for solving differential algebraic inequality (DAI) systems

This paper describes an algorithm for ``direct numerical integration'' of the initial value Differential-Algebraic Inequalities (DAI) in a time stepping fashion using a sequential quadratic programming (SQP) method solver for detecting and satisfying active path constraints at each time step. The activation of a path constraint generally increases the condition number of the active discretized differential algebraic equation's (DAE) Jacobian and this difficulty is addressed by a regularization property of the alpha method. The algorithm is locally stable when index 1 and index 2 active path constraints and bounds are active. Subject to available regularization it is seen to be stable for active index 3 active path constraints in the numerical examples. For the high index active path constraints, the algorithm uses a user-selectable parameter to perturb the smaller singular values of the Jacobian with a view to reducing the condition number so that the simulation can proceed. The algorithm can be used as a relatively cheaper estimation tool for trajectory and control planning and in the context of model predictive control solutions. It can also be used to generate initial guess values of optimization variables used as input to inequality path constrained dynamic optimization problems. The method is illustrated with examples from space vehicle trajectory and robot path planning.

Understanding solvent effects on structure and reactivity of organic intermediates: a Raman study

The effect of solvent on chemical reactivity has generally been explained on the basis of the dielectric constant and viscosity. However a number of spectroscopic studies, including UV-VIS, IR and Raman, has led to numerous empirical parameters to define solvent effect based on either solvating ability or polarity scale. These parameters include solvent polarizability, dipolarity, Lewis acidity and Lewis basicity, E-T(30), pi*, alpha, beta etc. However, from a structural point of view, we can separate solvation as static and dynamic processes. The static solvation basically relates to stabilization of the molecular structure by the solvent to attain the equilibrium structure, both in the intermediate and ground state. Dynamic solvation relates to solvent reorganization-induced dynamics prior to the structural reorganization to reach the equilibrium state. In this paper, we present (a) structural distortions induced by the solvent due to preferential solvation of the triplet excited state, and (b) the importance of dynamic solvation induced by vibronic coupling (pseudo-Jahn-Teller coupling). The examples include the effect of solvent on structure and reactivity of excited states of 2,2,2-trifluoroacetophenone (TFA). Based on the comparison of time resolved resonance Raman (TR3) data of TFA and other substituted acetophenone systems, it was found that change in solvent polarity indeed results in electronic state switching and structural changes in the excited state, which explains the trend in reactivity. Further, a TR3 study of fluoranil (FA) in the triplet excited state in solvents of varying polarities indicates that the structure of FA in the triplet excited state is determined by vibronic coupling effects and thus distorted structure. These experimental results have been well supported by density functional theoretical computational studies.

Metamodeling approach to predict friction factor of alluvial channel

The flow resistance of an alluvial channel flow is not only affected by the Reynolds number and the roughness conditions but also the Froude number. Froude number is the most basic parameter in the case of the alluvial channel, thus effect of Froude number on resistance to flow should be considered in the formulation of the friction factor, which is not in the case of present available resistance equations. At present, no generally acceptable quantitative description of the effects of the Froude number on hydraulic resistance has been developed. Metamodeling technique, which is particularly useful in modeling a complex processes or where knowledge of the physics is limited, is presented as a tool complimentary to modeling friction factor in alluvial channels. Present work uses, a radial basis metamodel, which is a type of neural network modeling, to find the effect of Froude number on the flow resistance. Based on the experimental data taken from different sources, it has been found that the predicting capability of the present model is on acceptable level. Present work also tries in formulating an empirical equation for resistance in alluvial channel comprising all the three majorm, parameters, namely, roughness parameter, Froude number and Reynolds number. (C) 2009 Elsevier B.V. All rights reserved.

A Model of Anomalous Chain Translocation Dynamics

A model of polymer translocation based on the stochastic dynamics of the number of monomers on one side of a pore-containing surface is formulated in terms of a one-dimensional generalized Langevin equation, in which the random force is assumed to be characterized by long-ranged temporal correlations. The model is introduced to rationalize anomalies in measured and simulated values of the average time of passage through the pore, which in general cannot be satisfactorily accounted for by simple Brownian diffusion mechanisms. Calculations are presented of the mean first passage time for barrier crossing and of the mean square displacement of a monomeric segment, in the limits of strong and weak diffusive bias. The calculations produce estimates of the exponents in various scaling relations that are in satisfactory agreement with available data.

A note on the Hadwiger number of circular arc graphs

The intention of this note is to motivate the researchers to study Hadwiger's conjecture for circular arc graphs. Let η(G) denote the largest clique minor of a graph G, and let χ(G) denote its chromatic number. Hadwiger's conjecture states that η(G)greater-or-equal, slantedχ(G) and is one of the most important and difficult open problems in graph theory. From the point of view of researchers who are sceptical of the validity of the conjecture, it is interesting to study the conjecture for graph classes where η(G) is guaranteed not to grow too fast with respect to χ(G), since such classes of graphs are indeed a reasonable place to look for possible counterexamples. We show that in any circular arc graph G, η(G)less-than-or-equals, slant2χ(G)−1, and there is a family with equality. So, it makes sense to study Hadwiger's conjecture for this family.

A new framework for solution of multidimensional population balance equations

A new framework is proposed in this work to solve multidimensional population balance equations (PBEs) using the method of discretization. A continuous PBE is considered as a statement of evolution of one evolving property of particles and conservation of their n internal attributes. Discretization must therefore preserve n + I properties of particles. Continuously distributed population is represented on discrete fixed pivots as in the fixed pivot technique of Kumar and Ramkrishna [1996a. On the solution of population balance equation by discretization-I A fixed pivot technique. Chemical Engineering Science 51(8), 1311-1332] for 1-d PBEs, but instead of the earlier extensions of this technique proposed in the literature which preserve 2(n) properties of non-pivot particles, the new framework requires n + I properties to be preserved. This opens up the use of triangular and tetrahedral elements to solve 2-d and 3-d PBEs, instead of the rectangles and cuboids that are suggested in the literature. Capabilities of computational fluid dynamics and other packages available for generating complex meshes can also be harnessed. The numerical results obtained indeed show the effectiveness of the new framework. It also brings out the hitherto unknown role of directionality of the grid in controlling the accuracy of the numerical solution of multidimensional PBEs. The numerical results obtained show that the quality of the numerical solution can be improved significantly just by altering the directionality of the grid, which does not require any increase in the number of points, or any refinement of the grid, or even redistribution of pivots in space. Directionality of a grid can be altered simply by regrouping of pivots.

Measurement of shock stand-off distance on a 120 degrees blunt cone model at hypersonic Mach number in Argon

The flow around a 120 degrees blunt cone model with a base radius of 60mm has been visualised at Mach 14.8 and 9.1 using argon as the test gas, at the newly established high speed schlieren facility in the IISc hypersonic shock tunnel HST2. The experimental shock stand off distance around the blunt cone is compared with that obtained using a commercial CFD package. The computed values of shock stand off distance of the blunt cone is found to agree reasonably well with the experimental data.

Constitutive equation for elevated temperature flow behavior of plain carbon steels using dimensional analysis

Dimensional analysis using π-theorem is applied to the variables associated with plastic deformation. The dimensionless groups thus obtained are then related and rewritten to obtain the constitutive equation. The constants in the constitutive equation are obtained using published flow stress data for carbon steels. The validity of the constitutive equation is tested for steels with up to 1.54 wt%C at temperatures: 850–1200 °C and strain rates: 6 × 10−6–2 × 10−2 s−1. The calculated flow stress agrees favorably with experimental data.

Dynamics of Barrierless and Activated Chemical Reactions in a Dispersive Medium within the Fractional Diffusion Equation Approach

Barrierless chemical reactions have often been modeled as a Brownian motion on a one-dimensional harmonic potential energy surface with a position-dependent reaction sink or window located near the minimum of the surface. This simple (but highly successful) description leads to a nonexponential survival probability only at small to intermediate times but exponential decay in the long-time limit. However, in several reactive events involving proteins and glasses, the reactions are found to exhibit a strongly nonexponential (power law) decay kinetics even in the long time. In order to address such reactions, here, we introduce a model of barrierless chemical reaction where the motion along the reaction coordinate sustains dispersive diffusion. A complete analytical solution of the model can be obtained only in the frequency domain, but an asymptotic solution is obtained in the limit of long time. In this case, the asymptotic long-time decay of the survival probability is a power law of the Mittag−Leffler functional form. When the barrier height is increased, the decay of the survival probability still remains nonexponential, in contrast to the ordinary Brownian motion case where the rate is given by the Smoluchowski limit of the well-known Kramers' expression. Interestingly, the reaction under dispersive diffusion is shown to exhibit strong dependence on the initial state of the system, thus predicting a strong dependence on the excitation wavelength for photoisomerization reactions in a dispersive medium. The theory also predicts a fractional viscosity dependence of the rate, which is often observed in the reactions occurring in complex environments.

Long-lived giant number fluctuations in a swarming granular nematic

Coherently moving flocks of birds, beasts, or bacteria are examples of living matter with spontaneous orientational order. How do these systems differ from thermal equilibrium systems with such liquid crystalline order? Working with a fluidized monolayer of macroscopic rods in the nematic liquid crystalline phase, we find giant number fluctuations consistent with a standard deviation growing linearly with the mean, in contrast to any situation where the central limit theorem applies. These fluctuations are long-lived, decaying only as a logarithmic function of time. This shows that flocking, coherent motion, and large-scale inhomogeneity can appear in a system in which particles do not communicate except by contact.

Hadwiger Number and the Cartesian Product of Graphs

The Hadwiger number eta(G) of a graph G is the largest integer n for which the complete graph K-n on n vertices is a minor of G. Hadwiger conjectured that for every graph G, eta(G) >= chi(G), where chi(G) is the chromatic number of G. In this paper, we study the Hadwiger number of the Cartesian product G square H of graphs. As the main result of this paper, we prove that eta(G(1) square G(2)) >= h root 1 (1 - o(1)) for any two graphs G(1) and G(2) with eta(G(1)) = h and eta(G(2)) = l. We show that the above lower bound is asymptotically best possible when h >= l. This asymptotically settles a question of Z. Miller (1978). As consequences of our main result, we show the following: 1. Let G be a connected graph. Let G = G(1) square G(2) square ... square G(k) be the ( unique) prime factorization of G. Then G satisfies Hadwiger's conjecture if k >= 2 log log chi(G) + c', where c' is a constant. This improves the 2 log chi(G) + 3 bound in [2] 2. Let G(1) and G(2) be two graphs such that chi(G1) >= chi(G2) >= clog(1.5)(chi(G(1))), where c is a constant. Then G1 square G2 satisfies Hadwiger's conjecture. 3. Hadwiger's conjecture is true for G(d) (Cartesian product of G taken d times) for every graph G and every d >= 2. This settles a question by Chandran and Sivadasan [2]. ( They had shown that the Hadiwger's conjecture is true for G(d) if d >= 3).

Design of plane sand-bed channels affected by seepage

The importance of seepage in the design of channels is discussed. Experimental investigations reveal that seepage, either in the downward direction (suction) or in the upward direction (injection), can significantly change the resistance as well as the mobility of the sand-bed particles. A resistance equation relating 'particle Reynolds number' and 'shear Reynolds number' under seepage conditions is developed for plane sediment beds. Finally, a detailed design procedure of the plane sediment beds affected by seepage is presented.