Modeling of precipitation in reverse micellar systems

A model of the precipitation process in reverse micelles has been developed to calculate the size of fine particles obtained therein. While the method shares several features of particle nucleation and growth common to precipitation in large systems, complexities arise in describing the processes of nucleation, due to the extremely small size of a micelle and of particle growth caused by fusion among the micelles. Occupancy of micelles by solubilized molecules is governed by Poisson statistics, implying most of them are empty and cannot nucleate of its own. The model therefore specifies the minimum number of solubilized molecules required to form a nucleus which is used to calculate the homogeneous nucleation rate. Simultaneously, interaction between micelles is assumed to occur by Brownian collision and instantaneous fusion. Analysis of time scales of various events shows growth of particles to be very fast compared to other phenomena occurring. This implies that nonempty micelles either are supersaturated or contain a single precipitated particle and allows application of deterministic population balance equations to describe the evolution of the system with time. The model successfully predicts the experimental measurements of Kandori ct al.(3) on the size of precipitated CaCO3 particles, obtained by carbonation of reverse micelles containing aqueous Ca(OH)(2) solution.

Physics based basis function for vibration analysis of high speed rotating beams

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

Automatic decomposition of unstructured meshes employing genetic algorithms for parallel FEM computations

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

Structural vibration control using resistively shunted piezoceramics

Application of piezoceramic materials in actuation and sensing of vibration is of current interest. Potential and more popular applications of piezoceramics are probably in the field of active vibration control. However, the objective of this work is to investigate the effect of shunted piezoceramics as passive vibration control devices when bonded to a host structure. Resistive shunting of a piezoceramic bonded to a cantilevered duralumin beam has been investigated. The piezoceramic is connected in parallel to an electrical network comprising of resistors and inductors. The piezoceramic is a capacitor that stores and discharges electrical energy that is transformed from the mechanical motion of the structure to which it is bonded. A resistor across the piezoceramic would be termed as a resistively shunted piezoceramic. Similarly, an inductor across the piezoceramic is termed as a resonantly shunted piezoceramic. In this study, the effect of resistive shunting on the nature of damping enhancement to the host structure has been investigated. Analytical studies are presented along with experimental results.

An improved Monte Carlo scheme for simulation of synthesis of nanoparticles in reverse micelles

An improved Monte Carlo technique is presented in this work to simulate nanoparticle formation through a micellar route. The technique builds on the simulation technique proposed by Bandyopadhyaya et al. (Langmuir 2000, 16, 7139) which is general and rigorous but at the same time very computation intensive, so much so that nanoparticle formation in low occupancy systems cannot be simulated in reasonable time. In view of this, several strategies, rationalized by simple mathematical analyses, are proposed to accelerate Monte Carlo simulations. These are elimination of infructuous events, removal of excess reactant postreaction, and use of smaller micelle population a large number of times. Infructuous events include collision of an empty micelle with another empty one or with another one containing only one molecule or only a solid particle. These strategies are incorporated in a new simulation technique which divides the entire micelle population in four classes and shifts micelles from one class to other as the simulation proceeds. The simulation results, throughly tested using chi-square and other tests, show that the predictions of the improved technique remain unchanged, but with more than an order of magnitude decrease in computational effort for some of the simulations reported in the literature. A post priori validation scheme for the correctness of the simulation results has been utilized to propose a new simulation strategy to arrive at converged simulation results with near minimum computational effort.

Closed form solutions for element equilibrium and flexibility matrices of eight node rectangular plate bending element using integrated force method

Closed form solutions for equilibrium and flexibility matrices of the Mindlin-Reissner theory based eight-node rectangular plate bending element (MRP8) using integrated Force Method (IFM) are presented in this paper. Though these closed form solutions of equilibrium and flexibility matrices are applicable to plate bending problems with square/rectangular boundaries, they reduce the computational time significantly and give more exact solutions. Presented closed form solutions are validated by solving large number of standard square/rectangular plate bending benchmark problems for deflections and moments and the results are compared with those of similar displacement-based eight-node quadrilateral plate bending elements available in the literature. The results are also compared with the exact solutions.

Filling Characteristics for an Activated Carbon Based Adsorbed Natural Gas Storage System

The storage capacity of an activated carbon bed is studied using a 2D transport model with constant inlet flow conditions. The predicted filling times and variation in bed pressure and temperature are in good agreement with experimental observations obtained using a 1.82 L prototype ANG storage cylinder. Storage efficiencies based on the maximum achievable V/V (volume of gas/volume of container) and filling times are used to quantify the performance of the charging process. For the high permeability beds used in the experiments, storage efficiencies are controlled by the rate of heat removal. Filling times, defined as the time at which the bed pressure reaches 3.5 MPa, range from 120 to 3.4 min for inlet flow rates of 1.0 L min(-1) and 30.0 L min(-1), respectively. The corresponding storage efficiencies, eta(s), vary from 90% to 76%, respectively. Simulations with L/D ratios ranging from 0.35 to 7.8 indicate that the storage efficiencies can be improved with an increase in the LID ratios and/or with water cooled convection. Thus for an inlet flow rate of 30.0 L min(-1), an eta(s) value of 90% can be obtained with water cooling for an L/D ratio of 7.8 and a filling time of a few minutes. In the absence of water cooling the eta(s) value reduces to 83% at the same L/D ratio. Our study suggests that with an appropriate choice of cylinder dimensions, solutions based on convective cooling during adsorptive storage are possible with some compromise in the storage capacity.

Queuing Network Models for Aircraft Availability and Spares Management

Interactions of major activities involved in airfleet operations, maintenance, and logistics are investigated in the framework of closed queuing networks with finite number of customers. The system is viewed at three levels, namely: operations at the flying-base, maintenance at the repair-depot, and logistics for subsystems and their interactions in achieving the system objectives. Several performance measures (eg, availability of aircraft at the flying-base, mean number of aircraft on ground at different stages of repair, use of repair facilities, and mean time an aircraft spends in various stages of repair) can easily be computed in this framework. At the subsystem level the quantities of interest are the unavailability (probability of stockout) of a spare and the duration of its unavailability. The repair-depot capability is affected by the unavailability of a spare which in turn, adversely affects the availability of aircraft at the flying-base level. Examples illustrate the utility of the proposed models.

The Goldman bracket characterizes homeomorphisms

We show that a homotopy equivalence between compact, connected, oriented surfaces with non-empty boundary is homotopic to a homeomorphism if and only if it commutes with the Goldman bracket. (C) 2011 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Males, but not females, mate with multiple partners: a laboratory study of a primitively eusocial wasp Ropalidia marginata

The intense interest in social Hymenoptera, on account of their elaborate sociality and the paradox of altruism, has often suffered from considerable gender imbalance. This is partly due to the fact that worker behaviour and altruism are restricted to the females and partly because males often live off the nest. Yet, understanding the males, especially in the context of mating biology is essential even for understanding the evolution of sociality. Mating patterns have a direct bearing on the levels of intra-colony genetic relatedness, which in turn, along with the associated costs and benefits of worker behaviour, are central to our understanding of the evolution of sociality. Although mating takes place away from the nest in natural colonies of the primitively eusocial wasp Ropalidia marginata, mating can be observed in the laboratory if a male and a female are placed in a transparent, aerated plastic container, and both wasps are in the range of 5-20 days of age. Here, we use this setup and show that males, but not females, mate serially with multiple partners. The multiple mating behaviour of the males is not surprising because in nature males have to mate with a number of females, only a few of whom will go on to lay eggs. The reluctance of R. marginata females to mate with multiple partners is consistent with the expectation of monogamy in primitively eusocial species with totipotent females, although the apparent discrepancy with a previous work with allozyme markers in natural colonies suggesting that females may sometimes mate with two or three different males remains to be resolved.

New nine-node Lagrangian quadrilateral plate element based on Mindlin-Reissner theory using IFM

This paper presents a new nine-node Lagrangian quadrilateral plate bending element (MQP9) using the Integrated Force Method (IFM) for the analysis of thin and moderately thick plate bending problems. Three degrees of freedom: transverse displacement w and two rotations theta(x) and theta(y) are considered at each node of the element. The Mindlin-Reissner theory has been employed in the formulation which accounts the effect of shear deformation. Many standard plate bending benchmark problems have been analyzed using the new element MQP9 for various grid sizes via Integrated Force Method to estimate defections and bending moments. These results of the new element MQP9 are compared with those of similar displacement-based plate bending elements available in the literature. The results are also compared with exact solutions. It is observed that the presented new element MQP9 is free from shear locking and produced, in general, excellent results in all plate bending benchmark problems considered.

Ultrasonic Wave Dispersion Characteristics of Fluid-Filled Single-Walled Carbon Nanotubes Based on Nonclassical Shell Model

In this paper, ultrasonic wave propagation analysis in fluid filled single-walled carbon nanotube (SWCNT) is studied using nonlocal elasticity theory. The SWCNT is modeled using Flugge's shell theory, with the wall having axial, circumferential and radial degrees of freedom and also including small scale effects. The fluid inside the SWCNT is assumed as water. Nonlocal governing equations for this system are derived and wave propagation analysis is also carried out. The presence of fluid in SWCNT alters the ultrasonic wave dispersion behavior. The wavenumber and wave velocity are smaller in presence of fluid as compared to the empty SWCNT. The nonlocal elasticity calculation shows that the wavenumber tends to reach the continuum limit at certain frequencies and the corresponding wave velocity tends to zero at those frequencies indicating localization and stationary behavior. It has been shown that the circumferential. waves will propagate non-dispersively at higher frequencies in nonlocality. The magnitudes of wave velocities of circumferential waves are smaller in nonlocal elasticity as compared to local elasticity. We also show that the cut-off frequency depend on the nonlocal scaling parameter and also on the density of the fluid inside the SWCNT, and the axial wavenumber, as the fluid becomes denser the cut-off frequency decreases. The effect of axial wavenumber on the ultrasonic wave behavior in SWCNTS filled with water is also discussed.

Influence of Bottom Bubbling Rate on Formation of Metal Emulsion in Al-Cu Alloy and Molten Salt System

In steel refining process, an increase of interfacial area between the metal and slag through the metal droplets emulsified into the slag, so-called ``metal emulsion'', is one prevailing view for improving the reaction rate. The formation of metal emulsion was experimentally evaluated using Al-Cu alloy as metal phase and chloride salt as slag phase under the bottom bubbling condition. Samples were collected from the center of the salt phase in the container. Large number of metal droplets were separated from the salt by dissolving it into water. The number, surface area, and weight of the droplets increased with the gas flow rate and have local maximum values. The formation and sedimentation rates of metal droplets were estimated using a mathematical model. The formation rate increased with the gas flow rate and has a local maximum value as a function of gas flow rate, while the sedimentation rate is independent of the gas flow rate under the bottom bubbling condition. Three types of formation mode of metal emulsion, which occurred by the rupture of metal film around the bubble, were observed using high speed camera. During the process, an elongated column covered with metal film was observed with the increasing gas flow rate. This elongated column sometimes reached to the top surface of the salt phase. In this case, it is considered that fine droplets were not formed and in consequence, the weight of metal emulsion decreased at higher gas flow rate.

Analysis of Parameter Values in the van der Waals and Platteeuw Theory for Methane Hydrates Using Monte Carlo Molecular Simulations

The van der Waals and Platteuw (vdVVP) theory has been successfully used to model the thermodynamics of gas hydrates. However, earlier studies have shown that this could be due to the presence of a large number of adjustable parameters whose values are obtained through regression with experimental data. To test this assertion, we carry out a systematic and rigorous study of the performance of various models of vdWP theory that have been proposed over the years. The hydrate phase equilibrium data used for this study is obtained from Monte Carlo molecular simulations of methane hydrates. The parameters of the vdWP theory are regressed from this equilibrium data and compared with their true values obtained directly from simulations. This comparison reveals that (i) methane-water interactions beyond the first cage and methane-methane interactions make a significant contribution to the partition function and thus cannot be neglected, (ii) the rigorous Monte Carlo integration should be used to evaluate the Langmuir constant instead of the spherical smoothed cell approximation, (iii) the parameter values describing the methane-water interactions cannot be correctly regressed from the equilibrium data using the vdVVP theory in its present form, (iv) the regressed empty hydrate property values closely match their true values irrespective of the level of rigor in the theory, and (v) the flexibility of the water lattice forming the hydrate phase needs to be incorporated in the vdWP theory. Since methane is among the simplest of hydrate forming molecules, the conclusions from this study should also hold true for more complicated hydrate guest molecules.

A comparative study of dragonfly inspired flapping wings actuated by single crystal piezoceramic

A dragonfly inspired flapping wing is investigated in this paper. The flapping wing is actuated from the root by a PZT-5H and PZN-7%PT single crystal unimorph in the piezofan configuration. The nonlinear governing equations of motion of the smart flapping wing are obtained using the Hamilton's principle. These equations are then discretized using the Galerkin method and solved using the method of multiple scales. Dynamic characteristics of smart flapping wings having the same size as the actual wings of three different dragonfly species Aeshna Multicolor, Anax Parthenope Julius and Sympetrum Frequens are analyzed using numerical simulations. An unsteady aerodynamic model is used to obtain the aerodynamic forces. Finally, a comparative study of performances of three piezoelectrically actuated flapping wings is performed. The numerical results in this paper show that use of PZN-7%PT single crystal piezoceramic can lead to considerable amount of wing weight reduction and increase of lift and thrust force compared to PZT-5H material. It is also shown that dragonfly inspired smart flapping wings actuated by single crystal piezoceramic are a viable contender for insect scale flapping wing micro air vehicles.