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.

DEVELOPMENT OF STATIC MIXERS FOR MISCIBLE FLUIDS IN LAMINAR FLOW WITH THE USE OF COMPUTATIONAL FLUID DYNAMICS (CFD)

Static mixers with improved performance were developed from CFD simulations in a stepwise approach. The relevant geometric features of simple mixer designs and the corresponding mixing mechanisms-laminar shear, elongational flow, and distributive mixing-were identified first. This information was used to formulate guidelines for the development of new geometries. The solid elements of the static mixer should: (a) provide restrictions to the flow; (b) deflect the flow; (c) be sequentially rotated around the flow direction to provide symmetry; (d) extend from the center of the pipe to the vicinity of the walls to avoid short-circuiting; and (e) distribute and remix the flow. Based on these guidelines, two improved mixer designs were developed: the DS A-I mixer has a good mixing efficiency and an acceptable pressure drop; the Fins 35 degrees mixer is more efficient and compact, but requires a larger pressure drop. Their performance indicates that their use is possible on industrial applications.

Computational fluid dynamics investigation of a centrifugal blood pump

In the development of a ventricular assist device, computational fluid dynamics (CFD) analysis is an efficient tool to obtain the best design before making the final prototype. In this study, different designs of a centrifugal blood pump were developed to investigate flow characteristics and performance. This study assumed the blood flow as being an incompressible homogeneous Newtonian fluid. A constant velocity was applied at the inlet; no slip boundary conditions were applied at device wall; and pressure boundary conditions were applied at the outlet. The CFD code used in this work was based on the finite volume method. In the future, the results of CFD analysis can be compared with flow visualization and hemolysis tests.

Characterization of Human Late Outgrowth Endothelial Progenitor-Derived Cells under Various Flow Conditions.

Background: Endothelial progenitor-derived cells (EPC) are a cell therapy tool in peripheral arterial disease and for re-endothelialization of bypasses and stents. Objective: To assess EPC behavior under flow conditions normally found in vivo. Results: EPC were isolated from human cord blood, cultured on compliant tubes and exposed in an in vitro flow system mimicking hemodynamic environments normally found in medium and large arteries. EPC exposed for 24 h to unidirectional (0.3 ± 0.1 or 6 ± 3 dynes/cm(2)) shear stress oriented along flow direction, while those exposed to bidirectional shear stress (0.3 ± 3 dynes/cm(2)) or static conditions had random orientation. Under bidirectional flow, tissue factor (TF) activity and mRNA expression were significantly increased (2.5- and 7.0-fold) compared to static conditions. Under low shear unidirectional flow TF mRNA increased 4.9 ± 0.5-fold. Similar flow-induced increases were observed for TF in mature umbilical vein-derived endothelial cells. Expression of tissue-type plasminogen activator (t-PA), urokinase (u-PA) and monocyte chemotactic protein 1 (MCP1) were reduced by 40-60% in late outgrowth endothelial progenitor-derived cells (LO-EPC) exposed to any flow environment, while MCP1, but not t-PA or u-PA, was decreased in HUVEC. Conclusions: Flow, in particular bidirectional, modifies the hemostatic balance in LO-EPC with increased TF and decreased plasminogen activator expression.

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.

Drag reduction by polyethylene glycol in the tail arterial bed of normotensive and hypertensive rats

This study was designed to evaluate the effect of drag reducer polymers (DRP) on arteries from normotensive (Wistar) and spontaneously hypertensive rats (SHR). Polyethylene glycol (PEG 4000 at 5000 ppm) was perfused in the tail arterial bed with (E+) and without endothelium (E-) from male, adult Wistar (N = 14) and SHR (N = 13) animals under basal conditions (constant flow at 2.5 mL/min). In these preparations, flow-pressure curves (1.5 to 10 mL/min) were constructed before and 1 h after PEG 4000 perfusion. Afterwards, the tail arterial bed was fixed and the internal diameters of the arteries were then measured by microscopy and drag reduction was assessed based on the values of wall shear stress (WSS) by computational simulation. In Wistar and SHR groups, perfusion of PEG 4000 significantly reduced pulsatile pressure (Wistar/E+: 17.5 ± 2.8; SHR/E+: 16.3 ± 2.7%), WSS (Wistar/E+: 36; SHR/E+: 40%) and the flow-pressure response. The E- reduced the effects of PEG 4000 on arteries from both groups, suggesting that endothelial damage decreased the effect of PEG 4000 as a DRP. Moreover, the effects of PEG 4000 were more pronounced in the tail arterial bed from SHR compared to Wistar rats. In conclusion, these data demonstrated for the first time that PEG 4000 was more effective in reducing the pressure-flow response as well as WSS in the tail arterial bed of hypertensive than of normotensive rats and these effects were amplified by, but not dependent on, endothelial integrity. Thus, these results show an additional mechanism of action of this polymer besides its mechanical effect through the release and/or bioavailability of endothelial factors.

Applying discrete element modelling to vertical and horizontal shaft impact crushers

The PFC3D (particle flow code) that models the movement and interaction of particles by the DEM techniques was employed to simulate the particle movement and to calculate the velocity and energy distribution of collision in two types of impact crusher: the Canica vertical shaft crusher and the BJD horizontal shaft swing hammer mill. The distribution of collision energies was then converted into a product size distribution for a particular ore type using JKMRC impact breakage test data. Experimental data of the Canica VSI crusher treating quarry and the BJD hammer mill treating coal were used to verify the DEM simulation results. Upon the DEM procedures being validated, a detailed simulation study was conducted to investigate the effects of the machine design and operational conditions on velocity and energy distributions of collision inside the milling chamber and on the particle breakage behaviour. (C) 2003 Elsevier Ltd. All rights reserved.

Ex vivo Pulsatile Perfusion of Human Saphenous Veins Induces Intimal Hyperplasia and Increased Levels of the Plasminogen Activator Inhibitor 1.

Vessel wall trauma induces vascular remodeling processes including the development of intimal hyperplasia (IH). To assess the development of IH in human veins, we have used an ex vivo vein support system (EVVSS) allowing the perfusion of freshly isolated segments of saphenous veins in the presence of a pulsatile flow which reproduced arterial conditions regarding shear stress, flow rate and pressure during a period of 7 and 14 days. Compared to the corresponding freshly harvested human veins, histomorphometric analysis showed a significant increase in the intimal thickness which was already maximal after 7 days of perfusion. Expression of the endothelial marker CD31 demonstrated the presence of endothelium up to 14 days of perfusion. In our EVVSS model, the activity as well as the mRNA and protein expression levels of plasminogen activator inhibitor 1, the inhibitor of urokinase-type plasminogen activator (uPA) and tissue-type plasminogen activator (tPA), were increased after 7 days of perfusion, whereas the expression levels of tPA and uPA were not altered. No major change was observed between 7 and 14 days of perfusion. These data show that our newly developed EVVSS is a valuable setting to study ex vivo remodeling of human veins submitted to a pulsatile flow.

Drag reduction by polyethylene glycol in the tail arterial bed of normotensive and hypertensive rats

This study was designed to evaluate the effect of drag reducer polymers (DRP) on arteries from normotensive (Wistar) and spontaneously hypertensive rats (SHR). Polyethylene glycol (PEG 4000 at 5000 ppm) was perfused in the tail arterial bed with (E+) and without endothelium (E-) from male, adult Wistar (N = 14) and SHR (N = 13) animals under basal conditions (constant flow at 2.5 mL/min). In these preparations, flow-pressure curves (1.5 to 10 mL/min) were constructed before and 1 h after PEG 4000 perfusion. Afterwards, the tail arterial bed was fixed and the internal diameters of the arteries were then measured by microscopy and drag reduction was assessed based on the values of wall shear stress (WSS) by computational simulation. In Wistar and SHR groups, perfusion of PEG 4000 significantly reduced pulsatile pressure (Wistar/E+: 17.5 ± 2.8; SHR/E+: 16.3 ± 2.7%), WSS (Wistar/E+: 36; SHR/E+: 40%) and the flow-pressure response. The E- reduced the effects of PEG 4000 on arteries from both groups, suggesting that endothelial damage decreased the effect of PEG 4000 as a DRP. Moreover, the effects of PEG 4000 were more pronounced in the tail arterial bed from SHR compared to Wistar rats. In conclusion, these data demonstrated for the first time that PEG 4000 was more effective in reducing the pressure-flow response as well as WSS in the tail arterial bed of hypertensive than of normotensive rats and these effects were amplified by, but not dependent on, endothelial integrity. Thus, these results show an additional mechanism of action of this polymer besides its mechanical effect through the release and/or bioavailability of endothelial factors.