On the solution of bivariate population balance equations for aggregation: Pivotwise expansion of solution domain

The solution of a bivariate population balance equation (PBE) for aggregation of particles necessitates a large 2-d domain to be covered. A correspondingly large number of discretized equations for particle populations on pivots (representative sizes for bins) are solved, although at the end only a relatively small number of pivots are found to participate in the evolution process. In the present work, we initiate solution of the governing PBE on a small set of pivots that can represent the initial size distribution. New pivots are added to expand the computational domain in directions in which the evolving size distribution advances. A self-sufficient set of rules is developed to automate the addition of pivots, taken from an underlying X-grid formed by intersection of the lines of constant composition and constant particle mass. In order to test the robustness of the rule-set, simulations carried out with pivotwise expansion of X-grid are compared with those obtained using sufficiently large fixed X-grids for a number of composition independent and composition dependent aggregation kernels and initial conditions. The two techniques lead to identical predictions, with the former requiring only a fraction of the computational effort. The rule-set automatically reduces aggregation of particles of same composition to a 1-d problem. A midway change in the direction of expansion of domain, effected by the addition of particles of different mean composition, is captured correctly by the rule-set. The evolving shape of a computational domain carries with it the signature of the aggregation process, which can be insightful in complex and time dependent aggregation conditions. (c) 2012 Elsevier Ltd. All rights reserved.

An operator-splitting Galerkin/SUPG finite element method for population balance equations : stability and convergence

We present a heterogeneous finite element method for the solution of a high-dimensional population balance equation, which depends both the physical and the internal property coordinates. The proposed scheme tackles the two main difficulties in the finite element solution of population balance equation: (i) spatial discretization with the standard finite elements, when the dimension of the equation is more than three, (ii) spurious oscillations in the solution induced by standard Galerkin approximation due to pure advection in the internal property coordinates. The key idea is to split the high-dimensional population balance equation into two low-dimensional equations, and discretize the low-dimensional equations separately. In the proposed splitting scheme, the shape of the physical domain can be arbitrary, and different discretizations can be applied to the low-dimensional equations. In particular, we discretize the physical and internal spaces with the standard Galerkin and Streamline Upwind Petrov Galerkin (SUPG) finite elements, respectively. The stability and error estimates of the Galerkin/SUPG finite element discretization of the population balance equation are derived. It is shown that a slightly more regularity, i.e. the mixed partial derivatives of the solution has to be bounded, is necessary for the optimal order of convergence. Numerical results are presented to support the analysis.

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.

Existence of Value in Stochastic Differential Games of Mixed Type

In this article, we address stochastic differential games of mixed type with both control and stopping times. Under standard assumptions, we show that the value of the game can be characterized as the unique viscosity solution of corresponding Hamilton-Jacobi-Isaacs (HJI) variational inequalities.

Perfect Static Balance of Linkages by Addition of Springs But Not Auxiliary Bodies

A linkage of rigid bodies under gravity loads can be statically counter-balanced by adding compensating gravity loads. Similarly, gravity loads or spring loads can be counterbalanced by adding springs. In the current literature, among the techniques that add springs, some achieve perfect static balance while others achieve only approximate balance. Further, all of them add auxiliary bodies to the linkage in addition to springs. We present a perfect static balancing technique that adds only springs but not auxiliary bodies, in contrast to the existing techniques. This technique can counter-balance both gravity loads and spring loads. The technique requires that every joint that connects two bodies in the linkage be either a revolute joint or a spherical joint. Apart from this, the linkage can have any number of bodies connected in any manner. In order to achieve perfect balance, this technique requires that all the spring loads have the feature of zero-free-length, as is the case with the existing techniques. This requirement is neither impractical nor restrictive since the feature can be practically incorporated into any normal spring either by modifying the spring or by adding another spring in parallel. DOI: 10.1115/1.4006521]

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.

Variation of Snowline and Mass Balance of Glaciers of Warwan and Bhut Basins of Western Himalaya Using Remote Sensing Technique

Glaciers are natural reservoirs of fresh water in frozen state and sensitive indicators of climate change. Among all the mountainous glaciated regions, glaciers of Himalayas form one of the largest concentrations of ice outside the Polar Regions. Almost all the major rivers of northern India originate from these glaciers and sustain perennial flow. Therefore, in view of the importance and role of the glaciers in sustaining the life on the Earth, monitoring the health of glaciers is necessary. Glacier's health is monitored in two ways (i) by mapping the change in extent of glaciers (ii) by finding variation in the annual mass balance. This paper has been discussed the later approach for monitoring the health of glaciers of Warwan and Bhut basins. Mass balance of glaciers of these two basins was determined based on the extraction of snow line at the end of ablation season. A series of satellite images of AWiFS sensor were analysed for extraction of snowline on the glaciers for the period of 2005, 2006 and 2007. The snow line at the end of ablation season is used to compute accumulation area ratio (AAR = Accumulation area/Glacier area) for each glacier of basins. An approach based on relationship of AAR to specific mass balance (computed in field) for glaciers of Basapa basin was employed in the present study. Mean of specific mass balance of individual glacier for the year 2005, 2006 and 2007 of Warwan basin was found to be -ve 0.19 m, -ve 0.27 m and -ve 0.2 m respectively. It is 0.05 m, -ve 0.11 m and -ve 0.19 m for Bhut basin. The analysis suggests a loss of 4.3 and 0.83 kmA(3) of glacier in the monitoring period of 3 years for Warwan and Bhut basins respectively. The overall results suggest that the glaciers of Warwan basin and Bhut basins have suffered more loss of ice than gain in the monitoring period of 3 years.

A Model for the Temperature Distribution in Resin Impregnated Paper Bushings

Resin impregnated paper (RIP) is a relatively new insulation system recommended for the use in transformer bushings. In the recent past, RIP has acquired prominence as insulation in bushings, over conventional oil impregnated paper (OIP), in view of its overwhelming advantages the more important among them being low dielectric loss and possibility for positioning the bushing at any desired angle over the transformer. In addition, the fact that such systems do not pose problems of fire hazard is counted as a very important consideration. The disadvantage of RIP compared to OIP, however, is its much higher cost and involved manufacturing process. The temperature rise in RIP bushings under normal operating conditions is seen to be a difficult parameter to control in view of the limited options for effective cooling. It is therefore essential to take serious note of this aspect, to arrest rapid deterioration of bushing. The degradation of dry-type insulation such as RIP is often due to thermal stress. The long time performance thereof, depends strongly, on the maximum operating temperature. With this in view, the Authors have developed a theoretical model and computational method to study the temperature distribution in the body of insulation. The Authors consider that the basis for the model as being the temperature and electric stress aided AC conductivity. The ensuing heat balance (continuity) equations in 2-D cylindrical geometry are treated as a Dirichelet-Neumann boundary value problem.

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.

Multivalent glycoliposomes and micelles to study carbohydrate-protein and carbohydrate-carbohydrate interactions

This tutorial review describes multivalent carbohydrate-protein and carbohydrate-carbohydrate interaction studies that utilize self-assembled aggregates of thermodynamically stable liposomes and micelles. Strategies to prepare multivalent glycoliposomes and micelles include: (i) insertion of synthetic glycolipids into matrix lipids; (ii) preparation of glycolipids that aggregate to liposomes and micelles and (iii) modification of the hydrophilic surfaces with desired sugars. Several design strategies have been developed in order to obtain constituent glycolipids, having multivalent sugar moieties and their subsequent interactions with proteins were assessed in relation to the type of linkers that connect the hydrophilic and lipophilic segments. Lipophilic segments other than alkyl chains have also been developed. Polymer based glycoliposomes and micelles form an emphasis. Further, glycoliposomes facilitate studies of carbohydrate-carbohydrate interactions. An overview of the various types of glycoliposomes and micelles used to study carbohydrate-protein and carbohydrate-carbohydrate recognition phenomena is presented.

Elasticity of DNA and the effect of dendrimer binding

Negatively charged DNA can be compacted by positively charged dendrimers and the degree of compaction is a delicate balance between the strength of the electrostatic interaction and the elasticity of DNA. We report various elastic properties of short double-stranded DNA (dsDNA) and the effect of dendrimer binding using fully atomistic molecular dynamics and numerical simulations. In equilibrium at room temperature, the contour length distribution P(L) and the end-to-end distance distribution P(R) are nearly Gaussian, the former gives an estimate of the stretch modulus gamma(1) of dsDNA in quantitative agreement with the literature value. The bend angle distribution P(.) of the dsDNA also has a Gaussian form and allows to extract a persistence length, L-p of 43 nm. When the dsDNA is compacted by positively charged dendrimer, the stretch modulus stays invariant but the effective bending rigidity estimated from the end-to-end distance distribution decreases dramatically due to backbone charge neutralization of dsDNA by dendrimer. We support our observations with numerical solutions of the worm-like-chain (WLC) model as well as using non-equilibrium dsDNA stretching simulations. These results are helpful in understanding the dsDNA elasticity at short length scales as well as how the elasticity is modulated when dsDNA binds to a charged object such as a dendrimer or protein.

Induction of Supramolecular Chirality in the Self-Assemblies of Lipophilic Pyrimidine Derivatives by Choice of the Amino Acid-Based Chiral Spacer

A new family of supramolecular organogelators, based on chiral amino acid derivatives of 2,4,6-trichloro-pyrimidine-5-carbaldehyde, has been synthesized. L-alanine was incorporated as a spacer between the pyrimidine core and long hydrocarbon tails to compare the effect of chirality and hydrogen bonding to that of the achiral analogue. The role of aromatic moiety on the chiral spacer was also investigated by introducing L-phenyl alanine moieties. The presence of intermolecular hydrogen-bonding leading to the chiral self-assembly was probed by concentration-dependent FTIR and UV/Vis spectroscopies, in addition to circular dichroism (CD) studies. Temperature and concentration-dependent CD spectroscopy ascribed to the formation of -sheet-type H-bonded networks. The morphology and the arrangements of the molecules in the freeze-dried gels were examined by scanning electron microscopy (SEM), transmission electron microscopy (TEM), atomic force microscopy (AFM), and X-ray diffraction (XRD) techniques. Calculation of the length of each molecular system by energy minimization in its extended conformation and comparison with the small-angle XRD pattern reveals that this class of gelator molecules adopts a lamellar organization. Polarized optical microscopy (POM) and differential scanning calorimetry (DSC) indicate that the solid state phase behavior of these molecules is totally dependent on the choice of their amino acid spacers. Structure-induced aggregation properties based on the H-bonding motifs and the packing of the molecule in three dimensions leading to gelation was elucidated by rheological studies. However, viscoelasticity was shown to depend only marginally on the H-bonding interactions; rather it depends on the packing of the gelators to a greater extent.

Evaluation of document binarization using eigen value decomposition

A necessary step for the recognition of scanned documents is binarization, which is essentially the segmentation of the document. In order to binarize a scanned document, we can find several algorithms in the literature. What is the best binarization result for a given document image? To answer this question, a user needs to check different binarization algorithms for suitability, since different algorithms may work better for different type of documents. Manually choosing the best from a set of binarized documents is time consuming. To automate the selection of the best segmented document, either we need to use ground-truth of the document or propose an evaluation metric. If ground-truth is available, then precision and recall can be used to choose the best binarized document. What is the case, when ground-truth is not available? Can we come up with a metric which evaluates these binarized documents? Hence, we propose a metric to evaluate binarized document images using eigen value decomposition. We have evaluated this measure on DIBCO and H-DIBCO datasets. The proposed method chooses the best binarized document that is close to the ground-truth of the document.

Glacial mass balance changes in the Karakoram and Himalaya based on CMIP5 multi-model climate projections

The impact of future climate change on the glaciers in the Karakoram and Himalaya (KH) is investigated using CMIP5 multi-model temperature and precipitation projections, and a relationship between glacial accumulation-area ratio and mass balance developed for the region based on the last 30 to 40 years of observational data. We estimate that the current glacial mass balance (year 2000) for the entire KH region is -6.6 +/- 1 Gta(-1), which decreases about sixfold to -35 +/- 2 Gta(-1) by the 2080s under the high emission scenario of RCP8.5. However, under the low emission scenario of RCP2.6 the glacial mass loss only doubles to -12 +/- 2 Gta(-1) by the 2080s. We also find that 10.6 and 27 % of the glaciers could face `eventual disappearance' by the end of the century under RCP2.6 and RCP8.5 respectively, underscoring the threat to water resources under high emission scenarios.