Coupled electrostatic-elastic analysis for topology optimization using material interpolation

In this paper, we present a novel analytical formulation for the coupled partial differential equations governing electrostatically actuated constrained elastic structures of inhomogeneous material composition. We also present a computationally efficient numerical framework for solving the coupled equations over a reference domain with a fixed finite-element mesh. This serves two purposes: (i) a series of problems with varying geometries and piece-wise homogeneous and/or inhomogeneous material distribution can be solved with a single pre-processing step, (ii) topology optimization methods can be easily implemented by interpolating the material at each point in the reference domain from a void to a dielectric or a conductor. This is attained by considering the steady-state electrical current conduction equation with a `leaky capacitor' model instead of the usual electrostatic equation. This formulation is amenable for both static and transient problems in the elastic domain coupled with the quasi-electrostatic electric field. The procedure is numerically implemented on the COMSOL Multiphysics (R) platform using the weak variational form of the governing equations. Examples have been presented to show the accuracy and versatility of the scheme. The accuracy of the scheme is validated for the special case of piece-wise homogeneous material in the limit of the leaky-capacitor model approaching the ideal case.

Performability of automated manufacturing systems with centralized material handling

We present a framework for performance evaluation of manufacturing systems subject to failure and repair. In particular, we determine the mean and variance of accumulated production over a specified time frame and show the usefulness of these results in system design and in evaluating operational policies for manufacturing systems. We extend this analysis for lead time as well. A detailed performability study is carried out for the generic model of a manufacturing system with centralized material handling. Several numerical results are presented, and the relevance of performability analysis in resolving system design issues is highlighted. Specific problems addressed include computing the distribution of total production over a shift period, determining the shift length necessary to deliver a given production target with a desired probability, and obtaining the distribution of Manufacturing Lead Time, all in the face of potential subsystem failures.

Composite material and piezoelectric coefficient uncertainty effects on structural health monitoring using feedback control gains as damage indicators

This article addresses uncertainty effect on the health monitoring of a smart structure using control gain shifts as damage indicators. A finite element model of the smart composite plate with surface-bonded piezoelectric sensors and actuators is formulated using first-order shear deformation theory and a matrix crack model is integrated into the finite element model. A constant gain velocity/position feedback control algorithm is used to provide active damping to the structure. Numerical results show that the response of the structure is changed due to matrix cracks and this change can be compensated by actively tuning the feedback controller. This change in control gain can be used as a damage indicator for structural health monitoring. Monte Carlo simulation is conducted to study the effect of material uncertainty on the damage indicator by considering composite material properties and piezoelectric coefficients as independent random variables. It is found that the change in position feedback control gain is a robust damage indicator.

Stabilized alpha-Ni(OH)2 as electrode material for for alkaline secondary cells

Hydrotalcite-like compounds of formula Ni1-xAl(x)(OH)2(CO3)x/2 . nH2O (x = 0.1 to 0.25), having the same structure as that of alpha-Ni(OH)2, have been synthesized by substituting nickel hydroxide with aluminum. Of these, the compounds of compositions x greater-than-or-equal-to 0.2 are found to have prolonged stability in strong alkaline medium. The electrodes comprising stabilized alpha-Ni(OH)2 of x = 0.2 composition are rechargeable with discharge-capacity values of 240 (+/- 15) mAh-g-1 and are attractive for applications in various alkaline secondary cells employing nickel-positive electrodes.

Influence of pH on the retention of heavy metal ions by fly ash

Clay liners have been widely used to contain toxic and hazardous wastes. Clays adsorb the contaminant cations due to their exchange capacity. To improve the performance of the clay liner, fly ash, a waste material arising out of combustion of coal has been studied as a pre-filter material. The results indicate that fly ash has the potential to retain heavy metal ions. This study concerns the retention of zinc by fly ash. The influence of pH on retention as well as leaching characteristics are examined. The results obtained from the retention experiments by permeameter method indicate that fly ash retains the zinc ions through precipitation in the pores as well as onto the surface when the ambient pH value is more than 6.9, and only through adsorption when the pH value is less than 6.9. It has been observed that fly ash did not release the retained zinc ions when the pH value is between 3.5 and 10.0. Hence, the retention of zinc ions by fly ash is likely to be permanent since the pH of most of the landfill leachates are between 3.7 to 8.8.

LaNiO3: a promising material for contact with YBa2 Cu3O7-x thin films

Epitaxial LaNiO3 metallic oxide thin films have been grown on c-axis oriented YBa2Cu3O7-delta thin films on LaAlO3 substrates by pulsed laser deposition technique and the interface formed between the two films has been examined by measuring the contact conductance of the same. The specific contact conductance of the interface measured using a modified four probe method was found to be 1.4 to 6 x 10(4) ohm(-1) cm(-2) at 77 K, There are indications that contact conductance can be brought closer to that obtained for noble metal-YBCO interface.

Effect of Melt Convection on The Optimum Thermal Design Of Heat Sinks With Phase Change Material

In this paper, the role of melt convection on the performance of heat sinks with phase change material (PCM) is investigated numerically. The heat sink consists of aluminum plate fins embedded in PCM, and is subjected to heat flux supplied from the bottom. A single-domain enthalpy-based CFD model is developed, which is capable of simulating the phase change process and the associated melt convection. The CFD model is coupled with a genetic algorithm for carrying out the optimization. Two cases are considered, namely, one without melt convection (i.e., conduction heat transfer analysis), and the other with convection. It is found that the geometrical optimizations of heat sinks are different for the two cases, indicating the importance of melt convection in the design of heat sinks with PCMs. In the case of conduction analysis, the optimum width of half fin (i.e., sum of half pitch and half fin thickness) is a constant, which is in good agreement with results reported in the literature. On the other hand, if melt convection is considered, the optimum half fin width depends on the effective thermal diffusivity due to conduction and convection. With melt convection, the optimized design results in a significant improvement of operational time.

The physics and technology of gallium antimonide: An emerging optoelectronic material

Recent advances in nonsilica fiber technology have prompted the development of suitable materials for devices operating beyond 1.55 mu m. The III-V ternaries and quaternaries (AlGaIn)(AsSb) lattice matched to GaSb seem to be the obvious choice and have turned out to be promising candidates for high speed electronic and long wavelength photonic devices. Consequently, there has been tremendous upthrust in research activities of GaSb-based systems. As a matter of fact, this compound has proved to be an interesting material for both basic and applied research. At present, GaSb technology is in its infancy and considerable research has to be carried out before it can be employed for large scale device fabrication. This article presents an up to date comprehensive account of research carried out hitherto. It explores in detail the material aspects of GaSb starting from crystal growth in bulk and epitaxial form, post growth material processing to device feasibility. An overview of the lattice, electronic, transport, optical and device related properties is presented. Some of the current areas of research and development have been critically reviewed and their significance for both understanding the basic physics as well as for device applications are addressed. These include the role of defects and impurities on the structural, optical and electrical properties of the material, various techniques employed for surface and bulk defect passivation and their effect on the device characteristics, development of novel device structures, etc. Several avenues where further work is required in order to upgrade this III-V compound for optoelectronic devices are listed. It is concluded that the present day knowledge in this material system is sufficient to understand the basic properties and what should be more vigorously pursued is their implementation for device fabrication. (C) 1997 American Institute of Physics.

Velocity distribution function for a dilute granular material in shear flow

The velocity distribution function for the steady shear flow of disks (in two dimensions) and spheres (in three dimensions) in a channel is determined in the limit where the frequency of particle-wall collisions is large compared to particle-particle collisions. An asymptotic analysis is used in the small parameter epsilon, which is naL in two dimensions and na(2)L in three dimensions, where; n is the number density of particles (per unit area in two dimensions and per unit volume in three dimensions), L is the separation of the walls of the channel and a is the particle diameter. The particle-wall collisions are inelastic, and are described by simple relations which involve coefficients of restitution e(t) and e(n) in the tangential and normal directions, and both elastic and inelastic binary collisions between particles are considered. In the absence of binary collisions between particles, it is found that the particle velocities converge to two constant values (u(x), u(y)) = (+/-V, O) after repeated collisions with the wall, where u(x) and u(y) are the velocities tangential and normal to the wall, V = (1 - e(t))V-w/(1 + e(t)), and V-w and -V-w, are the tangential velocities of the walls of the channel. The effect of binary collisions is included using a self-consistent calculation, and the distribution function is determined using the condition that the net collisional flux of particles at any point in velocity space is zero at steady state. Certain approximations are made regarding the velocities of particles undergoing binary collisions :in order to obtain analytical results for the distribution function, and these approximations are justified analytically by showing that the error incurred decreases proportional to epsilon(1/2) in the limit epsilon --> 0. A numerical calculation of the mean square of the difference between the exact flux and the approximate flux confirms that the error decreases proportional to epsilon(1/2) in the limit epsilon --> 0. The moments of the velocity distribution function are evaluated, and it is found that [u(x)(2)] --> V-2, [u(y)(2)] similar to V-2 epsilon and -[u(x)u(y)] similar to V-2 epsilon log(epsilon(-1)) in the limit epsilon --> 0. It is found that the distribution function and the scaling laws for the velocity moments are similar for both two- and three-dimensional systems.

Temperature scaling in a dense vibrofluidized granular material

The leading order "temperature" of a dense two-dimensional granular material fluidized by external vibrations is determined. The grain interactions are characterized by inelastic collisions, but the coefficient of restitution is considered to be close to 1, so that the dissipation of energy during a collision is small compared to the average energy of a particle. An asymptotic solution is obtained where the particles are considered to be elastic in the leading approximation. The velocity distribution is a Maxwell-Boltzmann distribution in the leading approximation,. The density profile is determined by solving the momentum balance equation in the vertical direction, where the relation between the pressure and density is provided by the virial equation of state. The temperature is determined by relating the source of energy due to the vibrating surface and the energy dissipation due to inelastic collisions. The predictions of the present analysis show good agreement with simulation results at higher densities where theories for a dilute vibrated granular material, with the pressure-density relation provided by the ideal gas law, sire in error. [:S1063-651X(99)04408-6].

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

Zinc-substituted alpha-nickel hydroxide as an electrode material for alkaline secondary cells

Double hydroxides of the formula, Ni1-xZn2x (OH)(2) (CO3)(x). nH(2)O (x = 0.1 to 0.25) having the same structure as that of alpha-nickel hydroxide have been synthesized by partial substitution of zinc for nickel. The hydroxide having the composition x = 0.25 exhibits prolonged stability in 6 M KOH. Pasted electrodes comprising this material are rechargeable with a stabilized reversible discharge capacity of 410 +/- 15 mAh g(-1) of nickel even under suboptimal conditions of electrode fabrication. This compares favorably with the capacity values achieved for beta-nickel hydroxide (221 mAh g(-1)', This work; 297 mAh g(-1), Delahaye-Vidal and Figlarz;(1) 456 mAh g(-1), theoretical). (C) 1999 The Electrochemical Society. S0013-4651(98)01-071-4. All rights reserved.

Lithium borate-strontium bismuth tantalate glass nanocomposite: a novel material for nonlinear optic and ferroelectric applications

Glass nanocomposites in the system (100 - x)Li2B4O7-xSrBi(2)Ta(2)O(9) (0 less than or equal to x less than or equal to 22.5, in molar ratio) were fabricated via a melt quenching technique followed by controlled heat-treatment. The as-quenched samples were confirmed to be glassy and amorphous by differential thermal analysis (DTA) and X-ray powder diffraction (XRD) techniques, respectively. The phase formation and crystallite size of the heat-treated samples (glass nanocomposites) were monitored by XRD and transmission electron microscopy (TEM). The relative permittivities (epsilon(tau)') of the glass nanocomposites for different compositions were found to lie in between that of the parent host glass (Li2B4O7) and strontium bismuth tantalate (SBT) ceramic in the frequency range 100 Hz-40 MHz at 300 K, whereas the dielectric loss (D) of the glass nanocomposite was less than that of both the parent phases. Among the various dielectric models employed to predict the effective relative permittivity of the glass nanocomposite, the one obtained using the Maxwell's model was in good agreement with the experimentally observed value. Impedance analysis was employed to rationalize the electrical behavior of the glasses and glass nanocomposites. The pyroelectric response of the glasses and glass nanocomposites was monitored as a function of temperature and the pyroelectric coefficient for glass and glass nanocomposite (x = 20) at 300 K were 27 muC m(-2) K-1 and 53 muC m(-2) K-1, respectively. The ferroelectric behavior of these glass nanocomposites was established by P vs. E hysteresis loop studies. The remnant polarization (P-r) of the glass nanocomposite increases with increase in SBT content. The coercive field (E-c) and P-r for the glass nanocomposite (x = 20) were 727 V cm(-1) and 0.527 muC cm(-2), respectively. The optical transmission properties of these glass nanocomposites were found to be composition dependent. The refractive index (n = 1.722), optical polarizability (am = 1.266 6 10 23 cm 3) and third-order nonlinear optical susceptibility (x(3) = 3.046 6 10(-21) cm(3)) of the glass nanocomposite (x = 15) were larger than those of the as-quenched glass. Second harmonic generation (SHG) was observed in transparent glass nanocomposites and the d(eff) for the glass nanocomposite (x = 20) was found to be 0.373 pm V-1.

Studies on compressive failure features in syntactic foam material

Syntactic foam made by mechanical mixing of glass hollow spheres in epoxy resin matrix is characterized for compressive properties in the present study. Volume fraction of hollow spheres in the syntactic foam under investigation is kept at 67.8%. Effect of specimen aspect ratio on failure behavior and stress-strain curve of the material is highlighted. Considerable differences are noted in the macroscopic fracture features of the specimen and the stress-strain curve with the variation in specimen aspect ratio, although compressive yield strength values were within a narrow range. Post compression test scanning electron microscopic observations coupled with the macroscopic observations taken during the test helped in explaining the deviation in specimen behavior and in gathering support for the proposed arguments.