Geochemistry of Miocene "blue tuff" samples

Distinctive, massive to stratified, pale blue volcaniclastics, initially referred to as the "blue tuff," were encountered at all four sites drilled during ODP Leg 127 in the Japan Sea. Detailed vertical sequence analysis, plagioclase chemistry, plagioclase 87Sr/86Sr isotopic composition, and 40Ar/39Ar age dating indicate that thick sequences of the blue tuff are not genetically related. Blue tuffs at Hole 794B were apparently deposited by density flows at ambient temperature. Deposition was penecontemporaneous with a large submarine phreatomagmatic eruption at 14.9 Ma in bathyal or deeper water depths. The blue tuffs at this location comprise mostly reworked hydroclastic glass shards and lesser amounts of plagioclase crystals. Pyrogenic plagioclase has an average An mole% of 18±3. Comparison of blue tuff plagioclase compositions with the composition of plagioclase from acoustic basement at Site 794 suggests that these rocks are not genetically related. As such, the extrapolation of sediment accumulation rate data in conjunction with this more precise age for the blue tuff corroborates previous minimum age estimates of 16.2 Ma for acoustic basement at Site 794. Blue tuffs at Hole 796B were probably deposited at ambient temperatures by downslope slumping and density flow of reworked pyrogenic debris. This debris includes abundant bubble wall glass shards and plagioclase crystals, with variable admixture of volcanic lithic and intrabasinal fragments. Pyrogenic fragments were produced by subaerial or shallow submarine, magmatic eruptions dated at 7.6 Ma. Blue tuffs contain a heterogeneous mixture of unrelated fragments including a mixed population of plagioclase crystals. The average An mole% of the predominant, probable comagmatic, plagioclase population is 30±4. The two sequences of blue tuff studied are distinct in age, mineral composition, and the eruptive origin of pyroclastic fragments. Preliminary 87Sr/86Sr isotopic compositions of plagioclase, however, indicates that blue tuffs at both locations are the product of typical, subduction-related island arc magmatism. Based on the results of this study, there is no justification for stratigraphic correlation of widespread, Miocene, blue to blue-gray bentonitic tuff and tuffaceous sandstones nor the interpretation that these strata are indicative of regional, explosive submarine volcanism genetically related to rifting and formation of the Japan Sea. Rather, these reworked pyroclastic strata of intermediate composition were deposited over a protracted 6-8 m.y. period in association with widespread, subduction-related submarine to subaerial volcanism in the Japan Sea backarc basin.

Description and chemical composition of manganese nodules during R/V "Tyro" cruise APNAP 1 to the Madeira Abyssal Plain in 1986

The porewater and sediment composition of two boxcores and of a small gravity core, taken on a manganese-nodule-covered hill and in the Madeira Abyssal Plain proper respectively, are compared. The pore-water study of the two boxcores indicates that oxic conditions prevail in both cores. In addition, it indicates that no detectable fluxes of Mn or Fe occur from the porewater to the ocean bottom water. Variations in the geochemical composition of the sediments can be explained by fluctuations in the amount of carbonate, which acts as a diluting agent. A clear carbonate minimum is observed at 20-22 cm depth in the two cores. This minimum is likely to be associated with the last glacial period (10-20 kyr B.P.). This association is supported by the sediment accumulation rate of 15 mm/kyr as found by extrapolation from the rate for pelagic sediments in the Madeira Abyssal Plain. The bulk composition of the manganese nodules recovered from the submarine hill is chemically almost identical to the average composition of Atlantic nodules. The trace metal and Rare Earth Elements composition indicate a hydrogenous origin for the manganese nodules of this study. On the basis of the chemical composition, and that of nodules relative to that of the adjacent sediments, an average nodule accretian rate of 2.8-3.3 mm/myr has been calculated. Although the analyses of the entire ferromanganese nodules that have been studied seem to indicate a homogenous composition, internal structures of the nodules reveal great inhomogeneity, both visually and chemically. These fluctuations may be related to variations in the fluxes of Mn and Fe, which in turn could be climate-related.

U-Th dating of marine sediments from the Bahamas slope

In order to understand the driving forces for Pleistocene climate change more fully we need to compare the timing of climate events with their possible forcing. In contrast to the last interglacial (marine isotope stage (MIS) 5) the timing of the penultimate interglacial (MIS 7) is poorly constrained. This study constrains its timing and structure by precise U-Th dating of high-resolution delta18O records from aragonite-rich Bahamian slope sediments of ODP Leg 166 (Sites 1008 and 1009). The major glacial-interglacial cycles in delta18O are distinct within these cores and some MIS 7 substages can be identified. These sediments are well suited for U-Th dating because they have uranium concentrations of up to 12 ppm and very low initial 230Th contributions with most samples showing 230Th/232Th activity ratio of >75. U and Th concentrations and isotope ratios were measured by thermal ionisation mass spectrometry and multiple collector inductively coupled plasma mass spectrometry, with the latter providing dramatically better precision. Twenty-nine of the 41 samples measured have a delta234U value close to modern seawater suggesting that they have experienced little diagenesis. Ages from 27 of the 41 samples were deemed reliable on the basis of both their U and their Th isotope ratios. Ages generally increase with depth, although we see a repeated section of stratigraphy in one core. Extrapolation of constant sedimentation rate through each substage suggests that the peak of MIS 7e lasted from ~237 to 228 ka and that 7c began at 215 ka. This timing is consistent with existing low precision radiometric dates from speleothem deposits. The beginning of both these substages appears to be slightly later than in orbitally tuned timescales. The end of MIS 7 is complex, but also appears to be somewhat later than is suggested by orbitally tuned timescales, although this event is not particularly well defined in these cores.

Ages and two-way travel times of ODP Sites 180-1109, 180-1115 and 180-1118

Synthetic seismograms are constructed from check shot-corrected velocity and density measurements collected during Ocean Drilling Program (ODP) Leg 180 at Sites 1109, 1115, and 1118. The synthetic seismograms facilitate direct correlation of a coincident multichannel seismic (MCS) profile with borehole data collected at the three sites. The MCS data and the synthetic seismograms correlate very well, with most major reflectors successfully reproduced in the synthetics. Our results enable a direct calibration of the MCS data in terms of age, paleoenvironment, and subsidence history. Seismic reflectors are time correlative within stratigraphic resolution but are often observed to result from different lithologies across strike. Our results facilitate the extrapolation of the sedimentation history into an unsampled section of Site 1118 and enable a full correlation between the three sites using all the data collected during ODP Leg 180. This study forms the foundation for regionalizing the site data to the northern margin of the Woodlark Basin, where the transition from continental rifting to seafloor spreading is taking place.

Isotope oxygen record for ODP sediments from the Exmouth Plateau

A new composite d18O record, generated from calcareous fine-fraction and bulk sediments from the Exmouth Plateau, details long-term Cretaceous climatic change at mid-latitudes in the Southern Hemisphere. Assessment of new and previously published d18O data indicates that a mid-Cretaceous global climatic optimum was achieved sometime between the time of the Cenomanian-Turonian boundary and the middle Turonian, when surface-ocean paleotemperatures were the highest of the past 115 m.y. Periods of cooling and warming that reversed the general patterns were superimposed on long-term Aptian-Turonian warming and Turonian-Maastrichtian cooling trends, respectively. Extrapolation of Southern Hemisphere paleotemperature trends to Maastrichtian paleotemperature data from a low-latitude Pacific guyot implies that maximum mid-Cretaceous low-latitude paleotemperatures could have been in excess of 33°C.

Methane concentrations in sediments of ODP sites 186-1150 and 186-1151

Substantial amounts of adsorbed methane were detected in authigenic carbonate concretions recovered from sedimentary layers from depths between 245 and 1,108 m below seafloor during Ocean Drilling Program Leg 186 to ODP sites 1150 and 1151 on the deep-sea terrace of the Japan Trench. Methane contents were almost two orders of magnitude higher in the concretions (291-4,528 nmol/g wet wt) than in the surrounding bulk sediments (5-93 nmol/g wet wt), whereas methane/ethane ratios and stable carbon isotopic compositions were very similar. Carbonate content of surrounding bulk sediments (0.02-3.2 wet wt%) and methane content of the surrounding bulk sediments correlated positively. Extrapolation of the carbonate contents of bulk sediments suggests that 100 wt% carbonate would correspond to 1,886±732 nmol methane per g bulk sediment, which is similar to the average value observed in the carbonate concretions (1,321±1,067 nmol/g wet wt, n = 13). These data support the hypothesis that, in sediments, adsorbed hydrocarbon gases are strongly associated with authigenic carbonates.

Age determination and carbonate content of sediment core KH-79-3_L3

Laminated sediments deposited under anoxic bottom waters in the Japan Sea during the last glacial maximum (LGM) contain extremely well preserved calcareous microfossils and eolian carbonates. The radiocarbon age-difference between bulk sediment and monospecific planktonic foraminifera in discrete laminae from a core in the southern Japan Sea implies that ~40% of the total carbonates in the sediments at the LGM are of eolian origin. Extrapolation of this result yields a rate of supply of eolian carbonates of ~2800 tons/d to the entire Japan Sea during the LGM. The climatic significance of this flux potentially lies in its broader geographic extension, particularly in the interaction of the carbonate-bearing dust with shallow, corrosive North Pacific waters and with rain in the atmosphere. By increasing the alkalinity of such waters and by enhancing the biological pump the dust flux could have increased CO2 absorption by both the ocean and rain during the LGM.

A magneto-mineralogical study of basalt at DSDP Leg 73 Holes

The basalts in Holes 519A, 522B, and 524 were studied for intensity of natural remanent magnetization, magnetic hysteresis, magnetic susceptibility, stability of isothermal remanence, and thermomagnetic behavior. Some of these properties are sensitive to both the composition and the microstructure of the magnetic minerals, others to composition only. Thus it is possible to separate the two effects and to trace the variation of effective magnetic grain size and degree of alteration within a lithologic unit or over a yet larger distance or time interval. The flow in Hole 519A is highly maghemitized at the top, the degree of maghemitization decreasing with depth in the flow. Effective grain size increases with increasing depth. Electron microprobe analysis of the titanomaghemite grains in these samples provides no support for the leaching out of iron during alteration. The pillows and flows in Hole 522B are distributed among a number of cooling units, and no systematic downhole variations are apparent. The inferred magneto-petrology is consistent with the cooling and alteration history that might be expected within the units. The upper and lower sills in Hole 524 are more uniform and have a larger concentration of well-developed magnetic mineral grains than the pillows and flows in Holes 519A and 522B. Maghemitization appears to have developed from the boundaries of the sills that are in contact with the sediments between the sills.

Pore pressure measurements on IODP Exp308 sites

Overpressures measured with pore pressure penetrometers during Integrated Ocean Drilling Program (IODP) Expedition 308 reach 70% and 60% of the hydrostatic effective stress (View the MathML source) in the first 200 meters below sea floor (mbsf) at Sites U1322 and U1324, respectively, in the deepwater Gulf of Mexico, offshore Louisiana. High overpressures are present within low permeability mudstones where there have been multiple, very large, submarine landslides during the Pleistocene. Beneath 200 mbsf at Site U1324, pore pressures drop significantly: there are no submarine landslides in this mixture of mudstone, siltstone, and sandstone. The penetrometer measurements did not reach the in situ pressure at the end of the deployment. We used a soil model to determine that an extrapolation approach based on the inverse of square route of time (View the MathML source) requires much less decay time to achieve a desirable accuracy than an inverse time (1/t) extrapolation. Expedition 308 examined how rapid and asymmetric sedimentation above a permeable aquifer drives lateral fluid flow, extreme pore pressures, and submarine landslides. We interpret that the high overpressures observed are driven by rapid sedimentation of low permeability material from the ancestral Mississippi River. Reduced overpressure at depth at Site U1324 suggests lateral flow (drainage) whereas high overpressure at Site U1322 requires inflow from below: lateral flow in the underlying permeable aquifer provides one mechanism for these observations. High overpressure near the seafloor reduces slope stability and provides a mechanism for the large submarine landslides and low regional gradient (2°) offshore from the Mississippi delta.

Accumulation rates, carbon composition and Eocene carbonate compensation depths of ODP Sites 199-1218 and 199-1219

CaCO3, Corg, and biogenic SiO2 were measured in Eocene equatorial Pacific sediments from Sites 1218 and 1219, and bulk oxygen and carbon isotopes were measured on selected intervals from Site 1219. These data delineate a series of CaCO3 events that first appeared at ~48 Ma and continued to the Eocene/Oligocene boundary. Each event lasted 1-2 m.y. and is separated from the next by a low CaCO3 interval of a similar time span. The largest of these carbonate accumulation events (CAE-3) is in Magnetochron 18. It began at ~42.2 Ma, lasted until ~40.3 Ma, and was marked by higher than average productivity. The end of CAE-3 was abrupt and was associated with a large-scale carbon transfer to the oceans prior to warming of high-latitude regions. Changes in carbonate compensation depth associated with CAE excursions were small in the early part of the middle Eocene but increased to as much as 800 m by the late middle Eocene before decreasing into the late Eocene. Oxygen isotope data indicate that the carbonate events are associated with cooling conditions and may mark small glaciations in the Eocene.

Hysteresis between coral reef calcification and the seawater aragonite saturation state

Some predictions of how ocean acidification (OA) will affect coral reefs assume a linear functional relationship between the ambient seawater aragonite saturation state (Omega a) and net ecosystem calcification (NEC). We quantified NEC in a healthy coral reef lagoon in the Great Barrier Reef during different times of the day. Our observations revealed a diel hysteresis pattern in the NEC versus Omega a relationship, with peak NEC rates occurring before the Omega a peak and relatively steady nighttime NEC in spite of variable Omega a. Net ecosystem production had stronger correlations with NEC than light, temperature, nutrients, pH, and Omega a. The observed hysteresis may represent an overlooked challenge for predicting the effects of OA on coral reefs. If widespread, the hysteresis could prevent the use of a linear extrapolation to determine critical Omega a threshold levels required to shift coral reefs from a net calcifying to a net dissolving state.

An overview on armor research for the laser fusion project HiPER

During the current preparatory phase of the European laser fusion project HiPER, an intensive effort has being placed to identify an armour material able to protect the internal walls of the chamber against the high thermal loads and high fluxes of x-rays and ions produced during the fusion explosions. This poster addresses the different threats and limitations of a poly-crystalline Tungsten armour. The analysis is carried out under the conditions of an experimental chamber hypothetically constructed to demonstrate laser fusion in a repetitive mode, subjected to a few thousand 48MJ shock ignition shots during its entire lifetime. If compared to the literature, an extrapolation of the thermomechanical and atomistic effects obtained from the simulations of the experimental chamber to the conditions of a Demo reactor (working 24/7 at hundreds of MW) or a future power plant (producing GW) suggests that “standard” tungsten will not be a suitable armour. Thus, new materials based on nano-structured W and C are being investigated as possible candidates. The research programme launched by the HiPER material team is introduced.

Independence in constraint logic programs

Studying independence of literals, variables, and substitutions has proven very useful in the context of logic programming (LP). Here we study independence in the broader context of constraint logic programming (CLP). We show that a naive extrapolation of the LP definitions of independence to CLP is unsatisfactory (in fact, wrong) for two reasons. First, because interaction between variables through constraints is more complex than in the case of logic programming. Second, in order to ensure the efUciency of several optimizations not only must independence of the search space be considered, but also an orthogonal issue - "independence of constraint solving." We clarify these issues by proposing various types of search independence and constraint solver independence, and show how they can be combined to allow different independence-related optimizations, from parallelism to intelligent backtracking. Sufficient conditions for independence which can be evaluated "a-priori" at run-time are also proposed. Our results suggest that independence, provided a suitable definition is chosen, is even more useful in CLP than in LP.

Modelización de costas rocosas acantiladas

El retroceso de las costas acantiladas es un fenómeno muy extendido sobre los litorales rocosos expuestos a la incidencia combinada de los procesos marinos y meteorológicos que se dan en la franja costera. Este fenómeno se revela violentamente como movimientos gravitacionales del terreno esporádicos, pudiendo causar pérdidas materiales y/o humanas. Aunque el conocimiento de estos riesgos de erosión resulta de vital importancia para la correcta gestión de la costa, el desarrollo de modelos predictivos se encuentra limitado desde el punto de vista geomorfológico debido a la complejidad e interacción de los procesos de desarrollo espacio-temporal que tienen lugar en la zona costera. Los modelos de predicción publicados son escasos y con importantes inconvenientes: a) extrapolación, extienden la información de registros históricos; b) empíricos, sobre registros históricos estudian la respuesta al cambio de un parámetro; c) estocásticos, determinan la cadencia y magnitud de los eventos futuros extrapolando las distribuciones de probabilidad extraídas de catálogos históricos; d) proceso-respuesta, de estabilidad y propagación del error inexplorada; e) en Ecuaciones en Derivadas Parciales, computacionalmente costosos y poco exactos. La primera parte de esta tesis detalla las principales características de los modelos más recientes de cada tipo y, para los más habitualmente utilizados, se indican sus rangos de aplicación, ventajas e inconvenientes. Finalmente como síntesis de los procesos más relevantes que contemplan los modelos revisados, se presenta un diagrama conceptual de la recesión costera, donde se recogen los procesos más influyentes que deben ser tenidos en cuenta, a la hora de utilizar o crear un modelo de recesión costera con el objetivo de evaluar la peligrosidad (tiempo/frecuencia) del fenómeno a medio-corto plazo. En esta tesis se desarrolla un modelo de proceso-respuesta de retroceso de acantilados costeros que incorpora el comportamiento geomecánico de materiales cuya resistencia a compresión no supere los 5 MPa. El modelo simula la evolución espaciotemporal de un perfil-2D del acantilado que puede estar formado por materiales heterogéneos. Para ello, se acoplan la dinámica marina: nivel medio del mar, cambios en el nivel medio del lago, mareas y oleaje; con la evolución del terreno: erosión, desprendimiento rocoso y formación de talud de derrubios. El modelo en sus diferentes variantes es capaz de incluir el análisis de la estabilidad geomecánica de los materiales, el efecto de los derrubios presentes al pie del acantilado, el efecto del agua subterránea, la playa, el run-up, cambios en el nivel medio del mar o cambios (estacionales o interanuales) en el nivel medio de la masa de agua (lagos). Se ha estudiado el error de discretización del modelo y su propagación en el tiempo a partir de las soluciones exactas para los dos primeros periodos de marea para diferentes aproximaciones numéricas tanto en tiempo como en espacio. Los resultados obtenidos han permitido justificar las elecciones que minimizan el error y los métodos de aproximación más adecuados para su posterior uso en la modelización. El modelo ha sido validado frente a datos reales en la costa de Holderness, Yorkshire, Reino Unido; y en la costa norte del lago Erie, Ontario, Canadá. Los resultados obtenidos presentan un importante avance en los modelos de recesión costera, especialmente en su relación con las condiciones geomecánicas del medio, la influencia del agua subterránea, la verticalización de los perfiles rocosos y su respuesta ante condiciones variables producidas por el cambio climático (por ejemplo, nivel medio del mar, cambios en los niveles de lago, etc.). The recession of coastal cliffs is a widespread phenomenon on the rocky shores that are exposed to the combined incidence of marine and meteorological processes that occur in the shoreline. This phenomenon is revealed violently and occasionally, as gravitational movements of the ground and can cause material or human losses. Although knowledge of the risks of erosion is vital for the proper management of the coast, the development of cliff erosion predictive models is limited by the complex interactions between environmental processes and material properties over a range of temporal and spatial scales. Published prediction models are scarce and present important drawbacks: extrapolation, that extend historical records to the future; empirical, that based on historical records studies the system response against the change in one parameter; stochastic, that represent of cliff behaviour based on assumptions regarding the magnitude and frequency of events in a probabilistic framework based on historical records; process-response, stability and error propagation unexplored; PDE´s, highly computationally expensive and not very accurate. The first part of this thesis describes the main features of the latest models of each type and, for the most commonly used, their ranges of application, advantages and disadvantages are given. Finally as a synthesis of the most relevant processes that include the revised models, a conceptual diagram of coastal recession is presented. This conceptual model includes the most influential processes that must be taken into account when using or creating a model of coastal recession to evaluate the dangerousness (time/frequency) of the phenomenon to medium-short term. A new process-response coastal recession model developed in this thesis has been designed to incorporate the behavioural and mechanical characteristics of coastal cliffs which are composed of with materials whose compressive strength is less than 5 MPa. The model simulates the spatial and temporal evolution of a cliff-2D profile that can consist of heterogeneous materials. To do so, marine dynamics: mean sea level, waves, tides, lake seasonal changes; is coupled with the evolution of land recession: erosion, cliff face failure and associated protective colluvial wedge. The model in its different variants can include analysis of material geomechanical stability, the effect of debris present at the cliff foot, groundwater effects, beach and run-up effects, changes in the mean sea level or changes (seasonal or inter-annual) in the mean lake level. Computational implementation and study of different numerical resolution techniques, in both time and space approximations, and the produced errors are exposed and analysed for the first two tidal periods. The results obtained in the errors analysis allow us to operate the model with a configuration that minimizes the error of the approximation methods. The model is validated through profile evolution assessment at various locations of coastline retreat on the Holderness Coast, Yorkshire, UK and on the north coast of Lake Erie, Ontario, Canada. The results represent an important stepforward in linking material properties to the processes of cliff recession, in considering the effect of groundwater charge and the slope oversteeping and their response to changing conditions caused by climate change (i.e. sea level, changes in lakes levels, etc.).

Numerical Error Prediction and its applications in CFD using tau-estimation

Nowadays, Computational Fluid Dynamics (CFD) solvers are widely used within the industry to model fluid flow phenomenons. Several fluid flow model equations have been employed in the last decades to simulate and predict forces acting, for example, on different aircraft configurations. Computational time and accuracy are strongly dependent on the fluid flow model equation and the spatial dimension of the problem considered. While simple models based on perfect flows, like panel methods or potential flow models can be very fast to solve, they usually suffer from a poor accuracy in order to simulate real flows (transonic, viscous). On the other hand, more complex models such as the full Navier- Stokes equations provide high fidelity predictions but at a much higher computational cost. Thus, a good compromise between accuracy and computational time has to be fixed for engineering applications. A discretisation technique widely used within the industry is the so-called Finite Volume approach on unstructured meshes. This technique spatially discretises the flow motion equations onto a set of elements which form a mesh, a discrete representation of the continuous domain. Using this approach, for a given flow model equation, the accuracy and computational time mainly depend on the distribution of nodes forming the mesh. Therefore, a good compromise between accuracy and computational time might be obtained by carefully defining the mesh. However, defining an optimal mesh for complex flows and geometries requires a very high level expertize in fluid mechanics and numerical analysis, and in most cases a simple guess of regions of the computational domain which might affect the most the accuracy is impossible. Thus, it is desirable to have an automatized remeshing tool, which is more flexible with unstructured meshes than its structured counterpart. However, adaptive methods currently in use still have an opened question: how to efficiently drive the adaptation ? Pioneering sensors based on flow features generally suffer from a lack of reliability, so in the last decade more effort has been made in developing numerical error-based sensors, like for instance the adjoint-based adaptation sensors. While very efficient at adapting meshes for a given functional output, the latter method is very expensive as it requires to solve a dual set of equations and computes the sensor on an embedded mesh. Therefore, it would be desirable to develop a more affordable numerical error estimation method. The current work aims at estimating the truncation error, which arises when discretising a partial differential equation. These are the higher order terms neglected in the construction of the numerical scheme. The truncation error provides very useful information as it is strongly related to the flow model equation and its discretisation. On one hand, it is a very reliable measure of the quality of the mesh, therefore very useful in order to drive a mesh adaptation procedure. On the other hand, it is strongly linked to the flow model equation, so that a careful estimation actually gives information on how well a given equation is solved, which may be useful in the context of _ -extrapolation or zonal modelling. The following work is organized as follows: Chap. 1 contains a short review of mesh adaptation techniques as well as numerical error prediction. In the first section, Sec. 1.1, the basic refinement strategies are reviewed and the main contribution to structured and unstructured mesh adaptation are presented. Sec. 1.2 introduces the definitions of errors encountered when solving Computational Fluid Dynamics problems and reviews the most common approaches to predict them. Chap. 2 is devoted to the mathematical formulation of truncation error estimation in the context of finite volume methodology, as well as a complete verification procedure. Several features are studied, such as the influence of grid non-uniformities, non-linearity, boundary conditions and non-converged numerical solutions. This verification part has been submitted and accepted for publication in the Journal of Computational Physics. Chap. 3 presents a mesh adaptation algorithm based on truncation error estimates and compares the results to a feature-based and an adjoint-based sensor (in collaboration with Jorge Ponsín, INTA). Two- and three-dimensional cases relevant for validation in the aeronautical industry are considered. This part has been submitted and accepted in the AIAA Journal. An extension to Reynolds Averaged Navier- Stokes equations is also included, where _ -estimation-based mesh adaptation and _ -extrapolation are applied to viscous wing profiles. The latter has been submitted in the Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering. Keywords: mesh adaptation, numerical error prediction, finite volume Hoy en día, la Dinámica de Fluidos Computacional (CFD) es ampliamente utilizada dentro de la industria para obtener información sobre fenómenos fluidos. La Dinámica de Fluidos Computacional considera distintas modelizaciones de las ecuaciones fluidas (Potencial, Euler, Navier-Stokes, etc) para simular y predecir las fuerzas que actúan, por ejemplo, sobre una configuración de aeronave. El tiempo de cálculo y la precisión en la solución depende en gran medida de los modelos utilizados, así como de la dimensión espacial del problema considerado. Mientras que modelos simples basados en flujos perfectos, como modelos de flujos potenciales, se pueden resolver rápidamente, por lo general aducen de una baja precisión a la hora de simular flujos reales (viscosos, transónicos, etc). Por otro lado, modelos más complejos tales como el conjunto de ecuaciones de Navier-Stokes proporcionan predicciones de alta fidelidad, a expensas de un coste computacional mucho más elevado. Por lo tanto, en términos de aplicaciones de ingeniería se debe fijar un buen compromiso entre precisión y tiempo de cálculo. Una técnica de discretización ampliamente utilizada en la industria es el método de los Volúmenes Finitos en mallas no estructuradas. Esta técnica discretiza espacialmente las ecuaciones del movimiento del flujo sobre un conjunto de elementos que forman una malla, una representación discreta del dominio continuo. Utilizando este enfoque, para una ecuación de flujo dado, la precisión y el tiempo computacional dependen principalmente de la distribución de los nodos que forman la malla. Por consiguiente, un buen compromiso entre precisión y tiempo de cálculo se podría obtener definiendo cuidadosamente la malla, concentrando sus elementos en aquellas zonas donde sea estrictamente necesario. Sin embargo, la definición de una malla óptima para corrientes y geometrías complejas requiere un nivel muy alto de experiencia en la mecánica de fluidos y el análisis numérico, así como un conocimiento previo de la solución. Aspecto que en la mayoría de los casos no está disponible. Por tanto, es deseable tener una herramienta que permita adaptar los elementos de malla de forma automática, acorde a la solución fluida (remallado). Esta herramienta es generalmente más flexible en mallas no estructuradas que con su homóloga estructurada. No obstante, los métodos de adaptación actualmente en uso todavía dejan una pregunta abierta: cómo conducir de manera eficiente la adaptación. Sensores pioneros basados en las características del flujo en general, adolecen de una falta de fiabilidad, por lo que en la última década se han realizado grandes esfuerzos en el desarrollo numérico de sensores basados en el error, como por ejemplo los sensores basados en el adjunto. A pesar de ser muy eficientes en la adaptación de mallas para un determinado funcional, este último método resulta muy costoso, pues requiere resolver un doble conjunto de ecuaciones: la solución y su adjunta. Por tanto, es deseable desarrollar un método numérico de estimación de error más asequible. El presente trabajo tiene como objetivo estimar el error local de truncación, que aparece cuando se discretiza una ecuación en derivadas parciales. Estos son los términos de orden superior olvidados en la construcción del esquema numérico. El error de truncación proporciona una información muy útil sobre la solución: es una medida muy fiable de la calidad de la malla, obteniendo información que permite llevar a cabo un procedimiento de adaptación de malla. Está fuertemente relacionado al modelo matemático fluido, de modo que una estimación precisa garantiza la idoneidad de dicho modelo en un campo fluido, lo que puede ser útil en el contexto de modelado zonal. Por último, permite mejorar la precisión de la solución resolviendo un nuevo sistema donde el error local actúa como término fuente (_ -extrapolación). El presenta trabajo se organiza de la siguiente manera: Cap. 1 contiene una breve reseña de las técnicas de adaptación de malla, así como de los métodos de predicción de los errores numéricos. En la primera sección, Sec. 1.1, se examinan las estrategias básicas de refinamiento y se presenta la principal contribución a la adaptación de malla estructurada y no estructurada. Sec 1.2 introduce las definiciones de los errores encontrados en la resolución de problemas de Dinámica Computacional de Fluidos y se examinan los enfoques más comunes para predecirlos. Cap. 2 está dedicado a la formulación matemática de la estimación del error de truncación en el contexto de la metodología de Volúmenes Finitos, así como a un procedimiento de verificación completo. Se estudian varias características que influyen en su estimación: la influencia de la falta de uniformidad de la malla, el efecto de las no linealidades del modelo matemático, diferentes condiciones de contorno y soluciones numéricas no convergidas. Esta parte de verificación ha sido presentada y aceptada para su publicación en el Journal of Computational Physics. Cap. 3 presenta un algoritmo de adaptación de malla basado en la estimación del error de truncación y compara los resultados con sensores de featured-based y adjointbased (en colaboración con Jorge Ponsín del INTA). Se consideran casos en dos y tres dimensiones, relevantes para la validación en la industria aeronáutica. Este trabajo ha sido presentado y aceptado en el AIAA Journal. También se incluye una extensión de estos métodos a las ecuaciones RANS (Reynolds Average Navier- Stokes), en donde adaptación de malla basada en _ y _ -extrapolación son aplicados a perfiles con viscosidad de alas. Este último trabajo se ha presentado en los Actas de la Institución de Ingenieros Mecánicos, Parte G: Journal of Aerospace Engineering. Palabras clave: adaptación de malla, predicción del error numérico, volúmenes finitos