Precipitation In Small Systems--II. Mean Field Equations More Effective Than Population Balance

Part I (Manjunath et al., 1994, Chem. Engng Sci. 49, 1451-1463) of this paper showed that the random particle numbers and size distributions in precipitation processes in very small drops obtained by stochastic simulation techniques deviate substantially from the predictions of conventional population balance. The foregoing problem is considered in this paper in terms of a mean field approximation obtained by applying a first-order closure to an unclosed set of mean field equations presented in Part I. The mean field approximation consists of two mutually coupled partial differential equations featuring (i) the probability distribution for residual supersaturation and (ii) the mean number density of particles for each size and supersaturation from which all average properties and fluctuations can be calculated. The mean field equations have been solved by finite difference methods for (i) crystallization and (ii) precipitation of a metal hydroxide both occurring in a single drop of specified initial supersaturation. The results for the average number of particles, average residual supersaturation, the average size distribution, and fluctuations about the average values have been compared with those obtained by stochastic simulation techniques and by population balance. This comparison shows that the mean field predictions are substantially superior to those of population balance as judged by the close proximity of results from the former to those from stochastic simulations. The agreement is excellent for broad initial supersaturations at short times but deteriorates progressively at larger times. For steep initial supersaturation distributions, predictions of the mean field theory are not satisfactory thus calling for higher-order approximations. The merit of the mean field approximation over stochastic simulation lies in its potential to reduce expensive computation times involved in simulation. More effective computational techniques could not only enhance this advantage of the mean field approximation but also make it possible to use higher-order approximations eliminating the constraints under which the stochastic dynamics of the process can be predicted accurately.

Approximate solution of u'=-A2u using a relation between exp(-tA2) and exp(itA)

Let A be a positive definite operator in a Hilbert space and consider the initial value problem for u(t) = -A(2)u. Using a representation of the semigroup exp(-A(2)t) in terms of the group exp(iAt) we express u in terms of the solution of the standard heat equation w(t) = W-yy, with initial values v solving the initial value problem for v(y) = iAv. This representation is used to construct a method for approximating u in terms of approximations of v. In the case that A is a 2(nd) order elliptic operator the method is combined with finite elements in the spatial variable and then reduces the solution of the 4(th) order equation for u to that of the 2(nd) order equation for v, followed by the solution of the heat equation in one space variable.

State space formulation of ammonia reactor dynamics

The specific objective of this paper is to develop a state space model of a tubular ammonia reactor which is the heart of an ammonia plant in a fertiliser complex. A ninth order model with three control inputs and two disturbance inputs is generated from the nonlinear distributed model using linearization and lumping approximations. The lumped model is chosen such that the steady state temperature at the exit of the catalyst bed computed from the simplified state space model is close enough to the one computed from the nonlinear steady state model. The model developed in this paper is very useful for the design of continuous/discrete versions of single variable/multivariable control algorithms.

Flow through a plateau border of cellular foam

Experimentally measured average velocities through plateau borders of stationary cellular foam, when compared with those calculated with the assumption of rigid Plateau Border walls, show that the assumption of rigid walls severely underestimates the velocities. An analysis of the situation wherein plateau border walls have velocities, as decided by the surface viscosity of the system, is presented here. The plateau border is idealized as a pipe of equilateral triangular cross-section with vertices of the triangle having zero velocity. The pertinent form of Navier-Stoke's equations with inhomogeneous boundary conditions and its solution through a procedure of successive approximations is presented in dimensionless form. The solution reduces to the known solution of slow steady flow through a triangular pipe, when surface viscosity is infinite. Results indicate that the assumption of rigid plateau border walls is valid only when value of the inverse of dimensionless surface viscosity is less than 0.044. Beyond that the assumption severely underestimates the flow and the effect of nonrigidity of the wall must be considered.

Non-linear flapping vibrations of rotating blades

The non-linear equations of motion of a rotating blade undergoing extensional and flapwise bending vibration are derived, including non-linearities up to O (ε3). The strain-displacement relationship derived is compared with expressions derived by earlier investigators and the errors and the approximations made in some of those are brought out. The equations of motion are solved under the inextensionality condition to obtain the influence of the amplitude on the fundamental flapwise natural frequency of the rotating blade. It is found that large finite amplitudes have a softening effect on the flapwise frequency and that this influence becomes stronger at higher speeds of rotation.

Boundary-layer convective heat transfer in a convergent channel

In this paper the numerical solution of the heat transfer problem in a convergent channel with uniform and non-uniform wall temperatures under boundary-layer approximations has been presented. Also, a semi-analytical solution for uniform wall temperature has been obtained.

An Analysis of a Boundary Layer Diffusion Flame Over a Porous Flat Plate in a Confined Flow

A numerical analysis of the gas dynamic structure of a two-dimensional laminar boundary layer diffusion flame over a porous flat plate in a confined flow is made on the basis of the familiar boundary layer and flame sheet approximations neglecting buoyancy effects. The governing equations of aerothermochemistry with the appropriate boundary conditions are solved using the Patankar-Spalding method. The analysis predicts the flame shape, profiles of temperature, concentrations of variousspecies, and the density of the mixture across the boundary layer. In addition, it also predicts the pressure gradient in the flow direction arising from the confinement ofthe flow and the consequent velocity overshoot near the flame surface. The results of thecomputation performed for an n-pentane-air system are compared with experimental data andthe agreement is found to be satisfactory.

Continuation of limit cycles near saddle homoclinic points using splines in phase space

We present a new algorithm for continuation of limit cycles of autonomous systems as a system parameter is varied. The algorithm works in phase space with an ordered set of points on the limit cycle, along with spline interpolation. Currently popular algorithms in bifurcation analysis packages compute time-domain approximations of limit cycles using either shooting or collocation. The present approach seems useful for continuation near saddle homoclinic points, where it encounters a corner while time-domain methods essentially encounter a discontinuity (a relatively short period of rapid variation). Other phase space-based algorithms use rescaled arclength in place of time, but subsequently resemble the time-domain methods. Compared to these, we introduce additional freedom through a variable stretching of arclength based on local curvature, through the use of an auxiliary index-based variable. Several numerical examples are presented. Comparisons with results from the popular package, MATCONT, are favorable close to saddle homoclinic points.

Vibrations of a Beam in Variable Contact With a Flat Surface

We study small vibrations of cantilever beams contacting a rigid surface. We study two cases: the first is a beam that sags onto the ground due to gravity, and the second is a beam that sticks to the ground through reversible adhesion. In both cases, the noncontacting length varies dynamically. We first obtain the governing equations and boundary conditions, including a transversality condition involving an end moment, using Hamilton's principle. Rescaling the variable length to a constant value, we obtain partial differential equations with time varying coefficients, which, upon linearization, give the natural frequencies of vibration. The natural frequencies for the first case (gravity without adhesion) match that of a clamped-clamped beam of the same nominal length; frequencies for the second case, however, show no such match. We develop simple, if atypical, single degree of freedom approximations for the first modes of these two systems, which provide insights into the role of the static deflection profile, as well as the end moment condition, in determining the first natural frequencies of these systems. Finally, we consider small transverse sinusoidal forcing of the first case and find that the governing equation contains both parametric and external forcing terms. For forcing at resonance, w find that either the internal or the external forcing may dominate.

Mass transfer effects on the generalised vortex flow over a stationary surface with or without magnetic field

The effect of injection and suction on the generalised vortex flow of a steady laminar incompressible fluid over a stationary infinite disc with or without magnetic field under boundary-layer approximations has been studied. The coupled nonlinear ordinary differential equations governing the self-similar flow have been numerically solved using the finite-difference scheme. The results indicate that the injection produces a deeper inflow layer and de-stabilises the motion while suction or magnetic field suppresses the inflow layer and produces stability. The effect of decreasingn, the parameter characterising the nature of vortex flow, is similar to that of increasing the injection rate.

State Dependent Attempt Rate Modeling of Single Cell IEEE 802.11 WLANs with Homogeneous Nodes and Poisson Arrivals

Analytical models of IEEE 802.11-based WLANs are invariably based on approximations, such as the well-known mean-field approximations proposed by Bianchi for saturated nodes. In this paper, we provide a new approach for modeling the situation when the nodes are not saturated. We study a State Dependent Attempt Rate (SDAR) approximation to model M queues (one queue per node) served by the CSMA/CA protocol as standardized in the IEEE 802.11 DCF. The approximation is that, when n of the M queues are non-empty, the attempt probability of the n non-empty nodes is given by the long-term attempt probability of n saturated nodes as provided by Bianchi's model. This yields a coupled queue system. When packets arrive to the M queues according to independent Poisson processes, we provide an exact model for the coupled queue system with SDAR service. The main contribution of this paper is to provide an analysis of the coupled queue process by studying a lower dimensional process and by introducing a certain conditional independence approximation. We show that the numerical results obtained from our finite buffer analysis are in excellent agreement with the corresponding results obtained from ns-2 simulations. We replace the CSMA/CA protocol as implemented in the ns-2 simulator with the SDAR service model to show that the SDAR approximation provides an accurate model for the CSMA/CA protocol. We also report the simulation speed-ups thus obtained by our model-based simulation.

Representation and properties of para-Bose oscillator operators. II. Coherent states and the minimum uncertainty states

The energy, position, and momentum eigenstates of a para-Bose oscillator system were considered in paper I. Here we consider the Bargmann or the analytic function description of the para-Bose system. This brings in, in a natural way, the coherent states ||z;alpha> defined as the eigenstates of the annihilation operator ?. The transformation functions relating this description to the energy, position, and momentum eigenstates are explicitly obtained. Possible resolution of the identity operator using coherent states is examined. A particular resolution contains two integrals, one containing the diagonal basis ||z;alpha><−z;alpha||. We briefly consider the normal and antinormal ordering of the operators and their diagonal and discrete diagonal coherent state approximations. The problem of constructing states with a minimum value of the product of the position and momentum uncertainties and the possible alpha dependence of this minimum value is considered. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Non-linear flapping vibrations of rotating blades

The non-linear equations of motion of a rotating blade undergoing extensional and flapwise bending vibration are derived, including non-linearities up to O (ε3). The strain-displacement relationship derived is compared with expressions derived by earlier investigators and the errors and the approximations made in some of those are brought out. The equations of motion are solved under the inextensionality condition to obtain the influence of the amplitude on the fundamental flapwise natural frequency of the rotating blade. It is found that large finite amplitudes have a softening effect on the flapwise frequency and that this influence becomes stronger at higher speeds of rotation.

Some direct numerical approaches to the solution of the singular nonlinear transonic integral equation

The nonlinear singular integral equation of transonic flow is examined, noting that standard numerical techniques are not applicable in solving it. The difficulties in approximating the integral term in this expression were solved by special methods mitigating the inaccuracies caused by standard approximations. It was shown how the infinite domain of integration can be reduced to a finite one; numerical results were plotted demonstrating that the methods proposed here improve accuracy and computational economy.

Structure of shock waves at re-entry speeds

The unified structure of steady, one-dimensional shock waves in argon, in the absence of an external electric or magnetic field, is investigated. The analysis is based on a two-temperature, three-fluid continuum approach, using the Navier—Stokes equations as a model and including non-equilibrium collisional as well as radiative ionization phenomena. Quasi charge neutrality and zero velocity slip are assumed. The integral nature of the radiative terms is reduced to analytical forms through suitable spectral and directional approximations. The analysis is based on the method of matched asymptotic expansions. With respect to a suitably chosen small parameter, which is the ratio of atom-atom elastic collisional mean free-path to photon mean free-path, the following shock morphology emerges: within the radiation and electron thermal conduction dominated outer layer occurs an optically transparent discontinuity which consists of a chemically frozen heavy particle (atoms and ions) shock and a collisional ionization relaxation layer. Solutions are obtained for the first order with respect to the small parameter of the problem for two cases: (i) including electron thermal conduction and (ii) neglecting it in the analysis of the outer layer. It has been found that the influence of electron thermal conduction on the shock structure is substantial. Results for various free-stream conditions are presented in the form of tables and figures.