Optimised form of acceleration correction algorithm within SPH-based simulations of impact mechanics

In the context of SPH-based simulations of impact dynamics, an optimised and automated form of the acceleration correction algorithm (Shaw and Reid, 2009a) is developed so as to remove spurious high frequency oscillations in computed responses whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. A rational framework for an insightful characterisation of the erstwhile acceleration correction method is first set up. This is followed by the proposal of an optimised version of the method, wherein the strength of the correction term in the momentum balance and energy equations is optimised. For the first time, this leads to an automated procedure to arrive at the artificial viscosity term. In particular, this is achieved by taking a spatially varying response-dependent support size for the kernel function through which the correction term is computed. The optimum value of the support size is deduced by minimising the (spatially localised) total variation of the high oscillation in the acceleration term with respect to its (local) mean. The derivation of the method, its advantages over the heuristic method and issues related to its numerical implementation are discussed in detail. (C) 2011 Elsevier Ltd. All rights reserved.

Decoupling of diffusion from viscosity: Difference scenario for translational and rotational motions

Recent experiments have indicated a dramatically different viscosity dependence of the translational and the rotational diffusion coefficients in a supercooled liquid as the glass transition temperature is approached from above. While the translational motion seems to be decoupled from the rising viscosity (eta), the rotational motion seems to remain firmly coupled to eta. In order to understand the microscopic origin of this behavior, we have carried nut detailed theoretical calculations of both the quantities by using a self-consistent mode-coupling theory (MCT). it is found that when the size of the solute is same as that of the solvent molecules, the conventional MCT fails to predict the observed decoupling. The solvent inhomogeneity is found to play a decisive role in determining the decoupling. The difference in the viscosity dependence between rotation and translational diffusion coefficient is discussed.

Implications of a viscosity bound on black hole accretion

Motivated by the viscosity bound in gauge/gravity duality, we consider the ratio of shear viscosity (eta) to entropy density (s) in black hole accretion flows. We use both an ideal gas equation of state and the QCD equation of state obtained from lattice for the fluid accreting onto a Kerr black hole. The QCD equation of state is considered since the temperature of accreting matter is expected to approach 10(12) K in certain hot flows. We find that in both the cases eta/s is small only for primordial black holes and several orders of magnitude larger than any known fluid for stellar and supermassive black holes. We show that a lower bound on the mass of primordial black holes leads to a lower bound on eta/s and vice versa. Finally we speculate that the Shakura-Sunyaev viscosity parameter should decrease with increasing density and/or temperatures. (C) 2012 Elsevier B.V. All rights reserved.

Stabilized SPH-based simulations of impact dynamics using acceleration-corrected artificial viscosity

Artificial viscosity in SPH-based computations of impact dynamics is a numerical artifice that helps stabilize spurious oscillations near the shock fronts and requires certain user-defined parameters. Improper choice of these parameters may lead to spurious entropy generation within the discretized system and make it over-dissipative. This is of particular concern in impact mechanics problems wherein the transient structural response may depend sensitively on the transfer of momentum and kinetic energy due to impact. In order to address this difficulty, an acceleration correction algorithm was proposed in Shaw and Reid (''Heuristic acceleration correction algorithm for use in SPH computations in impact mechanics'', Comput. Methods Appl. Mech. Engrg., 198, 3962-3974) and further rationalized in Shaw et al. (An Optimally Corrected Form of Acceleration Correction Algorithm within SPH-based Simulations of Solid Mechanics, submitted to Comput. Methods Appl. Mech. Engrg). It was shown that the acceleration correction algorithm removes spurious high frequency oscillations in the computed response whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. In this paper, we aim at gathering further insights into the acceleration correction algorithm by further exploring its application to problems related to impact dynamics. The numerical evidence in this work thus establishes that, together with the acceleration correction algorithm, SPH can be used as an accurate and efficient tool in dynamic, inelastic structural mechanics. (C) 2011 Elsevier Ltd. All rights reserved.

An experimental study on cracking evolution in concrete and cement mortar by the b-value analysis of acoustic emission technique

Notched three point bend (TPB) specimens made with plain concrete and cement mortar were tested under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/s and simultaneously acoustic emissions (AE) released were recorded during the experiments. Amplitude distribution analysis of AE released during concrete was carried out to study the development of fracture process in concrete and mortar specimens. The slope of the log-linear frequency-amplitude distribution of AE is known as the AE based b-value. The AE based b-value was computed in terms of physical process of time varying applied load using cumulative frequency distribution (Gutenberg-Richter relationship) and discrete frequency distribution (Aki's method) of AE released during concrete fracture. AE characteristics of plain concrete and cement mortar were studied and discussed and it was observed that the AE based b-value analysis serves as a tool to identify the damage in concrete structural members. (C) 2012 Elsevier Ltd. All rights reserved.

Ionic conductivity of bis(2-cyanoethyl) ether-lithium salt and poly(propylether imine)-lithium salt liquid electrolytes

We report here a multiple-nitrile based lithium-salt liquid electrolyte. The ionic conductivity of poly (propyl ether imine) (abbreviated as PETIM) lithium salt dendrimer liquid electrolyte was observed to be a function of dendrimer generation number, n=0 (monomer)-3. While the highest room temperature ionic conductivity value (similar to 10(-1) Sm-1) was recorded for the bis-2cyanoethyl ether monomer (i.e. zeroth generation; G(0)-CN), conductivity decreased progressively to lower values (similar to 10(-3) Sm-1) with increase in generation number (G(1)-CN -> G(3)-CN). The G(0)-CN and higher dendrimer generations showed high thermal stability (approximate to 150 to 200 degrees C), low moisture sensitivity and tunable viscosity (similar to 10(-2) (G(0)-CN) to 3 (G(3)-CN) Pa s). The linker ether group was found to be crucial for ion transport and also eliminated a large number of detrimental features, chiefly moisture sensitivity, chemical instability associated typically with prevalent molecular liquid solvents. Based on the combination of several beneficial physicochemical properties, we presently envisage that the PETIM dendrimers especially the G(0)-CN electrolytes hold promise as electrolytes in electrochemical devices such as lithium-ion batteries.

A QUADRATIC C degrees INTERIOR PENALTY METHOD FOR LINEAR FOURTH ORDER BOUNDARY VALUE PROBLEMS WITH BOUNDARY CONDITIONS OF THE CAHN-HILLIARD TYPE

We develop a quadratic C degrees interior penalty method for linear fourth order boundary value problems with essential and natural boundary conditions of the Cahn-Hilliard type. Both a priori and a posteriori error estimates are derived. The performance of the method is illustrated by numerical experiments.

Influence limitation in multi-campaign social networks: A Shapley value based approach

We investigate the problem of influence limitation in the presence of competing campaigns in a social network. Given a negative campaign which starts propagating from a specified source and a positive/counter campaign that is initiated, after a certain time delay, to limit the the influence or spread of misinformation by the negative campaign, we are interested in finding the top k influential nodes at which the positive campaign may be triggered. This problem has numerous applications in situations such as limiting the propagation of rumor, arresting the spread of virus through inoculation, initiating a counter-campaign against malicious propaganda, etc. The influence function for the generic influence limitation problem is non-submodular. Restricted versions of the influence limitation problem, reported in the literature, assume submodularity of the influence function and do not capture the problem in a realistic setting. In this paper, we propose a novel computational approach for the influence limitation problem based on Shapley value, a solution concept in cooperative game theory. Our approach works equally effectively for both submodular and non-submodular influence functions. Experiments on standard real world social network datasets reveal that the proposed approach outperforms existing heuristics in the literature. As a non-trivial extension, we also address the problem of influence limitation in the presence of multiple competing campaigns.

Seismic site classification and correlation between standard penetration test N value and shear wave velocity for Lucknow city in Indo-Gangetic basin

Subsurface lithology and seismic site classification of Lucknow urban center located in the central part of the Indo-Gangetic Basin (IGB) are presented based on detailed shallow subsurface investigations and borehole analysis. These are done by carrying out 47 seismic surface wave tests using multichannel analysis of surface waves (MASW) and 23 boreholes drilled up to 30 m with standard penetration test (SPT) N values. Subsurface lithology profiles drawn from the drilled boreholes show low- to medium-compressibility clay and silty to poorly graded sand available till depth of 30 m. In addition, deeper boreholes (depth >150 m) were collected from the Lucknow Jal Nigam (Water Corporation), Government of Uttar Pradesh to understand deeper subsoil stratification. Deeper boreholes in this paper refer to those with depth over 150 m. These reports show the presence of clay mix with sand and Kankar at some locations till a depth of 150 m, followed by layers of sand, clay, and Kankar up to 400 m. Based on the available details, shallow and deeper cross-sections through Lucknow are presented. Shear wave velocity (SWV) and N-SPT values were measured for the study area using MASW and SPT testing. Measured SWV and N-SPT values for the same locations were found to be comparable. These values were used to estimate 30 m average values of N-SPT (N-30) and SWV (V-s(30)) for seismic site classification of the study area as per the National Earthquake Hazards Reduction Program (NEHRP) soil classification system. Based on the NEHRP classification, the entire study area is classified into site class C and D based on V-s(30) and site class D and E based on N-30. The issue of larger amplification during future seismic events is highlighted for a major part of the study area which comes under site class D and E. Also, the mismatch of site classes based on N-30 and V-s(30) raises the question of the suitability of the NEHRP classification system for the study region. Further, 17 sets of SPT and SWV data are used to develop a correlation between N-SPT and SWV. This represents a first attempt of seismic site classification and correlation between N-SPT and SWV in the Indo-Gangetic Basin.

Synthesis, antimicrobial, DNA-binding and photonuclease studies of Cobalt(III) and Nickel(II) Schiff base complexes

New metal complexes of the type M(nih)(L)](PF6)(n)center dot xAH(2)O and M(nih)(2)](PF6)center dot xH(2)O (where M = Co(III) or Ni(II), L = 1,10-phenanthroline (phen)/or 2,2' bipyridine (bpy), nih = 2-hydroxy-1-naphthaldehyde isonicotinoyl hydrazone, n = 2 or 1 and x = 3 or 2) have been synthesized and characterized by elemental analysis, magnetic, IR and H-1 NMR spectral data. The electronic and magnetic moment 2.97-3.07 B.M. data infers octahedral geometry for all the complexes. The IR data reveals that Schiff base (nih) form coordination bond with the metal ion through azomethine-nitrogen, phenolic-oxygen and carbonyl-oxygen in a tridentate fashion. In addition, DNA-binding properties of these six metal complexes were investigated using absorption spectroscopy, viscosity measurements and thermal denaturation methods. The results indicated that the nickel(II) complex strongly bind with calf-thymus DNA with intrinsic DNA binding constant K-b value of 4.9 x 10(4) M-1 for (3), 4.2 x 10(4) M-1 for (4), presumably via an intercalation mechanism compared to cobalt(III) complex with K-b value of 4.6 x 10(4) M-1 (1) and 4.1 x 10(4) M-1 (2). The DNA Photoclevage experiment shows that, the complexes act as effective DNA cleavage agent. (C) 2012 Elsevier B.V. All rights reserved.

Can the viscosity in astrophysical black hole accretion disks be close to its string theory bound?

String theory and gauge/gravity duality suggest the lower bound of shear viscosity (eta) to entropy density (s) for any matter to be mu h/4 pi k(B), when h and k(B) are reduced Planck and Boltzmann constants respectively and mu <= 1. Motivated by this, we explore eta/s in black hole accretion flows, in order to understand if such exotic flows could be a natural site for the lowest eta/s. Accretion flow plays an important role in black hole physics in identifying the existence of the underlying black hole. This is a rotating shear flow with insignificant molecular viscosity, which could however have a significant turbulent viscosity, generating transport, heat and hence entropy in the flow. However, in presence of strong magnetic field, magnetic stresses can help in transporting matter independent of viscosity, via celebrated Blandford-Payne mechanism. In such cases, energy and then entropy produces via Ohmic dissipation. In,addition, certain optically thin, hot, accretion flows, of temperature greater than or similar to 10(9) K, may be favourable for nuclear burning which could generate/absorb huge energy, much higher than that in a star. We find that eta/s in accretion flows appears to be close to the lower bound suggested by theory, if they are embedded by strong magnetic field or producing nuclear energy, when the source of energy is not viscous effects. A lower bound on eta/s also leads to an upper bound on the Reynolds number of the flow.

Chapman-Enskog analysis for finite-volume formulation of lattice Boltzmann equation

The classical Chapman-Enskog expansion is performed for the recently proposed finite-volume formulation of lattice Boltzmann equation (LBE) method D.V. Patil, K.N. Lakshmisha, Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh, J. Comput. Phys. 228 (2009) 5262-5279]. First, a modified partial differential equation is derived from a numerical approximation of the discrete Boltzmann equation. Then, the multi-scale, small parameter expansion is followed to recover the continuity and the Navier-Stokes (NS) equations with additional error terms. The expression for apparent value of the kinematic viscosity is derived for finite-volume formulation under certain assumptions. The attenuation of a shear wave, Taylor-Green vortex flow and driven channel flow are studied to analyze the apparent viscosity relation.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.

Numerical simulations of CO2 migration during charnockite genesis

Charnockite is considered to be generated either through the dehydration of granitic magma by CO2 purging or by solid-state dehydration through CO2 metasomatism during granulite facies metamorphism. To understand the extent of dehydration, CO2 migration is quantitatively modeled in silicate melt and metasomatic fluid as a function of temperature, H2O wt%, pressure, basal CO2 flux and dynamic viscosity. Numerical simulations show that CO2 advection through porous and permeable high-grade metamorphic rocks can generate dehydrated patches close to the CO2 flow path, as illustrated by the occurrences of ``incipient charnockites.'' CO2 reaction-front velocity constrained by field observations is 0.69 km/m.y., a reasonable value, which matches well with other studies. On the other hand, temperature, rate of cooling, and basal CO2 flux are the critical parameters affecting CO2 diffusion through a silicate melt. CO2 diffusion through silicate melt can only occur at temperature greater than 840 degrees C and during slow cooling (<= 3.7 x 10(-5) degrees C/yr), features that are typical of magma emplacement in the lower crust. Stalling of CO2 fluxing at similar to 840 degrees C explains why some deep-level plutons contain both hydrous and anhydrous (charnockitic) mineral assemblages. CO2 diffusion through silicate melt is virtually insensitive to pressure. Addition of CO2 basal flux facilitates episodic dehydrated melt migration by generating fracture pathways.