Adaption and use of a compressible flow code for turbomachinery design

The present paper presents and discusses the use of dierent codes regarding the numerical simulation of a radial-in ow turbine. A radial-in ow turbine test case was selected from published literature [1] and commercial codes (Fluent and CFX) were used to perform the steady-state numerical simulations. An in-house compressible- ow simulation code, Eilmer3 [2] was also adapted in order to make it suitable to perform turbomachinery simulations and preliminary results are presented and discussed. The code itself as well as its adaptation, comprising the addition of terms for the rotating frame of reference, programmable boundary conditions for periodic boundaries and a mixing plane interface between the rotating and non-rotating blocks are also discussed. Several cases with dierent orders of complexity in terms of geometry were considered and the results were compared across the dierent codes. The agreement between these results and published data is also discussed.

The implementation of an aeronautical CFD flow code onto distributed memory parallel systems

The parallelization of an industrially important in-house computational fluid dynamics (CFD) code for calculating the airflow over complex aircraft configurations using the Euler or Navier–Stokes equations is presented. The code discussed is the flow solver module of the SAUNA CFD suite. This suite uses a novel grid system that may include block-structured hexahedral or pyramidal grids, unstructured tetrahedral grids or a hybrid combination of both. To assist in the rapid convergence to a solution, a number of convergence acceleration techniques are employed including implicit residual smoothing and a multigrid full approximation storage scheme (FAS). Key features of the parallelization approach are the use of domain decomposition and encapsulated message passing to enable the execution in parallel using a single programme multiple data (SPMD) paradigm. In the case where a hybrid grid is used, a unified grid partitioning scheme is employed to define the decomposition of the mesh. The parallel code has been tested using both structured and hybrid grids on a number of different distributed memory parallel systems and is now routinely used to perform industrial scale aeronautical simulations. Copyright © 2000 John Wiley & Sons, Ltd.

Parallelisation of ESAUNA within EUROPORT-1: A structured/unstructured aeronautical CFD flow code

General-purpose parallel processing for solving day-to-day industrial problems has been slow to develop, partly because of the lack of suitable hardware from well-established, mainstream computer manufacturers and suitably parallelized application software. The parallelization of a CFD-(computational fluid dynamics) flow solution code is known as ESAUNA. This code is part of SAUNA, a large CFD suite aimed at computing the flow around very complex aircraft configurations including complete aircraft. A novel feature of the SAUNA suite is that it is designed to use either block-structured hexahedral grids, unstructured tetrahedral grids, or a hybrid combination of both grid types. ESAUNA is designed to solve the Euler equations or the Navier-Stokes equations, the latter in conjunction with various turbulence models. Two fundamental parallelization concepts are used—namely, grid partitioning and encapsulation of communications. Grid partitioning is applied to both block-structured grid modules and unstructured grid modules. ESAUNA can also be coupled with other simulation codes for multidisciplinary computations such as flow simulations around an aircraft coupled with flutter prediction for transient flight simulations.

AN ASSESSMENT OF FLOW REGIME MAPS AND A NUMERICAL HEAT TRANSFER CORRELATION FOR THE STRATIFIED/ANNULAR REGIME OF SHEAR-DRIVEN INTERNAL CONDENSING FLOWS

Several modern-day cooling applications require the incorporation of mini/micro-channel shear-driven flow condensers. There are several design challenges that need to be overcome in order to meet those requirements. The difficulty in developing effective design tools for shear-driven flow condensers is exacerbated due to the lack of a bridge between the physics-based modelling of condensing flows and the current, popular approach based on semi-empirical heat transfer correlations. One of the primary contributors of this disconnect is a lack of understanding caused by the fact that typical heat transfer correlations eliminate the dependence of the heat transfer coefficient on the method of cooling employed on the condenser surface when it may very well not be the case. This is in direct contrast to direct physics-based modeling approaches where the thermal boundary conditions have a direct and huge impact on the heat transfer coefficient values. Typical heat transfer correlations instead introduce vapor quality as one of the variables on which the value of the heat transfer coefficient depends. This study shows how, under certain conditions, a heat transfer correlation from direct physics-based modeling can be equivalent to typical engineering heat transfer correlations without making the same apriori assumptions. Another huge factor that raises doubts on the validity of the heat-transfer correlations is the opacity associated with the application of flow regime maps for internal condensing flows. It is well known that flow regimes influence heat transfer rates strongly. However, several heat transfer correlations ignore flow regimes entirely and present a single heat transfer correlation for all flow regimes. This is believed to be inaccurate since one would expect significant differences in the heat transfer correlations for different flow regimes. Several other studies present a heat transfer correlation for a particular flow regime - however, they ignore the method by which extents of the flow regime is established. This thesis provides a definitive answer (in the context of stratified/annular flows) to: (i) whether a heat transfer correlation can always be independent of the thermal boundary condition and represented as a function of vapor quality, and (ii) whether a heat transfer correlation can be independently obtained for a flow regime without knowing the flow regime boundary (even if the flow regime boundary is represented through a separate and independent correlation). To obtain the results required to arrive at an answer to these questions, this study uses two numerical simulation tools - the approximate but highly efficient Quasi-1D simulation tool and the exact but more expensive 2D Steady Simulation tool. Using these tools and the approximate values of flow regime transitions, a deeper understanding of the current state of knowledge in flow regime maps and heat transfer correlations in shear-driven internal condensing flows is obtained. The ideas presented here can be extended for other flow regimes of shear-driven flows as well. Analogous correlations can also be obtained for internal condensers in the gravity-driven and mixed-driven configuration.

Yield stress, thixotropy and shear banding in a dilute aqueous suspension of few layer graphene oxide platelets

We demonstrate a rigidity percolation transition and the onset of yield stress in a dilute aqueous dispersion of graphene oxide platelets (aspect ratio similar to 5000) above a critical volume fraction of 3.75 x 10(-4) with a percolation exponent of 2.4 +/- 0.1. The viscoelastic moduli of the gel at rest measured as a function of time indicate the absence of structural evolution of the 3D percolated network of disks. However a shear-induced aging giving rise to a compact jammed state and shear rejuvenation indicating a homogenous flow is observed when a steady shear stress (sigma) is imposed in creep experiments. We construct a shear diagram (sigma vs. volume fraction phi) and the critical stress above which shear rejuvenation occurs is identified as the yield stress sigma(y) of the gel. The minimum steady state shear rate (gamma) over dot(m) obtained from creep experiments agrees well with the end of the plateau region in a controlled shear rate flow curve, indicating a shear localization below (gamma) over dot(m). A steady state shear banding in the plateau region of the flow curve observed in particle velocimetry measurements in a Couette geometry confirms that the dilute suspensions of GO platelets form a thixotropic yield stress fluid.

Studies on two-phase co-current air/non-Newtonian shear-thinning fluid flows in inclined smooth pipes

In this work. co-current flow characteristics of air/non-Newtonian liquid systems in inclined smooth pipes are studied experimentally and theoretically using transparent tubes of 20, 40 and 60 turn in diameter. Each tube includes two 10 m lone pipe branches connected by a U-bend that is capable of being inclined to any angle, from a completely horizontal to a fully vertical position. The flow rate of each phase is varied over a wide range. The studied flow phenomena are bubbly, plug flow, slug flow, churn flow and annular flow. These are observed and recorded by a high flow. stratified flow. -speed camera over a wide range of operating conditions. The effects of the liquid phase properties, the inclination angle and the pipe diameter on two-phase flow characteristics are systematically studied. The Heywood-Charles model for horizontal flow was modified to accommodate stratified flow in inclined pipes, taking into account the average void fraction and pressure drop of the mixture flow of a gas/non-Newtonian liquid. The pressure drop gradient model of Taitel and Barnea for a gas/Newtonian liquid slug flow was extended to include liquids possessing shear-thinning flow behaviour in inclined pipes. The comparison of the predicted values with the experimental data shows that the models presented here provide a reasonable estimate of the average void fraction and the corresponding pressure drop for the mixture flow of a gas/ non-Newtonian liquid. (C) 2007 Elsevier Ltd. All rights reserved.

Time development of small disturbances to plane Couette Flow

The question of linear sheared-disturbance evolution in constant-shear parallel flow is here reexamined with regard to the temporary-amplification phenomenon noted first by Orr in 1907. The results apply directly to Rossby waves on a beta-plane, and are also relevant to the Eady model of baroclinic instability. It is shown that an isotropic initial distribution of standing waves maintains a constant energy level throughout the shearing process, the amplification of some waves being precisely balanced by the decay of the others. An expression is obtained for the energy of a distribution of disturbances whose wavevectors lie within a given angular wedge and an upper bound derived. It is concluded that the case for ubiquitous amplification made in recent studies may have been somewhat overstated: while carefully-chosen individual Fourier components can amplify considerably before they decay. a general distribution will tend to exhibit little or no amplification.

DOCK2 is required for chemokine-promoted human T lymphocyte adhesion under shear stress mediated by the integrin alpha4beta1

The alpha4beta1 integrin is an essential adhesion molecule for recruitment of circulating lymphocytes into lymphoid organs and peripheral sites of inflammation. Chemokines stimulate alpha4beta1 adhesive activity allowing lymphocyte arrest on endothelium and subsequent diapedesis. Activation of the GTPase Rac by the guanine-nucleotide exchange factor Vav1 promoted by CXCL12 controls T lymphocyte adhesion mediated by alpha4beta1. In this study, we investigated the role of DOCK2, a lymphocyte guanine-nucleotide exchange factor also involved in Rac activation, in CXCL12-stimulated human T lymphocyte adhesion mediated by alpha4beta1. Using T cells transfected with DOCK2 mutant forms defective in Rac activation or with DOCK2 small interfering RNA, we demonstrate that DOCK2 is needed for efficient chemokine-stimulated lymphocyte attachment to VCAM-1 under shear stress. Flow chamber, soluble binding, and cell spreading assays identified the strengthening of alpha4beta1-VCAM-1 interaction, involving high affinity alpha4beta1 conformations, as the adhesion step mainly controlled by DOCK2 activity. The comparison of DOCK2 and Vav1 involvement in CXCL12-promoted Rac activation and alpha4beta1-dependent human T cell adhesion indicated a more prominent role of Vav1 than DOCK2. These results suggest that DOCK2-mediated signaling regulates chemokine-stimulated human T lymphocyte alpha4beta1 adhesive activity, and that cooperation with Vav1 might be required to induce sufficient Rac activation for efficient adhesion. In contrast, flow chamber experiments using lymph node and spleen T cells from DOCK2(-/-) mice revealed no significant alterations in CXCL12-promoted adhesion mediated by alpha4beta1, indicating that DOCK2 activity is dispensable for triggering of this adhesion in mouse T cells, and suggesting that Rac activation plays minor roles in this process.

(Table 1) Results from XRD analysis as well as sediment type and selected results from ring shear experiments

The location of the seaward tip of a subduction thrust controls material transfer at convergent plate margins, and hence global mass balances. At approximately half of those margins, the material of the subducting plate is completely underthrust so that no accretion or even subduction erosion takes place. Along the remaining margins, material is scraped off the subducting plate and added to the upper plate by frontal accretion. We here examine the physical properties of subducting sediments off Costa Rica and Nankai, type examples for an erosional and an accretionary margin, to investigate which parameters control the level where the frontal thrust cuts into the incoming sediment pile. A series of rotary-shear experiments to measure the frictional strength of the various lithologies entering the two subduction zones were carried out. Results include the following findings: (1) At Costa Rica, clay-rich strata at the top of the incoming succession have the lowest strength (µres = 0.19) while underlying calcareous ooze, chalk and diatomite are strong (up to µres = 0.43; µpeak = 0.56). Hence the entire sediment package is underthrust. (2) Off Japan, clay-rich deposits within the lower Shikoku Basin inventory are weakest (µres = 0.13-0.19) and favour the frontal proto-thrust to migrate into one particular horizon between sandy, competent turbidites below and ash-bearing mud above. (3) Taking in situ data and earlier geotechnical testing into account, it is suggested that mineralogical composition rather than pore-pressure defines the position of the frontal thrust, which locates in the weakest, clay mineral-rich (up to 85 wt.%) materials. (4) Smectite, the dominant clay mineral phase at either margin, shows rate strengthening and stable sliding in the frontal 50 km of the subduction thrust (0.0001-0.1 mm/s, 0.5-25 MPa effective normal stress). (5) Progressive illitization of smectite cannot explain seismogenesis, because illite-rich samples also show velocity strengthening at the conditions tested.

Determination of lifter design, speed and filling effects in AG mills by 3D DEM

The power required to operate large gyratory mills often exceeds 10 MW. Hence, optimisation of the power consumption will have a significant impact on the overall economic performance and environmental impact of the mineral processing plant. In most of the published models of tumbling mills (e.g. [Morrell, S., 1996. Power draw of wet tumbling mills and its relationship to charge dynamics, Part 2: An empirical approach to modelling of mill power draw. Trans. Inst. Mining Metall. (Section C: Mineral Processing Ext. Metall.) 105, C54-C62. Austin, L.G., 1990. A mill power equation for SAG mills. Miner. Metall. Process. 57-62]), the effect of lifter design and its interaction with mill speed and filling are not incorporated. Recent experience suggests that there is an opportunity for improving grinding efficiency by choosing the appropriate combination of these variables. However, it is difficult to experimentally determine the interactions of these variables in a full scale mill. Although some work has recently been published using DEM simulations, it was basically. limited to 2D. The discrete element code, Particle Flow Code 3D (PFC3D), has been used in this work to model the effects of lifter height (525 cm) and mill speed (50-90% of critical) on the power draw and frequency distribution of specific energy (J/kg) of normal impacts in a 5 m diameter autogenous (AG) mill. It was found that the distribution of the impact energy is affected by the number of lifters, lifter height, mill speed and mill filling. Interactions of lifter design, mill speed and mill filling are demonstrated through three dimensional distinct element methods (3D DEM) modelling. The intensity of the induced stresses (shear and normal) on lifters, and hence the lifter wear, is also simulated. (C) 2004 Elsevier Ltd. All rights reserved.

Punching shear of slab-column connection in flat plate construction

The flat plate system is currently widely used in construction. It permits architectural flexibility, more clear space, less building height, easier formwork, and shorter construction time. However, there remains the problem of brittle punching failure due to the transfer of shearing forces and unbalanced moments at the flat plate-column connection. It is the purpose of this paper to investigate the effects of various interdependent factors that govern the punching shear resistance and behaviour of the flat plate-column connection, as well as their inclusion in current Codes.

Relative effect of slope and equilibrium line altitude on the retreat of Himalayan glaciers

It has been observed that a majority of glaciers in the Himalayas have been retreating. In this paper, we show that there are two major factors which control the advance/retreat of the Himalayan glaciers. They are the slope of the glacier and changes in the equilibrium line altitude. While it is well known, that these factors are important, we propose a new way of combining them and use it to predict retreat. The functional form of this model has been derived from numerical simulations using an ice-flow code. The model has been successfully applied to the movement of eight Himalayan glaciers during the past 25 years. It explains why the Gangotri glacier is retreating while Zemu of nearly the same length is stationary, even if they are subject to similar environmental changes. The model has also been applied to a larger set of glaciers in the Parbati basin, for which retreat based on satellite data is available, though over a shorter time period.

Numerical Study Of Pile-Up In Bulk Metallic Glass During Spherical Indentation

Pile-up around indenter is usually observed during instrumented indentation tests on bulk metallic glass. Neglecting the pile-up effect may lead to errors in evaluating hardness, Young's modulus, stress-strain response, etc. Finite element analysis was employed to implement numerical simulation of spherical indentation tests on bulk metallic glass. A new model was proposed to describe the pile-up effect. By using this new model, the contact radius and hardness of Zr41.2Ti13.8Cu12.5Ni10Be22.5 bulk metallic glass were obtained under several different indenter loads with pile-up, and the results agree well with the data generated by numerical simulation.

TEMPERATURE AND STRAIN-RATE DEPENDENCE OF FRACTURE-TOUGHNESS OF PHENOLPHTHALEIN POLYETHER KETONE

A strong strain-rate and temperature dependence was observed for the fracture toughness of phenolphthalein polyether ketone (PEK-C). Two separate crack-blunting mechanisms have been proposed to account for the fracture-toughness data. The first mechanism involves thermal blunting due to adiabatic heating at the crack tip for the high temperatures studied. In the high-temperature range, thermal blunting increases the fracture toughness corresponding to an effectively higher test temperature. However, in the low-temperature range, the adiabatic temperature rise is insufficient to cause softening and Jic increases with increasing temperature owing to viscoelastic losses associated with the p-relaxation there. The second mechanism involves plastic blunting due to shear yield/flow processes at the crack tip and this takes place at slow strain testing of the single-edge notched bending (SENB) samples. The temperature and strain-rate dependence of the plastic zone size may also be responsible for the temperature and strain-rate dependence of fracture toughness.