Composite scheme using localized relaxation with non-standard finite difference method for hyperbolic conservation laws

Non-standard finite difference methods (NSFDM) introduced by Mickens [Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers–Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791–797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250–2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235–276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter (λ) is chosen locally on the three point stencil of grid which makes the proposed method more efficient. This composite scheme overcomes the problem of unphysical expansion shocks and captures the shock waves with an accuracy better than the upwind relaxation scheme, as demonstrated by the test cases, together with comparisons with popular numerical methods like Roe scheme and ENO schemes.

Nanoindentation behavior of nanotwinned Cu: Influence of indenter angle on hardness, strain rate sensitivity and activation volume

The influence of strain on the mechanical properties and deformation kinetic parameters of nanotwinned (at) copper is investigated by a series of nanoindentation experiments, which were performed by employing sharp indenters with five varying centerline-to-face angles (psi). Comparison experiments were also conducted on (1 1 0) single crystalline Cu. Experimental results indicate that, unlike coarsegrained materials, nt-Cu is prone to plastic flow softening with large material pile-up around the indentation impression at high levels of strains. Localized detwinning becomes more significant with decreasing psi, concomitant with reduced strain-rate sensitivity (m) and enhanced activation volume (V*). The m of nt-Cu is found to depend sensitively on psi with a variation of more than a factor of 3, whereas V* exhibits a much less sensitive trend. This paper discusses the validation of the experimental techniques and the implications of various deformation kinetic parameters on the underlying deformation mechanisms of nt-Ca. 2013 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Positive temperature coefficient and structural relaxations in selectively localized MWNTs in PE/PEO blends

The dispersion state of multiwall carbon nanotubes (MWNTs) in melt mixed polyethylene/polyethylene oxide (PE/PEO) blends has been assessed by both surface and volume electrical conductivity measurements and the structural relaxations have been assessed by broadband dielectric spectroscopy. The selective localization of MWNTs in the blends was controlled by the flow characteristics of the components, which led to their localization in the energetically less favored phase (PE). The electrical conductivity and positive temperature co-efficient (PTC) measurements were carried out on hot pressed samples. The neat blends exhibited only a negative temperature coefficient (NTC) effect while the blends with MWNTs exhibited both a PTC and a NTC at the melting temperatures of PE and PEO respectively. These phenomenal changes were corroborated with the different crystalline morphology in the blends. It was deduced that during compression molding, the more viscous PEO phase spreads less in contrast to the less viscous PE phase. This has further resulted in a gradient in morphology as well as the distribution state of the MWNTs in the samples and was supported by scanning electron and scanning acoustic microscopy (SAM) studies and contact angle measurements. SAM from different depths of the samples revealed a gradient in the microstructure in the PE/PEO blends which is contingent upon the flow characteristics of the components. Interestingly, the surface and volume electrical conductivity was different due to the different dispersion state of the MWNTs at the surface and bulk. The observed surface and volume electrical conductivity measurements were corroborated with the evolved morphology during processing. The structural relaxations in both PE and PEO were discerned from broadband dielectric spectroscopy. The segmental dynamics below and above the melting temperature of PEO were significantly different in the presence of MWNTs.

Diffusing-wave spectroscopy in an inhomogeneous object: Local viscoelastic spectra from ultrasound-assisted measurement of correlation decay arising from the ultrasound focal volume

We demonstrate diffusing-wave spectroscopy (DWS) in a localized region of a viscoelastically inhomogeneous object by measurement of the intensity autocorrelation g(2)(tau)] that captures only the decay introduced by the temperature-induced Brownian motion in the region. The region is roughly specified by the focal volume of an ultrasound transducer which introduces region specific mechanical vibration owing to insonification. Essential characteristics of the localized non-Markovian dynamics are contained in the decay of the modulation depth M(tau)], introduced by the ultrasound forcing in the focal volume selected, on g(2)(tau). The modulation depth M(tau(i)) at any delay time tau(i) can be measured by short-time Fourier transform of g(2)(tau) and measurement of the magnitude of the spectrum at the ultrasound drive frequency. By following the established theoretical framework of DWS, we are able to connect the decay in M(tau) to the mean-squared displacement (MSD) of scattering centers and the MSD to G*(omega), the complex viscoelastic spectrum. A two-region composite polyvinyl alcohol phantom with different viscoelastic properties is selected for demonstrating local DWS-based recovery of G*(omega) corresponding to these regions from the measured region specific M(tau(i))vs tau(i). The ultrasound-assisted measurement of MSD is verified by simulating, using a generalized Langevin equation (GLE), the dynamics of the particles in the region selected as well as by the usual DWS experiment without the ultrasound. It is shown that whereas the MSD obtained by solving the GLE without the ultrasound forcing agreed with its experimental counterpart covering small and large values of tau, the match was good only in the initial transients in regard to experimental measurements with ultrasound.

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.

Magnetic Field Induced Instabilities In Localized Two-Dimensional Electron Systems

The results of extensive transport studies in localized regime of mesoscopic two-dimensional electron systems (2DES) with varying disorder are presented. A quick overview of previously achieved result is given. The main focus is on the observation of density dependent instabilities manifested by strong resistance oscillations induced by high perpendicular magnetic fields B-perpendicular to. While the amplitude of the oscillations is strongly enhanced with increasing B-perpendicular to, their position in electron density remains unaffected. The temperature dependence of resistivity shows a transition from an activated behaviour at high temperature to a saturated behaviour at low T. In the positions of resistance minima, the T dependence can even become metal-like (d rho/dT > 0). The activation energies obtained from the high T behaviour exhibit a formation of plateaux in connection with the resistance oscillations when analyzed as a function of electron density. We suggest the interplay between a strongly interacting electron phase and the background disorder as a possible explanation for our observation.

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.

Mordenite-incorporated PVA-PSSA membranes as electrolytes for DMFCs

Composite membranes with mordenite (MOR) incorporated in poly vinyl alcohol (PVA)–polystyrene sulfonic acid (PSSA) blend tailored with varying degree of sulfonation are reported. Such a membrane comprises a dispersed phase of mordenite and a continuous phase of the polymer that help tuning the flow of methanol and water across it. The membranes on prolonged testing in a direct methanol fuel cell (DMFC) exhibit mitigated methanol cross-over from anode to the cathode. The membranes have been tested for their sorption behaviour, ion-exchange capacity, electrochemical selectivity and mechanical strength as also characterized by Fourier transform infrared spectroscopy and thermogravimetric analysis. Water release kinetics has been measured by magnetic resonance imaging (NMR imaging) and is found to be in agreement with the sorption data. Similarly, methanol release kinetics studied by volume-localized NMR spectroscopy (point resolved spectroscopy, PRESS) clearly demonstrates that the dispersion of mordenite in PVA–PSSA retards the methanol release kinetics considerably. A peak power-density of 74 mW/cm2 is achieved for the DMFC using a PVA–PSSA membrane electrolyte with 50% degree of sulfonation and 10 wt.% dispersed mordenite phase. A methanol cross-over current as low as 7.5 mA/cm2 with 2 M methanol feed at the DMFC anode is observed while using the optimized composite membrane as electrolyte in the DMFC, which is about 60% and 46% lower than Nafion-117 and PVA–PSSA membranes, respectively, when tested under identical conditions.

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.

Metal-insulator transition of naxwo3 studied by angle-resolved photoemission spectroscopy

The electronic structure of sodium tungsten bronzes NaxWO3 is investigated by high-resolution angle-resolved photoemission spectroscopy (ARPES). The ARPES spectra measured in both insulating and metallic phases of NaxWO3 reveals the origin of metal-insulator transition (MIT) in sodium tungsten bronze system. It is found that in insulating NaxWO3 the states near the Fermi level (E-F) are localized due to the strong disorder caused by the random distribution of Na+ ions in WO3 lattice. Due to the presence of disorder and long-range Coulomb interaction of conduction electrons, a soft Coulomb gap arises, where the density of states vanishes exactly at E-F. In the metallic regime the states near E-F are populated and the Fermi level shifts upward rigidly with increasing electron doping (x). Volume of electron-like Fermi surface (FS) at the Gamma(X) point of the Brillouin zone gradually increases with increasing Na concentration due to W 5d t(2g) band filling. A rigid shift of the Fermi energy is found to give a qualitatively good description of the Fermi surface evolution. As we move from bulk-sensitive to more surface sensitive photon energy, we found the emergence of Fermi surfaces at X(M) and M(R) point similar to the one at the Gamma(X) point in the metallic regime, suggesting that the reconstruction of surface was due to rotation/deformation of WO6 octahedra.