Sliding Modes for the Simulation of Mechanical and Electrical Systems Defined by Differential-Algebraic Equations

We offer a technique, motivated by feedback control and specifically sliding mode control, for the simulation of differential-algebraic equations (DAEs) that describe common engineering systems such as constrained multibody mechanical structures and electric networks. Our algorithm exploits the basic results from sliding mode control theory to establish a simulation environment that then requires only the most primitive of numerical solvers. We circumvent the most important requisite for the conventionalsimulation of DAEs: the calculation of a set of consistent initial conditions. Our algorithm, which relies on the enforcement and occurrence of sliding mode, will ensure that the algebraic equation is satisfied by the dynamic system even for inconsistent initial conditions and for all time thereafter. [DOI:10.1115/1.4001904]

Numerical simulation of solidification of liquid aluminum alloy flowing on cooling slope

Preparation of semisolid slurry using a cooling slope is increasingly becoming popular, primarily because of the simplicity in design and ease control of the process. In this process, liquid alloy is poured down an inclined surface which is cooled from underneath. The cooling enables partial solidification and the incline provides the necessary shear for producing semisolid slurry. However, the final microstructure of the ingot depends on several process parameters such as cooling rate, incline angle of the cooling slope, length of the slope and initial melt superheat. In this work, a CFD model using volume of fluid (VOF) method for simulating flow along the cooling slope was presented. Equations for conservation of mass, momentum, energy and species were solved to predict hydrodynamic and thermal behavior, in addition to predicting solid fraction distribution and macrosegregation. Solidification was modeled using an enthalpy approach and a volume averaged technique for the different phases. The mushy region was modeled as a multi-layered porous medium consisting of fixed columnar dendrites and mobile equiaxed/fragmented grains. The alloy chosen for the study was aluminum alloy A356, for which adequate experimental data were available in the literature. The effects of two key process parameters, namely the slope angle and the pouring temperature, on temperature distribution, velocity distribution and macrosegregation were also studied.

On the holographic simulation of non-circular zone plates

New methods involving the manipulation of fundamental wavefronts (e.g., plane and spherical) with simple optical components such as pinholes and spherical lenses have been developed for the fabrication of elliptic, hyperbolic and conical holographic zone plates. Also parabolic zone plates by holographic techniques have been obtained for the first time. The performance behaviour of these zone plates has been studied. Further a phenomenological explanation is offered for the observed improved fringe contrast obtained with a spherical reference wave.

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)].

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.

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.

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.

Digital simulation of galvanostatic current-potential data for gas-diffusion electrodes and estimation of electrodekinetic parameters

A computerized non-linear-least-squares regression procedure to analyse the galvanostatic current-potential data for kinetically hindered reactions on porous gas-diffusion electrodes is reported. The simulated data fit well with the corresponding measured values. The analytical estimates of electrode-kinetic parameters and uncompensated resistance are found to be in good agreement with their respective values obtained from Tafel plots and the current-interrupter method. The procedure circumvents the need to collect the data in the limiting-current region where the polarization values are usually prone to errors. The polarization data for two typical cases, namely, methanol oxidation on a carbon-supported platinum-tin electrode and oxygen reduction on a Nafion-coated platinized carbon electrode, are successfully analysed.

Brownian dynamics simulation of dense binary colloidal mixtures. I. Structural evolution and dynamics

We have carried out Brownian dynamics simulations of binary mixtures of charged colloidal suspensions of two different diameter particles with varying volume fractions phi and charged impurity concentrations n(i). For a given phi, the effective temperature is lowered in many steps by reducing n(i) to see how structure and dynamics evolve. The structural quantities studied are the partial and total pair distribution functions g(tau), the static structure factors, the time average g(<(tau)over bar>), and the Wendt-Abraham parameter. The dynamic quantity is the temporal evolution of the total meansquared displacement (MSD). All these parameters show that by lowering the effective temperature at phi = 0.2, liquid freezes into a body-centered-cubic crystal whereas at phi = 0.3, a glassy state is formed. The MSD at intermediate times shows significant subdiffusive behavior whose time span increases with a reduction in the effective temperature. The mean-squared displacements for the supercooled liquid with phi = 0.3 show staircase behavior indicating a strongly cooperative jump motion of the particles.

Brownian dynamics simulation of dense binary colloidal mixtures. II. Translational and bond-orientational order

We report the Brownian dynamics simulation results on the translational and bond-angle-orientational correlations for charged colloidal binary suspensions as the interparticle interactions are increased to form a crystalline (for a volume fraction phi = 0.2) or a glassy (phi = 0.3) state. The translational order is quantified in terms of the two- and four-point density autocorrelation functions whose comparisons show that there is no growing correlation length near the glass transition. The nearest-neighbor orientational order is determined in terms of the quadratic rotational invariant Q(l) and the bond-orientational correlation functions g(l)(t). The l dependence of Q(l) indicates that icosahedral (l = 6) order predominates at the cost of the cubic order (l = 4) near the glass as well as the crystal transition. The density and orientational correlation functions for a supercooled liquid freezing towards a glass fit well to the streched-exponential form exp[-(t/tau)(beta)]. The average relaxation times extracted from the fitted stretched-exponential functions as a function of effective temperatures T* obey the Arrhenius law for liquids freezing to a crystal whereas these obey the Vogel-Tamman-Fulcher law exp[AT(0)*/(T* - T-0*)] for supercooled Liquids tending towards a glassy state. The value of the parameter A suggests that the colloidal suspensions are ''fragile'' glass formers like the organic and molecular liquids.

Monte Carlo simulation of low temperature phase diagrams of YBa2Cu3O6+x

We report the results of Monte Carlo simulation of oxygen ordering in the oxygen deficient portion (x<0.5) of YBa2Cu3O6+x at low temperatures. We find qualitative agreement among cluster - variation, Monte Carlo and transfer matrix methods. However, low temperature and ground state simulations clearly indicate the presence of a tetragonal phase. There is also evidence for two second order phase transition lines separating the tetragonal and the �double cell� ortho II phase. The effect of decreasing the inter-chain repulsion on oxygen ordering has also been investigated.

Monte Carlo simulation of network formation based on structural fragments in epoxy-anhydride systems

A method combining the Monte Carlo technique and the simple fragment approach has been developed for simulating network formation in amine-catalysed epoxy-anhydride systems. The method affords a detailed insight into the nature and composition of the network, showing the distribution of various fragments. It has been used to characterize the network formation in the reaction of the diglycidyl ester of isophthalic acid with hexahydrophthalic anhydride, catalysed by benzyldimethylamine. Pre-gel properties like number and weight distributions and average molecular weights have been calculated as a function of epoxy conversion, leading to a prediction of the gel-point conversion. Analysis of the simulated network further yields other characteristic properties such as concentration of crosslink points, distribution and concentration of elastically active chains, average molecular weight between crosslinks, sol content and mass fraction of pendent chains. A comparison has been made of the properties obtained through simulation with those predicted by the fragment approach alone, which, however, gives only average properties. The Monte Carlo simulation results clearly show that loops and other cyclic structures occur in the gel. This may account for the differences observed between the results of the simulation and the fragment model in the post-gel phase. Copyright (C) 1996 Elsevier Science Ltd.