Diseño y Desarrollo de un Robot de Inspección de Tuberías

El siguiente trabajo presenta un desarrollo innovador de un Robot de Inspección de Tuberías, el cual fue diseñado para operación de Tuberías Verticales (RETOV Robot de operación en tuberías Verticales (siglas en Español)). Cuando RETOV es usado para inspeccionar pozos petroleros, puede ser controlado desde la superficie, esto permite al operador la inspección y monitoreo del pozo. Esto es efectuado a través de una interfaz de usuario que controla el Robot. RETOV fue diseñado con anillos articulados que abrazan la tubería, un nuevo sistema de amortiguación y ruedas móviles brinda la posibilidad de navegar y evitar obstáculos en la tubería realizando tres tipos de movimientos: rotación, traslación y helicoidal. Adicionalmente RETOV fue diseñado para realizar inspección entre dos tuberías(la tubería de producción y el pozo petrolero), con un sistema de seguridad y una estructura liviana, RETOV puede ser equipado con sensores para medir las variables de interés en el pozo. El diseño de la mecánica y el programa de control, la instrumentación, los modelos matemáticos y resultados de las pruebas han sido descritos en esta Tesis. The following Doctoral Thesis presents an innovative application of a Pipe Inspection Robot, which was designed especially as Operation Robot in Vertical Pipes (RETOV) (Robot for operation vertical pipes (Spanish acronyms)). When RETOV is used to inspect Oil Well, can be controlled from the surface, it´s allow the operator the well inspection and monitoring. This is done through our interface that controls the robot. The RETOV was designed with an articulated ring systems that embrace the pipe, a new damping systems and mobile wheels bring the possibility to navigate and avoid obstacles into the well, performing three types of movement: Rotational, translational and screw. Additional, RETOV was designed to perform the inspection between two pipes (pipe and oil well production), with a security braking systems and lightweight structure; RETOV can be equipment with sensors to measure the variables of interest in the well. The hardware and Software design, the instrumentation, the mathematical models, and the test results have been described in this thesis.

Mode II fracture energy in the adhesive bonding of dissimilar substrates: carbon fibre composite to aluminium joints

The end-notched flexure (ENF) test calculates the value of mode II fracture energy in adhesive bonding between the substrates of same nature. Traditional methods of calculating fracture energy in the ENF test are not suitable in cases where the thickness of the adhesive is non-negligible compared with adherent thicknesses. To address this issue, a specific methodology for calculating mode II fracture energy has been proposed in this paper. To illustrate the applicability of the proposed method, the fracture energy was calculated by the ENF test for adhesive bonds between aluminium and a composite material, which considered two different types of adhesive (epoxy and polyurethane) and various surface treatments. The proposed calculation model provides higher values of fracture energy than those obtained from the simplified models that consider the adhesive thickness to be zero, supporting the conclusion that the calculation of mode II fracture energy for adhesives with non-negligible thickness relative to their adherents should be based on mathematical models, such as the method proposed in this paper, that incorporate the influence of this thickness.

Computational Methods for Identification of Vibrating Structures

System identification deals with the problem of building mathematical models of dynamical systems based on observed data from the system" [1]. In the context of civil engineering, the system refers to a large scale structure such as a building, bridge, or an offshore structure, and identification mostly involves the determination of modal parameters (the natural frequencies, damping ratios, and mode shapes). This paper presents some modal identification results obtained using a state-of-the-art time domain system identification method (data-driven stochastic subspace algorithms [2]) applied to the output-only data measured in a steel arch bridge. First, a three dimensional finite element model was developed for the numerical analysis of the structure using ANSYS. Modal analysis was carried out and modal parameters were extracted in the frequency range of interest, 0-10 Hz. The results obtained from the finite element modal analysis were used to determine the location of the sensors. After that, ambient vibration tests were conducted during April 23-24, 2009. The response of the structure was measured using eight accelerometers. Two stations of three sensors were formed (triaxial stations). These sensors were held stationary for reference during the test. The two remaining sensors were placed at the different measurement points along the bridge deck, in which only vertical and transversal measurements were conducted (biaxial stations). Point estimate and interval estimate have been carried out in the state space model using these ambient vibration measurements. In the case of parametric models (like state space), the dynamic behaviour of a system is described using mathematical models. Then, mathematical relationships can be established between modal parameters and estimated point parameters (thus, it is common to use experimental modal analysis as a synonym for system identification). Stable modal parameters are found using a stabilization diagram. Furthermore, this paper proposes a method for assessing the precision of estimates of the parameters of state-space models (confidence interval). This approach employs the nonparametric bootstrap procedure [3] and is applied to subspace parameter estimation algorithm. Using bootstrap results, a plot similar to a stabilization diagram is developed. These graphics differentiate system modes from spurious noise modes for a given order system. Additionally, using the modal assurance criterion, the experimental modes obtained have been compared with those evaluated from a finite element analysis. A quite good agreement between numerical and experimental results is observed.

Fuzzy Optimal Control for Double Inverted Pendulum

In this paper a fuzzy optimal control for stabilizing an upright position a double inverted pendulum (DIP) is developed and compared. Modeling is based on Euler-Lagrange equations. This results in a complicated nonlinear fast reaction, unstable multivariable system. Firstly, the mathematical models of double pendulum system are presented. The weight variable fuzzy input is gained by combining the fuzzy control theory with the optimal control theory. Simulation results show that the controller, which the upper pendulum is considered as main control variable, has high accuracy, quick convergence speed and higher precision.

AUTOFLY-Aid: Flight Deck Automation Support with Dynamic 4D Trajectory Management for Responsive and Adaptive Airborne Collision Avoidance

AUTOFLY-Aid Project aims to develop and demonstrate novel automation support algorithms and tools to the flight crew for flight critical collision avoidance using “dynamic 4D trajectory management”. The automation support system is envisioned to improve the primary shortcomings of TCAS, and to aid the pilot through add-on avionics/head-up displays and reality augmentation devices in dynamically evolving collision avoidance scenarios. The main theoretical innovative and novel concepts to be developed by AUTOFLY-Aid project are a) design and development of the mathematical models of the full composite airspace picture from the flight deck’s perspective, as seen/measured/informed by the aircraft flying in SESAR 2020, b) design and development of a dynamic trajectory planning algorithm that can generate at real-time (on the order of seconds) flyable (i.e. dynamically and performance-wise feasible) alternative trajectories across the evolving stochastic composite airspace picture (which includes new conflicts, blunder risks, terrain and weather limitations) and c) development and testing of the Collision Avoidance Automation Support System on a Boeing 737 NG FNPT II Flight Simulator with synthetic vision and reality augmentation while providing the flight crew with quantified and visual understanding of collision risks in terms of time and directions and countermeasures.

Modeling Piezoelectric Harvesting Materials in Road Traffic Applications.

En esta comunicación se presenta el método para obtener modelos equivalentes eléctricos de materiales piezoeléctricos utilizados en entornos con tráfico vial para aplicaciones "Energy Harvesting". Los resultados experimentales se procesan para determinar la estructura topológica óptima y la tecnología de los elementos semiconductores utilizados en la etapa de entrada del sistema de alimentación "harvesting". Asimismo se presenta el modelo de la fuente de alimentación no regulada bajo demanda variable de corriente. Abstract: The method to obtain electrical equivalent models of piezoelectric materials used in energy harvesting road traffic environment is presented in this paper. The experimental results are processed in order to determine the optimal topological structure and technology of the semiconductor elements used in the input stage of the power harvesting system. The non regulated power supply model under variable current demand is also presented.

Morteros acumuladores con parafinas microencapsuladas para el aprovechamiento de la energía solar en suelos radiantes

Esta Tesis plantea la pregunta de si el uso de morteros con parafinas microencapsuladas combinado con colectores solares térmicos puede reducir el consumo de energías convencionales, en un sistema tradicional de suelo radiante. Se pretende contribuir al conocimiento acerca del efecto que produce en el edificio, el calor latente acumulado en suelos radiantes, utilizando morteros de cemento Portland con material de cambio de fase (PCM), en conjunto con la energía solar. Para cumplir con este propósito, la investigación se desarrolla considerando diversos aspectos. En primer lugar, se revisa y analiza la documentación disponible en la actualidad, de almacenamiento de energía mediante calor latente en la construcción, y en particular la aplicación de microcápsulas de PCM en morteros y suelos radiantes. También se revisa la documentación relacionada con la aplicación de la energía solar térmica y en suelo radiante. Se analiza la normativa vigente respecto al material, a los colectores solares y al suelo radiante. Se verifica que no hay normativa relacionada con mortero-PCM, debido a esto se aplica en la investigación una adaptación de la existente. La fase experimental desarrollada esta principalmente dirigida a la cuantificación, caracterización y evaluación de las propiedades físicas, mecánicas y térmicas del mortero de cemento Portland con parafinas microencapsuladas. Los resultados obtenidos y su análisis, permiten conocer el comportamiento de este tipo de morteros, con las diferentes variables aplicadas en la investigación. Además, permite disponer de la información necesaria, para crear una metodología para el diseño de morteros con parafina microencapsulada, tanto del punto de vista de su resistencia a la compresión y contenido de PCM, como de su comportamiento térmico como acumulador de calor. Esto se logra procesando la información obtenida y generando modelos matemáticos, para dosificar mezclas, y predecir la acumulación de calor en función de su composición. Se determinan los tipos y cantidades de PCM, y el cemento más adecuado. Se obtienen importantes conclusiones respecto a los aspectos constructivos a considerar en la aplicación de morteros con PCM, en suelo radiante. Se analiza y evalúa la demanda térmica que se puede cubrir con el suelo radiante, utilizando morteros con parafina microencapsulada, a través de la acumulación de energía solar producida por colectores solares, para condiciones climáticas, técnicas y tipologías constructivas específicas. Se determina que cuando los paneles cubren más de 60 % de la demanda por calefacción, se puede almacenar en los morteros con PCM, el excedente generado durante el día. Se puede cubrir la demanda de acumulación de energía con los morteros con PCM, en la mayoría de los casos analizados. Con esto, se determina que el uso de morteros con PCM, aporta a la eficiencia energética de los edificios, disminuyendo el consumo de energías convencionales, reemplazándola por energía solar térmica. En esta investigación, el énfasis está en las propiedades del material mortero de cemento-PCM y en poder generar metodologías que faciliten su uso. Se aborda el uso de la energía solar, para verificar que es posible su acumulación en morteros con PCM aplicados en suelo radiante, posibilitando el reemplazo de energías convencionales. Quedan algunos aspectos de la aplicación de energía solar a suelo radiante con morteros con PCM, que no han sido tratados con la profundidad que requieren, y que resultan interesantes de evaluar en este tipo de aplicaciones constructivas, como entre otros, los relacionados con la cuantificación de los ahorros de energía en las diferentes estaciones del año, de la estabilización de temperaturas internas, su análisis de costo y la optimización de este tipo de sistemas para utilización en verano, los que dan pie para otras Tesis o proyectos de investigación. ABSTRACT This Thesis proposes the question of whether the use of mortars with microencapsulated paraffin combined with solar thermal collectors can reduce conventional energy consumption in a traditional heating floor system. It aims to contribute to knowledge about the effect that it has on the building, the latent heat accumulated in heating floor, using Portland cement mortars with phase change material (PCM), in conjunction with solar energy. To fulfill this purpose, the research develops it considering various aspects. First, it reviews and analyzes the documentation available today, about energy storage by latent heat in the building, and in particular the application of PCM microcapsules in mortars and heating floors. It also reviews the documentation related to the application of solar thermal energy and heating floor. Additionally, it analyzes the current regulations regarding to material, solar collectors and heating floors. It verifies that there aren’t regulations related to PCM mortar, due to this, it applies an adaptation in the investigation. The experimental phase is aimed to the quantification, mainly, characterization and evaluation of physical, mechanical and thermal properties of Portland cement mortar with microencapsulated paraffin. The results and analysis, which allow us to know the behavior of this type of mortars with different variables applied in research. It also allows having the information necessary to create a methodology for designing mortars with microencapsulated paraffin, both from the standpoint of its resistance to compression and PCM content, and its thermal performance as a heat accumulator. This accomplishes by processing the information obtained, and generating mathematical models for dosing mixtures, and predicting heat accumulation depending on their composition. The research determines the kinds and amounts of PCM, and the most suitable cement. Relevant conclusions obtain it regarding constructive aspects to consider in the implementation of PCM mortars in heating floor. Also, it analyzes and evaluates the thermal demand that it can be covered in heating floor using microencapsulated paraffin mortars, through the accumulation of solar energy produced by solar collectors to weather conditions, technical and specific building typologies. It determines that if the panels cover more than 60% of the demand for heating, the surplus generated during the day can be stored in PCM mortars. It meets the demand of energy storage with PCM mortars, in most of the cases analyzed. With this, it determines that the use of PCM mortars contributes to building energy efficiency, reducing consumption of conventional energy, replacing it with solar thermal energy. In this research approaches the use of solar energy to determine that it’s possible to verify its accumulation in PCM mortars applied in heating floor, enabling the replacement of conventional energy. The emphasis is on material properties of PCM mortar and, in order to generate methodologies to facilitate their use. There are some aspects of solar energy application in PCM mortars in heating floor, which have not been discussed with the depth required, and that they are relevant to evaluate in this kind of construction applications, including among others: the applications related to the energy savings quantification in different seasons of the year, the stabilizing internal temperatures, its cost analysis and optimization of these systems for use in summer, which can give ideas for other thesis or research projects.

Robustness and Recoverability in Transport Logistics

Abstract Transport is the foundation of any economy: it boosts economic growth, creates wealth, enhances trade, geographical accessibility and the mobility of people. Transport is also a key ingredient for a high quality of life, making places accessible and bringing people together. The future prosperity of our world will depend on the ability of all of its regions to remain fully and competitively integrated in the world economy. Efficient transport is vital in making this happen. Operations research can help in efficiently planning the design and operating transport systems. Planning and operational processes are fields that are rich in combinatorial optimization problems. These problems can be analyzed and solved through the application of mathematical models and optimization techniques, which may lead to an improvement in the performance of the transport system, as well as to a reduction in the time required for solving these problems. The latter aspect is important, because it increases the flexibility of the system: the system can adapt in a faster way to changes in the environment (i.e.: weather conditions, crew illness, failures, etc.). These disturbing changes (called disruptions) often enforce the schedule to be adapted. The direct consequences are delays and cancellations, implying many schedule adjustments and huge costs. Consequently, robust schedules and recovery plans must be developed in order to fight against disruptions. This dissertation makes contributions to two different fields: rail and air applications. Robust planning and recovery methods are presented. In the field of railway transport we develop several mathematical models which answer to RENFE’s (the major railway operator in Spain) needs: 1. We study the rolling stock assignment problem: here, we introduce some robust aspects in order to ameliorate some operations which are likely to fail. Once the rolling stock assignment is known, we propose a robust routing model which aims at identifying the train units’ sequences while minimizing the expected delays and human resources needed to perform the sequences. 2. It is widely accepted that the sequential solving approach produces solutions that are not global optima. Therefore, we develop an integrated and robust model to determine the train schedule and rolling stock assignment. We also propose an integrated model to study the rolling stock circulations. Circulations are determined by the rolling stock assignment and routing of the train units. 3. Although our aim is to develop robust plans, disruptions will be likely to occur and recovery methods will be needed. Therefore, we propose a recovery method which aims to recover the train schedule and rolling stock assignment in an integrated fashion all while considering the passenger demand. In the field of air transport we develop several mathematical models which answer to IBERIA’s (the major airline in Spain) needs: 1. We look at the airline-scheduling problem and develop an integrated approach that optimizes schedule design, fleet assignment and passenger use so as to reduce costs and create fewer incompatibilities between decisions. Robust itineraries are created to ameliorate misconnected passengers. 2. Air transport operators are continuously facing competition from other air operators and different modes of transport (e.g., High Speed Rail). Consequently, airline profitability is critically influenced by the airline’s ability to estimate passenger demands and construct profitable flight schedules. We consider multi-modal competition including airline and rail, and develop a new approach that estimates the demand associated with a given schedule; and generates airline schedules and fleet assignments using an integrated schedule design and fleet assignment optimization model that captures the impacts of schedule decisions on passenger demand.

Medida de la conductividad térmica con el método de la aguja térmica, basado en la fuente lineal de calor transitorio, para su aplicación en los cerramientos de adobes y bloques de tierra comprimida

Existe un creciente interés internacional por el ahorro energético y la sostenibilidad en la edificación con importantes repercusiones en la Arquitectura. La inercia térmica es un parámetro fundamental para poder valorar energéticamente un edificio en condiciones reales. Para ello es necesario cambiar el enfoque tradicional de transmisión de calor en régimen estacionario por otro en régimen dinámico en el que se analizan las ondas térmicas y el flujo de calor oscilante que atraviesan los cerramientos. Los parámetros que definen la inercia térmica son: el espesor, la difusividad y el ciclo térmico. A su vez la difusividad está determinada por la conductividad térmica, la densidad y el calor específico del material. De estos parámetros la conductividad es el más complejo, variable y difícil de medir, especialmente en los cerramientos de tierra debido a su heterogeneidad y complejidad higrotérmica. En general, los métodos de medida de la conductividad o transmitancias en los paramentos presentan inconvenientes a la hora de medir un edificio construido con tierra: dificultades de implementación, el elevado coste o la fiabilidad de los resultados, principalmente. El Método de la Aguja Térmica (MAT) se basa en el principio de la evolución en el tiempo del calor emitido por una fuente lineal al insertarse en el seno de un material. Se ha escogido este método porque resulta práctico, de bajo coste y de fácil aplicación a gran escala pero tiene serios problemas de fiabilidad y exactitud. En esta tesis se desarrolla un método de medida de la conductividad térmica para Piezas de Albañilería de Tierra Cruda en laboratorio basado en el MAT, se mejora su fiabilidad, se analiza su incertidumbre, se compara con otros métodos de referencia y se aplica en adobes, Bloques de Tierra Comprimida y probetas de tierra estabilizada con distintas proporciones de paja. Este método servirá de base a una posterior aplicación in situ. Finalmente se proponen modelos matemáticos para mejorar la exactitud del dispositivo utilizado y para la estimación de la conductividad de cerramientos de tierra en función de su densidad. Con los resultados obtenidos se analizan las posibilidades de amortiguación y retardo de las ondas térmicas y capacidad de almacenaje de energía de los cerramientos en función de su densidad y humedad. There is growing international interest in energy saving and sustainability in buildings with significant impact on Architecture. Thermal inertia is a key parameter to assess energy in buildings in real conditions. This requires changing the traditional approach to heat transfer in steady state by another in dynamic regime which analyzes the thermal waves and oscillating heat flux passing through the external walls. The parameters defining the thermal inertia are: the thickness, the diffusivity and the thermal cycle. In turn, the diffusivity is determined by the thermal conductivity, density and specific heat of the material. Of these parameters, thermal conductivity is the most complex, variable and difficult to measure, especially in earth walls due to their heterogeneity and hygrothermal complexity. In general, the methods of measurement of conductivity and transmittance in walls have drawbacks when measuring a building with earth: implementation difficulties, high cost, or reliability of the results, mainly. The Thermal Needle Procedure (TNP) is based on the principle of evolution in time of heat from a line source when inserted within a material. This method was chosen because it is a practical, low cost and easy to implement on a large scale but has serious problems of reliability and accuracy. This thesis develops a laboratory method for measuring the thermal conductivity of Masonry Units Unfire Earth-based based on TNP, its uncertainty is analyzed, compared to other reference methods and applies in adobes, Compressed Earth Blocks and stabilized soil specimens with different proportions of straw. This method will form the basis of a subsequent application in situ. Finally, mathematical models are proposed to improve the accuracy of the device used, and to estimate the conductivity of earth enclosures depending on its density. With the results obtained earth enclosures are analyzed to estimate their possibilities of delay and buffer of termal waves and energy storage capacity according to their density and moisture.

Estudio de la aplicabilidad de las aproximaciones uniformes al análisis de series históricas : aplicación al caso de los precios del cobre

El objetivo de este proyecto de investigación es comparar dos técnicas matemáticas de aproximación polinómica, las aproximaciones según el criterio de mínimos cuadrados y las aproximaciones uniformes (“minimax”). Se describen tanto el mercado actual del cobre, con sus fluctuaciones a lo largo del tiempo, como los distintos modelos matemáticos y programas informáticos disponibles. Como herramienta informática se ha seleccionado Matlab®, cuya biblioteca matemática es muy amplia y de uso muy extendido y cuyo lenguaje de programación es suficientemente potente para desarrollar los programas que se necesiten. Se han obtenido diferentes polinomios de aproximación sobre una muestra (serie histórica) que recoge la variación del precio del cobre en los últimos años. Se ha analizado la serie histórica completa y dos tramos significativos de ella. Los resultados obtenidos incluyen valores de interés para otros proyectos. Abstract The aim of this research project is to compare two mathematical models for estimating polynomial approximation, the approximations according to the criterion of least squares approximations uniform (“Minimax”). Describes both the copper current market, fluctuating over time as different computer programs and mathematical models available. As a modeling tool is selected main Matlab® which math library is the largest and most widely used programming language and which is powerful enough to allow you to develop programs that are needed. We have obtained different approximating polynomials, applying mathematical methods chosen, a sample (historical series) which indicates the fluctuation in copper prices in last years. We analyzed the complete historical series and two significant sections of it. The results include values that we consider relevant to other projects

Validación del método de predicción de direcciones principales de drenaje subterráneo en macizos anisótropos, en los terrenos volcánicos del macizo de Anaga, Tenerife

Actualmente y desde hace ya más de 25 años, el Método de “Predicción de las Direcciones Principales de Drenaje Subterráneo en Macizos Anisótropos”, ha sido utilizado con éxito en diferentes terrenos Kársticos como: calizas, yesos, cuarcitas, pizarras, granitos y criokarst (karst en el hielo glaciar). Sin embargo hasta ahora, nunca se había validado en terrenos volcánicos donde está focalizada esta tesis que lleva por título, Validación de dicho Método en los Terrenos Volcánicos del Macizo de Anaga en Tenerife. Este Método matemático consiste esencialmente en “Predecir y Cuantificar” las direcciones principales de drenaje turbulento subterráneo en macizos anisótropos. Para ello se basa en el estudio realizado en campo de los tectoglifos o deformaciones permanentes del macizo, impresas éstas en la roca, como consecuencia de los esfuerzos tectónicos a los que ha estado sometido dicho macizo. Se consigue de esta manera cubrir el vacío para macizos anisótropos que existe con el modelo matemático de flujo subterráneo laminar (macizos isótropos) definido por Darcy (1856). Para validar el Método se ha elegido el macizo de Anaga, pues es la zona de mayor anisotropía existente en la isla de Tenerife, conformada por una gran y extensa red de diques de diversas formas y tamaños que pertenecen a la familia de diques del eje estructural NE de la isla. En dicho macizo se realizó un exhaustivo trabajo de campo con la toma 331 datos (diques basálticos) y se aplicó el Método, consiguiendo definir las direcciones preferentes de drenaje subterráneo en el macizo de Anaga. Esta predicción obtenida se contrastó con la realidad del drenaje en la zona, conocida gracias a la existencia de cinco galerías ubicadas en la zona trabajo, de las cuales se tiene información sobre sus alumbramientos. En todos los casos se demuestra la bondad de la predicción obtenida con el Método. Queda demostrado que a mayor caos geológico o geotectónico, se ha conseguido mejor predicción del Método, obteniéndose resultados muy satisfactorios para aquellas galerías de agua en las que su rumbo de avance fue coincidente con la dirección perpendicular a la obtenida con la predicción dada por el Método, como dirección preferente de drenaje en la zona en la que se encuentra ubicada cada galería. No cabe duda que la validación de Método en los terrenos volcánicos de Tenerife, supondrá un cambio considerable en el mundo de la hidrogeología en este tipo de terrenos. Es la única herramienta matemática que se dispone para predecir un rumbo acertado en el avance de la perforación de las galerías de aguas, lo que conlleva al mismo tiempo un ahorro importantísimo en la ejecución de las obras. Por otro lado, el Método deja un importante legado a la sociedad canaria, pues con él se abren numerosas vías de trabajo e investigación que generarán un importante desarrollo en el mundo de la hidrogeología volcánica. ABSTRACT Currently and for over 25 years now, the Method of "Prediction of Subsurface Drainage Main Directions in Anisotropic Massifs" has been successfully used in various karstic terrains such as: limestone, gypsum, quartzite, slate, granite and criokarst (karst in the glacier ice). However, until now, it had never been validated in volcanic terrains where is focused this thesis entitled Validation of such Method in the Anaga Massif Volcanic Terrains, in Tenerife. This mathematical method is essentially "predict and quantify" the main directions of groundwater turbulent drainage in anisotropic massifs. This is based on field study of tectoglifes or permanent deformation of the massif, printed on the rocks as a result of previous tectonic stresses. Therefore it is possible to use in anisotropic rock mathematical model instead of the isotropic laminar flow mathematical models defined by Darcy (1856). The Anaga Massif have been chosen to validate the method, because it presents the greatest anisotropy in Tenerife Island, shaped by a large and extensive network of dikes of various shapes and sizes that belong to the family of NE structural axis dikes of the island. An exhaustive field work was carried out in such massif, with 331 collected data (basaltic dikes) and the method was applied, in order to define the preferred direction of the underground drainage in the Anaga massif. This obtained prediction was contrasted to the reality of the drainage in the area, known thanks to the existence of five galleries located in the work area, from which information about their springs was available. In all cases it was possible to demonstrate the fitness of the prediction obtained by the method. It had been demonstrated that a greater geological or geotectonic chaos enhances a better prediction of the method, that predicted very satisfactory results for those water galleries which directions were perpendicular to that predicted by the Method as a drainage preferential direction, for the zone where was located each gallery. No doubt that the validation of the use of the Method in the volcanic terrain of Tenerife, means a considerable change in the world of hydrogeology in this type of terrain. It is the only mathematical tool available to predict a successful drilling direction in advancing water galleries, what also leads to major savings in execution of the drilling works. Furthermore, the method leaves an important legacy to the Canary Islands society, because it opens many lines of work and research to generate a significant development in the world of volcanic hydrogeology.