Geochemistry of rare earth elements in basalts from Walvis Ridge

Selected basalts from a suite of dredged and drilled samples (IPOD sites 525, 527, 528 and 530) from the Walvis Ridge have been analysed to determine their rare earth element (REE) contents in order to investigate the origin and evolution of this major structural feature in the South Atlantic Ocean. All of the samples show a high degree of light rare earth element (LREE) enrichment, quite unlike the flat or depleted patterns normally observed for normal mid-ocean ridge basalts (MORBs). Basalts from Sites 527, 528 and 530 show REE patterns characterised by an arcuate shape and relatively low (Ce/Yb)N ratios (1.46-5.22), and the ratios show a positive linear relationship to Nb content. A different trend is exhibited by the dredged basalts and the basalts from Site 525, and their REE patterns have a fairly constant slope, and higher (Ce/Yb)N ratios (4.31-8.50). These differences are further reflected in the ratios of incompatible trace elements, which also indicate considerable variations within the groups. Mixing hyperbolae for these ratios suggest that simple magma mixing between a 'hot spot' type of magma, similar to present-day volcanics of Tristan da Cunha, and a depleted source, possibly similar to that for magmas being erupted at the Mid-Atlantic Ridge, was an important process in the origin of parts of the Walvis Ridge, as exemplified by Sites 527, 528 and 530. Site 525 and dredged basalts cannot be explained by this mixing process, and their incompatible element ratios suggest either a mantle source of a different composition or some complexity to the mixing process. In addition, the occurrence of different types of basalt at the same location suggests there is vertical zonation within the volcanic pile, with the later erupted basalts becoming more alkaline arid more enriched in incompatible elements. The model proposed for the origin and evolution of the Walvis Ridge involves an initial stage of eruption in which the magma was essentially a mixture of enriched and depleted end-member sources, with the N-MORB component being small. The dredged basalts and Site 525, which represent either later-stage eruptives or those close to the hot spot plume, probably result from mixing of the enriched mantle source with variable amounts and variable low degrees of partial melting of the depleted mantle source. As the volcano leaves the hot spot, these late-stage eruptives continue for some time. The change from tholeiitic to alkalic volcanism is probably related either to evolution in the plumbing system and magma chamber of the individual volcano, or to changes in the depth of origin of the enriched mantle source melt, similar to processes in Hawaiian volcanoes.

Hydrate abundances at the high-flux Batumi seep area in the southeastern Black Sea

Submarine gas hydrates are a major global reservoir of the potent greenhouse gas methane. Since current assessments of worldwide hydrate-bound carbon vary by one order of magnitude, new technical efforts are required for improved and accurate hydrate quantifications. Here we present hydrate abundances determined for surface sediments at the high-flux Batumi seep area in the southeastern Black Sea at 840 m water depth using state-of-the art autoclave technology. Pressure sediment cores of up to 2.65 m in length were recovered with an autoclave piston corer backed by conventional gravity cores. Quantitative core degassing yielded volumetric gas/bulk sediment ratios of up to 20.3 proving hydrate presence. The cores represented late glacial to Holocene hemipelagic sediments with the shallowest hydrates found at 90 cmbsf. Calculated methane concentrations in the different cores surpassed methane equilibrium concentrations in the two lowermost lithological Black Sea units sampled. The results indicated hydrate fractions of 5.2% of pore volume in the sapropelic Unit 2 and mean values of 21% pore volume in the lacustrine Unit 3. We calculate that the studied area of ~ 0.5 km**2 currently contains about 11.3 kt of methane bound in shallow hydrates. Episodic detachment and rafting of such hydrates is suggested by a rugged seafloor topography along with variable thicknesses in lithologies. We propose that sealing by hydrate precipitation in coarse-grained deposits and gas accumulation beneath induces detachment of hydrate/sediment chunks. Floating hydrates will rapidly transport methane into shallower waters and potentially to the sea-atmosphere boundary. In contrast, persistent in situ dissociation of shallow hydrates appears unlikely in the near future as deep water warming by about 1.6 °C and/or decrease in hydrostatic pressure corresponding to a sea level drop of about 130 m would be required. Because hydrate detachment should be primarily controlled by internal factors in this area and in similar hydrated settings, it serves as source of methane in shallow waters and the atmosphere which is mainly decoupled from external forcing.

Petrology of ODP Leg 210 and Leg 173, and DSDP Leg 47 sites

Sand detrital modes of Albian-Eocene clastic gravity-flow deposits cored and recovered at Ocean Drilling Program Site 1276 reflect the postrift geologic evolution of the Newfoundland passive continental margin. Cretaceous sandstone compositions (average: Q57F23L20; Ls%Lsc = 35; total%bioclasts = 3) are consistent with a source on Grand Banks such as Avalon Uplift. Their relatively low potassium feldspar (Qm71K8P21) contents distinguish them from Iberian sandstones and appear to preclude an easterly source during the early history of the ocean basin. Isolated volcaniclastic input near the Paleocene/Eocene boundary (~60 Ma) at Site 1276 is also present in Iberian samples of this age, suggesting that magmatism was widespread across the North Atlantic during this time frame; the source(s) of this volcanic debris remains equivocal. In the Eocene, the development of carbonate bank facies on the shelf marks a profound compositional change to calcareous grainstones (average: Q27F11L62; Ls%Lsc = 82; total%bioclasts = 55) in basinal gravity-flow deposits at Site 1276. This calcareous petrofacies is present on the Iberian margin and in the Pyrenees, suggesting that it was a regional event. The production and downslope redistribution of carbonate debris, including bioclastic and lithic fragments, was likely eustatically controlled. The Newfoundland (Site 1276 and Jeanne d'Arc Basin) sandstones are mainly quartzolithic. Their composition and the contrast in composition between them and more quartzofeldspathic sandstones from the Iberian margin are likely a product of rifting along a Paleozoic suture zone separating distinct basement terranes. This prerift geologic setting contrasts with that of rifts developed within other cratonic settings with variable amounts of synrift volcanism. When synthesized, the spectrum of synrift and postrift sand compositions produces a general model of passive margin (rift-to-drift) sandstone provenance.

Lipid analyses of ODP Leg 117 samples

The solvent-extractable organic fractions of sediment samples from six Ocean Drilling Program Leg 117 sites were investigated by gas chromatography and gas chromatography-mass spectrometry. Sediments deposited in the Indus Fan (Site 720) as well as Miocene sediments from the Owen Ridge (Sites 722 and 731) contain almost exclusively organic matter of terrigenous origin. The organic matter in sediments from the Oman Margin (Sites 723, 725, and 728) and in the Pliocene/Pleistocene sections from the Owen Ridge is mainly of a marine origin with variable admixtures of terrigenous material. In these latter samples strong variations of the lipid composition and distribution are noted. However, the interpretation of the relation to potential biological sources is hampered by a lack of information on the possible lipid composition of appropriate organisms.

Chemistry of primary and secondary phases in intraplate basalts and volcaniclastic sediments at DSDP Leg 89 Holes

Electron microprobe data are presented for clinopyroxenes, plagioclases, palagonites, smectites, celadonites, and zeolites in Hole 462A sheet-flow basalts and Site 585 volcaniclastic sediments. Glomerocrystic clinopyroxenes in Hole 462A are predominantly Ti-poor augites with minor fractionation to ferroaugites in rim portions. Quenched plumose clinopyroxenes show considerable variation from Ca-rich to Ca-poor augites, although all are characterized by being Tirich and Cr-poor relative to the glomerocrysts. Two differentiated series of Site 585 pyroxene compositions, calcic augite and diopside-salite, demonstrate the coexistence, in the vitric and lithic clasts, of tholeiitic and alkali basalt types, respectively. Plagioclase compositions in all samples are mainly labradorites, although some zoned Hole 462A glomerocrysts range from An73 to An20 and are characterized by high Mg and Fe contents in the more calcic varieties. The K content of the plagioclases is highest in the more sodic crystals, although the overall higher orthoclase component of Site 585 plagioclases reflects the generally higher bulk-rock K content. The compositions of both secondary smectites and celadonites are similar irrespective of the alteration location (glass, matrix, vesicles, etc.), although brown smectites replacing interstitial glass have marginally higher total Fe contents than pale green and yellow smectites. Analyzed zeolites are mainly phillipsites with variable alkali content, and, together with associated celadonite, represent late-stage alteration repositories for K under mildly oxidizing conditions. The compositions of both early and late secondary minerals are typical of those formed by the submarine alteration of basaltic rocks at low temperatures.

Stable oxygen isotope, UK'37, productivity, CaCO3 and XRF records from sediment core MD06-3075

Proxy records of hydrologic variability in the West Pacific Warm Pool (WPWP) have revealed wide-scale changes in past convective activity in response to orbital and sub-orbital climate forcings. However, attributing proxy responses to regional changes in WPWP hydrology versus local variations in precipitation requires independent records linking the terrestrial and marine realms. We present high-resolution stable isotope, UK'37 sea-surface temperature, X-ray fluorescence (XRF) core scanning and coccolithophore-derived paleoproductivity records covering the past 120 ka from International Marine Global Change (IMAGES) Program Core MD06-3075 (6°29' N, 125°50' E, water depth 1878 m), situated in the Davao Gulf on the southern side of Mindanao. XRF-derived log(Fe/Ca) records provide a robust proxy for runoff-driven sedimentary discharge from Mindanao, whilst past changes in local productivity are associated with variable freshwater runoff and stratification of the surface layer. Significant precessional-scale variability in sedimentary discharge occurred during Marine Isotope Stage (MIS) 5, with peaks in discharge contemporaneous with Northern Hemisphere summer insolation minima. We attribute these changes to the latitudinal migration of the Intertropical Convergence Zone (ITCZ) over the WPWP together with variability in the strength of the Walker circulation acting on precessional timescales. Between 60 and 15 ka sedimentary discharge at Mindanao was muted, displaying little orbital- or millennial-scale variability, likely in response to weakened precessional insolation forcing and lower sea level driving increased subsidence of air masses over the exposed Sunda Shelf. These results highlight the high degree of local variability in the precipitation response to past climate changes in the WPWP.

Magnetostratigraphic synthesis of ODP Leg 123 sites

During ODP Leg 123, Sites 765 and 766 were drilled to examine the tectonic evolution, sedimentary history, and paleoceanography of the Argo Abyssal Plain and lower Exmouth Plateau. At each site, the quality of magnetostratigraphic and biostratigraphic records varies because of complicating factors, such as the predominance of turbidites, the presence of condensed horizons, or deposition beneath the CCD. Based primarily on the presence of nannofossils, the base of the sedimentary section at Site 765 was dated as Tithonian. A complete Cretaceous sequence was recovered at this site, although the sedimentation rate varies markedly through the section. The Cretaceous/Tertiary boundary is represented by a condensed horizon. The condensed Cenozoic sequence at Site 765 extends from the upper Paleocene to the lower Miocene. A dramatic increase in sedimentation rate was observed in the lower Miocene, and a 480-m-thick Neogene section is present. The Neogene section is continuous, except for a minor hiatus in the lower Pliocene. The base of the sedimentary section at Site 766 is Valanginian, in agreement with the site's position on marine magnetic anomaly Mil. Valanginian to Barremian sediments are terrigenous, with variable preservation of microfossils, and younger sediments are pelagic, with abundant well-preserved microfossils. Sedimentation rate is highest in the Lower Cretaceous and decreases continually upsection. Upper Cenozoic sediments are condensed, with several hiatuses.

The estimation of truncation error by tau-estimation revisited

The aim of this paper was to accurately estimate the local truncation error of partial differential equations, that are numerically solved using a finite difference or finite volume approach on structured and unstructured meshes. In this work, we approximated the local truncation error using the @t-estimation procedure, which aims to compare the residuals on a sequence of grids with different spacing. First, we focused the analysis on one-dimensional scalar linear and non-linear test cases to examine the accuracy of the estimation of the truncation error for both finite difference and finite volume approaches on different grid topologies. Then, we extended the analysis to two-dimensional problems: first on linear and non-linear scalar equations and finally on the Euler equations. We demonstrated that this approach yields a highly accurate estimation of the truncation error if some conditions are fulfilled. These conditions are related to the accuracy of the restriction operators, the choice of the boundary conditions, the distortion of the grids and the magnitude of the iteration error.

Simulation techniques for design of overhead conductor rail lines for speeds over 140 km/h

In overhead conductor rail lines, aluminium beams are usually mounted with support spacing between 8 and 12 meters, to limit the maximum vertical deflection in the center of the span. This small support spacing limits the use of overhead conductor rail to tunnels, therefore it has been used almost exclusively in metropolitan networks, with operation speeds below 110 km/h. Nevertheless, due to the lower cost of maintenance required for this electrification system, some railway administrations are beginning to install it in some tunnels on long-distance lines, requesting higher operation speeds [1]. Some examples are the Barcelona and Madrid suburban networks (Spain), and recent lines in Turkey and Malaysia. In order to adapt the design of the overhead conductor for higher speeds (V > 160 km/h), particular attention must be paid to the geometry of the conductor rail in critical zones as overlaps, crossings and, especially, transitions between conductor rail and conventional catenary, since the use of overhead conductor rail is limited to tunnels, as already mentioned. This paper describes simulation techniques developed in order to take into account these critical zones. Furthermore, some specific simulations results are presented that have been used to analyze and optimizes the geometry of this special zones to get a better current collection quality, in a real suburban network. This paper presents the work undertaken by the Railways Technology Research Centre (CITEF), having over 10 years of experience in railways research [1-4].

Optimal Design of Envelope Amplifier Based on Linear Assisted Buck Converter

Envelope Tracking (ET) and Envelope Elimination and Restoration (EER) are two techniques that have been used as a solution for highly efficient linear RF Power Amplifiers (PA). In both techniques the most important part is a dc-dc converter called envelope amplifier that has to supply the RF PA with variable voltage. Besides high efficiency, its bandwidth is very important as well. Envelope amplifier based on parallel combination of a switching dc-dc converter and a linear regulator is an architecture that is widely used due to its simplicity. In this paper we discuss about theoretical limitations of this architecture regarding its efficiency and we demonstrate two possible way of its implementation. In order to derive the presented conclusions, a theoretical model of envelope amplifier's efficiency has been presented. Additionally, the benefits of the new emerging GaN technology for this application have been shown as well.

Testing equipment for the improvement of mechanical devices to minimize damage to fruit in commercial packing lines.

Damages -reduced in fruit packing lines is a major cause of grace reduction and quality loos in fresh marks: fruit. Fruit must be treated gently during in sir handling to improve their qualityin order to get a good price in a competitive market. The correct post-hardvest handling in fruit packing lines is a prerequisite to cut down the heavy post-harvest losses. Fruit packing lines must be evaluated, studying their design, the impacts applied to the fruits, the characteristics of the materials, etc. This study establishes the possibility of carrying out modifications and tests in a packing line during a long period of time. For this purpose, an experimental fruit packing line has been designed and located in the Agricultural Engineering Department of the Polythecnic University of Madrid with the aim of improving mechanical devices and fruit handling conditions to minimize damage to fruit. The experimental line consists of several transporting belts, one rollers transporter, one sizer, one elevator, one singularizer, and three trays to receive the calibrated fruit. The line has a length of 6.15 m and a width cf 1.9 m. Movement of the different components is regulated by electric motors with variable velocity electronically controlled. The height of the transfer points is variable and can be easily modified. The experimental line has been calibrated using two instrumented spheres IS 100 (8.8 cm Ø and6.2cm Ø). Average acceleration values obtained in all the transfers of the experimental line lay under 80 g's, although there is big variation for some of them being some values above 100 g's.

State Convergence Theory Applied to a Delayed Teleoperation System

State convergence is a control strategy that was proposed in the early 2000s to ensure stability and transparency in a teleoperation system under specific control gains values. This control strategy has been implemented for a linear system with or without time delay. This paper represents the first attempt at demonstrating, theoretically and experimentantally, that this control strategy can also be applied to a nonlinear teleoperation system with n degrees of freedom and delay in the communication channel. It is assumed that the human operator applies a constant force on the local manipulator during the teleoperation. In addition, the interaction between the remote manipulator and the environment is considered passive. Communication between the local and remote sites is made by means of a communication channel with variable time delay. In this article the theory of Lyapunov-Krasovskii was used to demonstrate that the local-remote teleoperation system is asymptotically stable.

Towards cost reduction in concentrating solar power: innovative design for an efficient fresnel based solar field

Energía termosolar (de concentración) es uno de los nombres que hacen referencia en español al término inglés “concentrating solar power”. Se trata de una tecnología basada en la captura de la potencia térmica de la radiación solar, de forma que permita alcanzar temperaturas capaces de alimentar un ciclo termodinámico convencional (o avanzado); el futuro de esta tecnología depende principalmente de su capacidad para concentrar la radiación solar de manera eficiente y económica. La presente tesis está orientada hacia la resolución de ciertos problemas importantes relacionados con este objetivo. La mencionada necesidad de reducir costes en la concentración de radiación solar directa, asegurando el objetivo termodinámico de calentar un fluido hasta una determinada temperatura, es de vital importancia. Los colectores lineales Fresnel han sido identificados en la literatura científica como una tecnología con gran potencial para alcanzar esta reducción de costes. Dicha tecnología ha sido seleccionada por numerosas razones, entre las que destacan su gran libertad de diseño y su actual estado inmaduro. Con el objetivo de responder a este desafío se desarrollado un detallado estudio de las propiedades ópticas de los colectores lineales Fresnel, para lo cual se han utilizado métodos analíticos y numéricos de manera combinada. En primer lugar, se han usado unos modelos para la predicción de la localización y la irradiación normal directa del sol junto a unas relaciones analíticas desarrolladas para estudiar el efecto de múltiples variables de diseño en la energía incidente sobre los espejos. Del mismo modo, se han obtenido analíticamente los errores debidos al llamado “off-axis aberration”, a la apertura de los rayos reflejados en los espejos y a las sombras y bloqueos entre espejos. Esto ha permitido la comparación de diferentes formas de espejo –planos, circulares o parabólicos–, así como el diseño preliminar de la localización y anchura de los espejos y receptor sin necesidad de costosos métodos numéricos. En segundo lugar, se ha desarrollado un modelo de trazado de rayos de Monte Carlo con el objetivo de comprobar la validez del estudio analítico, pero sobre todo porque este no es preciso en el estudio de la reflexión en espejos. El código desarrollado está específicamente ideado para colectores lineales Fresnel, lo que ha permitido la reducción del tiempo de cálculo en varios órdenes de magnitud en comparación con un programa comercial más general. Esto justifica el desarrollo de un nuevo código en lugar de la compra de una licencia de otro programa. El modelo ha sido usado primeramente para comparar la intensidad de flujo térmico y rendimiento de colectores Fresnel, con y sin reflector secundario, con los colectores cilíndrico parabólicos. Finalmente, la conjunción de los resultados obtenidos en el estudio analítico con el programa numérico ha sido usada para optimizar el campo solar para diferentes orientaciones –Norte-Sur y Este-Oeste–, diferentes localizaciones –Almería y Aswan–, diferentes inclinaciones hacia el Trópico –desde 0 deg hasta 32 deg– y diferentes mínimos de intensidad del flujo en el centro del receptor –10 kW/m2 y 25 kW/m2–. La presente tesis ha conducido a importantes descubrimientos que deben ser considerados a la hora de diseñar un campo solar Fresnel. En primer lugar, los espejos utilizados no deben ser plano, sino cilíndricos o parabólicos, ya que los espejos curvos implican mayores concentraciones y rendimiento. Por otro lado, se ha llegado a la conclusión de que la orientación Este-Oeste es más propicia para localizaciones con altas latitudes, como Almería, mientras que en zonas más cercanas a los trópicos como Aswan los campos Norte-Sur conducen a mayores rendimientos. Es de destacar que la orientación Este-Oeste requiere aproximadamente la mitad de espejos que los campos Norte-Sur, puediendo estar inclinados hacia los Trópicos para mejorar el rendimiento, y que alcanzan parecidos valores de intensidad térmica en el receptor todos los días a mediodía. Sin embargo, los campos con orientación Norte-Sur permiten un flujo más constante a lo largo de un día. Por último, ha sido demostrado que el uso de diseños pre-optimizados analíticamente, con anchura de espejos y espaciado entre espejos variables a lo ancho del campo, pueden implicar aumentos de la energía generada por metro cuadrado de espejos de hasta el 6%. El rendimiento óptico anual de los colectores cilíndrico parabólicos es 23 % mayor que el rendimiento de los campos Fresnel en Almería, mientras que la diferencia es de solo 9 % en Aswan. Ello implica que, para alcanzar el mismo precio de electricidad que la tecnología de referencia, la reducción de costes de instalación por metro cuadrado de espejo debe estar entre el 10 % y el 25 %, y que los colectores lineales Fresnel tienen más posibilidades de ser desarrollados en zonas de bajas latitudes. Como consecuencia de los estudios desarrollados en esta tesis se ha patentado un sistema de almacenamiento que tiene en cuenta la variación del flujo térmico en el receptor a lo largo del día, especialmente para campos con orientación Este-Oeste. Este invento permitiría el aprovechamiento de la energía incidente durante más parte del año, aumentando de manera apreciable los rendimientos óptico y térmico. Abstract Concentrating solar power is the common name of a technology based on capturing the thermal power of solar radiation, in a suitable way to reach temperatures able to activate a conventional (or advanced) thermodynamic cycle to generate electricity; this quest mainly depends on our ability to concentrate solar radiation in a cheap and efficient way. The present thesis is focused to highlight and help solving some of the important issues related to this problem. The need of reducing costs in concentrating the direct solar radiation, but without jeopardizing the thermodynamic objective of heating a fluid up to the required temperature, is of prime importance. Linear Fresnel collectors have been identified in the scientific literature as a technology with high potential to reach this cost reduction. This technology has been selected because of a number of reasons, particularly the degrees of freedom of this type of concentrating configuration and its current immature state. In order to respond to this challenge, a very detailed exercise has been carried out on the optical properties of linear Fresnel collectors. This has been done combining analytic and numerical methods. First, the effect of the design variables on the ratio of energy impinging onto the reflecting surface has been studied using analytically developed equations, together with models that predict the location and direct normal irradiance of the sun at any moment. Similarly, errors due to off-axis aberration, to the aperture of the reflected energy beam and to shading and blocking effects have been obtained analytically. This has allowed the comparison of different shapes of mirrors –flat, cylindrical or parabolic–, as well as a preliminary optimization of the location and width of mirrors and receiver with no need of time-consuming numerical models. Second, in order to prove the validity of the analytic results, but also due to the fact that the study of the reflection process is not precise enough when using analytic equations, a Monte Carlo Ray Trace model has been developed. The developed code is designed specifically for linear Fresnel collectors, which has reduced the computing time by several orders of magnitude compared to a wider commercial software. This justifies the development of the new code. The model has been first used to compare radiation flux intensities and efficiencies of linear Fresnel collectors, both multitube receiver and secondary reflector receiver technologies, with parabolic trough collectors. Finally, the results obtained in the analytic study together with the numeric model have used in order to optimize the solar field for different orientations –North-South and East-West–, different locations –Almería and Aswan–, different tilts of the field towards the Tropic –from 0 deg to 32 deg– and different flux intensity minimum requirements –10 kW/m2 and 25 kW/m2. This thesis work has led to several important findings that should be considered in the design of Fresnel solar fields. First, flat mirrors should not be used in any case, as cylindrical and parabolic mirrors lead to higher flux intensities and efficiencies. Second, it has been concluded that, in locations relatively far from the Tropics such as Almería, East-West embodiments are more efficient, while in Aswan North- South orientation leads to a higher annual efficiency. It must be noted that East-West oriented solar fields require approximately half the number of mirrors than NS oriented fields, can be tilted towards the Equator in order to increase the efficiency and attain similar values of flux intensity at the receiver every day at midday. On the other hand, in NS embodiments the flux intensity is more even during each single day. Finally, it has been proved that the use of analytic designs with variable shift between mirrors and variable width of mirrors across the field can lead to improvements in the electricity generated per reflecting surface square meter up to 6%. The annual optical efficiency of parabolic troughs has been found to be 23% higher than the efficiency of Fresnel fields in Almería, but it is only around 9% higher in Aswan. This implies that, in order to attain the same levelized cost of electricity than parabolic troughs, the required reduction of installation costs per mirror square meter is in the range of 10-25%. Also, it is concluded that linear Fresnel collectors are more suitable for low latitude areas. As a consequence of the studies carried out in this thesis, an innovative storage system has been patented. This system takes into account the variation of the flux intensity along the day, especially for East-West oriented solar fields. As a result, the invention would allow to exploit the impinging radiation along longer time every day, increasing appreciably the optical and thermal efficiencies.

Truncation error estimation in the Discontinuous Galerkin Spectral Element Method

En esta tesis, el método de estimación de error de truncación conocido como restimation ha sido extendido de esquemas de bajo orden a esquemas de alto orden. La mayoría de los trabajos en la bibliografía utilizan soluciones convergidas en mallas de distinto refinamiento para realizar la estimación. En este trabajo se utiliza una solución en una única malla con distintos órdenes polinómicos. Además, no se requiere que esta solución esté completamente convergida, resultando en el método conocido como quasi-a priori T-estimation. La aproximación quasi-a priori estima el error mientras el residuo del método iterativo no es despreciable. En este trabajo se demuestra que algunas de las hipótesis fundamentales sobre el comportamiento del error, establecidas para métodos de bajo orden, dejan de ser válidas en esquemas de alto orden, haciendo necesaria una revisión completa del comportamiento del error antes de redefinir el algoritmo. Para facilitar esta tarea, en una primera etapa se considera el método conocido como Chebyshev Collocation, limitando la aplicación a geometrías simples. La extensión al método Discontinuouos Galerkin Spectral Element Method presenta dificultades adicionales para la definición precisa y la estimación del error, debidos a la formulación débil, la discretización multidominio y la formulación discontinua. En primer lugar, el análisis se enfoca en leyes de conservación escalares para examinar la precisión de la estimación del error de truncación. Después, la validez del análisis se demuestra para las ecuaciones incompresibles y compresibles de Euler y Navier Stokes. El método de aproximación quasi-a priori r-estimation permite desacoplar las contribuciones superficiales y volumétricas del error de truncación, proveyendo información sobre la anisotropía de las soluciones así como su ratio de convergencia con el orden polinómico. Se demuestra que esta aproximación quasi-a priori produce estimaciones del error de truncación con precisión espectral. ABSTRACT In this thesis, the τ-estimation method to estimate the truncation error is extended from low order to spectral methods. While most works in the literature rely on fully time-converged solutions on grids with different spacing to perform the estimation, only one grid with different polynomial orders is used in this work. Furthermore, a non timeconverged solution is used resulting in the quasi-a priori τ-estimation method. The quasi-a priori approach estimates the error when the residual of the time-iterative method is not negligible. It is shown in this work that some of the fundamental assumptions about error tendency, well established for low order methods, are no longer valid in high order schemes, making necessary a complete revision of the error behavior before redefining the algorithm. To facilitate this task, the Chebyshev Collocation Method is considered as a first step, limiting their application to simple geometries. The extension to the Discontinuous Galerkin Spectral Element Method introduces additional features to the accurate definition and estimation of the error due to the weak formulation, multidomain discretization and the discontinuous formulation. First, the analysis focuses on scalar conservation laws to examine the accuracy of the estimation of the truncation error. Then, the validity of the analysis is shown for the incompressible and compressible Euler and Navier Stokes equations. The developed quasi-a priori τ-estimation method permits one to decouple the interfacial and the interior contributions of the truncation error in the Discontinuous Galerkin Spectral Element Method, and provides information about the anisotropy of the solution, as well as its rate of convergence in polynomial order. It is demonstrated here that this quasi-a priori approach yields a spectrally accurate estimate of the truncation error.

Modelización del efecto de la pendiente de la vía sobre el consumo de vehículos turismo

El presente trabajo tiene por objetivo generar una metodología validada que permita predecir el consumo de vehículos turismo circulando en cualquier tramo de vía a partir del perfil orográfico y del diagrama velocidad-tiempo. Para la generación de la metodología, se ha realizado un modelo de simulación con el programa ADVISOR que permite calcular el consumo de combustible para un determinado recorrido en el que se tiene en cuenta el perfil orográfico. Este modelo fue validado con datos reales medidos con equipos on-board y se usó para calcular el consumo de combustible diferencial debido al efecto de la pendiente de la vía, al poderse simular con y sin pendiente. Se realizaron múltiples simulaciones de recorridos con velocidad máxima variable con el fin de obtener un número significativo de datos. Con los resultados de las diferentes simulaciones, se realizó un estudio estadístico, para determinar las variables influyentes y se generó una función estadística (Ecuación de Consumo Estimado – ECE) que permite calcular el consumo de combustible debido a la pendiente de la vía, conociendo el consumo del vehículo en carretera llana (sin pendiente). Esta función estadística generada (ECE), se validó con datos reales medidos en tráfico real. Con el fin de darle generalidad y aplicabilidad a la función generada, y teniendo en cuenta que el consumo de combustible en carretera llana no está siempre disponible, se ha calculado el consumo de combustible sin pendiente utilizando la metodología Copert 4, metodología oficial desarrollada por la Agencia de Medio Ambiente de Europa (EEA) para la estimación de emisiones y consumo de combustible que está basada en datos experimentales pero que no tiene en cuenta la pendiente de la vía. La Ecuación de Consumo Estimado (ECE) aplicada a los consumos calculados por la metodología Copert 4, se valida también usando datos reales medidos en tráfico real y se comprueba que esta función se ajusta considerablemente bien a la realidad, con un error en el consumo acumulado frente al del ensayo real de un 1% y una correlación con el consumo instantáneo del ensayo real de 0,93. Se concluye, que la Función de Consumo Estimado (ECE), permite predecir el efecto de la pendiente sobre el consumo de combustible de un vehículo turismo en tráfico real con un error menor del 1%. Abstract This projects aims to develop a validated methodology that enables to predict cars consumption while circulating at any kind of road section based on its orographic outline and the speed-time diagram. In order to develop this methodology, a simulation model has been performed with the programme ADVISOR, that allows fuel consumption calculation for an specific route in which the orographic outline is considered. This model was validated by real data measured with an on-board equipment and it was used to calculate the differential fuel consumption caused by the effect of the slope on the road, as it was possible to simulate with or without slope. Many simulations were run with routes with variable maximum speed, aiming to obtain a significant amount of data. An statistical study was carried out with the results of those simulations with the purpose to determine the influential variables and an statistical function ( Estimated Consumption Equation – ECE) that enables fuel consumption calculation due to the road’s slope when the consumption of a vehicle on horizontal road (without any slope) is known. This statistical function (ECE) was validated by real data measured in real traffic conditions. With the purpose to generalise the function and increase its applicability, considering that the consumption of a vehicle on horizontal road is not always available, the nonslope fuel consumption has been calculated through Copert 4 methodology, which is the official methodology developed by the European Environmental Agency (EEA) for emissions and fuel consumption calculation based on experimental data, but without taking into consideration the road’s slope. The Estimated Consumption Equation (ECE) applied to the consumption calculated through Copert 4 methodology is also validated using real data measured in real traffic conditions. It was verified that this function considerably adjusts to reality, with an error on the accumulated consumption compared to the real test of 1% and a correlation with the real test immediate fuel consumption of 0,93. It is concluded that the Estimated Consumption Equation (ECE) enables to predict the effect of the slope on the fuel consumption of a car in real traffic conditions with an error less than 1%.