Strain rate of the South American lithospheric plate by SIRGAS-CON geodetic observations

A multiyear solution of the SIRGAS-CON network was used to estimate the strain rates of the earth surface from the changing directions of the velocity vectors of 140 geodetic points located in the South American plate. The strain rate was determined by the finite element method using Delaunay triangulation points that formed sub-networks; each sub-network was considered a solid and homogeneous body. The results showed that strain rates vary along the South American plate and are more significant on the western portion of the plate, as expected, since this region is close to the subduction zone of the Nazca plate beneath the South American plate. After using Euler vectors to infer Nazca plate movement and to orient the velocity vectors of the South American plate, it was possible to estimate the convergence and accommodation rates of the Nazca and South American plates, respectively. Strain rate estimates permitted determination of predominant contraction and/or extension regions and to establish that contraction regions coincide with locations with most of the high magnitude seismic events. Some areas with extension and contraction strains were found to the east within the stable South American plate, which may result from different stresses associated with different geological characteristics. These results suggest that major movements detected on the surface near the Nazca plate occur in regions with more heterogeneous geological structures and multiple rupture events. Most seismic events in the South American plate are concentrated in areas with predominant contraction strain rates oriented northeast-southwest; significant amounts of elastic strain can be accumulated on geological structures away from the plate boundary faults; and, behavior of contractions and extensions is similar to what has been found in seismological studies. © 2013 Elsevier Ltd.

An analysis of some Eulerian methods for one dimensional, inviscid, compressible hydrodynamics,

"Supported in part by the Advanced Research Projects Agency as administered by the Rome Air Development Center under contract no. US AF 30(602) 4144."

Analysis of Non-Stationary Flows Using Eulerian and Lagrangian Frameworks

Recent developments have made researchers to reconsider Lagrangian measurement techniques as an alternative to their Eulerian counterpart when investigating non-stationary flows. This thesis advances the state-of-the-art of Lagrangian measurement techniques by pursuing three different objectives: (i) developing new Lagrangian measurement techniques for difficult-to-measure, in situ flow environments; (ii) developing new post-processing strategies designed for unstructured Lagrangian data, as well as providing guidelines towards their use; and (iii) presenting the advantages that the Lagrangian framework has over their Eulerian counterpart in various non-stationary flow problems. Towards the first objective, a large-scale particle tracking velocimetry apparatus is designed for atmospheric surface layer measurements. Towards the second objective, two techniques, one for identifying Lagrangian Coherent Structures (LCS) and the other for characterizing entrainment directly from unstructured Lagrangian data, are developed. Finally, towards the third objective, the advantages of Lagrangian-based measurements are showcased in two unsteady flow problems: the atmospheric surface layer, and entrainment in a non-stationary turbulent flow. Through developing new experimental and post-processing strategies for Lagrangian data, and through showcasing the advantages of Lagrangian data in various non-stationary flows, the thesis works to help investigators to more easily adopt Lagrangian-based measurement techniques.

Blast response and failure analysis of a segmented buried tunnel

Underground tunnels are vulnerable to terrorist attacks which can cause collapse of the tunnel structures or at least extensive damage, requiring lengthy repairs. This paper treats the blast impact on a reinforced concrete segmental tunnel buried in soil under a number of parametric conditions; soil properties, soil cover, distance of explosive from the tunnel centreline and explosive weight and analyses the possible failure patterns. A fully coupled Fluid Structure Interaction (FSI) technique incorporating the Arbitrary Lagrangian-Eulerian (ALE) method is used in this study. Results indicate that the tunnel in saturated soil is more vulnerable to severe damage than that buried in either partially saturated soil or dry soil. The tunnel is also more vulnerable to surface explosions which occur directly above the centre of the tunnel than those that occur at any equivalent distances in the ground away from the tunnel centre. The research findings provide useful information on modeling, analysis, overall tunnel response and failure patterns of segmented tunnels subjected to blast loads. This information will guide future development and application of research in this field.

Sand storms: CFD analysis of Reynolds stress and collision stress of particles near sand bed

Sand storm is a serious environmental threat to humans. Sand particles are transported by saltation and suspension, causing soil erosion in one place and deposition in another. In order to prevent and predict sand storms, the causes and the manners of particle motions must be studied in detail. In this paper a standard k-epsilon model is used for the gas phase simulation and the discrete element method (DEM) is used to predict the movements of particles using an in-house procedure. The data are summarized in an Eulerian-Eulerian regime after simulation to get the statistical particle Reynolds stress and particle collision stress. The results show that for the current case the Reynolds stress and the air shear stress predominate in the region 20-250 mm above the initial sand bed surface. However, in the region below 3 mm, the collision stress must be taken into account in predicting particle movement. (C) 2010 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.

Comparison of lagrangian and eulerian simulations of slurry flows in a sudden expansion

From a review of technical literature, it was not apparent if the Lagrangian or the Eulerian dispersed phase modeling approach was more valid to simulate dilute erosive slurry flow. In this study, both modeling approaches were employed and a comparative analysis of performances and accuracy between the two models was carried out. Due to an impossibility to define, for the Eulerian model already implemented in FLUENT, a set of boundary conditions consistent with the Lagrangian impulsive equations, an Eulerian dispersed phase model was integrated in the FLUENT code using subroutines and user-defined scalar equations. Numerical predictions obtained from the two different approaches for two-phase flow in a sudden expansion were compared with the measured data. Excellent agreement was attained between the predicted and observed fluid and particle velocity in the axial direction and for the kinetic energy. Erosion profiles in a sudden expansion computed using the Lagrangian scheme yielded good qualitative agreement with measured data and predicted a maximum impact angle of 29 deg at the fluid reattachment point. The Eulerian model was adversely affected by the reattachment of the fluid phase to the wall and the simulated erosion profiles were not in agreement with the Lagrangian or measured data. Furthermore, the Eulerian model under-predicted the Lagrangian impact angle at all locations except the reattachment point. © 2010 American Society of Mechanical Engineers.

Investigating Ambrosia Pollen Episodes in Poland Using Back-Trajectory Analysis

Background: The pollen grains of Ambrosia spp. are considered to be important aeroallergens. Previous studies have shown that the long-range transport of Ambrosia pollen to Poland is intermittent and mainly related to the passage of air masses over the Carpathian and Sudetes mountains from sources to the south, e.g. the Czech Republic, Slovakia and Hungary. In this study, Ambrosia pollen counts and back-trajectories from specific episodes in 1999 and 2002 have been analysed with the aim of identifying possible new sources of Ambrosia pollen arriving at three sites in Poland. Method: The combination of Ambrosia pollen measurements (daily average and bi-hourly concentrations) and air mass trajectory calculations were used to investigate two Ambrosia pollen episodes recorded at Rzeszow, Krakow and Poznań on the 4th and 5th September 1999 and 3rd September 2002. Ambrosia pollen counts were recorded by volumetric spore traps of the Hirst design. Trajectories were calculated using the transport model within the Lagrangian air pollution model, ACDEP (Atmospheric Chemistry and Deposition). Results: The collective results of pollen measurements and back-trajectory analysis indicate plumes of Ambrosia pollen travelling up through Poland from the southeast during the investigated episodes. In 1999, the plume was first recorded at Rzeszow in Southeastern Poland during the morning of the 4th September. Its route can be followed as it passed Krakow during the afternoon of the 4th, and later on the 4th and 5th September at Poznań. Similarly, back-trajectories calculated during the morning and afternoon from Krakow and Rzeszow on the 3rd September 2002 indicates that the air masses arrived at these sites from the East or Southeast. Conclusion: This study shows the progress of Ambrosia plumes into Poland from the southeast. Ambrosia pollen release occurs mainly during the day and so a midday peak in Ambrosia pollen concentrations may indicate a local source. However, if the plume of Ambrosia pollen tracked along its northwesterly path over Poland during investigated episodes did not originate from inside Poland, then it is likely that it came from the Ukraine. This identifies a possible new source of ragweed pollen for Poland. Trajectory analysis can only show the path along which an air mass travels, not the specific source area. Further investigation could therefore include source based transport models such as 3D Eulerian atmospheric transport models.

Parallel performance and scalability experiments with the Danish Eulerian Model on the EPCC supercomputers

The Danish Eulerian Model (DEM) is a powerful air pollution model, designed to calculate the concentrations of various dangerous species over a large geographical region (e.g. Europe). It takes into account the main physical and chemical processes between these species, the actual meteorological conditions, emissions, etc.. This is a huge computational task and requires significant resources of storage and CPU time. Parallel computing is essential for the efficient practical use of the model. Some efficient parallel versions of the model were created over the past several years. A suitable parallel version of DEM by using the Message Passing Interface library (AIPI) was implemented on two powerful supercomputers of the EPCC - Edinburgh, available via the HPC-Europa programme for transnational access to research infrastructures in EC: a Sun Fire E15K and an IBM HPCx cluster. Although the implementation is in principal, the same for both supercomputers, few modifications had to be done for successful porting of the code on the IBM HPCx cluster. Performance analysis and parallel optimization was done next. Results from bench marking experiments will be presented in this paper. Another set of experiments was carried out in order to investigate the sensitivity of the model to variation of some chemical rate constants in the chemical submodel. Certain modifications of the code were necessary to be done in accordance with this task. The obtained results will be used for further sensitivity analysis Studies by using Monte Carlo simulation.

Fluid dynamic analysis and modelling of industrial chemical equipment

The research is aimed at contributing to the identification of reliable fully predictive Computational Fluid Dynamics (CFD) methods for the numerical simulation of equipment typically adopted in the chemical and process industries. The apparatuses selected for the investigation, specifically membrane modules, stirred vessels and fluidized beds, were characterized by a different and often complex fluid dynamic behaviour and in some cases the momentum transfer phenomena were coupled with mass transfer or multiphase interactions. Firs of all, a novel modelling approach based on CFD for the prediction of the gas separation process in membrane modules for hydrogen purification is developed. The reliability of the gas velocity field calculated numerically is assessed by comparison of the predictions with experimental velocity data collected by Particle Image Velocimetry, while the applicability of the model to properly predict the separation process under a wide range of operating conditions is assessed through a strict comparison with permeation experimental data. Then, the effect of numerical issues on the RANS-based predictions of single phase stirred tanks is analysed. The homogenisation process of a scalar tracer is also investigated and simulation results are compared to original passive tracer homogenisation curves determined with Planar Laser Induced Fluorescence. The capability of a CFD approach based on the solution of RANS equations is also investigated for describing the fluid dynamic characteristics of the dispersion of organics in water. Finally, an Eulerian-Eulerian fluid-dynamic model is used to simulate mono-disperse suspensions of Geldart A Group particles fluidized by a Newtonian incompressible fluid as well as binary segregating fluidized beds of particles differing in size and density. The results obtained under a number of different operating conditions are compared with literature experimental data and the effect of numerical uncertainties on axial segregation is also discussed.

Multi-Scale Modeling for Image Analysis of Brain Tumor Studies

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multi-scale, multi-physics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlasbased segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

A computational model of hydrogen production by steam reforming of dimethyl ether in a large scale CFB reactor. Part II:parametric analysis

This study presents a computational parametric analysis of DME steam reforming in a large scale Circulating Fluidized Bed (CFB) reactor. The Computational Fluid Dynamic (CFD) model used, which is based on Eulerian-Eulerian dispersed flow, has been developed and validated in Part I of this study [1]. The effect of the reactor inlet configuration, gas residence time, inlet temperature and steam to DME ratio on the overall reactor performance and products have all been investigated. The results have shown that the use of double sided solid feeding system remarkable improvement in the flow uniformity, but with limited effect on the reactions and products. The temperature has been found to play a dominant role in increasing the DME conversion and the hydrogen yield. According to the parametric analysis, it is recommended to run the CFB reactor at around 300 °C inlet temperature, 5.5 steam to DME molar ratio, 4 s gas residence time and 37,104 ml gcat -1 h-1 space velocity. At these conditions, the DME conversion and hydrogen molar concentration in the product gas were both found to be around 80%.

Eulerian-Lagrangian definition of coarse bed-load transport: Theory and verification with low-cost inertial measurement units

Fluvial sediment transport is controlled by hydraulics, sediment properties and arrangement, and flow history across a range of time scales. This physical complexity has led to ambiguous definition of the reference frame (Lagrangian or Eulerian) in which sediment transport is analysed. A general Eulerian-Lagrangian approach accounts for inertial characteristics of particles in a Lagrangian (particle fixed) frame, and for the hydrodynamics in an independent Eulerian frame. The necessary Eulerian-Lagrangian transformations are simplified under the assumption of an ideal Inertial Measurement Unit (IMU), rigidly attached at the centre of the mass of a sediment particle. Real, commercially available IMU sensors can provide high frequency data on accelerations and angular velocities (hence forces and energy) experienced by grains during entrainment and motion, if adequately customized. IMUs are subjected to significant error accu- mulation but they can be used for statistical parametrisation of an Eulerian-Lagrangian model, for coarse sediment particles and over the temporal scale of individual entrainment events. In this thesis an Eulerian-Lagrangian model is introduced and evaluated experimentally. Absolute inertial accelerations were recorded at a 4 Hz frequency from a spherical instrumented particle (111 mm diameter and 2383 kg/m3 density) in a series of entrainment threshold experiments on a fixed idealised bed. The grain-top inertial acceleration entrainment threshold was approximated at 44 and 51 mg for slopes 0.026 and 0.037 respectively. The saddle inertial acceleration entrainment threshold was at 32 and 25 mg for slopes 0.044 and 0.057 respectively. For the evaluation of the complete Eulerian-Lagrangian model two prototype sensors are presented: an idealised (spherical) with a diameter of 90 mm and an ellipsoidal with axes 100, 70 and 30 mm. Both are instrumented with a complete IMU, capable of sampling 3D inertial accelerations and 3D angular velocities at 50 Hz. After signal analysis, the results can be used to parametrize sediment movement but they do not contain positional information. The two sensors (spherical and ellipsoidal) were tested in a series of entrainment experiments, similar to the evaluation of the 111 mm prototype, for a slope of 0.02. The spherical sensor entrained at discharges of 24.8 ± 1.8 l/s while the same threshold for the ellipsoidal sensor was 45.2 ± 2.2 l/s. Kinetic energy calculations were used to quantify the particle-bed energy exchange under fluvial (discharge at 30 l/s) and non-fluvial conditions. All the experiments suggest that the effect of the inertial characteristics of coarse sediments on their motion is comparable to the effect hydrodynamic forces. The coupling of IMU sensors with advanced telemetric systems can lead to the tracking of Lagrangian particle trajectories, at a frequency and accuracy that will permit the testing of diffusion/dispersion models across the range of particle diameters.

Raman microscopy of formamide-intercalated kaolinites treated by controlled-rate thermal analysis technology

Raman spectroscopy of formamide-intercalated kaolinites treated using controlled-rate thermal analysis technology (CRTA), allowing the separation of adsorbed formamide from intercalated formamide in formamide-intercalated kaolinites, is reported. The Raman spectra of the CRTA-treated formamide-intercalated kaolinites are significantly different from those of the intercalated kaolinites, which display a combination of both intercalated and adsorbed formamide. An intense band is observed at 3629 cm-1, attributed to the inner surface hydroxyls hydrogen bonded to the formamide. Broad bands are observed at 3600 and 3639 cm-1, assigned to the inner surface hydroxyls, which are hydrogen bonded to the adsorbed water molecules. The hydroxyl-stretching band of the inner hydroxyl is observed at 3621 cm-1 in the Raman spectra of the CRTA-treated formamide-intercalated kaolinites. The results of thermal analysis show that the amount of intercalated formamide between the kaolinite layers is independent of the presence of water. Significant differences are observed in the CO stretching region between the adsorbed and intercalated formamide.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.