GIS-based debris flow source and runout susceptibility assessment from DEM data - a case study in NW Nicaragua

In October 1998, Hurricane Mitch triggered numerous landslides (mainly debris flows) in Honduras and Nicaragua, resulting in a high death toll and in considerable damage to property. The potential application of relatively simple and affordable spatial prediction models for landslide hazard mapping in developing countries was studied. Our attention was focused on a region in NW Nicaragua, one of the most severely hit places during the Mitch event. A landslide map was obtained at 1:10 000 scale in a Geographic Information System (GIS) environment from the interpretation of aerial photographs and detailed field work. In this map the terrain failure zones were distinguished from the areas within the reach of the mobilized materials. A Digital Elevation Model (DEM) with 20 m×20 m of pixel size was also employed in the study area. A comparative analysis of the terrain failures caused by Hurricane Mitch and a selection of 4 terrain factors extracted from the DEM which, contributed to the terrain instability, was carried out. Land propensity to failure was determined with the aid of a bivariate analysis and GIS tools in a terrain failure susceptibility map. In order to estimate the areas that could be affected by the path or deposition of the mobilized materials, we considered the fact that under intense rainfall events debris flows tend to travel long distances following the maximum slope and merging with the drainage network. Using the TauDEM extension for ArcGIS software we generated automatically flow lines following the maximum slope in the DEM starting from the areas prone to failure in the terrain failure susceptibility map. The areas crossed by the flow lines from each terrain failure susceptibility class correspond to the runout susceptibility classes represented in a runout susceptibility map. The study of terrain failure and runout susceptibility enabled us to obtain a spatial prediction for landslides, which could contribute to landslide risk mitigation.

Evolution of a polymineralic mantle shear zone and the role of second phases in the localization of deformation

The influence of second phases (e.g., pyroxenes) on olivine grain size was studied by quantitative microfabric analyses of samples of the Hilti massif mantle shear zone (Semail ophiolite, Oman). The microstructures range from porphyroclastic tectonites to ultramylonites, from outside to the center of the shear zone. Starting at conditions of ridge-related flow, they formed under continuous cooling leading to progressive strain localization. The dependence of the average olivine grain size on the second-phase content can be split into a second-phase controlled and a dynamic recrystallization-controlled field. In the former, the olivine grain size is related to the ratio between the second-phase grain size and volume fraction (Zener parameter). In the latter, dynamic recrystallization manifested by a balance between grain growth and grain size reduction processes yields a stable olivine grain size. In both fields the average olivine and second-phase grain size decreases with decreasing temperature. Combining the microstructural information with deformation mechanism maps suggests that the porphyroclastic tectonites (similar to 1100 degrees C) and mylonites (similar to 800 degrees C) formed under the predominance of dislocation creep. Since olivine-rich layers are intercalated with layer parallel, polymineralic bands in the mylonites, nearly equiviscous conditions can be assumed. In the ultramylonites, diffusion creep represents the major deformation mechanism in the polymineralic layers. It is this switch in deformation mechanism from dislocation creep to diffusion creep that forces strain to localize in the fine-grained polymineralic domains at low temperatures (<similar to 700 degrees C), underlining the role of the second phases on strain localization in cooling mantle rocks.

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.

Buprenorphine and norbuprenorphine quantification in human plasma by simple protein precipitation and ultra-high performance liquid chromatography tandem mass spectrometry.

A highly sensitive ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) method was developed for the quantification of buprenorphine and its major metabolite norbuprenorphine in human plasma. In order to speed up the process and decrease costs, sample preparation was performed by simple protein precipitation with acetonitrile. To the best of our knowledge, this is the first application of this extraction technique for the quantification of buprenorphine in plasma. Matrix effects were strongly reduced and selectivity increased by using an efficient chromatographic separation on a sub-2μm column (Acquity UPLC BEH C18 1.7μm, 2.1×50mm) in 5min with a gradient of ammonium formate 20mM pH 3.05 and acetonitrile as mobile phase at a flow rate of 0.4ml/min. Detection was made using a tandem quadrupole mass spectrometer operating in positive electrospray ionization mode, using multiple reaction monitoring. The procedure was fully validated according to the latest Food and Drug Administration guidelines and the Société Française des Sciences et Techniques Pharmaceutiques. Very good results were obtained by using a stable isotope-labeled internal standard for each analyte, to compensate for the variability due to the extraction and ionization steps. The method was very sensitive with lower limits of quantification of 0.1ng/ml for buprenorphine and 0.25ng/ml for norbuprenorphine. The upper limit of quantification was 250ng/ml for both drugs. Trueness (98.4-113.7%), repeatability (1.9-7.7%), intermediate precision (2.6-7.9%) and internal standard-normalized matrix effects (94-101%) were in accordance with international recommendations. The procedure was successfully used to quantify plasma samples from patients included in a clinical pharmacogenetic study and can be transferred for routine therapeutic drug monitoring in clinical laboratories without further development.

Dynamic drag modeling of submerged aquatic vegetation canopy flows

Vegetation has a profound effect on flow and sediment transport processes in natural rivers, by increasing both skin friction and form drag. The increase in drag introduces a drag discontinuity between the in-canopy flow and the flow above, which leads to the development of an inflection point in the velocity profile, resembling a free shear layer. Therefore, drag acts as the primary driver for the entire canopy system. Most current numerical hydraulic models which incorporate vegetation rely either on simple, static plant forms, or canopy-scaled drag terms. However, it is suggested that these are insufficient as vegetation canopies represent complex, dynamic, porous blockages within the flow, which are subject to spatially and temporally dynamic drag forces. Here we present a dynamic drag methodology within a CFD framework. Preliminary results for a benchmark cylinder case highlight the accuracy of the method, and suggest its applicability to more complex cases.

Extending the canopy flow model for natural, highly flexible macrophyte canopies

Flow structures above vegetation canopies have received much attention within terrestrial and aquatic literature. This research has led to a good process understanding of mean and turbulent canopy flow structure. However, much of this research has focused on rigid or semi-rigid vegetation with relatively simple morphology. Aquatic macrophytes differ from this form, exhibiting more complex morphologies, predominantly horizontal posture in the flow and a different force balance. While some recent studies have investigated such canopies, there is still the need to examine the relevance and applicability of general canopy layer theory to these types of vegetation. Here, we report on a range of numerical experiments, using both semi-rigid and highly flexible canopies. The results for the semi-rigid canopies support existing canopy layer theory. However, for the highly flexible vegetation, the flow pattern is much more complex and suggests that a new canopy model may be required.

Analysis of flow dynamics in right ventricular outflow tract.

BACKGROUND: The mechanism behind early graft failure after right ventricular outflow tract (RVOT) reconstruction is not fully understood. Our aim was to establish a three-dimensional computational fluid dynamics (CFD) model of RVOT to investigate the hemodynamic conditions that may trigger the development of intimal hyperplasia and arteriosclerosis. METHODS: Pressure, flow, and diameter at the RVOT, pulmonary artery (PA), bifurcation of the PA, and left and right PAs were measured in 10 normal pigs with a mean weight of 24.8 ± 0.78 kg. Data obtained from the experimental scenario were used for CFD simulation of pressure, flow, and shear stress profile from the RVOT to the left and right PAs. RESULTS: Using experimental data, a CFD model was obtained for 2.0 and 2.5-L/min pulsatile inflow profiles. In both velocity profiles, time and space averaged in the low-shear stress profile range from 0-6.0 Pa at the pulmonary trunk, its bifurcation, and at the openings of both PAs. These low-shear stress areas were accompanied to high-pressure regions 14.0-20.0 mm Hg (1866.2-2666 Pa). Flow analysis revealed a turbulent flow at the PA bifurcation and ostia of both PAs. CONCLUSIONS: Identified local low-shear stress, high pressure, and turbulent flow correspond to a well-defined trigger pattern for the development of intimal hyperplasia and arteriosclerosis. As such, this real-time three-dimensional CFD model may in the future serve as a tool for the planning of RVOT reconstruction, its analysis, and prediction of outcome.

A simple model to describe the thixotropic behavior of paints.

We propose a simple rheological model to describe the thixotropic behavior of paints, since the classical hysteresis area, which is usually used, is not enough to evaluate thixotropy. The model is based on the assumption that viscosity is a direct measure of the structural level of the paint. The model depends on two equations: the Cross-Carreau equation to describe the equilibrium viscosity and a second order kinetic equation to express the time dependence of viscosity. Two characteristic thixotropic times are differentiated: one for the net structure breakdown, which is defined as a power law function of shear rate, and an other for the net structure buildup, which is not dependent on the shear rate. The knowledge of both kinetic processes can be used to improve the quality and applicability of paints. Five representative commercial protective marine paints are tested. They are based on chlorinated rubber, acrylic, alkyd, vinyl, and epoxy resins. The temperature dependence of the rheological behavior is also studied with the temperature ranging from 5 ºC to 35 ºC. It is found that the paints exhibit both shear thinning and thixotropic behavior. The model fits satisfactorily the thixotropy of the studied paints. It is also able to predict the thixotropy dependence on temperature. Both viscosity and the degree of thixotropy increase as the temperature decreases.

Calculations of Boiling Two-Phase Flow Using a Porous Media Model

Boiling two-phase flow and the equations governing the motion of fluid in two-phase flows are discussed in this thesis. Disposition of the governing equations in three-dimensional complex geometries is considered from the perspective of the porous medium concept. The equations governing motion in two-phase flows were formulated, discretized and implemented in a subroutine for pressure-velocity solution utilizing the SIMPLE algorithm modified for two-phase flow. The subroutine was included in PORFLO, which is a three-dimensional 5-equation porous media model developed at VTT by Jaakko Miettinen. The development of two-phase flow and the resulting void fraction distribution was predicted in a geometry resembling a section of BWR fuel bundle in a couple of test cases using PORFLO.

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.

Bernoulli correction to viscous losses: Radial flow between two parallel discs

For a massless fluid (density = 0), the steady flow along a duct is governed exclusively by viscous losses. In this paper, we show that the velocity profile obtained in this limit can be used to calculate the pressure drop up to the first order in density. This method has been applied to the particular case of a duct, defined by two plane-parallel discs. For this case, the first-order approximation results in a simple analytical solution which has been favorably checked against numerical simulations. Finally, an experiment has been carried out with water flowing between the discs. The experimental results show good agreement with the approximate solution

Identification of new transformation products during enzymatic treatment of tetracycline and erythromycin antibiotics at laboratory scale by an on-line turbulent flow liquid-chromatography coupled to a high resolution mass spectrometer LTQ-Orbitrap

This work describes the formation of transformation products (TPs) by the enzymatic degradation at laboratory scale of two highly consumed antibiotics: tetracycline (Tc) and erythromycin (ERY). The analysis of the samples was carried out by a fast and simple method based on the novel configuration of the on-line turbulent flow system coupled to a hybrid linear ion trap – high resolution mass spectrometer. The method was optimized and validated for the complete analysis of ERY, Tc and their transformation products within 10 min without any other sample manipulation. Furthermore, the applicability of the on-line procedure was evaluated for 25 additional antibiotics, covering a wide range of chemical classes in different environmental waters with satisfactory quality parameters. Degradation rates obtained for Tc by laccase enzyme and ERY by EreB esterase enzyme without the presence of mediators were ∼78% and ∼50%, respectively. Concerning the identification of TPs, three suspected compounds for Tc and five of ERY have been proposed. In the case of Tc, the tentative molecular formulas with errors mass within 2 ppm have been based on the hypothesis of dehydroxylation, (bi)demethylation and oxidation of the rings A and C as major reactions. In contrast, the major TP detected for ERY has been identified as the “dehydration ERY-A”, with the same molecular formula of its parent compound. In addition, the evaluation of the antibiotic activity of the samples along the enzymatic treatments showed a decrease around 100% in both cases

Process intensification: From optimised flow patterns to microprocess technology

This thesis is focused on process intensification. Several significant problems and applications of this theme are covered. Process intensification is nowadays one of the most popular trends in chemical engineering and attempts have been made to develop a general, systematic methodology for intensification. This seems, however, to be very difficult, because intensified processes are often based on creativity and novel ideas. Monolith reactors and microreactors are successful examples of process intensification. They are usually multichannel devices in which a proper feed technique is important for creating even fluid distribution into the channels. Two different feed techniques were tested for monoliths. In the first technique a shower method was implemented by means of perforated plates. The second technique was a dispersion method using static mixers. Both techniques offered stable operation and uniform fluid distribution. The dispersion method enabled a wider operational range in terms of liquid superficial velocity. Using dispersion method, a volumetric gas-liquid mass transfer coefficient of 2 s-1 was reached. Flow patterns play a significant role in terms of the mixing performance of micromixers. Although the geometry of a T-mixer is simple, channel configurations and dimensions had a clear effect on mixing efficiency. The flow in the microchannel was laminar, but the formation of vortices promoted mixing in micro T-mixers. The generation of vortices was dependent on the channel dimensions, configurations and flow rate. Microreactors offer a high ratio of surface area to volume. Surface forces and interactions between fluids and surfaces are, therefore, often dominant factors. In certain cases, the interactions can be effectively utilised. Different wetting properties of solid materials (PTFE and stainless steel) were applied in the separation of immiscible liquid phases. A micro-scale plate coalescer with hydrophilic and hydrophobic surfaces was used for the continuous separation of organic and aqueous phases. Complete phase separation occurred in less than 20 seconds, whereas the separation time by settling exceeded 30 min. Fluid flows can be also intensified in suitable conditions. By adding certain additives into turbulent fluid flow, it was possible to reduce friction (drag) by 40 %. Drag reduction decreases frictional pressure drop in pipelines which leads to remarkable energy savings and decreases the size or number of pumping facilities required, e.g., in oil transport pipes. Process intensification enables operation often under more optimal conditions. The consequent cost savings from reduced use of raw materials and reduced waste lead to greater economic benefits in processing.

Flow analysis by using solenoid valves for As(III) determination in natural waters by an on-line separation and pre-concentration system coupled to a tungsten coil atomizer

A flow system coupled to a tungsten coil atomizer in an atomic absorption spectrometer (TCA-AAS) was developed for As(III) determination in waters, by extraction with sodium diethyldithiocarbamate (NaDDTC) as complexing agent, and by sorption of the As(III)-DDTC complex in a micro-column filled with 5 mg C18 reversed phase (10 µL dry sorbent), followed by elution with ethanol. A complete pre-concentration/elution cycle took 208 s, with 30 s sample load time (1.7 mL) and 4 s elution time (71 µL). The interface and software for the synchronous control of two peristaltic pumps (RUN/ STOP), an autosampler arm, seven solenoid valves, one injection valve, the electrothermal atomizer and the spectrometer Read function were constructed. The system was characterized and validated by analytical recovery studies performed both in synthetic solutions and in natural waters. Using a 30 s pre-concentration period, the working curve was linear between 0.25 and 6.0 µg L-1 (r = 0.9976), the retention efficiency was 94±1% (6.0 µg L-1), and the pre-concentration coefficient was 28.9. The characteristic mass was 58 pg, the mean repeatability (expressed as the variation coefficient) was 3.4% (n=5), the detection limit was 0.058 µg L-1 (4.1 pg in 71 µL of eluate injected into the coil), and the mean analytical recovery in natural waters was 92.6 ± 9.5 % (n=15). The procedure is simple, economic, less prone to sample loss and contamination and the useful lifetime of the micro-column was between 200-300 pre-concentration cycles.

Pulsatile Blood Flow Simulations in Computed Tomography(CT) Scan-Based and Idealized Geometries of Human Aorta

Blood flow in human aorta is an unsteady and complex phenomenon. The complex patterns are related to the geometrical features like curvature, bends, and branching and pulsatile nature of flow from left ventricle of heart. The aim of this work was to understand the effect of aorta geometry on the flow dynamics. To achieve this, 3D realistic and idealized models of descending aorta were reconstructed from Computed Tomography (CT) images of a female patient. The geometries were reconstructed using medical image processing code. The blood flow in aorta was assumed to be laminar and incompressible and the blood was assumed to be Newtonian fluid. A time dependent pulsatile and parabolic boundary condition was deployed at inlet. Steady and unsteady blood flow simulations were performed in real and idealized geometries of descending aorta using a Finite Volume Method (FVM) code. Analysis of Wall Shear Stress (WSS) distribution, pressure distribution, and axial velocity profiles were carried out in both geometries at steady and unsteady state conditions. The results obtained in thesis work reveal that the idealization of geometry underestimates the values of WSS especially near the region with sudden change of diameter. However, the resultant pressure and velocity in idealized geometry are close to those in real geometry