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.

Rotational cut-off and anharmonicity effects on thermodynamic properties of gases at high temperatures

The exact expressions for the partition function (Q) and the coefficient of specific heat at constant volume (Cv) for a rotating-anharmonic oscillator molecule, including coupling and rotational cut-off, have been formulated and values of Q and Cv have been computed in the temperature range of 100 to 100,000 K for O2, N2 and H2 gases. The exact Q and Cv values are also compared with the corresponding rigid-rotator harmonic-oscillator (infinite rotational and vibrational levels) and rigid-rotator anharmonic-oscillator (infinite rotational levels) values. The rigid-rotator harmonic-oscillator approximation can be accepted for temperatures up to about 5000 K for O2 and N2. Beyond these temperatures the error in Cv will be significant, because of anharmonicity and rotational cut-off effects. For H2, the rigid-rotator harmonic-oscillator approximation becomes unacceptable even for temperatures as low as 2000 K.

Structure and optical properties of Cd-substituted ZnO (Zn1-xCdxO) nanostructures synthesized by the high-pressure solution route

We report the synthesis of Cd-substituted ZnO nanostructures (Zn1-xCdxO with x up to approximate to 0.09) by the high-pressure solution growth method. The synthesized nanostructures comprise nanocrystals that are both particles (similar to 10-15 nm) and rods which grow along the [002] direction as established by transmission electron microscope (TEM) and x-ray diffraction (XRD) analysis. Rietveld analysis of the XRD data shows a monotonic increase of the unit cell volume with the increase of Cd concentration. The optical absorption, as well as the photoluminescence (PL), shows a red shift on Cd substitution. The line width of the PL spectrum is related to the strain inhomogeneity and it peaks in the region where the CdO phase separates from the Zn1-xCdxO nanostructures. The time-resolved photoemission showed a long-lived (similar to 10 ns) component. We propose that the PL behaviour of the Zn1-xCdxO is dominated by strain in the sample with the red shift of the PL linked to the expansion of the unit cell volume on Cd substitution.

Structure and optical properties of Cd-substituted ZnO (Zn1-xCdxO) nanostructures synthesized by the high-pressure solution route

We report the synthesis of Cd-substituted ZnO nanostructures (Zn1-xCdxO with x up to approximate to 0.09) by the high-pressure solution growth method. The synthesized nanostructures comprise nanocrystals that are both particles (similar to 10-15 nm) and rods which grow along the [002] direction as established by transmission electron microscope (TEM) and x-ray diffraction (XRD) analysis. Rietveld analysis of the XRD data shows a monotonic increase of the unit cell volume with the increase of Cd concentration. The optical absorption, as well as the photoluminescence (PL), shows a red shift on Cd substitution. The line width of the PL spectrum is related to the strain inhomogeneity and it peaks in the region where the CdO phase separates from the Zn1-xCdxO nanostructures. The time-resolved photoemission showed a long-lived (similar to 10 ns) component. We propose that the PL behaviour of the Zn1-xCdxO is dominated by strain in the sample with the red shift of the PL linked to the expansion of the unit cell volume on Cd substitution.

Parallelization of an unstructured data based cell centre finite volume code, HIFUN-3D

This work describes the parallelization of High Resolution flow solver on unstructured meshes, HIFUN-3D, an unstructured data based finite volume solver for 3-D Euler equations. For mesh partitioning, we use METIS, a software based on multilevel graph partitioning. The unstructured graph used for partitioning is associated with weights both on its vertices and edges. The data residing on every processor is split into four layers. Such a novel procedure of handling data helps in maintaining the effectiveness of the serial code. The communication of data across the processors is achieved by explicit message passing using the standard blocking mode feature of Message Passing Interface (MPI). The parallel code is tested on PACE++128 available in CFD Center

High-pressure NMR investigations of the protonic conductor (NH4)4Fe(CN)6.1.5H2O

The effect of high hydrostatic pressure up to 1.5 GPa on ionic motion in (NH4)4Fe(CN)6.1.5H2O has been studied by wide-line 1H NMR experiments performed in the temperature range from room temperature to 77 K. The experiments at room temperature have shown a large increase in the second moment at 0.45 GPa as a result of a pressure-induced phase transition. The temperature dependence study up to 0.425 GPa has shown a gradual increase in the values of activation energy and attempt frequency with increase in pressure. The activation volume for motion at 300 K has been estimated to be 6% of molar volume. Vacancy-assisted ionic jumps are concluded to be the mode of charge transport. Second moments estimated at 77 K show evidence for tunnelling reorientation of at least one of the two NH4+ groups in the compound.

Reactive hot pressing of ZrB2-ZrCx ultra-high temperature ceramic composites with the addition of SiC particulate

Starting with non-stoichiometric Zr-B4C powder mixture ZrB2-ZrC matrix composites with SiC particulate addition have been made. It was found that variable amounts (5-25 vol%) of SiC could be incorporated and reactively hot pressed (RHPed) to relative densities of 97-99% at 1400-1500 degrees C. This technique has the potential to fabricate ZrB2-based matrices at low temperatures with a variety of reinforcements whose composition and volume fraction are not limited by stoichiometric considerations. The hardness of the composites is in the range of 17-22 GPa. (C) 2010 Elsevier Ltd. All rights reserved.

A numerical investigation of ductile fracture initiation in a high-strength low-alloy steel

In this work, static and drop-weight impact experiments, which have been conducted using three-point bend fracture specimens of a high-strength low-alloy steel, are analysed by performing finite-element simulations. The Gurson constitutive model that accounts for the ductile failure mechanisms of microvoid nucleation, growth and is employed within the framework of a finite deformation plasticity theory. Two populations of second-phase particles are considered, including large inclusions which initiate voids at an early stage and small particles which require large strains to nucleate voids. The most important objective of the work is to assess quantitatively the effects of material inertia, strain rate sensitivity and local adiabatic temperature rise (due to conversion of plastic work into heat) on dynamic ductile crack initiation. This is accomplished by comparing the evolution histories of void volume fraction near the notch tip in the static analysis with the dynamic analyses. The results indicate that increased strain hardening caused by strain rate sensitivity, which becomes important under dynamic loading, plays a benign role in considerably slowing down the void growth rate near the notch tip. This is partially opposed by thermal softening caused by adiabatic heating near the notch tip.

High swirl-inducing piston bowls in small diesel engines for emission reduction

Detailed three-dimensional CFD simulations involving flow and combustion chemistry are used to study the effect of swirl induced by re-entrant piston bowl geometries on pollutant emissions from a single-cylinder diesel engine. The baseline engine configuration consists of a hemispherical piston bowl and an injector with finite sac volume. The first iteration involved using a torroidal, slightly re-entrant bowl geometry, and a sac-less injector. Pollutant emission measurements indicated a reduction in emissions with this modification. Simulations on both configurations were then conducted to understand the effect of the changes. The simulation results indicate that the selected piston bowl geometry could actually be reducing the in-cylinder swirl and turbulence and the emission reduction may be entirely due to the introduction of the sac-less injector. In-cylinder air motion was then studied in a number of combustion chamber geometries, and a geometry which produced the highest in-cylinder swirl and Turbulence Kinetic Energy (TKE) around the compression top dead centre (TDC) was identified. The optimal nature of this re-entrant piston bowl geometry is confirmed by detailed combustion simulations and emission predictions. (C) 2010 Elsevier Ltd. All rights reserved.

Physics based basis function for vibration analysis of high speed rotating beams

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

Comparative study on flexural response of full and partial depth fiber-reinforced high-strength concrete

This study presents the results of an experimental and analytical comparison of the flexural behavior of a high-strength concrete specimen (no conventional reinforcement) with an average plain concrete cube strength of nearly 65 MPa and containing trough shape steel fibers. Trough shape steel fibers with a volume fraction ranging from 0 to 1.5% and having a constant aspect ratio of 80 have been used in this study. Increased toughness and a more ductile stress-strain response were observed with an increase in fiber content, when the fibers were distributed over the full/partial depth of the beam cross section. Based on the tests, a robust analytical procedure has been proposed to establish the required partial depth to contain fiber-reinforced concrete (FRC) so as to obtain the flexural capacity of a member with FRC over the full depth. It is expected that this procedure will help designers in properly estimating the required partial depth of fibers in composite sections for specific structural applications. Empirical and mechanistic relations have also been proposed in this study to establish the load-deflection behavior of high-strength FRC.

High-pressure kinetics of oxidative copolymerization of styrene with alpha-methyl styrene

This is the first report on the study carried out on high-pressure free-radical initiated oxidated copolymerization of styrene (STY) with alpha-methylstyrene (AMS) at various temperatures (45-65degreesC) at constant pressure (100 psi) and then at various pressures (50-300 psi) keeping the temperature (50degreesC) constant. The compositions of the copolyperoxides obtained from the H-1 NMT spectra were utilized to determine the reactivity ratios of the monomers. The reactivity ratios indicate that STY forms an ideal copolyperoxide with AMS and the copolyperoxide is richer in AMS. The effect of temperature and oxygen pressure in the reactivity ratios of the monomers was studied. The rates of copolymerization (R-p) were used to determine the overall activation energies (E-a) and activation volume (DeltaV(#)) of copolymerization. The unusually higher values of the DeltaV(#) may be due to the pressurizing fluid oxygen which itself is a reactant in the copolymerization, the side reactions, and the chain-transfer reactions occuring during copolymerizations.

Influence of martensite content and morphology on the toughness and fatigue behavior of high-martensite dual-phase steels

A series of high-martensite dual-phase (HMDP) steels exhibiting a 0.3 to 0.8 volume fraction of martensite (V m ), produced by intermediate quenching (IQ) of a vanadium and boron-containing microalloyed steel, have been studied for toughness and fatigue behavior to supplement the contents of a recent report by the present authors on the unusual tensile behavior of these steels. The studies included assessment of the quasi-static and dynamic fracture toughness and fatigue-crack growth (FCG) behavior of the developed steels. The experimental results show that the quasi-static fracturetoughness (K ICV ) increases with increasing V m in the range between V m =0.3 and 0.6 and then decreases, whereas the dynamic fracture-toughness parameters (K ID , K D , and J ID ) exhibit a significant increase in their magnitudes for steels containing 0.45 to 0.60 V m before achieving a saturation plateau. Both the quasi-static and dynamic fracture-toughness values exhibit the best range of toughnesses for specimens containing approximately equal amounts of precipitate-free ferrite and martensite in a refined microstructural state. The magnitudes of the fatigue threshold in HMDP steels, for V m between 0.55 and 0.60, appear to be superior to those of structural steels of a similar strength level. The Paris-law exponents (m) for the developed HMDP steels increase with increasing V m , with an attendant decrease in the pre-exponential factor (C).

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 2. The fluctuating-force model

A fluctuating-force model is developed for representing the effect of the turbulent fluid velocity fluctuations on the particle phase in a turbulent gas–solid suspension in the limit of high Stokes number, where the particle relaxation time is large compared with the correlation time for the fluid velocity fluctuations. In the model, a fluctuating force is incorporated in the equation of motion for the particles, and the force distribution is assumed to be an anisotropic Gaussian white noise. It is shown that this is equivalent to incorporating a diffusion term in the Boltzmann equation for the particle velocity distribution functions. The variance of the force distribution, or equivalently the diffusion coefficient in the Boltzmann equation, is related to the time correlation functions for the fluid velocity fluctuations. The fluctuating-force model is applied to the specific case of a Couette flow of a turbulent particle–gas suspension, for which both the fluid and particle velocity distributions were evaluated using direct numerical simulations by Goswami & Kumaran (2010). It is found that the fluctuating-force simulation is able to quantitatively predict the concentration, mean velocity profiles and the mean square velocities, both at relatively low volume fractions, where the viscous relaxation time is small compared with the time between collisions, and at higher volume fractions, where the time between collisions is small compared with the viscous relaxation time. The simulations are also able to predict the velocity distributions in the centre of the Couette, even in cases in which the velocity distribution is very different from a Gaussian distribution.

Stream-Close: Fast mining of Closed Frequent Itemsets in high speed data streams

With the emergence of large-volume and high-speed streaming data, the recent techniques for stream mining of CFIpsilas (closed frequent itemsets) will become inefficient. When concept drift occurs at a slow rate in high speed data streams, the rate of change of information across different sliding windows will be negligible. So, the user wonpsilat be devoid of change in information if we slide window by multiple transactions at a time. Therefore, we propose a novel approach for mining CFIpsilas cumulatively by making sliding width(ges1) over high speed data streams. However, it is nontrivial to mine CFIpsilas cumulatively over stream, because such growth may lead to the generation of exponential number of candidates for closure checking. In this study, we develop an efficient algorithm, stream-close, for mining CFIpsilas over stream by exploring some interesting properties. Our performance study reveals that stream-close achieves good scalability and has promising results.