A Numerical-Simulation Of The Gneiss-Charnockite Transformations In Southern India

This paper presents a numerical simulation of the well-documented, fluid-controlled Kabbal and Ponmudi type gneiss-chamockite transformations in southern India using a free energy minimization method. The computations have considered all the major solid phases and important fluid species in the rock - C-O-H and rock - C-O-H-N systems. Appropriate activity-composition relations for the solid solutions and equations of state for the fluids have been included in order to evaluate the mineral-fluid equilibria attending the incipient chamockite development in the gneisses. The C-O-H fluid speciation pattern in both the Kabbal and Ponmudi type systems indicates that CO2 and H2O make up the bulk of the fluid phase with CO, CH4, H-2 and O2 as minor constituents. In the graphite-buffered Ponmudi-system, the abundance of CO, CH4 and H-2 is orders of magnitude higher than that in the graphite-free Kabbal system. Simulation with C-O-H-N fluids of varying composition demonstrates the complementary role of CO2 and N2 as rather inert dilutants of H2O in the fluid phase. The simulation, carried out on available whole-rock data, has demonstrated the dependence of the transformation X(H2O) on P,T, and phase and chemical composition of the precursor gneiss.

Cell-dynamical simulation of magnetic hysteresis in the two-dimensional Ising system

We present results from numerical simulations using a ‘‘cell-dynamical system’’ to obtain solutions to the time-dependent Ginzburg-Landau equation for a scalar, two-dimensional (2D), (Φ2)2 model in the presence of a sinusoidal external magnetic field. Our results confirm a recent scaling law proposed by Rao, Krishnamurthy, and Pandit [Phys. Rev. B 42, 856 (1990)], and are also in excellent agreement with recent Monte Carlo simulations of hysteretic behavior of 2D Ising spins by Lo and Pelcovits [Phys. Rev. A 42, 7471 (1990)].

Influence of Die Surface Textures during Metal Forming-A Study Using Experiments and Simulation

Friction plays an important role in metal forming processes, and the surface texture of the die is a major factor that influences friction. In the present investigation, experiments were conducted to understand the role of surface texture of the harder die surface and load on coefficient of friction. The data analysis showed that the coefficient of friction is highly dependent on the surface texture of the die surface. Assigning different magnitude of coefficients of friction, obtained in the experiments, at different regions between the die and the workpiece, Finite element (FE) simulation of a compression test was carried out to understand the effect of friction on deformation and stress/strain-rate distribution. Simulation results revealed that, owing to the difference in coefficient of friction, there is a change in metal flow pattern. Both experimental and simulation results confirmed that the surface texture of the die surface and thus coefficient of friction directly affects the strain rate and flow pattern of the workpiece.

Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (epsilon) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.

Crystal structures, spectral and magnetic properties of (µ-hydroxo)(µ-acetato)dicopper(II) complexes containing chelating amines

The reaction of [Cu2(O2CMe)4(H2O)2] with N, N, N′, N′-tetramethylethane- 1,2-diamine (tmen) in ethanol yielded the dicopper(II) complex [Cu2(OH)(O2CMe)(tmen)2][ClO4]21. A similar reaction with N, N- dimethylethane- 1,2-diamine (dmen) afforded a crystalline product 2 in which two dicopper(II) complexes, [Cu2(OH)(O2CMe)(dmen)2][ClO4]22a and [Cu2(OH)(O2CMe)(H2O)2(dmen)2][ClO4]22b, are cocrystallized in a 1 : 1 molar ratio along with 2NaClO4. The crystal structures of 1 and 2 have been determined. The complexes have an asymmetrically dibridged [Cu2(µ-OH)(µ-O2CMe)]2+ core. The co-ordination geometry of the metal is square planar (CuO2N2). The copper atoms in 2b have a square-pyramidal CuO3N2 co-ordination sphere. The Cu Cu distances and Cu–O–Cu angles in 1, 2a and 2b are 3.339(2), 3.368(3), 3.395(7)Å, 120.1(2), 116.4(1) and 123.6(2)°, respectively. Complex 1 exhibits an axial ESR spectrum in a methanol glass giving g∥= 2.26 (A∥= 164 × 10–4 cm–1) and g⊥= 2.04. The ESR spectra obtained from the bulk material of the dmen product are indicative of the presence of two dimers, viz. complex 2a(g∥= 2.25, A∥= 165 × 10–4 cm–1; g⊥= 2.03) and 2b(g∥= 2.19, A∥= 184 × 10–4 cm–1; g⊥= 2.0). Variable-temperature magnetic susceptibility measurements on these complexes show an intramolecular antiferromagnetic coupling in the dimeric core. The fitting parameters are J=–27.8 cm–1, g= 2.1 for complex 1 and J=–10.1 cm–1, g= 2.0 for 2. The magnetostructural properties of the complexes are discussed. There is a linear correlation of the –2J values with the Cu Cu distances among dibridged complexes having square-planar copper(II) centres.

Raman spectral studies of borate glasses

Raman spectroscopic measurements in borate glasses have been reviewe. The review shows that the technique is useful in identifying the structural groups present in the borate on the basis of the Krogh-Moe hypothesis. Vitreous B2O3 and alkali borates are extensvvely studied and a satisfactory assignment of bands is possible by a careful consideration of the literature. A cation effect on the borate netwoork is observed. Availaable measurements on binary borates other than alkali borates and on ternary borates are limited and more work is required to identify the structural modifications that take place with composition. Mixed alkali effect is reported only lithium-caesium borade and shows the formation of non-bridging oxygens, destroying the six-membered rings when Li2O is replaced by Cs2O. Fast ionic glasses (alkali borates containing alkali halides) yield the same Raman spectra as the alkali borates, except when the alkali is a fluoride.

Solvation dynamics in a Brownian dipolar lattice. Comparison between computer simulation and various molecular theories of solvation dynamics

Several recent theoretical and computer simulation studies have considered solvation dynamics in a Brownian dipolar lattice which provides a simple model solvent for which detailed calculations can be carried out. In this article a fully microscopic calculation of the solvation dynamics of an ion in a Brownian dipolar lattice is presented. The calculation is based on the non‐Markovian molecular hydrodynamic theory developed recently. The main assumption of the present calculation is that the two‐particle orientational correlation functions of the solid can be replaced by those of the liquid state. It is shown that such a calculation provides an excellent agreement with the computer simulation results. More importantly, the present calculations clearly demonstrate that the frequency‐dependent dielectric friction plays an important role in the long time decay of the solvation time correlation function. We also find that the present calculation provides somewhat better agreement than either the dynamic mean spherical approximation (DMSA) or the Fried–Mukamel theory which use the simulated frequency‐dependent dielectric function. It is found that the dissipative kernels used in the molecular hydrodynamic approach and in the Fried–Mukamel theory are vastly different, especially at short times. However, in spite of this disagreement, the two theories still lead to comparable results in good agreement with computer simulation, which suggests that even a semiquantitatively accurate dissipative kernel may be sufficient to obtain a reliable solvation time correlation function. A new wave vector and frequency‐dependent dissipative kernel (or memory function) is proposed which correctly goes over to the appropriate expressions in both the single particle and the collective limits. This form is expected to lead to better results than all the existing descriptions.

Spectral approach for the soliton and periodic solutions of the nonlinear wave equation

A spectral method that obtains the soliton and periodic solutions to the nonlinear wave equation is presented. The results show that the nonlinear group velocity is a function of the frequency shift as well as of the soliton power. When the frequency shift is a function of time, a solution in terms of the Jacobian elliptic function is obtained. This solution is periodic in nature, and, to generate such an optical pulse train, one must simultaneously amplitude- and frequency-modulate the optical carrier. Finally, we extend the method to include the effect of self-steepening.

Synthesis, spectral characterisation and single-crystal structure of a bicyclic tetraphosphapentazane tetraoxide, N5P4ET5O4(OC6H3Me2-2,6)2

Reaction of the bicyclic phosphazane N5P4Et5Cl2 with 2,6-dimethylphenol and subsequent oxidation of the product by aqueous hydrogen peroxide yields N5P4Et5O4(OC6H3Me2-2,6)2 in 85% yield. Its structure has been established by NMR spectroscopy and single-crystal X-ray diffraction. The compound crystallises in the monoclinic space group C2/c with a= 21.245(5), b= 10.879(2), c= 16.450(6)Å, ?= 123.94(2)°, Z= 4, R= 0.066. The structural features are compared with those of bicyclic ?5-phosphazenes of type N5P4R3(NR1R2)5(NHR3)(R1,R3= Me or Et, R2= H or Me). The observed conformation of the N3P3 rings in the present compound is mainly dictated by the maximisation of the stabilising influence of �negative hyperconjugative interactions� between the nitrogen lone pairs and the adjacent P�X ?* orbitals.

A simple method for potential flow simulation of cascades

A simple method using a combination of conformal mapping and vortex panel method to simulate potential flow in cascades is presented. The cascade is first transformed to a single body using a conformal mapping, and the potential flow over this body is solved using a simple higher order vortex panel method. The advantage of this method over existing methodologies is that it enables the use of higher order panel methods, as are used to solve flow past an isolated airfoil, to solve the cascade problem without the need for any numerical integrations or iterations. The fluid loading on the blades, such as the normal force and pitching moment, may be easily calculated from the resultant velocity field. The coefficient of pressure on cascade blades calculated with this methodology shows good agreement with previous numerical and experimental results.

Analysis of Hydrogen Bonding and Stability of Protein Secondary Structures in Molecular Dynamics Simulation

Molecular Dynamics (MD) simulations provide an atomic level account of the molecular motions and have proven to be immensely useful in the investigation of the dynamical structure of proteins. Once an MD trajectory is obtained, specific interactions at the molecular level can be directly studied by setting up appropriate combinations of distance and angle monitors. However, if a study of the dynamical behavior of secondary structures in proteins becomes important, this approach can become unwieldy. We present herein a method to study the dynamical stability of secondary structures in proteins, based on a relatively simple analysis of backbone hydrogen bonds. The method was developed for studying the thermal unfolding of beta-lactamases, but can be extended to other systems and adapted to study relevant properties.

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.

Stochastic Approximation Algorithms for Constrained Optimization via Simulation

We develop four algorithms for simulation-based optimization under multiple inequality constraints. Both the cost and the constraint functions are considered to be long-run averages of certain state-dependent single-stage functions. We pose the problem in the simulation optimization framework by using the Lagrange multiplier method. Two of our algorithms estimate only the gradient of the Lagrangian, while the other two estimate both the gradient and the Hessian of it. In the process, we also develop various new estimators for the gradient and Hessian. All our algorithms use two simulations each. Two of these algorithms are based on the smoothed functional (SF) technique, while the other two are based on the simultaneous perturbation stochastic approximation (SPSA) method. We prove the convergence of our algorithms and show numerical experiments on a setting involving an open Jackson network. The Newton-based SF algorithm is seen to show the best overall performance.

Autoregressive spectral array for graphical display of EEG data

A graphical display of the frequency content of background,electroencephalogram (EEG) activity is obtained by calculating the spectral estimates using autocorrelation autoregressive method and the classical Fourier transform method, Display of spectral content of consecutive data segments is made using hidden-line suppression technique so as to get a spectral array, The autoregressive spectral array (ASA) is found to be sensitive to baseline drift, Following baseline correction the autoregressive technique is found to be superior to the Fourier method of compressed spectral array (CSA) in detecting the transitions in the frequencies of the signal. The smoothed ASA gives a better picture of transitions and changes in the background activity, The ASA can be made to adapt to specific changes of dominant frequencies while eliminating unnecessary peaks in the spectrum. The utility,of the ASA for background EEG analysis is discussed,

High Yield Room Temperature Synthesis And Spectral Studies Of Tri(Amino)Silanes: (R2n)3SIH

Tri(amino)silanes were prepared by the condensation of trichlorosilane with secondary amines in 1:6 molar ratio. Reactions of trichlorosilane with pyrrolidine, piperidine, hexamethyleneimine, morpholine, N-methylpiperazine and diethylamine afford the tri(amino)silanes in nearly quantitative yields. Their physical and spectroscopic properties are discussed. All these compounds are highly sensitive to moisture and hydrolyse to silica and the respective amine with the evolution of hydrogen. The compounds have been characterised by IR, 1H NMR, [1H]29Si NMR spectroscopic methods and CHN elemental analysis.