Accurate Parameter Estimation of Contaminant Transport Inverse Problem using Particle Swarm Optimization

Swarm Intelligence techniques such as particle swarm optimization (PSO) are shown to be incompetent for an accurate estimation of global solutions in several engineering applications. This problem is more severe in case of inverse optimization problems where fitness calculations are computationally expensive. In this work, a novel strategy is introduced to alleviate this problem. The proposed inverse model based on modified particle swarm optimization algorithm is applied for a contaminant transport inverse model. The inverse models based on standard-PSO and proposed-PSO are validated to estimate the accuracy of the models. The proposed model is shown to be out performing the standard one in terms of accuracy in parameter estimation. The preliminary results obtained using the proposed model is presented in this work.

Oxide Particle Surface Chemistry and Ion Transport in "Soggy Sand" Electrolytes

The crucial role of oxide surface chemical composition on ion transport in "soggy sand" electrolytes is discussed in a systematic manner. A prototype soggy sand electrolytic system comprising aerosil silica functionalized with various hydrophilic and hydrophobic moieties dispersed in lithium perchlorate-ethylene glycol solution was used for the study. Detailed rheology studies show that the attractive particle network in the case of the composite with unmodified aerosil silica (with surface silanol groups) is most favorable for percolation in ionic conductivity, as well as rendering the composite with beneficial elastic mechanical properties: Though weaker in strength compared to the composite with unmodified aerosil particles, attractive particle networks are also observed in composites of aerosil particles with surfaces partially substituted with hydrophobic groups. The percolation in ionic conductivity is, however, dependent on the size of the hydrophobic moiety. No spanning attractive particle network was formed for aerosil particles with surfaces modified with stronger hydrophilic groups (than silanol), and as a result, no percolation in ionic conductivity was observed. The composite with hydrophilic particles was a sol, contrary to gels obtained in the case of unmodified aerosil, and partially substituted with hydrophobic groups.

Performance evaluation of elitist-mutated multi-objective particle swarm optimization for integrated water resources management

Optimal allocation of water resources for various stakeholders often involves considerable complexity with several conflicting goals, which often leads to multi-objective optimization. In aid of effective decision-making to the water managers, apart from developing effective multi-objective mathematical models, there is a greater necessity of providing efficient Pareto optimal solutions to the real world problems. This study proposes a swarm-intelligence-based multi-objective technique, namely the elitist-mutated multi-objective particle swarm optimization technique (EM-MOPSO), for arriving at efficient Pareto optimal solutions to the multi-objective water resource management problems. The EM-MOPSO technique is applied to a case study of the multi-objective reservoir operation problem. The model performance is evaluated by comparing with results of a non-dominated sorting genetic algorithm (NSGA-II) model, and it is found that the EM-MOPSO method results in better performance. The developed method can be used as an effective aid for multi-objective decision-making in integrated water resource management.

Nonlinear structural dynamical system identification using adaptive particle filters

The problem of identifying parameters of nonlinear vibrating systems using spatially incomplete, noisy, time-domain measurements is considered. The problem is formulated within the framework of dynamic state estimation formalisms that employ particle filters. The parameters of the system, which are to be identified, are treated as a set of random variables with finite number of discrete states. The study develops a procedure that combines a bank of self-learning particle filters with a global iteration strategy to estimate the probability distribution of the system parameters to be identified. Individual particle filters are based on the sequential importance sampling filter algorithm that is readily available in the existing literature. The paper develops the requisite recursive formulary for evaluating the evolution of weights associated with system parameter states. The correctness of the formulations developed is demonstrated first by applying the proposed procedure to a few linear vibrating systems for which an alternative solution using adaptive Kalman filter method is possible. Subsequently, illustrative examples on three nonlinear vibrating systems, using synthetic vibration data, are presented to reveal the correct functioning of the method. (c) 2007 Elsevier Ltd. All rights reserved.

Vector evaluated particle swarm optimization (VEPSO) for multi-objective design optimization of composite structures

We present a generic method/model for multi-objective design optimization of laminated composite components, based on vector evaluated particle swarm optimization (VEPSO) algorithm. VEPSO is a novel, co-evolutionary multi-objective variant of the popular particle swarm optimization algorithm (PSO). In the current work a modified version of VEPSO algorithm for discrete variables has been developed and implemented successfully for the, multi-objective design optimization of composites. The problem is formulated with multiple objectives of minimizing weight and the total cost of the composite component to achieve a specified strength. The primary optimization variables are - the number of layers, its stacking sequence (the orientation of the layers) and thickness of each layer. The classical lamination theory is utilized to determine the stresses in the component and the design is evaluated based on three failure criteria; failure mechanism based failure criteria, Maximum stress failure criteria and the Tsai-Wu failure criteria. The optimization method is validated for a number of different loading configurations - uniaxial, biaxial and bending loads. The design optimization has been carried for both variable stacking sequences, as well fixed standard stacking schemes and a comparative study of the different design configurations evolved has been presented. (C) 2007 Elsevier Ltd. All rights reserved.

Modeling of formation of gold nanoparticles by citrate method

Properties of nanoparticles are size dependent, and a model to predict particle size is of importance. Gold nanoparticles are commonly synthesized by reducing tetrachloroauric acid with trisodium citrate, a method pioneered by Turkevich et al (Discuss. Faraday Soc. 1951, 11, 55). Data from several investigators that used this method show that when the ratio of initial concentrations of citrate to gold is varied from 0.4 to similar to 2, the final mean size of the particles formed varies by a factor of 7, while subsequent increases in the ratio hardly have any effect on the size. In this paper, a model is developed to explain this widely varying dependence. The steps that lead to the formation of particles are as follows: reduction of Au3+ in solution, disproportionation of Au+ to gold atoms and their nucleation, growth by disproportionation on particle surface, and coagulation. Oxidation of citrate results in the formation of dicarboxy acetone, which aids nucleation but also decomposes into side products. A detailed kinetic model is developed on the basis of these steps and is combined with population balance to predict particle-size distribution. The model shows that, unlike the usual balance between nucleation and growth that determines the particle size, it is the balance between rate of nucleation and degradation of dicarboxy acetone that determines the particle size in the citrate process. It is this feature that is able to explain the unusual dependence of the mean particle size on the ratio of citrate to gold salt concentration. It is also found that coagulation plays an important role in determining the particle size at high concentrations of citrate.

Vector Evaluated Particle Swarm Optimization (VEPSO) of Supersonic Ejector for Hydrogen Fuel Cells

Fuel cells are emerging as alternate green power producers for both large power production and for use in automobiles. Hydrogen is seen as the best option as a fuel; however, hydrogen fuel cells require recirculation of unspent hydrogen. A supersonic ejector is an apt device for recirculation in the operating regimes of a hydrogen fuel cell. Optimal ejectors have to be designed to achieve best performances. The use of the vector evaluated particle swarm optimization technique to optimize supersonic ejectors with a focus on its application for hydrogen recirculation in fuel cells is presented here. Two parameters, compression ratio and efficiency, have been identified as the objective functions to be optimized. Their relation to operating and design parameters of ejector is obtained by control volume based analysis using a constant area mixing approximation. The independent parameters considered are the area ratio and the exit Mach number of the nozzle. The optimization is carried out at a particularentrainment ratio and results in a set of nondominated solutions, the Pareto front. A set of such curves can be used for choosing the optimal design parameters of the ejector.

Strength design of composite beam using gradient and particle swarm optimization

This work addresses the optimum design of a composite box-beam structure subject to strength constraints. Such box-beams are used as the main load carrying members of helicopter rotor blades. A computationally efficient analytical model for box-beam is used. Optimal ply orientation angles are sought which maximize the failure margins with respect to the applied loading. The Tsai-Wu-Hahn failure criterion is used to calculate the reserve factor for each wall and ply and the minimum reserve factor is maximized. Ply angles are used as design variables and various cases of initial starting design and loadings are investigated. Both gradient-based and particle swarm optimization (PSO) methods are used. It is found that the optimization approach leads to the design of a box-beam with greatly improved reserve factors which can be useful for helicopter rotor structures. While the PSO yields globally best designs, the gradient-based method can also be used with appropriate starting designs to obtain useful designs efficiently. (C) 2006 Elsevier Ltd. All rights reserved.

Mechanism of crystal structure transformations. Part 3. - Factors affecting the anatase-rutile transformation

The particle size and crystallite size of anatase increase markedly in the region of the crystal structure transformation. The unit cell of anatase seems to expand prior to the transformation to rutile. This expansion has been attributed to a displacive transformation of the type defined by Buerger. Smaller particle size and larger surface area seem to favour the transformation. The kinetics of the transformation of anatase prepared by the hydrolysis of titanium sulphate have been studied at different temperatures and are found to be considerably different from the kinetics of the transformation of pure anatase. The transformation becomes immeasurably slow below ∼695 ± 10°C compared to ∼610°C for pure anatase. An induction period is observed in the transformation of anatase obtained from sulphate hydrolysis and the duration decreases with increase in temperature. The activation energy is ∼120 kcal/mole, a value higher than that for the pure anatase-rutile transformation. The results have been interpreted in terms of the relative rates of nucleation and propagation processes. The activation energy for the nucleation process seems to be much larger than for the propagation process. The kinetics of the transformation of anatase samples doped with different amounts of sulphate ion impurity have also been studied and the transformation is found to be progressively decelerated with increase in the impurity concentration. The energy of activation for the transformation appears to increase progressively with increase in impurity concentration.

Influence of Oxide Particle Network Morphology on Ion Solvation and Transport in ``Soggy Sand''Electrolytes

The role of oxide surface chemical composition and solvent on ion solvation and ion transport of ``soggy sand'' electrolytes are discussed here. A ``soggy sand'' electrolyte system comprising dispersions of hydrophilic/hydrophobic functionalized aerosil silica in lithium perchlorate methoxy polyethylene glycol solution was employed for the study. Static and dynamic rheology measurements show formation of an attractive particle network in the case of the composite with unmodified aerosil silica (i.e., with surface silanol groups) as well as composites with hydrophobic alkane groups. While particle network in the composite with hydrophilic aerosil silica (unmodified) were due to hydrogen bonding, hydrophobic aerosil silica particles were held together via van der Waals forces. The network strength in the latter case (i.e., for hydrophobic composites) were weaker compared with the composite with unmodified aerosil silica. Both unmodified silica as well as hydrophobic silica composites displayed solid-like mechanical strength. No enhancement in ionic conductivity compared to the liquid electrolyte was observed in the case of the unmodified silica. This was attributed to the existence of a very strong particle network, which led to the ``expulsion'' of all conducting entities from the interfacial region between adjacent particles. The ionic conductivity for composites with hydrophobic aerosil particles displayed ionic conductivity dependent on the size of the hydrophobic chemical moiety. No spanning attractive particle network was observed for aerosil particles with surfaces modified with stronger hydrophilic groups (than silanol). The composite resembled a sol, and no percolation in ionic conductivity was observed.

A particle filtering approach for structural system identification in vehicle–structure interaction problems

The problem of identification of parameters of a beam-moving oscillator system based on measurement of time histories of beam strains and displacements is considered. The governing equations of motion here have time varying coefficients. The parameters to be identified are however time invariant and consist of mass, stiffness and damping characteristics of the beam and oscillator subsystems. A strategy based on dynamic state estimation method, that employs particle filtering algorithms, is proposed to tackle the identification problem. The method can take into account measurement noise, guideway unevenness, spatially incomplete measurements, finite element models for supporting structure and moving vehicle, and imperfections in the formulation of the mathematical models. Numerical illustrations based on synthetic data on beam-oscillator system are presented to demonstrate the satisfactory performance of the proposed procedure.

Growth, self-assembly and dynamics of nano-scale films at fluid interfaces

Ultrathin films at fluid interfaces are important not only from a fundamental point of view as 2D complex fluids but have also become increasingly relevant in the development of novel functional materials. There has been an explosion in the synthesis work in this area over the last decade, giving rise to many exotic nanostructures at fluid interfaces. However, the factors controlling particle nucleation, growth and self-assembly at interfaces are poorly understood on a quantitative level. We will outline some of the recent attempts in this direction. Some of the selected investigations examining the macroscopic mechanical properties of molecular and particulate films at fluid interfaces will be reviewed. We conclude with a discussion of the electronic properties of these films that have potential technological and biological applications.

Simulation of precipitation reactions in reverse micelles

Precipitation involving mixing of two sets of reverse micellar solutions-containing a reactant and precipitant respectively-has been analyzed. Particle formation in such systems has been simulated by a Monte Carlo (MC) scheme (Li, Y.; Park, C. W. Langmuir 1999, 15, 952), which however is very restrictive in its approach. We have simulated particle formation by developing a general Monte Carlo scheme, using the interval of quiescence technique (IQ). It uses Poisson distribution with realistic, low micellar occupancies of reactants, Brownian collision of micelles with coalescence efficiency, fission of dimers with binomial redispersion of solutes, finite nucleation rate of particles with critical number of molecules, and instantaneous particle growth. With the incorporation of these features, the previous work becomes a special case of our simulation. The present scheme was then used to predict experimental data on two systems. The first is the experimental results of Lianos and Thomas (Chem. Phys. Lett. 1986, 125, 299, J. Colloid Interface Sci. 1987, 117, 505) on formation of CdS nanoparticles. They reported the number of molecules in a particle as a function of micellar size and reactant concentrations, which have been predicted very well. The second is on the formation of Fe(OH)(3) nanoparticles, reported by Li and Park. Our simulation in this case provides a better prediction of the experimental particle size range than the prediction of the authors. The present simulation scheme is general and can be applied to explain nanoparticle formation in other systems.

Peptide hairpin nucleation with the obligatory Type I' beta-turn Aib-(D)Pro segment

The alpha-aminoisobutyric acid-D-proline (Aib-(D)Pro) dipeptide is an obligatory Type I' beta-turn forming segment that nucleates hairpin formation.

The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.