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.

CFD analysis of high frequency miniature pulse tube refrigerators for space applications with thermal non-equilibrium model

High frequency, miniature, pulse tube cryocoolers are extensively used in space applications because of their simplicity. Parametric studies of inertance type pulse tube cooler are performed with different length-to-diameter ratios of the pulse tube with the help of the FLUENT (R) package. The local thermal non-equilibrium of the gas and the matrix is taken into account for the modeling of porous zones, in addition to the wall thickness of the components. Dynamic characteristics and the actual mechanism of energy transfer in pulse are examined with the help of the pulse tube wall time constant. The heat interaction between pulse tube wall and the oscillating gas, leading to surface heat pumping, is quantified. The axial heat conduction is found to reduce the performance of the pulse tube refrigerator. The thermal non-equilibrium predicts a higher cold heat exchanger temperature compared to thermal equilibrium. The pressure drop through the porous medium has a strong non-linear effect due to the dominating influence of Forchheimer term over that of the linear Darcy term at high operating frequencies. The phase angle relationships among the pressure, temperature and the mass flow rate in the porous zones are also important in determining the performance of pulse tuberefrigerator.

Analyzing information flow in brain networks with nonparametric Granger causality

Multielectrode neurophysiological recording and high-resolution neuroimaging generate multivariate data that are the basis for understanding the patterns of neural interactions. How to extract directions of information flow in brain networks from these data remains a key challenge. Research over the last few years has identified Granger causality as a statistically principled technique to furnish this capability. The estimation of Granger causality currently requires autoregressive modeling of neural data. Here, we propose a nonparametric approach based on widely used Fourier and wavelet transforms to estimate both pairwise and conditional measures of Granger causality, eliminating the need of explicit autoregressive data modeling. We demonstrate the effectiveness of this approach by applying it to synthetic data generated by network models with known connectivity and to local field potentials recorded from monkeys performing a sensorimotor task.

Influence of spatial variability of permeability property on steady state seepage flow and slope stability analysis

In recent years, spatial variability modeling of soil parameters using random field theory has gained distinct importance in geotechnical analysis. In the present Study, commercially available finite difference numerical code FLAC 5.0 is used for modeling the permeability parameter as spatially correlated log-normally distributed random variable and its influence on the steady state seepage flow and on the slope stability analysis are studied. Considering the case of a 5.0 m high cohesive-frictional soil slope of 30 degrees, a range of coefficients of variation (CoV%) from 60 to 90% in the permeability Values, and taking different values of correlation distance in the range of 0.5-15 m, parametric studies, using Monte Carlo simulations, are performed to study the following three aspects, i.e., (i) effect ostochastic soil permeability on the statistics of seepage flow in comparison to the analytic (Dupuit's) solution available for the uniformly constant permeability property; (ii) strain and deformation pattern, and (iii) stability of the given slope assessed in terms of factor of safety (FS). The results obtained in this study are useful to understand the role of permeability variations in slope stability analysis under different slope conditions and material properties. (C) 2009 Elsevier B.V. All rights reserved.

Constitutive equation for elevated temperature flow behavior of plain carbon steels using dimensional analysis

Dimensional analysis using π-theorem is applied to the variables associated with plastic deformation. The dimensionless groups thus obtained are then related and rewritten to obtain the constitutive equation. The constants in the constitutive equation are obtained using published flow stress data for carbon steels. The validity of the constitutive equation is tested for steels with up to 1.54 wt%C at temperatures: 850–1200 °C and strain rates: 6 × 10−6–2 × 10−2 s−1. The calculated flow stress agrees favorably with experimental data.

Response time measurement in flow induced signal generation on semiconductors

Measurable electrical signal is generated when a gas flows over a variety of solids, including doped semiconductors, even at the modest speed of a few meters per second. The underlying mechanism is an interesting interplay of Bernoulli's principle and the Seebeck effect. The electrical signal depends on the square of Mach number (M) and is proportional to the Seebeck coefficient (S) of the solids. Here we present experimental estimate of the response time of the signal rise and fall process, i.e. how fast the semiconductor materials respond to a steady flow as soon as it is set on or off. A theoretical model is also presented to understand the process and the dependence of the response time on the nature and physical dimensions of the semiconductor material used and they are compared with the experimental observations. (c) 2007 Elsevier B.V. All rights reserved.

Correlation between terrestrial heat flow and thermal conductivity

The significant correlation coefficient between the terrestial heat flow and thermal conductivity computed from the continental heat flow data by Horai and Nur [1]2) may be explained as a natural consequence of terrestrial heat flow through a random medium. The theory predicts a value of 0.40 for the correlation coefficient. A simple statistical test shows that the majority of the computed coefficients belong to the statistical population whose mean is equal to the theoretical correlation coefficient. There are, however, a few observations of unsually high correlation coefficient which cannot be explained by the above hypothesis.

Deformation and crystallization of Zr-based amorphous alloys in homogeneous flow regime

The purpose of this study is to experimentally investigate the interaction of inelastic deformation and microstructural changes of two Zr-based bulk metallic glasses (BMGs): Zr41.25Ti13.75Cu12.5Ni10Be22.5 (commercially designated as Vitreloy 1 or Vit1) and Zr46.75Ti8.25Cu7.5Ni10Be27.5 (Vitreloy 4, Vit4). High-temperature uniaxial compression tests were performed on the two Zr alloys at various strain rates, followed by structural characterization using differential scanning calorimetry (DSC) and transmission electron microscopy (TEM). Two distinct modes of mechanically induced atomic disordering in the two alloys were observed, with Vit1 featuring clear phase separation and crystallization after deformation as observed with TEM, while Vit4 showing only structural relaxation with no crystallization. The influence of the structural changes on the mechanical behaviors of the two materials was further investigated by jump-in-strain-rate tests, and flow softening was observed in Vit4. A free volume theory was applied to explain the deformation behaviors, and the activation volumes were calculated for both alloys.

Polarization switching in near-stoichiometric Zn:LiNbO3 at high temperatures

The effect of temperature and stoichiometry on the polarization switching rate in lithium niobate is presented. An increased polarization switching rate in congruent and near-stoichiometric lithium niobate (CLN and SLN) and SLN doped with 1.6 mol% Zn (SLN:Zn(1.6)) is observed using a pulsed field switching technique near the transition temperature (TO. Compared to CLN, the observed switching rate and domain wall mobility for SLN and SLN:Zn(1.6) are higher. The extra charge flow was observed during switching at high temperatures,and is attributed to the creation of defect dipoles and increase in ionic conductivity. Forward domain motion is expected to be the mechanism involved in switching. (C) 2010 Elsevier Ltd. All rights reserved.

Transient Momentum and Enthalpy Transfer in Packed Beds at High Reynolds Numbers

An experimental study for transient temperature response and pressure drop in a randomly packed bed at high Reynolds numbers is presented.The packed bed is used as a compact heat exchanger along with a solid-propellant gas generator, to generate room-temperature gases for use in control actuation, air bottle pressurization, etc. Packed beds of lengths 200 and 300 mm were characterized for packing-sphere-based Reynolds numbers ranging from 0.8 x 10(4) to 8.5 x 10(4).The solid packing used in the bed consisted of phi 9.5 mm steel spheres. The bed-to-particle diameter ratio was with the average packed-bed porosity around 0.43. The inlet flow temperature was unsteady and a mesh of spheres was used at either end to eliminate flow entrance and exit effects. Gas temperature and pressure were measured at the entry, exit,and at three axial locations along centerline in the packed beds. The solid packing temperature was measured at three axial locations in the packed bed. A correlation based on the ratio of pressure drop and inlet-flow momentum (Euler number) exhibited an asymptotically decreasing trend with increasing Reynolds number. Axial conduction across the packed bed was found to he negligible in the investigated Reynolds number range. The enthalpy absorption rate to solid packing from hot gases is plotted as a function of a nondimensional time constant for different Reynolds numbers. A longer packed bed had high enthalpy absorption rate at Reynolds number similar to 10(4), which decreased at Reynolds number similar to 10(5). The enthalpy absorption plots can be used for estimating enthalpy drop across packed bed with different material, but for a geometrically similar packing.

The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.

The high speed water tunnel facility at the Indian Institute of Science

The present article about the high speed water tunnel facility at the Indian Institute of Science, Bangalore, provides a general description of the tunnel circuit, and brief reports on the performance of the facility and some typical results from investigations carried out in it. A unique aspect of the facility is that it has a horizontal resorber in the form of a large cylindrical tank located in the lower leg of the circuit. The facility has been used, among other things, for flow visualization studies, and investigations on marine propeller hydrodynamics and “synthetic cavitation”. The last topic has been primarily developed at the Indian Institute of Science and shows considerable promise for basic work in cavitation inception and noise.

Accelerometer-Based Force Balance for High Enthalpy Facilities

A single component accelerometer-based force balance is developed, calibrated, and used for high enthalpy applications. Functionality of this force balance, for such applications, is demonstrated for the first time during high enthalpy tests in a newly established free piston-driven shock tunnel, HST3, using a 60 degrees apex angle blunt cone model at 0 degrees angle of incidence. Usefulness of this force balance is also confirmed, for much complicated high enthalpy flow situations, during the drag reduction studies with counterflow supersonic jet from the stagnation point.

High-Temperature (1023 K to 1273 K 750 degrees C to 1000 degrees C]) Plastic Deformation Behavior of B-Modified Ti-6Al-4V Alloys: Temperature and Strain Rate Effects

A minor addition of B to the Ti-6Al-4V alloy, by similar to 0.1 wt pct, reduces its as-cast prior beta grain size by an order of magnitude, whereas higher B content leads to the presence of in situ formed TiB needles in significant amounts. An experimental investigation into the role played by these microstructural modifications on the high-temperature deformation behavior of Ti-6Al-4V-xB alloys, with x varying between 0 wt pct and 0.55 wt pct, was conducted. Uniaxial compression tests were performed in the temperature range of 1023 K to 1273 K (750 degrees C to 1000 degrees C) and in the strain rate range of 10(-3) to 10(+1) s(-1). True stress-true strain responses of all alloys exhibit flow softening at lower strain rates and oscillations at higher strain rates. The flow softening is aided by the occurrence of dynamic recrystallization through lath globularization in high temperature (1173 K to 1273 K 900 degrees C to 1000 degrees C]) and a lower strain rate (10(-2) to 10(-3) s(-1)) regime. The grain size refinement with the B addition to Ti64, despite being marked, had no significant effect on this. Oscillations in the flow curve at a higher strain rate (10(0) to 10(+1) s(-1)), however, are associated with microstructural instabilities such as bending of laths, breaking of lath boundaries, generation of cavities, and breakage of TiB needles. The presence of TiB needles affected the instability regime. Microstructural evidence suggests that the matrix cavitation is aided by the easy fracture of TiB needles.

Structural effects on the vaporization of high molecular weight esters

To understand the effect of molecular weight and branching on the heats of vaporization (AH,) and their flow behavior, AH, and viscosity (7) were measured at different temperatures in the high molecular weight ester series: linear flexible di-n-alkyl sebacates and compact branched triglycerides with molecular weight ranging from 300 to 900. AHv" values (AHv corrected to 298 K) have been obtained with experimental AH, and also computed according to the group additivity method; a smaller-CH,- group value of 3.8 kJ mol-' compared to the normal value of 5.0 kJ mol-' is found to give good agreement with the experimental data (within 2-5% error). Both ester series have the same AH," irrespective of their molecular features, namely,shape, flexibility, and polarity, suggesting the coiling of the molecules during vaporization. The segmental motion of these ester series during their flow and its dependence on their molecular features unlike AH,' are demonstrated by the correlation of the enthalpy of activation for viscous flow (AH*) and the ratio AE,/AH* = n (AE, is the energy of vaporization) with molecular weight.