Focal laminar cortical infarcts following aneurysmal subarachnoid haemorrhage

INTRODUCTION: The aim of this prospective study was to analyse small band-like cortical infarcts after subarachnoid haemorrhage (SAH) using magnetic resonance imaging (MRI) with reference to additional digital subtraction angiography (DSA). METHODS: In a 5-year period between January 2002 and January 2007 10 out of 188 patients with aneurysmal SAH were evaluated (one patient Hunt and Hess grade I, one patient grade II, four patients grade III, two patients grade IV, and two patients grade V). The imaging protocol included serially performed MRI with diffusion- and perfusion-weighted images (DWI/PWI) at three time points after aneurysm treatment, and cerebral vasospasm (CVS) was analysed on follow-up DSA on day 7+/-3 after SAH. RESULTS: The lesions were located in the frontal lobe (n=10), in the insular cortex (n=3) and in the parietal lobe (n=1). The band-like infarcts occurred after a mean time interval of 5.8 days (range 3-10 days) and showed unexceptional adjacent thick sulcal clots. Seven out of ten patients with cortical infarcts had no or mild CVS, and in the remaining three patients DSA disclosed moderate (n=2) or severe (n=1) CVS. CONCLUSION: The infarct pattern after aneurysmal SAH includes cortical band-like lesions. In contrast to territorial infarcts or lacunar infarcts in the white matter which develop as a result of moderate or severe proximal and/or distal vasospasm visible on angiography, the cortical band-like lesions adjacent to sulcal clots may also develop without evidence of macroscopic vasospasm, implying a vasospastic reaction of the most distal superficial and intraparenchymal vessels.

Modeling of ore-forming and geoenvironmental systems: Roles of fluid flow and chemical reaction processes

Ore-forming and geoenviromental systems commonly involve coupled fluid flowand chemical reaction processes. The advanced numerical methods and computational modeling have become indispensable tools for simulating such processes in recent years. This enables many hitherto unsolvable geoscience problems to be addressed using numerical methods and computational modeling approaches. For example, computational modeling has been successfully used to solve ore-forming and mine site contamination/remediation problems, in which fluid flow and geochemical processes play important roles in the controlling dynamic mechanisms. The main purpose of this paper is to present a generalized overview of: (1) the various classes and models associated with fluid flow/chemically reacting systems in order to highlight possible opportunities and developments for the future; (2) some more general issues that need attention in the development of computational models and codes for simulating ore-forming and geoenviromental systems; (3) the related progresses achieved on the geochemical modeling over the past 50 years or so; (4) the general methodology for modeling of oreforming and geoenvironmental systems; and (5) the future development directions associated with modeling of ore-forming and geoenviromental systems.

Experimental validation of a nonlinear μ FE model based on cohesive-frictional plasticity for trabecular bone

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of μFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

Modeling asteroid collisions and impact processes

As a complement to experimental and theoretical approaches, numerical modeling has become an important component to study asteroid collisions and impact processes. In the last decade, there have been significant advances in both computational resources and numerical methods. We discuss the present state-of-the-art numerical methods and material models used in "shock physics codes" to simulate impacts and collisions and give some examples of those codes. Finally, recent modeling studies are presented, focussing on the effects of various material properties and target structures on the outcome of a collision.

Pollen records and lithology of profiles from Lobsigensee, Switzerland

Lobsigensee is a small kettle hole lake 15 km north-west of Bern on the Swiss Plateau, at an altitude of 514 m asl. Its surface is 2ha today, its maximum depth 2.7 m; it has no inlet and the overflow functions mainly during snow melting. The area was covered by Rhone ice during the Last Glaciation (map in Fig.2). Local geology, climate and vegetation are summarized in Figure 3A-C, the history of settlement in Figures 5-7. In order to reconstruct the vegetational and environmental history of the lake and its surroundings pollen analysis and other bio- and isotope stratigraphies were applied to twelve profiles cored across the basin with modified Livingstone corers (Fig.3 D). (1) The standard diagram: The central core LQ-90 is described as the standard pollen diagram (Chapter 3) with 10 local pollen assemblage zones of the Late-Glacial (local PAZ Ll to Ll0, from about 16'000(7) to 10'000 years BP) and 20 PAZ of the Holocene (local PAZ L11 to L30), see Figs. 8-10 and 20-24. Local PAZ L 1 to L3 are in the Late-Glacial clay and record the vegetational development after the ice retreat: L1 shows very low pollen concentration and high Pinus percentages due to long-distance transport and reworking; the latter mechanism is corroborated by the findings of thermophilous and pre-Quaternary taxa. Local PAZ L2 has a high di versi ty of non-arboreal pollen (NAP) and reflects the Late-Glacial steppe rich in heliophilous species. Local PAZ L3 is similar but additionally rich in Betula nana and Sal1x, thus reflecting a "shrub tundra". The PAZ L1 to L3 belong to the Oldest Dryas biozone. Local PAZ L4 to L 10 are found in the gyttja of the profundal or in the lake marl of the littoral and record the Late-Glacial forests. L4 is the shrub phase of reforestation with very high Junlperus and rapidly increasing Betula percentages. L5 is the PAZ with a first, L7 with a second dominance of tree-birches, separated by L6 showing a depression in the Betula curve. L4 to L7 can be assigned to the Balling biozone. Possible correlation of the Betula depression to the Older Dryas biozone is discussed. In local PAZ L8 Plnus immigrates and expands. L9 shows a facies difference in that Plnus dominates over Betula in littoral but not in profundal spectra. L8 and L9 belong to the Allerod biozone. In its youngest part the volcanic ash from Laach/Eifel is regularly found (11,000 BP). The local PAZ Ll0 corresponds to the Younger Dryas blozone. The merely slight increase of the NAP indicates that the pine forests of the lowland were not strongly affected by a cooler climate. In order to evaluate the significance of the littoral accumulation of coniferous pollen the littoral profile LQ-150 is compared to the profundal. Radiocarbon stratigraphies derived from different materials are presented in Figures 13 and 14 and in Tables 2 and 3. The hard-water errors in the gyttja samples and the carbonate samples are similar. The samples of terrestrial plant macrofossils are not affected by hard-water errors. Two plateaux of constant age appear in the age-depth relationship; their consequence for biostratigraphy as well as pollen concentration and influx diagrams are discussed. Radiocarbon ages of the Late-Glacial pollen zones are shown in Table 10. The Holocene vegetational history is recorded in the local PAZ L 11 to L30. After a Preboreal (PAZ L11) dominated by pine and birch the expansions of Corylus, Ulmus and Quercus are very rapid. Among these taxa Corylus dominates dur ing the Boreal (PAZ L 12 and L 1 3), whereas the components of the mixed oak forest dominate in the Older Atlantic (PAZ L14 to L16). In the Younger Atlantic (PAZ L 17 to L 19) Fagus and Alnus play an increasing, the mixed oak forest a decreasing role. During the period of local PAZ L19 Neolithic settlers lived on the shore of Lobsigensee. During the Subboreal (PAZ L20 and L21) and the Older Subatlantic (L22 to L25) strong fluctuations of Fagus and often antagonistic peaks of NAP, Alnus, Betula and Corylus can be interpreted as signs of human impact on vegetation. L23 is characterized not only by high values of NAP (especially apophytes and anthropochorous species) but also by the appearance of Juglans, Castanea and Secale which point to the Roman colonization of the area. For a certain period during the Younger Subatlantic (PAZ L26 to L30) the lake was used for retting hemp (Cannabis). Later the dominance of Quercus pollen indicates the importance of wood pastures. The youngest sediments reflect the wide-spread agricultural grass lands and the plantation of Pinus and Picea. Radiocarbon dates for the Holocene are given in Figure 23 and Table 4, the extrapolated ages of the Holocene pollen zones in Table 15. (2) The cross sections: Figures 25 and 26 give a summary of the litho- and palynostratigraphy of the two cross sections. Based on 11 Late-Glacial and 9 Holocene pollen diagrams (in addition to the standard ones), the consistency of the criteria for the definition of the pollen zones is examined in Tables 7 and 8 for the Late-Glacial and in Tables 11 to 14 for the Holocene. Sediment thicknesses across the basin for each pollen zone are presented in these tables as well as in Figures 43 to 45 for the Late-Glacial and in Figures 59 to 65 for the Holocene. Sediment focusing can explain differences between the gyttja cores of the profundal. Focusing is more than compensated for through "stretching" by carbonate precipitation on the littoral terrace. Pollen influx to the cross section are discussed (Chapters 4.1.5. and 4.2.3.). (3) The regional pollen zones: Based on some selected sites between Lake Geneva and Lake Constance regional pollen zones are proposed (Table 16, 17 and 19). (4) Paleoecology: Climatic change in the Late-Glacial can be inferred from Coleoptera, Trichoptera, Chironomidae and d18O of carbonates: a distinct warming is recorded around 12' 600 BP and around 10' 000 BP. The Younger Dryas biozone (10'700-10'000 BP) was the only cooling found in the Late-Glacial. The Betula depression often correlated wi th the Older Dryas biozone was possibl not colder but dryer than the previous period. During the Holocene the lowland site is not very sensitive to the minor climatic changes. Table 22 summarizes climatic and trophic changes before 8'000 BP as deduced from various biostratigraphies studied by a number of authors. Ostracods, Chironomids and fossil pigments indicate that anoxic conditions prevailed during the BoIling (possibly meromixis). Changes in the lake level are illustrated in Figure 74. A first lake-level lowering occurred in the early Holocene (10'000 to 9'000 BP), a second during the Atlantic (about 6'800 to 5'200 BP). The first "shrinking" of the lake volume resulted in a eutrophication recorded by laminations in the profundal and by pigments of Cyanophyceae. The second fall in water level corresponds to an increase of Nymphaeaceae. Human impact can be inferred in three ways: eutrophication of the lake (since the Neolithic), changes of terrestrial vegetation by deforestations (cyclicity of Fagus, see Figures 78 to 80), and enhanced erosion (increasing sedimentation rates by inwashed clay, particularly since the Roman Colonization, see Figures 49 and 81). Summary: This paper was planned as the final report on Lobsigensee. However, a number of issues are not answered but can only be asked more precisely, for example: (1) For the two periods with the highest rates of change, Le. the Bolling and the Preboreal biozones, pollen influx may reflect vegetation dynamics. Detailed investigations of these periods in annually laminated sediments are planned. (2) Biostratigraphies other than palynostratigraphy are needed to estimate the degree of linkage or independence in the development of terrestrial and lacustrine ecosystems. Often our sampling intervals were not identical, thus influencing our temporal resolution. (3) 6180- and 14C-stratigraPhies with high resolution will elucidate the leads and lags of these dynamic periods. Plateaux of constant age in the age-depth relationship have a strong bearing on both biological and geophysical understanding of Late-Glacial and early Holocene developments. (4) Numerical methods applied to the pollen diagrams of the cross section will help to quantify the significance of similari ties and dissimilarities across a single basin (with Prof. Birks). (5) Numerical methods applied to different sites on the Swiss Plateau and on the transect across the Alps will be helpful in evaluating the influence of different environmental factors (with Prof. Birks). (6) A new map 1: 1000 with 50cm-contour lines prov ided by Prof. Zurbuchen will be combined with a grid of cores sampling the transition from lake marl to peat enabling us to calculate paleo-volumes of the lake. This is interesting for the two "shrinking periods" (in Fig. 74A numbers 2-6 and 7-10), both accompanied by eutrophication. The pal eo-volume during the Neoli thic set tlement of the Cortaillod culture linked wi th an est l.mate of trophic change derived from diatoms (Prof. Smol in prep.) could possibly give an indication of the size of the human population of this period. (7) For the period with the antagonism between Fagus peaks and ABC-peaks close collaboration between palynologists, geochemists and archeologists should enable us to determine the influence of prehistoric and historic people on vegetation (collaboration with Prof. Stockli and Prof. Herzig). (8) The core LL-75 taken with a "cold letter box" will be analysed for major and trace elements by Dr. Sturm for 210pb and 137Cs by Prof.von Gunten and for pollen. We will see if our local PAZ L30 really corresponds to the surface sediment and if the small seepage lake reflects modern pollution.

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

Cargas explosivas sobre pilares y forjados modelados con elementos tipo viga y lámina

Los fenómenos de impacto y explosión sobre estructuras de hormigón tienen efectos en muchos casos catastróficos a pesar de su reducida probabilidad. Las estructuras de hormigón no suelen estar diseñadas para resistir este tipo de solicitaciones dinámicas. El análisis numérico mediante elementos finitos con integración explícita permite una aproximación suficiente a los efectos de la onda explosiva sobre pilares y forjados de estructuras reticuladas de hormigón. Los materiales recientemente implementados en LS-Dyna para hormigón como el CSCM [1], para elementos de continuo 3D, y la formulación que proporciona la debida compatibilidad con los elementos viga de acero dispuestos de forma segregada, permite estudiar de forma realista modelos detallados de pilares y forjados. Pero las limitaciones computacionales hacen inviable emplear estos métodos en estructuras completas. Como alternativa es posible usar modelos de elementos estructurales de vigas y láminas para el análisis de estas estructuras. Sin embargo es necesario un adecuado ajuste de parámetros y propiedades en estos modelos. Este trabajo muestra un método con en el que obtener modelos de elementos estructurales, elementos viga y lámina, usando modelos de material [2] adecuados para ellos, junto a un procedimiento para incluir la armadura de forma adecuada. Utilizando este método es posible representar con suficiente aproximación el comportamiento de modelos detallados realistas de forjados y pilares de estructuras reticuladas de hormigón frente a acciones explosivas, posibilitando el análisis de una estructura completa frente a explosión.

Análisis de los parámetros que definen un modelo de daño aplicado al deterioro de la capacidad resistente de paredes arteriales

El comportamiento mecánico de las paredes arteriales constituye un ingrediente fundamental para entender enfermedades cardiovasculares como las disecciones aórticas o la fisuración de placas de ateroma. El fenómeno mecánico que subyace en estas enfermedades es el fallo constitutivo del material, que puede deberse a la degradación de las propiedades mecánicas del tejido, a la aplicación de una carga excesiva o la concentración de tensiones que tiene lugar en traumatizaciones de la íntima. Un modelo que busque reproducir este comportamiento biomecánico debe incorporar además la información histológica de las paredes arteriales, formadas en una primera aproximación por una matriz acuosa incompresible y dos familias de fibras dispuestas helicoidalmente. De entre los modelos mecánicos existentes cabe destacar los modelos de daño, que con un conjunto pequeño de parámetros son capaces de reproducir la degradación de la capacidad resistente tanto en la matriz como en las fibras. Esta simplicidad convierte a esta clase de modelos en una herramienta muy útil para estudiar este fenómeno. En este trabajo se presenta un estudio para un modelo de daño de los parámetros que lo definen, asociados a la imposición de la incompresibilidad de la matriz acuosa y a las evoluciones del daño en dicha matriz y en las dos familias de fibras. Se incluye una revisión de las técnicas existentes para imponer la incompresibilidad, un resumen bibliográfico de la caracterización experimental de los parámetros materiales y un análisis de su influencia en el comportamiento estructural de la aorta, analizada mediante ensayos de tracción

Dynamic Effects On Culverts for High Speed Trains

Culverts are very common in recent railway lines. Wild life corridors and drainage conducts often fall in this category of partially buried structures. Their dynamic behavior has received far less attention than other structures such as bridges but its large number makes that study an interesting challenge from the point of view of safety and savings. In this paper a complete study of a culvert, including on-site measurements as well as numerical modelling, will be presented. The structure belongs to the high speed railway line linking Segovia and Valladolid, in Spain. The line was opened to traffic in 2004. Its dimensions (3x3m) are the most frequent along the line. Other factors such as reduced overburden (0.6m) and an almost right angle with the track axis make it an interesting example to extract generalized conclusions. On site measurements have been performed in the structure recording the dynamic response at selected points of the structure during the passage of high speed trains at speeds ranging between 200 and 300km/h. The measurements by themselves provide a good insight into the main features of the dynamic behaviour of the structure. A 3D finite element model of the structure, representing its key features was also studied as it allows further understanding of the dynamic response to the train loads . In the paper the discrepancies between predicted and measured vibration levels will be analyzed and some advices on numerical modelling will be proposed

Towards the foundation of a global modes concept

A contribution is presented, intended to provide theoretical foundations for the ongoing efforts to employ global instability theory for the analysis of the classic boundary-layer flow, and address the associated issue of appropriate inflow/outflow boundary conditions to close the PDE-based global eigenvalue problem in open flows. Starting from a theoretically clean and numerically simple application, in which results are also known analytically and thus serve as a guidance for the assessment of the performance of the numerical methods employed herein, a sequence of issues is systematically built into the target application, until we arrive at one representative of open systems whose instability is presently addressed by global linear theory applied to open flows, the latter application being neither tractable theoretically nor straightforward to solve by numerical means. Experience gained along the way is documented. It regards quantification of the depar- ture of the numerical solution from the analytical one in the simple problem, the generation of numerical boundary layers at artificially truncated boundaries, no matter how far the latter are placed from the region of highest flow gradients and, ultimately the impracti- cally large number of (direct and adjoint) modes necessary to project an arbitrary initial perturbation and follow its temporal evolution by a global analysis approach, a finding which may question the purported robustness reported in the literature of the recovery of optimal perturbations as part of global analyses yielding under-resolved eigenspectra.

Modelling Crowd Load for Floor Vibration Analysis

The dynamic floor loads induced by crowds in gymnasium or stadium structures are commonly modelled by superposition of the individual contributions using reduction factors for the different Fourier coefficients. These Fourier coefficients and the reduction factors are calculated using full scale measurements. Generally the testing is performed on platforms or structures that can be considered rigid, such that the natural frequencies are higher than the frequencies of the spectator movement. In this paper we shall present the testing done on a structure that used to be a gymnasium as well as the procedure used to identify its dynamic properties and a first evaluation of the socalled “group effect”.

Statistical mechanics approaches to granular media: between micromechanics and macromechanics

The mechanical behavior of granular materials has been traditionally approached through two theoretical and computational frameworks: macromechanics and micromechanics. Macromechanics focuses on continuum based models. In consequence it is assumed that the matter in the granular material is homogeneous and continuously distributed over its volume so that the smallest element cut from the body possesses the same physical properties as the body. In particular, it has some equivalent mechanical properties, represented by complex and non-linear constitutive relationships. Engineering problems are usually solved using computational methods such as FEM or FDM. On the other hand, micromechanics is the analysis of heterogeneous materials on the level of their individual constituents. In granular materials, if the properties of particles are known, a micromechanical approach can lead to a predictive response of the whole heterogeneous material. Two classes of numerical techniques can be differentiated: computational micromechanics, which consists on applying continuum mechanics on each of the phases of a representative volume element and then solving numerically the equations, and atomistic methods (DEM), which consist on applying rigid body dynamics together with interaction potentials to the particles. Statistical mechanics approaches arise between micro and macromechanics. It tries to state which the expected macroscopic properties of a granular system are, by starting from a micromechanical analysis of the features of the particles and the interactions. The main objective of this paper is to introduce this approach.

A C1 finite element for Kirchhoff plate bending

After a short introduction the possibilities and limitations of polynomial simple elements with C1 continuity are discussed with reference to plate bending analysis. A family of this kind of elements is presented.. These elements are applied to simple cases in order to assess their computational efficiency. Finally some conclusions are shown, and future research is also proposed.

Aplicación de la regresión logística en la predicción empírica de fenómenos complejos en obras subterráneas : squeezing y rotura de pilares de carbón

La heterogeneidad del medio geológico introduce en el proyecto de obra subterránea un alto grado de incertidumbre que debe ser debidamente gestionado a fin de reducir los riesgos asociados, que son fundamentalmente de tipo geotécnico. Entre los principales problemas a los que se enfrenta la Mecánica de Rocas moderna en el ámbito de la construcción subterránea, se encuentran la fluencia de roca en túneles (squeezing) y la rotura de pilares de carbón. Es ampliamente conocido que su aparición causa importantes perjuicios en el coste y la seguridad de los proyectos por lo que su estudio, ha estado tradicionalmente vinculado a la predicción de su ocurrencia. Entre las soluciones existentes para la determinación de estos problemas se encuentran las que se basan en métodos analíticos y numéricos. Estas metodologías son capaces de proporcionar un alto nivel de representatividad respecto del comportamiento geotécnico real, sin embargo, su utilización solo es posible cuando se dispone de una suficiente caracterización geotécnica y por tanto de una detallada definición de los parámetros que alimentan los complejos modelos constitutivos y criterios de rotura que los fenómenos estudiados requieren. Como es lógico, este nivel de definición solo es posible cuando se alcanzan etapas avanzadas de proyecto, incluso durante la propia construcción, a fin de calibrar adecuadamente los parámetros introducidos en los modelos, lo que supone una limitación de uso en etapas iniciales, cuando su predicción tiene verdadero sentido. Por su parte, los métodos empíricos permiten proporcionar soluciones a estos complejos problemas de un modo sencillo, con una baja parametrización y, dado su eminente enfoque observacional, de gran fiabilidad cuando se implementan sobre condiciones de contorno similares a las originales. La sencillez y escasez de los parámetros utilizados permiten a estas metodologías ser utilizadas desde las fases preliminares del proyecto, ya que estos constituyen en general, información habitual de fácil y económica adquisición. Este aspecto permite por tanto incorporar la predicción desde el principio del proceso de diseño, anticipando el riesgo en origen. En esta tesis doctoral, se presenta una nueva metodología empírica que sirve para proporcionar predicciones para la ocurrencia de squeezing y el fallo de pilares de carbón basada en una extensa recopilación de información de casos reales de túneles y minas en las que ambos fenómenos fueron evaluados. Esta información, recogida de referencias bibliográficas de prestigio, ha permitido recopilar una de las más extensas bases de datos existentes hasta la fecha relativa a estos fenómenos, lo que supone en sí mismo una importante contribución sobre el estado del arte. Con toda esta información, y con la ayuda de la teoría de clasificadores estadísticos, se ha implementado sobre las bases de datos un clasificador lineal de tipo regresión logística que permite hacer predicciones sobre la ocurrencia de ambos fenómenos en términos de probabilidad, y por tanto ponderar la incertidumbre asociada a la heterogeneidad incorporada por el medio geológico. Este aspecto del desarrollo es el verdadero valor añadido proporcionado por la tesis y la principal ventaja de la solución propuesta respecto de otras metodologías empíricas. Esta capacidad de ponderación probabilística permite al clasificador constituir una solución muy interesante como metodología para la evaluación de riesgo geotécnico y la toma de decisiones. De hecho, y como ejercicio de validación práctica, se ha implementado la solución desarrollada en un modelo coste-beneficio asociado a la optimización del diseño de pilares involucrados en una de mina “virtual” explotada por tajos largos. La capacidad del clasificador para cuantificar la probabilidad de fallo del diseño, junto con una adecuada cuantificación de las consecuencias de ese fallo, ha permitido definir una ley de riesgo que se ha incorporado al balance de costes y beneficios, que es capaz, a partir del redimensionamiento iterativo del sistema de pilares y de la propia configuración de la mina, maximizar el resultado económico del proyecto minero bajo unas condiciones de seguridad aceptables, fijadas de antemano. Geological media variability introduces to the subterranean project a high grade of uncertainty that should be properly managed with the aim to reduce the associated risks, which are mainly geotechnical. Among the major problems facing the modern Rock Mechanics in the field of underground construction are both, the rock squeezing while tunneling and the failure of coal pillars. Given their harmfulness to the cost and safety of the projects, their study has been traditionally linked to the determination of its occurrence. Among the existing solutions for the determination of these problems are those that are based on analytical and numerical methods. Those methodologies allow providing a high level of reliability of the geotechnical behavior, and therefore a detailed definition of the parameters that feed the complex constitutive models and failure criteria that require the studied phenomena. Obviously, this level of definition is only possible when advanced stages of the project are achieved and even during construction in order to properly calibrate the parameters entered in the models, which suppose a limited use in early stages, when the prediction has true sense. Meanwhile, empirical methods provide solutions to these complex problems in a simple way, with low parameterization and, given his observational scope, with highly reliability when implemented on similar conditions to the original context. The simplicity and scarcity of the parameters used allow these methodologies be applied in the early stages of the project, since that information should be commonly easy and cheaply to get. This aspect can therefore incorporate the prediction from the beginning of the design process, anticipating the risk beforehand. This thesis, based on the extensive data collection of case histories of tunnels and underground mines, presents a novel empirical approach used to provide predictions for the occurrence of both, squeezing and coal pillars failures. The information has been collected from prestigious references, providing one of the largest databases to date concerning phenomena, a fact which provides an important contribution to the state of the art. With all this information, and with the aid of the theory of statistical classifiers, it has been implemented on both databases, a type linear logistic regression classifier that allows predictions about the occurrence of these phenomena in terms of probability, and therefore weighting the uncertainty associated with geological variability. This aspect of the development is the real added value provided by the thesis and the main advantage of the proposed solution over other empirical methodologies. This probabilistic weighting capacity, allows being the classifier a very interesting methodology for the evaluation of geotechnical risk and decision making. In fact, in order to provide a practical validation, we have implemented the developed solution within a cost-benefit analysis associated with the optimization of the design of coal pillar systems involved in a "virtual" longwall mine. The ability of the classifier to quantify the probability of failure of the design along with proper quantification of the consequences of that failure, has allowed defining a risk law which is introduced into the cost-benefits model, which is able, from iterative resizing of the pillar system and the configuration of the mine, maximize the economic performance of the mining project under acceptable safety conditions established beforehand.

Towards cost reduction in concentrating solar power: innovative design for an efficient fresnel based solar field

Energía termosolar (de concentración) es uno de los nombres que hacen referencia en español al término inglés “concentrating solar power”. Se trata de una tecnología basada en la captura de la potencia térmica de la radiación solar, de forma que permita alcanzar temperaturas capaces de alimentar un ciclo termodinámico convencional (o avanzado); el futuro de esta tecnología depende principalmente de su capacidad para concentrar la radiación solar de manera eficiente y económica. La presente tesis está orientada hacia la resolución de ciertos problemas importantes relacionados con este objetivo. La mencionada necesidad de reducir costes en la concentración de radiación solar directa, asegurando el objetivo termodinámico de calentar un fluido hasta una determinada temperatura, es de vital importancia. Los colectores lineales Fresnel han sido identificados en la literatura científica como una tecnología con gran potencial para alcanzar esta reducción de costes. Dicha tecnología ha sido seleccionada por numerosas razones, entre las que destacan su gran libertad de diseño y su actual estado inmaduro. Con el objetivo de responder a este desafío se desarrollado un detallado estudio de las propiedades ópticas de los colectores lineales Fresnel, para lo cual se han utilizado métodos analíticos y numéricos de manera combinada. En primer lugar, se han usado unos modelos para la predicción de la localización y la irradiación normal directa del sol junto a unas relaciones analíticas desarrolladas para estudiar el efecto de múltiples variables de diseño en la energía incidente sobre los espejos. Del mismo modo, se han obtenido analíticamente los errores debidos al llamado “off-axis aberration”, a la apertura de los rayos reflejados en los espejos y a las sombras y bloqueos entre espejos. Esto ha permitido la comparación de diferentes formas de espejo –planos, circulares o parabólicos–, así como el diseño preliminar de la localización y anchura de los espejos y receptor sin necesidad de costosos métodos numéricos. En segundo lugar, se ha desarrollado un modelo de trazado de rayos de Monte Carlo con el objetivo de comprobar la validez del estudio analítico, pero sobre todo porque este no es preciso en el estudio de la reflexión en espejos. El código desarrollado está específicamente ideado para colectores lineales Fresnel, lo que ha permitido la reducción del tiempo de cálculo en varios órdenes de magnitud en comparación con un programa comercial más general. Esto justifica el desarrollo de un nuevo código en lugar de la compra de una licencia de otro programa. El modelo ha sido usado primeramente para comparar la intensidad de flujo térmico y rendimiento de colectores Fresnel, con y sin reflector secundario, con los colectores cilíndrico parabólicos. Finalmente, la conjunción de los resultados obtenidos en el estudio analítico con el programa numérico ha sido usada para optimizar el campo solar para diferentes orientaciones –Norte-Sur y Este-Oeste–, diferentes localizaciones –Almería y Aswan–, diferentes inclinaciones hacia el Trópico –desde 0 deg hasta 32 deg– y diferentes mínimos de intensidad del flujo en el centro del receptor –10 kW/m2 y 25 kW/m2–. La presente tesis ha conducido a importantes descubrimientos que deben ser considerados a la hora de diseñar un campo solar Fresnel. En primer lugar, los espejos utilizados no deben ser plano, sino cilíndricos o parabólicos, ya que los espejos curvos implican mayores concentraciones y rendimiento. Por otro lado, se ha llegado a la conclusión de que la orientación Este-Oeste es más propicia para localizaciones con altas latitudes, como Almería, mientras que en zonas más cercanas a los trópicos como Aswan los campos Norte-Sur conducen a mayores rendimientos. Es de destacar que la orientación Este-Oeste requiere aproximadamente la mitad de espejos que los campos Norte-Sur, puediendo estar inclinados hacia los Trópicos para mejorar el rendimiento, y que alcanzan parecidos valores de intensidad térmica en el receptor todos los días a mediodía. Sin embargo, los campos con orientación Norte-Sur permiten un flujo más constante a lo largo de un día. Por último, ha sido demostrado que el uso de diseños pre-optimizados analíticamente, con anchura de espejos y espaciado entre espejos variables a lo ancho del campo, pueden implicar aumentos de la energía generada por metro cuadrado de espejos de hasta el 6%. El rendimiento óptico anual de los colectores cilíndrico parabólicos es 23 % mayor que el rendimiento de los campos Fresnel en Almería, mientras que la diferencia es de solo 9 % en Aswan. Ello implica que, para alcanzar el mismo precio de electricidad que la tecnología de referencia, la reducción de costes de instalación por metro cuadrado de espejo debe estar entre el 10 % y el 25 %, y que los colectores lineales Fresnel tienen más posibilidades de ser desarrollados en zonas de bajas latitudes. Como consecuencia de los estudios desarrollados en esta tesis se ha patentado un sistema de almacenamiento que tiene en cuenta la variación del flujo térmico en el receptor a lo largo del día, especialmente para campos con orientación Este-Oeste. Este invento permitiría el aprovechamiento de la energía incidente durante más parte del año, aumentando de manera apreciable los rendimientos óptico y térmico. Abstract Concentrating solar power is the common name of a technology based on capturing the thermal power of solar radiation, in a suitable way to reach temperatures able to activate a conventional (or advanced) thermodynamic cycle to generate electricity; this quest mainly depends on our ability to concentrate solar radiation in a cheap and efficient way. The present thesis is focused to highlight and help solving some of the important issues related to this problem. The need of reducing costs in concentrating the direct solar radiation, but without jeopardizing the thermodynamic objective of heating a fluid up to the required temperature, is of prime importance. Linear Fresnel collectors have been identified in the scientific literature as a technology with high potential to reach this cost reduction. This technology has been selected because of a number of reasons, particularly the degrees of freedom of this type of concentrating configuration and its current immature state. In order to respond to this challenge, a very detailed exercise has been carried out on the optical properties of linear Fresnel collectors. This has been done combining analytic and numerical methods. First, the effect of the design variables on the ratio of energy impinging onto the reflecting surface has been studied using analytically developed equations, together with models that predict the location and direct normal irradiance of the sun at any moment. Similarly, errors due to off-axis aberration, to the aperture of the reflected energy beam and to shading and blocking effects have been obtained analytically. This has allowed the comparison of different shapes of mirrors –flat, cylindrical or parabolic–, as well as a preliminary optimization of the location and width of mirrors and receiver with no need of time-consuming numerical models. Second, in order to prove the validity of the analytic results, but also due to the fact that the study of the reflection process is not precise enough when using analytic equations, a Monte Carlo Ray Trace model has been developed. The developed code is designed specifically for linear Fresnel collectors, which has reduced the computing time by several orders of magnitude compared to a wider commercial software. This justifies the development of the new code. The model has been first used to compare radiation flux intensities and efficiencies of linear Fresnel collectors, both multitube receiver and secondary reflector receiver technologies, with parabolic trough collectors. Finally, the results obtained in the analytic study together with the numeric model have used in order to optimize the solar field for different orientations –North-South and East-West–, different locations –Almería and Aswan–, different tilts of the field towards the Tropic –from 0 deg to 32 deg– and different flux intensity minimum requirements –10 kW/m2 and 25 kW/m2. This thesis work has led to several important findings that should be considered in the design of Fresnel solar fields. First, flat mirrors should not be used in any case, as cylindrical and parabolic mirrors lead to higher flux intensities and efficiencies. Second, it has been concluded that, in locations relatively far from the Tropics such as Almería, East-West embodiments are more efficient, while in Aswan North- South orientation leads to a higher annual efficiency. It must be noted that East-West oriented solar fields require approximately half the number of mirrors than NS oriented fields, can be tilted towards the Equator in order to increase the efficiency and attain similar values of flux intensity at the receiver every day at midday. On the other hand, in NS embodiments the flux intensity is more even during each single day. Finally, it has been proved that the use of analytic designs with variable shift between mirrors and variable width of mirrors across the field can lead to improvements in the electricity generated per reflecting surface square meter up to 6%. The annual optical efficiency of parabolic troughs has been found to be 23% higher than the efficiency of Fresnel fields in Almería, but it is only around 9% higher in Aswan. This implies that, in order to attain the same levelized cost of electricity than parabolic troughs, the required reduction of installation costs per mirror square meter is in the range of 10-25%. Also, it is concluded that linear Fresnel collectors are more suitable for low latitude areas. As a consequence of the studies carried out in this thesis, an innovative storage system has been patented. This system takes into account the variation of the flux intensity along the day, especially for East-West oriented solar fields. As a result, the invention would allow to exploit the impinging radiation along longer time every day, increasing appreciably the optical and thermal efficiencies.