A numerical study of the effect of loading conditions on tubular hydroforming

In this work, the mechanics of tubular hydroforming under various types of loading conditions is investigated. The main objective is to contrast the effects of prescribing fluid pressure or volume flow rate, in conjunction with axial displacement, on the stress and strain histories experienced by the tube and the process of bulging. To this end, axisymmetric finite element simulations of free hydroforming (without external die contact) of aluminium alloy tubes are carried out. Hill’s normally anisotropic yield theory along with material properties determined in a previous experimental study [A. Kulkarni, P. Biswas, R. Narasimhan, A. Luo, T. Stoughton, R. Mishra, A.K. Sachdev, An experimental and numerical study of necking initiation in aluminium alloy tubes during hydroforming, Int. J. Mech. Sci. 46 (2004) 1727–1746] are employed in the computations. It is found that while prescribed fluid pressure leads to highly non-proportional strain paths, specified fluid volume flow rate may result in almost proportional ones for the predominant portion of loading. The peak pressure increases with axial compression for the former, while the reverse trend applies under the latter. The implication of these results on failure by localized necking of the tube wall is addressed in a subsequent investigation.

A one-step technique to prepare aligned arrays of carbon nanotubes

A simple effective pyrolysis technique has been developed to synthesize aligned arrays of multi-walled carbon nanotubes (MWCNTs) without using any carrier gas in a single-stage furnace at 700 °C. This technique eliminates nearly the entire complex and expensive machinery associated with other extensively used methods for preparation of CNTs such as chemical vapour deposition (CVD) and pyrolysis. Carbon source materials such as xylene, cyclohexane, camphor, hexane, toluene, pyridine and benzene have been pyrolyzed separately with the catalyst source material ferrocene to obtain aligned arrays of MWCNTs. The synthesized CNTs have been characterized by scanning electron microscopy (SEM), x-ray diffraction (XRD), transmission electron microscopy (TEM), thermogravimetric analysis (TGA) and Raman spectroscopy. In this technique, the need for the tedious and time-consuming preparation of metal catalysts and continuously fed carbon source material containing carrier gas can be avoided. This method is a single-step process where not many parameters are required to be monitored in order to prepare aligned MWCNTs. For the production of CNTs, the technique has great advantages such as low cost and easy operation.

A new framework for solution of multidimensional population balance equations

A new framework is proposed in this work to solve multidimensional population balance equations (PBEs) using the method of discretization. A continuous PBE is considered as a statement of evolution of one evolving property of particles and conservation of their n internal attributes. Discretization must therefore preserve n + I properties of particles. Continuously distributed population is represented on discrete fixed pivots as in the fixed pivot technique of Kumar and Ramkrishna [1996a. On the solution of population balance equation by discretization-I A fixed pivot technique. Chemical Engineering Science 51(8), 1311-1332] for 1-d PBEs, but instead of the earlier extensions of this technique proposed in the literature which preserve 2(n) properties of non-pivot particles, the new framework requires n + I properties to be preserved. This opens up the use of triangular and tetrahedral elements to solve 2-d and 3-d PBEs, instead of the rectangles and cuboids that are suggested in the literature. Capabilities of computational fluid dynamics and other packages available for generating complex meshes can also be harnessed. The numerical results obtained indeed show the effectiveness of the new framework. It also brings out the hitherto unknown role of directionality of the grid in controlling the accuracy of the numerical solution of multidimensional PBEs. The numerical results obtained show that the quality of the numerical solution can be improved significantly just by altering the directionality of the grid, which does not require any increase in the number of points, or any refinement of the grid, or even redistribution of pivots in space. Directionality of a grid can be altered simply by regrouping of pivots.

Testing and Modeling of MR Damper and Its Application to SDOF Systems Using Integral Backstepping Technique

Magnetorheological dampers are intrinsically nonlinear devices, which make the modeling and design of a suitable control algorithm an interesting and challenging task. To evaluate the potential of magnetorheological (MR) dampers in control applications and to take full advantages of its unique features, a mathematical model to accurately reproduce its dynamic behavior has to be developed and then a proper control strategy has to be taken that is implementable and can fully utilize their capabilities as a semi-active control device. The present paper focuses on both the aspects. First, the paper reports the testing of a magnetorheological damper with an universal testing machine, for a set of frequency, amplitude, and current. A modified Bouc-Wen model considering the amplitude and input current dependence of the damper parameters has been proposed. It has been shown that the damper response can be satisfactorily predicted with this model. Second, a backstepping based nonlinear current monitoring of magnetorheological dampers for semi-active control of structures under earthquakes has been developed. It provides a stable nonlinear magnetorheological damper current monitoring directly based on system feedback such that current change in magnetorheological damper is gradual. Unlike other MR damper control techniques available in literature, the main advantage of the proposed technique lies in its current input prediction directly based on system feedback and smooth update of input current. Furthermore, while developing the proposed semi-active algorithm, the dynamics of the supplied and commanded current to the damper has been considered. The efficiency of the proposed technique has been shown taking a base isolated three story building under a set of seismic excitation. Comparison with widely used clipped-optimal strategy has also been shown.

Low power pipelined TCAM employing mismatch dependent power allocation technique

This paper presents design of a Low power 256x72 bit TCAM in 0.13um CMOS technology. In contrast to conventional Match line (ML) sensing scheme in which equal power is consumed irrespective of match or mismatch, the ML scheme employed in this design allocates less power to match decisions involving a large number of mismatched bits. Typically, the probability of mismatch is high so this scheme results in significant CAM power reduction. We propose to use this technique along with pipelining of search operation in which the MLs are broken into several segments. Since most words fail to match in first segment, the search operation for subsequent segments is discontinued, resulting in further reduction in power consumption. The above architecture provides 70% power reduction while performing search in 3ns.

Landau-Zener problem with waiting at the minimum gap and related quench dynamics of a many-body system

We discuss a technique for solving the Landau-Zener (LZ) problem of finding the probability of excitation in a two-level system. The idea of time reversal for the Schrodinger equation is employed to obtain the state reached at the final time and hence the excitation probability. Using this method, which can reproduce the well-known expression for the LZ transition probability, we solve a variant of the LZ problem, which involves waiting at the minimum gap for a time t(w); we find an exact expression for the excitation probability as a function of t(w). We provide numerical results to support our analytical expressions. We then discuss the problem of waiting at the quantum critical point of a many-body system and calculate the residual energy generated by the time-dependent Hamiltonian. Finally, we discuss possible experimental realizations of this work.

Numerical Simulation of Groundwater Flow with Gradually Increasing Heterogeneity due to Clogging

Well injection replenishes depleting water levels in a well field. Observation well water levels some distance away from the injection well are the indicators of the success of a well injection program. Simulation of the observation well response, located a few tens of meters from the injection well, is likely to be affected by the effects of nonhomogeneous medium, inclined initial water table, and aquifer clogging. Existing algorithms, such as the U.S. Geological Survey groundwater flow software MODFLOW, are capable of handling the first two conditions, whereas time-dependent clogging effects are yet to be introduced in the groundwater flow models. Elsewhere, aquifer clogging is extensively researched in theory of filtration; scope for its application in a well field is a potential research problem. In the present paper, coupling of one such filtration theory to MODFLOW is introduced. Simulation of clogging effects during “Hansol” well recharge in the parts of western India is found to be encouraging.

Numerical studies on channel formation and growth during solidification: Effect of process parameters

In the present work, solidification of a hyper-eutectic ammonium chloride solution in a bottom-cooled cavity (i.e. with stable thermal gradient) is numerically studied. A Rayleigh number based criterion is developed, which determines the conditions favorable for freckles formation. This criterion, when expressed in terms of physical properties and process parameters, yields the condition for plume formation as a function of concentration, liquid fraction, permeability, growth rate of a mushy layer and thermophysical properties. Subsequently, numerical simulations are performed for cases with initial and boundary conditions favoring freckle formation. The effects of parameters, such as cooling rate and initial concentration, on the formation and growth of freckles are investigated. It was found that a high cooling rate produced larger and more defined channels which are retained for a longer durations. Similarly, a lower initial concentration of solute resulted in fewer but more pronounced channels. The number and size of channels are also found to be related to the mushy zone thickness. The trends predicted with regard to the variation of number of channels with time under different process conditions are in accordance with the experimental observations reported in the literature.

Novel fluorescence detection technique for non-contact temperature sensing in microchip PCR

DNA amplification using Polymerase Chain Reaction (PCR) in a small volume is used in Lab-on-a-chip systems involving DNA manipulation. For few microliters of volume of liquid, it becomes difficult to measure and monitor the thermal profile accurately and reproducibly, which is an essential requirement for successful amplification. Conventional temperature sensors are either not biocompatible or too large and hence positioned away from the liquid leading to calibration errors. In this work we present a fluorescence based detection technique that is completely biocompatible and measures directly the liquid temperature. PCR is demonstrated in a 3 ILL silicon-glass microfabricated device using non-contact induction heating whose temperature is controlled using fluorescence feedback from SYBR green I dye molecules intercalated within sensor DNA. The performance is compared with temperature feedback using a thermocouple sensor. Melting curve followed by gel electrophoresis is used to confirm product specificity after the PCR cycles. (c) 2007 Elsevier B.V. All rights reserved.

Optimization of the distance between source and substrate for device-grade SnS films grown by the thermal evaporation technique

Tin monosulfide (SnS) films with varying distance between the source and substrate (DSS) were prepared by the thermal evaporation technique at a temperature of 300 degrees C to investigate the effect of the DSS on the physical properties. The physical properties of the as-deposited films are strongly influenced by the variation of DSS. The thickness, Sn to S at.% ratio, grain size, and root mean square (rms) roughness of the films decreased with the increase of DSS. The films grown at DSS = 10 and 15 cm exhibited nearly single-crystalline nature with low electrical resistivity. From Hall-effect measurements, it is observed that the films grown at DSS <= 15 cm have p-type conduction and the films grown at higher distances have n-type conduction due to the variation of the Sn/S ratio. The films grown at DSS = 15 cm showed higher optical band gap of 1.36 eV as compared with the films grown at other distances. The effect of the DSS on the physical properties of SnS films is discussed and reported.

The partition of unity finite element method for elastic wave propagation in Reissner-Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four-noded element for the Reissner-indlin plate is derived in this paper, which is free of shear locking. Such a locking-free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency-dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright (C)006 John Wiley & Sons, Ltd.

Effect of process parameters on transport phenomena during solidification in the presence of electromagnetic stirring

In the present work, a numerical study is performed to predict the effect of process parameters on transport phenomena during solidification of aluminium alloy A356 in the presence of electromagnetic stirring. A set of single-phase governing equations of mass, momentum, energy and species conservation is used to represent the solidification process and the associated fluid flow, heat and mass transfer. In the model, the electromagnetic forces are incorporated using an analytical solution of Maxwell equation in the momentum conservation equations and the slurry rheology during solidification is represented using an experimentally determined variable viscosity function. Finally, the set of governing equations is solved for various process conditions using a pressure based finite volume technique, along with an enthalpy based phase change algorithm. In present work, the effect of stirring intensity and cooling rate are considered. It is found that increasing stirring intensity results in increase of slurry velocity and corresponding increase in the fraction of solid in the slurry. In addition, the increasing stirring intensity results uniform distribution of species and fraction of solid in the slurry. It is also found from the simulation that the distribution of solid fraction and species is dependent on cooling rate conditions. At low cooling rate, the fragmentation of dendrites from the solid/liquid interface is more.

A Simple Instrumentation Calibration Technique for Electrical Impedance Tomography (EIT) Using A 16-Electrode Phantom

A simple analog instrumentation for Electrical Impedance Tomography is developed and calibrated using the practical phantoms. A constant current injector consisting of a modified Howland voltage controlled current source fed by a voltage controlled oscillator is developed to inject a constant current to the phantom boundary. An instrumentation amplifier, 50 Hz notch filter and a narrow band pass filter are developed and used for signal conditioning. Practical biological phantoms are developed and the forward problem is studied to calibrate the EIT-instrumentation. An array of sixteen stainless steel electrodes is developed and placed inside the phantom tank filled with KCl solution. 1 mA, 50 kHz sinusoidal current is injected at the phantom boundary using adjacent current injection protocol. The differential potentials developed at the voltage electrodes are measured for sixteen current injections. Differential voltage signal is passed through an instrumentation amplifier and a filtering block and measured by a digital multimeter. A forward solver is developed using Finite Element Method in MATLAB7.0 for solving the EIT governing equation. Differential potentials are numerically calculated using the forward solver with a simulated current and bathing solution conductivity. Measured potential data is compared with the differential potentials calculated for calibrating the instrumentation to acquire the voltage data suitable for better image reconstruction.

Effect of boundary condition description on convergence of solution in a boundary-value problem

It is well known that the numerical accuracy of a series solution to a boundary-value problem by the direct method depends on the technique of approximate satisfaction of the boundary conditions and on the stage of truncation of the series. On the other hand, it does not appear to be generally recognized that, when the boundary conditions can be described in alternative equivalent forms, the convergence of the solution is significantly affected by the actual form in which they are stated. The importance of the last aspect is studied for three different techniques of computing the deflections of simply supported regular polygonal plates under uniform pressure. It is also shown that it is sometimes possible to modify the technique of analysis to make the accuracy independent of the description of the boundary conditions.

Numerical procedure for second order non-linear ordinary differential equations and application to heat transfer problem

In this paper, we have first given a numerical procedure for the solution of second order non-linear ordinary differential equations of the type y″ = f (x;y, y′) with given initial conditions. The method is based on geometrical interpretation of the equation, which suggests a simple geometrical construction of the integral curve. We then translate this geometrical method to the numerical procedure adaptable to desk calculators and digital computers. We have studied the efficacy of this method with the help of an illustrative example with known exact solution. We have also compared it with Runge-Kutta method. We have then applied this method to a physical problem, namely, the study of the temperature distribution in a semi-infinite solid homogeneous medium for temperature-dependent conductivity coefficient.