Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

Effect of hydrogen on the yielding behavior and shear transformation zone volume in metallic glass ribbons

The effect of hydrogen (H) charging on the shear yield strength (tau(max)) and shear transformation zone volume (Omega) of Ni-Nb-Zr metallic glass ribbons, with varying Zr content, were studied through the first pop-in loads during nanoindentation. Weight gain measurements after H charging and desorption studies were utilized to identify how the total H absorbed during charging is partitioned into mobile and immobile (or trapped) parts. These, in turn, indicate the significant role of H mobility in the amorphous structure on the yielding behavior. In high-Zr alloys, tau(max) increases significantly whereas Omega decreases. In low-Zr alloys, a slight decrease in tau(max) and increase in Omega were noted. These experimental observations are rationalized in terms of the mobility of the absorbed H in the amorphous structure and the possible role of it in the shear transformation zone dynamics during deformation of the metallic glass. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Estimating the volume of glaciers in the Himalayan-Karakoram region using different methods

Ice volume estimates are crucial for assessing water reserves stored in glaciers. Due to its large glacier coverage, such estimates are of particular interest for the Himalayan-Karakoram (HK) region. In this study, different existing methodologies are used to estimate the ice reserves: three area-volume relations, one slope-dependent volume estimation method, and two ice-thickness distribution models are applied to a recent, detailed, and complete glacier inventory of the HK region, spanning over the period 2000-2010 and revealing an ice coverage of 40 775 km(2). An uncertainty and sensitivity assessment is performed to investigate the influence of the observed glacier area and important model parameters on the resulting total ice volume. Results of the two ice-thickness distribution models are validated with local ice-thickness measurements at six glaciers. The resulting ice volumes for the entire HK region range from 2955 to 4737 km(3), depending on the approach. This range is lower than most previous estimates. Results from the ice thickness distribution models and the slope-dependent thickness estimations agree well with measured local ice thicknesses. However, total volume estimates from area-related relations are larger than those from other approaches. The study provides evidence on the significant effect of the selected method on results and underlines the importance of a careful and critical evaluation.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Role of the Molar Volume on Estimated Diffusion Coefficients

The role of the molar volume on the estimated diffusion parameters has been speculated for decades. The Matano-Boltzmann method was the first to be developed for the estimation of the variation of the interdiffusion coefficients with composition. However, this could be used only when the molar volume varies ideally or remains constant. Although there are no such systems, this method is still being used to consider the ideal variation. More efficient methods were developed by Sauer-Freise, Den Broeder, and Wagner to tackle this problem. However, there is a lack of research indicating the most efficient method. We have shown that Wagner's method is the most suitable one when the molar volume deviates from the ideal value. Similarly, there are two methods for the estimation of the ratio of intrinsic diffusion coefficients at the Kirkendall marker plane proposed by Heumann and van Loo. The Heumann method, like the Matano-Boltzmann method, is suitable to use only when the molar volume varies more or less ideally or remains constant. In most of the real systems, where molar volume deviates from the ideality, it is safe to use the van Loo method. We have shown that the Heumann method introduces large errors even for a very small deviation of the molar volume from the ideal value. On the other hand, the van Loo method is relatively less sensitive to it. Overall, the estimation of the intrinsic diffusion coefficient is more sensitive than the interdiffusion coefficient.

Mechanism of improvement in the strength and volume change behavior of lime stabilized soil

An attempt has been made to bring out the influence on strength and volume change behavior of fabric changes and new cementitious compound formation in a soil upon addition of various lime contents and with curing periods. The effects of changes in fabric of treatment with various lime contents (0, 2,4 and 6%) and with curing periods (0, 7, 14 and 28 days) have been evaluated by one-dimensional consolidation tests, in terms of void ratio changes and compressibility. The strength of soil treated with different lime contents with curing periods up to 28 days, and with the optimum lime content of 6% up to one year has been determined by unconfined compression tests. Comparison of effects of lime on the strength and volume change behavior of the soil brings out that the formation of flocculated fabric and cation exchange significantly reduces the compressibility of soil but marginally increases the strength. Cementation of soil particles and filling with cementitious compounds of the voids of flocculated fabric in the soil marginally reduces the compressibility but significantly increases the strength. Thus, the mechanism of volume change behavior of soil treated with lower lime content at short curing periods is distinctly different from that of the soil treated with optimum lime content at longer curing periods. This is consistent with the increase in the permeability caused by the addition from 2 to 4% lime and the decrease following the addition of 6% lime. Changes consistent with mechanical behavior have been determined by scanning electron microscope, X-ray diffraction and thermal analyses, energy dispersive X-ray spectrometer and pH value in microstructure, mineralogy, chemical composition and alkalinity, respectively. (C) 2015 Published by Elsevier B.V.

Volume change behavior of lime treated gypseous soil - influence of mineralogy and microstructure

Series of oedometer tests and micro-analytical studies (XRD, SEM and EDAX) have been carried out to investigate the influence of varying gypsum content on swell, compressibility and permeability of lime treated montmorillonitic soil after curing for different period. Immediate swell is observed on inundation of compacted samples with water and continuously increased with gypsum content. However, changes in swell are found to be marginal with curing. This is attributed to the formation and growth of ettringite crystals by ionic reactions of aluminum calcium-sulfate in the presence of water which is confirmed through detailed micro-analysis. The higher swell in uncured specimens and gradual reduction in swell with increase in curing periods are due to relative dominance of formation and growth of ettringite and cementitious compounds, respectively. Also, the ionic reaction products are found to bear a significant influence on the compressibility and permeability behavior. (C) 2015 Elsevier B.V. All rights reserved.

Limited-view light sheet fluorescence microscopy for three dimensional volume imaging

We propose and demonstrate a limited-view light sheet microscopy (LV-LSM) for three dimensional (3D) volume imaging. Realizing that longer and frequent image acquisition results in significant photo-bleaching, we have taken limited angular views (18 views) of the macroscopic specimen and integrated with maximum likelihood (ML) technique for reconstructing high quality 3D volume images. Existing variants of light-sheet microscopy require both rotation and translation with a total of approximately 10-fold more views to render a 3D volume image. Comparatively, LV-LSM technique reduces data acquisition time and consequently minimizes light-exposure by many-folds. Since ML is a post-processing technique and highly parallelizable, this does not cost precious imaging time. Results show noise-free and high contrast volume images when compared to the state-of-the-art selective plane illumination microscopy. (C) 2015 AIP Publishing LLC.

Perturbational Finite Volume Method For The Solution of 2-D Navier-Stokes Equations on Unstructured and Structured Colocated Meshes

Based on the first-order upwind and second-order central type of finite volume( UFV and CFV) scheme, upwind and central type of perturbation finite volume ( UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from 0.3 to 0. 8 and convergence perform excellent with Reynolds number variation from 102 to 104.

Simulation of Free-Surface Flow in a Tank Using the Navier-Stokes Model and Unstructured Finite Volume Method

Modelling free-surface flow has very important applications in many engineering areas such as oil transportation and offshore structures. Current research focuses on the modelling of free surface flow in a tank by solving the Navier-Stokes equation. An unstructured finite volume method is used to discretize the governing equations. The free surface is tracked by dynamically adapting the mesh and making it always surface conforming. A mesh-smoothing scheme based on the spring analogy is also implemented to ensure mesh quality throughout the computaiton. Studies are performed on the sloshing response of a liquid in an elastic container subjected to various excitation frequencies. Further investigations are also carried out on the critical frequency that leads to large deformation of the tank walls. Another numerical simulation involves the free-surface flow past as submerged obstacle placed in the tank to show the flow separation and vortices. All these cases demonstrate the capability of this numerical method in modelling complicated practical problems.

Effect of liquid bridge volume on the instability in small-Prandtl-number half zones

A perturbation method is used to examine the linear instability of thermocapillary convection in a liquid bridge of floating half-zone filled with a small Prandtl number fluid. The influence of liquid bridge volume on critical Marangoni number and flow features is analyzed. The neutral modes show that the instability is mainly caused by the bulk flow that is driven by the nonuniform thermocapillary forces acting on the free surface. The hydrodynamic instability is dominant in the case of small Prandtl number fluid and the first instability mode is a stationary bifurcation. The azimuthal wave number for the most dangerous mode depends on the liquid bridge volume, and is not always two as in the case of a cylindrical liquid bridge with aspect ratio near 0.6. Its value may be equal to unity when the liquid bridge is relatively slender.

Constant volume combustion of aluminum and cornstarch dust in microgravity

The subject of the present work is to report an experimental comparative study of the effect of dispersion-induced turbulence on dust combustion in constant volume vessel, carried out both in normal gravity and in microgravity environment. Dispersion system with small scale of turbulence, creating uniform homogeneous mixture, was used in experiments. To improve reproducibility of the explosion data an ignitor of small energy, with local soft ignition was developed. Both factors contributed to acquisition of more reproducible experimental data. In experiments under microgravity conditions a dust suspension during combustion remains constant. This makes possible to study dust explosion under stationary dust suspension without influence of turbulence.

Influence of liquid bridge volume on instability of floating half zone convection

Thermocapillary instabilities on floating half zone convection in microgravity environment were investigated by linear instability analysis method. The critical Marangoni numbers were obtained and compared with the experimental ones. The influences of the liquid bridge volume and the aspect ratio on the critical Marangoni number were analyzed. It is found that the liquid bridge volume and the aspect ratio have great influence on the critical Marangoni number. There was a gap region where the oscillatory convection will not be observed in present analyses and in experiments in the curve of the critical Marangoni number vs the liquid bridge volume for the case of large Prandtl number and small aspect ratio.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection III. Three-dimensional model

An unsteady and three-dimensional model of the floating-half-zone convection on the ground is studied by the direct numerical simulation for the medium of 10 cSt silicon oil, and the influence of the liquid bridge volume on the critical applied temperature difference is especially discussed. The marginal curves for the onset of oscillation are separated into two branches related, respectively, to the slender liquid bridge and the fat liquid bridge. The oscillatory features of the floating-half-zone convection are also discussed.