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.

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.

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.

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.

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.

Modified governing equations for gas-particle nozzle flows

A modified set of governing equations for gas-particle flows in nozzles is suggested to include the inertial forces acting on the particle phase. The problem of gas-particle flow through a nozzle is solved using a first order finite difference scheme. A suitable stability condition for the numerical scheme for gas-particle flows is defined. Results obtained from the present set of equations are compared with those of the previous set of equations. It is also found that present set of equations give results which are in good agreement with the experimental observation.

Investigations on the melting and bending modulus of polymer grafted bilayers using dissipative particle dynamics

Understanding the influence of polymer grafted bilayers on the physicomechanical properties of lipid membranes is important while developing liposomal based drug delivery systems. The melting characteristics and bending moduli of polymer grafted bilayers are investigated using dissipative particle dynamics simulations as a function of the amount of grafted polymer and lipid tail length. Simulations are carried out using a modified Andersen barostat, whereby the membrane is maintained in a tensionless state. For lipids made up of four to six tail beads, the transition from the low temperature L-beta phase to the L-alpha phase is lowered only above a grafting fraction of G(f)=0.12 for polymers made up of 20 beads. Below G(f)=0.12 small changes are observed only for the HT4 bilayer. The bending modulus of the bilayers is obtained as a function of G(f) from a Fourier analysis of the height fluctuations. Using the theory developed by Marsh Biochim. Biophys. Acta 1615, 33 (2003)] for polymer grafted membranes, the contributions to the bending modulus due to changes arising from the grafted polymer and bilayer thinning are partitioned. The contributions to the changes in kappa from bilayer thinning were found to lie within 11% for the lipids with four to six tail beads, increasing to 15% for the lipids containing nine tail beads. The changes in the area stretch modulus were also assessed and were found to have a small influence on the overall contribution from membrane thinning. The increase in the area per head group of the lipids was found to be consistent with the scalings predicted by self-consistent mean field results. (C) 2010 American Institute of Physics.

Biomimetic synthesis of silver nanoparticles by Citrus limon (lemon) aqueous extract and theoretical prediction of particle size

In the present study silver nanoparticles were rapidly synthesized at room temperature by treating silver ions with the Citrus limon (lemon) extract The effect of various process parameters like the reductant con centration mixing ratio of the reactants and the concentration of silver nitrate were studied in detail In the standardized process 10(-2) M silver nitrate solution was interacted for 411 with lemon Juice (2% citric acid concentration and 0 5% ascorbic acid concentration) in the ratio of 1 4(vol vol) The formation of silver nanoparticles was confirmed by Surface Plasmon Resonance as determined by UV-Visible spectra in the range of 400-500 nm X ray diffraction analysis revealed the distinctive facets (1 1 1 200 220 2 2 2 and 3 1 1 planes) of silver nanoparticles We found that citric acid was the principal reducing agent for the nanosynthesis process FT IR spectral studies demonstrated citric acid as the probable stabilizing agent Silver nanoparticles below 50 nm with spherical and spheroidal shape were observed from transmission electron microscopy The correlation between absorption maxima and particle sizes were derived for different UV-Visible absorption maxima (corresponding to different citric acid concentrations) employing MiePlot v 3 4 The theoretical particle size corresponding to 2% citric acid concentration was corn pared to those obtained by various experimental techniques like X ray diffraction analysis atomic force microscopy and transmission electron microscopy (C) 2010 Elsevier B V All rights reserved

Dense granular flow down an inclined plane A comparison between the hard particle model and soft particle simulations

The granular flow down an inclined plane is simulated using the discrete element (DE) technique to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model First, we examine the average coordination number for the particles in the flow down an inclined plane using the DE technique using the linear contact model with and without friction, and the Hertzian contact model with friction The simulations show that the average coordination number decreases below 1 for values of the spring stiffness corresponding to real materials, such as sand and glass, even when the angle of inclination is only 10 larger than the angle of repose Additional measures of correlations in the system, such as the fraction of particles with multibody contact, the force ratio (average ratio of the magnitudes of the largest and the second largest force on a particle), and the angle between the two largest forces on the particle, show no evidence of force chains or other correlated motions in the system An analysis of the bond-orientational order parameter indicates that the flow is in the random state, as in event-driven (ED) simulations V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are compared with the theory V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] and are found to be in quantitative agreement In addition, we also conduct a comparison of the collision frequency and the distribution of the precollisional relative velocities of particles in contact The strong correlation effects exhibited by these two quantities in event-driven simulations V Kumaran, J Fluid Mech 632, 145 (2009)] are also found in the DE simulations (C) 2010 American Institute of Physics doi 10 1063/1 3504660]

Grinding abrasive wear and associated particle size effect

The type of abrasion that the grinding medium experiences inside a ball mill is classified as high stress or grinding abrasion, because the stress levels at the surface of the medium exceed the yield stress of the metal when hard abrasives are crushed. During dry grinding of ores the medium undergoes not only abrasion but also erosion and impact. As all three mechanisms of wear occur simultaneously, it is difficult to follow the individual components of wear. However, it is possible to show that the overall kinetics of wear follows a simple power law of the type w = at(b), where w is the weight loss of the grinding medium for a specified grinding time t and a and b are constants. Experimental data, obtained from dry grinding of quartz for a wide range of times using AISI 52100 steel balls having various microstructures in a laboratory scale batch mill, are fitted to the proposed equation and the wear rate w is calculated from the first derivative of the equation. The mean particle sizes of the quartz charge DBAR corresponding to 50 and 80% retained size are determined by mechanical sieving of the ground product after a grinding time t and thus the relationship between wear rate and particle size of the abrasive is established. It is found that w increases rapidly with DBAR up to some critical size and then increases at a much lower rate.