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.

Computer modelling studies of ribonuclease A-pyrimidine nucleotide complexes

Different modes of binding of pyrimidine monophosphates 2'-UMP, 3'-UMP, 2'-CMP and 3'-CMP to ribonuclease (RNase) A are studied by energy minimization in torsion angle and subsequently in Cartesian coordinate space. The results are analysed in the light of primary binding sites. The hydrogen bonding pattern brings out roles for amino acids such as Asn44 and Ser123 apart from the well known active site residues viz., His12,Lys41,Thr45 and His119. Amino acid segments 43-45 and 119-121 seem to be guiding the ligand binding by forming a pocket. Many of the active site charged residues display considerable movement upon nucleotide binding.

Modelling and analysis of the variance in parallelism in parallel computations

In this paper, we introduce an analytical technique based on queueing networks and Petri nets for making a performance analysis of dataflow computations when executed on the Manchester machine. This technique is also applicable for the analysis of parallel computations on multiprocessors. We characterize the parallelism in dataflow computations through a four-parameter characterization, namely, the minimum parallelism, the maximum parallelism, the average parallelism and the variance in parallelism. We observe through detailed investigation of our analytical models that the average parallelism is a good characterization of the dataflow computations only as long as the variance in parallelism is small. However, significant difference in performance measures will result when the variance in parallelism is comparable to or higher than the average parallelism.

Integrated modelling approaches to analysis of climate change impacts on forests and forest management

This paper reviews integrated economic and ecological models that address impacts and adaptation to climate change in the forest sector. Early economic model studies considered forests as one out of many possible impacts of climate change, while ecological model studies tended to limit the economic impacts to fixed price-assumptions. More recent studies include broader representations of both systems, but there are still few studies which can be regarded fully integrated. Full integration of ecological and economic models is needed to address forest management under climate change appropriately. The conclusion so far is that there are vast uncertainties about how climate change affects forests. This is partly due to the limited knowledge about the global implications of the social and economical adaptation to the effects of climate change on forests.

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.

Characterization of substrate UpA binding to RNase A--computer modelling and energetics approach.

In the past two decades RNase A has been the focus of diverse investigations in order to understand the nature of substrate binding and to know the mechanism of enzyme action. Although this system is reasonably well characterized from the view point of some of the binding sites, the details of interactions in the second base binding (B2) site is insufficient. Further, the nature of ligand-protein interaction is elucidated generally by studies on RNase A-substrate analog complexes (mainly with the help of X-ray crystallography). Hence, the details of interactions at atomic level arising due to substrates are inferred indirectly. In the present paper, the dinucleotide substrate UpA is fitted into the active site of RNase A Several possible substrate conformations are investigated and the binding modes have been selected based on Contact Criteria. Thus identified RNase A-UpA complexes are energy minimized in coordinate space and are analysed in terms of conformations, energetics and interactions. The best possible ligand conformations for binding to RNase A are identified by experimentally known interactions and by the energetics. Upon binding of UpA to RNase A the changes associated,with protein back bone, Side chains in general and at the binding sites in particular are described. Further, the detailed interactions between UpA and RNase A are characterized in terms of hydrogen bonds and energetics. An extensive study has helped in interpreting the diverse results obtained from a number of experiments and also in evaluating the extent of changes the protein and the substrate undergo in order to maximize their interactions.

Modelling of a batch sonochemical reactor

Ultrasonication of aqueous KI solution is known to yield I2 due to reaction of iodide ions with hydroxyl radicals, which in turn are generated due to cavitation. Based on this conceptual framework, a model has been developed to predict the rate of iodine formation for KI solutions of various concentrations under different gas atmospheres. The model follows the growth and collapse of a gas—vapour cavity using the Rayleigh—Plesset bubble dynamics equation. The bubble is assumed to behave isothermally during its growth phase and a part of the collapse phase. Thereafter it is assumed to collapse adiabatically, yielding high temperatures and pressures. Thermodynamic equilibrium is assumed in the bubble at the end of collapse phase. The contents of the bubble are assumed to mix with the liquid, and the reactor contents are assumed to be well stirred. The model has been verified by conducting experiments with KI solutions of different concentrations and using different gas atmospheres. The model not only explains these results but also the existence of a maximum when Ar---O2 mixtures of different compositions are employed.

Modelling of a batch sonochemical reactor

Ultrasonication of aqueous KI solution is known to yield I2 due to reaction of iodide ions with hydroxyl radicals, which in turn are generated due to cavitation. Based on this conceptual framework, a model has been developed to predict the rate of iodine formation for KI solutions of various concentrations under different gas atmospheres. The model follows the growth and collapse of a gas-vapour cavity using the Rayleigh-Plesset bubble dynamics equation. The bubble is assumed to behave isothermally during its growth phase and a part of the collapse phase. Thereafter it is assumed to collapse adiabatically, yielding high temperatures and pressures. Thermodynamic equilibrium is assumed in the bubble at the end of collapse phase. The contents of the bubble are assumed to mix with the liquid, and the reactor contents are assumed to be well stirred. The model has been verified by conducting experiments with KI solutions of different concentrations and using different gas atmospheres. The model not only explains these results but also the existence of a maximum when Ar-O2 mixtures of different compositions are employed.

A note on the choice of interfacial friction angle

Estimates of interfacial friction angle (delta) are necessary for the design of retaining structures and deep foundations, Recommendations in the literature regarding delta values are often contradictory and are therefore not easy to apply in geotechnical design, A critical examination of past studies in terms of data generation techniques used and conclusions drawn indicates that two distinctly different test procedures/techniques have been evolved. The interfacial situation in practice can also be categorized into two broad types, These two types of interface problems in geotechnical engineering are (a) the structure is placed on the free surface of prepared fill (type A situation) and (b) the fill is placed against the material surface which functions as a confined boundary (type B situation), The friction angle delta depends on the surface roughness of the construction material, But in the type A situation, it is independent of density and its limiting maximum value (delta(lim)) is the critical state friction angle phi(cv). In the type B situation, it is dependent on density of the fill and its limiting maximum value is the peak angle of internal friction phi(p) of the fill.

Machine learning modelling for predicting soil liquefaction susceptibility

This study describes two machine learning techniques applied to predict liquefaction susceptibility of soil based on the standard penetration test (SPT) data from the 1999 Chi-Chi, Taiwan earthquake. The first machine learning technique which uses Artificial Neural Network (ANN) based on multi-layer perceptions (MLP) that are trained with Levenberg-Marquardt backpropagation algorithm. The second machine learning technique uses the Support Vector machine (SVM) that is firmly based on the theory of statistical learning theory, uses classification technique. ANN and SVM have been developed to predict liquefaction susceptibility using corrected SPT (N-1)(60)] and cyclic stress ratio (CSR). Further, an attempt has been made to simplify the models, requiring only the two parameters (N-1)(60) and peck ground acceleration (a(max)/g)], for the prediction of liquefaction susceptibility. The developed ANN and SVM models have also been applied to different case histories available globally. The paper also highlights the capability of the SVM over the ANN models.

Joining of dissimilar metals: issues and modelling techniques

Joining of dissimilar metals involves a number of scientific issues, the modelling of which offers unique challenges. This review discusses the complexities in different joining processes and dissimilar combinations, and the corresponding computational techniques that have the potential to address the same. Future directions in modelling at both macroscopic and microscopic scales are also suggested.

Stochastic modelling of flexible manufacturing systems

Mathematical modelling plays a vital role in the design, planning and operation of flexible manufacturing systems (FMSs). In this paper, attention is focused on stochastic modelling of FMSs using Markov chains, queueing networks, and stochastic Petri nets. We bring out the role of these modelling tools in FMS performance evaluation through several illustrative examples and provide a critical comparative evaluation. We also include a discussion on the modelling of deadlocks which constitute an important source of performance degradation in fully automated FMSs.

Choice of Classifiers in Hierarchical Recognition of Online Handwritten Kannada and Tamil Aksharas

In this paper, we propose a novel dexterous technique for fast and accurate recognition of online handwritten Kannada and Tamil characters. Based on the primary classifier output and prior knowledge, the best classifier is chosen from set of three classifiers for second stage classification. Prior knowledge is obtained through analysis of the confusion matrix of primary classifier which helped in identifying the multiple sets of confused characters. Further, studies were carried out to check the performance of secondary classifiers in disambiguating among the confusion sets. Using this technique we have achieved an average accuracy of 92.6% for Kannada characters on the MILE lab dataset and 90.2% for Tamil characters on the HP Labs dataset.

Thermodynamic modelling of absorption-resorption heating cycles with some new working pairs

This paper presents the results of a thermodynamic cycle analysis of single stage resorption heat pump (RHP) and resorption heat transformer (RHT) cycles with the new working pairs R22-NMP and R22-DMA. The coefficients of performance (COP) are correlated with the low grade source temperature, temperature at which useful heat is obtained and ambient temperature. The COPs are in the range 1.20–1.60 for the RHP mode and 0.25–0.45 for the RHT mode. Absorber temperatures (useful temperatures) as high as 50°C in the RHP mode and 87°C in the RHT mode have been obtained. It is observed that absorption-resorption systems are inflexible in their range of operating temperature and necessitate a higher pump work as compared with simple single-stage absorption heating systems. However, single stage RHTs show higher temperature boosts than simple absorption heat transformers.