Composition of the argillite sequence, continental slope of the Sea of Japan

The argillite sequence located at the base of the sedimentary cover on the continental slope of the Sea of Japan was studied by petrographic, palynological, and X-ray diffraction methods. Two spores-pollen complexes were distinguished in it: the Late Oligocene reflecting cooling and the Early Miocene corresponding to initiated warming. Data obtained indicate that the sequence is composed of terrigenous silty-clayey sediments that accumulated in shallow coastal-marine settings. The global sea-level rise at the Early-Middle Miocene transition, combined with regional tectonic processes, determined basin's deepening, owing to which the argillite sequence was overlain by a thick layer of Middle Miocene diatomaceous-clayey sediments. Due to tectonic movement along existing faults in the terminal Late Miocene, the argillite sequence occurring initially at depths of at least 400-500 m was locally exhumed to the basin bottom.

Lithology and pollen record of 4 profiles from the Bavaria, southern Germany

Palynological investigations of sediments from northern Bavaria (Rhön, Grabfeld, Lange Berge) reveal the Late Glacial and Postglacial history of the regional vegetation. The older sedirnents were found in the Rhön (Schwarzes Moor) and date back into the Bölling Period. At the end of that period pine spread into the Grabfeld. In both areas Lacher Tuff has been found. A radiocarbon date of 10,300 BP was found for the Late Glacial - Postglacial transition and one of 9300 BP for the Preboreal - Boreal transition. Hazel reached its highest values in the Rhön around 7,400 BP. During the Atlanticum a deciduous mixed oak forest covered the Rhön and Grabfeld regions. Beech dominated since the Subatlanticum. In the Lange Berge region, however, a mixed forest with Fagus, Picea, Pinus and Abies developed. In the Rhön first anthropogenic influence was found during the Latene Period. The boundary between zone IX and X has been dated at 820 A.D., and the start of extensive forest clearances at 1000 A. D. A culmination of landuse was found for the Medieval Period. At the end of that period however the Rhön was deserted. New forest clearances started around 1500 A.D., but were interrupted by the 'Thirty Years War'. Afterwards the Rhön got its present appearance.

Pollen analysis of soil samples from the A.D. 79 level. Boscoreale, Oplontis, and Pompeii

Fifty samples of Roman time soil preserved under the thick ash layer of the A.D.79 eruption of Mt Vesuvius were studied by pollen analysis: 33 samples from a former vineyard surrounding a Villa Rustica at Boscoreale (excavation site 40 x 50 m), 13 samples taken along the 60 m long swimming pool in the sculpture garden of the Villa of Poppaea at Oplontis, and four samples from the formal garden (12.4 x 17.5 m) of the House of the Gold Bracelet in Pompeii. To avoid contamination with modern pollen all samples were taken immediately after uncovering a new portion of the A.D. 79 soil. For comparison also samples of modern Italian soils were studied. Using standard methods for pollen preparation the pollen content of 15 of the archaeological samples proved to be too little to reach a pollen sum of more than 100 grains. The pollen spectra of these samples are not shown in the pollen tables. (Flotation with a sodium tungstate solution, Na2WO4, D = 2.05, following treatment with HCl and NaOH would probably have given a somewhat better result. This method was, however, not available as too expensive at that time.) Although the archaeological samples were taken a few meters apart their pollen values differ very much from one sample to the other. E.g., at Boscoreale (SW quarter). the pollen values of Pinus range from 1.5 to 54.5% resp. from 1 to 244 pine pollen grains per 1 gram of soil, the extremes even found under pine trees. Vitis pollen was present in 7 of the 11 vineyard samples from Boscoreale (NE quarter) only. Although a maximum of 21.7% is reached, the values of Vitis are mostly below 1.5%. Even the values of common weeds differ very much, not only at Boscoreale, but also at the other two sites. The pollen concentration values show similar variations: 3 to 3053 grains and spores were found in 1 g of soil. The mean value (290) is much less than the number of pollen grains, which would fall on 1 cm2 of soil surface during one year. In contrast, the pollen and spore concentrations of the recent soil samples, treated in exactly the same manner, range from 9313 to almost 80000 grains per 1 g of soil. Evidently most of the Roman time pollen has disappeared since its deposition, the reasons not being clear. Not even species which are known to have been cultivated in the garden of Oplontis, like Citrus and Nerium, plant species with easily distinguishable pollen grains, could be traced by pollen analysis. The loss of most of the pollen grains originally contained in the soil prohibits any detailed interpretation of the Pompeian pollen data. The pollen counts merely name plant species which grew in the region, but not necessarily on the excavated plots.

Results of n-alkane biomarker analyses and paleoclimate reconstruction

The current study presents quantitative reconstructions of tree cover, annual precipitation and mean July temperature derived from the pollen record from Lake Billyakh (65°17'N, 126°47'E, 340 m above sea level) spanning the last ca. 50 kyr. The reconstruction of tree cover suggests presence of woody plants through the entire analyzed time interval, although trees played only a minor role in the vegetation around Lake Billyakh prior to 14 kyr BP (<5%). This result corroborates low percentages of tree pollen and low scores of the cold deciduous forest biome in the PG1755 record from Lake Billyakh. The reconstructed values of the mean temperature of the warmest month ~8-10 °C do not support larch forest or woodland around Lake Billyakh during the coldest phase of the last glacial between ~32 and ~15 kyr BP. However, modern cases from northern Siberia, ca. 750 km north of Lake Billyakh, demonstrate that individual larch plants can grow within shrub and grass tundra landscape in very low mean July temperatures of about 8 °C. This makes plausible our hypothesis that the western and southern foreland of the Verkhoyansk Mountains could provide enough moist and warm microhabitats and allow individual larch specimens to survive climatic extremes of the last glacial. Reconstructed mean values of precipitation are about 270 mm/yr during the last glacial interval. This value is almost 100 mm higher than modern averages reported for the extreme-continental north-eastern Siberia east of Lake Billyakh, where larch-dominated cold deciduous forest grows at present. This suggests that last glacial environments around Lake Billyakh were never too dry for larch to grow and that the summer warmth was the main factor, which limited tree growth during the last glacial interval. The n-alkane analysis of the Siberian plants presented in this study demonstrates rather complex alkane distribution patterns, which challenge the interpretation of the fossil records. In particular, extremely low n-alkane concentrations in the leaves of local coniferous trees and shrubs suggest that their contribution to the litter and therefore to the fossil lake sediments might be not high enough for tracing the Quaternary history of the needleleaved taxa using the n-alkane biomarker method.

WCSPH viscosity diffusion processes in vortex flows

The aim of this paper is to clarify the role played by the most commonly used viscous terms in simulating viscous laminar flows using the weakly compressible approach in the context of smooth particle hydrodynamics (WCSPH). To achieve this, Takeda et al. (Prog. Theor. Phys. 1994; 92(5):939–960), Morris et al. (J. Comput. Phys. 1997; 136:214–226) and Monaghan–Cleary–Gingold's (Appl. Math. Model. 1998; 22(12):981–993; Monthly Notices of the Royal Astronomical Society 2005; 365:199–213) viscous terms will be analysed, discussing their origins, structures and conservation properties. Their performance will be monitored with canonical flows of which related viscosity phenomena are well understood, and in which boundary effects are not relevant. Following the validation process of three previously published examples, two vortex flows of engineering importance have been studied. First, an isolated Lamb–Oseen vortex evolution where viscous effects are dominant and second, a pair of co-rotating vortices in which viscous effects are combined with transport phenomena. The corresponding SPH solutions have been compared to finite-element numerical solutions. The SPH viscosity model's behaviour in modelling the viscosity related effects for these canonical flows is adequate

Theoretical analysis and numerical verification of the consistency of viscous smoothed-particle-hydrodynamics formulations in simulating free-surface flows

The theoretical formulation of the smoothed particle hydrodynamics (SPH) method deserves great care because of some inconsistencies occurring when considering free-surface inviscid flows. Actually, in SPH formulations one usually assumes that (i) surface integral terms on the boundary of the interpolation kernel support are neglected, (ii) free-surface conditions are implicitly verified. These assumptions are studied in detail in the present work for free-surface Newtonian viscous flow. The consistency of classical viscous weakly compressible SPH formulations is investigated. In particular, the principle of virtual work is used to study the verification of the free-surface boundary conditions in a weak sense. The latter can be related to the global energy dissipation induced by the viscous term formulations and their consistency. Numerical verification of this theoretical analysis is provided on three free-surface test cases including a standing wave, with the three viscous term formulations investigated.

Particle methods parallel implementations by GP-GPU strategies

This paper outlines the problems found in the parallelization of SPH (Smoothed Particle Hydrodynamics) algorithms using Graphics Processing Units. Different results of some parallel GPU implementations in terms of the speed-up and the scalability compared to the CPU sequential codes are shown. The most problematic stage in the GPU-SPH algorithms is the one responsible for locating neighboring particles and building the vectors where this information is stored, since these specific algorithms raise many dificulties for a data-level parallelization. Because of the fact that the neighbor location using linked lists does not show enough data-level parallelism, two new approaches have been pro- posed to minimize bank conflicts in the writing and subsequent reading of the neighbor lists. The first strategy proposes an efficient coordination between CPU-GPU, using GPU algorithms for those stages that allow a straight forward parallelization, and sequential CPU algorithms for those instructions that involve some kind of vector reduction. This coordination provides a relatively orderly reading of the neighbor lists in the interactions stage, achieving a speed-up factor of x47 in this stage. However, since the construction of the neighbor lists is quite expensive, it is achieved an overall speed-up of x41. The second strategy seeks to maximize the use of the GPU in the neighbor's location process by executing a specific vector sorting algorithm that allows some data-level parallelism. Al- though this strategy has succeeded in improving the speed-up on the stage of neighboring location, the global speed-up on the interactions stage falls, due to inefficient reading of the neighbor vectors. Some changes to these strategies are proposed, aimed at maximizing the computational load of the GPU and using the GPU texture-units, in order to reach the maximum speed-up for such codes. Different practical applications have been added to the mentioned GPU codes. First, the classical dam-break problem is studied. Second, the wave impact of the sloshing fluid contained in LNG vessel tanks is also simulated as a practical example of particle methods

Benefits of using a Wendland Kernel for free-surface flows

The aim of this paper Is lo discuss the influence of the selection of the interpolation kernel in the accuracy of the modeling of the internal viscous dissipation in Tree surface Hows, Simulations corresponding to a standing wave* for which an analytic solution available, are presented. Wendland and renormalized Gaussian kernels are considered. The differences in the flow pattern* and Internal dissipation mechanisms are documented for a range of Reynolds numbers. It is shown that the simulations with Wendland kernels replicate the dissipation mechanisms more accurately than those with a renormalized Gaussian kernel. Although some explanations are hinted we have Tailed to clarify which the core structural reasons for Mich differences are*

Calibración de Modelos Constitutivos de Hormigón

En este artículo se describen diversos modelos constitutivos del comportamiento dilatante del hormigón. Se describe una metodología para caracterizar el comportamiento mecánico del hormigón. En particular se calibra el modelo de plasticidad y daño para hormigón de Abaqus/Standard y Abaqus/Explicit a partir de ensayos uniaxiales y tria xiales en probetas cilíndricas. El modelo se enriquece con una subrutina de usuario para capturar mejor la ductilidad del material sometido a moderadas presiones de confinamiento. Además se define una variable de estado para controlar la eliminación de elementos. Finalmente se muestran algunos ejemplos de aplicación de impactos de proyectiles deformables y rígidos sobre muros de hormigón armado, incluyendo técnicas especiales como la conversión de elementos sólidos a partículas (SPH).

Modelización de la rotura y fluidificación en geomateriales y geoestructuras

La rotura de las geoestructuras puede ser causada por cambios en las tensiones efectivas debidos bien por cargas externas (terremotos, por ejemplo), bien por variación de las presiones intersticiales (lluvia), o cambios en la geometría (erosión), así como por una disminución de las propiedades resistentes de los materiales (meteorización, ataque químico, etc). El caso particular de los deslizamientos es interesante, existiendo diversas clasificaciones que tienen en cuenta su forma, velocidad de propagación, etc. Dos de estos casos son los deslizamientos propiamente dichos y los flujos. En el primer caso, la deformación se concentra en zonas de pequeño espesor, que se idealiza como una superficie (superficie de rotura). La cinemática de esta rotura se puede considerar como el movimiento relativo de dos masas cuyas deformaciones no son grandes. Este mecanismo está usualmente asociado a materiales sobreconsolidados que presentan reblandecimiento. Los flujos se producen generalmente en materiales de baja densidad y estructura metaestable, con tendencia a compactar, de forma que se generan presiones intersticiales que aumentan el ángulo de rozamiento movilizado, pudiéndose llegar en algunos casos a la licuefacción. Este mecanismo de rotura se conoce como rotura difusa, y no ha sido tan estudiado como el de localización a pesar de se trata frecuentemente de roturas de tipo catastrófico. De hecho, el suelo pasa de un estado sólido a un estado fluidificado, con una gran movilidad. Es importante para el ingeniero predecir tanto el comportamiento de las geoestructuras bajo las cargas de cálculo como las condiciones en las que se producirá la rotura. De esta manera, en algunos casos, se podrán reforzar las zonas más débiles aumentando así su seguridad. En otros casos, no se podrá realizar este refuerzo (grandes deslizamientos como avalanchas, lahares, etc), pero sí se podrán conocer las consecuencias de la rotura: velocidad de propagación, alcance, espesores, etc. La modelización de estos problemas es compleja, ya que aparecen dificultades en los modelos matemáticos, constitutivos o reológicos y numéricos. Dado que en los geomateriales aparece una interacción fuerte entre el esqueleto sólido y los fluidos intersticiales, esto debe ser tenido en cuenta en los modelos matemáticos. En este trabajo se describirán, pues, el desarrollo y aplicación de técnicas avanzadas de modelización; matemática, constitutiva/reológica y numérica. Se dedicará especial atención a los problemas de transición entre suelos y suelos fluidificados, que hoy en día se estudian en una gran mayoría de los casos con modelos diferentes. Así por ejemplo, se emplean modelos constitutivos para el comportamiento previo a la rotura, y reológicos para los materiales fluidificados. En lo que respecta a los modelos matemáticos, existen formulaciones nuevas en velocidades (o desplazamientos), tensiones, y presiones de aire y agua intersticial, de los que se pueden obtener modelos simplificados integrados en profundidad para deslizamientos rápidos. Respecto de los modelos constitutivos, es interesante la teoría de la Plasticidad Generalizada (modelo básico de Pastor-Zienkiewicz y extensiones a suelos no saturados). Se estudiará la extensión de estos modelos elastoplásticos a la viscoplasticidad (Perzyna y consistente), explorando la posibilidad de emplearlos tanto antes como después de la rotura. Finalmente, en lo que a modelos numéricos se refiere, se describirá la implementación de los modelos matemáticos y constitutivos desarrollados en (i) modelos clásicos de elementos finitos, como el GeHoMadrid desarrollado en los grupo investigador M2i al que pertenece el autor de este trabajo, (ii) Métodos de tipo Taylor Galerkin, y (iii) métodos sin malla como el SPH y el Material Point Model. Estos modelos se aplicarán, principalmente a (i) Licuefacción de estructuras cimentadas en el fondo marino (ii) presas de residuos mineros (iii) deslizamientos rápidos de laderas.

Inception of debris avalanches: remarks on geomechanical modelling

Debris avalanches are complex phenomena due to the variety of mechanisms that control the failure stage and the avalanche formation. Regarding these issues, in the literature, either field evidence or qualitative interpretations can be found while few experimental laboratory tests and rare examples of geomechanical modelling are available for technical and/or scientific purposes. As a contribution to the topic, the paper firstly highlights as the problem can be analysed referring to a unique mathematical framework from which different modelling approaches can be derived based on limit equilibrium method (LEM), finite element method (FEM), or smooth particle hydrodynamics (SPH). Potentialities and limitations of these approaches are then tested for a large study area where huge debris avalanches affected shallow deposits of pyroclastic soils (Sarno-Quindici, Southern Italy). The numerical results show that LEM as well as uncoupled and coupled stress–strain FEM analyses are able to individuate the major triggering mechanisms. On the other hand, coupled SPH analyses outline the relevance of erosion phenomena, which can modify the kinematic features of debris avalanches in their source areas, i.e. velocity, propagation patterns and later spreading of the unstable mass. As a whole, the obtained results encourage the application of the introduced approaches to further analyse real cases in order to enhance the current capability to forecast the inception of these dangerous phenomena.

Dynamic interactions between solids and viscous liquids with free surface = Interacciones dinámicas entre sólidos y líquidos viscosos con superficie libre

En esta tesis se investiga la interacción entre un fluido viscoso y un cuerpo sólido en presencia de una superficie libre. El problema se expresa teóricamente poniendo especial atención a los aspectos de conservación de energía y de la interacción del fluido con el cuerpo. El problema se considera 2D y monofásico, y un desarrollo matemático permite una descomposición de los términos disipativos en términos relacionados con la superficie libre y términos relacionados con la enstrofía. El modelo numérico utilizado en la tesis se basa en el método sin malla Smoothed Particle Hydrodynamics (SPH). De manera análoga a lo que se hace a nivel continuo, las propiedades de conservación se estudian en la tesis con el sistema discreto de partículas. Se tratan también las condiciones de contorno de un cuerpo que se mueve en un flujo viscoso, implementadas con el método ghost-fluid. Se ha desarrollado un algoritmo explícito de interacción fluido / cuerpo. Se han documentado algunos casos de modo detallado con el objetivo de comprobar la capacidad del modelo para reproducir correctamente la disipación de energía y el movimiento del cuerpo. En particular se ha investigado la atenuación de una onda estacionaria, comparando la simulación numérica con predicciones teóricas. Se han realizado otras pruebas para monitorizar la disipación de energía para flujos más violentos que implican la fragmentación de la superficie libre. La cantidad de energía disipada con los diferentes términos se ha evaluado en los casos estudiados con el modelo numérico. Se han realizado otras pruebas numéricas para verificar la técnica de modelización de la interacción fluido / cuerpo, concretamente las fuerzas ejercidas por las olas en cuerpos con formas simples, y el equilibrio de un cuerpo flotante con una forma compleja. Una vez que el modelo numérico ha sido validado, se han realizado simulaciones numéricas para obtener una comprensión más completa de la física implicada en casos (casi) realistas sobre los había aspectos que no se conocían suficientemente. En primer lugar se ha estudiado el el flujo alrededor de un cilindro bajo la superficie libre. El estudio se ha realizado con un número de Reynolds moderado, para un rango de inmersiones del cilindro y números de Froude. La solución numérica permite una investigación de los patrones complejos que se producen. La estela del cilindro interactúa con la superficie libre. Se han identificado algunos inestabilidades características. El segundo estudio se ha realizado sobre el problema de sloshing, tanto experimentalmente como numéricamente. El análisis se restringe a aguas poco profundas y con oscilación horizontal, pero se ha estudiado un gran número de condiciones, lo que lleva a una comprensión bastante completa de los sistemas de onda involucradas. La última parte de la tesis trata también sobre un problema de sloshing pero esta vez el tanque está oscilando con rotación y hay acoplamiento con un sistema mecánico. El sistema se llama pendulum-TLD (Tuned Liquid Damper - con líquido amortiguador). Este tipo de sistema se utiliza normalmente para la amortiguación de las estructuras civiles. El análisis se ha realizado analíticamente, numéricamente y experimentalmente utilizando líquidos con viscosidades diferentes, centrándose en características no lineales y mecanismos de disipación. ABSTRA C T The subject of the present thesis is the interaction between a viscous fluid and a solid body in the presence of a free surface. The problem is expressed first theoretically with a particular focus on the energy conservation and the fluid-body interaction. The problem is considered 2D and monophasic, and some mathematical development allows for a decomposition of the energy dissipation into terms related to the Free Surface and others related to the enstrophy. The numerical model used on the thesis is based on Smoothed Particle Hydrodynamics (SPH): a computational method that works by dividing the fluid into particles. Analogously to what is done at continuum level, the conservation properties are studied on the discrete system of particles. Additionally the boundary conditions for a moving body in a viscous flow are treated and discussed using the ghost-fluid method. An explicit algorithm for handling fluid-body coupling is also developed. Following these theoretical developments on the numerical model, some test cases are devised in order to test the ability of the model to correctly reproduce the energy dissipation and the motion of the body. The attenuation of a standing wave is used to compare what is numerically simulated to what is theoretically predicted. Further tests are done in order to monitor the energy dissipation in case of more violent flows involving the fragmentation of the free-surface. The amount of energy dissipated with the different terms is assessed with the numerical model. Other numerical tests are performed in order to test the fluid/body interaction method: forces exerted by waves on simple shapes, and equilibrium of a floating body with a complex shape. Once the numerical model has been validated, numerical tests are performed in order to get a more complete understanding of the physics involved in (almost) realistic cases. First a study is performed on the flow passing a cylinder under the free surface. The study is performed at moderate Reynolds numbers, for various cylinder submergences, and various Froude numbers. The capacity of the numerical solver allows for an investigation of the complex patterns which occur. The wake from the cylinder interacts with the free surface, and some characteristical flow mechanisms are identified. The second study is done on the sloshing problem, both experimentally and numerically. The analysis is restrained to shallow water and horizontal excitation, but a large number of conditions are studied, leading to quite a complete understanding of the wave systems involved. The last part of the thesis still involves a sloshing problem but this time the tank is rolling and there is coupling with a mechanical system. The system is named pendulum-TLD (Tuned Liquid Damper). This kind of system is normally used for damping of civil structures. The analysis is then performed analytically, numerically and experimentally for using liquids with different viscosities, focusing on non-linear features and dissipation mechanisms.

Depth Averaged Models for Fast Landslide Propagation: Mathematical, Rheological and Numerical Aspects

This paper presents an overview of depth averaged modelling of fast catastrophic landslides where coupling of solid skeleton and pore fluid (air and water) is important. The first goal is to show how Biot-Zienkiewicz models can be applied to develop depth integrated, coupled models. The second objective of the paper is to consider a link which can be established between rheological and constitutive models. Perzyna´s viscoplasticity can be considered a general framework within which rheological models such as Bingham and cohesive frictional fluids can be derived. Among the several alternative numerical models, we will focus here on SPH which has not been widely applied by engineers to model landslide propagation. We propose an improvement, based on combining Finite Difference meshes associated to SPH nodes to describe pore pressure evolution inside the landslide mass. We devote a Section to analyze the performance of the models, considering three sets of tests and examples which allows to assess the model performance and limitations: (i) Problems having an analytical solution, (ii) Small scale laboratory tests, and (iii) Real cases for which we have had access to reliable information

Smooth particle hydrodynamics study of surface defect machining for diamond turning of silicon

Acknowledgments The authors would like to thank EPSRC (EP/ K018345/1) and Royal Society-NSFC International Exchange Scheme for providing financial support to this research.