Reynolds-stress turbulence model in the KIVA code for engine simulation

A Reynolds-Stress Turbulence Model has been incorporated with success into the KIVA code, a computational fluid dynamics hydrocode for three-dimensional simulation of fluid flow in engines. The newly implemented Reynolds-stress turbulence model greatly improves the robustness of KIVA, which in its original version has only eddy-viscosity turbulence models. Validation of the Reynolds-stress turbulence model is accomplished by conducting pipe-flow and channel-flow simulations, and comparing the computed results with experimental and direct numerical simulation data. Flows in engines of various geometry and operating conditions are calculated using the model, to study the complex flow fields as well as confirm the model’s validity. Results show that the Reynolds-stress turbulence model is able to resolve flow details such as swirl and recirculation bubbles. The model is proven to be an appropriate choice for engine simulations, with consistency and robustness, while requiring relatively low computational effort.

Experimental results on gravity driven fully condensing flows in vertical tubes, their agreement with theory, and their differences with shear driven flows' boundary-condition sensitivities

This doctoral thesis presents the experimental results along with a suitable synthesis with computational/theoretical results towards development of a reliable heat transfer correlation for a specific annular condensation flow regime inside a vertical tube. For fully condensing flows of pure vapor (FC-72) inside a vertical cylindrical tube of 6.6 mm diameter and 0.7 m length, the experimental measurements are shown to yield values of average heat transfer co-efficient, and approximate length of full condensation. The experimental conditions cover: mass flux G over a range of 2.9 kg/m2-s ≤ G ≤ 87.7 kg/m2-s, temperature difference ∆T (saturation temperature at the inlet pressure minus the mean condensing surface temperature) of 5 ºC to 45 ºC, and cases for which the length of full condensation xFC is in the range of 0 < xFC < 0.7 m. The range of flow conditions over which there is good agreement (within 15%) with the theory and its modeling assumptions has been identified. Additionally, the ranges of flow conditions for which there are significant discrepancies (between 15 -30% and greater than 30%) with theory have also been identified. The paper also refers to a brief set of key experimental results with regard to sensitivity of the flow to time-varying or quasi-steady (i.e. steady in the mean) impositions of pressure at both the inlet and the outlet. The experimental results support the updated theoretical/computational results that gravity dominated condensing flows do not allow such elliptic impositions.

Bislug flow in circular and noncircular channels and the role of interface stretching on energy dissipation

The area of microfluidics has increased in popularity with such fields as MEMS, microreactors, microscaleheat exchangers, etc. A comprehensive understanding of dissipation mechanisms for fluid flow in microchannels is required to accurately predict the behavior in these small systems. Tests were performed using a constant pressure potential created by two immiscible fluids juxtaposed in a microchannel. This study focused on the flow and dissipation mechanisms in round and square microchannels. There are four major dissipation mechanisms in slug flow; wall shear, dissipation at the contact line, menisci interaction and the stretching of the interface. A force balance between the internal driving potential, viscous drag and interface stretching was used to develop a model for the prediction of the velocity of a bislug in a microchannel. Interface stretching is a dissipation mechanism that has been included due to the unique system properties and becomes increasingly more important as the bislug decreases in length.

Two dimensional Lattice Boltzmann simulation of fluid flow through an idealized micro-structure of disordered media

This technical report discusses the application of Lattice Boltzmann Method (LBM) in the fluid flow simulation through porous filter-wall of disordered media. The diesel particulate filter (DPF) is an example of disordered media. DPF is developed as a cutting edge technology to reduce harmful particulate matter in the engine exhaust. Porous filter-wall of DPF traps these soot particles in the after-treatment of the exhaust gas. To examine the phenomena inside the DPF, researchers are looking forward to use the Lattice Boltzmann Method as a promising alternative simulation tool. The lattice Boltzmann method is comparatively a newer numerical scheme and can be used to simulate fluid flow for single-component single-phase, single-component multi-phase. It is also an excellent method for modelling flow through disordered media. The current work focuses on a single-phase fluid flow simulation inside the porous micro-structure using LBM. Firstly, the theory concerning the development of LBM is discussed. LBM evolution is always related to Lattice gas Cellular Automata (LGCA), but it is also shown that this method is a special discretized form of the continuous Boltzmann equation. Since all the simulations are conducted in two-dimensions, the equations developed are in reference with D2Q9 (two-dimensional 9-velocity) model. The artificially created porous micro-structure is used in this study. The flow simulations are conducted by considering air and CO2 gas as fluids. The numerical model used in this study is explained with a flowchart and the coding steps. The numerical code is constructed in MATLAB. Different types of boundary conditions and their importance is discussed separately. Also the equations specific to boundary conditions are derived. The pressure and velocity contours over the porous domain are studied and recorded. The results are compared with the published work. The permeability values obtained in this study can be fitted to the relation proposed by Nabovati [8], and the results are in excellent agreement within porosity range of 0.4 to 0.8.

Flow-Centric, Back-In-Time Debugging

Conventional debugging tools present developers with means to explore the run-time context in which an error has occurred. In many cases this is enough to help the developer discover the faulty source code and correct it. However, rather often errors occur due to code that has executed in the past, leaving certain objects in an inconsistent state. The actual run-time error only occurs when these inconsistent objects are used later in the program. So-called back-in-time debuggers help developers step back through earlier states of the program and explore execution contexts not available to conventional debuggers. Nevertheless, even back-in-time debuggers do not help answer the question, ``Where did this object come from?'' The Object-Flow Virtual Machine, which we have proposed in previous work, tracks the flow of objects to answer precisely such questions, but this VM does not provide dedicated debugging support to explore faulty programs. In this paper we present a novel debugger, called Compass, to navigate between conventional run-time stack-oriented control flow views and object flows. Compass enables a developer to effectively navigate from an object contributing to an error back-in-time through all the code that has touched the object. We present the design and implementation of Compass, and we demonstrate how flow-centric, back-in-time debugging can be used to effectively locate the source of hard-to-find bugs.

Hemodynamics and Flow Characteristics of a New Dialysis Port

Renal replacement therapy by hemodialysis requires a permanent vascular access. Implantable ports offer a potential alternative to standard vascular access strategies although their development is limited both in number and extent. We explored the fluid dynamics within two new percutaneous bone-anchored dialysis port prototypes, both by in vitro experiments and computer simulation. The new port is to be fixed to bone and allows the connection of a dialysis machine to a central venous catheter via a built-in valve. We found that the pressure drop induced by the two ports was between 20 and 50 mmHg at 500 ml/min, which is comparable with commercial catheter connectors (15–80 mmHg). We observed the formation of vortices in both geometries, and a shear rate in the physiological range (<10,000s-1), which is lower than maximal shear rates reported in commercial catheters (up to 13,000s-1). A difference in surface shear rate of 15% between the two ports was obtained.

Impact of normal and shear stresses due to wheel slip on hydrological properties of an agricultural clay loam: Experimental and new computerized approach

The main purpose of this study was to evaluate the effect that mechanical stresses acting under the slipping driving wheels of agricultural equipment have on the soil’s pore system and water flow process (surface runoff generation during extreme event). The field experiment simulated low slip (1%) and high slip (27%) on a clay loam. The stress on the soil surface and changes in the amounts of water flowing from macropores were simulated using the Tires/tracks And Soil Compaction (TASC) tool and the MACRO model, respectively. Taking a 65 kW tractor on a clay loam as a reference, results showed that an increase in slip of the rear wheels from 1% to 27% caused normal stress to increase from 90.6 kPa to 104.4 kPa at the topsoil level, and the maximum shear contact stress to rise drastically from 6.0 kPa to 61.6 kPa. At 27% slip, topsoil was sheared and displaced over a distance of 0.35 m. Excessive normal and shear stress values with high slip caused severe reductions of the soil’s macroporosity, saturated hydraulic conductivity, and water quantities flowing from topsoil macropores. Assuming that, under conditions of intense rainfall on sloping land, a loss in vertical water flow would mean an increase in surface runoff, we calculated that a rainfall intensity of 100 mm h-1 and a rainfall duration of 1 h would increase the runoff coefficient to 0.79 at low slip and to 1.00 at high slip, indicating that 100% of rainwater would be transformed into surface runoff at high slip. We expect that these effects have a significant impact on soil erosion and floods in steeper terrain (slope > 15°) and across larger surface areas (> 16 m2) than those included in our study.

K-t GRAPPA-accelerated 4D flow MRI of liver hemodynamics: influence of different acceleration factors on qualitative and quantitative assessment of blood flow.

OBJECTIVE We sought to evaluate the feasibility of k-t parallel imaging for accelerated 4D flow MRI in the hepatic vascular system by investigating the impact of different acceleration factors. MATERIALS AND METHODS k-t GRAPPA accelerated 4D flow MRI of the liver vasculature was evaluated in 16 healthy volunteers at 3T with acceleration factors R = 3, R = 5, and R = 8 (2.0 × 2.5 × 2.4 mm(3), TR = 82 ms), and R = 5 (TR = 41 ms); GRAPPA R = 2 was used as the reference standard. Qualitative flow analysis included grading of 3D streamlines and time-resolved particle traces. Quantitative evaluation assessed velocities, net flow, and wall shear stress (WSS). RESULTS Significant scan time savings were realized for all acceleration factors compared to standard GRAPPA R = 2 (21-71 %) (p < 0.001). Quantification of velocities and net flow offered similar results between k-t GRAPPA R = 3 and R = 5 compared to standard GRAPPA R = 2. Significantly increased leakage artifacts and noise were seen between standard GRAPPA R = 2 and k-t GRAPPA R = 8 (p < 0.001) with significant underestimation of peak velocities and WSS of up to 31 % in the hepatic arterial system (p <0.05). WSS was significantly underestimated up to 13 % in all vessels of the portal venous system for k-t GRAPPA R = 5, while significantly higher values were observed for the same acceleration with higher temporal resolution in two veins (p < 0.05). CONCLUSION k-t acceleration of 4D flow MRI is feasible for liver hemodynamic assessment with acceleration factors R = 3 and R = 5 resulting in a scan time reduction of at least 40 % with similar quantitation of liver hemodynamics compared with GRAPPA R = 2.

Impact of local endothelial shear stress on neointima and plaque following stent implantation in patients with ST-elevation myocardial infarction: A subgroup-analysis of the COMFORTABLE AMI-IBIS 4 trial.

BACKGROUND Numerous studies have demonstrated an association between endothelial shear stress (ESS) and neointimal formation after stent implantation. However, the role of ESS on the composition of neointima and underlying plaque remains unclear. METHODS Patients recruited in the Comfortable AMI-IBIS 4 study implanted with bare metal stents (BMS) or biolimus eluting stents (BES) that had biplane coronary angiography at 13month follow-up were included in the analysis. The intravascular ultrasound virtual-histology (IVUS-VH) and the angiographic data were used to reconstruct the luminal surface, and the stent in the stented segments. Blood flow simulation was performed in the stent surface, which was assumed to represent the luminal surface at baseline, to assess the association between ESS and neointima thickness. The predominant ESS was estimated in 3-mm segments and was correlated with the amount of neointima, neointimal tissue composition, and with the changes in the underlying plaque burden and composition. RESULTS Forty three patients (18 implanted with BMS and 25 with BES) were studied. In both stent groups negative correlations were noted between ESS and neointima thickness in BMS (P<0.001) and BES (P=0.002). In BMS there was a negative correlation between predominant ESS and the percentage of the neointimal necrotic core component (P=0.015). In BES group, the limited neointima formation did not allow evaluation of the effect of ESS on its tissue characteristics. ESS did not affect vessel wall remodeling and the plaque burden and composition behind BMS (P>0.10) and BES (P>0.45). CONCLUSIONS ESS determines neointimal formation in both BMS and BES and affects the composition of the neointima in BMS. Conversely, ESS does not impact the plaque behind struts irrespective of stent type throughout 13months of follow-up.

Structural, metamorphic and geochronological relations between the Zanskar Shear Zone and the Miyar Shear Zone (NW Indian Himalaya): Evidence for two distinct tectonic structures and implications for the evolution of the High Himalayan Crystalline of Zanskar

The geologic structures and metamorphic zonation of the northwestern Indian Himalaya contrast significantly with those in the central and eastern parts of the range, where the high-grade metamorphic rocks of the High Himalayan Crystalline (HHC) thrust southward over the weakly metamorphosed sediments of the Lesser Himalaya along the Main Central Thrust (MCT). Indeed, the hanging wall of the MCT in the NW Himalaya mainly consists of the greenschist facies metasediments of the Chamba zone, whereas HHC high-grade rocks are exposed more internally in the range as a large-scale dome called the Gianbul dome. This Gianbul dome is bounded by two oppositely directed shear zones, the NE-dipping Zanskar Shear Zone (ZSZ) on the northern flank and the SW-dipping Miyar Shear Zone (MSZ) on the southern limb. Current models for the emplacement of the HHC in NW India as a dome structure differ mainly in terms of the roles played by both the ZSZ and the MSZ during the tectonothermal evolution of the HHC. In both the channel flow model and wedge extrusion model, the ZSZ acts as a backstop normal fault along which the high-grade metamorphic rocks of the HHC of Zanskar are exhumed. In contrast, the recently proposed tectonic wedging model argues that the ZSZ and the MSZ correspond to one single detachment system that operates as a subhorizontal backthrust off of the MCT. Thus, the kinematic evolution of the two shear zones, the ZSZ and the MSZ, and their structural, metamorphic and chronological relations appear to be diagnostic features for discriminating the different models. In this paper, structural, metamorphic and geochronological data demonstrate that the MSZ and the ZSZ experienced two distinct kinematic evolutions. As such, the data presented in this paper rule out the hypothesis that the MSZ and the ZSZ constitute one single detachment system, as postulated by the tectonic wedging model. Structural, metamorphic and geochronological data are used to present an alternative tectonic model for the large-scale doming in the NW Indian Himalaya involving early NE-directed tectonics, weakness in the upper crust, reduced erosion at the orogenic front and rapid exhumation along both the ZSZ and the MSZ.

Granular flow in polymineralic rocks bearing sheet silicates: New evidence from natural examples

Natural deformation in carbonate mylonites bearing sheet silicates occurs via a complex interaction of granular flow and solution transfer processes and involves continuous cycles of dissolution, grain boundary diffusion, nucleation and growth. In this way, new sheet silicates (a) nucleate within voids formed by grain boundary sliding of calcite grains. (b) grow, and (c) rotate towards the shear plane. As a consequence, small mica grains show a wide range of orientations with respect to the shear plane, but moderate to large grains are subparallel both to each other and to the shear plane. Increases of average grain sizes with increasing temperature of sheet silicates in mica-rich layers is more pronounced than in mica-poor layers. In the calcitic matrix however, sheet silicates can only grow via solution-precipitation and mass transfer processes. Therefore, the observed grain size variability indicates drastic differences in mass transfer behavior between the individual layers, which might be related to differences in the fluid flux. Based on these observations, a conceptual model for the microfabric evolution in sheet silicate bearing mylonites is presented. © 2001 Elsevier Science B.V. All rights reserved.

FOXC2 and fluid shear stress stabilize postnatal lymphatic vasculature

Biomechanical forces, such as fluid shear stress, govern multiple aspects of endothelial cell biology. In blood vessels, disturbed flow is associated with vascular diseases, such as atherosclerosis, and promotes endothelial cell proliferation and apoptosis. Here, we identified an important role for disturbed flow in lymphatic vessels, in which it cooperates with the transcription factor FOXC2 to ensure lifelong stability of the lymphatic vasculature. In cultured lymphatic endothelial cells, FOXC2 inactivation conferred abnormal shear stress sensing, promoting junction disassembly and entry into the cell cycle. Loss of FOXC2-dependent quiescence was mediated by the Hippo pathway transcriptional coactivator TAZ and, ultimately, led to cell death. In murine models, inducible deletion of Foxc2 within the lymphatic vasculature led to cell-cell junction defects, regression of valves, and focal vascular lumen collapse, which triggered generalized lymphatic vascular dysfunction and lethality. Together, our work describes a fundamental mechanism by which FOXC2 and oscillatory shear stress maintain lymphatic endothelial cell quiescence through intercellular junction and cytoskeleton stabilization and provides an essential link between biomechanical forces and endothelial cell identity that is necessary for postnatal vessel homeostasis. As FOXC2 is mutated in lymphedema-distichiasis syndrome, our data also underscore the role of impaired mechanotransduction in the pathology of this hereditary human disease.

Chemistry of minerals from amphibole-chlorite veins of ODP Site 209-1270

We examined small-scale shear zones in drillcore samples of abyssal peridotites from the Mid-Atlantic Ridge. These shear zones are associated with veins consisting of chlorite + actinolite/tremolite assemblages, with accessory phases zircon and apatite, and they are interpreted as altered plagiogranite melt impregnations, which originate from hydrous partial melting of gabbroic intrusion in an oceanic detachment fault. Ti-in-zircon thermometry yields temperatures around 820°C for the crystallization of the evolved melt. Reaction path modeling indicates that the alteration assemblage includes serpentine of the adjacent altered peridotites. Based on the model results, we propose that formation of chlorite occurred at higher temperatures than serpentinization, thus leading to strain localization around former plagiogranites during alteration. The detachment fault represents a major pathway for fluids through the oceanic crust, as evidenced by extremely low d18O of altered plagiogranite veins (+3.0-4.2 per mil) and adjacent serpentinites (+ 2.6-3.7 per mil). The uniform oxygen isotope data indicate that fluid flow in the detachment fault system affected veins and adjacent host serpentinites likewise.

(Table 2) Compressional and shear wave velocites from basalts and sediments of the Atlantic and Pacific Ocean

Compressional (Vp) and shear (Vs) wave velocities have been measured to 1.0 kbar for 14 cores of well-consolidated sedimentary rock from Atlantic and Pacific sites of the Deep Sea Drilling Project. The range of VP (2.05-5.38 km/sec at 0.5 kbar) shows significant overlap with the range of oceanic layer-2 seismic velocities determined by marine refraction surveys, suggesting that sedimentary rocks may, in some regions, constitute the upper portion of layer 2. Differing linear relationships between VP and Vs for basalts and sedimentary rocks, however, may provide a method of resolving layer-2 composition. This is illustra ted for a refraction survey site on the flank of the Mid-Atlantic Ridge where layer-2 velocities agree with basalt, and two sites on the Saya de Malha Bank in the Indian Ocean where layer-2 velocities appear to represent sedimentary rock.