Developing discontinuous Galerkin methods for the resolution of the incompressible Navier-Stokes equations

Informe de investigación elaborado a partir de una estancia en el Laboratorio de Diseño Computacional en Aeroespacial en el Massachusetts Institute of Technology (MIT), Estados Unidos, entre noviembre de 2006 y agosto de 2007. La aerodinámica es una rama de la dinámica de fluidos referida al estudio de los movimientos de los líquidos o gases, cuya meta principal es predecir las fuerzas aerodinámicas en un avión o cualquier tipo de vehículo, incluyendo los automóviles. Las ecuaciones de Navier-Stokes representan un estado dinámico del equilibrio de las fuerzas que actúan en cualquier región dada del fluido. Son uno de los sistemas de ecuaciones más útiles porque describen la física de una gran cantidad de fenómenos como corrientes del océano, flujos alrededor de una superficie de sustentación, etc. En el contexto de una tesis doctoral, se está estudiando un flujo viscoso e incompresible, solucionando las ecuaciones de Navier- Stokes incompresibles de una manera eficiente. Durante la estancia en el MIT, se ha utilizado un método de Galerkin discontinuo para solucionar las ecuaciones de Navier-Stokes incompresibles usando, o bien un parámetro de penalti para asegurar la continuidad de los flujos entre elementos, o bien un método de Galerkin discontinuo compacto. Ambos métodos han dado buenos resultados y varios ejemplos numéricos se han simulado para validar el buen comportamiento de los métodos desarrollados. También se han estudiado elementos particulares, los elementos de Raviart y Thomas, que se podrían utilizar en una formulación mixta para obtener un algoritmo eficiente para solucionar problemas numéricos complejos.

Increased flow of fatty acids toward beta-oxidation in developing seeds of Arabidopsis deficient in diacylglycerol acyltransferase activity or synthesizing medium-chain-length fatty acids.

Synthesis of polyhydroxyalkanoates (PHAs) from intermediates of fatty acid beta-oxidation was used as a tool to study fatty acid degradation in developing seeds of Arabidopsis. Transgenic plants expressing a peroxisomal PHA synthase under the control of a napin promoter accumulated PHA in developing seeds to a final level of 0. 06 mg g(-1) dry weight. In plants co-expressing a plastidial acyl-acyl carrier protein thioesterase from Cuphea lanceolata and a peroxisomal PHA synthase, approximately 18-fold more PHA accumulated in developing seeds. The proportion of 3-hydroxydecanoic acid monomer in the PHA was strongly increased, indicating a large flow of capric acid toward beta-oxidation. Furthermore, expression of the peroxisomal PHA synthase in an Arabidopsis mutant deficient in the enzyme diacylglycerol acyltransferase resulted in a 10-fold increase in PHA accumulation in developing seeds. These data indicate that plants can respond to the inadequate incorporation of fatty acids into triacylglycerides by recycling the fatty acids via beta-oxidation and that a considerable flow toward beta-oxidation can occur even in a plant tissue primarily devoted to the accumulation of storage lipids.

Metabolic syndrome according to different definitions in a rapidly developing country of the African region

AIMS: We examined, in a country of the African region, i) the prevalence of the metabolic syndrome (MetS) according to three definitions (ATP, WHO and IDF); ii) the distribution of the MetS criteria; iii) the level of agreement between these three definitions and iv) we also examined these issues upon exclusion of people with diabetes. METHODS: We conducted an examination survey on a sample representative of the general population aged 25-64 years in the Seychelles (Indian Ocean, African region), attended by 1255 participants (participation rate of 80.3%). RESULTS: The prevalence of MetS increased markedly with age. According to the ATP, WHO and IDF definitions, the prevalence of MetS was, respectively, 24.0%, 25.0%, 25.1% in men and 32.2%, 24.6%, 35.4% in women. Approximately 80% of participants with diabetes also had MetS and the prevalence of MetS was approximately 7% lower upon exclusion of diabetic individuals. High blood pressure and adiposity were the criteria found most frequently among MetS holders irrespective of the MetS definitions. Among people with MetS based on any of the three definitions, 78% met both ATP and IDF criteria, 67% both WHO and IDF criteria, 54% both WHO and ATP criteria and only 37% met all three definitions. CONCLUSION: We identified a high prevalence of MetS in this population in epidemiological transition. The prevalence of MetS decreased by approximately 32% upon exclusion of persons with diabetes. Because of limited agreement between the MetS definitions, the fairly similar proportions of MetS based on any of the three MetS definitions classified, to a substantial extent, different subjects as having MetS.

Metabolic syndrome according to different definitions in a rapidly developing country of the African region

Rapport de synthèseLe syndrome métabolique représente un ensemble de facteurs de risque métaboliques souvent présents simultanément et il est associé à un risque accru de développer des maladies cardiovasculaires. La prevalence du syndrome métabolique est à la hausse au niveau mondial comme cela a souvent été documenté, en particulier dans les pays développés. Pourtant, les données concernant le syndrome métabolique dans les pays de la région sub-saharienne restent rares.Au cours des dernières années, plusieurs définitions du syndrome métabolique ont été formulées, dont celle du 'National Cholesterol Education Program Adult Treatment Panel III', celle de 1 Organisation Mondiale de la Santé et celle du 'International Diabetes Federation'. Parmi les controverses au sujet du syndrome métabolique persiste la question de l'utilité de rechercher la présence du syndrome métabolique chez les patients diabétiques, étant donné que la présence d'un diabète en soit suffit pour identifier un individu à haut risque de faire un événement cardiovasculaire.L'objectif de ce travail de thèse a été de déterminer la prévalence du syndrome métabolique selon les trois définitions majeures mentionnées ci-dessus, grâce à une étude de population transversale, réalisée aux Seychelles en 2004 dans un échantillon représentatif de la population âgée de 24-65 ans (n=1255, taux de participation de 80.3%). L'intérêt d'examiner cette question dans ce pays était d'obtenir des informations dans un pays en transition épidémiologique.Les résultats de ce travail montrent que la prévalence du syndrome métabolique aux Seychelles est élevée, quelque soit la définition utilisée. Selon la définition utilisée, cette prévalence était d'environ 25% chez les hommes et variant entre 25 et 35% chez les femmes.Cependant, malgré des prévalences semblables selon ces trois définitions, la concordance entre ces définitions n'était pas bonne, impliquant que ces différentes définitions classifient, à un certain degré, des individus différents comme étant porteurs du syndrome métabolique.En outre, la plupart (environ 80%) des individus diabétiques avaient un syndrome métabolique. Après exclusion des individus diabétiques, la prévalence du syndrome métabolique dans la population est réduite d'environ un tiers, à environ 20-25%.Ces résultats montrent que, d'une part, le fardeau de maladie dû au syndrome métabolique aux Seychelles, un pays en voie de développement, est considérable. Cette observation peut potentiellement s'appliquer à d'autres pays à un stade de développement semblable. Cela renforce le besoin de mettre en oeuvre des stratégies de santé publique afin de cibler les causes de ces désordres métaboliques, tels que le surpoids et la sédentarité. D'un point de vue du diagnostic, les trois définitions du syndrome métabolique semblent classifier un nombre semblable de personnes atteints du syndrome métabolique dans cette population. Par contre, la relativement mauvaise concordance entre ces définitions - certaines personnes identifiés comme porteurs du syndrome métabolique selon une définition ne le sont pas selon une autre - confirme la nécessité de clarifier la signification de ces différentes définitions et/ou éventuellement de développer une définition unifiée et fiable du syndrome métabolique.

GIS-based debris flow source and runout susceptibility assessment from DEM data - a case study in NW Nicaragua

In October 1998, Hurricane Mitch triggered numerous landslides (mainly debris flows) in Honduras and Nicaragua, resulting in a high death toll and in considerable damage to property. The potential application of relatively simple and affordable spatial prediction models for landslide hazard mapping in developing countries was studied. Our attention was focused on a region in NW Nicaragua, one of the most severely hit places during the Mitch event. A landslide map was obtained at 1:10 000 scale in a Geographic Information System (GIS) environment from the interpretation of aerial photographs and detailed field work. In this map the terrain failure zones were distinguished from the areas within the reach of the mobilized materials. A Digital Elevation Model (DEM) with 20 m×20 m of pixel size was also employed in the study area. A comparative analysis of the terrain failures caused by Hurricane Mitch and a selection of 4 terrain factors extracted from the DEM which, contributed to the terrain instability, was carried out. Land propensity to failure was determined with the aid of a bivariate analysis and GIS tools in a terrain failure susceptibility map. In order to estimate the areas that could be affected by the path or deposition of the mobilized materials, we considered the fact that under intense rainfall events debris flows tend to travel long distances following the maximum slope and merging with the drainage network. Using the TauDEM extension for ArcGIS software we generated automatically flow lines following the maximum slope in the DEM starting from the areas prone to failure in the terrain failure susceptibility map. The areas crossed by the flow lines from each terrain failure susceptibility class correspond to the runout susceptibility classes represented in a runout susceptibility map. The study of terrain failure and runout susceptibility enabled us to obtain a spatial prediction for landslides, which could contribute to landslide risk mitigation.

Trust region versus line search for computing the optical flow

We consider the numerical treatment of the optical flow problem by evaluating the performance of the trust region method versus the line search method. To the best of our knowledge, the trust region method is studied here for the first time for variational optical flow computation. Four different optical flow models are used to test the performance of the proposed algorithm combining linear and nonlinear data terms with quadratic and TV regularization. We show that trust region often performs better than line search; especially in the presence of non-linearity and non-convexity in the model.

Trust region versus line search for computing the optical flow

We consider the numerical treatment of the optical flow problem by evaluating the performance of the trust region method versus the line search method. To the best of our knowledge, the trust region method is studied here for the first time for variational optical flow computation. Four different optical flow models are used to test the performance of the proposed algorithm combining linear and nonlinear data terms with quadratic and TV regularization. We show that trust region often performs better than line search; especially in the presence of non-linearity and non-convexity in the model.

Developing a desertification indicator system for a small Mediterranean catchment: a case study from the Serra de Rodes, Alt Empordà, Catalunya, NE Spain

La desertificació és un problema de degradació de sòls de gran importància en regions àrides, semi-àrides i sub-humides, amb serioses conseqüències ambientals, socials i econòmiques com a resultat de l'impacte d'activitats humanes en combinació amb condicions físiques i medi ambientals desfavorables (UNEP, 1994). L'objectiu principal d'aquesta tesi va ser el desenvolupament d'una metodologia simple per tal de poder avaluar de forma precisa l'estat i l'evolució de la desertificació a escala local, a través de la creació d'un model anomenat sistema d'indicators de desertificació (DIS). En aquest mateix context, un dels dos objectius específics d'aquesta recerca es va centrar en l'estudi dels factors més importants de degradació de sòls a escala de parcel.la, comportant un extens treball de camp, analisi de laboratori i la corresponent interpretació i discussió dels resultats obtinguts. El segon objectiu específic es va basar en el desenvolupament i aplicació del DIS. L'àrea d'estudi seleccionada va ser la conca de la Serra de Rodes, un ambient típic Mediterràni inclòs en el Parc Natural del Cap de Creus, NE Espanya, el qual ha estat progressivament abandonat pels agricultors durant el segle passat. Actualment, els incendis forestals així com el canvi d'ús del sòl i especialment l'abandonament de terres són considerats els problemes ambientals més importants a l'àrea d'estudi (Dunjó et al., 2003). En primer lloc, es va realitzar l'estudi dels processos i causes de la degradació dels sòls a l'àrea d'interés. En base a aquest coneixement, es va dur a terme la identificació i selecció dels indicadors de desertificació més rellevants. Finalment, els indicadors de desertificació seleccionats a escala de conca, incloent l'erosió del sòl i l'escolament superficial, es van integrar en un model espaial de procés. Ja que el sòl és considerat el principal indicador dels processos d'erosió, segons la FAO/UNEP/UNESCO (1979), tant el paisatge original així com els dos escenaris d'ús del sòl desenvolupats, un centrat en el cas hipotétic del pas d'un incendi forestal, i l'altre un paisatge completament cultivat, poden ser ambients classificats sota baixa o moderada degradació. En comparació amb l'escenari original, els dos escenaris creats van revelar uns valors més elevats d'erosió i escolament superficial, i en particular l'escenari cultivat. Per tant, aquests dos hipotètic escenaris no semblen ser una alternativa sostenible vàlida als processos de degradació que es donen a l'àrea d'estudi. No obstant, un ampli ventall d'escenaris alternatius poden ser desenvolupats amb el DIS, tinguent en compte les polítiques d'especial interés per la regió de manera que puguin contribuir a determinar les conseqüències potencials de desertificació derivades d'aquestes polítiques aplicades en aquest escenari tan complexe espaialment. En conclusió, el model desenvolupat sembla ser un sistema força acurat per la identificació de riscs presents i futurs, així com per programar efectivament mesures per combatre la desertificació a escala de conca. No obstant, aquesta primera versió del model presenta varies limitacions i la necessitat de realitzar més recerca en cas de voler desenvolupar una versió futura i millor del DIS.

The impact of targeted observations made during the Greenland Flow Distortion Experiment

The impact of targeted sonde observations on the 1-3 day forecasts for northern Europe is evaluated using the Met Office four-dimensional variational data assimilation scheme and a 24 km gridlength limited-area version of the Unified Model (MetUM). The targeted observations were carried out during February and March 2007 as part of the Greenland Flow Distortion Experiment, using a research aircraft based in Iceland. Sensitive area predictions using either total energy singular vectors or an ensemble transform Kalman filter were used to predict where additional observations should be made to reduce errors in the initial conditions of forecasts for northern Europe. Targeted sonde data was assimilated operationally into the MetUM. Hindcasts show that the impact of the sondes was mixed. Only two out of the five cases showed clear forecast improvement; the maximum forecast improvement seen over the verifying region was approximately 5% of the forecast error 24 hours into the forecast. These two cases are presented in more detail: in the first the improvement propagates into the verification region with a developing polar low; and in the second the improvement is associated with an upper-level trough. The impact of cycling targeted data in the background of the forecast (including the memory of previous targeted observations) is investigated. This is shown to cause a greater forecast impact, but does not necessarily lead to a greater forecast improvement. Finally, the robustness of the results is assessed using a small ensemble of forecasts.

Developing new components for variable flow distribution system modelling in TRNSYS

This study attempts to fill the existing gap in the simulation of variable flow distribution systems through developing new pressure governing components. These components are able to capture the actual ever-changing system performance curve in variable flow distribution systems together with the prediction of controversial issues such as starving, over-flow and the lack of controllability on the flow rate of different branches in a hydronic system. The performance of the proposed components is verified using a case study under design and off-design circumstances. Full integration of the new components within the TRNSYS simulation package is another advantage of this study, which makes it more applicable for designers in both the design and commissioning of hydronic systems.

Heat flow values in the Sea of Okhotsk region

Tectonic structure and anomalous distributions of geophysical fields of the Sea of Okhotsk region are considered; the lack of reliable data on age of the lithosphere beneath basins of various origin in the Sea of Okhotsk is noted. Model calculations based on geological and geophysical data yielded 65 Ma (Cretaceous-Paleocene boundary) age for the Central Okhotsk rise underlain by the continental lithosphere. This estimate agrees with the age (the end of Cretaceous) derived from seismostratigraphic data. A comparative analysis of theoretical and measured heat flows in the Akademii Nauk Rise, underlain by the thinned continental crust, is performed. The analysis points to a higher (by 20%) value of the measured thermal background of the rise, which is consistent with high negative gradient of gravity anomalies in this area. Calculations yielded 36 Ma (Early Oligocene) age and lithosphere thickness of 50 km for the South Okhotsk depression, whose seafloor was formed by processes of back-arc spreading. The estimated age of the depression is supported by kinematic data on the region; the calculated thickness of the lithosphere coincides with the value estimated from data of magnetotelluric sounding here. This indicates that formation time (36 Ma) of the South Okhotsk depression was estimated correctly. Numerical modeling performed for determination of the basement age of rifting basins in the Sea of Okhotsk gave the following estimates: 18 Ma (Early Miocene) for the Deryugin Basin, 12 Ma (Middle Miocene) for the TINRO Basin, and 23 Ma (Late Oligocene) for the West Kamchatka Trough. These estimates agree with formation time (Oligocene-Quaternary) of the sedimentary cover in rifting basins of the Sea of Okhotsk derived from geological and geophysical data. Model temperature estimates are obtained for lithologic and stratigraphic boundaries of the sedimentary cover in the Deryugin and TINRO Basins and the West Kamchatka Trough; the temperature analysis indicates that the latter two structures are promising for oil and hydrocarbon gas generation; the West Kamchatka Trough possesses better reservoir properties compared to the TINRO and Deryugin Basins. The latter is promising for generation of hydrocarbon gas. Paleogeodynamic reconstructions of the Sea of Okhotsk region evolution are obtained for times of 90, 66, and 36 Ma on the base of kinematic, geomagnetic, structural, tectonic, geothermal, and other geological and geophysical data.

Background heat flow values in the region of the Sea of Okhotsk

Particular features of tectonic structure and anomalous distribution of geothermal, geomagnetic, and gravity fields in the region of the Sea of Okhotsk are considered. On the basis of heat flow data, ages of large-scale structures in the Sea of Okhotsk are estimated at 65 Ma for the Central Okhotsk Rise and 36 Ma for the South Okhotsk Basin. Age of the South Okhotsk Basin is confirmed by data on kinematics and corresponds to 50 km thickness of the lithosphere. This is in accordance with thickness value obtained by magnetotelluric soundings. Comparative analysis of model geothermal background and measured heat flow values on the Akademii Nauk Rise is performed. Analysis points to abnormally high (~20%) measured heat flow agrees with high negative gradient of gravity anomalies. Estimates of deep heat flow and basement age of riftogenic basins in the Sea of Okhotsk were carried out in the following areas: Deryugin Basin (18 Ma, Early Miocene), TINRO Basin (12 Ma, Middle Miocene), and West Kamchatka Basin (23 Ma, Late Oligocene). Temperatures at boundaries of the main lithological complexes of the sedimentary cover are calculated and zones of oil and gas generation are defined. On the basis of geothermal, magnetic, structural, and other geological-geophysical data a kinematic model of the region of the Sea of Okhotsk for period of 36 Ma was calculated and constructed.