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.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.

An Optimizing Code Generator for a Class of Lattice-Boltzmann Computations

The Lattice-Boltzmann method (LBM), a promising new particle-based simulation technique for complex and multiscale fluid flows, has seen tremendous adoption in recent years in computational fluid dynamics. Even with a state-of-the-art LBM solver such as Palabos, a user has to still manually write the program using library-supplied primitives. We propose an automated code generator for a class of LBM computations with the objective to achieve high performance on modern architectures. Few studies have looked at time tiling for LBM codes. We exploit a key similarity between stencils and LBM to enable polyhedral optimizations and in turn time tiling for LBM. We also characterize the performance of LBM with the Roofline performance model. Experimental results for standard LBM simulations like Lid Driven Cavity, Flow Past Cylinder, and Poiseuille Flow show that our scheme consistently outperforms Palabos-on average by up to 3x while running on 16 cores of an Intel Xeon (Sandybridge). We also obtain an improvement of 2.47x on the SPEC LBM benchmark.

Void Geometry Driven Spin Crossover in Zeolite-Encapsulated Cobalt Tris(bipyridyl) Complex Ion

The cobalt(II) tris(bipyridyl) complex ion encapsulated in zeolite-Y supercages exhibits a thermally driven interconversion between a low-spin and a high-spin state-a phenomenon not observed for this ion either in solid state or in solution. From a comparative study of the magnetism and optical spectroscopy of the encapsulated and unencapsulated complex ion, supported by molecular modeling, such spin behavior is shown to be intramolecular in origin. In the unencapsulated or free state, the [Co(bipy)(3)](2+) ion exhibits a marked trigonal prismatic distortion, but on encapsulation, the topology of the supercage forces it to adopt a near-octahedral geometry. An analysis using the angular overlap ligand field model with spectroscopically derived parameters shows that the geometry does indeed give rise to a low-spin ground state, and suggests a possible scenario for the spin state interconversion.

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

Turbulent mixed convection in a shallow enclosure with a series of heat generating components

Turbulent mixed convection flow and heat transfer in a shallow enclosure with and without partitions and with a series of block-like heat generating components is studied numerically for a range of Reynolds and Grashof numbers with a time-dependent formulation. The flow and temperature distributions are taken to be two-dimensional. Regions with the same velocity and temperature distributions can be identified assuming repeated placement of the blocks and fluid entry and exit openings at regular distances, neglecting the end wall effects. One half of such module is chosen as the computational domain taking into account the symmetry about the vertical centreline. The mixed convection inlet velocity is treated as the sum of forced and natural convection components, with the individual components delineated based on pressure drop across the enclosure. The Reynolds number is based on forced convection velocity. Turbulence computations are performed using the standard k– model and the Launder–Sharma low-Reynolds number k– model. The results show that higher Reynolds numbers tend to create a recirculation region of increasing strength in the core region and that the effect of buoyancy becomes insignificant beyond a Reynolds number of typically 5×105. The Euler number in turbulent flows is higher by about 30 per cent than that in the laminar regime. The dimensionless inlet velocity in pure natural convection varies as Gr1/3. Results are also presented for a number of quantities of interest such as the flow and temperature distributions, Nusselt number, pressure drop and the maximum dimensionless temperature in the block, along with correlations.

Local structure of Sr2FeMoxW1-xO6 double perovskites across the composition-driven metal to insulator transition

Sr2FeMoO6 oxides exhibit a half-metallic ferromagnetic (HM-FM) ground state and peculiar magnetic and magnetotransport properties, which are interesting for applications in the emerging field of spintronics and attractive for fundamental research in the field of heavily correlated electron systems. Sr2FeWO6 is an insulator with an antiferromagnetic (I-AFM) ground state. The solid solutions Sr2FeMoxW1-xO6 also have peculiar properties-W doping enhances chemical order which allows stabilization of the HM-FM state; as the W content exceeds a certain value a metal to insulator transition (MIT) occurs. The role of W in determining the physical properties of Sr2FeMoxW1-xO6 systems has been a matter of intense investigation. This work deals with the problem of the structural and electronic changes related to the MIT from a local perspective by means of x-ray absorption spectroscopy (XAS). This technique allows one to probe in detail the local structure and electronic modifications around selected absorber ions (W, Mo, Fe and Sr in our case). The results of XAS analysis in the whole composition range (0 <= x <= 1), in the near edge (XANES) and extended (EXAFS) regions, demonstrate an abrupt change of the local structure around the Fe and Mo sites at the critical composition, x(c). This change represents the microstructural counterpart associated with the MIT. Conversely, the local structure and electronic configuration of W ions remain unaltered in the whole composition range, suggesting indirect participation of W in the MIT.

Competing A-site and B-site driven ferroelectric instabilities in the (1-x)PbTiO3-(x)BiAlO3 system

The system (1-x)PbTiO3-(x)BiAlO3 has been investigated with regard to its solid solubility, crystal structure, microstructure, and ferroelectric transition. The unit cell volume and the tetragonality exhibit anomalous behavior near x=0.10. The Curie point (T-C) of PbTiO3 was however found to be nearly unchanged. The study seems to suggest that the decrease in the stability of the ferroelectric state due to dilution of the Ti-sublattice by smaller sized Al+3 ions is compensated by the increase in the ferroelectric stability by the Bi+3 ions.

Stress and strain-driven algorithmic formulations for finite strain viscoplasticity for hybrid and standard finite elements

This work deals with the formulation and implementation of finite deformation viscoplasticity within the framework of stress-based hybrid finite element methods. Hybrid elements, which are based on a two-field variational formulation, are much less susceptible to locking than conventional displacement-based elements. The conventional return-mapping scheme cannot be used in the context of hybrid stress methods since the stress is known, and the strain and the internal plastic variables have to be recovered using this known stress field.We discuss the formulation and implementation of the consistent tangent tensor, and the return-mapping algorithm within the context of the hybrid method. We demonstrate the efficacy of the algorithm on a wide range of problems.

Studies on micro explosive driven blast wave propagation in confined domains using NONEL tubes

Explosive driven micro blast waves are generated in the laboratory using NONEL tubes. The explosive mixture coated to the inner walls of the plastic Nonel tube comprises of HMX and Aluminum ( 18mg/m). The detonation is triggered electrically to generate micro blast waves from the open end of the tube. Flow visualization and over pressure measurements have been carried out to understand the propagation dynamics of these micro-blast waves in both confined and unconfined domains. The classical cubic root law used for large scale blast correlation appears to hold good even for these micro-blasts generated in the laboratory.

Material processing and surface reaction studies in free piston driven shock tube

Recently established moderate size free piston driven hypersonic shock tunnel HST3 along with its calibration is described here. The extreme thermodynamic conditions prevalent behind the reflected shock wave have been utilized to study the catalytic and non-catalytic reactions of shock heated test gases like Ar, N2 or O2 with different material like C60 carbon, zirconia and ceria substituted zirconia. The exposed test samples are investigated using different experimental methods. These studies show the formation of carbon nitride due to the non-catalytic interaction of shock heated nitrogen gas with C60 carbon film. On the other hand, the ZrO2 undergoes only phase transformation from cubic to monoclinic structure and Ce0.5Zr0.5O2 in fluorite cubic phase changes to pyrochlore (Ce2Zr2O7±δ) phase by releasing oxygen from the lattice due to heterogeneous catalytic surface reaction.

Low Complexity LID using Pruned Pattern Tables of LZW

We present two discriminative language modelling techniques for Lempel-Ziv-Welch (LZW) based LID system. The previous approach to LID using LZW algorithm was to directly use the LZW pattern tables forlanguage modelling. But, since the patterns in a language pattern table are shared by other language pattern tables, confusability prevailed in the LID task. For overcoming this, we present two pruning techniques (i) Language Specific (LS-LZW)-in which patterns common to more than one pattern table are removed. (ii) Length-Frequency product based (LF-LZW)-in which patterns having their length-frequency product below a threshold are removed. These approaches reduce the classification score (Compression Ratio [LZW-CR] or the weighted discriminant score [LZW-WDS]) for non native languages and increases the LID performance considerably. Also the memory and computational requirements of these techniques are much less compared to basic LZW techniques.

Large Low-Frequency Fluctuations in the Velocity of a Driven Vortex Lattice in a Single Crystal of 2H-NbSe2 Superconductor

The driven state of a well-ordered flux line lattice in a single crystal of 2H-NbSe2 in the time domain has revealed the presence of substantial fluctuations in velocity, with large and distinct time periods (similar to seconds). A superposition of a periodic drive in the driven vortex lattice causes distinct changes in these fluctuations. We propose that prior to the onset of the peak effect there exists a heretofore unexplored regime of coherent dynamics, with unexpected behavior in velocity fluctuations.

Effect of surface radiation on conjugate natural convection in a horizontal annulus driven by inner heat generating solid cylinder

This paper reports a numerical study of the laminar conjugate natural convection heat transfer with and without the interaction of the surface radiation in a horizontal cylindrical annulus formed between an inner heat generating solid circular cylinder and an outer isothermal circular boundary. Numerical solutions are obtained by solving the governing equations with a pressure correction method on a collocated (non-staggered) mesh. Steady-state results are presented for the flow and temperature distributions and Nusselt numbers for the heat generation based Grashof number ranging from 10(7) to 10(10), solid-to-fluid thermal conductivity ratios of 1, 5, 10, 50 and 100, radius ratios of 0.226 and 0.452 and surface emissivities of 0-0.8 with air as the working medium. It is observed that surface radiation reduces the convective heat transfer in the annulus compared to the pure natural convection case and enhances the overall Nusselt number.