Relationship between microscopic and macroscopic orientational relaxation times in polar liquids

Microscopic relations between single-particle orientational relaxation time (T, ) , dielectric relaxation time ( T ~ )a,n d many-body orientational relaxation time ( T ~o)f a dipolar liquid are derived. We show that both T~ and T~ are influenced significantly by many-body effects. In the present theory, these many-body effects enter through the anisotropic part of the two-particle direct correlation function of the polar liquid. We use mean-spherical approximation (MSA) for dipolar hard spheres for explicit numerical evaluation of the relaxation times. We find that, although the dipolar correlation function is biexponential, the frequency-dependent dielectric constant is of simple Debye form, with T~ equal to the transverse polarization relaxation time. The microscopic T~ falls in between Debye and Onsager-Glarum expressions at large values of the static dielectric constant.

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.

Chemical toughening of epoxies. II. Mechanical, thermal, and microscopic studies of epoxies toughened with hydroxyl-terminated poly(butadiene-co-acrylonitrile)

Diglycidyl ether–bisphenol-A-based epoxies toughened with various levels (0–12%) of chemically reacted liquid rubber, hydroxyl-terminated poly(butadiene-co-acrylonitrile) (HTBN) were studied for some of the mechanical and thermal properties. Although the ultimate tensile strength showed a continuous decrease with increasing rubber content, the toughness as measured by the area under the stress-vs.-strain curve and flexural strength reach a maximum around an optimum rubber concentration of 3% before decreasing. Tensile modulus was found to increase for concentrations below 6%. The glass transition temperature Tg as measured by DTA showed no variation for the toughened formulations. The TGA showed no variations in the pattern of decomposition. The weight losses for the toughened epoxies at elevated temperatures compare well with that of the neat epoxy. Scanning electron microscopy revealed the presence of a dual phase morphology with the spherical rubber particles precipitating out in the cured resin with diameter varying between 0.33 and 6.3 μm. In contrast, a physically blended rubber–epoxy showed much less effect towards toughening with the precipitated rubber particles of much bigger diameter (0.6–21.3 μm).

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.

Dielectric friction and solvation dynamics: Novel results on relaxation in dipolar liquids

In this article we present a new, general but simple, microscopic expression for time-dependent solvation energy of an ion. This expression is surprisingly similar to the expression for the time-dependent dielectric friction on a moving ion. We show that both the Chandra-Bagchi and the Fried-Mukamel formulations of solvation dynamics can be easily derived from this expression. This expression leads to an almost perfect agreement of the theory with all the available computer simulation results. Second, we show here for the first time that the mobility of a light solute ion can significantly accelerate its own solvation, specially in the underdamped limit. The latter result is also in excellent agreement with the computer simulations.

Microscopic theory of solvation of an ion in a binary dipolar liquid

A microscopic theory of the statics and the dynamics of solvation of an ion in a binary dipolar liquid is presented. The theory properly includes the different intermolecular correlations that are present in a binary mixture. As a result, the theory can explain several important aspects of both the statics and the dynamics of solvation that are observed in experiments. It provides a microscopic explanation of the preferential solvation of the more polar species by the solute ion. The dynamics of solvation is predicted to be highly non-exponential, in general. The average relaxation time is found to change nonlinearly with the composition of the mixture. These predictions are in qualitative agreement with the experimental results.

Molecular expression for dielectric friction on a rotating dipole: Reduction to the continuum theory

Recently we presented a microscopic expression for dielectric friction on a rotating dipole. This expression has a rather curious structure, involving the contributions of the transverse polarization modes of the solvent and also of the molecular length scale processes. It is shown here that under proper limiting conditions, this expression reduces exactly to the classical continuum model expression of Nee and Zwanzig [J. Chem. Phys. 52, 6353 (1970)]. The derivation requires the use of the asymptotic form of the orientation‐dependent total pair correlation function, the neglect of the contributions of translational modes of the solvent, and also the use of the limit that the size of the solvent molecules goes to zero. Thus, the derivation can be important in understanding the validity of the continuum model and can also help in explaining the results of a recent computer simulation study of dielectric relaxation in a Brownian dipolar lattice.

Molecular theory of underdamped dielectric relaxation: understanding collective effects in dipolar liquids

A molecular theory of underdamped dielectric relaxation of a dense dipolar liquid is presented. This theory properly takes into account the collective effects that are present (due to strong intermolecular correlations) in a dipolar liquid. For small rigid molecules, the theory again leads to a three-variable description which, however, is somewhat different from the traditional version. In particular, two of the three parameters are collective in nature and are determined by the orientational pair correlation function. A detailed comparison between the theory and the computer simulation results of Neria and Nitzan is performed and an excellent agreement is obtained without the use of any adjustable or free parameter - the calculation is fully microscopic. The theory can also provide a systematic description of the Poley absorption often observed in dipolar liquids in the high-frequency regime.

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.

A Microscopic Theory of the f.c.c.-b.c.c. Interface .

We present a first-principles theory of the equilibrium b.c.c.-f.c.c. interface at coexistence using the density functional method. We assume that the interfacial region has local body-centred tetragonal (b.c.t.) symmetry and predict typical interfacial widths to be of order 2 to 3 lattice spacings with typical energies close to 0.05 J/m2. These quantities are in good agreement with laboratory measurements on coherent interfaces.

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.