Numerical studies on quasilinear and linear elliptic equations

We explore here the acceleration of convergence of iterative methods for the solution of a class of quasilinear and linear algebraic equations. The specific systems are the finite difference form of the Navier-Stokes equations and the energy equation for recirculating flows. The acceleration procedures considered are: the successive over relaxation scheme; several implicit methods; and a second-order procedure. A new implicit method—the alternating direction line iterative method—is proposed in this paper. The method combines the advantages of the line successive over relaxation and alternating direction implicit methods. The various methods are tested for their computational economy and accuracy on a typical recirculating flow situation. The numerical experiments show that the alternating direction line iterative method is the most economical method of solving the Navier-Stokes equations for all Reynolds numbers in the laminar regime. The usual ADI method is shown to be not so attractive for large Reynolds numbers because of the loss of diagonal dominance. This loss can however be restored by a suitable choice of the relaxation parameter, but at the cost of accuracy. The accuracy of the new procedure is comparable to that of the well-tested successive overrelaxation method and to the available results in the literature. The second-order procedure turns out to be the most efficient method for the solution of the linear energy equation.

Structural studies of analgesics and their interactions. Part 4. Crystal structures of phenylbutazone and a 2 : 1 complex between phenylbutazone and piperazine

The X-ray crystal structures of 4-butyl-1,2-diphenylpyrazolidine-3,5-dione (phenylbutazone)(I). and its 2 : 1 complex (II) with piperazine have been determined by direct methods and the structures refined to R 0.096 (2 300 observed reflections measured by diffractometer) and 0.074 (2 494 observed reflections visuallyestimated). Crystals are monoclinic, space group P21/c; for (I)a= 21.695(4), b= 5.823(2), c= 27.881(4)Å, = 108.06 (10)°, Z= 8, and for (II)a= 8.048(4), b= 15.081(4), c= 15.583(7)Å, = 95.9(3)°, Z= 2. The two crystallographically independant molecules in the structure of (I) are similar except for the conformation of the butyl group, which is disordered in one of the molecules. In the pyrazolidinedione group, the two C–C bonds are single and the two C–O bonds double. The two nitrogen atoms in the five-membered ring are pyramidal with the attached phenyl groups lying on the opposite sides of the mean plane of the ring. The phenylbutazone molecule in (II) exists as a negative ion owing to deprotonation of C-4. C-4 is therefore trigonal and the orientation of the Bu group with respect to the pyrazolidinedione group is considerably different from that in (I); there is also considerable electron delocalization along the C–O and C–C bonds. These changes in geometry and electronic structure may relate to biological activity. The doubly charged cationic piperazine molecule exists in the chair form with the nitrogen atoms at the apices. The crystal structure of (II) is stabilized by ionic interactions and N–H O hydrogen bonds.

NUPARM and NUCGEN: software for analysis and generation of sequence dependent nucleic acid structures

Software packages NUPARM and NUCGEN, are described, which can be used to understand sequence directed structural variations in nucleic acids, by analysis and generation of non-uniform structures. A set of local inter basepair parameters (viz. tilt, roll, twist, shift, slide and rise) have been defined, which use geometry and coordinates of two successive basepairs only and can be used to generate polymeric structures with varying geometries for each of the 16 possible dinucleotide steps. Intra basepair parameters, propeller, buckle, opening and the C6...C8 distance can also be varied, if required, while the sugar phosphate backbone atoms are fixed in some standard conformation ill each of the nucleotides. NUPARM can be used to analyse both DNA and RNA structures, with single as well as double stranded helices. The NUCGEN software generates double helical models with the backbone fixed in B-form DNA, but with appropriate modifications in the input data, it can also generate A-form DNA ar rd RNA duplex structures.

Computational systems approach for drug target discovery

Importance of the field: The shift in focus from ligand based design approaches to target based discovery over the last two to three decades has been a major milestone in drug discovery research. Currently, it is witnessing another major paradigm shift by leaning towards the holistic systems based approaches rather the reductionist single molecule based methods. The effect of this new trend is likely to be felt strongly in terms of new strategies for therapeutic intervention, new targets individually and in combinations, and design of specific and safer drugs. Computational modeling and simulation form important constituents of new-age biology because they are essential to comprehend the large-scale data generated by high-throughput experiments and to generate hypotheses, which are typically iterated with experimental validation. Areas covered in this review: This review focuses on the repertoire of systems-level computational approaches currently available for target identification. The review starts with a discussion on levels of abstraction of biological systems and describes different modeling methodologies that are available for this purpose. The review then focuses on how such modeling and simulations can be applied for drug target discovery. Finally, it discusses methods for studying other important issues such as understanding targetability, identifying target combinations and predicting drug resistance, and considering them during the target identification stage itself. What the reader will gain: The reader will get an account of the various approaches for target discovery and the need for systems approaches, followed by an overview of the different modeling and simulation approaches that have been developed. An idea of the promise and limitations of the various approaches and perspectives for future development will also be obtained. Take home message: Systems thinking has now come of age enabling a `bird's eye view' of the biological systems under study, at the same time allowing us to `zoom in', where necessary, for a detailed description of individual components. A number of different methods available for computational modeling and simulation of biological systems can be used effectively for drug target discovery.

Crack Tip Fields in a Single Edge Notched Aluminum Single Crystal Specimen

We report a combined experimental and computational study of a low constraint aluminum single crystal fracture geometry and investigate the near-tip stress and strain fields. To this end, a single edge notched tensile (SENT) specimen is considered. A notch, with a radius of 50 µm, is taken to lie in the (010) plane and its front is aligned along the [101] direction. Experiments are conducted by subjecting the specimen to tensile loading using a special fixture inside a scanning electron microscope chamber. Both SEM micrographs and electron back-scattered diffraction (EBSD) maps are obtained from the near-tip region. The experiments are complemented by performing 3D and 2D plane strain finite element simulations within a continuum crystal plasticity framework assuming an isotropic hardening response characterized by the Pierce–Asaro–Needleman model. The simulations show a distinct slip band forming at about 55 deg with respect to the notch line corresponding to slip on (11-bar 1)[011] system, which corroborates well with experimental data. Furthermore, two kink bands occur at about 45 deg and 90 deg with respect to the notch line within which large rotations in the crystal orientation take place. These predictions are in good agreement with the EBSD observations. Finally, the near-tip angular variations of the 3D stress and plastic strain fields in the low constraint SENT fracture geometry are examined in detail.

On the Cubicity of AT-Free Graphs and Circular-Arc Graphs

A unit cube in k dimensions (k-cube) is defined as the Cartesian product R-1 x R-2 x ... x R-k where R-i (for 1 <= i <= k) is a closed interval of the form [a(i), a(i) + 1] on the real line. A graph G on n nodes is said to be representable as the intersection of k-cubes (cube representation in k dimensions) if each vertex of C can be mapped to a k-cube such that two vertices are adjacent in G if and only if their corresponding k-cubes have a non-empty intersection. The cubicity of G denoted as cub(G) is the minimum k for which G can be represented as the intersection of k-cubes. An interesting aspect about cubicity is that many problems known to be NP-complete for general graphs have polynomial time deterministic algorithms or have good approximation ratios in graphs of low cubicity. In most of these algorithms, computing a low dimensional cube representation of the given graph is usually the first step. We give an O(bw . n) algorithm to compute the cube representation of a general graph G in bw + 1 dimensions given a bandwidth ordering of the vertices of G, where bw is the bandwidth of G. As a consequence, we get O(Delta) upper bounds on the cubicity of many well-known graph classes such as AT-free graphs, circular-arc graphs and cocomparability graphs which have O(Delta) bandwidth. Thus we have: 1. cub(G) <= 3 Delta - 1, if G is an AT-free graph. 2. cub(G) <= 2 Delta + 1, if G is a circular-arc graph. 3. cub(G) <= 2 Delta, if G is a cocomparability graph. Also for these graph classes, there axe constant factor approximation algorithms for bandwidth computation that generate orderings of vertices with O(Delta) width. We can thus generate the cube representation of such graphs in O(Delta) dimensions in polynomial time.

DNA cleavage in red light promoted by copper(II) complexes of -amino acids and photoactive phenanthroline bases

Ternary copper(II) complexes [Cu(L-trp)(B)(H2O)](NO3) ( 1–3) and [Cu(L-phe)(B)(H2O)](NO3) ( 4–6) of L-tryptophan (L-trp) and L-phenylalanine (L-phe) having phenanthroline bases (B), viz. 1,10-phenanthroline (phen, 1 and 4), dipyrido[3,2-d:2,3-f]quinoxaline (dpq, 2 and 5) and dipyrido[3,2-a:2,3-c]phenazine (dppz, 3 and 6), were prepared and characterized by physico-chemical techniques. Complexes 3 and 6 were structurally characterized by X-ray crystallography and show the presence of a square pyramidal (4 + 1) CuN3O2 coordination geometry in which the N,O-donor amino acid (L-trp or L-phe) and N,N-donor phenanthroline base bind at the equatorial plane with an aqua ligand coordinated at the elongated axial site. Complex 3 shows significant distortion from the square pyramidal geometry and a strong intramolecular – stacking interaction between the pendant indole ring of L-trp and the planar dppz aromatic moiety. All the complexes display good binding propensity to the calf thymus DNA giving an order: 3, 6 (dppz) > 2, 5 (dpq) > 1, 4 (phen). The binding constant (Kb) values are in the range of 2.1 × 104–1.1 × 106 mol-1 with the binding site size (s) values of 0.17–0.63. The phen and dpq complexes are minor groove binders while the dppz analogues bind at the DNA major groove. Theoretical DNA docking studies on 2 and 3 show the close proximity of two photosensitizers, viz. the indole moiety of L-trp and the quinoxaline/phenazine of the dpq/dppz bases, to the complementary DNA strands. Complexes 2 and 3 show oxidative DNA double strand breaks (dsb) of supercoiled (SC) DNA forming a significant quantity of linear DNA along with the nicked circular (NC) form on photoexposure to UV-A light of 365 nm and red light of 647.1 nm (Ar–Kr laser). Complexes 1, 5 and 6 show only single strand breaks (ssb) forming NC DNA. The red light induced DNA cleavage involves metal-assisted photosensitization of L-trp and dpq/dppz base resulting in the formation of a reactive singlet oxygen (1O2) species.

Electromechanical interactions in a carbon nanotube based thin film field emitting diode

Carbon nanotubes (CNTs) have emerged as promising candidates for biomedical x-ray devices and other applications of field emission. CNTs grown/deposited in a thin film are used as cathodes for field emission. In spite of the good performance of such cathodes, the procedure to estimate the device current is not straightforward and the required insight towards design optimization is not well developed. In this paper, we report an analysis aided by a computational model and experiments by which the process of evolution and self-assembly (reorientation) of CNTs is characterized and the device current is estimated. The modeling approach involves two steps: (i) a phenomenological description of the degradation and fragmentation of CNTs and (ii) a mechanics based modeling of electromechanical interaction among CNTs during field emission. A computational scheme is developed by which the states of CNTs are updated in a time incremental manner. Finally, the device current is obtained by using the Fowler–Nordheim equation for field emission and by integrating the current density over computational cells. A detailed analysis of the results reveals the deflected shapes of the CNTs in an ensemble and the extent to which the initial state of geometry and orientation angles affect the device current. Experimental results confirm these effects.

Localization of handwritten text in documents using moment invariants and Delaunay triangulation

This paper describes an approach based on Zernike moments and Delaunay triangulation for localization of hand-written text in machine printed text documents. The Zernike moments of the image are first evaluated and we classify the text as hand-written using the nearest neighbor classifier. These features are independent of size, slant, orientation, translation and other variations in handwritten text. We then use Delaunay triangulation to reclassify the misclassified text regions. When imposing Delaunay triangulation on the centroid points of the connected components, we extract features based on the triangles and reclassify the text. We remove the noise components in the document as part of the preprocessing step so this method works well on noisy documents. The success rate of the method is found to be 86%. Also for specific hand-written elements such as signatures or similar text the accuracy is found to be even higher at 93%.

Line removal and restoration of handwritten strokes

In document images, we often find printed lines over-lapping with hand written elements especially in case of signatures. Typical examples of such images are bank cheques and payment slips. Although the detection and removal of the horizontal lines has been addressed, the restoration of the handwritten area after removal of lines, persists to be a problem of interest. lit this paper, we propose a method for line removal and restoration of the erased areas of the handwritten elements. Subjective evaluation of the results have been conducted to analyze the effectiveness of the proposed method. The results are promising with an accuracy of 86.33%. The entire Process takes less than half a second for completion on a 2.4 GHz 512 MB RAM Pentium IV PC for a document image.

Visual Analysis of Biomolecular Surfaces

Surface models of biomolecules have become crucially important for the study and understanding of interaction between biomolecules and their environment. We argue for the need for a detailed understanding of biomolecular surfaces by describing several applications in computational and structural biology. We review methods used to model, represent, characterize, and visualize biomolecular surfaces focusing on the role that geometry and topology play in identifying features on the surface. These methods enable the development of efficient computational and visualization tools for studying the function of biomolecules.

Control aspects in proportionality Markov manpower systems

In this paper, the control aspects of a hierarchical organization under the influence of "proportionality" policies are analyzed. Proportionality policies are those that restrict the recruitment to every level of the hierarchy (except the bottom most level or base level) to be in strict proportion to the promotions into that level. Both long term and short term control analysis have been discussed. In long term control the specific roles of the parameters of the system with regard to control of the shape and size of the system have been analyzed and yield suitable control strategies. In short term control, the attainability of a target or goal structure within a specific time from a given initial structure has been analyzed and yields the required recruitment strategies. The theoretical analyses have been illustrated with computational examples and also with real world data.

Selective dissemination of XML documents using GAs and SVM

XML has emerged as a medium for interoperability over the Internet. As the number of documents published in the form of XML is increasing there is a need for selective dissemination of XML documents based on user interests. In the proposed technique, a combination of Self Adaptive Migration Model Genetic Algorithm (SAMCA)[5] and multi class Support Vector Machine (SVM) are used to learn a user model. Based on the feedback from the users the system automatically adapts to the user's preference and interests. The user model and a similarity metric are used for selective dissemination of a continuous stream of XML documents. Experimental evaluations performed over a wide range of XML documents indicate that the proposed approach significantly improves the performance of the selective dissemination task, with respect to accuracy and efficiency.

Control aspects in proportionality Markov manpower systems

In this paper, the control aspects of a hierarchical organization under the influence of "proportionality" policies are analyzed. Proportionality policies are those that restrict the recruitment to every level of the hierarchy (except the bottom most level or base level) to be in strict proportion to the promotions into that level. Both long term and short term control analysis have been discussed. In long term control the specific roles of the parameters of the system with regard to control of the shape and size of the system have been analyzed and yield suitable control strategies. In short term control, the attainability of a target or goal structure within a specific time from a given initial structure has been analyzed and yields the required recruitment strategies. The theoretical analyses have been illustrated with computational examples and also with real world data. The control of such proportionality systems is then compared with that of the general systems (which do not follow such policies) with some significant conclusions. The control relations of such proportionality systems are found to be simpler and more practically feasible than those of general Markov systems, which do not have such restrictions. Such proportionality systems thus not only retain and match the flexibility of general Markov systems but also have the added advantage of simpler and more practically feasible controls. The proportionality policies hence act as an alternative and more practicably feasible means of control. (C) 2004 Elsevier Inc. All rights reserved.

Constructing Reeb Graphs using Cylinder Maps

The Reeb graph of a scalar function represents the evolution of the topology of its level sets. In this video, we describe a near-optimal output-sensitive algorithm for computing the Reeb graph of scalar functions defined over manifolds. Key to the simplicity and efficiency of the algorithm is an alternate definition of the Reeb graph that considers equivalence classes of level sets instead of individual level sets. The algorithm works in two steps. The first step locates all critical points of the function in the domain. Arcs in the Reeb graph are computed in the second step using a simple search procedure that works on a small subset of the domain that corresponds to a pair of critical points. The algorithm is also able to handle non-manifold domains.