What mainly controls recession flows in river basins?

The ubiquity of the power law relationship between dQ/dt and Q for recession periods (-dQ/dt kQ(alpha); Q being discharge at the basin outlet at time t) clearly hints at the existence of a dominant recession flow process that is common to all real basins. It is commonly assumed that a basin, during recession events, functions as a single phreatic aquifer resting on a impermeable horizontal bed or the Dupuit-Boussinesq (DB) aquifer, and with time different aquifer geometric conditions arise that give different values of alpha and k. The recently proposed alternative model, geomorphological recession flow model, however, suggests that recession flows are controlled primarily by the dynamics of the active drainage network (ADN). In this study we use data for several basins and compare the above two contrasting recession flow models in order to understand which of the above two factors dominates during recession periods in steep basins. Particularly, we do the comparison by selecting three key recession flow properties: (1) power law exponent alpha, (2) dynamic dQ/dt-Q relationship (characterized by k) and (3) recession timescale (time period for which a recession event lasts). Our observations suggest that neither drainage from phreatic aquifers nor evapotranspiration significantly controls recession flows. Results show that the value of a and recession timescale are not modeled well by DB aquifer model. However, the above mentioned three recession curve properties can be captured satisfactorily by considering the dynamics of the ADN as described by geomorphological recession flow model, possibly indicating that the ADN represents not just phreatic aquifers but the organization of various sub-surface storage systems within the basin. (C) 2014 Elsevier Ltd. All rights reserved.

`Universal' recession curves and their geomorphological interpretation

The study of recession flows offers fundamental insights into basin hydrological processes and, in particular, into the collective behavior of the governing dominant subsurface flows and properties. We use here an existing geomorphological interpretation of recession dynamics, which links the exponent in the classic recession curve -dQ/dt - kQ(alpha) to the geometric properties of the time-varying drainage network to study the general properties of recession curves across a wide variety of river basins. In particular, we show how the parameter k depends on the initial soil moisture state of the basin and can be made to explicitly depend on an index discharge, representative of initial sub-subsurface storage. Through this framework we obtain a non-dimensional, event-independent, recession curve. We subsequently quantify the variability of k across different basins on the basis of their geometry, and, by rescaling, collapse curves from different events and basins to obtain a generalized, or `universal', recession curve. Finally, we analyze the resulting normalized recession curves and explain their universal characteristics, lending further support to the notion that the statistical properties of observed recession curves bear the signature of the geomorphological structure of the networks producing them. (C) 2014 Elsevier Ltd. All rights reserved.

Analysis of composite single lap joints using numerical and experimental approach

Adhesives are widely used to execute the assembly of aerospace and automotive structures due to their ability to join dissimilar materials, reduced stress concentration, and improved fatigue resistance. The mechanical behavior of adhesive joints can be studied either using analytical models or by conducting mechanical tests. However, the complexity owing to multiple interfaces, layers with different properties, material and geometric nonlinearity and its three-dimensional nature combine to increase the difficulty in obtaining an overall system of governing equations to predict the joint behavior. On the other hand, experiments are often time consuming and expensive due to a number of parameters involved. Finite element analysis (FEA) is profoundly used in recent years to overcome these limitations. The work presented in this paper involves the finite element modeling and analysis of a composite single lap joint where the adhesive-adherend interface region was modeled using connector elements. The computed stresses were compared with the experimental stresses obtained using digital image correlation technique. The results showed an agreement. Further, the failure load predicted using FEA was found to be closer to the actual failure load obtained by mechanical tests.

Dynamics of a driven magneto-martensitic ribbon

We investigate the dynamics of a sinusoidally driven ferromagnetic martensitic ribbon by adopting a recently introduced model that involves strain and magnetization as order parameters. Retaining only the dominant mode of excitation we reduce the coupled set of partial differential equations for strain and magnetization to a set of coupled ordinary nonlinear equations for the strain and magnetization amplitudes. The equation for the strain amplitude takes the form of parametrically driven oscillator. Finite strain amplitude can only be induced beyond a critical value of the strength of the magnetic field. Chaotic response is seen for a range of values of all the physically interesting parameters. The nature of the bifurcations depends on the choice of temperature relative to the ordering of the Curie and the martensite transformation temperatures. We have studied the nature of response as a function of the strength and frequency of the magnetic field, and magneto-elastic coupling. In general, the bifurcation diagrams with respect to these parameters do not follow any standard route. The rich dynamics exhibited by the model is further illustrated by the presence of mixed mode oscillations seen for low frequencies. The geometric structure of the mixed mode oscillations in the phase space has an unusual deep crater structure with an outer and inner cone on which the orbits circulate. We suggest that these features should be seen in experiments on driven magneto-martensitic ribbons. (C) 2014 Elsevier B. V. All rights reserved.

Effect of quantum nuclear motion on hydrogen bonding

This work considers how the properties of hydrogen bonded complexes, X-H center dot center dot center dot Y, are modified by the quantum motion of the shared proton. Using a simple two-diabatic state model Hamiltonian, the analysis of the symmetric case, where the donor (X) and acceptor (Y) have the same proton affinity, is carried out. For quantitative comparisons, a parametrization specific to the O-H center dot center dot center dot O complexes is used. The vibrational energy levels of the one-dimensional ground state adiabatic potential of the model are used to make quantitative comparisons with a vast body of condensed phase data, spanning a donor-acceptor separation (R) range of about 2.4-3.0 angstrom, i.e., from strong to weak hydrogen bonds. The position of the proton (which determines the X-H bond length) and its longitudinal vibrational frequency, along with the isotope effects in both are described quantitatively. An analysis of the secondary geometric isotope effect, using a simple extension of the two-state model, yields an improved agreement of the predicted variation with R of frequency isotope effects. The role of bending modes is also considered: their quantum effects compete with those of the stretching mode for weak to moderate H-bond strengths. In spite of the economy in the parametrization of the model used, it offers key insights into the defining features of H-bonds, and semi-quantitatively captures several trends. (C) 2014 AIP Publishing LLC.

Transition in vortex breakdown modes in a coaxial isothermal unconfined swirling jet

This paper reports first observations of transition in recirculation pattern from an open-bubble type axisymmetric vortex breakdown to partially open bubble mode through an intermediate, critical regime of conical sheet formation in an unconfined, co-axial isothermal swirling flow. This time-mean transition is studied for two distinct flow modes which are characterized based on the modified Rossby number (Ro(m)), i.e., Ro(m) <= 1 and Ro(m) > 1. Flow modes with Ro(m) <= 1 are observed to first undergo cone-type breakdown and then to partially open bubble state as the geometric swirl number (S-G) is increased by similar to 20% and similar to 40%, respectively, from the baseline open-bubble state. However, the flow modes with Ro(m) > 1 fail to undergo such sequential transition. This distinct behavior is explained based on the physical significance associated with Ro(m) and the swirl momentum factor (xi). In essence, xi represents the ratio of angular momentum distributed across the flow structure to that distributed from central axis to the edge of the vortex core. It is observed that xi increases by similar to 100% in the critical swirl number band where conical breakdown occurs as compared to its magnitude in the S-G regime where open bubble state is seen. This results from the fact that flow modes with Ro(m) <= 1 are dominated by radial pressure gradient due to swirl/rotational effect when compared to radial pressure deficit arising from entrainment (due to the presence of co-stream). Consequently, the imparted swirl tends to penetrate easily towards the central axis causing it to spread laterally and finally undergo conical sheet breakdown. However, the flow modes with Ro(m) > 1 are dominated by pressure deficit due to entrainment effect. This blocks the radial inward penetration of imparted angular momentum thus preventing the lateral spread of these flow modes. As such these structures fail to undergo cone mode of vortex breakdown which is substantiated by a mere 30%-40% rise in xi in the critical swirl number range. (C) 2014 AIP Publishing LLC.

Time-Optimal Convergence to a Rectilinear Path in the Presence of Wind

This paper considers the problem of determining the time-optimal path of a fixed-wing Miniature Air Vehicle (MAV), in the presence of wind. The MAV, which is subject to a bounded turn rate, is required to eventually converge to a straight line starting from a known initial position and orientation. Earlier work in the literature uses Pontryagin's Minimum Principle (PMP) to solve this problem only for the no-wind case. In contrast, the present work uses a geometric approach to solve the problem completely in the presence of wind. In addition, it also shows how PMP can be used to partially solve the problem. Using a 6-DOF model of a MAV the generated optimal path is tracked by an autopilot consisting of proportional-integral-derivative (PID) controllers. The simulation results show the path generation and tracking for cases with steady and time-varying wind. Some issues on real-time path planning are also addressed.

Extended liaison as an interface between product and process model in assembly

This paper describes the use of liaison to better integrate product model and assembly process model so as to enable sharing of design and assembly process information in a common integrated form and reason about them. Liaison can be viewed as a set, usually a pair, of features in proximity with which process information can be associated. A liaison is defined as a set of geometric entities on the parts being assembled and relations between these geometric entities. Liaisons have been defined for riveting, welding, bolt fastening, screw fastening, adhesive bonding (gluing) and blind fastening processes. The liaison captures process specific information through attributes associated with it. The attributes are associated with process details at varying levels of abstraction. A data structure for liaison has been developed to cluster the attributes of the liaison based on the level of abstraction. As information about the liaisons is not explicitly available in either the part model or the assembly model, algorithms have been developed for extracting liaisons from the assembly model. The use of liaison is proposed to enable both the construction of process model as the product model is fleshed out, as well as maintaining integrity of both product and process models as the inevitable changes happen to both design and the manufacturing environment during the product lifecycle. Results from aerospace and automotive domains have been provided to illustrate and validate the use of liaisons. (C) 2014 Elsevier Ltd. All rights reserved.

Toward a Scalable Working Set Size Estimation Method and Its Application for Chip Multiprocessors

It is essential to accurately estimate the working set size (WSS) of an application for various optimizations such as to partition cache among virtual machines or reduce leakage power dissipated in an over-allocated cache by switching it OFF. However, the state-of-the-art heuristics such as average memory access latency (AMAL) or cache miss ratio (CMR) are poorly correlated to the WSS of an application due to 1) over-sized caches and 2) their dispersed nature. Past studies focus on estimating WSS of an application executing on a uniprocessor platform. Estimating the same for a chip multiprocessor (CMP) with a large dispersed cache is challenging due to the presence of concurrently executing threads/processes. Hence, we propose a scalable, highly accurate method to estimate WSS of an application. We call this method ``tagged WSS (TWSS)'' estimation method. We demonstrate the use of TWSS to switch-OFF the over-allocated cache ways in Static and Dynamic NonUniform Cache Architectures (SNUCA, DNUCA) on a tiled CMP. In our implementation of adaptable way SNUCA and DNUCA caches, decision of altering associativity is taken by each L2 controller. Hence, this approach scales better with the number of cores present on a CMP. It gives overall (geometric mean) 26% and 19% higher energy-delay product savings compared to AMAL and CMR heuristics on SNUCA, respectively.

Crack stability in edge-notched clamped beam specimens: modeling and experiments

Stability of a fracture toughness testing geometry is important to determine the crack trajectory and R-curve behavior of the specimen. Few configurations provide for inherent geometric stability, especially when the specimen being tested is brittle. We propose a new geometrical construction called the single edge notched clamped bend specimen (SENCB), a modified form of three point bending, yielding stable cracking under load control. It is shown to be particularly suitable for small-scale structures which cannot be made free-standing, (e.g., thin films, coatings). The SENCB is elastically clamped at the two ends to its parent material. A notch is inserted at the bottom center and loaded in bending, to fracture. Numerical simulations are carried out through extended finite element method to derive the geometrical factor f(a/W) and for different beam dimensions. Experimental corroborations of the FEM results are carried out on both micro-scale and macro-scale brittle specimens. A plot of vs a/W, is shown to rise initially and fall off, beyond a critical a/W ratio. The difference between conventional SENB and SENCB is highlighted in terms of and FEM simulated stress contours across the beam cross-section. The `s of bulk NiAl and Si determined experimentally are shown to match closely with literature values. Crack stability and R-curve effect is demonstrated in a PtNiAl bond coat sample and compared with predicted crack trajectories from the simulations. The stability of SENCB is shown for a critical range of a/W ratios, proving that it can be used to get controlled crack growth even in brittle samples under load control.

Halogen Bonds in Crystal Engineering: Like Hydrogen Bonds yet Different

CONSPECTUS: The halogen bond is an attractive interaction in which an electrophilic halogen atom approaches a negatively polarized species. Short halogen atom contacts in crystals have been known for around 50 years. Such contacts are found in two varieties: type I, which is symmetrical, and type II, which is bent. Both are influenced by geometric and chemical considerations. Our research group has been using halogen atom interactions as design elements in crystal engineering, for nearly 30 years. These interactions include halogen center dot center dot center dot halogen interactions (X center dot center dot center dot X) and halogen center dot center dot center dot heteroatom interactions (X center dot center dot center dot B). Many X center dot center dot center dot X and almost all X center dot center dot center dot B contacts can be classified as halogen bonds. In this Account, we illustrate examples of crystal engineering where one can build up from previous knowledge with a focus that is provided by the modern definition of the halogen bond. We also comment on the similarities and differences between halogen bonds and hydrogen bonds. These interactions are similar because the protagonist atoms halogen and hydrogen are both electrophilic in nature. The interactions are distinctive because the size of a halogen atom is of consequence when compared with the atomic sizes of, for example, C, N, and O, unlike that of a hydrogen atom. Conclusions may be drawn pertaining to the nature of X center dot center dot center dot X interactions from the Cambridge Structural Database (CSD). There is a clear geometric and chemical distinction between type I and type II, with only type II being halogen bonds. Cl/Br isostructurality is explained based on a geometric model. In parallel, experimental studies on 3,4-dichlorophenol and its congeners shed light on the nature of halogen center dot center dot center dot halogen interactions and reveal the chemical difference between Cl and Br. Variable temperature studies also show differences between type I and type II contacts. In terms of crystal design, halogen bonds offer a unique opportunity in the strength, atom size and interaction gradation; this may be used in the design of ternary cocrystals. Structural modularity in which an entire crystal structure is defined as a combination of modules is rationalized on the basis of the intermediate strength of a halogen bond. The specific directionality of the halogen bond makes it a good tool to achieve orthogonality in molecular crystals. Mechanical properties can be tuned systematically by varying these orthogonally oriented halogen center dot center dot center dot halogen interactions. In a further development, halogen bonds are shown to play a systematic role in organization of LSAMs (long range synthon aufbau module), which are bigger structural units containing multiple synthons. With a formal definition in place, this may be the right time to look at differences between halogen bonds and hydrogen bonds and exploit them in more subtle ways in crystal engineering.

MolBridge: a program for identifying nonbonded interactions in small molecules and biomolecular structures

Identification and analysis of nonbonded interactions within a molecule and with the surrounding molecules are an essential part of structural studies, given the importance of these interactions in defining the structure and function of any supramolecular entity. MolBridge is an easy to use algorithm based purely on geometric criteria that can identify all possible nonbonded interactions, such as hydrogen bond, halogen bond, cation-pi, pi-pi and van der Waals, in small molecules as well as biomolecules. The user can either upload three-dimensional coordinate files or enter the molecular ID corresponding to the relevant database. The program is available in a standalone form and as an interactive web server with Jmol and JME incorporated into it. The program is freely downloadable and the web server version is also available at http://nucleix.mbu.iisc.ernet.in/molbridge/index.php.

Turbulent electrical activity at sharp-edged inexcitable obstacles in a model for human cardiac tissue

Wave propagation around various geometric expansions, structures, and obstacles in cardiac tissue may result in the formation of unidirectional block of wave propagation and the onset of reentrant arrhythmias in the heart. Therefore, we investigated the conditions under which reentrant spiral waves can be generated by high-frequency stimulation at sharp-edged obstacles in the ten Tusscher-Noble-Noble-Panfilov (TNNP) ionic model for human cardiac tissue. We show that, in a large range of parameters that account for the conductance of major inward and outward ionic currents of the model fast inward Na+ current (INa), L-type slow inward Ca2+ current (I-CaL), slow delayed-rectifier current (I-Ks), rapid delayed-rectifier current (I-Kr), inward rectifier K+ current (I-K1)], the critical period necessary for spiral formation is close to the period of a spiral wave rotating in the same tissue. We also show that there is a minimal size of the obstacle for which formation of spirals is possible; this size is similar to 2.5 cm and decreases with a decrease in the excitability of cardiac tissue. We show that other factors, such as the obstacle thickness and direction of wave propagation in relation to the obstacle, are of secondary importance and affect the conditions for spiral wave initiation only slightly. We also perform studies for obstacle shapes derived from experimental measurements of infarction scars and show that the formation of spiral waves there is facilitated by tissue remodeling around it. Overall, we demonstrate that the formation of reentrant sources around inexcitable obstacles is a potential mechanism for the onset of cardiac arrhythmias in the presence of a fast heart rate.

Fixed-orientation equilateral triangle matching of point sets

Given a point set P and a class C of geometric objects, G(C)(P) is a geometric graph with vertex set P such that any two vertices p and q are adjacent if and only if there is some C is an element of C containing both p and q but no other points from P. We study G(del)(P) graphs where del is the class of downward equilateral triangles (i.e., equilateral triangles with one of their sides parallel to the x-axis and the corner opposite to this side below that side). For point sets in general position, these graphs have been shown to be equivalent to half-Theta(6) graphs and TD-Delaunay graphs. The main result in our paper is that for point sets P in general position, G(del)(P) always contains a matching of size at least vertical bar P vertical bar-1/3] and this bound is tight. We also give some structural properties of G(star)(P) graphs, where is the class which contains both upward and downward equilateral triangles. We show that for point sets in general position, the block cut point graph of G(star)(P) is simply a path. Through the equivalence of G(star)(P) graphs with Theta(6) graphs, we also derive that any Theta(6) graph can have at most 5n-11 edges, for point sets in general position. (C) 2013 Elsevier B.V. All rights reserved.

Transverse orthotropic elastic moduli of bundled coated thin elliptical tubes

Light weight structures with tailored mechanical properties have evolved beyond regular hexagonal/circular honeycomb topology. For applications which demand anisotropic mechanical properties, elliptical-celled structures offer interesting features. This paper characterizes the anisotropic in-plane elastic response of coated thin elliptical tubes in different array patterns viz, close-packed, diagonal and rectangular patterns under compression. This paper also extends earlier works on elliptical close-packed structure to a more general case of coated tubes. Theoretical framework using thin ring theory provides formulae in terms of geometric and material parameters. These are compared with a series of FE simulations using contact elements. The FE results are presented as graphs to aid in design. (C) 2014 Elsevier Ltd. All rights reserved.