DNA Free Energy-Based Promoter Prediction and Comparative Analysis of Arabidopsis and Rice Genomes

The cis-regulatory regions on DNA serve as binding sites for proteins such as transcription factors and RNA polymerase. The combinatorial interaction of these proteins plays a crucial role in transcription initiation, which is an important point of control in the regulation of gene expression. We present here an analysis of the performance of an in silico method for predicting cis-regulatory regions in the plant genomes of Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) on the basis of free energy of DNA melting. For protein-coding genes, we achieve recall and precision of 96% and 42% for Arabidopsis and 97% and 31% for rice, respectively. For noncoding RNA genes, the program gives recall and precision of 94% and 75% for Arabidopsis and 95% and 90% for rice, respectively. Moreover, 96% of the false-positive predictions were located in noncoding regions of primary transcripts, out of which 20% were found in the first intron alone, indicating possible regulatory roles. The predictions for orthologous genes from the two genomes showed a good correlation with respect to prediction scores and promoter organization. Comparison of our results with an existing program for promoter prediction in plant genomes indicates that our method shows improved prediction capability.

Generalized model for infrared perception from an engine exhaust

A comprehensive scheme has been developed for the prediction of radiation from engine exhaust and its incidence on an arbitrarily located sensor. Existing codes have been modified for the simulation of flows inside nozzles and jets. A novel view factor computation scheme has been applied for the determination of the radiosities of the discrete panels of a diffuse and gray nozzle surface. The narrowband model has been used to model the radiation from the gas inside the nozzle and the nonhomogeneous jet. The gas radiation from the nozzle inclusive of nozzle surface radiosities have been used as boundary conditions on the jet radiation. Geometric modeling techniques have been developed to identify and isolate nozzle surface panels and gas columns of the nozzle and jet to determine the radiation signals incident on the sensor. The scheme has been validated for intensity and heat flux predictions, and some useful results of practical importance have been generated to establish its viability for infrared signature analysis of jets.

Modelling, simulation, and graphical user interface for industrial gas carburising process

A mathematical model has been developed for the gas carburising (diffusion) process using finite volume method. The computer simulation has been carried out for an industrial gas carburising process. The model's predictions are in good agreement with industrial experimental data and with data collected from the literature. A study of various mass transfer and diffusion coefficients has been carried out in order to suggest which correlations should be used for the gas carburising process. The model has been interfaced in a Windows environment using a graphical user interface. In this way, the model is extremely user friendly. The sensitivity analysis of various parameters such as initial carbon concentration in the specimen, carbon potential of the atmosphere, temperature of the process, etc. has been carried out using the model.

Nonlinear dynamics of eucaryotic pyruvate dehydrogenase multienzyme complex: Decarboxylation rate, oscillations, and multiplicity

Pyruvate conversion to acetyl-CoA by the pyruvate dehydrogenase (PDH) multienzyme complex is known as a key node in affecting the metabolic fluxes of animal cell culture. However, its possible role in causing possible nonlinear dynamic behavior such as oscillations and multiplicity of animal cells has received little attention. In this work, the kinetic and dynamic behavior of PDH of eucaryotic cells has been analyzed by using both in vitro and simplified in vivo models. With the in vitro model the overall reaction rate (v(1)) of PDH is shown to be a nonlinear function of pyruvate concentration, leading to oscillations under certain conditions. All enzyme components affect v, and the nonlinearity of PDH significantly, the protein X and the core enzyme dihydrolipoamide acyltransferase (E2) being mostly predominant. By considering the synthesis rates of pyruvate and PDH components the in vitro model is expanded to emulate in vivo conditions. Analysis using the in vivo model reveals another interesting kinetic feature of the PDH system, namely, multiple steady states. Depending on the pyruvate and enzyme levels or the operation mode, either a steady state with high pyruvate decarboxylation rate or a steady state with significantly lower decarboxylation rate can be achieved under otherwise identical conditions. In general, the more efficient steady state is associated with a lower pyruvate concentration. A possible time delay in the substrate supply and enzyme synthesis can also affect the steady state to be achieved and lead's to oscillations under certain conditions. Overall, the predictions of multiplicity for the PDH system agree qualitatively well with recent experimental observations in animal cell cultures. The model analysis gives some hints for improving pyruavte metabolism in animal cell culture.

Sandwich propellant combustion: Modeling and experimental comparison

This paper reports reacting fluid dynamics calculations for an ammonium percholrate binder sandwich and extracts experimentally observed features including surface profiles and maximum regression rates as a function of pressure and binder thickness. These studies have been carried out by solving the two-dimensional unsteady Navier-Stokes equations with energy and species conservation equations and a kinetic model of three reaction steps (ammonium perchlorate decomposition flame, primary diffusion flame, and final diffusion flame) in the gas phase. The unsteady two-dimensional conduction equation is solved in the condensed phase. The regressing surface is unsteady and two dimensional. Computations have been carried out for a binder thickness range of 25-125 mum and a pressure range of 1.4 to 6.9 MPa. Good comparisons at several levels of detail are used to demonstrate the need for condensed-phase two-dimensional unsteady conduction and three-step gas-phase reactions. The choice of kinetic and thermodynamic parameters is crucial to good comparison with experiments. The choice of activation energy parameters for ammonium percholrate combustion has been made with stability of combustion in addition to experimentally determined values reported in literature. The choice of gas-phase parameters for the diffusion flames are made considering that (a) primary diffusion flame affects the low-pressure behavior and (b) final diffusion flame affects high-pressure behavior. The predictions include the low-pressure deflagration limit of the sandwich apart from others noted above. Finally, this study demonstrates the possibility of making meaningful comparisons with experimental observations on sandwich propellant combustion.

Scaling analysis of momentum and heat transport in gas tungsten arc weld pools

A systematic procedure is outlined for scaling analysis of momentum and heat transfer in gas tungsten arc weld pools. With suitable selections of non-dimentionalised parameters, the governing equations coupled with appropriate boundary conditions are first scaled, and the relative significance of various terms appearing in them is analysed accordingly. The analysis is then used to predict the orders of magnitude of some important quantities, such as the velocity scene lit the top surface, velocity boundary layer thickness, maximum temperature increase in the pool, and time required for initiation of melting. Some of the quantities predicted from the scaling analysis can also be used for optimised selection of appropriate grid size and time steps for full numerical simulation of the process. The scaling predictions are finally assessed by comparison with numerical results quoted in the literature, and a good qualitative agreement is observed.

Scaling analysis of momentum, heat, and mass transfer in binary alloy solidification problems

A systematic approach is developed for scaling analysis of momentum, heat and species conservation equations pertaining to the case of solidification of a binary mixture. The problem formulation and description of boundary conditions are kept fairly general, so that a large class of problems can be addressed. Analysis of the momentum equations coupled with phase change considerations leads to the establishment of an advection velocity scale. Analysis of the energy equation leads to an estimation of the solid layer thickness. Different regimes corresponding to different dominant modes of transport are simultaneously identified. A comparative study involving several cases of possible thermal boundary conditions is also performed. Finally, a scaling analysis of the species conservation equation is carried out, revealing the effect of a non-equilibrium solidification model on solute segregation and species distribution. It is shown that non-equilibrium effects result in an enhanced macrosegregation compared with the case of an equilibrium model. For the sake of assessment of the scaling analysis, the predictions are validated against corresponding computational results.

Study of translational and rotational mobility and orientational preference of ethane in one-dimensional channels

Structural and dynamical properties of ethane in one-dimensional channels of AlPO4-5 and carbon nanotube have been investigated at dilute concentration with the help of molecular dynamics simulation. Density distributions and orientational structure of ethane have been analyzed. Repulsive interactions seem to play an important role when ethane is located in the narrow part of the AlPO4-5 channel. In AlPO4-5, parallel orientation is predominant over perpendicular orientation except when ethane is located in the broader part of the channel. Unlike in the case of single-file diffusion, our results in carbon nanotube show that at dilute concentrations the mean squared displacement, mu(2)(t) approximate to t(alpha), alpha = 1.8. The autocorrelation function for the z-component of angular velocity of ethane in space-fixed frame of reference shows a pronounced negative correlation. This is attributed to the restriction in the movement of ethane along the x- and y- directions. It is seen that the ratio of reorientational correlation times does not follow the Debye model for confined ethane but it is closer to the predictions of the Debye model for bulk ethane.

A scaling analysis of momentum and heat transport in laser surface melting

In this paper, we outline a systematic procedure for scaling analysis of momentum and heat transfer in laser melted pools. With suitable choices of non-dimensionalising parameters, the governing equations coupled with appropriate boundary conditions are first scaled, and the relative significance of various terms appearing in them are accordingly analysed. The analysis is then utilised to predict the orders of magnitude of some important quantities, such as the velocity scale at the top surface, velocity boundary layer thickness, maximum temperature rise in the pool, fully developed pool-depth, and time required for initiation of melting. Using the scaling predictions, the influence of various processing parameters on the system variables can be well recognised, which enables us to develop a deeper insight into the physical problem of interest. Moreover, some of the quantities predicted from the scaling analysis can be utilised for optimised selection of appropriate grid-size and time-steps for full numerical simulation of the process. The scaling predictions are finally assessed by comparison with experimental and numerical results quoted in the literature, and an excellent qualitative agreement is observed.

Grain Floatation During Equiaxed Solidification of an Al-Cu Alloy in a Side-Cooled Cavity: Part II-Numerical Studies

In this article, a single-phase, one-domain macroscopic model is developed for studying binary alloy solidification with moving equiaxed solid phase, along with the associated transport phenomena. In this model, issues such as thermosolutal convection, motion of solid phase relative to liquid and viscosity variations of the solid-liquid mixture with solid fraction in the mobile zone are taken into account. Using the model, the associated transport phenomena during solidification of Al-Cu alloys in a rectangular cavity are predicted. The results for temperature variation, segregation patterns, and eutectic fraction distribution are compared with data from in-house experiments. The model predictions compare well with the experimental results. To highlight the influence of solid phase movement on convection and final macrosegregation, the results of the current model are also compared with those obtained from the conventional solidification model with stationary solid phase. By including the independent movement of the solid phase into the fluid transport model, better predictions of macrosegregation, microstructure, and even shrinkage locations were obtained. Mechanical property prediction models based on microstructure will benefit from the improved accuracy of this model.

Roy equation analysis of ππ scattering

We analyse the Roy equations for the lowest partial waves of elastic ππ scattering. In the first part of the paper, we review the mathematical properties of these equations as well as their phenomenological applications. In particular, the experimental situation concerning the contributions from intermediate energies and the evaluation of the driving terms are discussed in detail. We then demonstrate that the two S-wave scattering lengths a00 and a02 are the essential parameters in the low energy region: Once these are known, the available experimental information determines the behaviour near threshold to within remarkably small uncertainties. An explicit numerical representation for the energy dependence of the S- and P-waves is given and it is shown that the threshold parameters of the D- and F-waves are also fixed very sharply in terms of a00 and a20. In agreement with earlier work, which is reviewed in some detail, we find that the Roy equations admit physically acceptable solutions only within a band of the (a00,a02) plane. We show that the data on the reactions e+e−→ππ and τ→ππν reduce the width of this band quite significantly. Furthermore, we discuss the relevance of the decay K→ππeν in restricting the allowed range of a00, preparing the grounds for an analysis of the forthcoming precision data on this decay and on pionic atoms. We expect these to reduce the uncertainties in the two basic low energy parameters very substantially, so that a meaningful test of the chiral perturbation theory predictions will become possible.

An improved Monte Carlo scheme for simulation of synthesis of nanoparticles in reverse micelles

An improved Monte Carlo technique is presented in this work to simulate nanoparticle formation through a micellar route. The technique builds on the simulation technique proposed by Bandyopadhyaya et al. (Langmuir 2000, 16, 7139) which is general and rigorous but at the same time very computation intensive, so much so that nanoparticle formation in low occupancy systems cannot be simulated in reasonable time. In view of this, several strategies, rationalized by simple mathematical analyses, are proposed to accelerate Monte Carlo simulations. These are elimination of infructuous events, removal of excess reactant postreaction, and use of smaller micelle population a large number of times. Infructuous events include collision of an empty micelle with another empty one or with another one containing only one molecule or only a solid particle. These strategies are incorporated in a new simulation technique which divides the entire micelle population in four classes and shifts micelles from one class to other as the simulation proceeds. The simulation results, throughly tested using chi-square and other tests, show that the predictions of the improved technique remain unchanged, but with more than an order of magnitude decrease in computational effort for some of the simulations reported in the literature. A post priori validation scheme for the correctness of the simulation results has been utilized to propose a new simulation strategy to arrive at converged simulation results with near minimum computational effort.

Turbulent non-equlibrium wakes

We consider here the detailed application of a model Reynolds stress equation (Narasimha 1969) to plane turbulent wakes subjected to pressure gradients. The model, which is a transport equation for the stress exhibiting relaxation and diffusion, is found to be consistent with the observed response of a wake to a nearly impulsive pressure gradient (Narasimha & Prabhu 1971). It implies in particular that a wake can be in equilibrium only if the longitudinal strain rate is appreciably less than the wake shear. We then describe a further series of experiments, undertaken to investigate the range of validity of the model. It is found that, with an appropriate convergence correction when necessary, the model provides excellent predictions of wake development under favourable, adverse and mixed pressure gradients. Furthermore, the behaviour of constant-pressure distorted wakes, as reported by Keffer (1965, 1967), is also explained very well by the model when account is taken of the effective flow convergence produced by the distortion. In all these calculations, only a simple version of the model is used, involving two non-dimensional constants both of which have been estimated from a single relaxation experiment.

Analysis of shear bands in slow granular flows using a frictional Cosserat model

The slow flow of granular materials is often marked by the existence of narrow shear layers, adjacent to large regions that suffer little or no deformation. This behaviour, in the regime where shear stress is generated primarily by the frictional interactions between grains, has so far eluded theoretical description. In this paper, we present a rigid-plastic frictional Cosserat model that captures thin shear layers by incorporating a microscopic length scale. We treat the granular medium as a Cosserat continuum, which allows the existence of localised couple stresses and, therefore, the possibility of an asymmetric stress tensor. In addition, the local rotation is an independent field variable and is not necessarily equal to the vorticity. The angular momentum balance, which is implicitly satisfied for a classical continuum, must now be solved in conjunction with the linear momentum balances. We extend the critical state model, used in soil plasticity, for a Cosserat continuum and obtain predictions for flow in plane and cylindrical Couette devices. The velocity profile predicted by our model is in qualitative agreement with available experimental data. In addition, our model can predict scaling laws for the shear layer thickness as a function of the Couette gap, which must be verified in future experiments. Most significantly, our model can determine the velocity field in viscometric flows, which classical plasticity-based model cannot.

Plume structure in high-Rayleigh-number convection

Near-wall structures in turbulent natural convection at Rayleigh numbers of $10^{10}$ to $10^{11}$ at A Schmidt number of 602 are visualized by a new method of driving the convection across a fine membrane using concentration differences of sodium chloride. The visualizations show the near-wall flow to consist of sheet plumes. A wide variety of large-scale flow cells, scaling with the cross-section dimension, are observed. Multiple large-scale flow cells are seen at aspect ratio (AR)= 0.65, while only a single circulation cell is detected at AR= 0.435. The cells (or the mean wind) are driven by plumes coming together to form columns of rising lighter fluid. The wind in turn aligns the sheet plumes along the direction of shear. the mean wind direction is seen to change with time. The near-wall dynamics show plumes initiated at points, which elongate to form sheets and then merge. Increase in rayleigh number results in a larger number of closely and regularly spaced plumes. The plume spacings show a common log–normal probability distribution function, independent of the rayleigh number and the aspect ratio. We propose that the near-wall structure is made of laminar natural-convection boundary layers, which become unstable to give rise to sheet plumes, and show that the predictions of a model constructed on this hypothesis match the experiments. Based on these findings, we conclude that in the presence of a mean wind, the local near-wall boundary layers associated with each sheet plume in high-rayleigh-number turbulent natural convection are likely to be laminar mixed convection type.