Remarks on methods for the computation of boundary-element integrals by co-ordinate transformation

It is well known that the evaluation of the influence matrices in the boundary-element method requires the computation of singular integrals. Quadrature formulae exist which are especially tailored to the specific nature of the singularity, i.e. log(*- x0)9 Ijx- JC0), etc. Clearly the nodes and weights of these formulae vary with the location Xo of the singular point. A drawback of this approach is that a given problem usually includes different types of singularities, and therefore a general-purpose code would have to include many alternative formulae to cater for all possible cases. Recently, several authors1"3 have suggested a type independent alternative technique based on the combination of standard Gaussian rules with non-linear co-ordinate transformations. The transformation approach is particularly appealing in connection with the p.adaptive version, where the location of the collocation points varies at each step of the refinement process. The purpose of this paper is to analyse the technique in eference 3. We show that this technique is asymptotically correct as the number of Gauss points increases. However, the method possesses a 'hidden' source of error that is analysed and can easily be removed.

A bi-cubic transformation for the numerical evaluation of the Cauchy principal value integrals in boundary methods

The numerical strategies employed in the evaluation of singular integrals existing in the Cauchy principal value (CPV) sense are, undoubtedly, one of the key aspects which remarkably affect the performance and accuracy of the boundary element method (BEM). Thus, a new procedure, based upon a bi-cubic co-ordinate transformation and oriented towards the numerical evaluation of both the CPV integrals and some others which contain different types of singularity is developed. Both the ideas and some details involved in the proposed formulae are presented, obtaining rather simple and-attractive expressions for the numerical quadrature which are also easily embodied into existing BEM codes. Some illustrative examples which assess the stability and accuracy of the new formulae are included.

Absorption coefficient for the intraband transitions in quantum dot materials

In this paper, we present calculations of the absorption coefficient for transitions between the bound states of quantum dots grown within a semiconductor and the extended states of the conduction band. For completeness, transitions among bound states are also presented. In the separation of variables, single band k·p model is used in which most elements may be expressed analytically. The analytical formulae are collected in the appendix of this paper. It is concluded that the transitions are strong enough to provide a quick path to the conduction band for electrons pumped from the valence to the intermediate band

Caracterización de los sistemas de protección basados en materiales naturales destinados al control de la socavación en obras marítimas presentes en instalaciones eólicas marinas

El fenómeno de la socavación puede poner en peligro la operatividad o la seguridad de estructuras offshore comprometiendo su estabilidad global. Hasta el momento, la gran mayoría de las investigaciones destinadas a estudiar el origen y el desarrollo de este fenómeno han estado centradas en entornos fluviales, bajo condiciones de corriente continua. En la última década, debido al crecimiento del mercado de la eólica marina, han surgido numerosos estudios para la caracterización de este fenómeno en el entorno marino, teniendo en cuenta que estas estructuras se encuentran sometidas de forma simultánea a los efectos de la corriente y el oleaje, y las corrientes provocadas por las mareas. Ante la observada carencia de criterios existentes para el diseño de protecciones frente a la socavación, la presente Tesis Doctoral surge con el objetivo principal de desarrollar una serie de recomendaciones que permitan mejorar y optimizar el diseño de estas estructuras, teniendo en cuenta no sólo los parámetros geométricos de las cimentaciones de los aerogeneradores, o la propia erosión, sino también, variables características del oleaje como la altura de ola, el periodo o la longitud de onda, así como la profundidad o la batimétrica de la cimentación. La caracterización de los sistemas de protección basados en materiales naturales destinados al control de la socavación en obras marítimas presentes en instalaciones eólicas marinas no es el único objetivo alcanzado en la presente Tesis Doctoral. A través de la calibración del parámetro de altura de ola adimensional (H0) en diferentes parques eólicos europeos, y de acuerdo al criterio propuesto por Van der Meer (1988), se propone la clasificación de este tipo de estructuras de protección, basadas en material granular o escollera. La información recopilada sobre la socavación registrada en numerosos parques eólicos con protección también ha permitido estudiar la funcionalidad de este tipo de protecciones instaladas hasta la fecha. Asímismo, gracias al conocimiento adquirido sobre el fenómeno de la socavación, se plantea una propuesta para la mejora de la caracterización de este fenómeno en ambiente marino basada en el estudio de la influencia del campo de aceleraciones. La presente Tesis Doctoral permite mejorar el diseño de las protecciones frente a la socavación que se utilizan en parques eólicos marinos teniendo en cuenta las acciones del clima marítimo, supliendo de este modo las carencias que hasta el momento presentan las formulaciones existentes, las cuales únicamente toman en consideración la geometría de las cimentaciones, el ángulo de rozamiento interno del terreno y la estimación de la máxima socavación que puede llegar a producirse. Scour phenomenon jeopardizes the stability and functionality of offshore structures compromising its overall stability. So far, most studies about the origin and the development of this phenomenon have been focused on river environments (under steady current conditions). In the last decade a lot of research projects about the characterization of this phenomenon have been carried out due to the growth of offshore wind industry. These projects take into account that these structures are subjected simultaneously to current, waves and tidal effects. This PhD Thesis arises due to the current lack of criteria for the design of scour protections. Its main objective is to develop some recommendations to improve and to optimize the design of scour protection structures. For that it is necessary to take into account not only the geometrical parameters of foundations or the erosion forecasted, but also wave variables such as wave height, wave period or wavelength. Characterization of protection systems based on natural materials for the control of the scour in offshore wind farms was not the only goal achieved in this PhD Thesis. Through the calibration of the dimensionless wave height parameter (H0) in different European offshore wind farms, and according to the criteria proposed by Van der Meer (1988), a classification of these protection structures based on natural elements (rocks or riprap) has been proposed. Scour data registered in numerous offshore wind farms with scour protection systems also allowed to study the functionality of this type of protection installed up to now. Thanks to the knowledge acquired about the scour development, a proposal for the improvement of the characterization of this phenomenon in marine environment is proposed. This has been based on the study of the influence of the acceleration parameters. This PhD Thesis improves the design of scour protections used in offshore wind facilities taking into account maritime climate actions. To solve the current formulae deficiencies only considering the foundation geometry, the internal friction angle of the seabed and the maximum scour depth forecasted.

Order 10 4 speedup in global linear instability analysis using matrix formation

A unified solution framework is presented for one-, two- or three-dimensional complex non-symmetric eigenvalue problems, respectively governing linear modal instability of incompressible fluid flows in rectangular domains having two, one or no homogeneous spatial directions. The solution algorithm is based on subspace iteration in which the spatial discretization matrix is formed, stored and inverted serially. Results delivered by spectral collocation based on the Chebyshev-Gauss-Lobatto (CGL) points and a suite of high-order finite-difference methods comprising the previously employed for this type of work Dispersion-Relation-Preserving (DRP) and Padé finite-difference schemes, as well as the Summationby- parts (SBP) and the new high-order finite-difference scheme of order q (FD-q) have been compared from the point of view of accuracy and efficiency in standard validation cases of temporal local and BiGlobal linear instability. The FD-q method has been found to significantly outperform all other finite difference schemes in solving classic linear local, BiGlobal, and TriGlobal eigenvalue problems, as regards both memory and CPU time requirements. Results shown in the present study disprove the paradigm that spectral methods are superior to finite difference methods in terms of computational cost, at equal accuracy, FD-q spatial discretization delivering a speedup of ð (10 4). Consequently, accurate solutions of the three-dimensional (TriGlobal) eigenvalue problems may be solved on typical desktop computers with modest computational effort.

Influence of the Second Flexural Mode on the Response of High-Speed Bridges

This paper deals with the assessment of the contribution of the second flexural mode to the dynamic behaviour of simply supported railway bridges. Alluding to the works of other authors, it is suggested in some references that the dynamic behaviour of simply supported bridges could be adequately represented taking into account only the contribution of the fundamental flexural mode. On the other hand, the European Rail Research Institute (ERRI) proposes that the second mode should also be included whenever the associated natural frequency is lower than 30 Hz]. This investigation endeavours to clarify the question as much as possible by establishing whether the maximum response of the bridge, in terms of displacements, accelerations and bending moments, can be computed accurately not taking account of the contribution of the second mode. To this end, a dimensionless formulation of the equations of motion of a simply supported beam traversed by a series of equally spaced moving loads is presented. This formulation brings to light the fundamental parameters governing the behaviour of the beam: damping ratio, dimensionless speed $ \alpha$=VT/L, and L/d ratio (L stands for the span of the beam, V for the speed of the train, T represents the fundamental period of the bridge and d symbolises the distance between consecutive loads). Assuming a damping ratio equal to 1%, which is a usual value for prestressed high-speed bridges, a parametric analysis is conducted over realistic ranges of values of $ \alpha$ and L/d. The results can be extended to any simply supported bridge subjected to a train of equally spaced loads in virtue of the so-called Similarity Formulae. The validity of these formulae can be derived from the dimensionless formulation mentioned above. In the parametric analysis the maximum response of the bridge is obtained for one thousand values of speed that cover the range from the fourth resonance of the first mode to the first resonance of the second mode. The response at twenty-one different locations along the span of the beam is compared in order to decide if the maximum can be accurately computed with the sole contribution of the fundamental mode.

Modelling of Topographic Irregularities for Seismic Site Response

Seismic evaluation methodology is applied to an existing viaduct in the south of Spain, near Granada, which is a medium seismicity region. The influence of both geology and topography in the spatial variability of ground motion are studied as well as seismic hazard analysis and ground motion characterization. Artificial hazard-consistent ground motion records are synthesised applying seismic hazard analysis and site effects are estimated through a diffraction study. Direct BEM is used to calculate the valley displacement response to vertically propagating SV waves and transfer functions are generated allowing the transformation of free field motion to motion at each support. A closed formulae is used to estimate these transfer function. Finally, the results obtained are compared.

Geometric characteristics of conics in Bézier form

Synthetic derivation of closed formulae of the geometric characteristic of a conic given in Bézier form in terms of its control polygon, (P; Q; R) and weights, (1; w; 1g)

Efficient algorithms for dynamic probabilistic safety assessment. An application to convex damage domain.

The design of nuclear power plant has to follow a number of regulations aimed at limiting the risks inherent in this type of installation. The goal is to prevent and to limit the consequences of any possible incident that might threaten the public or the environment. To verify that the safety requirements are met a safety assessment process is followed. Safety analysis is as key component of a safety assessment, which incorporates both probabilistic and deterministic approaches. The deterministic approach attempts to ensure that the various situations, and in particular accidents, that are considered to be plausible, have been taken into account, and that the monitoring systems and engineered safety and safeguard systems will be capable of ensuring the safety goals. On the other hand, probabilistic safety analysis tries to demonstrate that the safety requirements are met for potential accidents both within and beyond the design basis, thus identifying vulnerabilities not necessarily accessible through deterministic safety analysis alone. Probabilistic safety assessment (PSA) methodology is widely used in the nuclear industry and is especially effective in comprehensive assessment of the measures needed to prevent accidents with small probability but severe consequences. Still, the trend towards a risk informed regulation (RIR) demanded a more extended use of risk assessment techniques with a significant need to further extend PSA’s scope and quality. Here is where the theory of stimulated dynamics (TSD) intervenes, as it is the mathematical foundation of the integrated safety assessment (ISA) methodology developed by the CSN(Consejo de Seguridad Nuclear) branch of Modelling and Simulation (MOSI). Such methodology attempts to extend classical PSA including accident dynamic analysis, an assessment of the damage associated to the transients and a computation of the damage frequency. The application of this ISA methodology requires a computational framework called SCAIS (Simulation Code System for Integrated Safety Assessment). SCAIS provides accident dynamic analysis support through simulation of nuclear accident sequences and operating procedures. Furthermore, it includes probabilistic quantification of fault trees and sequences; and integration and statistic treatment of risk metrics. SCAIS comprehensively implies an intensive use of code coupling techniques to join typical thermal hydraulic analysis, severe accident and probability calculation codes. The integration of accident simulation in the risk assessment process and thus requiring the use of complex nuclear plant models is what makes it so powerful, yet at the cost of an enormous increase in complexity. As the complexity of the process is primarily focused on such accident simulation codes, the question of whether it is possible to reduce the number of required simulation arises, which will be the focus of the present work. This document presents the work done on the investigation of more efficient techniques applied to the process of risk assessment inside the mentioned ISA methodology. Therefore such techniques will have the primary goal of decreasing the number of simulation needed for an adequate estimation of the damage probability. As the methodology and tools are relatively recent, there is not much work done inside this line of investigation, making it a quite difficult but necessary task, and because of time limitations the scope of the work had to be reduced. Therefore, some assumptions were made to work in simplified scenarios best suited for an initial approximation to the problem. The following section tries to explain in detail the process followed to design and test the developed techniques. Then, the next section introduces the general concepts and formulae of the TSD theory which are at the core of the risk assessment process. Afterwards a description of the simulation framework requirements and design is given. Followed by an introduction to the developed techniques, giving full detail of its mathematical background and its procedures. Later, the test case used is described and result from the application of the techniques is shown. Finally the conclusions are presented and future lines of work are exposed.

Estudio del comportamiento del hormigón de alta resistencia reforzado con fibras de acero frente al impacto de proyectiles

En la presente investigación se buscó estudiar el efecto de la adición de fibras metálicas como refuerzo en hormigones de alta resistencia, y en especial su comportamiento frente al impacto de proyectiles. Se efectuó el estudio sobre un hormigón de alta resistencia (HAR), analizando los aspectos mecánicos, durabilidad y trabajabilidad para su colocación en obra. Las pruebas de laboratorio se llevaron a cabo en el Laboratorio de Materiales de Construcción de la Escuela Técnica Superior de Caminos Canales y Puertos de la UPM y los ensayos balísticos en la galería de tiro cubierta del Polígono de Experiencia de Carabanchel, adscrito a la Dirección General de Infraestructura del Ministerio de la Defensa. La caracterización del HAR empleado en el estudio se centró en los aspectos de resistencias mecánicas a compresión, tracción, flexotracción, tenacidad a flexotracción, punzonamiento, retracción, fluencia, temperatura interna y resistencia al impacto de proyectiles, siempre buscando de manera primordial analizar el efecto de la adición de fibras en el hormigón de alta resistencia. El programa de ensayos balísticos comprendió la fabricación de 47 placas de hormigón de diferentes espesores, desde 5 a 40 cm., 26 de dichas placas eran de HAR con una adición de fibras metálicas de 80 kg/m3, 11 de ellas eran de HAR sin fibras y 10 de un hormigón de resistencia convencional con y sin fibras; sobre dichas placas se efectuaron diversos impactos con proyectiles de los cuatro calibres siguientes: 7.62 AP, 12.70 M8, 20 mm APDS y 25 mm APDS. Las pruebas mostraron que el HAR presenta una mayor resistencia a los impactos de proyectiles, aunque sin la adición de fibras su fragilidad es un serio inconveniente para su utilización como barrera protectora, la adición de fibras reduce considerablemente la fragmentación en la cara posterior “scabbing” y en menor medida en la cara anterior “spalling”. También se incrementa la capacidad del hormigón a la resistencia de múltiples impactos. Se efectuó un estudio de las diferentes formulas y modelos, en especial el modelo desarrollado por Moreno [60], que se vienen utilizando para el diseño de barreras protectoras de hormigón contra impacto de proyectiles, analizando su viabilidad en el caso del hormigón de alta resistencia, hormigón para el cual no fueron desarrolladas y para el que no existen bases de cálculo específicas. In this research we have tried to study the effect of adding metallic fibres as a means of reinforcing high strength concrete, and especially its behaviour when impacted upon by projectiles. The study was carried out using high strength concrete (HSC), analysing its mechanical facets, durability and malleability when used in construction. The laboratory tests took place in the Laboratorio de Materiales de Construcción of the Escuela Técnica Superior de Caminos Canales y Puertos of the Universidad Politécnica de Madrid, and the ballistic tests were carried out in the covered shooting gallery of the Polígono de Experiencias in Carabanchel (Madrid), belongs to the Departamento de Infraestructura of the Ministerio de Defensa. The aspects of the HSC studied are its mechanical strength to compression, traction, flexotraction, resilience to flexo-traction, shear strength, creep, shrinkage, internal temperature and strength to the impact of projectiles, always looking to analyse the effect of adding fibres to HSC. The ballistic testing process required the construction of 47 concrete plates of different thicknesses, from 5 to 40 cm, 26 made which HSC containing of 80 kg/m3 metallic fibres of, 11 made of HSC without fibres, and 10 made with concrete of normal strength with and without fibres. These plates were subjected to a variety of impacts by four projectile, 7.62 AP, 12.70 M8, 20 mm APDS and 25 mm APDS. The results showed that HSC has a greater resistance to the impact of projectiles, although without the addition of fibres, its fragility makes it much less suitable for use as a protective barrier. The addition of fibres reduces considerably frontal fragmentation, known as “scabbing”, and to a lesser extent causes fragmentation of the reverse side, known as “spalling”. In addition, the concrete’s capacity to resist multiple impacts is improved by its letter ductility. A study was carried out on the various formulae and models used to design protective concrete barriers impacted on by projectiles, analysing their viability in the case of HSC for which they were not developed and for which no specific calculations exist.

A theory of electron collection by a grid sphere

Use of a spherical grid as electron collector at the anodic end of a tether, as recently proposed, is considered. The standard analysis of space-charge limited current to a solid sphere (with neither magnetic nor plasma-motion effects), which has been shown to best fit TSS1R in-orbit results at very high bias, is used to determine effects from grid transparency on current collected; the analysis is first reformulated in the formalism recently introduced in the two-dimensional analysis of bare-tethers. A discussion of the electric potential created by a spherical grid in vacuum is then carried out; it is shown that each grid-wire collects current well below its maximum OML current, the effective grid transparency being close to its optical value. Formulae for the current to a spherical grid, showing the effects of grid transparency, is determined. A fully consistent analysis of electric potential and electron density, outside and inside the grid, is completed.

SMS freeforms for illumination

SMS 3D (simultaneous multiple surfaces in their three-dimensional version) is a well-known design method comprising two freeform surfaces that allow the perfect coupling of two wavefronts with another two. The design algorithm provides a collection of line pairs on both surfaces (called SMS spines), whose three-dimensional shape seems arbitrary at first sight. This paper shows that the shapes of the spines are partially governed by applying the étendue conservation theorem to the biparametric bundle of rays linking the paired spines, which is one lesser known étendue invariants found by Poincaré. The resulting formulae for the spines in three-dimensional space happen to coincide with the conventional étendue formulas of two-dimensional geometry, like for instance, the Hottel formula.

Un método de diferencias finitas generalizadas para simular la actividad eléctrica de un tejido cardiaco

Una evolución del método de diferencias finitas ha sido el desarrollo del método de diferencias finitas generalizadas (MDFG) que se puede aplicar a mallas irregulares o nubes de puntos. En este método se emplea una expansión en serie de Taylor junto con una aproximación por mínimos cuadrados móviles (MCM). De ese modo, las fórmulas explícitas de diferencias para nubes irregulares de puntos se pueden obtener fácilmente usando el método de Cholesky. El MDFG-MCM es un método sin malla que emplea únicamente puntos. Una contribución de esta Tesis es la aplicación del MDFG-MCM al caso de la modelización de problemas anisótropos elípticos de conductividad eléctrica incluyendo el caso de tejidos reales cuando la dirección de las fibras no es fija, sino que varía a lo largo del tejido. En esta Tesis también se muestra la extensión del método de diferencias finitas generalizadas a la solución explícita de ecuaciones parabólicas anisótropas. El método explícito incluye la formulación de un límite de estabilidad para el caso de nubes irregulares de nodos que es fácilmente calculable. Además se presenta una nueva solución analítica para una ecuación parabólica anisótropa y el MDFG-MCM explícito se aplica al caso de problemas parabólicos anisótropos de conductividad eléctrica. La evidente dificultad de realizar mediciones directas en electrocardiología ha motivado un gran interés en la simulación numérica de modelos cardiacos. La contribución más importante de esta Tesis es la aplicación de un esquema explícito con el MDFG-MCM al caso de la modelización monodominio de problemas de conductividad eléctrica. En esta Tesis presentamos un algoritmo altamente eficiente, exacto y condicionalmente estable para resolver el modelo monodominio, que describe la actividad eléctrica del corazón. El modelo consiste en una ecuación en derivadas parciales parabólica anisótropa (EDP) que está acoplada con un sistema de ecuaciones diferenciales ordinarias (EDOs) que describen las reacciones electroquímicas en las células cardiacas. El sistema resultante es difícil de resolver numéricamente debido a su complejidad. Proponemos un método basado en una separación de operadores y un método sin malla para resolver la EDP junto a un método de Runge-Kutta para resolver el sistema de EDOs de la membrana y las corrientes iónicas. ABSTRACT An evolution of the method of finite differences has been the development of generalized finite difference (GFD) method that can be applied to irregular grids or clouds of points. In this method a Taylor series expansion is used together with a moving least squares (MLS) approximation. Then, the explicit difference formulae for irregular clouds of points can be easily obtained using a simple Cholesky method. The MLS-GFD is a mesh-free method using only points. A contribution of this Thesis is the application of the MLS-GFDM to the case of modelling elliptic anisotropic electrical conductivity problems including the case of real tissues when the fiber direction is not fixed, but varies throughout the tissue. In this Thesis the extension of the generalized finite difference method to the explicit solution of parabolic anisotropic equations is also given. The explicit method includes a stability limit formulated for the case of irregular clouds of nodes that can be easily calculated. Also a new analytical solution for homogeneous parabolic anisotropic equation has been presented and an explicit MLS- GFDM has been applied to the case of parabolic anisotropic electrical conductivity problems. The obvious difficulty of performing direct measurements in electrocardiology has motivated wide interest in the numerical simulation of cardiac models. The main contribution of this Thesis is the application of an explicit scheme based in the MLS-GFDM to the case of modelling monodomain electrical conductivity problems using operator splitting including the case of anisotropic real tissues. In this Thesis we present a highly efficient, accurate and conditionally stable algorithm to solve a monodomain model, which describes the electrical activity in the heart. The model consists of a parabolic anisotropic partial differential equation (PDE), which is coupled to systems of ordinary differential equations (ODEs) describing electrochemical reactions in the cardiac cells. The resulting system is challenging to solve numerically, because of its complexity. We propose a method based on operator splitting and a meshless method for solving the PDE together with a Runge-Kutta method for solving the system of ODE’s for the membrane and ionic currents.

Ice volume estimates from ground-penetrating radar surveys, western Nordenskiöld Land glaciers, Svalbard

As part of ongoing work within the SvalGlac project aimed to obtain a reliable estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea level rise, in this contribution we present volume calculations, with detailed error estimates, for ten glaciers on western Nordenskiöld Land, central Spitsbergen, Svalbard. The volume estimates are based upon a dense net of GPR-retrieved ice thickness data collected over several field campaigns spanning the period 1999-2012, all of them except one within 2010-2012. The total area and volume of the ensemble are 113.38±0.09 km2 and 10.439±0.185 km3, respectively, while the individual areas, volumes and average ice thickness lie within 2.5-49.1 km2, 0.08-5.48 km3 and 29-108 m, respectively. The maximum recorded ice thickness, 265±15 m, corresponds to Fridtjovbreen, which has also the largest average thickness (108±1m). Available empirical formulae for Svalbard glaciers overestimate the total volume of these glaciers by 24% with respect to our calculation. On the basis of the pattern of scattering in the radargrams, we also analyse the hydrothermal structure of these glaciers. Nine out of ten are polythermal, while only one is entirely cold.

Analysis of the geometric altimetry to support aircraft optimal trajectories within the future 4d trajectory management

Over the past few years, the common practice within air traffic management has been that commercial aircraft fly by following a set of predefined routes to reach their destination. Currently, aircraft operators are requesting more flexibility to fly according to their prefer- ences, in order to achieve their business objectives. Due to this reason, much research effort is being invested in developing different techniques which evaluate aircraft optimal trajectory and traffic synchronisation. Also, the inefficient use of the airspace using barometric altitude overall in the landing and takeoff phases or in Continuous Descent Approach (CDA) trajectories where currently it is necessary introduce the necessary reference setting (QNH or QFE). To solve this problem and to permit a better airspace management born the interest of this research. Where the main goals will be to evaluate the impact, weakness and strength of the use of geometrical altitude instead of the use of barometric altitude. Moreover, this dissertation propose the design a simplified trajectory simulator which is able to predict aircraft trajectories. The model is based on a three degrees of freedom aircraft point mass model that can adapt aircraft performance data from Base of Aircraft Data, and meteorological information. A feature of this trajectory simulator is to support the improvement of the strategic and pre-tactical trajectory planning in the future Air Traffic Management. To this end, the error of the tool (aircraft Trajectory Simulator) is measured by comparing its performance variables with actual flown trajectories obtained from Flight Data Recorder information. The trajectory simulator is validated by analysing the performance of different type of aircraft and considering different routes. A fuel consumption estimation error was identified and a correction is proposed for each type of aircraft model. In the future Air Traffic Management (ATM) system, the trajectory becomes the fundamental element of a new set of operating procedures collectively referred to as Trajectory-Based Operations (TBO). Thus, governmental institutions, academia, and industry have shown a renewed interest for the application of trajectory optimisation techniques in com- mercial aviation. The trajectory optimisation problem can be solved using optimal control methods. In this research we present and discuss the existing methods for solving optimal control problems focusing on direct collocation, which has received recent attention by the scientific community. In particular, two families of collocation methods are analysed, i.e., Hermite-Legendre-Gauss-Lobatto collocation and the pseudospectral collocation. They are first compared based on a benchmark case study: the minimum fuel trajectory problem with fixed arrival time. For the sake of scalability to more realistic problems, the different meth- ods are also tested based on a real Airbus 319 El Cairo-Madrid flight. Results show that pseudospectral collocation, which has shown to be numerically more accurate and computa- tionally much faster, is suitable for the type of problems arising in trajectory optimisation with application to ATM. Fast and accurate optimal trajectory can contribute properly to achieve the new challenges of the future ATM. As atmosphere uncertainties are one of the most important issues in the trajectory plan- ning, the final objective of this dissertation is to have a magnitude order of how different is the fuel consumption under different atmosphere condition. Is important to note that in the strategic phase planning the optimal trajectories are determined by meteorological predictions which differ from the moment of the flight. The optimal trajectories have shown savings of at least 500 [kg] in the majority of the atmosphere condition (different pressure, and temperature at Mean Sea Level, and different lapse rate temperature) with respect to the conventional procedure simulated at the same atmosphere condition.This results show that the implementation of optimal profiles are beneficial under the current Air traffic Management (ATM).