Modelling of landslides propagation with SPH : effects of rheology and pore water pressure

El estudio desarrollado en este trabajo de tesis se centra en la modelización numérica de la fase de propagación de los deslizamientos rápidos de ladera a través del método sin malla Smoothed Particle Hydrodynamics (SPH). Este método tiene la gran ventaja de permitir el análisis de problemas de grandes deformaciones evitando operaciones costosas de remallado como en el caso de métodos numéricos con mallas tal como el método de los Elementos Finitos. En esta tesis, particular atención viene dada al rol que la reología y la presión de poros desempeñan durante estos eventos. El modelo matemático utilizado se basa en la formulación de Biot-Zienkiewicz v - pw, que representa el comportamiento, expresado en términos de velocidad del esqueleto sólido y presiones de poros, de la mezcla de partículas sólidas en un medio saturado. Las ecuaciones que gobiernan el problema son: • la ecuación de balance de masa de la fase del fluido intersticial, • la ecuación de balance de momento de la fase del fluido intersticial y de la mezcla, • la ecuación constitutiva y • una ecuación cinemática. Debido a sus propiedades geométricas, los deslizamientos de ladera se caracterizan por tener una profundidad muy pequeña frente a su longitud y a su anchura, y, consecuentemente, el modelo matemático mencionado anteriormente se puede simplificar integrando en profundidad las ecuaciones, pasando de un modelo 3D a 2D, el cual presenta una combinación excelente de precisión, sencillez y costes computacionales. El modelo propuesto en este trabajo se diferencia de los modelos integrados en profundidad existentes por incorporar un ulterior modelo capaz de proveer información sobre la presión del fluido intersticial a cada paso computacional de la propagación del deslizamiento. En una manera muy eficaz, la evolución de los perfiles de la presión de poros está numéricamente resuelta a través de un esquema explicito de Diferencias Finitas a cada nodo SPH. Este nuevo enfoque es capaz de tener en cuenta la variación de presión de poros debida a cambios de altura, de consolidación vertical o de cambios en las tensiones totales. Con respecto al comportamiento constitutivo, uno de los problemas principales al modelizar numéricamente deslizamientos rápidos de ladera está en la dificultad de simular con la misma ley constitutiva o reológica la transición de la fase de iniciación, donde el material se comporta como un sólido, a la fase de propagación donde el material se comporta como un fluido. En este trabajo de tesis, se propone un nuevo modelo reológico basado en el modelo viscoplástico de Perzyna, pensando a la viscoplasticidad como a la llave para poder simular tanto la fase de iniciación como la de propagación con el mismo modelo constitutivo. Con el fin de validar el modelo matemático y numérico se reproducen tanto ejemplos de referencia con solución analítica como experimentos de laboratorio. Finalmente, el modelo se aplica a casos reales, con especial atención al caso del deslizamiento de 1966 en Aberfan, mostrando como los resultados obtenidos simulan con éxito estos tipos de riesgos naturales. The study developed in this thesis focuses on the modelling of landslides propagation with the Smoothed Particle Hydrodynamics (SPH) meshless method which has the great advantage of allowing to deal with large deformation problems by avoiding expensive remeshing operations as happens for mesh methods such as, for example, the Finite Element Method. In this thesis, special attention is given to the role played by rheology and pore water pressure during these natural hazards. The mathematical framework used is based on the v - pw Biot-Zienkiewicz formulation, which represents the behaviour, formulated in terms of soil skeleton velocity and pore water pressure, of the mixture of solid particles and pore water in a saturated media. The governing equations are: • the mass balance equation for the pore water phase, • the momentum balance equation for the pore water phase and the mixture, • the constitutive equation and • a kinematic equation. Landslides, due to their shape and geometrical properties, have small depths in comparison with their length or width, therefore, the mathematical model aforementioned can then be simplified by depth integrating the equations, switching from a 3D to a 2D model, which presents an excellent combination of accuracy, computational costs and simplicity. The proposed model differs from previous depth integrated models by including a sub-model able to provide information on pore water pressure profiles at each computational step of the landslide's propagation. In an effective way, the evolution of the pore water pressure profiles is numerically solved through a set of 1D Finite Differences explicit scheme at each SPH node. This new approach is able to take into account the variation of the pore water pressure due to changes of height, vertical consolidation or changes of total stress. Concerning the constitutive behaviour, one of the main issues when modelling fast landslides is the difficulty to simulate with the same constitutive or rheological model the transition from the triggering phase, where the landslide behaves like a solid, to the propagation phase, where the landslide behaves in a fluid-like manner. In this work thesis, a new rheological model is proposed, based on the Perzyna viscoplastic model, thinking of viscoplasticity as the key to close the gap between the triggering and the propagation phase. In order to validate the mathematical model and the numerical approach, benchmarks and laboratory experiments are reproduced and compared to analytical solutions when possible. Finally, applications to real cases are studied, with particular attention paid to the Aberfan flowslide of 1966, showing how the mathematical model accurately and successfully simulate these kind of natural hazards.

Pulsatile blood flow interaction with arterial walls of aorta : autoregulation and impedance pressure boundary condition and its biomedical applications

Para las decisiones urgentes sobre intervenciones quirúrgicas en el sistema cardiovascular se necesitan simulaciones computacionales con resultados fiables y que consuman un tiempo de cálculo razonable. Durante años los investigadores han trabajado en diversos métodos numéricos de cálculo que resulten atractivos para los cirujanos. Estos métodos, precisos pero costosos desde el punto de vista del coste computacional, crean un desajuste entre la oferta de los ingenieros que realizan las simulaciones y los médicos que operan en el quirófano. Por otra parte, los métodos de cálculo más simplificados reducen el tiempo de cálculo pero pueden proporcionar resultados no realistas. El objetivo de esta tesis es combinar los conceptos de autorregulación e impedancia del sistema circulatorio, la interacción flujo sanguíneo-pared arterial y modelos geométricos idealizados tridimensionales de las arterias pero sin pérdida de realismo, con objeto de proponer una metodología de simulación que proporcione resultados correctos y completos, con tiempos de cálculo moderados. En las simulaciones numéricas, las condiciones de contorno basadas en historias de presión presentan inconvenientes por ser difícil conocerlas con detalle, y porque los resultados son muy sensibles ante pequeñas variaciones de dichas historias. La metodología propuesta se basa en los conceptos de autorregulación, para imponer la demanda de flujo aguas abajo del modelo en el ciclo cardiaco, y la impedancia, para representar el efecto que ejerce el flujo en el resto del sistema circulatorio sobre las arterias modeladas. De este modo las historias de presión en el contorno son resultados del cálculo, que se obtienen de manera iterativa. El método propuesto se aplica en una geometría idealizada del arco aórtico sin patologías y en otra geometría correspondiente a una disección Stanford de tipo A, considerando la interacción del flujo pulsátil con las paredes arteriales. El efecto de los tejidos circundantes también se incorpora en los modelos. También se hacen aplicaciones considerando la interacción en una geometría especifica de un paciente anciano que proviene de una tomografía computarizada. Finalmente se analiza una disección Stanford tipo B con tres modelos que incluyen la fenestración del saco. Clinicians demand fast and reliable numerical results of cardiovascular biomechanic simulations for their urgent pre-surgery decissions. Researchers during many years have work on different numerical methods in order to attract the clinicians' confidence to their colorful contours. Though precise but expensive and time-consuming methodologies create a gap between numerical biomechanics and hospital personnel. On the other hand, simulation simplifications with the aim of reduction in computational time may cause in production of unrealistic outcomes. The main objective of the current investigation is to combine ideas such as autoregulation, impedance, fluid-solid interaction and idealized geometries in order to propose a computationally cheap methodology without excessive or unrealistic simplifications. The pressure boundary conditions are critical and polemic in numerical simulations of cardiovascular system, in which a specific arterial site is of interest and the rest of the netwrok is neglected but represented by a boundary condition. The proposed methodology is a pressure boundary condition which takes advantage of numerical simplicity of application of an imposed pressure boundary condition on outlets, while it includes more sophisticated concepts such as autoregulation and impedance to gain more realistic results. Incorporation of autoregulation and impedance converts the pressure boundary conditions to an active and dynamic boundary conditions, receiving feedback from the results during the numerical calculations and comparing them with the physiological requirements. On the other hand, the impedance boundary condition defines the shapes of the pressure history curves applied at outlets. The applications of the proposed method are seen on idealized geometry of the healthy arotic arch as well as idealized Stanford type A dissection, considering the interaction of the arterial walls with the pulsatile blood flow. The effect of surrounding tissues is incorporated and studied in the models. The simulations continue with FSI analysis of a patient-specific CT scanned geometry of an old individual. Finally, inspiring of the statistic results of mortality rates in Stanford type B dissection, three models of fenestrated dissection sac is studied and discussed. Applying the developed boundary condition, an alternative hypothesis is proposed by the author with respect to the decrease in mortality rates in patients with fenestrations.

Role of Dopamine D2 Receptors in the Consolidation of Amphetamine-Cue Memory in Conditioned Activity in Rats

The neurotransmitter dopamine (DA) plays an essential role in reward-related incentive learning, whereby neutral stimuli gain the ability to elicit approach and other responses. In an incentive learning paradigm called conditioned activity, animals receive a stimulant drug in a specific environment over the course of several days. When then placed in that environment drug-free, they generally display a conditioned hyperactive response. Modulating DA transmission at different time points during the paradigm has been shown to disrupt or enhance conditioning effects. For instance, blocking DA D2 receptors before sessions generally impedes the acquisition of conditioned activity. To date, no studies have examined the role of D2 receptors in the consolidation phase of conditioned activity; this phase occurs immediately after acquisition and involves the stabilization of memories for long-term storage. To investigate this possible role, I trained Wistar rats (N = 108) in the conditioned activity paradigm produced by amphetamine (2.0 mg/kg, intraperitoneally) to examine the effects of the D2 antagonist haloperidol (doses 0.10, 0.25, 0.50, 0.75, 1.0, & 2.0 mg/kg, intraperitoneally) administered 5 min after conditioning sessions. Two positive control groups received haloperidol 1 h before conditioning sessions (doses 1.0 mg/kg and 2.0 mg/kg). The results revealed that post-session haloperidol at all doses tested did not disrupt the consolidation of conditioned activity, while pre-session haloperidol at 2.0 mg/kg prevented acquisition, with the 1.0 mg/kg group trending toward a block. Additionally, post-session haloperidol did not diminish activity during conditioning days, unlike pre-session haloperidol. One possible reason for these findings is that the consolidation phase may have begun earlier than when haloperidol was administered, since the conditioned activity paradigm uses longer learning sessions than those generally used in consolidation studies. Future studies may test if conditioned activity can be achieved with shorter sessions; if so, haloperidol would then be re-tested at an earlier time point. D2 receptor second messenger systems may also be investigated in consolidation. Since drug-related incentive stimuli can evoke cravings in those with drug addiction, a better understanding of the mechanisms of incentive learning may lead to the development of solutions for these individuals.

Stepwise consolidation of glacigenic sediments of the Prydz Bay

A series of samples from the five sites drilled across the continental shelf and upper slope in Prydz Bay during ODP Leg 119 were consolidation tested in an oedometer. Preconsolidation stresses increase downcore at Sites 739 and 742 in a stepwise manner, and the steps are interpreted to represent periods of increased action of grounded glaciers covering the entire shelf. By the use of theoretical ice sheet surface profiles giving the range of possible ice thicknesses, sediment loading and subsequent erosion seem to be the most important factor for increasing the overconsolidation ratios, and a total glacial erosion exceeding 1 km is possible. Four separate steps in consolidation, here termed "load events" have been identified. The lowermost load event, 1, is correlated to the onset of glaciations reaching the shelf edge and an early period of extensive glaciations, starting in early Oligocene or possibly earlier. Glacial activity related to the buildup of ice in West Antarctica in the late Miocene is tentatively correlated to load event 2. Event 3 is the trace of relatively extensive glacial erosion probably in the Pliocene, whereas the upper step in preconsolidation stress, load event 4, results from the last glaciation reaching the shelf edge, possibly during the late Weichselian. Correlations to other data related to Antarctic glacial history are, however, hampered by the poor age control of the cored diamictites. Consolidation tests may provide a tool for finding the position for hiatuses and unconformities formed subglacially and obscured by subglacial reworking.

(Table 1) Petrological description and consolidation properties of ODP Site 186-1150 and Hole 186-1151A

In order to determine the shear parameters of the forearc sedimentary strata drilled during Ocean Drilling Program Leg 186, West Pacific Seismic Network, Japan Trench, eight whole-round samples were selected from different depths in the drilled sections of Sites 1150 and 1151. Whereas Site 1150 lays above the seismically active part of the subduction zone, Site 1151 is situated in an aseismic zone. The aim of the triaxial tests was, apart from determination of the static stress strain behavior of the sediments, to test the hypothesis that the static stress strain parameter could differ for each sites. In order to simulate undrained deformation conditions according to the high clay mineral content of the strata, consolidated undrained shear tests were performed in a triaxial testing setup. Measurements of water content, grain density, organic content, and microtextural investigations under the scanning electron microscope (SEM) accompanied the compression experiments. After the saturation and consolidation stages were completed, failure occurred in the compression stage of the experiments at peak strengths of 280-7278 kPa. The stiffness moduli calculated for each sample from differential stress vs. strain curves show a linear relationship with depth and range between 181 and 5827 kPa. Under the SEM, the artificial fault planes of the tested specimen only show partial alignment of clay minerals because of the high content of microfossils.

(Table 1) Consolidation test for sediments of ODP Leg 130 sites

Consolidation tests were performed on 19 samples of calcareous ooze from the Ontong Java Plateau, obtained during Ocean Drilling Program Leg 130. Rebound curves from consolidation tests on Ontong Java Plateau samples yield porosity rebounds of 1%-4% for these sediments at equivalent depths up to 1200 mbsf. The exception is a radiolarian-rich sample that has 6% rebound. A rebound correction derived from the porosity rebound vs. depth data has been combined with a correction for pore-water expansion to correct the shipboard laboratory porosity data to in-situ values. Comparison of the laboratory porosity data corrected in this manner with the downhole log data shows good agreement.

Consolidation characteristics of outer Bengal Fan deposits (Table 1)

The state of consolidation of outer Bengal Fan deposits is closely related to their depositional history. Early Pleistocene sediments are underconsolidated, i.e., fabric strength is lower than under equilibrium with the present overburden stress, at depths greater than 60-80 mbsf. This may be due to rapid accumulation (>15-20 cm/k.y.) and overlying low-permeability (<10**-7 cm/s) deposits. The underlying Pliocene sediments are overconsolidated, i.e., fabric strength is higher than in equilibrium with the present overburden stress; at Sites 717 and 719 the sediments are slightly overconsolidated, whereas at Site 718 the overconsolidation is significant. Overconsolidation is explained by erosion that occurred during the early Pleistocene for which a stratigraphic gap was recorded in the drilled cores. The eroded section is estimated to be in the range of 10-40 m at Sites 717 and 719, respectively, and 130-150 m at Site 718. Below 250-300 mbsf the most sediments are normally consolidated. The amount of erosion seems to be related to block rotation and uplift due to intraplate deformation.

Who's really hypertensive? - Quality control issues in the assessment of blood pressure for randomized trials

The characterization of blood pressure in treatment trials assessing the benefits of blood pressure lowering regimens is a critical factor for the appropriate interpretation of study results. With numerous operators involved in the measurement of blood pressure in many thousands of patients being screened for entry into clinical trials, it is essential that operators follow pre-defined measurement protocols involving multiple measurements and standardized techniques. Blood pressure measurement protocols have been developed by international societies and emphasize the importance of appropriate choice of cuff size, identification of Korotkoff sounds, and digit preference. Training of operators and auditing of blood pressure measurement may assist in reducing the operator-related errors in measurement. This paper describes the quality control activities adopted for the screening stage of the 2nd Australian National Blood Pressure Study (ANBP2). ANBP2 is cardiovascular outcome trial of the treatment of hypertension in the elderly that was conducted entirely in general practices in Australia. A total of 54 288 subjects were screened; 3688 previously untreated subjects were identified as having blood pressure >140/90 mmHg at the initial screening visit, 898 (24%) were not eligible for study entry after two further visits due to the elevated reading not being sustained. For both systolic and diastolic blood pressure recording, observed digit preference fell within 7 percentage points of the expected frequency. Protocol adherence, in terms of the required minimum blood pressure difference between the last two successive recordings, was 99.8%. These data suggest that adherence to blood pressure recording protocols and elimination of digit preferences can be achieved through appropriate training programs and quality control activities in large multi-centre community-based trials in general practice. Repeated blood pressure measurement prior to initial diagnosis and study entry is essential to appropriately characterize hypertension in these elderly patients.

Gas and particle dynamics of a contoured shock tube for pre-clinical microparticle drug delivery

We investigate the gas-particle dynamics of a device designed for biological pre-clinical experiments. The device uses transonic/supersonic gas flow to accelerate microparticles such that they penetrate the outer skin layers. By using a shock tube coupled to a correctly expanded nozzle, a quasi-one-dimensional, quasi-steady flow (QSF) is produced to uniformly accelerate the microparticles. The system utilises a microparticle cassette (a diaphragm sealed container) that incorporates a jet mixing mechanism to stir the particles prior to diaphragm rupture. Pressure measurements reveal that a QSF exit period - suitable for uniformly accelerating microparticles - exists between 155 and 220 mus after diaphragm rupture. Immediately preceding the QSF period, a starting process secondary shock was shown to form with its (x,t) trajectory comparing well to theoretical estimates. To characterise the microparticle, flow particle image velocimetry experiments were conducted at the nozzle exit, using particle payloads with varying diameter (2.7-48 mu m), density (600-16,800 kg/m(3)) and mass (0.25-10 mg). The resultant microparticle velocities were temporally uniform. The experiments also show that the starting process does not significantly influence the microparticle nozzle exit velocities. The velocity distribution across the nozzle exit was also uniform for the majority of microparticle types tested. For payload masses typically used in pre-clinical drug and vaccine applications (

Predicting changes in measured introcular pressure with non-contact tonometry after LASIK

Purpose: Optometrists are becoming more integrally involved in the diagnosis of and care for glaucoma patients in the UK. The correlation of apparent change in non contact tonometry (NCT) IOP measurement and change in other ocular parameters such as refractive error, corneal curvature, corneal thickness and treatment zone size (data available to optometrists after LASIK) would facilitate care of these patients. Setting: A UK Laser Eye Clinic. Methods: This is a retrospective study study of 200 sequential eyes with myopia with or without astigmatism which underwent LASIK using a Hansatome and an Alcon LADARvision 4000 excimer laser. Refraction keratometry, pachymetry and NCT IOP mesurements were taken before treatmebnt and agian 3 months after treatment. The relationship between these variables anfd teh treatment zones were studied using stepwise multiple regression analysis. Results: There was a mean difference of 5.54mmHg comnparing pre and postoperative NCT IOP. IOP change correlates with refractive error change (P < 0.001), preoperative corneal thickness (P < 0.001) and treatment zone size (P = 0.047). Preoperative corneal thickness correlates with preoperative IOP (P < 0.001) and postoperative IOP (P < 0.001). Using these correlations, the measured difference in NCT IIOP can be predicted preoperatively or postoperatively using derived equations.Conclusion: There is a significant reduction in measured NCT IOP after LASIK. The amount of reduction can be calculated using data acquired by optometrists. This is helpful for opthalmologists and optometrists who co-manage glaucoma patients who have had LASIK or with glaucoma pateints who are consideraing having LASIK.

Working with the public sector:SME medical device innovation and pre-procurement

Healthcare providers are under ever increasing pressure to deliver more technologically advanced care without increasing costs. Innovation is essential (Darzi, 2008), and for this healthcare providers rely on innovation within commercial companies. SMEs have an important part to play in this sector (NHS Supply Chain Parliamentary Brief, 2013). Collaboration between SME suppliers and the NHS for innovation forms the focus of this paper. We examine the academic literature on interorganizational innovation including academic insights from the areas of forward commitment procurement (Environmental Innovation Advisory Group, 2003-2008), pre-commercial procurement (Bos & Corvers, 2007), innovation and SMEs. We then explore practice, first from a policy and business sector perspective. Second, we present evidence from fifteen cases of interorganizational innovation projects involving SMEs and UK healthcare providers. Our findings show much effort is being put into creating opportunities for more interorganizational innovation of medical devices. Working across organizational boundaries presents added complexity to the innovation environment and process, and the challenge of developing high-quality cross-boundary group interaction.

Maternal low-protein diet during mouse pre-implantation development induces vascular dysfunction and altered renin-angiotensin-system homeostasis in the offspring

Environmental perturbations during early mammalian development can affect aspects of offspring growth and cardiovascular health. We have demonstrated previously that maternal gestational dietary protein restriction in mice significantly elevated adult offspring systolic blood pressure. Therefore, the present study investigates the key mechanisms of blood pressure regulation in these mice. Following mating, female MF-1 mice were assigned to either a normal-protein diet (NPD; 18% casein) or an isocaloric low-protein diet throughout gestation (LPD; 9% casein), or fed the LPD exclusively during the pre-implantation period (3.5d) before returning to the NPD for the remainder of gestation (Emb-LPD). All offspring received standard chow. At 22 weeks, isolated mesenteric arteries from LPD and Emb-LPD males displayed significantly attenuated vasodilatation to isoprenaline (P=0.04 and P=0.025, respectively), when compared with NPD arteries. At 28 weeks, stereological analysis of glomerular number in female left kidneys revealed no significant difference between the groups. Real-time RT-PCR analysis of type 1a angiotensin II receptor, Na /K ATPase transporter subunits and glucocorticoid receptor expression in male and female left kidneys revealed no significant differences between the groups. LPD females displayed elevated serum angiotensin-converting enzyme (ACE) activity (P=0.044), whilst Emb-LPD males had elevated lung ACE activity (P=0.001), when compared with NPD offspring. These data demonstrate that elevated offspring systolic blood pressure following maternal gestational protein undernutrition is associated with impaired arterial vasodilatation in male offspring, elevated serum and lung ACE activity in female and male offspring, respectively, but kidney glomerular number in females and kidney gene expression in male and female offspring appear unaffected. © 2010 The Authors.

Performance evaluation of reverse osmosis (RO) pre-treatment technologies for in-land brackish water treatment

Integration of renewable energy with desalination technologies has emerged as an attractive solution to augment fresh water supply sustainably. Fouling and scaling are still considered as limiting factors in membrane desalination processes. For brackish water treatment, pre-treatment of reverse osmosis (RO) feed water is a key step in designing RO plants avoiding membrane fouling. This study aims to compare at pilot scale the rejection efficiency of RO membranes with multiple pre-treatment options at different water recoveries (30, 35, 40, 45 and 50%) and TDS concentrations (3500, 4000, and 4500mg/L). Synthetic brackish water was prepared and performance evaluation were carried out using brackish water reverse osmosis (BWRO) membranes (Filmtec LC-LE-4040 and Hydranautics CPA5-LD-4040) preceded by 5 and 1μm cartridge filters, 0.02μm ultra-filtration (UF) membrane, and forward osmosis (FO) membrane using 0.25M NaCl and MgCl2 as draw solutions (DS). It was revealed that FO membrane with 0.25M MgCl2 used as a draw solution (DS) and Ultra-filtration (UF) membrane followed by Filmtec membrane gave overall 98% rejection but UF facing high fouling potential due to high applied pressure. Use of 5 and 1μm cartridge filter prior to Filmtec membrane also showed effective results with 95% salt rejection.