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.

Size effect on the lattice parameter of KCl during mechanical milling

The size effect on the lattice parameter of ionic KCl nanocrystals was studied systematically during mechanical milling of Pure KCl powder under vacuum. The results suggest anomalous lattice expansion, with the lattice parameter increasing from 6.278 angstrom at d = 6 mu m to 6.30307 angstrom at d = 85 mn. The defects generated during ball milling of KCl and surface stress are deemed to be responsible for this lattice parameter expansion. (C) 2009 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Simulation of laminar flow in a three-dimensional lid-driven cavity by lattice Boltzmann method

Purpose - The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid-driven flow in a three-dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two-dimensional (2D) cavity. Design/methodology/approach - The second-order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25-4) and Reynolds numbers (0.01-1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier-Stokes (NS) equations and with available experimental data. Findings - The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side-walls in transverse direction (in 3D), due to which the vorticity in the core-region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry-plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid-motion are observed; these are weak at low. Reynolds numbers, but become quite strong at high Reynolds numbers. Originality/value - The findings will be useful in the study of variety of enclosed fluid flows.

Modelling heterogeneity of concrete using 2D lattice network for concrete fracture and comparison with AE study

In this paper, numerical modelling of fracture in concrete using two-dimensional lattice model is presented and also a few issues related to lattice modelling technique applicable to concrete fracture are reviewed. A comparison is made with acoustic emission (AE) events with the number of fractured elements. To implement the heterogeneity of the plain concrete, two methods namely, by generating grain structure of the concrete using Fuller's distribution and the concrete material properties are randomly distributed following Gaussian distribution are used. In the first method, the modelling of the concrete at meso level is carried out following the existing methods available in literature. The shape of the aggregates present in the concrete are assumed as perfect spheres and shape of the same in two-dimensional lattice network is circular. A three-point bend (TPB) specimen is tested in the experiment under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/sec and the fracture process in the same TPB specimen is modelled using regular triangular 2D lattice network. Load versus crack mouth opening isplacement (CMOD) plots thus obtained by using both the methods are compared with experimental results. It was observed that the number of fractured elements increases near the peak load and beyond the peak load. That is once the crack starts to propagate. AE hits also increase rapidly beyond the peak load. It is compulsory here to mention that although the lattice modelling of concrete fracture used in this present study is very similar to those already available in literature, the present work brings out certain finer details which are not available explicitly in the earlier works.

Verification of the applicability of lattice model to concrete fracture by AE study

Notched three-point bend specimens (TPB) were tested under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/s and the entire fracture process was simulated using a regular triangular two-dimensional lattice network only over the expected fracture proces zone width. The rest of the beam specimen was discretised by a coarse triangular finite element mesh. The discrete grain structure of the concrete was generated assuming the grains to be spherical. The load versus CMOD plots thus simulated agreed reasonably well with the experimental results. Moreover, acoustic emission (AE) hits were recorded during the test and compared with the number of fractured lattice elements. It was found that the cumulative AE hits correlated well with the cumulative fractured lattice elements at all load levels thus providing a useful means for predicting when the micro-cracks form during the fracturing process, both in the pre-peak and in the post-peak regimes.

Engineering behavior of uplifted Smectite-rich Cochin and Mangalore marine clays

The present study aims to assess whether the smectite-rich Cochin and Mangalore clays, which were deposited in a marine medium and subsequently uplifted, exhibit consistency limits response typical of expanding lattice or nonexpanding (fixed) lattice-type clays on artificially changing the chemical environment. The chemical and engineering behaviors of Cochin and Mangalore marine clays are also compared with those of the smectite-rich Ariake Bay marine clay from Japan. Although Cochin, Mangalore, and Ariake clays contain comparable amounts of smectite (32-45%), Ariake clay exhibits lower consistency limits and much higher ranges of liquidity indices than the Indian marine clays. The lower consistency limits of the Ariake clay are attributed to the absence of well-developed, long-range, interparticle forces associated with the clay. Also, Ariake clay exhibits a significantly large (48-714 times) decrease in undrained strength on remolding in comparison to Cochin and Mangalore clays (sensitivity ranges between 1 and 4). A preponderance of long-range, interparticle forces reflected in the high consistency limits of Cochin and Mangalore clays (wL range from 75 to 180%) combined with low natural water contents yield low liquidity indices (typically <1) and high, remolded, undrained strengths and are considered to be responsible for the low sensitivity of the Indian marine clays.

Role of lattice defects and crystallite morphology in the UV and visible-light-induced photo-catalytic properties of combustion-prepared TiO2

The physico-chemical, photo-physical and micro-structural properties responsible for the strikingly different photocatalytic behavior of combustion-prepared TiO2 (c.TiO2) and Degussa P25 (d.TiO2) samples are elucidated in this study. Electron microscopy and selected area electron diffraction micrographs revealed that the two samples exhibited different morphologies. The grains of c.TiO2 were spherical and comprised of 5-6 nm size primary particle. On the other hand, d.TiO2 consisted of large (0.5-3.0 mu m) size and irregular shape aggregates having primary particles of 15-40 nm cross-sectional diameter. The ESR study revealed that the presence of certain defect states in c.TiO2 helped in stabilization of O-. and Ti3+-OH type species during room-temperature UV-irradiation. No such paramagnetic species were however formed over d.TiO2 under similar conditions. C1s and Ti 2p XPS spectra provide evidence for the presence of some lattice vacancies in c.TiO2 and also for the bulk Ti4+ -> Ti3+ conversion during its UV-irradiation. Compared to d.TiO2, c.TiO2 displayed considerably higher activity for discoloration of methyl orange but very poor activity for splitting of water, both under UV and visible light radiations. This is attributed to enhanced surface adsorption of dye molecules over c.TiO2, because of its textural features and also the presence of photo-active ion-radicals. On the other hand, the poor activity of c.TiO2 for water splitting is related to certain defect-induced inter-band charge trapping states in the close vicinity of valence and conduction bands of c.TiO2, as revealed by thermoluminescence spectroscopy. Further, the dispersion of nanosize gold particles gave rise to augmented activity of both the catalysts, particularly for water splitting. This is explained by the promotional role of Au-0 or Au-0/TiO2 interfacial sites in the adsorption and charge-adsorbate interaction processes. (C) 2011 Elsevier B.V. All rights reserved.

Role of organic fluorine in crystal engineering

The design of compounds with novel and improved physico-chemical properties as advanced functional materials with a specific application spectrum requires the knowledge about possible supramolecular packing motifs and their experimental control in crystalline lattice. Besides the structure of the individual molecule, non-covalent interactions play a significant role in the determination of molecular conformation, along with the formation of three-dimensional supramolecular architecture in a crystal as a requirement for molecular recognition processes, and the related bioactivity. Involvement of functional groups will contribute to the formation of a predefined packing motif due to their well-defined interactions. The strength and directionality of these interactions create characteristic packing motifs, which can be used for the design of supramolecular arrangements by the development of appropriate strategies for the precise control of their topology. Most relevant of these non-covalent interactions are stacking interactions and hydrogen bonds, which have been subjects of extensive study in the last two decades. In recent literature, substantial efforts have been put in by various researchers towards the understanding of interactions involving organic fluorine and the role they play in generating different packing motifs which guides assembling of molecules in the crystal lattice.

Neutron diffraction study of the coupling between spin, lattice, and structural degrees of freedom in 0.8BiFeO(3)-0.2PbTiO(3)

A powder neutron diffraction study was carried out on 0.8BiFeO(3)-0.2PbTiO(3) in the temperature range 27-1000 degrees C. The system exhibits magnetic transition at similar to 300 degrees C and a rhombohedral (R3c)-cubic (Pm3m) ferroelectric phase transition at similar to 650 degrees C. Anomalous variation in the lattice parameters and the octahedral tilt angle is observed across the magnetic transition temperature. In the magnetic phase, the c parameter is contracted and the octahedral tilt angle is slightly increased. The results suggest coupling between the spin, lattice and structural degrees of freedom. (C) 2011 American Institute of Physics. doi:10.1063/1.3555093]

Mean field lattice model for adsorption isotherms in zeolite NaA

Using a lattice model for adsorption in microporous materials, pure component adsorption isotherms are obtained within a mean field approximation for methane at 300 K and xenon at 300 and 360 K in zeolite NaA. It is argued that the increased repulsive adsorbate-adsorbate interactions at high coverages must play an important role in determining the adsorption behavior. Therefore, this feature is incorporated through a "coverage-dependent interaction'' model, which introduces a free, adjustable parameter. Another important feature, the site volume reduction, has been treated in two ways: a van der Waal model and a 1D hard-rod theory [van Tassel et al., AIChE J. 40, 925 (1994)]; we have also generalized the latter to include all possible adsorbate overlap scenarios. In particular, the 1D hard-rod model, with our coverage-dependent interaction model, is shown to be in best quantitative agreement with the previous grand canonical Monte Carlo isotherms. The expressions for the isosteric heats of adsorption indicate that attractive and repulsive adsorbate-adsorbate interactions increase and decrease the heats of adsorption, respectively. It is concluded that within the mean field approximation, our simple model for repulsive interactions and the 1D hard-rod model for site volume reduction are able to capture most of the important features of adsorption in confined regions. (C) 1999 American Institute of Physics. [S0021-9606(99)70515-5].

Statistical thermodynamics of lattice models in zeolites: Implications of local versus global mean field interactions

The statistical thermodynamics of adsorption in caged zeolites is developed by treating the zeolite as an ensemble of M identical cages or subsystems. Within each cage adsorption is assumed to occur onto a lattice of n identical sites. Expressions for the average occupancy per cage are obtained by minimizing the Helmholtz free energy in the canonical ensemble subject to the constraints of constant M and constant number of adsorbates N. Adsorbate-adsorbate interactions in the Brag-Williams or mean field approximation are treated in two ways. The local mean field approximation (LMFA) is based on the local cage occupancy and the global mean field approximation (GMFA) is based on the average coverage of the ensemble. The GMFA is shown to be equivalent in formulation to treating the zeolite as a collection of interacting single site subsystems. In contrast, the treatment in the LMFA retains the description of the zeolite as an ensemble of identical cages, whose thermodynamic properties are conveniently derived in the grand canonical ensemble. For a z coordinated lattice within the zeolite cage, with epsilon(aa) as the adsorbate-adsorbate interaction parameter, the comparisons for different values of epsilon(aa)(*)=epsilon(aa)z/2kT, and number of sites per cage, n, illustrate that for -1 0. We compare the isotherms predicted with the LMFA with previous GMFA predictions [K. G. Ayappa, C. R. Kamala, and T. A. Abinandanan, J. Chem. Phys. 110, 8714 (1999)] (which incorporates both the site volume reduction and a coverage-dependent epsilon(aa)) for xenon and methane in zeolite NaA. In all cases the predicted isotherms are very similar, with the exception of a small steplike feature present in the LMFA for xenon at higher coverages. (C) 1999 American Institute of Physics. [S0021-9606(99)70333-8].

Studies in crystal engineering: Steering ability of fluorine in styrylcoumarins

In continuation of our studies on crystal engineering using fluorine as a steering group, the photobehaviour of di and tri fluoro 4-styrylcoumarins has been examined. It is found that out of the five derivatives, four crystallize into P-packing mode producing syn-HH photodimer upon irradiation whereas the parent hydrocarbon produces an anti K-T dimer. The packing features of the photolabile crystals of 4-(4-fluorostyryl)-6-fluorocoumarin (1), 4-(2,6-difluorostyryl) 6-fluorocoumarin (2) and the photodimer (3a) of 4-(2,6-fluorostyryl)-7-fluorocoumarin (3) have been determined by single crystal X-ray diffraction studies. The stereochemistry of the photodimer of 4-(2-fluorostyryl)-6-fluorocoumarin (4) is deduced based on preliminary X-ray crystallographic data. However, 4-(2,6-difluorostyryl) coumarin (5) is photoinert. The remarkable steering ability of fluorine is established with the molecular packing in the crystal lattice leading to the formation of syn H-H dimer in the above four examples. (C) 1999 Elsevier Science Ltd. All rights reserved.

Effect of lattice orientation on crack tip constraint in ductile single crystals

In this work, the effect of lattice orientation on the fields prevailing near a notch tip is investigated pertaining to various constraint levels in FCC single crystals. A modified boundary layer formulation is employed and numerical solutions under mode I, plane strain conditions are generated by assuming an elastic-perfectly plastic FCC single crystal. The analysis is carried out corresponding to different lattice orientations with respect to the notch line. It is found that the near-tip deformation field, especially the development of kink or slip shear bands is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the notch tip are also strongly influenced by the level of T-stress. The present results clearly establish that ductile single crystal fracture geometries would progressively lose crack tip constraint as the T-stress becomes more negative irrespective of lattice orientation. Also, the near-tip field for a range of constraint levels can be characterized by two-parameters such as K-T or J-Q as in isotropic plastic solids.

Theory of Gas Hydrates: Effect of the Approximation of Rigid Water Lattice

One of the assumptions of the van der Waals and Platteeuw theory for gas hydrates is that the host water lattice is rigid and not distorted by the presence of guest molecules. In this work, we study the effect of this approximation on the triple-point lines of the gas hydrates. We calculate the triple-point lines of methane and ethane hydrates via Monte Carlo molecular simulations and compare the simulation results with the predictions of van der Waals and Platteeuw theory. Our study shows that even if the exact intermolecular potential between the guest molecules and water is known, the dissociation temperatures predicted by the theory are significantly higher. This has serious implications to the modeling of gas hydrate thermodynamics, and in spite of the several impressive efforts made toward obtaining an accurate description of intermolecular interactions in gas hydrates, the theory will suffer from the problem of robustness if the issue of movement of water molecules is not adequately addressed.