Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Lukacova-Medvid'ova,J. Saibertov'a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533-562; M. Luacova-Medvid'ova, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1-30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor-corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.

Prediction of compressive strength of SCC and HPC with high volume fly ash using ANN

An artificial neural network (ANN) is presented to predict a 28-day compressive strength of a normal and high strength self compacting concrete (SCC) and high performance concrete (HPC) with high volume fly ash. The ANN is trained by the data available in literature on normal volume fly ash because data on SCC with high volume fly ash is not available in sufficient quantity. Further, while predicting the strength of HPC the same data meant for SCC has been used to train in order to economise on computational effort. The compressive strengths of SCC and HPC as well as slump flow of SCC estimated by the proposed neural network are validated by experimental results.

The Milky Way's Dark Matter Halo Appears To Be Lopsided

The atomic hydrogen gas (H I) disk in the outer region (beyond similar to 10 kpc from the center) of Milky Way can provide valuable information about the structure of the dark matter halo. The recent three-dimensional thickness map of the outer H I disk from the all sky 21 cm line Leiden/Argentine/Bonn survey, gives us a unique opportunity to investigate the structure of the dark matter halo of Milky Way in great detail. A striking feature of this new survey is the north-south (N-S) asymmetry in the thickness map of the atomic hydrogen gas. Assuming vertical hydrostatic equilibrium under the total potential of the Galaxy, we derive the model thickness map of the H I gas. We show that simple axisymmetric halo models, such as softened isothermal halo (producing a flat rotation curve with V-c similar to 220 km s(-1)) or any halo with density falling faster than the isothermal one, are not able to explain the observed radial variation of the gas thickness. We also show that such axisymmetric halos along with different H I velocity dispersion in the two halves, cannot explain the observed asymmetry in the thickness map. Amongst the nonaxisymmetric models, it is shown that a purely lopsided (m = 1, first harmonic) dark matter halo with reasonable H I velocity dispersion fails to explain the N-S asymmetry satisfactorily. However, we show that by superposing a second harmonic (m = 2) out of phase onto a purely lopsided halo, e. g., our best fit and more acceptable model A (with parameters epsilon(1)(h) = 0.2, epsilon(2)(h) = 0.18, and sigma(H I) = 8.5 km s(-1)) can provide an excellent fit to the observation and reproduce the N-S asymmetry naturally. The emerging picture of the asymmetric dark matter halo is supported by the. cold dark matter halos formed in the cosmological N-body simulation.

Role of free volume in strain softening of as-cast and annealed bulk metallic glass

Plasticity in amorphous alloys is associated with strain softening, induced by the creation of additional free volume during deformation. In this paper, the role of free volume, which was a priori in the material, on work softening was investigated. For this, an as-cast Zr-based bulk metallic glass (BMG) was systematically annealed below its glass transition temperature, so as to reduce the free volume content. The bonded-interface indentation technique is used to generate extensively deformed and well defined plastic zones. Nanoindentation was utilized to estimate the hardness of the deformed as well as undeformed regions. The results show that the structural relaxation annealing enhances the hardness and that both the subsurface shear band number density and the plastic zone size decrease with annealing time. The serrations in the nanoindentation load-displacement curves become smoother with structural relaxation. Regardless of the annealing condition, the nanohardness of the deformed regions is similar to 12-15% lower, implying that the prior free volume only changes the yield stress (or hardness) but not the relative flow stress (or the extent of strain softening). Statistical distributions of the nanohardness obtained from deformed and undeformed regions have no overlap, suggesting that shear band number density has no influence on the plastic characteristics of the deformed region.

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Critical Evaluation of Determining Swelling Pressure by Swell-Load Method and Constant Volume Method

For any construction activity in expansive soils, determination of swelling pressure/heave is an essential step. Though many attempts have been made to develop laboratory procedures by using the laboratory one-dimensional oedometer to determine swelling pressure of expansive soils, they are reported to yield varying results. The main reason for these variations could be heterogeneous moisture distribution of the sample over its thickness. To overcome this variation the experimental procedure should be such that the soil gets fully saturated. Attempts were made to introduce vertical sand drains in addition to the top and bottom drains. In this study five and nine vertical sand drains were introduced to experimentally find out the variations in the swell and swelling pressure. The variations in the moisture content at middle, top, and bottom of the sample in the oedometer test are also reported. It is found that swell-load method is better as compared to zero-swell method. Further, five number of vertical sand drains are found to be sufficient to obtain uniform moisture content distribution.

Magneto-structural study of a Cr-doped CaRuO3

Temperature-dependent neutron powder diffraction, magnetization and XPS studies were carried out on an optimally Cr-doped CaRuO3, i.e. CaRu0.85Cr0.15O3 (CRC-15). XPS data revealed that Cr exist in 3+ and 6+ oxidation states. The charge dissociation preserves the overall 4+ nominal charge of the Ru site. Although ferromagnetic correlations develop around 100 K, the system exhibits a large coercive field below 50 K. The unit cell volume exhibits negative thermal expansion below 50 K since the lattice expansion due to the magnetostrictive effect outweighs the thermal contraction due to the phonon-driven mechanism.

Effectiveness factor for zeroth-order catalytic reactions with volume change and nonuniform catalyst activity profiles

Analytical solutions are presented for the effectiveness factor of a zeroth-order reaction with volume change and nonuniform catalyst activity profile in slab, cylinder and spherical pellets. The possibility of shape normalization is considered for a variety of activity profiles and pellet shapes. When the catalyst activity at the external surface of the pellet is non-zero, shape normalization is obtained, which makes the asymptotic behavior of the effectiveness factor identical for small and large values of Thiele modulus, however, the normalization can lead to significant errors, particularly for the case of activity profiles decreasing towards the outer surface of the catalyst.

R-parameter: A local truncation error based adaptive framework for finite volume compressible flow solvers

A residual-based strategy to estimate the local truncation error in a finite volume framework for steady compressible flows is proposed. This estimator, referred to as the -parameter, is derived from the imbalance arising from the use of an exact operator on the numerical solution for conservation laws. The behaviour of the residual estimator for linear and non-linear hyperbolic problems is systematically analysed. The relationship of the residual to the global error is also studied. The -parameter is used to derive a target length scale and consequently devise a suitable criterion for refinement/derefinement. This strategy, devoid of any user-defined parameters, is validated using two standard test cases involving smooth flows. A hybrid adaptive strategy based on both the error indicators and the -parameter, for flows involving shocks is also developed. Numerical studies on several compressible flow cases show that the adaptive algorithm performs excellently well in both two and three dimensions.

Rocket Propellant Combustion Studies in a Constant Volume Bomb

A constant volume window bomb has been used to measure the characteristic velocity (c*) of rocket propellants. Analysis of the combustion process inside the bomb including heat losses has been made. The experiments on double base and composite propellants have revealed some (i) basic heat transfer aspects inside the bomb and (ii) combustion characteristics of Ammonium Perchlorate-Polyester propellants. It has been found that combustion continues even beyond the peak pressure and temperature points. Lithium Fluoride mixed propellants do not seem to indicate significant differences in c*) though the low pressure deflagration limit is increased with percentage of Lithium Fluoride.

Phase diagram of vortex matter in layered superconductors with tilted columnar pinning centers

We study the vortex matter phase diagram of a layered superconductor in the presence of columnar pinning defects, tilted with respect to the normal to the layers. We use numerical minimization of the free energy written as a functional of the time-averaged vortex density of the Ramakrishnan-Yussouff form, supplemented by the appropriate pinning potential. We study the case where the pin density is smaller than the areal vortex density. At lower pin concentrations, we find, for temperatures of the order of the melting temperature of the unpinned lattice, a Bose glass type phase which at lower temperatures converts, via a first-order transition, to a Bragg glass, while, at higher temperatures, it crosses over to an interstitial liquid. At somewhat higher concentrations, no transition to a Bragg glass is found even at the lowest temperatures studied. While qualitatively the behavior we find is similar to that obtained using the same procedures for columnar pins normal to the layers, there are important and observable quantitative differences, which we discuss.

Conserving Biodiversity as If People Matter: A Case Study from India

India has rich traditions of nature conservation as well as a vigorous official program of protection of nature reserves developed over the last 40 years. However, the officialp rograms uffers fromt otal relianceo n authoritarianm anagement arrangements in which decisions are made centrally and coercion is used to implement them. At the same time, the state apparatus organises subsidized resource flows to the urbanindustrial- intensivea griculturalc omplex which promote inefficient,n on-sustainable resource-use patterns that are inimical to conservation of biodiversity. These processes are illustrated within the concrete setting of the district of Uttara Kannada in southern India. It is suggested that the interests of conservation would be served far better by an approach that withdraws the subsidies to the elite so that a much more efficient, sustainable and equitable pattern of resource use, compatible with conservation of biodiversity, can be instituted. In conjunction with this, the larger society should involve local people in working out detailed plans for conservation of biodiversity and offer them adequate authority as well as appropriate financial incentives to implement these plans. The paper goes on to illustrate how such an approach may be implemented in the case of Uttara Kannada.

Dark matter dominance at all radii in the superthin galaxy UGC 7321

We model the shape and density profile of the dark matter halo of the low surface brightness, superthin galaxy UGC 7321, using the observed rotation curve and the H i scale height data as simultaneous constraints. We treat the galaxy as a gravitationally coupled system of stars and gas, responding to the gravitational potential of the dark matter halo. An isothermal halo of spherical shape with a core density in the range of View the MathML source and a core radius between 2.5 and 2.9 kpc, gives the best fit to the observations for a range of realistic gas parameters assumed. We find that the best-fit core radius is only slightly higher than the stellar disc scale length (2.1 kpc), unlike the case of the high surface brightness galaxies where the halo core radius is typically 3–4 times the disc scale length of the stars. Thus our model shows that the dark matter halo dominates the dynamics of the low surface brightness, superthin galaxy UGC 7321 at all radii, including the inner parts of the galaxy.

Soft matter lithium salt electrolytes: ion conduction and application to rechargeable batteries

Soft matter provides diverse opportunities for the development of electrolytes for all solid state lithium batteries. Here we review soft matter solid electrolytes for lithium batteriesthat are primarily obtained starting from liquid electrolytic systems. This concept of solid electrolyte synthesis from liquid is significantly different from prevalent approaches. The novelty of our approach is discussed in the light of various fundamental issues and in relation to its application to rechargeable lithium batteries.