On the solution of bivariate population balance equations for aggregation: X-discretization of space for expansion and contraction of computational domain

A new structured discretization of 2D space, named X-discretization, is proposed to solve bivariate population balance equations using the framework of minimal internal consistency of discretization of Chakraborty and Kumar [2007, A new framework for solution of multidimensional population balance equations. Chem. Eng. Sci. 62, 4112-4125] for breakup and aggregation of particles. The 2D space of particle constituents (internal attributes) is discretized into bins by using arbitrarily spaced constant composition radial lines and constant mass lines of slope -1. The quadrilaterals are triangulated by using straight lines pointing towards the mean composition line. The monotonicity of the new discretization makes is quite easy to implement, like a rectangular grid but with significantly reduced numerical dispersion. We use the new discretization of space to automate the expansion and contraction of the computational domain for the aggregation process, corresponding to the formation of larger particles and the disappearance of smaller particles by adding and removing the constant mass lines at the boundaries. The results show that the predictions of particle size distribution on fixed X-grid are in better agreement with the analytical solution than those obtained with the earlier techniques. The simulations carried out with expansion and/or contraction of the computational domain as population evolves show that the proposed strategy of evolving the computational domain with the aggregation process brings down the computational effort quite substantially; larger the extent of evolution, greater is the reduction in computational effort. (C) 2011 Elsevier Ltd. All rights reserved.

An operator-splitting finite element method for the efficient parallel solution of multidimensional population balance systems

A finite element method for solving multidimensional population balance systems is proposed where the balance of fluid velocity, temperature and solute partial density is considered as a two-dimensional system and the balance of particle size distribution as a three-dimensional one. The method is based on a dimensional splitting into physical space and internal property variables. In addition, the operator splitting allows to decouple the equations for temperature, solute partial density and particle size distribution. Further, a nodal point based parallel finite element algorithm for multi-dimensional population balance systems is presented. The method is applied to study a crystallization process assuming, for simplicity, a size independent growth rate and neglecting agglomeration and breakage of particles. Simulations for different wall temperatures are performed to show the effect of cooling on the crystal growth. Although the method is described in detail only for the case of d=2 space and s=1 internal property variables it has the potential to be extendable to d+s variables, d=2, 3 and s >= 1. (C) 2011 Elsevier Ltd. All rights reserved.

A new 6-component accelerometer force balance for short duration ground testing facilities

A new six-component accelerometer force balance is developed and used in the HST2 shock tunnel of Indian Institute of Science. Aerodynamic forces and moments for a hypersonic slender body measured using this balance system at a free stream Mach number of 5.75 and Reynolds number of 1.5 million and stagnation enthalpy of 1.5 and 2 MJ/kg are presented. These measured values compare well with the theoretical values estimated using modified Newtonian theory.

A novel free floating accelerometer force balance system for shock tunnel applications

In order to overcome the interference of the model mounting system with the external aerodynamics of the body during shock tunnel testing, a new free floating internally mountable balance system that ensures unrestrained model motion during testing has been designed, fabricated and tested. Minimal friction ball bearings are used for ensuring the free floating condition of the model during tunnel testing. The drag force acting on a blunt leading edge flat plate at hypersonic Mach number has been measured using the new balance system. Finite element model (FEM) and CFD are exhaustively used in the design as well as for calibrating the new balance system. The experimentally measured drag force on the blunt leading edge flat plate at stagnation enthalpy of 0.7 and 1.2 MJ/kg and nominal Mach number of 5.75 matches well with FEM results. The concept can also be extended for measuring all the three fundamental aerodynamic forces in short duration test facilities like free piston driven shock tunnels.

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.

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]

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.

Low-thermal-budget solution processing of thin films of zinc ferrite and other complex oxides

Further miniaturization of magnetic and electronic devices demands thin films of advanced nanomaterials with unique properties. Spinel ferrites have been studied extensively owing to their interesting magnetic and electrical properties coupled with stability against oxidation. Being an important ferrospinel, zinc ferrite has wide applications in the biological (MRI) and electronics (RF-CMOS) arenas. The performance of an oxide like ZnFe2O4 depends on stoichiometry (defect structure), and technological applications require thin films of high density, low porosity and controlled microstructure, which depend on the preparation process. While there are many methods for the synthesis of polycrystalline ZnFe2O4 powder, few methods exist for the deposition of its thin films, where prolonged processing at elevated temperature is not required. We report a novel, microwave-assisted, low temperature (<100°C) deposition process that is conducted in the liquid medium, developed for obtaining high quality, polycrystalline ZnFe2O4 thin films on technologically important substrates like Si(100). An environment-friendly solvent (ethanol) and non-hazardous oxide precursors (β-diketonates of Zn and Fe in 1:2 molar ratio), forming a solution together, is subjected to irradiation in a domestic microwave oven (2.45 GHz) for a few minutes, leading to reactions which result in the deposition of ZnFe2O4 films on Si (100) substrates suspended in the solution. Selected surfactants added to the reactant solution in optimum concentration can be used to control film microstructure. The nominal temperature of the irradiated solution, i.e., film deposition temperature, seldom exceeds 100°C, thus sharply lowering the thermal budget. Surface roughness and uniformity of large area depositions (50x50 mm2) are controlled by tweaking the concentration of the mother solution. Thickness of the films thus grown on Si (100) within 5 min of microwave irradiation can be as high as several microns. The present process, not requiring a vacuum system, carries a very low thermal budget and, together with a proper choice of solvents, is compatible with CMOS integration. This novel solution-based process for depositing highly resistive, adherent, smooth ferrimagnetic films on Si (100) is promising to RF engineers for the fabrication of passive circuit components. It is readily extended to a wide variety of functional oxide films.

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.

Vegetation impact on stream chemical fluxes: Mule Hole watershed (South India)

The proportion of chemical elements passing through vegetation prior to being exported in a stream was quantified for a forested tropical watershed(Mule Hole, South India) using an extensive hydrological and geochemical monitoring at several scales. First, a solute annual mass balance was established at the scale of the soil-plant profile for assessing the contribution of canopy interaction and litter decay to the solute fluxes of soil inputs (overland flow) and soil outputs (pore water flow as seepages). Second, based on the respective contributions of overland flow and seepages to the stream flow as estimated by a hydrological lumped model, we assigned the proportion of chemical elements in the stream that transited through the vegetation at both flood event (End Member Mixing Analysis) and seasonal scales. At the scale of the 1D soil-plant profile, leaching from the canopy constituted the main source of K above the ground surface. Litter decay was the main source of Si, whereas alkalinity, Ca and Mg originated in the same proportions from both sources. The contribution of vegetation was negligible for Na. Within the soil, all elements but Na were removed from the pore water in proportions varying from 20% for Cl to 95% for K: The soil output fluxes corresponded to a residual fraction of the infiltration fluxes. The behavior of K, Cl, Ca and Mg in the soil-plant profile can be explained by internal cycling, as their soil output fluxes were similar to the atmospheric inputs. Na was released from soils as a result of Na-plagioclase weathering and accompanied by additional release of Si. Concentration of soil pore water by evapotranspiration might limit the chemical weathering in the soil. Overall, the solute K, Ca, Mg, alkalinity and Si fluxes associated with the vegetation turnover within the small experimental watershed represented 10-15 times the solute fluxes exported by the stream, of which 83-97% transited through the vegetation. One important finding is that alkalinity and Si fluxes at the outlet were not linked to the ``current weathering'' of silicates in this watershed. These results highlight the dual effect of the vegetation cover on the solute fluxes exported from the watershed: On one hand the runoff was limited by evapotranspiration and represented only 10% of the annual rainfall, while on the other hand, 80-90% of the overall solute flux exported by the stream transited through the vegetation. The approach combining geochemical monitoring and accurate knowledge of the watershed hydrological budget provided detailed understanding of several effects of vegetation on stream fluxes: (1) evapotranspiration (limiting), (2) vertical transfer through vegetation from vadose zone to ground surface (enhancing) and (3) redistribution by throughfalls and litter decay. It provides a good basis for calibrating geochemical models and more precisely assessing the role of vegetation on soil processes. (C) 2014 Elsevier Ltd. All rights reserved.

High-conductance states and A-type K+ channels are potential regulators of the conductance-current balance triggered by HCN channels

An increase in the hyperpolarization-activated cyclic nucleotide-gated (HCN) channel conductance reduces input resistance, whereas the consequent increase in the inward h current depolarizes the membrane. This results in a delicate and unique conductance-current balance triggered by the expression of HCN channels. In this study, we employ experimentally constrained, morphologically realistic, conductance-based models of hippocampal neurons to explore certain aspects of this conductance-current balance. First, we found that the inclusion of an experimentally determined gradient in A-type K+ conductance, but not in M-type K+ conductance, tilts the HCN conductance-current balance heavily in favor of conductance, thereby exerting an overall restorative influence on neural excitability. Next, motivated by the well-established modulation of neuronal excitability by synaptically driven high-conductance states observed under in vivo conditions, we inserted thousands of excitatory and inhibitory synapses with different somatodendritic distributions. We measured the efficacy of HCN channels, independently and in conjunction with other channels, in altering resting membrane potential (RMP) and input resistance (R-in) when the neuron received randomized or rhythmic synaptic bombardments through variable numbers of synaptic inputs. We found that the impact of HCN channels on average RMP, R in, firing frequency, and peak-to-peak voltage response was severely weakened under high-conductance states, with the impinging synaptic drive playing a dominant role in regulating these measurements. Our results suggest that the debate on the role of HCN channels in altering excitability should encompass physiological and pathophysiological neuronal states under in vivo conditions and the spatiotemporal interactions of HCN channels with other channels.

A comparative study of a direct discretization and an operator-splitting solver for population balance systems

A direct discretization approach and an operator-splitting scheme are applied for the numerical simulation of a population balance system which models the synthesis of urea with a uni-variate population. The problem is formulated in axisymmetric form and the setup is chosen such that a steady state is reached. Both solvers are assessed with respect to the accuracy of the results, where experimental data are used for comparison, and the efficiency of the simulations. Depending on the goal of simulations, to track the evolution of the process accurately or to reach the steady state fast, recommendations for the choice of the solver are given. (C) 2015 Elsevier Ltd. All rights reserved.

Sources of black carbon aerosols in South Asia and surrounding regions during the Integrated Campaign for Aerosols, Gases and Radiation Budget (ICARB)

This study examines differences in the surface black carbon (BC) aerosol loading between the Bay of Bengal (BoB) and the Arabian Sea (AS) and identifies dominant sources of BC in South Asia and surrounding regions during March-May 2006 (Integrated Campaign for Aerosols, Gases and Radiation Budget, ICARB) period. A total of 13 BC tracers are introduced in the Weather Research and Forecasting Model coupled with Chemistry to address these objectives. The model reproduced the temporal and spatial variability of BC distribution observed over the AS and the BoB during the ICARB ship cruise and captured spatial variability at the inland sites. In general, the model underestimates the observed BC mass concentrations. However, the model-observation discrepancy in this study is smaller compared to previous studies. Model results show that ICARB measurements were fairly well representative of the AS and the BoB during the pre-monsoon season. Elevated BC mass concentrations in the BoB are due to 5 times stronger influence of anthropogenic emissions on the BoB compared to the AS. Biomass burning in Burma also affects the BoB much more strongly than the AS. Results show that anthropogenic and biomass burning emissions, respectively, accounted for 60 and 37% of the average +/- standard deviation (representing spatial and temporal variability) BC mass concentration (1341 +/- 2353 ng m(-3)) in South Asia. BC emissions from residential (61 %) and industrial (23 %) sectors are the major anthropogenic sources, except in the Himalayas where vehicular emissions dominate. We find that regional-scale transport of anthropogenic emissions contributes up to 25% of BC mass concentrations in western and eastern India, suggesting that surface BC mass concentrations cannot be linked directly to the local emissions in different regions of South Asia.