Water-modeling study of the surface disturbances in continuous slab caster

The present work is based on four static molds using nozzles of different port diameter, port angle, and immersion depth. It has been observed that the meniscus is wavy. The wave amplitude shows a parabolic variation with the nozzle exit velocity. The dimensionless amplitude is found to vary linearly with the Froude number. Vortex formation and bubble entrainment by the wave occurs at the meniscus beyond a critical flow rate, depending upon the nozzle configuration, immersion depth, and the mold aspect ratio.

Boson Inverse Operators and a New Family of Two-photon Annihilation Operators

We introduce the inverse of the Hermitian operator (acircacirc†) and express the Boson inverse operators acirc-1 and acirc†-1 in terms of the operators acirc, acirc† and (acircacirc†)-1. We show that these Boson inverse operators may be realized by Susskind-Glogower phase operators. In this way, we find a new two-photon annihilation operator and denote it as acirc2(acircacirc†)-1. We show that the eigenstates of this operator have interesting non-classical properties. We find that the eigenstates of the operators (acircacirc†)-1 acirc2, acirc(acircacirc†)-1 acirc and acirc2(acircacirc†)-1 have many similar properties and thus they constitute a family of two-photon annihilation operators.

Modeling Of Hydrolytic Polymerization In A Semibatch Nylon 6 Reactor

Hydrolytic polymerization of caprolactam to Nylon 6 in a semibatch reactor is carried out by heating a mixture of water and caprolactam. Evaporation of volatiles caused by heating results in a pressure build-up. After the pressure reaches a predetermined value, vapors are vented to keep the pressure constant for some time, and thereafter, to lower the pressure to a value slightly above atmospheric in a preprogrammed manner. The characteristics of the polymer are determined by the chemical reactions and the vaporization of water and caprolactam. The semibatch operation has been simulated and the predictions have been compared with industria data. The observed temperature and pressure histories were predicted with a fair degree of accuracy. It was found that the predictions of the degree of polymerization however are sensitive to the vapor-liquid equilibrium relations. A comparison with an earlier model, which neglected mass transfer resistance, indicates that simulation using the VLE data of Giori and Hayes and accounting for mass transfer resistance is more reliable.

Modes of binding of 2'-AMP to RNase T1. A computer modeling study

The modes of binding of adenosine 2'-monophosphate (2'-AMP) to the enzyme ribonuclease (RNase) T1 were determined by computer modelling studies. The phosphate moiety of 2'-AMP binds at the primary phosphate binding site. However, adenine can occupy two distinct sites--(1) The primary base binding site where the guanine of 2'-GMP binds and (2) The subsite close to the N1 subsite for the base on the 3'-side of guanine in a guanyl dinucleotide. The minimum energy conformers corresponding to the two modes of binding of 2'-AMP to RNase T1 were found to be of nearly the same energy implying that in solution 2'-AMP binds to the enzyme in both modes. The conformation of the inhibitor and the predicted hydrogen bonding scheme for the RNase T1-2'-AMP complex in the second binding mode (S) agrees well with the reported x-ray crystallographic study. The existence of the first mode of binding explains the experimental observations that RNase T1 catalyses the hydrolysis of phosphodiester bonds adjacent to adenosine at high enzyme concentrations. A comparison of the interactions of 2'-AMP and 2'-GMP with RNase T1 reveals that Glu58 and Asn98 at the phosphate binding site and Glu46 at the base binding site preferentially stabilise the enzyme-2'-GMP complex.

Neural network modeling of associative memory: Beyond the Hopfield model

A number of neural network models, in which fixed-point and limit-cycle attractors of the underlying dynamics are used to store and associatively recall information, are described. In the first class of models, a hierarchical structure is used to store an exponentially large number of strongly correlated memories. The second class of models uses limit cycles to store and retrieve individual memories. A neurobiologically plausible network that generates low-amplitude periodic variations of activity, similar to the oscillations observed in electroencephalographic recordings, is also described. Results obtained from analytic and numerical studies of the properties of these networks are discussed.

DockYard-a repository to assist modeling of protein-protein docking

In the absence of interlogs, building docking models is a time intensive task, involving generation of a large pool of docking decoys followed by refinement and screening to identify near native docking solutions. This limits the researcher interested in building docking methods with the choice of benchmarking only a limited number of protein complexes. We have created a repository called dockYard (http://pallab.serc.iisc.ernet.in/dockYard), that allows modelers interested in protein-protein interaction to access large volume of information on protein dimers and their interlogs, and also download decoys for their work if they are interested in building modeling methods. dockYard currently offers four categories of docking decoys derived from: Bound (native dimer co-crystallized), Unbound (individual subunits are crystallized, as well as the target dimer), Variants (match the previous two categories in at least one subunit with 100% sequence identity), and Interlogs (match the previous categories in at least one subunit with >= 90% or >= 50% sequence identity). The web service offers options for full or selective download based on search parameters. Our portal also serves as a repository to modelers who may want to share their decoy sets with the community.

Modeling of a reacting nanofilm on a composite substrate

This article provides a detailed computational analysis of the reaction of dense nanofilms and the heat transfer characteristics on a composite substrate. Although traditional energetic compounds based on organic materials have similar energy per unit weight, non-organic material in nanofilm configuration offers much higher energy density and higher flame speed. The reaction of a multilayer thin film of aluminum and copper oxide has been studied by varying the substrate material and thicknesses. The numerical analysis of the thermal transport of the reacting film deposited on the substrate combined a hybrid approach in which a traditional two-dimensional black box theory was used in conjunction with the sandwich model to estimate the appropriate heat flux on the substrate accounting for the heat loss to the surroundings. A procedure to estimate this heat flux using stoichiometric calculations is provided. This work highlights two important findings. One is that there is very little difference in the temperature profiles between a single substrate of silica and a composite substrate of silicon silica. Secondly, with increase in substrate thickness, the quenching effect is progressively diminished at a given speed. These findings show that the composite substrate is effective and that the average speed and quenching of flames depend on the thickness of the silica substrate, and can be controlled by a careful choice of the substrate configuration. (C) 2011 Elsevier Ltd. All rights reserved.

The structures of bacteriorhodopsin with different retinal-Schiff base orientations--computer modeling and energy minimization studies

Bacteriorhodopsin has been the subject of intense study in order to understand its photochemical function. The recent atomic model proposed by Henderson and coworkers based on electron cryo-microscopic studies has helped in understanding many of the structural and functional aspects of bacteriorhodopsin. However, the accuracy of the positions of the side chains is not very high since the model is based on low-resolution data. In this study, we have minimized the energy of this structure of bacteriorhodopsin and analyzed various types of interactions such as - intrahelical and interhelical hydrogen bonds and retinal environment. In order to understand the photochemical action, it is necessary to obtain information on the structures adopted at the intermediate states. In this direction, we have generated some intermediate structures taking into account certain experimental data, by computer modeling studies. Various isomers of retinal with 13-cis and/or 15-cis conformations and all possible staggered orientations of Lys-216 side chain were generated. The resultant structures were examined for the distance between Lys-216-schiff base nitrogen and the carboxylate oxygen atoms of Asp-96 - a residue which is known to reprotonate the schiff base at later stages of photocycle. Some of the structures were selected on the basis of suitable retinal orientation and the stability of these structures were tested by energy minimization studies. Further, the minimized structures are analyzed for the hydrogen bond interactions and retinal environment and the results are compared with those of the minimized rest state structure. The importance of functional groups in stabilizing the structure of bacteriorhodopsin and in participating dynamically during the photocycle have been discussed.

Water Modeling Study Of The Jet Characteristics In A Continuous-Casting Mold

The jet characteristics and the fluid flow pattern in a continuous slab caster have been studied using a water model. The fluid jet is studied under free fall and submerged discharge conditions. In the latter case, the jet was followed by dye-injection technique and image analyser was used to find out the effect of nozzle parameters on jet-spread angle, jet-discharge angle and the volume entrainment by the jet. All free-fall jets with nozzle port angle zero and upward are found to be spinning. Some of the free-fall jets with downward nozzle-port angle are found to be spinning and rest are smooth. The spinning direction of the jets are found to change with time. The well depth, port diameter and the inner diameter of the nozzle have a clear effect on the free-fall jets with downward port angle. The jet-spread angle is found to be about 17-degrees for smooth jets. The spread angle for spinning jet increases as the nozzle-port angle is increased from downward 25 to upward 15-degrees. The jet-discharge angle is always downward even when the nozzle-discharge ports are angled upward. The extent of volume entrainment by the spinning jet is higher and it increases as the nozzle-port angle is increased from 25 downward to 15-degrees upward.

Improving the convergence rate of the inverse iteration method

Solution of generalized eigenproblem, K phi = lambda M phi, by the classical inverse iteration method exhibits slow convergence for some eigenproblems. In this paper, a modified inverse iteration algorithm is presented for improving the convergence rate. At every iteration, an optimal linear combination of the latest and the preceding iteration vectors is used as the input vector for the next iteration. The effectiveness of the proposed algorithm is demonstrated for three typical eigenproblems, i.e. eigenproblems with distinct, close and repeated eigenvalues. The algorithm yields 29, 96 and 23% savings in computational time, respectively, for these problems. The algorithm is simple and easy to implement, and this renders the algorithm even more attractive.

Modeling of angiogenin-3 '-NMP complex

Angiogenin belongs to the Ribonuclease superfamily and has a weak enzymatic activity that is crucial for its biological function of stimulating blood vessel growth. Structural studies on ligand bound Angiogenin will go a long way in understanding the mechanism of the protein as well as help in designing drugs against it. In this study we present the first available structure of nucleotide ligand bound Angiogenin obtained by computer modeling. The importance of this study in itself notwithstanding, is a precursor to modeling a full dinucleotide substrate onto Angiogenin. Bovine Angiogenin, the structure of which has been solved at a high resolution, was earlier subjected to Molecular Dynamics simulations for a nanosecond. The MD structures offer better starting points for docking as they offer lesser obstruction than the crystal structure to ligand binding. The MD structure with the least serious short contacts was modeled to obtain a steric free Angiogenin - 3' mononucleotide complex structure. The structures were energetically minimized and subjected to a brief spell of Molecular Dynamics. The results of the simulation show that all the li,ligand-Angiogenin interactions and hydrogen bonds are retained, redeeming the structure and docking procedure. Further, following ligand - protein interactions in the case of the ligands 3'-CMP and 3'-UMP we were able to speculate on how Angiogenin, a predominantly prymidine specific ribonuclease prefers Cytosine to Uracil in the first base position.

Inverse surface melting in confined clusters: Ar-13 in zeolite L

Monte Carlo and molecular dynamics simulations on an Ar-13 cluster in zeolite L have been carried out at a series of temperatures to understand the rigid-nonrigid transition corresponding to the solid-liquid transition exhibited by the free Ar-13 cluster. The icosahedral geometry of the free cluster is no longer preferred when the cluster is confined in the zeolite. The root-mean-squared pair distance fluctuation, delta, exhibits a sharp, well-defined rigid-nonrigid transition at 17 K as compared to 27 K for the free cluster. Multiple peaks in the distribution of short-time averages of the guest-host interaction energy indicate coexistence of two phases.; It is shown that this transition is associated with the inner atoms becoming mobile at 17 K even while the outer layer atoms, which are in close proximity to the zeolitic wall, continue to be comparatively immobile. This may be contrasted with the melting of large free clusters of 40 or more atoms which exhibit surface melting. Guest-host interactions seem to play a predominant role in determining the properties of confined clusters. We demonstrate that the volume of the cluster increases rather sharply at 17 and 27 K respectively for the confined and the free cluster. Power spectra suggest that the motion of the inner atoms is generally parallel to the atoms which form the cage wall.

Modeling of precipitation in reverse micellar systems

A model of the precipitation process in reverse micelles has been developed to calculate the size of fine particles obtained therein. While the method shares several features of particle nucleation and growth common to precipitation in large systems, complexities arise in describing the processes of nucleation, due to the extremely small size of a micelle and of particle growth caused by fusion among the micelles. Occupancy of micelles by solubilized molecules is governed by Poisson statistics, implying most of them are empty and cannot nucleate of its own. The model therefore specifies the minimum number of solubilized molecules required to form a nucleus which is used to calculate the homogeneous nucleation rate. Simultaneously, interaction between micelles is assumed to occur by Brownian collision and instantaneous fusion. Analysis of time scales of various events shows growth of particles to be very fast compared to other phenomena occurring. This implies that nonempty micelles either are supersaturated or contain a single precipitated particle and allows application of deterministic population balance equations to describe the evolution of the system with time. The model successfully predicts the experimental measurements of Kandori ct al.(3) on the size of precipitated CaCO3 particles, obtained by carbonation of reverse micelles containing aqueous Ca(OH)(2) solution.

Modeling of the cooling process on the runout table of a hot strip mill - A parallel approach

This paper deals with the development of a new model for the cooling process on the runout table of hot strip mills, The suitability of different numerical methods for the solution of the proposed model equation from the point of view of accuracy and computation time are studied, Parallel solutions for the model equation are proposed.