Pinus uncinata Ramond ex DC.: Conservation and proposal of management at its southwestern limit of distribution (Castillo de Vinuesa, Soria, Spain)

Pinus uncinata forms forests in the centre and southwest of the Alps and in the subalpine Pyrenees (at around 1700 – 2600 m) (Costa Tenorio et al., 1997). The species reaches the southwestern limit of its distribution at the top of Mount Castillo de Vinuesa (Soria, Spain). The small population on this mountain occupies just 66 ha, but is very important from a geobotanical viewpoint since it is just one of two populations (the other being in the Sierra de Gúdar range in Teruel, Spain) isolated from the main area where the species is found in the Iberian Peninsula (The Pyrenees)

El problema del rozamiento en el análisis de estructuras de fábrica mediante modelos de sólidos rígidos

La evaluación de la seguridad de estructuras antiguas de fábrica es un problema abierto.El material es heterogéneo y anisótropo, el estado previo de tensiones difícil de conocer y las condiciones de contorno inciertas. A comienzos de los años 50 se demostró que el análisis límite era aplicable a este tipo de estructuras, considerándose desde entonces como una herramienta adecuada. En los casos en los que no se produce deslizamiento la aplicación de los teoremas del análisis límite estándar constituye una herramienta formidable por su simplicidad y robustez. No es necesario conocer el estado real de tensiones. Basta con encontrar cualquier solución de equilibrio, y que satisfaga las condiciones de límite del material, en la seguridad de que su carga será igual o inferior a la carga real de inicio de colapso. Además esta carga de inicio de colapso es única (teorema de la unicidad) y se puede obtener como el óptimo de uno cualquiera entre un par de programas matemáticos convexos duales. Sin embargo, cuando puedan existir mecanismos de inicio de colapso que impliquen deslizamientos, cualquier solución debe satisfacer tanto las restricciones estáticas como las cinemáticas, así como un tipo especial de restricciones disyuntivas que ligan las anteriores y que pueden plantearse como de complementariedad. En este último caso no está asegurada la existencia de una solución única, por lo que es necesaria la búsqueda de otros métodos para tratar la incertidumbre asociada a su multiplicidad. En los últimos años, la investigación se ha centrado en la búsqueda de un mínimo absoluto por debajo del cual el colapso sea imposible. Este método es fácil de plantear desde el punto de vista matemático, pero intratable computacionalmente, debido a las restricciones de complementariedad 0 y z 0 que no son ni convexas ni suaves. El problema de decisión resultante es de complejidad computacional No determinista Polinomial (NP)- completo y el problema de optimización global NP-difícil. A pesar de ello, obtener una solución (sin garantía de exito) es un problema asequible. La presente tesis propone resolver el problema mediante Programación Lineal Secuencial, aprovechando las especiales características de las restricciones de complementariedad, que escritas en forma bilineal son del tipo y z = 0; y 0; z 0 , y aprovechando que el error de complementariedad (en forma bilineal) es una función de penalización exacta. Pero cuando se trata de encontrar la peor solución, el problema de optimización global equivalente es intratable (NP-difícil). Además, en tanto no se demuestre la existencia de un principio de máximo o mínimo, existe la duda de que el esfuerzo empleado en aproximar este mínimo esté justificado. En el capítulo 5, se propone hallar la distribución de frecuencias del factor de carga, para todas las soluciones de inicio de colapso posibles, sobre un sencillo ejemplo. Para ello, se realiza un muestreo de soluciones mediante el método de Monte Carlo, utilizando como contraste un método exacto de computación de politopos. El objetivo final es plantear hasta que punto está justificada la busqueda del mínimo absoluto y proponer un método alternativo de evaluación de la seguridad basado en probabilidades. Las distribuciones de frecuencias, de los factores de carga correspondientes a las soluciones de inicio de colapso obtenidas para el caso estudiado, muestran que tanto el valor máximo como el mínimo de los factores de carga son muy infrecuentes, y tanto más, cuanto más perfecto y contínuo es el contacto. Los resultados obtenidos confirman el interés de desarrollar nuevos métodos probabilistas. En el capítulo 6, se propone un método de este tipo basado en la obtención de múltiples soluciones, desde puntos de partida aleatorios y calificando los resultados mediante la Estadística de Orden. El propósito es determinar la probabilidad de inicio de colapso para cada solución.El método se aplica (de acuerdo a la reducción de expectativas propuesta por la Optimización Ordinal) para obtener una solución que se encuentre en un porcentaje determinado de las peores. Finalmente, en el capítulo 7, se proponen métodos híbridos, incorporando metaheurísticas, para los casos en que la búsqueda del mínimo global esté justificada. Abstract Safety assessment of the historic masonry structures is an open problem. The material is heterogeneous and anisotropic, the previous state of stress is hard to know and the boundary conditions are uncertain. In the early 50's it was proven that limit analysis was applicable to this kind of structures, being considered a suitable tool since then. In cases where no slip occurs, the application of the standard limit analysis theorems constitutes an excellent tool due to its simplicity and robustness. It is enough find any equilibrium solution which satisfy the limit constraints of the material. As we are certain that this load will be equal to or less than the actual load of the onset of collapse, it is not necessary to know the actual stresses state. Furthermore this load for the onset of collapse is unique (uniqueness theorem), and it can be obtained as the optimal from any of two mathematical convex duals programs However, if the mechanisms of the onset of collapse involve sliding, any solution must satisfy both static and kinematic constraints, and also a special kind of disjunctive constraints linking the previous ones, which can be formulated as complementarity constraints. In the latter case, it is not guaranted the existence of a single solution, so it is necessary to look for other ways to treat the uncertainty associated with its multiplicity. In recent years, research has been focused on finding an absolute minimum below which collapse is impossible. This method is easy to set from a mathematical point of view, but computationally intractable. This is due to the complementarity constraints 0 y z 0 , which are neither convex nor smooth. The computational complexity of the resulting decision problem is "Not-deterministic Polynomialcomplete" (NP-complete), and the corresponding global optimization problem is NP-hard. However, obtaining a solution (success is not guaranteed) is an affordable problem. This thesis proposes solve that problem through Successive Linear Programming: taking advantage of the special characteristics of complementarity constraints, which written in bilinear form are y z = 0; y 0; z 0 ; and taking advantage of the fact that the complementarity error (bilinear form) is an exact penalty function. But when it comes to finding the worst solution, the (equivalent) global optimization problem is intractable (NP-hard). Furthermore, until a minimum or maximum principle is not demonstrated, it is questionable that the effort expended in approximating this minimum is justified. XIV In chapter 5, it is proposed find the frequency distribution of the load factor, for all possible solutions of the onset of collapse, on a simple example. For this purpose, a Monte Carlo sampling of solutions is performed using a contrast method "exact computation of polytopes". The ultimate goal is to determine to which extent the search of the global minimum is justified, and to propose an alternative approach to safety assessment based on probabilities. The frequency distributions for the case study show that both the maximum and the minimum load factors are very infrequent, especially when the contact gets more perfect and more continuous. The results indicates the interest of developing new probabilistic methods. In Chapter 6, is proposed a method based on multiple solutions obtained from random starting points, and qualifying the results through Order Statistics. The purpose is to determine the probability for each solution of the onset of collapse. The method is applied (according to expectations reduction given by the Ordinal Optimization) to obtain a solution that is in a certain percentage of the worst. Finally, in Chapter 7, hybrid methods incorporating metaheuristics are proposed for cases in which the search for the global minimum is justified.

Holocene distribution of woody taxa at the westernmost limit of the Circumboreal/Mediterranean boundary: evidence from wood remains

Información paleoecológica relativa a la distribución de táxones arbóreos durante el holoceno, basado en restos fósiles leñosos.

Critical comparative evaluation of design methodologies for the ultimate limit state of fatigue in concrete structures

La aparición de la fatiga ha sido ampliamente investigada en el acero y en otros materiales metálicos, sin embargo no se conoce en tanta profundidad en el hormigón estructural. Esto crea falta de uniformidad y enfoque en el proceso de verificación de estructuras de hormigón para el estado límite último de la fatiga. A medida que se llevan a cabo más investigaciones, la información sobre los parámetros que afectan a la fatiga en el hormigón comienzan a ser difundidos e incluso los que les afectan de forma indirecta. Esto conlleva a que se estén incorporando en las guías de diseño de todo el mundo, a pesar de que la comprobación del estado límite último no se trata por igual entre los distintos órganos de diseño. Este trabajo presentará un conocimiento básico del fenómeno de la fatiga, qué lo causa y qué condiciones de carga o propiedades materiales amplían o reducen la probabilidad de fallo por fatiga. Cuatro distintos códigos de diseño serán expuestos y su proceso de verificación ha sido examinado, comparados y valorados cualitativa y cuantitativamente. Una torre eólica, como ejemplo, fue analizada usando los procedimientos de verificación como se indica en sus respectivos códigos de referencia. The occurrence of fatigue has been extensively researched in steel and other metallic materials it is however, not as broadly understood in concrete. This produces a lack of uniformity in the approach and process in the verification of concrete structures for the ultimate limit state of fatigue. As more research is conducted and more information is known about the parameters which cause, propagate, and indirectly affect fatigue in concrete, they are incorporated in design guides around the world. Nevertheless, this ultimate limit state verification is not addressed equally by various design governing bodies. This report presents a baseline understanding of what the phenomenon of fatigue is, what causes it, and what loading or material conditions amplify or reduce the likelihood of fatigue failure. Four different design codes are exposed and their verification process has been examined, compared and evaluated both qualitatively and quantitatively. Using a wind turbine tower structure as case study, this report presents calculated results following the verification processes as instructed in the respective reference codes.

Precipitated iron: a limit on gettering efficacy in multicrystalline silicon

A phosphorus diffusion gettering model is used to examine the efficacy of a standard gettering process on interstitial and precipitated iron in multicrystalline silicon. The model predicts a large concentration of precipitated iron remaining after standard gettering for most as-grown iron distributions. Although changes in the precipitated iron distribution are predicted to be small, the simulated post-processing interstitial iron concentration is predicted to depend strongly on the as-grown distribution of precipitates, indicating that precipitates must be considered as internal sources of contamination during processing. To inform and validate the model, the iron distributions before and after a standard phosphorus diffusion step are studied in samples from the bottom, middle, and top of an intentionally Fe-contaminated laboratory ingot. A census of iron-silicide precipitates taken by synchrotron-based X-ray fluorescence microscopy confirms the presence of a high density of iron-silicide precipitates both before and after phosphorus diffusion. A comparable precipitated iron distribution was measured in a sister wafer after hydrogenation during a firing step. The similar distributions of precipitated iron seen after each step in the solar cell process confirm that the effect of standard gettering on precipitated iron is strongly limited as predicted by simulation. Good agreement between the experimental and simulated data supports the hypothesis that gettering kinetics is governed by not only the total iron concentration but also by the distribution of precipitated iron. Finally, future directions based on the modeling are suggested for the improvement of effective minority carrier lifetime in multicrystalline silicon solar cells.

Can seed production and restricted dispersal limit recruitment in Pinus pinaster Aiton from the Spanish Northern Plateau?

Natural regeneration faces increasing difficulties in dry forests from the Mediterranean basin, including for normally well-regenerating species such as maritime pine (Pinus pinaster Aiton). In this paper, we studied female fertility, seed dispersal and spread rates in P. pinaster from the Spanish Northern Plateau, where natural regeneration failure is a main concern for forest managers. For this purpose we periodically collected data from seed traps and trees located at two core locations across several years. We found significant variation in interannual cone production, with the best seed trees being the same across years. In addition, we found highly skewed distributions of female reproductive effort and large fertility differences across stands located few kilometres away. Annual seed dispersal kernels fitted lognormal or 2Dt models depending on the stand analysed, with median dispersal distances between 14 and 25 m. Kernels fitted for maximum dispersal periods showed an outstanding intraseasonal variation of median dispersal distances, from 10 to 54 m, in association to variable patterns of rainfall and maximum wind speed. The amount of seed produced appeared to be enough to guarantee the natural regeneration of the stands during the typical 20-year regeneration period. Colonisation simulations concluded that Mediterranean maritime pine has a notable dispersion capacity, which is strongly influenced by levels of fecundity and, especially, by the number and frequency of long-distance dispersal events. The latter play a key role in tree dispersion processes through enlarging the occupied area and fostering the invasion of abandoned crop land.

Asymptotic analysis of vertical geothermal boreholes in the limit of slowly varyng heat injection rates

Theoretical models for the thermal response of vertical geothermal boreholes often assume that the characteristic time of variation of the heat injection rate is much larger than the characteristic diffusion time across the borehole. In this case, heat transfer inside the borehole and in its immediate surroundings is quasi-steady in the first approximation, while unsteady effects enter only in the far field. Previous studies have exploited this disparity of time scales, incorporating approximate matching conditions to couple the near-borehole region with the outer unsteady temperatura field. In the present work matched asymptotic expansion techniques are used to analyze the heat transfer problem, delivering a rigorous derivation of the true matching condition between the two regions and of the correct definition of the network of thermal resistances that represents the quasi-steady solution near the borehole. Additionally, an apparent temperature due to the unsteady far field is identified that needs to be taken into account by the near-borehole region for the correct computation of the heat injection rate. This temperature differs from the usual mean borehole temperature employed in the literatura.

Monte Carlo limit cycle characterization

The fixed point implementation of IIR digital filters usually leads to the appearance of zero-input limit cycles, which degrade the performance of the system. In this paper, we develop an efficient Monte Carlo algorithm to detect and characterize limit cycles in fixed-point IIR digital filters. The proposed approach considers filters formulated in the state space and is valid for any fixed point representation and quantization function. Numerical simulations on several high-order filters, where an exhaustive search is unfeasible, show the effectiveness of the proposed approach.

Limit load instabilities in structural elements

Este trabajo analiza distintas inestabilidades en estructuras formadas por distintos materiales. En particular, se capturan y se modelan las inestabilidades usando el método de Riks. Inicialmente, se analiza la bifurcación en depósitos cilíndricos formados por material anisótropo sometidos a carga axial y presión interna. El análisis de bifurcación y post-bifurcación asociados con cilindros de pared gruesa se formula para un material incompresible reforzado con dos fibras que son mecánicamente equivalentes y están dispuestas simétricamente. Consideramos dos casos en la naturaleza de la anisotropía: (i) Fibras refuerzo que tienen una influencia particular sobre la respuesta a cortante del material y (ii) Fibras refuerzo que influyen sólo si la fibra cambia de longitud con la deformación. Se analiza la propagación de las inestabilidades. En concreto, se diferencia en el abultamiento (bulging) entre la propagación axial y la propagación radial de la inestabilidad. Distintos modelos sufren una u otra propagación. Por último, distintas inestabilidades asociadas al mecanismo de ablandamiento del material (material softening) en contraposición al de endurecimiento (hardening) en una estructura (viga) de a: hormigón y b: hormigón reforzado son modeladas utilizando una metodología paralela a la desarrollada en el análisis de inestabilidades en tubos sometidos a presión interna. This present work deals with the instability of structures made of various materials. It captures and models different types of instabilities using numerical analysis. Firstly, we consider bifurcation for anisotropic cylindrical shells subject to axial loading and internal pressure. Analysis of bifurcation and post bifurcation of inflated hyperelastic thick-walled cylinder is formulated using a numerical procedure based on the modified Riks method for an incompressible material with two preferred directions which are mechanically equivalent and are symmetrically disposed. Secondly, bulging/necking motion in doubly fiber-reinforced incompressible nonlinearly elastic cylindrical shells is captured and we consider two cases for the nature of the anisotropy: (i) reinforcing models that have a particular influence on the shear response of the material and (ii) reinforcing models that depend only on the stretch in the fiber direction. The different instability motions are considered. Axial propagation of the bulging instability mode in thin-walled cylinders under inflation is analyzed. We present the analytical solution for this particular motion as well as for radial expansion during bulging evolution. For illustration, cylinders that are made of either isotropic incompressible non-linearly elastic materials or doubly fiber reinforced incompressible non-linearly elastic materials are considered. Finally, strain-softening constitutive models are considered to analyze two concrete structures: a reinforced concrete beam and an unreinforced notch beam. The bifurcation point is captured using the Riks method used previously to analyze bifurcation of a pressurized cylinder.

An experimental study of the breakage of liquid bridges at stability limit of minimum volume

An experimental apparatus to study the breaking process of axisymmetric liquid bridges has been developed, and the breaking sequences of a large number of liquid bridge configurations at minimum-volume stability limit have been analyzed. Experimental results show that very close to the breaking moment the neck radius of the liquid bridge varies as t1/3, where t is the time to breakage, irrespective of the value of the distance between the solid disks that support the liquid column.

The extended Betz–Lanchester limit

A proposal for an extended formulation of the power coefficient of a wind turbine is presented. This new formulation is a generalization of the Betz–Lanchester expression for the power coefficient as function of the axial deceleration of the wind speed provoked by the wind turbine in operation. The extended power coefficient takes into account the benefits of the power produced and the cost associated to the production of this energy. By the simple model proposed is evidenced that the purely energetic optimum operation condition giving rise to the Betz–Lanchester limit (maximum energy produced) does not coincide with the global optimum operational condition (maximum benefit generated) if cost of energy and degradation of the wind turbine during operation is considered. The new extended power coefficient is a general parameter useful to define global optimum operation conditions for wind turbines, considering not only the energy production but also the maintenance cost and the economic cost associated to the life reduction of the machine.