Resilience analysis of networked systems-of-systems based on structural and dynamic interdependencies

Critical infrastructures support everyday activities in modern societies, facilitating the exchange of services and quantities of various nature. Their functioning is the result of the integration of diverse technologies, systems and organizations into a complex network of interconnections. Benefits from networking are accompanied by new threats and risks. In particular, because of the increased interdependency, disturbances and failures may propagate and render unstable the whole infrastructure network. This paper presents a methodology of resilience analysis of networked systems of systems. Resilience generalizes the concept of stability of a system around a state of equilibrium, with respect to a disturbance and its ability of preventing, resisting and recovery. The methodology provides a tool for the analysis of off-equilibrium conditions that may occur in a single system and propagate through the network of dependencies. The analysis is conducted in two stages. The first stage of the analysis is qualitative. It identifies the resilience scenarios, i.e. the sequence of events, triggered by an initial disturbance, which include failures and the system response. The second stage is quantitative. The most critical scenarios can be simulated, for the desired parameter settings, in order to check if they are successfully handled, i.e recovered to nominal conditions, or they end into the network failure. The proposed methodology aims at providing an effective support to resilience-informed design.

Intermediate Band to Conduction Band optical absorption in ZnTe:O

ZnTe doped with high concentrations of oxygen has been proposed in previous works as intermediate band (IB) material for photovoltaic applications. The existence of extra optical transitions related to the presence of an IB has already been demonstrated in this material and it has been possible to measure the absorption coefficient of the transitions from the valence band (VB) to the IB. In this work we present the first measurement of the absorption coefficient associated to transitions from the IB to the conduction band (CB) in ZnTe:O. The samples used are 4 ?m thick ZnTe layers with or without O in a concentration ~ 1019 cm-3, which have been grown on semi-insulating GaAs substrates by molecular beam epitaxy (MBE). The IB-CB absorption coefficient peaks for photon energies ~ 0.4 eV. It is extracted from reflectance and transmittance spectra measured using Fourier Transform Infrared (FTIR) spectroscopy. Under typical FTIR measurement conditions (low light intensity, broadband spectrum) the absorption coefficient in IB-to-CB transitions reaches 700 cm-1. This is much weaker than the one observed for VB-IB absorption. This result is consistent with the fact that the IB is expected to be nearly empty of electrons under equilibrium conditions in ZnTe(:O). The absorption for VB to IB transitions is also observed in the same samples through reflectance measurements performed in the visible range using a monochromator. These measurements are compared with the quantum efficiency (QE) from solar cells fabricated under similar conditions.

The breaking of axisymmetric slender liquid bridges

Liquids held by surface tension forces can bridge the gap between two solid bodies placed not too far apart from each other. The equilibrium conditions and stability criteria for static, cylindrical liquid bridges are well known. However, the behaviour of an unstable liquid bridge, regarding both its transition toward breaking and the resulting configuration, is a matter for discussion. The dynamical problem of axisymmetric rupture of a long liquid bridge anchored at two equal coaxial disks is treated in this paper through the adoption of one-dimensional theories which are widely used in capillary jet problems

Intermediate band to conduction band optical absorption in ZnTeO

ZnTe doped with high concentrations of oxygen has been proposed in previous works as an intermediate band (IB) material for photovoltaic applications. The existence of extra optical transitions related to the presence of an IB has already been demonstrated in this material and it has been possible to measure the absorption coefficient of the transitions from the valence band (VB) to the IB. In this study, we present the first measurement of the absorption coefficient associated with transitions from the IB to the conduction band (CB) in ZnTeO. The samples used are 4-mum-thick ZnTe layers with or without O in a concentration ~10 19 cm -3, which have been grown on semiinsulating GaAs substrates by molecular beam epitaxy (MBE). The IB-CB absorption coefficient peaks for photon energies ~0.4 eV. It is extracted from reflectance and transmittance spectra measured using Fourier transform infrared (FTIR) spectroscopy. Under typical FTIR measurement conditions (low light intensity, broadband spectrum), the absorption coefficient in IB-to-CB transitions reaches 700 cm -1. This is much weaker than the one observed for VB-IB absorption. This result is consistent with the fact that the IB is expected to be nearly empty of electrons under equilibrium conditions in ZnTe(O).

Human Factor in Resilience Systems

This paper shows the importance of a holistic comprehension of the Earth as a living planet, where man inhabits and is exposed to environmental incidences of different nature. The aim of the paper here summarized is a reflection on all these concepts and scientific considerations related to the important role of men in the handling of natural hazards. Our Planet is an unstable and dynamical system highly sensitive to initial conditions, as proposed by Chaos theory (González-Miranda 2004); it is a complex organic whole, which responds to minimal variations which can affect several natural phenomena such as plate tectonics, solar flares, fluid turbulences, landscape formation, forest fires, growth and migration of populations and biological evolution. This is known as the “butterfly effect” (Lorenz 1972), which means that a small change of the system causes a chain of events leading to large-scale unpredictable consequences. The aim of this work is dwelling on the importance of the knowledge of these natural and catastrophic geological, biological and human systems so much sensible to equilibrium conditions, to prevent, avoid and mend their effects, and to face them in a resilient way

Evaluación de la seguridad al deslizamiento de una presa de fábrica en un cimiento con una familia de discontinuidades y con un criterio de rotura con ley de fluencia no asociada

En España existen del orden de 1,300 grandes presas, de las cuales un 20% fueron construidas antes de los años 60. El hecho de que existan actualmente una gran cantidad de presas antiguas aún en operación, ha producido un creciente interés en reevaluar su seguridad empleando herramientas nuevas o modificadas que incorporan modelos de fallo teóricos más completos, conceptos geotécnicos más complejos y nuevas técnicas de evaluación de la seguridad. Una manera muy común de abordar el análisis de estabilidad de presas de gravedad es, por ejemplo, considerar el deslizamiento a través de la interfase presa-cimiento empleando el criterio de rotura lineal de Mohr-Coulomb, en donde la cohesión y el ángulo de rozamiento son los parámetros que definen la resistencia al corte de la superficie de contacto. Sin embargo la influencia de aspectos como la presencia de planos de debilidad en el macizo rocoso de cimentación; la influencia de otros criterios de rotura para la junta y para el macizo rocoso (ej. el criterio de rotura de Hoek-Brown); las deformaciones volumétricas que ocurren durante la deformación plástica en el fallo del macizo rocoso (i.e., influencia de la dilatancia) no son usualmente consideradas durante el diseño original de la presa. En este contexto, en la presente tesis doctoral se propone una metodología analítica para el análisis de la estabilidad al deslizamiento de presas de hormigón, considerando un mecanismo de fallo en la cimentación caracterizado por la presencia de una familia de discontinuidades. En particular, se considera la posibilidad de que exista una junta sub-horizontal, preexistente y persistente en el macizo rocoso de la cimentación, con una superficie potencial de fallo que se extiende a través del macizo rocoso. El coeficiente de seguridad es entonces estimado usando una combinación de las resistencias a lo largo de los planos de rotura, cuyas resistencias son evaluadas empleando los criterios de rotura no lineales de Barton y Choubey (1977) y Barton y Bandis (1990), a lo largo del plano de deslizamiento de la junta; y el criterio de rotura de Hoek y Brown (1980) en su versión generalizada (Hoek et al. 2002), a lo largo del macizo rocoso. La metodología propuesta también considera la influencia del comportamiento del macizo rocoso cuando este sigue una ley de flujo no asociada con ángulo de dilatancia constante (Hoek y Brown 1997). La nueva metodología analítica propuesta es usada para evaluar las condiciones de estabilidad empleando dos modelos: un modelo determinista y un modelo probabilista, cuyos resultados son el valor del coeficiente de seguridad y la probabilidad de fallo al deslizamiento, respectivamente. El modelo determinista, implementado en MATLAB, es validado usando soluciones numéricas calculadas mediante el método de las diferencias finitas, empleando el código FLAC 6.0. El modelo propuesto proporciona resultados que son bastante similares a aquellos calculados con FLAC; sin embargo, los costos computacionales de la formulación propuesta son significativamente menores, facilitando el análisis de sensibilidad de la influencia de los diferentes parámetros de entrada sobre la seguridad de la presa, de cuyos resultados se obtienen los parámetros que más peso tienen en la estabilidad al deslizamiento de la estructura, manifestándose además la influencia de la ley de flujo en la rotura del macizo rocoso. La probabilidad de fallo es obtenida empleando el método de fiabilidad de primer orden (First Order Reliability Method; FORM), y los resultados de FORM son posteriormente validados mediante simulaciones de Monte Carlo. Los resultados obtenidos mediante ambas metodologías demuestran que, para el caso no asociado, los valores de probabilidad de fallo se ajustan de manera satisfactoria a los obtenidos mediante las simulaciones de Monte Carlo. Los resultados del caso asociado no son tan buenos, ya que producen resultados con errores del 0.7% al 66%, en los que no obstante se obtiene una buena concordancia cuando los casos se encuentran en, o cerca de, la situación de equilibrio límite. La eficiencia computacional es la principal ventaja que ofrece el método FORM para el análisis de la estabilidad de presas de hormigón, a diferencia de las simulaciones de Monte Carlo (que requiere de al menos 4 horas por cada ejecución) FORM requiere tan solo de 1 a 3 minutos en cada ejecución. There are 1,300 large dams in Spain, 20% of which were built before 1960. The fact that there are still many old dams in operation has produced an interest of reevaluate their safety using new or updated tools that incorporate state-of-the-art failure modes, geotechnical concepts and new safety assessment techniques. For instance, for gravity dams one common design approach considers the sliding through the dam-foundation interface, using a simple linear Mohr-Coulomb failure criterion with constant friction angle and cohesion parameters. But the influence of aspects such as the persistence of joint sets in the rock mass below the dam foundation; of the influence of others failure criteria proposed for rock joint and rock masses (e.g. the Hoek-Brown criterion); or the volumetric strains that occur during plastic failure of rock masses (i.e., the influence of dilatancy) are often no considered during the original dam design. In this context, an analytical methodology is proposed herein to assess the sliding stability of concrete dams, considering an extended failure mechanism in its rock foundation, which is characterized by the presence of an inclined, and impersistent joint set. In particular, the possibility of a preexisting sub-horizontal and impersistent joint set is considered, with a potential failure surface that could extend through the rock mass; the safety factor is therefore computed using a combination of strength along the rock joint (using the nonlinear Barton and Choubey (1977) and Barton and Bandis (1990) failure criteria) and along the rock mass (using the nonlinear failure criterion of Hoek and Brown (1980) in its generalized expression from Hoek et al. (2002)). The proposed methodology also considers the influence of a non-associative flow rule that has been incorporated using a (constant) dilation angle (Hoek and Brown 1997). The newly proposed analytical methodology is used to assess the dam stability conditions, employing for this purpose the deterministic and probabilistic models, resulting in the sliding safety factor and the probability of failure respectively. The deterministic model, implemented in MATLAB, is validated using numerical solution computed with the finite difference code FLAC 6.0. The proposed deterministic model provides results that are very similar to those computed with FLAC; however, since the new formulation can be implemented in a spreadsheet, the computational cost of the proposed model is significantly smaller, hence allowing to more easily conduct parametric analyses of the influence of the different input parameters on the dam’s safety. Once the model is validated, parametric analyses are conducting using the main parameters that describe the dam’s foundation. From this study, the impact of the more influential parameters on the sliding stability analysis is obtained and the error of considering the flow rule is assessed. The probability of failure is obtained employing the First Order Reliability Method (FORM). The probabilistic model is then validated using the Monte Carlo simulation method. Results obtained using both methodologies show good agreement for cases in which the rock mass has a nonassociate flow rule. For cases with an associated flow rule errors between 0.70% and 66% are obtained, so that the better adjustments are obtained for cases with, or close to, limit equilibrium conditions. The main advantage of FORM on sliding stability analyses of gravity dams is its computational efficiency, so that Monte Carlo simulations require at least 4 hours on each execution, whereas FORM requires only 1 to 3 minutes on each execution.

Determinación del riesgo de incendio en una instalación de secado de alperujo = Determination of the risk of fire in an olive-oil mill waste dryer

For the energy valorization of alperujo, residue of the olive oil two phases extraction process, it is necessary to perform a drying process to reduce moisture content from over 60% to less than 10%. In order to reduce primary energy consumption and get an economic return, usually in this kind of drying facilities Gas Turbine CHP is used as a heat source. There have been recently in Spain some fires in this kind of GT-CHP facilities, which have caused high material losses. In some of these fires it has been suggested that the fire was caused by the output of incandescent alperujo in the flue gasesof the drying system. Therefore, the aim of this study is to determine experimentally and analytically under which operational conditions a process of alperujo self-ignition in the drying process can begin, and determine the actual fire hazard in this type of TG-CHP system. For analytical study, the temperature and initial composition of the combustion gases of the Gas Turbine at the entrance of the drying process was calculated and the gas equilibrium conditions reached in contact with the biomass were calculated and, therefore, the temperature of the biomass during the drying process. Moreover, the layer and dust ignition temperature of alperujo has been experimentally determined, according to EN 50281-2-1: 2000. With these results, the operating conditions of the drying process, in which there are real risk of auto-ignition of alperujo have been established.Para la valorización energética del alperujo, residuo del proceso de extracción en dos fases del aceite de oliva, es necesario realizar un proceso de secado para reducir su contenido de humedad de más del 60% al 10% m/m en b.h. Con el fin de reducir el consumo de energía primaria y obtener una rentabilidad económica, normalmente en este tipo de instalaciones de secado se usa la cogeneración con turbina de gas (TG) como fuente de calor. En España en los últimos años han ocurrido algunos casos de incendio en este tipo de instalaciones de cogeneración, que han supuesto pérdidas materiales muy elevadas. Por esta razón, el objetivo de este trabajo es determinar analítica y experimentalmente las condiciones operativas del secadero bajo las cuales podría comenzar un proceso de autoinflamación del alperujo y determinar el riesgo real de incendio en este tipo de instalaciones. Para el estudio analítico, se ha planteado y validado el modelo matemático que permite calcular la temperatura y la composición de los gases de combustión a la entrada y a la salida del secadero, en función de las curvas características de la TG, de las condiciones atmosféricas, del caudal y del grado de humedad de la biomasa tratada. El modelo permite además calcular la temperatura de bulbo húmedo, que es la máxima temperatura que podría alcanzar la biomasa durante el proceso de secado y determinar la cantidad de biomasa que se puede secar completamente en función del caudal y de las condiciones de entrada de los gases de combustión. Con estos resultados y la temperatura mínima de autoinflamación del alperujo determinada experimentalmente siguiendo la norma EN 50281- 2-1:2000, se demuestra que en un proceso de secado de alperujo en condiciones normales de operación no existe riesgo de autoencendido que pueda dar origen a un incendio.

Contribution to Theme A of the Benchmark Workshop: Effect of Concrete Swelling on the Equilibrium and Displacements of an Arch Dam

The Kariba dam is undergoing concrete expansion as a result of an alkali-aggregate reaction. The model adopted to simulate the process is explained in the paper; it is based on the model first proposed by Ulm et al, as later modified by Saouma and Perotti. It has been implemented in the commercial finite element code Abaqus and applied to solve the benchmark problem. The parameters of the model were calibrated using the data recorded up to 1995. The calibrated model was then used for predicting the evolution of the dam up to the present date. Apart from this prediction the paper offers a number of conclusions, such as the fact that the stress level appears to have a major influence on the expansion process; and it presents some suggestions to improve the formulation of the benchmark, such as providing temperature data and widening the locations and conditions of the data employed in the calibration

Effect of Concrete Swelling on the Equilibrium and Displacements of an Arch Dam

The Kariba dam is undergoing concrete expansion as a result of an alkali-aggregate reaction. The model adopted to simulate the process is explained in the paper; it is based on the model first proposed by Ulm et al, as later modified by Saouma and Perotti. It has been implemented in the commercial finite element code Abaqus and applied to solve the benchmark problem. The parameters of the model were calibrated using the data recorded up to 1995. The calibrated model was then used for predicting the evolution of the dam up to the present date. Apart from this prediction the paper offers a number of conclusions, such as the fact that the stress level appears to have a major influence on the expansion process; and it presents some suggestions to improve the formulation of the benchmark, such as providing temperature data and widening the locations and conditions of the data employed in the calibration

A characterization of conserved quantities in non-equilibrium thermodynamics

The well-known Noether theorem in Lagrangian and Hamiltonian mechanics associates symmetries in the evolution equations of a mechanical system with conserved quantities. In this work, we extend this classical idea to problems of non-equilibrium thermodynamics formulated within the GENERIC (General Equations for Non-Equilibrium Reversible-Irreversible Coupling) framework. The geometric meaning of symmetry is reviewed in this formal setting and then utilized to identify possible conserved quantities and the conditions that guarantee their strict conservation. Examples are provided that demonstrate the validity of the proposed definition in the context of finite and infinite dimensional thermoelastic problems.

Enhanced thermionic currents by non equilibrium electron population of metals

An analytical expression is derived for the electron thermionic current from heated metals by using a non equilibrium, modified Kappa energy distribution for electrons. This isotropic distribution characterizes the long high energy tails in the electron energy spectrum for low values of the index ? and also accounts for the Fermi energy for the metal electrons. The limit for large ? recovers the classical equilibrium Fermi-Dirac distribution. The predicted electron thermionic current for low ? increases between four and five orders of magnitude with respect to the predictions of the equilibrium Richardson-Dushmann current. The observed departures from this classical expression, also recovered for large ?, would correspond to moderate values of this index. The strong increments predicted by the thermionic emission currents suggest that, under appropriate conditions, materials with non equilibrium electron populations would become more efficient electron emitters at low temperatures.

Flight dynamics and stability of kites in steady and unsteady wind conditions

The flight dynamics and stability of a kite with a single main line flying in steady and unsteady wind conditions are discussed. A simple dynamic model with five degrees of freedom is derived with the aid of Lagrangian formulation, which explicitly avoids any constraint force in the equations of motion. The longitudinal and lateral–directional modes and stability of the steady flight under constant wind conditions are analyzed by using both numerical and analytical methods. Taking advantage of the appearance of small dimensionless parameters in the model, useful analytical formulas for stable-designed kites are found. Under nonsteady wind-velocity conditions, the equilibrium state disappears and periodic orbits occur. The kite stability and an interesting resonance phenomenon are explored with the aid of a numerical method based on Floquet theory.