A perturbed martingale approach to global optimization

A new global stochastic search, guided mainly through derivative-free directional information computable from the sample statistical moments of the design variables within a Monte Carlo setup, is proposed. The search is aided by imparting to the directional update term additional layers of random perturbations referred to as `coalescence' and `scrambling'. A selection step, constituting yet another avenue for random perturbation, completes the global search. The direction-driven nature of the search is manifest in the local extremization and coalescence components, which are posed as martingale problems that yield gain-like update terms upon discretization. As anticipated and numerically demonstrated, to a limited extent, against the problem of parameter recovery given the chaotic response histories of a couple of nonlinear oscillators, the proposed method appears to offer a more rational, more accurate and faster alternative to most available evolutionary schemes, prominently the particle swarm optimization. (C) 2014 Elsevier B.V. All rights reserved.

Phthalate mediated hydrogelation of a pyrene based system: a novel scaffold for shape-persistent, self-healing luminescent soft material

Two-component super-hydrogelation triggered by the acid-base interaction of a L-histidine appended pyrenyl derivative (PyHis) and phthalic acid (PA) was reported. The use of isomeric isophthalic or terephthalic acid or other comparable acids in place of PA does not lead to salt formation and therefore hydrogelation is not observed. Excimer formation of the pyrenyl unit has not been detected although the PyHis : PA = 1: 1 system undergoes extensive self-assembly in aqueous solution. The synergistic effect of intermolecular H-bonding forces, pi-pi stacking, electrostatic interactions, etc. is found to be responsible for robust hydrogel formation. Development of chiral supramotecular assemblies has been verified through circular dichroism spectroscopy. Morphological investigations involving the PyHis : PA = 1: 1 system show vesicular nano-structures with a definite bilayer width at relatively low concentrations. The latter fuses to construct coiled-coil left-handed helical fibers upon increase in the concentrations of the gelators. The intertwining of the resultant helical fibers eventually results in hydrogel formation. The probable bilayer packing in the self-assembled structures has been probed using X-ray diffraction (XRD) studies and lanthanide sensitization, which suggests that the polar imidazolium hydrogen phthalate unit of the gelator forms the head group and faces the hydrophilic water environment while the hydrophobic pyrenyl units sit inside the hydrophobic core of the bilayer. The hydrogel exhibits multi-stimuli responsiveness including thixotropic behavior. In addition, shape-persistent as well as rapid self-healing behaviour of the hydrogel was established. Furthermore load-bearing characteristics of the hydrogel have also been demonstrated.

Newton-based stochastic optimization using q-Gaussian smoothed functional algorithms

We present the first q-Gaussian smoothed functional (SF) estimator of the Hessian and the first Newton-based stochastic optimization algorithm that estimates both the Hessian and the gradient of the objective function using q-Gaussian perturbations. Our algorithm requires only two system simulations (regardless of the parameter dimension) and estimates both the gradient and the Hessian at each update epoch using these. We also present a proof of convergence of the proposed algorithm. In a related recent work (Ghoshdastidar, Dukkipati, & Bhatnagar, 2014), we presented gradient SF algorithms based on the q-Gaussian perturbations. Our work extends prior work on SF algorithms by generalizing the class of perturbation distributions as most distributions reported in the literature for which SF algorithms are known to work turn out to be special cases of the q-Gaussian distribution. Besides studying the convergence properties of our algorithm analytically, we also show the results of numerical simulations on a model of a queuing network, that illustrate the significance of the proposed method. In particular, we observe that our algorithm performs better in most cases, over a wide range of q-values, in comparison to Newton SF algorithms with the Gaussian and Cauchy perturbations, as well as the gradient q-Gaussian SF algorithms. (C) 2014 Elsevier Ltd. All rights reserved.

Tailoring the second mode of Euler-Bernoulli beams: an analytical approach

In this paper, we study the inverse mode shape problem for an Euler-Bernoulli beam, using an analytical approach. The mass and stiffness variations are determined for a beam, having various boundary conditions, which has a prescribed polynomial second mode shape with an internal node. It is found that physically feasible rectangular cross-section beams which satisfy the inverse problem exist for a variety of boundary conditions. The effect of the location of the internal node on the mass and stiffness variations and on the deflection of the beam is studied. The derived functions are used to verify the p-version finite element code, for the cantilever boundary condition. The paper also presents the bounds on the location of the internal node, for a valid mass and stiffness variation, for any given boundary condition. The derived property variations, corresponding to a given mode shape and boundary condition, also provides a simple closed-form solution for a class of non-uniform Euler-Bernoulli beams. These closed-form solutions can also be used to check optimization algorithms proposed for modal tailoring.

A bottom-up approach for optimization of friction stir processing parameters; a study on aluminium 2024-T3 alloy

Friction stir processing (FSP) is emerging as one of the most competent severe plastic deformation (SPD) method for producing bulk ultra-fine grained materials with improved properties. Optimizing the process parameters for a defect free process is one of the challenging aspects of FSP to mark its commercial use. For the commercial aluminium alloy 2024-T3 plate of 6 mm thickness, a bottom-up approach has been attempted to optimize major independent parameters of the process such as plunge depth, tool rotation speed and traverse speed. Tensile properties of the optimum friction stir processed sample were correlated with the microstructural characterization done using Scanning Electron Microscope (SEM) and Electron Back-Scattered Diffraction (EBSD). Optimum parameters from the bottom-up approach have led to a defect free FSP having a maximum strength of 93% the base material strength. Micro tensile testing of the samples taken from the center of processed zone has shown an increased strength of 1.3 times the base material. Measured maximum longitudinal residual stress on the processed surface was only 30 MPa which was attributed to the solid state nature of FSP. Microstructural observation reveals significant grain refinement with less variation in the grain size across the thickness and a large amount of grain boundary precipitation compared to the base metal. The proposed experimental bottom-up approach can be applied as an effective method for optimizing parameters during FSP of aluminium alloys, which is otherwise difficult through analytical methods due to the complex interactions between work-piece, tool and process parameters. Precipitation mechanisms during FSP were responsible for the fine grained microstructure in the nugget zone that provided better mechanical properties than the base metal. (C) 2014 Elsevier Ltd. All rights reserved.

Simultaneous Perturbation Newton Algorithms for Simulation Optimization

We present a new Hessian estimator based on the simultaneous perturbation procedure, that requires three system simulations regardless of the parameter dimension. We then present two Newton-based simulation optimization algorithms that incorporate this Hessian estimator. The two algorithms differ primarily in the manner in which the Hessian estimate is used. Both our algorithms do not compute the inverse Hessian explicitly, thereby saving on computational effort. While our first algorithm directly obtains the product of the inverse Hessian with the gradient of the objective, our second algorithm makes use of the Sherman-Morrison matrix inversion lemma to recursively estimate the inverse Hessian. We provide proofs of convergence for both our algorithms. Next, we consider an interesting application of our algorithms on a problem of road traffic control. Our algorithms are seen to exhibit better performance than two Newton algorithms from a recent prior work.

Manipulating Crystallographic Texture of Sn Coatings by Optimization of Electrodeposition Process Conditions to Suppress Growth of Whiskers

The effects of two major electrodeposition process conditions, electrolyte bath temperature and current density, on the microstructure and crystallographic texture of pure tin coatings on brass and, ultimately, on the extent of whisker formation have been examined. The grain size of the deposited coatings increased with increasing electrolyte bath temperature and current density, which significantly affected the dominant texture: (211) or (420) was the dominant texture at low current densities whereas, depending on deposition temperature, (200) or (220) became the dominant texture at high current densities. After deposition, coatings were subjected to different environmental conditions, for example isothermal aging (room temperature, 50A degrees C, or 150A degrees C) for up to 90 days and thermal cycling between -25A degrees C and 85A degrees C for 100 cycles, and whisker growth was studied. The Sn coatings with low Miller index planes, for example (200) and (220), and with moderate aging temperature were more prone to whiskering than coating with high Miller index planes, for example (420), and high aging temperature. A processing route involving the optimum combination of current density and deposition temperature is proposed for suppressing whisker growth.

Microstructure dependent elastic modulus variation in NiTi shape memory alloy

Modulus variation of NiTi shape memory alloy has been investigated at microstructural level through nano dynamical mechanical analysis and compared with bulk experimental measurements. The differences between the modulus values at the macro and micro level as well as within the micro level are discussed and the corresponding variations have been explained based on the crystal structure, orientation and misorientation. The experimental results confirm a higher modulus value for the martensite phase that is in agreement with the theoretical predictions. (C) 2015 Elsevier B. V. All rights reserved.

Shape and size directed self-selection in organic cage formation

The selective formation of a single isomer of a 3+2] self-assembled organic cage from a reaction mixture of an unsymmetrical aldehyde and a flexible amine is discussed. The experimental and theoretical findings suggest that in such a process, the geometric features of the aldehyde play a key role.

Effects of Interface Treatment on the Fatigue Behavior of Shape Memory Alloy Reinforced Polymer Composites

Interfacial properties of Shape Memory Alloy (SMA) reinforced polymer matrix composites can be enhanced by improving the interfacial bonding. This paper focuses on studying the interfacial stresses developed in the SMA-epoxy interface due to various laser shot penning conditions. Fiber-pull test-setup is designed to understand the role of mechanical bias stress cycling and thermal actuation cycling. Phase transformation is tracked over mechanical and thermal fatigue cycles. A micromechanics based model developed earlier based on shear lag in SMA and energy based consistent homogenization is extended here to incorporate the stress-temperature phase diagram parameters for modeling fatigue.

A Polymatroid Approach to Separable Convex Optimization with Linear Ascending Constraints

We revisit a problem studied by Padakandla and Sundaresan SIAM J. Optim., August 2009] on the minimization of a separable convex function subject to linear ascending constraints. The problem arises as the core optimization in several resource allocation problems in wireless communication settings. It is also a special case of an optimization of a separable convex function over the bases of a specially structured polymatroid. We give an alternative proof of the correctness of the algorithm of Padakandla and Sundaresan. In the process we relax some of their restrictions placed on the objective function.

Drift parameters optimization of a TPC polarimeter: a simulation study

Time Projection Chamber (TPC) based X-ray polarimeters using Gas Electron Multiplier (GEM) are currently being developed to make sensitive measurement of polarization in 2-10 keV energy range. The emission direction of the photoelectron ejected via photoelectric effect carries the information of the polarization of the incident X-ray photon. Performance of a gas based polarimeter is affected by the operating drift parameters such as gas pressure, drift field and drift-gap. We present simulation studies carried out in order to understand the effect of these operating parameters on the modulation factor of a TPC polarimeter. Models of Garfield are used to study photoelectron interaction in gas and drift of electron cloud towards GEM. Our study is aimed at achieving higher modulation factors by optimizing drift parameters. Study has shown that Ne/DME (50/50) at lower pressure and drift field can lead to desired performance of a TPC polarimeter.

Optimization of clamped beam geometry for fracture toughness testing of micron-scale samples

Fracture toughness measurements at the small scale have gained prominence over the years due to the continuing miniaturization of structural systems. Measurements carried out on bulk materials cannot be extrapolated to smaller length scales either due to the complexity of the microstructure or due to the size and geometric effect. Many new geometries have been proposed for fracture property measurements at small-length scales depending on the material behaviour and the type of device used in service. In situ testing provides the necessary environment to observe fracture at these length scales so as to determine the actual failure mechanism in these systems. In this paper, several improvements are incorporated to a previously proposed geometry of bending a doubly clamped beam for fracture toughness measurements. Both monotonic and cyclic loading conditions have been imposed on the beam to study R-curve and fatigue effects. In addition to the advantages that in situ SEM-based testing offers in such tests, FEM has been used as a simulation tool to replace cumbersome and expensive experiments to optimize the geometry. A description of all the improvements made to this specific geometry of clamped beam bending to make a variety of fracture property measurements is given in this paper.

Multiobjective Optimization of Piezoelectric Bimorph Actuator with Rigid Extension

Production of high tip deflection in a piezoelectric bimorph laminar actuator by applying high voltage is limited by many physical constraints. Therefore, piezoelectric bimorph actuator with a rigid extension of non-piezoelectric material at its tip is used to increase the tip deflection of such an actuator. Research on this type of piezoelectric bending actuator is either limited to first order constitutive relations, which do not include non-linear behavior of piezoelectric element at high electric field, or limited to curve fitting techniques. Therefore, this paper considers high electric field, and analytically models tapered piezoelectric bimorph actuator with a rigid extension of non-piezoelectric material at its tip. The stiffness, capacitance, effective tip deflection, block force, output strain energy, output energy density, input electrical energy and energy efficiency of the actuator are calculated analytically. The paper also discusses the multi-objective optimization of this type of actuator subjected to the mechanical and electrical constraints.

Exploring isospectral cantilever beams using electromagnetism inspired optimization technique

Isospectral beams have identical free vibration frequency spectrum for a specific boundary condition. The problem of finding non-uniform beams which are isospectral to a given uniform beam, with fixed-free boundary condition, leads to a multimodal optimization problem. The first Q natural frequencies of the given uniform Euler-Bernoulli beam are determined using analytical solution. The first Q natural frequencies of a non-uniform beam are obtained with the help of finite element modeling. In order to obtain the non-uniform beams isospectral to a given uniform beam, an error function is designed, which calculates the difference between the spectra of the given uniform beam and the non-uniform beam. In our study, this error function is minimized using electromagnetism inspired optimization technique, a population based iterative algorithm inspired by the attraction-repulsion physics of electromagnetism. Numerical results show the existence of the isospectral non-uniform beams for a given uniform beam, which occur as local minima. Non-uniform beams isospectral to a damaged beam, are also explored using the proposed methodology to illustrate the fact that accurate structural damage identification is difficult by just frequency measurements. (C) 2012 Elsevier B.V. All rights reserved.