Search for resonant diboson production in the lvjj decay channels with the ATLAS detector at 7 TeV

A search for resonant diboson production using a data sample corresponding to 4.7 fb(-1) of integrated luminosity collected by the ATLAS experiment at the Large Hadron Collider in pp collisions at root s = 7 TeV is presented. The search for a narrow resonance in the WW or WZ mass distribution is conducted in a final state with an electron or a muon, missing transverse momentum, and at least two jets. No significant excess is observed and limits are set using three benchmark models: WW resonance masses below 940 and 710 GeV are excluded at 95% confidence level for spin-2 Randall-Sundrum and bulk Randall-Sundrum gravitons, respectively; WZ resonance masses below 950 GeV are excluded at 95% confidence level for a spin-1 extended gauge model W' boson.

Search for WZ resonances in the fully leptonic channel using pp collisions at √s = 8 TeV with the ATLAS detector

A search for resonant WZ production in the ℓνℓ′ℓ′ℓνℓ′ℓ′ (ℓ,ℓ′=e,μℓ,ℓ′=e,μ) decay channel using 20.3 fb−1 of View the MathML sources=8 TeVpp collision data collected by the ATLAS experiment at LHC is presented. No significant deviation from the Standard Model prediction is observed and upper limits on the production cross sections of WZ resonances from an extended gauge model W′W′ and from a simplified model of heavy vector triplets are derived. A corresponding observed (expected) lower mass limit of 1.52 (1.49) TeV is derived for the W′W′ at the 95% confidence level.

Bayesian Adaptive Designs for Early Phase Clinical Trials

My dissertation focuses mainly on Bayesian adaptive designs for phase I and phase II clinical trials. It includes three specific topics: (1) proposing a novel two-dimensional dose-finding algorithm for biological agents, (2) developing Bayesian adaptive screening designs to provide more efficient and ethical clinical trials, and (3) incorporating missing late-onset responses to make an early stopping decision. Treating patients with novel biological agents is becoming a leading trend in oncology. Unlike cytotoxic agents, for which toxicity and efficacy monotonically increase with dose, biological agents may exhibit non-monotonic patterns in their dose-response relationships. Using a trial with two biological agents as an example, we propose a phase I/II trial design to identify the biologically optimal dose combination (BODC), which is defined as the dose combination of the two agents with the highest efficacy and tolerable toxicity. A change-point model is used to reflect the fact that the dose-toxicity surface of the combinational agents may plateau at higher dose levels, and a flexible logistic model is proposed to accommodate the possible non-monotonic pattern for the dose-efficacy relationship. During the trial, we continuously update the posterior estimates of toxicity and efficacy and assign patients to the most appropriate dose combination. We propose a novel dose-finding algorithm to encourage sufficient exploration of untried dose combinations in the two-dimensional space. Extensive simulation studies show that the proposed design has desirable operating characteristics in identifying the BODC under various patterns of dose-toxicity and dose-efficacy relationships. Trials of combination therapies for the treatment of cancer are playing an increasingly important role in the battle against this disease. To more efficiently handle the large number of combination therapies that must be tested, we propose a novel Bayesian phase II adaptive screening design to simultaneously select among possible treatment combinations involving multiple agents. Our design is based on formulating the selection procedure as a Bayesian hypothesis testing problem in which the superiority of each treatment combination is equated to a single hypothesis. During the trial conduct, we use the current values of the posterior probabilities of all hypotheses to adaptively allocate patients to treatment combinations. Simulation studies show that the proposed design substantially outperforms the conventional multi-arm balanced factorial trial design. The proposed design yields a significantly higher probability for selecting the best treatment while at the same time allocating substantially more patients to efficacious treatments. The proposed design is most appropriate for the trials combining multiple agents and screening out the efficacious combination to be further investigated. The proposed Bayesian adaptive phase II screening design substantially outperformed the conventional complete factorial design. Our design allocates more patients to better treatments while at the same time providing higher power to identify the best treatment at the end of the trial. Phase II trial studies usually are single-arm trials which are conducted to test the efficacy of experimental agents and decide whether agents are promising to be sent to phase III trials. Interim monitoring is employed to stop the trial early for futility to avoid assigning unacceptable number of patients to inferior treatments. We propose a Bayesian single-arm phase II design with continuous monitoring for estimating the response rate of the experimental drug. To address the issue of late-onset responses, we use a piece-wise exponential model to estimate the hazard function of time to response data and handle the missing responses using the multiple imputation approach. We evaluate the operating characteristics of the proposed method through extensive simulation studies. We show that the proposed method reduces the total length of the trial duration and yields desirable operating characteristics for different physician-specified lower bounds of response rate with different true response rates.

Growth rates and stable carbon and oxygen isotope record of corals from the Red Sea

Monthly delta18O records of 2 coral colonies (Porites cf. lutea and P. cf. nodifera) from different localities (Aqaba and Eilat) from the northern Gulf of Aqaba, Red Sea, were calibrated with recorded sea surface temperatures (SST) between 1988 and 2000. The results show high correlation coefficients between SST and delta18O. Seasonal variations of coral delta18O in both locations could explain 91% of the recorded SST. Different delta18O/SST relations from both colonies and from the same colonies were obtained, indicating that delta18O from coral skeletons were subject to an extension rate effect. Significant delta18O depletions are associated with high extension rates and higher values with low extension rates. The relation between coral skeletal delta18O and extension rate is not linear and can be described by a simple exponential model. An inverse relationship extends over extension rates from 1 to 5 mm/yr, while for more rapidly growing corals and portions of colonies the relation is constant and the extension rate does not appear to have a significant effect. We recommend that delta18O values be obtained from fast-growing corals or from portions in which the isotopic disequilibrium is fairly constant (extension rate >5 mm/yr). The results show that interspecific differences in corals may produce a significant delta18O profile offset between 2 colonies that is independent of environmental and extension-rate effects. We conclude that the rate of skeletal extension and the species of coral involved have an important influence on coral delta18O and must be considered when using delta18O records for paleoclimatic reconstructions.

Modelización y generación de mapas del comportamiento espacial de Sporisorium reilianum ((Kühn) Langdon y Fullerton) en Estado de México, México

El presente estudio se realizó con la finalidad de modelizar la distribución espacial del carbón de la espiga del maíz causada por Sporisorium reilianum durante 2006 en el Estado de México y su visualización a través de la generación mapas de densidad. El muestreo se realizó en 100 parcelas georreferenciadas por cada localidad analizada. La incidencia de la enfermedad (porcentaje de plantas enfermas) se determinó al establecer cinco puntos parcela, en cada punto se contabilizaron 100 plantas. Se realizó el análisis geoestadístico para estimar el semivariograma experimental, una vez obtenido, se ajustó a un modelo teórico (esférico, exponencial o gaussiano) a través de los programas Variowin 2.2., su ajuste se validó a través de la validación cruzada. Posteriormente, se elaboraron mapas de agregación de la enfermedad con el método de interpolación geoestadística o krigeado. Los resultados indicaron que la enfermedad se presentó en 20 localidades de 19 municipios del Estado de México; todas las localidades presentaron un comportamiento espacial agregado de la enfermedad, 16 localidades se ajustaron al modelo esférico, dos al modelo exponencial y dos localidades se ajustaron al modelo gaussiano. En todos los modelos se lograron establecer mapas de agregación que permitirá adecuar las acciones de manejo en términos de puntos o sitios específicos.

Production sharing, demand spillovers and CO2 emissions : the case of Chinese regions in GVCs

This study adopts the perspective of demand spillovers to provide new insights regarding Chinese domestic-regions' production position in global value chains and their associated CO2 emissions. To this end, we constructed a new type of World Input-Output Database in which China's domestic interregional input-output table for 2007 is endogenously embedded. Then, the pattern of China's regional demand spillovers across both domestic regions and countries are revealed by employing this new database. These results were further connected to endowments theory, which help to make sense of the empirical results. It is found that China's regions locate relatively upstream in GVCs, and had CO2 emissions in net exports, which were entirely predicted by the environmental extended HOV model. Our study points to micro policy instruments to combat climate change, for example, the tax reform for energy inputs that helps to change the production pattern thus has impact on trade pattern and so forth.

Amplificador de señal analógico basado en conceptos biomiméticos. Desarrollo de modelos computacionales y validación experimental

La presente tesis doctoral estudia las características de un transductor acústico bioinspirado en la estructura del maxilar inferior de un Zifio de Couvier (Ziphius cavirostris). El mecanismo de funcionamiento del sensor se basa en las características de un sistema acoplado formado por los distintos componentes acústicos identificados en el maxilar. Para analizar las características del sensor se propone un modelo simplificado 2D que consta de una cavidad cerrada con forma de bocina acoplada a una lengüeta. Una parte de la lengüeta se encuentra dentro de la cavidad y otra en el exterior. Dicha lengüeta detecta los cambios de presión acústica y las vibraciones generadas por el sonido y las transmite con ondas de flexión al interior de la cavidad. La excitación prolongada sobre la placa puede provocar la activación de los modos propios del sistema acoplado. Dichos modos se caracterizan porque presentan un máximo de presión en el cono de la bocina la cual a su vez actúa como un amplificador acústico. Mediante el Método de los elementos Finitos se analizan las características acústicas del sensor y se construye un prototipo experimental para validar los resultados evaluados en el modelo numérico. Se propone una metodología numérica que permite desarrollar y validar un elemento tetraédrico para caracterizar el comportamiento isotrópico de los medios porosos. La metodología permite construir elementos de línea y bidimensionales. A partir de esta metodología se desarrolla un elemento plano ortotrópico. Se realiza un modelo de la cavidad en el que una de las paredes de la bocina está constituida por material poroso y se une la parte exterior de la lengüeta a dicho material para que constituya una viga sobre un apoyo elástico. Se calcula la respuesta modal y se discuten los efectos del material poroso en la eficiencia del transductor y las posibles mejoras a introducir en el mismo. SUMMARY This Thesis studies the characteristics of an acoustic transducer bioinspired by the structure of the lower maxilla of an odontoceto. In this case a Cuvier’s beaked whale (Ziphius cavirostris). The transducer working mechanism is based in a coupled system, with components identified in the maxilla. To analyze the transducer a simplified 2D model composed by a horn shaped closed cavity is modeled. The cavity is coupled with a flat belt. The belt has one part inside the cavity and the other part outside of it. The belt traverses the cavity wall and it is in charge to pick the vibrations from outside and introduce it inside the cavity. The transmission is obtained through the belt bending. A sustained external load with the right frequency contents will allow the system to reach a stationary pressure intensity distribution inside the cavity. Frequencies with modal shapes that show an important intensity increase at the horn tip are of special interest because of the signal amplification. A finite element model is constructed to study the transducer coupled modes and a prototype is constructed to validate the numerical results. A numerical methodology to construct and validate a tetrahedral finite element for isotropic porous materials is presented. The methodology allows constructing linear and 2D elements. It is extended to model orthotropic porous materials behavior. At the end, one of the horn walls is made of an orthotropic material and the external belt is glued to it in order to configure a belt over an elastic foundation. Modal response is evaluated and the porous material effects in the transducer efficiency and further improvements are discussed.

Modelización de un proceso de biorremediación de suelos contaminados con gasoil

El objeto de esta tesis doctoral es encontrar, mediante herramientas estadísticas, un modelo matemático que prediga la evolución temporal de las concentraciones de los hidrocarburos totales del petróleo (TPH) en un proceso de biorremediación de suelos contaminados con gasoleo. Obtenido el modelo matemático realizar las comparaciones pertinentes para determinar la eficiencia de distintos tratamientos, la influencia tanto del suelo como del nivel de concentración. Para cumplir este objetivo se realizó un diseño de experimentos que tomaba en cuenta dos tipos de suelos, dos niveles de concentración de hidrocarburo y seis tratamientos. Se realizaron en total 1824 ensayos en laboratorio repartidos en 8 campañas de campo durante un año de experiencia tanto en parcelas de 4 x 4 m como en microcosmos. Los resultados obtenidos muestran que el proceso de biorremediación se puede ajustar a un comportamiento exponencial, que es una ecuación de cinética de primer orden, y que las comparaciones realizadas han mostrado resultados satisfactorios de la eficiencia de algunos tratamientos. ABSTRACT The purpose of this thesis is using statistical tools find a mathematical model to predict the time evolution of the total petroleum hydrocarbons (TPH) concentrations in the bioremediation process of diesel contaminated soil. Retrieved mathematical model make relevant comparisons to determine the efficiency of different treatments and the influence of both soil and concentration levels. To achieve this goal a design of experiments was developed, it took into account two types of soil, two concentration levels of hydrocarbon and six treatments. There were a total of 1824 laboratory tests spread over 8 field campaigns during a year of experience in both plots of 4 x 4 m in microcosm. The results obtained show that the bioremediation process can be adjusted to an exponential model, it is an equation of kinetic of first order and that the comparisons have shown satisfactory results of the efficiency of some treatments.

Service innovation: Inward and outward related activities and cooperation mode

Knowledge of how customers co-create value, the way that suppliers and providers co-produce services, and how research and development centers and universities transfer technologies is becoming increasingly important to scholars' understanding of service innovation. This paper presents an analysis of the relationship between inward and outward innovation activities in service organizations and their modes of innovation, using network innovation premises and an extended innovation model. Empirical data from retail, health and education sector service organizations show the existence of a relationship between the degree of development of the inward innovation process and the degree of development of outward innovation activities. The majority of service organizations have innovation processes with an orientation toward customers and suppliers rather than other service network members, and leading service organizations follow a path that the literature defines as oriented toward the service value network. Findings lead to implications of how innovation managers could develop their internal innovation capacity to balance inward and outward activities properly.

Modelo gráfico para la datación de vías romanas empedradas a partir del estudio de sus estados de frecuentación y del análisis superficial de roderas

La presente tesis doctoral se orienta al estudio y análisis de los caminos empedrados antiguos, desde la época prerromana, tanto desde el punto de vista histórico como desde el técnico. La cuantificación de la romanidad de un camino representa un objetivo importante para la mayoría de los estudiosos de la caminería antigua, así como para los arqueólogos, por los datos que ofrece acerca del uso del territorio, los trazados de caminos en la antigüedad y los tráficos asociados. Cuantificar la romanidad de un camino no es tarea sencilla debido a que intervienen multitud de condicionantes que están vivos y son cambiantes como consecuencia del dinamismo inherente al propio camino. En cuanto al aspecto histórico, se realiza una descripción y análisis de la evolución del camino en la Península Ibérica desde sus orígenes hasta mediados del siglo XX, que permite diferenciar la red itineraria según su momento histórico. Así mismo, se describen y analizan: las ruedas y los carros desde sus orígenes, especialmente en la época romana -incluyendo una toma de medidas de distintos tipos de carro, existentes en instituciones y colecciones particulares-; las técnicas de transporte en la antigüedad y las características de la infraestructura viaria de época romana, detallando aspectos generales de sus técnicas de ingeniería y construcción. Desde el punto de vista técnico, el enfoque metodológico ha sido definir un Índice de Romanidad del Camino (IRC) para la datación de vías romanas empedradas, basado en un análisis multicriterio, a partir de los distintos factores que caracterizan su romanidad. Se ha realizado un exhaustivo estudio de campo, con la correspondiente toma de datos en las vías. Se han realizado una serie de ensayos de laboratorio con un prototipo creado exprofeso para simular el desgaste de la piedra producido por el traqueteo del carro al circular por el camino empedrado y dar una hipótesis de datación del camino. Se ha realizado un tratamiento estadístico con la muestra de datos medidos en campo. Se ha definido además el concepto de elasticidad de rodera usando la noción de derivada elástica. En cuanto a los resultados obtenidos: se ha calculado el Índice de Romanidad del Camino (IRC) en una serie de vías empedradas, para cuantificar su romanidad, obteniéndose un resultado coherente con la hipótesis previa sobre la datación de dichas vías; y se ha formulado un modelo exponencial para el número de frecuentaciones de carga que lo relaciona con la elasticidad de rodera y con su esbeltez y que se ha utilizado para relacionar la elasticidad de la rodera con la geología de la roca. Se ha iniciado una línea de investigación sobre la estimación de tráficos históricos en la caminería antigua, considerando que el volumen de tráfico a lo largo del tiempo en un tramo de vía está relacionado con los valores de elasticidad de rodera de dicho tramo a través de la tipología de la roca. En resumen, la presente tesis doctoral proporciona un método para sistematizar el estudio de los caminos antiguos, así como para datarlos y estimar la evolución de sus tráficos. The present Ph. D. Thesis aims to study and to analyze ancient cobbled ways, since pre-roman times, both from the historical and technical points of view. The quantification of the Roman character of a way represents an important target for most of the researchers of ancient ways, as well as for the archaeologists, due to information that it offers about the use of the territory, the tracings of ways in the antiquity and the associate flows. To quantify the Roman character of a way is not a simple task because it involves multitude of influent factors that are alive and variable as a result of the dynamism inherent to the way. As for the historical aspect, a description and analysis of the evolution of the way in the Iberian Peninsula from its origins until the middle of the twentieth century has been done. This allows us to distinguish between elements of the network according to its historical moment. Likewise, a description and analysis is given about: the wheels and the cars since their origins, especially in the Roman time - including a capture of measurements of different types of car, belonging to institutions and to particular collections-; the transport techniques on the antiquity and the characteristics of the road infrastructure of Roman epoch, detailing general technical engineering and constructive aspects. From the technical point of view, the methodological approach has been to define an Index of the Roman Character of the Way (IRC) for the dating of cobbled Roman routes, based on a multi-criterion analysis, involving different factors typical of Roman ways. An exhaustive field study has been realized, with the corresponding capture of information in the routes. A series of laboratory essays has been realized with an ad hoc prototype created to simulate the wear of the stone produced by cars circulating along the cobbled way, and to give a dating hypothesis of the way. A statistical treatment has been realized with the sample of information measured in field. There has been defined also the concept of elasticity of rolling trace using the notion of elastic derivative. As for the obtained results: there has been calculated the Index of Roman Character of the Way (IRC) in a series of cobbled routes, to quantify its Roman character, obtaining a coherent result with the previous dating hypothesis of the above mentioned routes; and an exponential model has been formulated for the number of frequent attendances of load that relates this number to the elasticity of rolling trace and to its slenderness and that has been used to relate the elasticity of the rolling trace to the geology of the rock. An investigation line has been opened about the estimation of historical flows in ancient ways, considering that the traffic volume over the course of time in a route stretch is related to the values of elasticity of rolling trace of the above mentioned stretch by means of the typology of the rock. In short, the present Ph. D. Thesis provides a method to systematize the study of ancient ways, as well as to date them and to estimate the evolution of their flows.

Carbon dioxide emissions from semi-arid soils amended with biochar alone or combined with mineral and organic fertilizers

Semi-arid soils cover a significant area of Earth s land surface and typically contain large amounts of inorganic C. Determining the effects of biochar additions on CO2 emissions fromsemi-arid soils is therefore essential for evaluating the potential of biochar as a climate change mitigation strategy. Here, we measured the CO2 that evolved from semi-arid calcareous soils amended with biochar at rates of 0 and 20 t ha?1 in a full factorial combination with three different fertilizers (mineral fertilizer, municipal solid waste compost, and sewage sludge) applied at four rates (equivalent to 0, 75, 150, and 225 kg potentially available N ha?1) during 182 days of aerobic incubation. A double exponential model, which describes cumulative CO2 emissions from two active soil C compartments with different turnover rates (one relatively stable and the other more labile), was found to fit verywell all the experimental datasets. In general, the organic fertilizers increased the size and decomposition rate of the stable and labile soil C pools. In contrast, biochar addition had no effects on any of the double exponential model parameters and did not interact with the effects ascribed to the type and rate of fertilizer. After 182 days of incubation, soil organic and microbial biomass C contents tended to increase with increasing the application rates of organic fertilizer, especially of compost, whereas increasing the rate of mineral fertilizer tended to suppress microbial biomass. Biochar was found to increase both organic and inorganic C contents in soil and not to interactwith the effects of type and rate of fertilizer on C fractions. As a whole, our results suggest that the use of biochar as enhancer of semi-arid soils, either alone or combined with mineral and organic fertilizers, is unlikely to increase abiotic and biotic soil CO2 emissions.

Optimization of the Integrated Guidance and Control for a Dual Aerodynamic Control Missile

Las futuras misiones para misiles aire-aire operando dentro de la atmósfera requieren la interceptación de blancos a mayores velocidades y más maniobrables, incluyendo los esperados vehículos aéreos de combate no tripulados. La intercepción tiene que lograrse desde cualquier ángulo de lanzamiento. Una de las principales discusiones en la tecnología de misiles en la actualidad es cómo satisfacer estos nuevos requisitos incrementando la capacidad de maniobra del misil y en paralelo, a través de mejoras en los métodos de guiado y control modernos. Esta Tesis aborda estos dos objetivos simultáneamente, al proponer un diseño integrando el guiado y el control de vuelo (autopiloto) y aplicarlo a misiles con control aerodinámico simultáneo en canard y cola. Un primer avance de los resultados obtenidos ha sido publicado recientemente en el Journal of Aerospace Engineering, en Abril de 2015, [Ibarrondo y Sanz-Aranguez, 2015]. El valor del diseño integrado obtenido es que permite al misil cumplir con los requisitos operacionales mencionados empleando únicamente control aerodinámico. El diseño propuesto se compara favorablemente con esquemas más tradicionales, consiguiendo menores distancias de paso al blanco y necesitando de menores esfuerzos de control incluso en presencia de ruidos. En esta Tesis se demostrará cómo la introducción del doble mando, donde tanto el canard como las aletas de cola son móviles, puede mejorar las actuaciones de un misil existente. Comparado con un misil con control en cola, el doble control requiere sólo introducir dos servos adicionales para accionar los canards también en guiñada y cabeceo. La sección de cola será responsable de controlar el misil en balanceo mediante deflexiones diferenciales de los controles. En el caso del doble mando, la complicación añadida es que los vórtices desprendidos de los canards se propagan corriente abajo y pueden incidir sobre las superficies de cola, alterando sus características de control. Como un primer aporte, se ha desarrollado un modelo analítico completo para la aerodinámica no lineal de un misil con doble control, incluyendo la caracterización de este efecto de acoplamiento aerodinámico. Hay dos modos de funcionamiento en picado y guiñada para un misil de doble mando: ”desviación” y ”opuesto”. En modo ”desviación”, los controles actúan en la misma dirección, generando un cambio inmediato en la sustentación y produciendo un movimiento de translación en el misil. La respuesta es rápida, pero en el modo ”desviación” los misiles con doble control pueden tener dificultades para alcanzar grandes ángulos de ataque y altas aceleraciones laterales. Cuando los controles actúan en direcciones opuestas, el misil rota y el ángulo de ataque del fuselaje se incrementa para generar mayores aceleraciones en estado estacionario, aunque el tiempo de respuesta es mayor. Con el modelo aerodinámico completo, es posible obtener una parametrización dependiente de los estados de la dinámica de corto periodo del misil. Debido al efecto de acoplamiento entre los controles, la respuesta en bucle abierto no depende linealmente de los controles. El autopiloto se optimiza para obtener la maniobra requerida por la ley de guiado sin exceder ninguno de los límites aerodinámicos o mecánicos del misil. Una segunda contribución de la tesis es el desarrollo de un autopiloto con múltiples entradas de control y que integra la aerodinámica no lineal, controlando los tres canales de picado, guiñada y cabeceo de forma simultánea. Las ganancias del autopiloto dependen de los estados del misil y se calculan a cada paso de integración mediante la resolución de una ecuación de Riccati de orden 21x21. Las ganancias obtenidas son sub-óptimas, debido a que una solución completa de la ecuación de Hamilton-Jacobi-Bellman no puede obtenerse de manera práctica, y se asumen ciertas simplificaciones. Se incorpora asimismo un mecanismo que permite acelerar la respuesta en caso necesario. Como parte del autopiloto, se define una estrategia para repartir el esfuerzo de control entre el canard y la cola. Esto se consigue mediante un controlador aumentado situado antes del bucle de optimización, que minimiza el esfuerzo total de control para maniobrar. Esta ley de alimentación directa mantiene al misil cerca de sus condiciones de equilibrio, garantizando una respuesta transitoria adecuada. El controlador no lineal elimina la respuesta de fase no-mínima característica de la cola. En esta Tesis se consideran dos diseños para el guiado y control, el control en Doble-Lazo y el control Integrado. En la aproximación de Doble-Lazo, el autopiloto se sitúa dentro de un bucle interior y se diseña independientemente del guiado, que conforma el bucle más exterior del control. Esta estructura asume que existe separación espectral entre los dos, esto es, que los tiempos de respuesta del autopiloto son mucho mayores que los tiempos característicos del guiado. En el estudio se combina el autopiloto desarrollado con una ley de guiado óptimo. Los resultados obtenidos demuestran que se consiguen aumentos muy importantes en las actuaciones frente a misiles con control canard o control en cola, y que la interceptación, cuando se lanza cerca del curso de colisión, se consigue desde cualquier ángulo alrededor del blanco. Para el misil de doble mando, la estrategia óptima resulta en utilizar el modo de control opuesto en la aproximación al blanco y utilizar el modo de desviación justo antes del impacto. Sin embargo la lógica de doble bucle no consigue el impacto cuando hay desviaciones importantes con respecto al curso de colisión. Una de las razones es que parte de la demanda de guiado se pierde, ya que el misil solo es capaz de modificar su aceleración lateral, y no tiene control sobre su aceleración axial, a no ser que incorpore un motor de empuje regulable. La hipótesis de separación mencionada, y que constituye la base del Doble-Bucle, puede no ser aplicable cuando la dinámica del misil es muy alta en las proximidades del blanco. Si se combinan el guiado y el autopiloto en un único bucle, la información de los estados del misil está disponible para el cálculo de la ley de guiado, y puede calcularse la estrategia optima de guiado considerando las capacidades y la actitud del misil. Una tercera contribución de la Tesis es la resolución de este segundo diseño, la integración no lineal del guiado y del autopiloto (IGA) para el misil de doble control. Aproximaciones anteriores en la literatura han planteado este sistema en ejes cuerpo, resultando en un sistema muy inestable debido al bajo amortiguamiento del misil en cabeceo y guiñada. Las simplificaciones que se tomaron también causan que el misil se deslice alrededor del blanco y no consiga la intercepción. En nuestra aproximación el problema se plantea en ejes inerciales y se recurre a la dinámica de los cuaterniones, eliminado estos inconvenientes. No se limita a la dinámica de corto periodo del misil, porque se construye incluyendo de modo explícito la velocidad dentro del bucle de optimización. La formulación resultante en el IGA es independiente de la maniobra del blanco, que sin embargo se ha de incluir en el cálculo del modelo en Doble-bucle. Un típico inconveniente de los sistemas integrados con controlador proporcional, es el problema de las escalas. Los errores de guiado dominan sobre los errores de posición del misil y saturan el controlador, provocando la pérdida del misil. Este problema se ha tratado aquí con un controlador aumentado previo al bucle de optimización, que define un estado de equilibrio local para el sistema integrado, que pasa a actuar como un regulador. Los criterios de actuaciones para el IGA son los mismos que para el sistema de Doble-Bucle. Sin embargo el problema matemático resultante es muy complejo. El problema óptimo para tiempo finito resulta en una ecuación diferencial de Riccati con condiciones terminales, que no puede resolverse. Mediante un cambio de variable y la introducción de una matriz de transición, este problema se transforma en una ecuación diferencial de Lyapunov que puede resolverse mediante métodos numéricos. La solución resultante solo es aplicable en un entorno cercano del blanco. Cuando la distancia entre misil y blanco es mayor, se desarrolla una solución aproximada basada en la solución de una ecuación algebraica de Riccati para cada paso de integración. Los resultados que se han obtenido demuestran, a través de análisis numéricos en distintos escenarios, que la solución integrada es mejor que el sistema de Doble-Bucle. Las trayectorias resultantes son muy distintas. El IGA preserva el guiado del misil y consigue maximizar el uso de la propulsión, consiguiendo la interceptación del blanco en menores tiempos de vuelo. El sistema es capaz de lograr el impacto donde el Doble-Bucle falla, y además requiere un orden menos de magnitud en la cantidad de cálculos necesarios. El efecto de los ruidos radar, datos discretos y errores del radomo se investigan. El IGA es más robusto, resultando menos afectado por perturbaciones que el Doble- Bucle, especialmente porque el núcleo de optimización en el IGA es independiente de la maniobra del blanco. La estimación de la maniobra del blanco es siempre imprecisa y contaminada por ruido, y degrada la precisión de la solución de Doble-Bucle. Finalmente, como una cuarta contribución, se demuestra que el misil con guiado IGA es capaz de realizar una maniobra de defensa contra un blanco que ataque por su cola, sólo con control aerodinámico. Las trayectorias estudiadas consideran una fase pre-programada de alta velocidad de giro, manteniendo siempre el misil dentro de su envuelta de vuelo. Este procedimiento no necesita recurrir a soluciones técnicamente más complejas como el control vectorial del empuje o control por chorro para ejecutar esta maniobra. En todas las demostraciones matemáticas se utiliza el producto de Kronecker como una herramienta practica para manejar las parametrizaciones dependientes de variables, que resultan en matrices de grandes dimensiones. ABSTRACT Future missions for air to air endo-atmospheric missiles require the interception of targets with higher speeds and more maneuverable, including forthcoming unmanned supersonic combat vehicles. The interception will need to be achieved from any angle and off-boresight launch conditions. One of the most significant discussions in missile technology today is how to satisfy these new operational requirements by increasing missile maneuvering capabilities and in parallel, through the development of more advanced guidance and control methods. This Thesis addresses these two objectives by proposing a novel optimal integrated guidance and autopilot design scheme, applicable to more maneuverable missiles with forward and rearward aerodynamic controls. A first insight of these results have been recently published in the Journal of Aerospace Engineering in April 2015, [Ibarrondo and Sanz-Aránguez, 2015]. The value of this integrated solution is that it allows the missile to comply with the aforementioned requirements only by applying aerodynamic control. The proposed design is compared against more traditional guidance and control approaches with positive results, achieving reduced control efforts and lower miss distances with the integrated logic even in the presence of noises. In this Thesis it will be demonstrated how the dual control missile, where canard and tail fins are both movable, can enhance the capabilities of an existing missile airframe. Compared to a tail missile, dual control only requires two additional servos to actuate the canards in pitch and yaw. The tail section will be responsible to maintain the missile stabilized in roll, like in a classic tail missile. The additional complexity is that the vortices shed from the canard propagate downstream where they interact with the tail surfaces, altering the tail expected control characteristics. These aerodynamic phenomena must be properly described, as a preliminary step, with high enough precision for advanced guidance and control studies. As a first contribution we have developed a full analytical model of the nonlinear aerodynamics of a missile with dual control, including the characterization of this cross-control coupling effect. This development has been produced from a theoretical model validated with reliable practical data obtained from wind tunnel experiments available in the scientific literature, complement with computer fluid dynamics and semi-experimental methods. There are two modes of operating a missile with forward and rear controls, ”divert” and ”opposite” modes. In divert mode, controls are deflected in the same direction, generating an increment in direct lift and missile translation. Response is fast, but in this mode, dual control missiles may have difficulties in achieving large angles of attack and high level of lateral accelerations. When controls are deflected in opposite directions (opposite mode) the missile airframe rotates and the body angle of attack is increased to generate greater accelerations in steady-state, although the response time is larger. With the aero-model, a state dependent parametrization of the dual control missile short term dynamics can be obtained. Due to the cross-coupling effect, the open loop dynamics for the dual control missile is not linearly dependent of the fin positions. The short term missile dynamics are blended with the servo system to obtain an extended autopilot model, where the response is linear with the control fins turning rates, that will be the control variables. The flight control loop is optimized to achieve the maneuver required by the guidance law without exceeding any of the missile aerodynamic or mechanical limitations. The specific aero-limitations and relevant performance indicators for the dual control are set as part of the analysis. A second contribution of this Thesis is the development of a step-tracking multi-input autopilot that integrates non-linear aerodynamics. The designed dual control missile autopilot is a full three dimensional autopilot, where roll, pitch and yaw are integrated, calculating command inputs simultaneously. The autopilot control gains are state dependent, and calculated at each integration step solving a matrix Riccati equation of order 21x21. The resulting gains are sub-optimal as a full solution for the Hamilton-Jacobi-Bellman equation cannot be resolved in practical terms and some simplifications are taken. Acceleration mechanisms with an λ-shift is incorporated in the design. As part of the autopilot, a strategy is defined for proper allocation of control effort between canard and tail channels. This is achieved with an augmented feed forward controller that minimizes the total control effort of the missile to maneuver. The feedforward law also maintains the missile near trim conditions, obtaining a well manner response of the missile. The nonlinear controller proves to eliminate the non-minimum phase effect of the tail. Two guidance and control designs have been considered in this Thesis: the Two- Loop and the Integrated approaches. In the Two-Loop approach, the autopilot is placed in an inner loop and designed separately from an outer guidance loop. This structure assumes that spectral separation holds, meaning that the autopilot response times are much higher than the guidance command updates. The developed nonlinear autopilot is linked in the study to an optimal guidance law. Simulations are carried on launching close to collision course against supersonic and highly maneuver targets. Results demonstrate a large boost in performance provided by the dual control versus more traditional canard and tail missiles, where interception with the dual control close to collision course is achieved form 365deg all around the target. It is shown that for the dual control missile the optimal flight strategy results in using opposite control in its approach to target and quick corrections with divert just before impact. However the Two-Loop logic fails to achieve target interception when there are large deviations initially from collision course. One of the reasons is that part of the guidance command is not followed, because the missile is not able to control its axial acceleration without a throttleable engine. Also the separation hypothesis may not be applicable for a high dynamic vehicle like a dual control missile approaching a maneuvering target. If the guidance and autopilot are combined into a single loop, the guidance law will have information of the missile states and could calculate the most optimal approach to the target considering the actual capabilities and attitude of the missile. A third contribution of this Thesis is the resolution of the mentioned second design, the non-linear integrated guidance and autopilot (IGA) problem for the dual control missile. Previous approaches in the literature have posed the problem in body axes, resulting in high unstable behavior due to the low damping of the missile, and have also caused the missile to slide around the target and not actually hitting it. The IGA system is posed here in inertial axes and quaternion dynamics, eliminating these inconveniences. It is not restricted to the missile short term dynamic, and we have explicitly included the missile speed as a state variable. The IGA formulation is also independent of the target maneuver model that is explicitly included in the Two-loop optimal guidance law model. A typical problem of the integrated systems with a proportional control law is the problem of scales. The guidance errors are larger than missile state errors during most of the flight and result in high gains, control saturation and loss of control. It has been addressed here with an integrated feedforward controller that defines a local equilibrium state at each flight point and the controller acts as a regulator to minimize the IGA states excursions versus the defined feedforward state. The performance criteria for the IGA are the same as in the Two-Loop case. However the resulting optimization problem is mathematically very complex. The optimal problem in a finite-time horizon results in an irresoluble state dependent differential Riccati equation with terminal conditions. With a change of variable and the introduction of a transition matrix, the equation is transformed into a time differential Lyapunov equation that can be solved with known numerical methods in real time. This solution results range limited, and applicable when the missile is in a close neighborhood of the target. For larger ranges, an approximate solution is used, obtained from solution of an algebraic matrix Riccati equation at each integration step. The results obtained show, by mean of several comparative numerical tests in diverse homing scenarios, than the integrated approach is a better solution that the Two- Loop scheme. Trajectories obtained are very different in the two cases. The IGA fully preserves the guidance command and it is able to maximize the utilization of the missile propulsion system, achieving interception with lower miss distances and in lower flight times. The IGA can achieve interception against off-boresight targets where the Two- Loop was not able to success. As an additional advantage, the IGA also requires one order of magnitude less calculations than the Two-Loop solution. The effects of radar noises, discrete radar data and radome errors are investigated. IGA solution is robust, and less affected by radar than the Two-Loop, especially because the target maneuvers are not part of the IGA core optimization loop. Estimation of target acceleration is always imprecise and noisy and degrade the performance of the two-Loop solution. The IGA trajectories are such that minimize the impact of radome errors in the guidance loop. Finally, as a fourth contribution, it is demonstrated that the missile with IGA guidance is capable of performing a defense against attacks from its rear hemisphere, as a tail attack, only with aerodynamic control. The studied trajectories have a preprogrammed high rate turn maneuver, maintaining the missile within its controllable envelope. This solution does not recur to more complex features in service today, like vector control of the missile thrust or side thrusters. In all the mathematical treatments and demonstrations, the Kronecker product has been introduced as a practical tool to handle the state dependent parametrizations that have resulted in very high order matrix equations.

Estrategias de fijación de precios de las aerolíneas de bajo coste: una aproximación al modelo de la rivalidad ampliada

En la presente tesis doctoral se realiza una aproximación a las estrategias de fijación de precios de las aerolíneas de bajo coste europeas a partir del estudio de su entorno (demanda, proveedores, rivales, productos sustitutos, etcétera). Para ello se emplea el marco del modelo de las cinco fuerzas descrito por Michael Porter. Para su análisis, se utiliza un método combinado de Análisis Factorial (ACP) y datos de panel.

Modeling of Distributed Generation under Next Generation Interconnection Requirements

Thesis (Master's)--University of Washington, 2016-06

Learning Models of Gene Expression from Synthetic DNA Sequences

Thesis (Ph.D.)--University of Washington, 2016-06