Effect of capillarity and heterogeneity in the numerical modelling of multiphase flow of fluids in unsaturated porous medium

The multiphase flow of fluids in the unsaturated porous medium is considered as a three phase flow of water, NAPL, and air simultaneously in the porous medium. The adaptive solution fully implicit modified sequential method is used for the numerical modelling. The effect of capillarity and heterogeneity effect at the interface between the media is studied and it is observed that the interface criteria has to be taken into account for the correct prediction of NAPL migration especially in heterogeneous media. The modified Newton Raphson method is used for the linearization and Hestines and Steifel Conjugate Gradient method is used as the solver.

Design, modelling and simulation of vibratory micromachined gyroscopes

Among various MEMS sensors, a rate gyroscope is one of the most complex sensors from the design point of view. The gyro normally consists of a proof mass suspended by an elaborate assembly of beams that allow the system to vibrate in two transverse modes. The structure is normally analysed and designed using commercial FEM packages such as ANSYS or MEMS specific commercial tools such as Coventor or Intellisuite. In either case, the complexity in analysis rises manyfolds when one considers the etch hole topography and the associated fluid flow calculation for damping. In most cases, the FEM analysis becomes prohibitive and one resorts to equivalent electrical circuit simulations using tools like SABER in Coventor. Here, we present a simplified lumped parameter model of the tuning fork gyro and show how easily it can be implemented using a generic tool like SIMULINK. The results obtained are compared with those obtained from more elaborate and intense simulations in Coventor. The comparison shows that lumped parameter SIMULINK model gives equally good results with fractional effort in modelling and computation. Next, the performance of a symmetric and decoupled vibratory gyroscope structure is also evaluated using this approach and a few modifications are made in this design to enhance the sensitivity of the device.

Numerical Modelling of the Observed Diurnal Cycle of Indian Monsoon Rainfall

We have analysed the diurnal cycle of rainfall over the Indian region (10S-35N, 60E-100E) using both satellite and in-situ data, and found many interesting features associated with this fundamental, yet under-explored, mode of variability. Since there is a distinct and strong diurnal mode of variability associated with the Indian summer monsoon rainfall, we evaluate the ability of the Weather Research and Forecasting Model (WRF) to simulate the observed diurnal rainfall characteristics. The model (at 54km grid-spacing) is integrated for the month of July, 2006, since this period was particularly favourable for the study of diurnal cycle. We first evaluate the sensitivity of the model to the prescribed sea surface temperature (SST), by using two different SST datasets, namely, Final Analyses (FNL) and Real-time Global (RTG). It was found that with RTG SST the rainfall simulation over central India (CI) was significantly better than that with FNL. On the other hand, over the Bay of Bengal (BoB), rainfall simulated with FNL was marginally better than with RTG. However, the overall performance of RTG SST was found to be better than FNL, and hence it was used for further model simulations. Next, we investigated the role of the convective parameterization scheme on the simulation of diurnal cycle of rainfall. We found that the Kain-Fritsch (KF) scheme performs significantly better than Betts-Miller-Janjić (BMJ) and Grell-Devenyi schemes. We also studied the impact of other physical parameterizations, namely, microphysics, boundary layer, land surface, and the radiation parameterization, on the simulation of diurnal cycle of rainfall, and identified the “best” model configuration. We used this configuration of the “best” model to perform a sensitivity study on the role of various convective components used in the KF scheme. In particular, we studied the role of convective downdrafts, convective timescale, and feedback fraction, on the simulated diurnal cycle of rainfall. The “best” model simulations, in general, show a good agreement with observations. Specifically, (i) Over CI, the simulated diurnal rainfall peak is at 1430 IST, in comparison to the observed 1430-1730 IST peak; (ii) Over Western Ghats and Burmese mountains, the model simulates a diurnal rainfall peak at 1430 IST, as opposed to the observed peak of 1430-1730 IST; (iii) Over Sumatra, both model and observations show a diurnal peak at 1730 IST; (iv) The observed southward propagating diurnal rainfall bands over BoB are weakly simulated by WRF. Besides the diurnal cycle of rainfall, the mean spatial pattern of total rainfall and its partitioning between the convective and stratiform components, are also well simulated. The “best” model configuration was used to conduct two nested simulations with one-way, three-level nesting (54-18-6km) over CI and BoB. While, the 54km and 18km simulations were conducted for the whole of July, 2006, the 6km simulation was carried out for the period 18 - 24 July, 2006. The results of our coarse- and fine-scale numerical simulations of the diurnal cycle of monsoon rainfall will be discussed.

Wavelet based spectral finite element modelling and detection of de-lamination in composite beams

In this paper, a model for composite beam with embedded de-lamination is developed using the wavelet based spectral finite element (WSFE) method particularly for damage detection using wave propagation analysis. The simulated responses are used as surrogate experimental results for the inverse problem of detection of damage using wavelet filtering. The WSFE technique is very similar to the fast fourier transform (FFT) based spectral finite element (FSFE) except that it uses compactly supported Daubechies scaling function approximation in time. Unlike FSFE formulation with periodicity assumption, the wavelet-based method allows imposition of initial values and thus is free from wrap around problems. This helps in analysis of finite length undamped structures, where the FSFE method fails to simulate accurate response. First, numerical experiments are performed to study the effect of de-lamination on the wave propagation characteristics. The responses are simulated for different de-lamination configurations for both broad-band and narrow-band excitations. Next, simulated responses are used for damage detection using wavelet analysis.

Modelling and analysis of an open-loop induction motor drive incorporating the effect of inverter dead-time

The objective of this paper is to study the influence of inverter dead-time on steady as well as dynamic operation of an open-loop induction motor drive fed from a voltage source inverter (VSI). Towards this goal, this paper presents a systematic derivation of a dynamic model for an inverter-fed induction motor, incorporating the effect of inverter dead-time, in the synchronously revolving dq reference frame. Simulation results based on this dynamic model bring out the impact of inverter dead-time on both the transient response and steady-state operation of the motor drive. For the purpose of steady-state analysis, the dynamic model of the motor drive is used to derive a steady-state model, which is found to be non-linear. The steady-state model shows that the impact of dead-time can be seen as an additional resistance in the stator circuit, whose value depends on the stator current. Towards precise evaluation of this dead-time equivalent resistance, an analytical expression is proposed for the same in terms of inverter dead-time, switching frequency, modulation index and load impedance. The notion of dead-time equivalent resistance is shown to simplify the solution of the non-linear steady-state model. The analytically evaluated steady-state solutions are validated through numerical simulations and experiments.

Increasing ionic conductivity and mechanical strength of a plastic electrolyte by inclusion of a polymer

in this contribution we present a soft matter solid electrolyte which was obtained by inclusion of a polymer (polyacrylonitrile, PAN) in LiClO4/LiTFSI-succinonitrile (SN), a semi-solid organic plastic electrolyte. Addition of the polymer resulted in considerable enhancement in ionic conductivity as well as mechanical strength of LiX-SN (X=ClO4, TFSI) plastic electrolyte. Ionic conductivity of 92.5%-[1 M LiClO4-SN]:7.5%-PAN (PAN amount as per SN weight) composite at 25 degrees C recorded a remarkably high value of 7 x 10(-3) Omega(-1) cm(-1), higher by few tens of order in magnitude compared to 1 M LiClO4-SN. Composite conductivity at sub-ambient temperature is also quite high. At -20 degrees C, the ionic conductivity of (100 -x)%-[1 M LiClO4-SN]:x%-PAN composites are in the range 3 x 10(-5)-4.5 x 10(-4) Omega(-1) cm(-1), approximately one to two orders of magnitude higher with respect to 1 M LiClO4-SN electrolyte conductivity. Addition of PAN resulted in an increase of the Young's modulus (Y) from Y -> 0 for LiClO4-SN to a maximum of 0.4MPa for the composites. Microstructural studies based on X-ray diffraction, differential scanning calorimetry and Fourier transform infrared spectroscopy suggest that enhancement in composite ionic conductivity is a combined effect of decrease in crystallinity and enhanced trans conformer concentration. (c) 2008 Elsevier Ltd. All rights reserved.

Synthesis, reactivity and structural characterisation of bicyclic tetraphosphapentazanes

The diphenoxy bicyclic tetraphosphapentazane derivatives (EtN)(5)P-4(OPh)(2) 2 and its monoxide (EtN)(5)P-4(O)(OPh)(2) 3 have been prepared. Both 2 and 3 exist as a mixture of two isomers. One isomer of (EtN)(5)P-4(O)(OPh)(2) 3a has been isolated and its reaction with tetrachloro-1,2-benzoquinone yielded (EtN)(5)P-4(O)(OPh)(2)(O2C6Cl4) 5 in which the junction phosphorus atom becomes five-co-ordinated. Treatment of 2 or 3a with [Mo(CO)(4)(nbd)] (nbd = norbornadiene, bicyclo[2.2.1]hepta-2,5-diene), on the other hand, yielded the chelate complex [Mo(CO)(4){(EtN)(5)P-4(O)(n)(OPh)(2)}] (n = 0 or 1; 6 or 7) in which the peripheral phosphorus atoms are bonded to the metal. The structures of 3a and 5-7 have been confirmed by single-crystal X-ray diffraction studies. The two P3N3 rings in 3a and 5 adopt twist/twist and irregular/twist conformations respectively; the phenoxy substituents occupy the 'pseudo axial' positions. However, an ideal chair conformation is observed for the P3N3 rings in 6 and 7 with the phenoxy substituents taking up the 'pseudo equatorial' positions. The NMR spectroscopic data for the compounds are discussed.

Microstructural evolution and mechanical characteristics in nanocrystalline nickel with a bimodal grain-size distribution

An attractive microstructural possibility for enhancing the ductility of high-strength nanocrystals is to develop a bimodal grain-size distribution, in which the fine grains provide strength, and the coarser grains enable strain hardening. Annealing of nanocrystalline Ni over a range of temperatures and times led to microstructures with varying volume fractions of coarse grains and a change in texture. Tensile tests revealed a drastic reduction in ductility with increasing volume fraction of coarse grains. The reduction in ductility may be related to the segregation of sulphur to grain boundaries.

Electrical and mechanical properties of diluted magnetic semiconductor Zn1-xMnxS nanocrystalline films

Nanostructured Zn1-xMnxS films (0 less-than-or-equals, slant x less-than-or-equals, slant 0.25) were deposited on glass substrates by simple resistive thermal evaporation technique. All the films were deposited at 300 K in a vacuum of 2*10-6 m bar. All the films temperature dependence of resistivity revealed semiconducting behaviour of the samples. Hot probe test revealed that all the samples exhibited n-type conductivity. The nanohardness of the films ranges from 4.7 to 9.9 GPa, Young's modulus value ranging 69.7-94.2 GPa.

X-ray and IR characterisation of Nasicon solid solution, Na1+xZr2SixP3−xO12

X-ray and IR studies on Nasicon solid solutions, Na1+xZr2SixP3−xO12, are carried out as a function of composition x. X-ray diffraction studies show that the unit cell volume increases as x increases and exhibits a maximum at x ≈ 2. On further increase in x the unit cell volume decreases. The infrared absorption peak positions and the splitting of these absorption peaks suggest a distortion of the PO4 and SiO4 tetrahedra. But the distortion is not large enough to change the local symmetry around the phosphorus or silicon ion from Td to C3v.

Effect of microstructure and mechanical properties on the seizure resistance of cast aluminium alloys

Seizure resistance of several cast aluminium base alloys has been examined using a standard Hohman Wear Tester. Disks of aluminium base alloys were run against a standard aluminium 12% silicon base alloy. The seizure resistance of the alloys (as measured by the lowest bearing parameter reached before seizure) increased with hardness, yield and tensile strength. In Al-Si-Ni alloys where silicon and nickel have little solid solubility in α-aluminium and Si and Ni Al3 hard phases are formed, the minimum bearing parameter decreased with the parameter V (The product of vol. % of hard phases in the disk and the shoe). Apparently the silicon and NiAl3 particles provided discontinuities in the matrix and reduced the probability (1 − V) of the α-aluminium phase in the disk coming into contact with the α-aluminium phase in the shoe. The copper and magnesium containing Al-Si-Ni alloys with lesser volumes of hard phases exhibit considerably better seizure resistance indicating that a slight increase in the solute content or the hardness of the primary α-phase leads to a considerable increase in seizure resistance. Deformation during wear and seizure leads to fragmentation of the original hard particles into considerably smaller particles uniformly dispersed in the deformed α-aluminium matrix.

Microstructure and mechanical properties of unidirectionally solidified Al-Si-Ni ternary eutectic

Al-10.98 pct Si-4.9 pct Ni ternary eutectic alloy was unidirectionally solidified at growth rates from 1.39μm/sec to 6.95μm/sec. Binary Al-Ni and Al-Si eutectics prepared from the same purity metals were also solidified under similar conditions to characterize the growth conditions under the conditions of present study. NiAl3 phase appeared as fibers in the binary Al-Ni eutectic and silicon appeared as irregular plates in the binary Al-Si eutectic. However, in the ternary Al-Si-Ni eutectic alloy both NiAl3 and silicon phases appeared as irregular plates dispersed in α-Al phase, without any regular repctitive arrangement. The size and spacing of NiAl3 and Si platelets in cone shaped colonies decreased with an increase in the growth rate of the ternary eutectic. Examination of specimen quenched during unidirectional solidification indicated that the ternary eutectic grows with a non-planar interface with both Si and NiAl3 phases protruding into the liquid. It is concluded that it will be difficult to grow regular ternary eutectic structures even if only one phase has a high entropy of melting. The tensile strength and modulus of unidirectionally solidified Al-Si-Ni eutectic was lower than the chill cast alloys of the same composition, and decreased with a decrease in growth rate. Tensile modulus and strength of ternary Al-Si-Ni eutectic alloys was greater than binary Al-Si eutectic alloy under similar growth conditions, both in the chill cast and in unidirectionally solidified conditions.