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.

Fuzzy-logic-based health monitoring and residual-life prediction for composite helicopter rotor

A health-monitoring and life-estimation strategy for composite rotor blades is developed in this work. The cross-sectional stiffness reduction obtained by physics-based models is expressed as a function of the life of the structure using a recent phenomenological damage model. This stiffness reduction is further used to study the behavior of measurable system parameters such as blade deflections, loads, and strains of a composite rotor blade in static analysis and forward flight. The simulated measurements are obtained using an aeroelastic analysis of the composite rotor blade based on the finite element in space and time with physics-based damage modes that are then linked to the life consumption of the blade. The model-based measurements are contaminated with noise to simulate real data. Genetic fuzzy systems are developed for global online prediction of physical damage and life consumption using displacement- and force-based measurement deviations between damaged and undamaged conditions. Furthermore, local online prediction of physical damage and life consumption is done using strains measured along the blade length. It is observed that the life consumption in the matrix-cracking zone is about 12-15% and life consumption in debonding/delamination zone is about 45-55% of the total life of the blade. It is also observed that the success rate of the genetic fuzzy systems depends upon the number of measurements, type of measurements and training, and the testing noise level. The genetic fuzzy systems work quite well with noisy data and are recommended for online structural health monitoring of composite helicopter rotor blades.

Collapse behaviour of a residual soil

Residually derived red soils occur in Bangalore District of Karnataka State, India. The porous and unsaturated nature of the red soils makes them susceptible to collapse on wetting under load. The present study analyses the collapse behaviour of an unsaturated bonded (undisturbed) red soil from Bangalore referenced to tests on samples in an unbonded (remoulded) state. A filter paper method was used to determine the matric suction of the bonded and unbonded specimens, and mercury intrusion porosimetry (MIP) was used to determine their soil structure. Analysis of the experimental results shows that bonding plays an important role in the collapse behaviour of the unsaturated residual soil. The results of the study also provide insight into the volume change behaviour of unsaturated bonded soils on wetting within and beyond the yield locus.

A h-Adaptive Algorithm Using Residual Error Estimates for Fluid Flows

Algorithms for adaptive mesh refinement using a residual error estimator are proposed for fluid flow problems in a finite volume framework. The residual error estimator, referred to as the R-parameter is used to derive refinement and coarsening criteria for the adaptive algorithms. An adaptive strategy based on the R-parameter is proposed for continuous flows, while a hybrid adaptive algorithm employing a combination of error indicators and the R-parameter is developed for discontinuous flows. Numerical experiments for inviscid and viscous flows on different grid topologies demonstrate the effectiveness of the proposed algorithms on arbitrary polygonal grids.

Influence of surface area to volume ratio of fuel particles on gasification process in a fixed bed

The paper addresses the effect of particle size on tar generation in a fixed bed gasification system. Pyrolysis, a diffusion limited process, depends on the heating rate and the surface area of the particle influencing the release of the volatile fraction leaving behind residual char. The flaming time has been estimated for different biomass samples. It is found that the flaming time for wood flakes is almost one fourth than that of coconut shells for same equivalent diameter fuel samples. The particle density of the coconut shell is more than twice that of wood spheres, and almost four times compared with wood flakes; having a significant influence on the flaming time. The ratio of the particle surface area to that of an equivalent diameter is nearly two times higher for flakes compared with wood pieces. Accounting for the density effect, on normalizing with density of the particle, the flaming rate is double in the case of wood flakes or coconut shells compared with the wood sphere for an equivalent diameter. This is due to increased surface area per unit volume of the particle. Experiments are conducted on estimation of tar content in the raw gas for wood flakes and standard wood pieces. It is observed that the tar level in the raw gas is about 80% higher in the case of wood flakes compared with wood pieces. The analysis suggests that the time for pyrolysis is lower with a higher surface area particle and is subjected to fast pyrolysis process resulting in higher tar fraction with low char yield. Increased residence time with staged air flow has a better control on residence time and lower tar in the raw gas. (C) 2014 International Energy Initiative. Published by Elsevier Inc. All rights reserved.

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.

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.

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.

Effect of carbon equivalent and inoculation on residual stresses in grey iron castings

It is virtually impossible to produce castings free from internal stresses using conventional methods of founding. Castings with appreciable stresses distort during storage, transportation, machining and service. Though composition and melt treatment are known to affect the magnitude of residual stress in castings, the data on the effect of carbon equivalent and inoculation on the magnitude of residual stress in castings are limited. In the present investigation, an attempt is made to study (i) the effect of carbon equivalent on residual stress in cast iron castings, and (ii) the effect of inoculants such as calcium silicide and ferrosilicon on residual stress in iron castings in the carbon equivalent range 3.0–4.0%. The results of the investigation indicate the following: (i) the residual strains decrease linearly with increase in carbon equivalent in the uninoculated and inoculated irons; (ii) the tensile residual stresses decrease linearly with increase in carbon equivalent value of the uninoculated, calcium silicide-inoculated and ferrosilicon-inoculated cast iron castings; (iii) the ratio of UTS to residual stress increased on inoculating the grid castings. This increase is higher for calcium silicide-inoculated grids than for ferrosilicon-inoculated grid castings. This implies that from the residual stress point of view, inoculation of the iron with calcium silicide is beneficial.

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.

Effects of mould filling on evolution of the solid-liquid interface during solidification

This work presents a numerical analysis of simultaneous mould filling and phase change for solidification in a two-dimensional rectangular cavity. The role of residual flow strength and temperature gradients within the solidifying domain, caused by the filling process, on the evolution of solidification interface are investigated. An implicit volume of fluid (VOF)-based algorithm has been employed for simulating the free surface flows during the filling process, while the model for solidification is based on a fixed-grid enthalpy-based control volume approach. Solidification modeling is coupled with VOF through User Defined Functions developed in the commercial computational fluid dynamics (CFD) code FLUENT 6.3.26. Comparison between results of the conventional analysis without filling effect and those of the present analysis shows that the residual flow resulting from the filling process significantly influences the progress of the solidification interface. A parametric study is also performed with variables such as cooling rate, filling velocity and filling configuration, in order to investigate the coupled effects of the buoyancy-driven flow and the residual flow on the solidification behavior.

Effect of metal coating and residual stress on the resonant frequency of MEMS resonators

MEMS resonators are designed for a fixed resonant frequency. Therefore, any shift in the resonant frequency of the final fabricated structure can be a denting factor for its suitability towards a desired application. There are numerous factors which alter the designed resonant frequency of the fabricated resonator such as the metal layer deposited on top of the beam and the residual stresses present in the fabricated structure. While the metal coating, which acts as electrode, increases the stiffness and the effective mass of the composite structure, the residual stress increases or decreases the net stiffness if it is a tensile or compressive type respectively. In this paper, we investigate both these cases by taking two different structures, namely, the micro cantilever beam with gold layer deposited on its top surface and the MEMS gyroscope with residual stresses. First, we carry out experiments to characterize both these structures to find their resonant frequencies. Later, we analytically model those effects and compare them with the experimentally obtained values. Finally, it is found that the analytical models give an error of less than 10% with respect to the experimental results in both the cases.