Avaliação do comportamento de uma ponte metálica centenária sujeita à ação de tráfego ferroviário ligeiro

A presente dissertação tem como objetivo principal a análise numérica do comportamento dinâmico de uma ponte ferroviária, sob ação de tráfego ligeiro ferroviário. Neste contexto são apresentados alguns fundamentos teóricos a ter em conta nestes domínios, visando uma melhor compreensão dos fenómenos existentes no comportamento dinâmico de pontes ferroviárias quando sujeitas ao tráfego. O caso de estudo teve como foco a ponte Luiz I, uma ponte metálica situada sobre o rio Douro, que liga as cidades do Porto e Vila Nova de Gaia, sob ação de tráfego ligeiro ferroviário no seu tabuleiro superior para a condição anterior aos trabalhos de reabilitação e reforço realizados entre 2004 e 2005. Para o efeito foi desenvolvido um modelo numérico de elementos finitos da ponte realizado com recurso ao programa ANSYS, assim como um modelo numérico do veículo do Metro de Lisboa. Com base nestes modelos foram obtidos os parâmetros modais, nomeadamente as frequências naturais e os modos de vibração de toda a estrutura e do veículo. O estudo do comportamento dinâmico da ponte foi realizado por intermédio de uma metodologia de cargas móveis e de interação veículo-estrutura, através da ferramenta computacional Train-Bridge Interaction (TBI). As análises dinâmicas foram efetuadas para a passagem dos veículos de passageiros das redes de Metros do Porto e Lisboa. Nestas análises é estudada a resposta da estrutura em função da variabilidade ao nível da secção transversal, dependência do tramo, influência do veículo, da sua velocidade de circulação e impacto das frequências de vibração estimadas pelo modelo numérico.

Structural assessment of a masonry vault in Portugal

This paper reports on a structural safety assessment and performance evaluation of the upper choir of the Santa Maria de Belém Church in the Jerónimos monastery, Lisbon, one of the most important cultural heritage buildings in Portugal. The possibility of adding a new 20 t organ to the upper choir and its effects on the church structure's response are presented. A refined and a simplified finite-element model is developed to investigate the structure's performance under self-weight and seismic actions. A sensitivity analysis is performed to investigate the effect of masonry mechanical properties and rib cross-sections on the structural response, given the difficulty in accurately obtaining this information. The results show that the safety level of the structure is acceptable, even in the case of adding a heavy new organ.

Seismic evaluation of environmentally degraded FRP-strengthened masonry walls

Fiber Reinforced Polymers (FRPs) have been extensively used for externally bonded reinforcement of masonry structures during the last years. Available information shows that FRPs can significantly improve the seismic performance of masonry elements without altering their structural mass. However, the durability and long-term performance of the strengthened elements are not clearly known yet. Recent experimental results show that environmental conditions can lead to degradation of the bond between FRP and masonry and FRP delaminations. But the effect of these local degradation mechanisms on the global structural response is not studied yet. This paper is therefore aimed at numerically investigating the effect of environmental degradation on the global performance of strengthened masonry walls. The nonlinear behavior of masonry walls strengthened with FRP composites is initially simulated with the aim of a FE package. The adopted numerical modeling strategy is verified by comparison of numerical and experimental results. The model, once validated, is used for investigating the effect of materials and bond degradation on the global behavior and failure modes of strengthened walls. The effect of strengthening scheme on the long-term performance of strengthened walls is also investigated. The degradation data are taken from experimental tests previously performed by the authors. The numerical results show that the effect of local material degradation on the global response of strengthened walls depends on the strengthening schemes, and severity of the environmental conditions. Moreover, environmental induced degradations and FRP delaminations can lead to change of expected failure modes in the strengthened elements. These observations, that are usually neglected at the design stage, can be critical in the long-term performance of strengthened structures.

Methodology for real-time adaptation of tunnels support using the observational method

The observational method in tunnel engineering allows the evaluation in real time of the actual conditions of the ground and to take measures if its behavior deviates considerably from predictions. However, it lacks a consistent and structured methodology to use the monitoring data to adapt the support system in real time. The definition of limit criteria above which adaptation is required are not defined and complex inverse analysis procedures (Rechea et al. 2008, Levasseur et al. 2010, Zentar et al. 2001, Lecampion et al. 2002, Finno and Calvello 2005, Goh 1999, Cui and Pan 2012, Deng et al. 2010, Mathew and Lehane 2013, Sharifzadeh et al. 2012, 2013) may be needed to consistently analyze the problem. In this paper a methodology for the real time adaptation of the support systems during tunneling is presented. In a first step limit criteria for displacements and stresses are proposed. The methodology uses graphics that are constructed during the project stage based on parametric calculations to assist in the process and when these graphics are not available, since it is not possible to predict every possible scenario, inverse analysis calculations are carried out. The methodology is applied to the “Bois de Peu” tunnel which is composed by two tubes with over 500 m long. High uncertainty levels existed concerning the heterogeneity of the soil and consequently in the geomechanical design parameters. The methodology was applied in four sections and the results focus on two of them. It is shown that the methodology has potential to be applied in real cases contributing for a consistent approach of a real time adaptation of the support system and highlight the importance of the existence of good quality and specific monitoring data to improve the inverse analysis procedure.

Evaluation of the DPMFoam solver

[Extrat] Multiphase flows are relevant in several industrial processes, thus the availability of accurate numerical modeling tools, able to support the design of products and processes, is of much significance. OpenFOAM version 2.3.x comprises a multiphase flow solver able to couple Eulerian and Lagrangian phases using the discrete particles method (DPM), the DPMFoam. In this work the DPMFoam solver is assessed by comparing its predictions with analytical results and experimental and simulated data available in the literature. They are results from Goldschmidt’s [1] and Hoomans’s [2] theses and the analytical Ergun equation. The goal was to define accuracy and performance of DPMFoam in general scientific or commercial applications. Obtained results demonstrate a good agreement with the reference simulation data and is within reasonable deviations from the experimental values. (...)

An interactive web-based tool for breast reduction surgery simulation

In this work, we present a 3D web-based interactive tool for numerical modeling and simulation approach to breast reduction surgery simulation, to assist surgeons in planning all aspects related to breast reduction surgery before the actual procedure takes place, thereby avoiding unnecessary risks. In particular, it allows the modeling of the initial breast geometry, the definition of all aspects related to the surgery and the visualization of the post-surgery breast shape in a realistic environment.

Hazard assessment of storm events for the battery, New York

The environmental and socio-economic importance of coastal areas is widely recognized, but at present these areas face severe weaknesses and high-risk situations. The increased demand and growing human occupation of coastal zones have greatly contributed to exacerbating such weaknesses. Today, throughout the world, in all countries with coastal regions, episodes of waves overtopping and coastal flooding are frequent. These episodes are usually responsible for property losses and often put human lives at risk. The floods are caused by coastal storms primarily due to the action of very strong winds. The propagation of these storms towards the coast induces high water levels. It is expected that climate change phenomena will contribute to the intensification of coastal storms. In this context, an estimation of coastal flooding hazards is of paramount importance for the planning and management of coastal zones. Consequently, carrying out a series of storm scenarios and analyzing their impacts through numerical modeling is of prime interest to coastal decision-makers. Firstly, throughout this work, historical storm tracks and intensities are characterized for the northeastern region of United States coast, in terms of probability of occurrence. Secondly, several storm events with high potential of occurrence are generated using a specific tool of DelftDashboard interface for Delft3D software. Hydrodynamic models are then used to generate ensemble simulations to assess storms' effects on coastal water levels. For the United States’ northeastern coast, a highly refined regional domain is considered surrounding the area of The Battery, New York, situated in New York Harbor. Based on statistical data of numerical modeling results, a review of the impact of coastal storms to different locations within the study area is performed.

Численное моделирование уменьшения солёности солонцовой почвы

ნიადაგში წყლისა და გახსნილი მარილის ინფილტრაციის განტოლებების რიცხვითი ინტეგრირების საშუალებით მოდელირდება ბიცობიანი ნიადაგის ზედა ფენის მარილიანობის შემცირების პრაქტიკაში გამოყენებული ერთ-ერთი მეთოდი. ნაჩვენებია, რომ მეთოდი ნიადაგის ფილტრაციის და დიფუზიის კოეფიციენტების გარკვეული მნიშვნელობებისა და 90 დღეღამის განმავლობაში სამჯერადი მორწყვის შემთხვევებში იწვევს ნიადაგის ზედა 70 სმ ფენის მარილიანობის შემცირებას.

Perfeccionamiento de las metodologías docentes en la enseñanza de la construcción mixta en grado y postgrado en la Universidad Politécnica de Cataluña (UPC) y la Universidad de Camaguey (UC)

Estudio elaborado a partir de una estancia en la Universitat Polictecnica de Cataluña entre marzo i junio del 2007.El proyecto tiene como objetivo general el perfeccionamiento de las metodologías docentes en la enseñanza de la construcción mixta en grado y postgrado en la Universidad Politécnica de Cataluña (UPC) y la Universidad de Camagüey (UC), en correspondencia con la situación actual y perspectiva de la normativa europea, española y cubana. En el trabajo se realiza una valoración sobre programas experimentales en estructuras mixtas y sus resultados, así como su interrelación con la modelación numérica. Se evalúa la correspondencia entre las normativas europeas, española y cubana y se definen los aspectos específicos que deben ser perfeccionados o complementados en la normatividad vigente en Cuba, así como la estrategia general para acometerlo. Para caracterizar la enseñanza de la construcción mixta en la UPC y la UC, se parte de la valoración de información sobre asignaturas similares en un grupo seleccionado de universidades y sobre el contexto que sirve de marco al perfeccionamiento. Se realiza el diagnostico del sistema de objetivos, conocimientos y organización metodológica en los programas de grado y postgrado y se ofrecen recomendaciones generales para el perfeccionamiento de la enseñanza de la construcción mixta en la UPC y la UC. Finalmente, se realizan propuestas detalladas para ambas instituciones, para el perfeccionamiento de las metodologías docentes, con enfoques generales similares, respetando los criterios específicos y formatos de cada entorno. Se concibe un sistema integrado de actividades para cada asignatura, lográndose que, a pesar de las particularidades de cada contexto, el enfoque de trabajo sea coincidente en gran medida, lo que favorece la continuidad en las relaciones de colaboración entre ambas instituciones.

Metamorphism and kinetics in the Torres del Paine contact aureole

Contact aureoles provide an excellent geologic environment to study the mechanisms of metamorphic reactions in a natural system. The Torres del Paine (TP) intrusion is one of the most spectacular natural laboratories because of its excellent outcrop conditions. It formed in a period from 12.59 to 12.43 Ma and consists of three large granite and four smaller mafic batches. The oldest granite is on top, the youngest at the bottom of the granitic complex, and the granites overly the mafic laccolith. The TP intruded at a depth of 2-3 km into regional metamorphic anchizone to greenschist facies pelites, sandstones, and conglomerates of the Cerro Toro and Punta Barrosa formations. It formed a thin contact aureole of 150-400 m width. This thesis focuses on the reaction kinetics of the mineral cordierite in the contact aureole using quantitative textural analysis methods. First cordierite was formed from chlorite break¬down (zone I, ca. 480 °C, 750 bar). The second cordierite forming reaction was the muscovite break-down, which is accompanied by a modal decrease in biotite and the appearance of k- feldspar (zone II, 540-550 °C, 750 bar). Crystal sizes of the roundish, poikiloblastic cordierites were determined from microscope thin section images by manually marking each crystal. Images were then automatically processed with Matlab. The correction for the intersection probability of each crystal radius yields the crystal size distribution in the rock. Samples from zone I below the laccolith have the largest crystals (0.09 mm). Cordierites from zone II are smaller, with a maximum crystal radius of 0.057 mm. Rocks from zone II have a larger number of small cordierite crystals than rocks from zone I. A combination of these quantitative analysis with numerical modeling of nucleation and growth, is used to infer nucleation and growth parameters which are responsible for the observed mineral textures. For this, the temperature-time paths of the samples need to be known. The thermal history is complex because the main body of the intrusion was formed by several intrusive batches. The emplacement mechanism and duration of each batch can influence the thermal structure in the aureole. A possible subdivision of batches in smaller increments, so called pulses, will focus heat at the side of the intrusion. Focusing all pulses on one side increases the contact aureole size on that side, but decreases it on the other side. It forms a strongly asymmetric contact aureole. Detailed modeling shows that the relative thicknesses of the TP contact aureole above and below the intrusion (150 and 400 m) are best explained by a rapid emplacement of at least the oldest granite batch. Nevertheless, temperatures are significantly too low in all models, compared to observed mineral assemblages in the hornfelses. Hence, an other important thermal mechanisms needs to take place in the host rock. Clastic minerals in the immature sediments outside the contact aureole are hydrated due to small amounts of expelled fluids during contact metamorphism. This leads to a temperature increase of up to 50 °C. The origin of fluids can be traced by stable isotopes. Whole rock stable isotope data (6D and δ180) and chlorine concentrations in biotite document that the TP intrusion induced only very small amounts of fluid flow. Oxygen whole rock data show δ180 values between 9.0 and 10.0 %o within the first 5 m of the contact. Values increase to 13.0 - 15.0 %o further away from the intrusion. Whole rock 6D values display a more complex zoning. First, host rock values (-90 to -70 %o) smoothly decrease towards the contact by ca. 20 %o, up to a distance of ca. 150 m. This is followed by an increase of ca. 20 %o within the innermost 150 m of the aureole (-97.0 to -78 %o at the contact). The initial decrease in 6D values is interpreted to be due to Rayleigh fractionation accompanying the dehydration reactions forming cordierite, while the final increase reflects infiltration of water-rich fluids from the intrusion. An over-estimate on the quantity and the corresponding thermal effect yields a temperature increase of less than 30 °C. This suggests that fluid flow might have contributed only for a small amount to the thermal evolution of the system. A combination of the numerical growth model with the thermal model, including the hydration reaction enthalpies but neglecting fluid flow and incremental growth, can be used to numerically reproduce the observed cordierite textures in the contact aureole. This yields kinetic parameters which indicate fast cordierite crystallization before the thermal peak in the inner aureole, and continued reaction after the thermal peak in the outermost aureole. Only small temperature dependencies of the kinetic parameters seem to be needed to explain the obtained crystal size data. - Les auréoles de contact offrent un cadre géologique privilégié pour l'étude des mécanismes de réactions métamorphiques associés à la mise en place de magmas dans la croûte terrestre. Par ses conditions d'affleurements excellentes, l'intrusion de Torres del Paine représente un site exceptionnel pour améliorer nos connaissances de ces processus. La formation de cette intrusion composée de trois injections granitiques principales et de quatre injections mafiques de volume inférieur couvre une période allant de 12.50 à 12.43 Ma. Le plus vieux granite forme la partie sommitale de l'intrusion alors que l'injection la plus jeune s'observe à la base du complexe granitique; les granites recouvrent la partie mafique du laccolite. L'intrusion du Torres del Paine s'est mise en place a 2-3 km de profondeur dans un encaissant métamorphique. Cet encaissant est caractérisé par un métamorphisme régional de faciès anchizonal à schiste vert et est composé de pélites, de grès, et des conglomérats des formations du Cerro Toro et Punta Barrosa. La mise en place des différentes injections granitiques a généré une auréole de contact de 150-400 m d'épaisseur autour de l'intrusion. Cette thèse se concentre sur la cinétique de réaction associée à la formation de la cordiérite dans les auréoles de contact en utilisant des méthodes quantitatives d'analyses de texture. On observe plusieurs générations de cordiérite dans l'auréole de contact. La première cordiérite est formée par la décomposition de la chlorite (zone I, environ 480 °C, 750 bar), alors qu'une seconde génération de cordiérite est associée à la décomposition de la muscovite, laquelle est accompagnée par une diminution modale de la teneur en biotite et l'apparition de feldspath potassique (zone II, 540-550 °C, 750 bar). Les tailles des cristaux de cordiérites arrondies et blastic ont été déterminées en utilisant des images digitalisées des lames minces et en marquant individuellement chaque cristal. Les images sont ensuite traitées automatiquement à l'aide du programme Matlab. La correction de la probabilité d'intersection en fonction du rayon des cristaux permet de déterminer la distribution de la taille des cristaux dans la roche. Les échantillons de la zone I, en dessous du lacolite, sont caractérisés par de relativement grands cristaux (0.09 mm). Les cristaux de cordiérite de la zone II sont plus petits, avec un rayon maximal de 0.057 mm. Les roches de la zone II présentent un plus grand nombre de petits cristaux de cordiérite que les roches de la zone I. Une combinaison de ces analyses quantitatives avec un modèle numérique de nucléation et croissance a été utilisée pour déduire les paramètres de nucléation et croissance contrôlant les différentes textures minérales observées. Pour développer le modèle de nucléation et de croissance, il est nécessaire de connaître le chemin température - temps des échantillons. L'histoire thermique est complexe parce que l'intrusion est produite par plusieurs injections successives. En effet, le mécanisme d'emplace¬ment et la durée de chaque injection peuvent influencer la structure thermique dans l'auréole. Une subdivision des injections en plus petits incréments, appelés puises, permet de concentrer la chaleur dans les bords de l'intrusion. Une mise en place préférentielle de ces puises sur un côté de l'intrusion modifie l'apport thermique et influence la taille de l'auréole de contact produite, auréole qui devient asymétrique. Dans le cas de la première injection de granite, une modélisation détaillée montre que l'épaisseur relative de l'auréole de contact de Torres del Paine au-dessus et en dessous de l'intrusion (150 et 400 m) est mieux expliquée par un emplacement rapide du granite. Néanmoins, les températures calculées dans l'auréole de con¬tact sont trop basses pour que les modèles thermiques soient cohérants par rapport à la taille de cette auréole. Ainsi, un autre mecanisme exothermique est nécessaire pour permettre à la roche encais¬sante de produire les assemblages observés. L'observation des roches encaissantes entourant les granites montre que les minéraux clastiques dans les sédiments immatures au-dehors de l'auréole sont hydratés suite à la petite quantité de fluide expulsée durant le métamorphisme de contact et/ou la mise en place des granites. Les réactions d'hydratation peuvent permettre une augmentation de la température jusqu'à 50 °C. Afin de déterminer l'origine des fluides, une étude isotopique de roches de l'auréole de contact a été entreprise. Les isotopes stables d'oxygène et d'hydrogène sur la roche totale ainsi que la concentration en chlore dans la biotite indiquent que la mise en place des granites du Torres del Paine n'induit qu'une circulation de fluide limitée. Les données d'oxygène sur roche totale montrent des valeurs δ180 entre 9.0 et 10.0%o au sein des cinq premiers mètres du contact. Les valeurs augmentent jusqu'à 13.0 - 15.0 plus on s'éloigne de l'intrusion. Les valeurs 5D sur roche totale montrent une zonation plus complexe. Les valeurs de la roche encaissante (-90 à -70%o) diminuent progressivement d'environ 20%o depuis l'extérieur de l'auréole jusqu'à une distance d'environ 150 m du granite. Cette diminution est suivie par une augmentation d'environ 20%o au sein des 150 mètres les plus proches du contact (-97.0 à -78%o au contact). La diminution initiale des valeurs de 6D est interprétée comme la conséquence du fractionnement de Rayleigh qui accompagne les réactions de déshydratation formant la cordiérite, alors que l'augmentation finale reflète l'infiltration de fluide riche en eau venant de l'intrusion. A partir de ces résultats, le volume du fluide issu du granite ainsi que son effet thermique a pu être estimé. Ces résultats montrent que l'augmentation de température associée à ces fluides est limitée à un maximum de 30 °C. La contribution de ces fluides dans le bilan thermique est donc faible. Ces différents résultats nous ont permis de créer un modèle thermique associé à la for¬mation de l'auréole de contact qui intègre la mise en place rapide du granite et les réactions d'hydratation lors du métamorphisme. L'intégration de ce modèle thermique dans le modèle numérique de croissance minérale nous permet de calculer les textures des cordiérites. Cepen¬dant, ce modèle est dépendant de la vitesse de croissance et de nucléation de ces cordiérites. Nous avons obtenu ces paramètres en comparant les textures prédites par le modèle et les textures observées dans les roches de l'auréole de contact du Torres del Paine. Les paramètres cinétiques extraits du modèle optimisé indiquent une cristallisation rapide de la cordiérite avant le pic thermique dans la partie interne de l'auréole, et une réaction continue après le pic thermique dans la partie la plus externe de l'auréole. Seules de petites dépendances de température des paramètres de cinétique semblent être nécessaires pour expliquer les don¬nées obtenues sur la distribution des tailles de cristaux. Ces résultats apportent un éclairage nouveau sur la cinétique qui contrôle les réactions métamorphiques.

Structure et évolution thermique de la croûte supérieure des chaînes de subduction : approches thermologique et modélisation numérique, dans les North Cascades (U.S.A) et de la zone de faille Motagua (Guatemala)

Cette thèse cible l'étude de la structure thermique de la croûte supérieure (<10km) dans les arcs magmatiques continentaux, et son influence sur l'enregistrement thermochronologique de leur exhumation et de leur évolution topographique. Nous portons notre regard sur deux chaînes de montagne appartenant aux Cordillères Américaines : Les Cascades Nord (USA) et la zone de faille Motagua (Guatemala). L'approche utilisée est axée sur la thermochronologie (U-Th-Sm)/He sur apatite et zircon, couplée avec la modélisation numérique de la structure thermique de la croûte. Nous mettons en évidence la variabilité à la fois spatiale et temporelle du gradient géothermique, et attirons l'attention du lecteur sur l'importance de prendre en compte la multitude des processus géologiques perturbant la structure thermique dans les chaînes de type cordillère, c'est à dire formées lors de la subduction océanique sous un continent.Une nouvelle approche est ainsi développée pour étudier et contraindre la perturbation thermique autour des chambres magmatiques. Deux profiles âge-elevation (U-Th-Sm)/He sur apatite et zircon, ont été collectées 7 km au sud du batholithe de Chilliwack, Cascades Nord. Les résultats montrent une variabilité spatiale et temporelle du gradient géothermique lors de l'emplacement magmatique qui peut être contrainte et séparé de l'exhumation. Durant l'emplacement de l'intrusion, la perturbation thermique y atteint un état d'équilibre (-80-100 °C/km) qui est fonction du flux de magma et de ia distance à la source du magma, puis rejoint 40 °C/km à la fin du processus d'emplacement magmatique.Quelques nouvelles données (U-Th)/He, replacées dans une compilation des données existantes dans les Cascades Nord, indiquent une vitesse d'exhumation constante (-100 m/Ma) dans le temps et l'espace entre 35 Ma et 2 Ma, associée à un soulèvement uniforme de la chaîne contrôlé par l'emplacement de magma dans la croûte durant toute l'activité de l'arc. Par contre, après ~2 Ma, le versant humide de la chaîne est affecté par une accélération des taux d'exhumation, jusqu'à 3 km de croûte y sont érodés. Les glaciations ont un triple effet sur l'érosion de cette chaîne: (1) augmentation des vitesses d'érosion, d'exhumation et de soulèvement la où les précipitations sont suffisantes, (2) limitation de l'altitude contrôlé par la position de Γ Ε LA, (3) élargissement du versant humide et contraction du versant aride de la chaîne.Les modifications des réseaux de drainage sont des processus de surface souvent sous-estimés au profil d'événements climatiques ou tectoniques. Nous proposons une nouvelle approche couplant une analyse géomorphologique, des données thermochronologiques de basse température ((U-Th-Sm)/He sur apatite et zircon), et l'utilisation de modélisation numérique thermo-cinématique pour les mettre en évidence et les dater; nous testons cette approche sur la gorge de la Skagit river dans les North Cascades.De nouvelles données (U-Th)/He sur zircons, complétant les données existantes, montrent que le déplacement horizontal le long de la faille transformante continentale Motagua, la limite des plaques Caraïbe/Amérique du Nord, a juxtaposé un bloc froid, le bloc Maya (s.s.), contre un bloque chaud, le bloc Chortis (s.s.) originellement en position d'arc. En plus de donner des gammes d'âges thermochronologiques très différents des deux côtés de la faille, le déplacement horizontal rapide (~2 cm/a) a produit un fort échange thermique latéral, résultant en un réchauffement du côté froid et un refroidissement du côté chaud de la zone de faille de Motagua.Enfin des données (U-Th-Sm)/He sur apatite témoignent d'un refroidissement Oligocène enregistré uniquement dans la croûte supérieure de la bordure nord de la zone de faille Motagua. Nous tenterons ultérieurement de reproduire ce découplage vertical de la structure thermique par la modélisation de la formation d'un bassin transtensif et de circulation de fluides le long de la faille de Motagua. - This thesis focuses on the influence of the dynamic thermal structure of the upper crust (<10km) on the thermochronologic record of the exhumational and topographic history of magmatic continental arcs. Two mountain belts from the American Cordillera are studied: the North Cascades (USA) and the Motagua fault zone (Guatemala). I use a combined approach coupling apatite and zircon (U-Th-Sm}/He thermochronology and thermo- kinematic numerical modelling. This study highlights the temporal and spatial variability of the geothermal gradient and the importance to take into account the different geological processes that perturb the thermal structure of Cordilleran-type mountain belts (i.e. mountain belts related to oceanic subduction underneath a continent}.We integrate apatite and zircon (U-Th)/He data with numerical thermo-kinematic models to study the relative effects of magmatic and surface processes on the thermal evolution of the crust and cooling patterns in the Cenozoic North Cascades arc (Washington State, USA). Two age-elevation profiles that are located 7 km south of the well-studied Chiliiwack intrusions shows that spatial and temporal variability in geothermal gradients linked to magma emplacement can be contrained and separated from exhumation processes. During Chiliiwack batholith emplacement at -35-20 Ma, the geothermal gradient of the country rocks increased to a very high steady-state value (80-100°C/km), which is likely a function of magma flux and the distance from the magma source area. Including temporally varying geothermal gradients in the analysis allows quantifying the thermal perturbation around magmatic intrusions and retrieving a relatively simple denudation history from the data.The synthesis of new and previously published (U-Th)/He data reveals that denudation of the Northern Cascades is spatially and temporally constant at -100 m/Ma between ~32 and ~2 Ma, which likely reflects uplift due to magmatic crustal thickening since the initiation of the Cenozoic stage of the continental magmatic arc. In contrast, the humid flank of the North Cascades is affected by a ten-fold acceleration in exhumation rate at ~2 Ma, which we interpret as forced by the initiation of glaciations; around 3 km of crust have been eroded since that time. Glaciations have three distinct effects on the dynamics of this mountain range: (1) they increase erosion, exhumation and uplift rates where precipitation rates are sufficient to drive efficient glacial erosion; (2) they efficiently limit the elevation of the range; (3) they lead to widening of the humid flank and contraction of the arid flank of the belt.Drainage reorganizations constitute an important agent of landscape evolution that is often underestimated to the benefit of tectonic or climatic events. We propose a new method that integrates geomorphology, low-temperature thermochronometry (apatite and zircon {U-Th-Sm)/He), and 3D numerical thermal-kinematic modelling to detect and date drainage instability producing recent gorge incision, and apply this approach to the Skagit River Gorge, North Cascades.Two zircon (U-Th)/He age-elevation profiles sampled on both sides of the Motagua Fault Zone (MFZ), the boundary between the North American and the Caribbean plates, combined with published thermochronological data show that strike-slip displacement has juxtaposed the cold Maya block (s.s.) against the hot, arc derived, Chortis block (s.s ), producing different age patterns on both sides of the fault and short-wavelength lateral thermal exchange, resulting in recent heating of the cool side and cooling of the hot side of the MFZ.Finally, an apatite (U-Th-Sm)/He age-elevation profile records rapid cooling at -35 Ma localized only in the upper crust along the northern side of the Motagua fault zone. We will try to reproduce these data by modeling the thermal perturbation resulting from the formation of a transtensional basin and of fluid flow activity along a crustal- scale strike-slip fault.

Dynamic compact thermal models with multiple power sources: application to an ultrathin chip stacking technology

Whereas numerical modeling using finite-element methods (FEM) can provide transient temperature distribution in the component with enough accuracy, it is of the most importance the development of compact dynamic thermal models that can be used for electrothermal simulation. While in most cases single power sources are considered, here we focus on the simultaneous presence of multiple sources. The thermal model will be in the form of a thermal impedance matrix containing the thermal impedance transfer functions between two arbitrary ports. Eachindividual transfer function element ( ) is obtained from the analysis of the thermal temperature transient at node ¿ ¿ after a power step at node ¿ .¿ Different options for multiexponential transient analysis are detailed and compared. Among the options explored, small thermal models can be obtained by constrained nonlinear least squares (NLSQ) methods if the order is selected properly using validation signals. The methods are applied to the extraction of dynamic compact thermal models for a new ultrathin chip stack technology (UTCS).

Design Procedures and Field Monitoring of Submerged Barbs for Streambank Protection, June 2007

The main objective of this study was to evaluate the hydraulic performance of riprap spurs and weirs in controlling bank erosion at the Southern part of the Raccoon River upstream U.S. Highway 169 Bridge utilizing the commercially available model FESWMS and field monitoring. It was found based on a 2 year monitoring and numerical modeling that the design of structures was overall successful, including their spacing and stability. The riprap material incorporated into the structures was directly and favorably correlated to the flow transmission through the structure, or in other words, dictated the permeable nature of the structure. It was found that the permeable dikes and weirs chosen in this study created less volume of scour in the vicinity of the structure toes and thus have less risk comparatively to other impermeable structures to collapse. The fact that the structures permitted the transmission of flow through them it allowed fine sand particles to fill in the gaps of the rock interstices and thus cement and better stabilize the structures. During bank-full flows the maximum scour hole was recorded away from the structures toe and the scourhole size was directly related to the protrusion angle of the structure to the flow. It was concluded that the proposed structure inclination with respect to the main flow direction was appropriate since it provides maximum bank protection while creating the largest volume of local scour away from the structure and towards the center of the channel. Furthermore, the lowest potential for bank erosion also occurs with the present set-up design chosen by the IDOT. About 2 ft of new material was deposited in the area located between the structures for the period extending from the construction day to May 2007. Surveys obtained by sonar and the presence of vegetation indicate that new material has been added at the bank toes. Finally, the structures provided higher variability in bed topography forming resting pools, creating flow shade on the leeward side of the structure, and separation of bed substrate due to different flow conditions. Another notable environmental benefit to rock riprap weirs and dikes is the creation of resting pools, especially in year 2007 (2nd year of the project). The magnitude of these benefits to aquatic habitat has been found in the literature that is directly related to the induced scour-hole volume.

Patterns of gravity induced aggregate migration during casting of fluid concretes

In this paper, aggregate migration patterns during fluid concrete castings are studied through experiments, dimensionless approach and numerical modeling. The experimental results obtained on two beams show that gravity induced migration is primarily affecting the coarsest aggregates resulting in a decrease of coarse aggregates volume fraction with the horizontal distance from the pouring point and in a puzzling vertical multi-layer structure. The origin of this multi layer structure is discussed and analyzed with the help of numerical simulations of free surface flow. Our results suggest that it finds its origin in the non Newtonian nature of fresh concrete and that increasing casting rate shall decrease the magnitude of gravity induced particle migration. (C) 2012 Elsevier Ltd. All rights reserved.

Development of Self-Cleaning Box Culvert Design – Phase II, TR-619, 2013

Culverts are common means to convey flow through the roadway system for small streams. In general, larger flows and road embankment heights entail the use of multibarrel culverts (a.k.a. multi-box) culverts. Box culverts are generally designed to handle events with a 50-year return period, and therefore convey considerably lower flows much of the time. While there are no issues with conveying high flows, many multi-box culverts in Iowa pose a significant problem related to sedimentation. The highly erosive Iowa soils can easily lead to the situation that some of the barrels can silt-in early after their construction, becoming partially filled with sediment in few years. Silting can reduce considerably the capacity of the culvert to handle larger flow events. Phase I of this Iowa Highway Research Board project (TR-545) led to an innovative solution for preventing sedimentation. The solution was comprehensively investigated through laboratory experiments and numerical modeling aimed at screening design alternatives and testing their hydraulic and sediment conveyance performance. Following this study phase, the Technical Advisory Committee suggested to implement the recommended sediment mitigation design to a field site. The site selected for implementation was a 3-box culvert crossing Willow Creek on IA Hwy 1W in Iowa City. The culvert was constructed in 1981 and the first cleanup was needed in 2000. Phase II of the TR 545 entailed the monitoring of the site with and without the selfcleaning sedimentation structure in place (similarly with the study conducted in laboratory). The first monitoring stage (Sept 2010 to December 2012) was aimed at providing a baseline for the operation of the as-designed culvert. In order to support Phase II research, a cleanup of the IA Hwy 1W culvert was conducted in September 2011. Subsequently, a monitoring program was initiated to document the sedimentation produced by individual and multiple storms propagating through the culvert. The first two years of monitoring showed inception of the sedimentation in the first spring following the cleanup. Sedimentation continued to increase throughout the monitoring program following the depositional patterns observed in the laboratory tests and those documented in the pre-cleaning surveys. The second part of Phase II of the study was aimed at monitoring the constructed self-cleaning structure. Since its construction in December 2012, the culvert site was continuously monitored through systematic observations. The evidence garnered in this phase of the study demonstrates the good performance of the self-cleaning structure in mitigating the sediment deposition at culverts. Besides their beneficial role in sediment mitigation, the designed self-cleaning structures maintain a clean and clear area upstream the culvert, keep a healthy flow through the central barrel offering hydraulic and aquatic habitat similar with that in the undisturbed stream reaches upstream and downstream the culvert. It can be concluded that the proposed self-cleaning structural solution “streamlines” the area upstream the culvert in a way that secures the safety of the culvert structure at high flows while producing much less disturbance in the stream behavior compared with the current constructive approaches.