CFD modeling of single-phase flow in a packed bed with MRI validation

Computational fluid dynamics (CFD) simulations are becoming increasingly widespread with the advent of more powerful computers and more sophisticated software. The aim of these developments is to facilitate more accurate reactor design and optimization methods compared to traditional lumped-parameter models. However, in order for CFD to be a trusted method, it must be validated using experimental data acquired at sufficiently high spatial resolution. This article validates an in-house CFD code by comparison with flow-field data obtained using magnetic resonance imaging (MRI) for a packed bed with a particle-to-column diameter ratio of 2. Flows characterized by inlet Reynolds numbers, based on particle diameter, of 27, 55, 111, and 216 are considered. The code used employs preconditioning to directly solve for pressure in low-velocity flow regimes. Excellent agreement was found between the MRI and CFD data with relative error between the experimentally determined and numerically predicted flow-fields being in the range of 3-9%. © 2012 American Institute of Chemical Engineers (AIChE).

Numerical computation of laminar flow in a heated rotating cavity with an axial throughflow of air

A heated rotating cavity with an axial throughflow of cooling air is used as a model for the flow in the cylindrical cavities between adjacent discs of a high-pressure gas-turbine compressor. In an engine the flow is expected to be turbulent, the limitations of this laminar study are fully realised but it is considered an essential step to understand the fundamental nature of the flow. The three-dimensional, time-dependent governing equations are solved using a code based on the finite volume technique and a multigrid algorithm. The computed flow structure shows that flow enters the cavity in one or more radial arms and then forms regions of cyclonic and anticyclonic circulation. This basic flow structure is consistent with existing experimental evidence obtained from flow visualization. The flow structure also undergoes cyclic changes with time. For example, a single radial arm, and pair of recirculation regions can commute to two radial arms and two pairs of recirculation regions and then revert back to one. The flow structure inside the cavity is found to be heavily influenced by the radial distribution of surface temperature imposed on the discs. As the radial location of the maximum disc temperature moves radially outward, this appears to increase the number of radial arms and pairs of recirculation regions (from one to three for the distributions considered here). If the peripheral shroud is also heated there appear to be many radial arms which exchange fluid with a strong cyclonic flow adjacent to the shroud. One surface temperature distribution is studied in detail and profiles of the relative tangential and radial velocities are presented. The disc heat transfer is also found to be influenced by the disc surface temperature distribution. It is also found that the computed Nusselt numbers are in reasonable accord over most of the disc surface with a correlation found from previous experimental measurements. © 1994, MCB UP Limited.

Scaling of transitionally-rough flow structures in DNSes of channels with riblets up to Reτ=550

The linear, drag-reducing effect of vanishingly small riblets breaks down once their size is in the transitionally-rough regime. We have previously reported that this breakdown is caused by the additional Reynolds stresses produced by the appearance of elongated spanwise rollers just above the riblet surface. These rollers are related with the Kelvin--Helmholtz instability of free shear layers, and to similar structures appearing over other rough and porous surfaces. However, because of the limited Reτ=180 in our previous DNSes, it could not be determined whether those structures scaled in inner or outer units. Furthermore, it is questionable if results in the transitionally-rough regime at Reτ=180 can be extrapolated to configurations of practical interest. At such small Reynolds numbers, roughness of transitional size can perturb a large portion of the boundary layer, which is not the case in most industrial and atmospheric applications. To clarify these issues we have conducted a set of DNSes at Reτ=550. Our results indicate that the spanwise rollers scale in wall units, and support the validity of the extrapolation to configurations of practical interest.

The acoustic analogy in an annular duct with swirling mean flow

This paper is concerned with modelling the effects of swirling flow on turbomachinery noise. We develop an acoustic analogy to predict sound generation in a swirling and sheared base flow in an annular duct, including the presence of moving solid surfaces to account for blade rows. In so doing we have extended a number of classical earlier results, including Ffowcs Williams & Hawkings' equation in a medium at rest with moving surfaces, and Lilley's equation for a sheared but non-swirling jet. By rearranging the Navier-Stokes equations we find a single equation, in the form of a sixth-order differential operator acting on the fluctuating pressure field on the left-hand side and a series of volume and surface source terms on the right-hand side; the form of these source terms depends strongly on the presence of swirl and radial shear. The integral form of this equation is then derived, using the Green's function tailored to the base flow in the (rigid) duct. As is often the case in duct acoustics, it is then convenient to move into temporal, axial and azimuthal Fourier space, where the Green's function is computed numerically. This formulation can then be applied to a number of turbomachinery noise sources. For definiteness here we consider the noise produced downstream when a steady distortion flow is incident on the fan from upstream, and compare our results with those obtained using a simplistic but commonly used Doppler correction method. We show that in all but the simplest case the full inclusion of swirl within an acoustic analogy, as described in this paper, is required. © 2013 Cambridge University Press.

A thermopile based SOI CMOS MEMS wall shear stress sensor

In this paper we present for the first time, a novel silicon on insulator (SOI) complementary metal oxide semiconductor (CMOS) MEMS thermal wall shear stress sensor based on a tungsten hot-film and three thermopiles. These devices have been fabricated using a commercial 1 μm SOI-CMOS process followed by a deep reactive ion etch (DRIE) back-etch step to create silicon oxide membranes under the hot-film for effective thermal isolation. The sensors show an excellent repeatability of electro-thermal characteristics and can be used to measure wall shear stress in both constant current anemometric as well as calorimetric modes. The sensors have been calibrated for wall shear stress measurement of air in the range of 0-0.48 Pa using a suction type, 2-D flow wind tunnel. The calibration results show that the sensors have a higher sensitivity (up to four times) in calorimetric mode compared to anemometric mode for wall shear stress lower than 0.3 Pa. © 2013 IEEE.

Modelling and validation of wet steam flow in a low pressure steam turbine

Results of numerical investigations of the wet steam flow in a three stage low pressure steam turbine test rig are presented. The test rig is a scale model of a modern steam turbine design and provides flow measurements over a range of operating conditions which are used for detailed comparisons with the numerical results. For the numerical analysis a modern CFD code with user defined models for specific wet steam modelling is used. The effect of different theoretical models for nucleation and droplet growth are examined. It is shown that heterogeneous condensation is highly dependent on steam quality and, in this model turbine with high quality steam, a homogeneous theory appears to be the best choice. The homogeneous theory gives good agreement between the test rig traverse measurements and the numerical results. The differences in the droplet size distribution of the three stage turbine are shown for different loads and modelling assumptions. The different droplet growth models can influence the droplet size by a factor of two. An estimate of the influence of unsteady effects is made by means of an unsteady two-dimensional simulation. The unsteady modelling leads to a shift of nucleation into the next blade row. For the investigated three stage turbine the influence due to wake chopping on the condensation process is weak but to confirm this conclusion further investigations are needed in complete three dimensions and on turbines with more stages. Copyright © 2011 by ASME.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Strain rate-dependent shear band behavior in bulk metallic glasses

In this article, we review our recent advances in understanding the deformation behavior of a typical tough Zr41.2Ti13.8Cu12.5Ni10Be22.5 (Vit 1) bulk metallic glass (BMG), as a model material, under various loading modes and strain rates, focusing particularly on the rate-dependence and formation mechanism of shear-banding. Dynamic and quasi-static mechanical experiments, including plate shear, shear punch and spherical indentation, and continuum as well as atomistic modeling on shear-banding are discussed. The results demonstrate that higher strain rate slows down the annihilation process of free volume, but promotes the free-volume coalescence, which is responsible for the rate-dependent shear banding. The physical origin of shear bands, that is the free volume softening underpinned by irreversible rearrangements of atoms, is unveiled. Finally, some concluding remarks are given.

Quasi-3D Cytoskeletal Dynamics of Osteocytes under Fluid Flow

Osteocytes respond to dynamic fluid shear loading by activating various biochemical pathways, mediating a dynamic process of bone formation and resorption. Whole-cell deformation and regional deformation of the cytoskeleton may be able to directly regulate this process. Attempts to image cellular deformation by conventional microscopy techniques have been hindered by low temporal or spatial resolution. In this study, we developed a quasi-three-dimensional microscopy technique that enabled us to simultaneously visualize an osteocyte's traditional bottom-view profile and a side-view profile at high temporal resolution. Quantitative analysis of the plasma membrane and either the intracellular actin or microtubule (MT) cytoskeletal networks provided characterization of their deformations over time. Although no volumetric dilatation of the whole cell was observed under flow, both the actin and MT networks experienced primarily tensile strains in all measured strain components. Regional heterogeneity in the strain field of normal strains was observed in the actin networks, especially in the leading edge to flow, but not in the MT networks. In contrast, side-view shear strains exhibited similar subcellular distribution patterns in both networks. Disruption of MT networks caused actin normal strains to decrease, whereas actin disruption had little effect on the MT network strains, highlighting the networks' mechanical interactions in osteocytes.

Non-modal instabilities of two-dimensional disturbances in plane Couette flow of a power-law fluid

Instabilities of fluid flows have traditionally been investigated by normal mode analysis, i.e. by linearizing the equations of flow and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In this paper we study the instabilities of two-dimensional Couette flow of a polymeric fluid in the framework of non-modal stability theory rather than normal mode analysis. A power-law model is used to describe the polymeric liquid. We focus on the response to external excitations and initial conditions by examining the pseudospectra structures and the transient energy growths. For both Newtonian and non-Newtonian flows, the results show that there can be a rather large transient growth even though the linear operator of Couette flow has no unstable eigenvalue. The effects of non-Newtonian viscosity on the transient behaviors are examined in this study. The results show that the "shear-thinning/shear-thickening" effect increases/decreases the amplitude of responses to external excitations and initial conditions. (C) 2010 Elsevier B.V. All rights reserved.

Shear-enhanced crystallization of a thermoplastic polyimide derived from 3,3 ',4,4 '-oxydiphthalic dianhydride and 4,4 '-oxydianiline

The aim of this work is to investigate the effect of consecutive shear on the crystallization of an amorphous aromatic polyimide (PI) derived from 3,3',4,4'oxydiphthalic dianhydride (3,3',4,4'-ODPA) and 4,4-oxydianiline (ODA). At 260 degrees C, the increase of shear rate or shear time leads to the increase of crystallinity. Indeed, increasing shear rate can also accelerate the crystallization behavior. Moreover, it was found that a new melting peak appeared at higher temperature for long time or high rate sheared sample. The enhancement of crystallization behavior appears directly linked to the increase of crystal thickness. Particularly, the effect of shear temperature was investigated, and the results revealed that the crystallization of the PI was more sensitive to shear at 260 degrees C, which was 10 degrees above the glass transition temperature (250 degrees C) of the PI. Possible mechanism was proposed to illustrate the effect of consecutive shear on the crystallization of the PI polymer.

Influence of shear on polypropylene crystallization kinetics

Isothermal crystallization kinetics under shear in the melt of iPP was investigated by optical microscopy. It appears that shearing from 200 to the crystallization temperatures enhanced the kinetics, but the shear effect was not obvious if the melt of iPP was sheared only at 200. The experiment results show that relaxation plays an important role during crystallization, and that spherulite growth rates increased with shear rates and were governed by relaxation. The effect of flow on the crystallization kinetics can be understood by considering that the increase of the degree of order due to flow results is an effective change of the melt free energy. The Laurizen-Hoffman theory and the DE-IAA model were used to describe the shear-induced crystallization kinetics of iPP excellently.

Influence of shear on crystallization behavior of the beta phase in isotactic polypropylene with beta-nucleating agent

Wide-angle X-ray diffraction (WAXD) was used to investigate the effects of shear on the crystallization behavior of polypropylene (PP) with beta-nucleating agent. The melt was subjected to shear at the shear rate from 0.5 to 60 s(-1) for 5 s with a CSS450 shear stage. For the PP with low content of the additive, the formation mechanism of the beta crystals is almost the same as that of pure isotactic polypropylene (iPP), viz., shear induces. Otherwise, for the samples with high content of the additive, the formation mechanism of the beta form are nucleating agent induces. The results clearly show that shear restrains the formation of high beta phase for the melt with additive.

Polymer concentration, shear and stretch field effects on the surface morphology evolution during the spin-coating

Polymer concentration and shear and stretch field effects on the surface morphology evolution of three different kinds of polymers (polystyrene (PS), polybutadiene (PB) and polystyrene-b-polybutadiene-b-polystyrene (SBS)) during the spin-coating were investigated by means of atomic force microscopy (AFM). For PS and SBS, continuous film, net-like structure and particle structure were observed at different concentrations. For PB, net-like structures were not observed and continuous films and radial array of droplets emerged. Moreover, we compared surface morphology transitions on different substrate locations from the center to the edge. For PS, net-like structure, broken net-like structure and irregular array of particles were observed. For SBS, net-like structure, periodically orientated string-like structure and broken-line structure appeared. But for PB, flower-like holes in the continuous film, distorted stream-like structure and irregular distributions of droplets emerged. These different transitions of surface morphologies were discussed in terms of individual material property.

Theoretical representation of shear and pressure influence on the phase separation behavior of PVME/PS mixture

In the framework of lattice fluid model, the Gibbs energy and equation of state are derived by introducing the energy (E-s) stored during flow for polymer blends under shear. From the calculation of the spinodal of poly(vinyl methyl ether) (PVME) and polystyrene (PS) mixtures, we have found the influence of E., an equation of state in pure component is inappreciable, but it is appreciable in the mixture. However, the effect of E, on phase separation behavior is extremely striking. In the calculation of spinodal for the PVME/PS system, a thin, long and banana miscibility gap generated by shear is seen beside the miscibility gap with lower critical solution temperature. Meanwhile, a binodal coalescence of upper and lower miscibility gaps is occurred. The three points of the three-phase equilibrium are forecasted. The shear rate dependence of cloud point temperature at a certain composition is discussed. The calculated results are acceptable compared with the experiment values obtained by Higgins et at. However, the maximum positive shift and the minimum negative shift of cloud point temperature guessed by Higgins are not obtained, Furthermore, the combining effects of pressure and shear on spinodal shift are predicted.