On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

Tight-Binding Approximations in 1D and 2D Coupled-Cavity Photonic Crystal Structures

Light confinement and controlling an optical field has numerous applications in the field of telecommunications for optical signals processing. When the wavelength of the electromagnetic field is on the order of the period of a photonic microstructure, the field undergoes reflection, refraction, and coherent scattering. This produces photonic bandgaps, forbidden frequency regions or spectral stop bands where light cannot exist. Dielectric perturbations that break the perfect periodicity of these structures produce what is analogous to an impurity state in the bandgap of a semiconductor. The defect modes that exist at discrete frequencies within the photonic bandgap are spatially localized about the cavity-defects in the photonic crystal. In this thesis the properties of two tight-binding approximations (TBAs) are investigated in one-dimensional and two-dimensional coupled-cavity photonic crystal structures We require an efficient and simple approach that ensures the continuity of the electromagnetic field across dielectric interfaces in complex structures. In this thesis we develop \textrm{E} -- and \textrm{D} --TBAs to calculate the modes in finite 1D and 2D two-defect coupled-cavity photonic crystal structures. In the \textrm{E} -- and \textrm{D} --TBAs we expand the coupled-cavity \overrightarrow{E} --modes in terms of the individual \overrightarrow{E} -- and \overrightarrow{D} --modes, respectively. We investigate the dependence of the defect modes, their frequencies and quality factors on the relative placement of the defects in the photonic crystal structures. We then elucidate the differences between the two TBA formulations, and describe the conditions under which these formulations may be more robust when encountering a dielectric perturbation. Our 1D analysis showed that the 1D modes were sensitive to the structure geometry. The antisymmetric \textrm{D} mode amplitudes show that the \textrm{D} --TBA did not capture the correct (tangential \overrightarrow{E} --field) boundary conditions. However, the \textrm{D} --TBA did not yield significantly poorer results compared to the \textrm{E} --TBA. Our 2D analysis reveals that the \textrm{E} -- and \textrm{D} --TBAs produced nearly identical mode profiles for every structure. Plots of the relative difference between the \textrm{E} and \textrm{D} mode amplitudes show that the \textrm{D} --TBA did capture the correct (normal \overrightarrow{E} --field) boundary conditions. We found that the 2D TBA CC mode calculations were 125-150 times faster than an FDTD calculation for the same two-defect PCS. Notwithstanding this efficiency, the appropriateness of either TBA was found to depend on the geometry of the structure and the mode(s), i.e. whether or not the mode has a large normal or tangential component.

Análise comparativa dos métodos da ISO 13790 e sua adequabilidade na estimativa das necessidades de energia para aquecimento e arrefecimento e da temperatura do ar interior

Tese de mestrado integrado em Engenharia da Energia e do Ambiente, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2016

Benchmarking numerical models of brittle thrust wedges

We report quantitative results from three brittle thrust wedge experiments, comparing numerical results directly with each other and with corresponding analogue results. We first test whether the participating codes reproduce predictions from analytical critical taper theory. Eleven codes pass the stable wedge test, showing negligible internal deformation and maintaining the initial surface slope upon horizontal translation over a frictional interface. Eight codes participated in the unstable wedge test that examines the evolution of a wedge by thrust formation from a subcritical state to the critical taper geometry. The critical taper is recovered, but the models show two deformation modes characterised by either mainly forward dipping thrusts or a series of thrust pop-ups. We speculate that the two modes are caused by differences in effective basal boundary friction related to different algorithms for modelling boundary friction. The third experiment examines stacking of forward thrusts that are translated upward along a backward thrust. The results of the seven codes that run this experiment show variability in deformation style, number of thrusts, thrust dip angles and surface slope. Overall, our experiments show that numerical models run with different numerical techniques can successfully simulate laboratory brittle thrust wedge models at the cm-scale. In more detail, however, we find that it is challenging to reproduce sandbox-type setups numerically, because of frictional boundary conditions and velocity discontinuities. We recommend that future numerical-analogue comparisons use simple boundary conditions and that the numerical Earth Science community defines a plasticity test to resolve the variability in model shear zones.

Tropically-driven climate shifts in southwestern Europe during MIS 19, a low eccentricity interglacial

The relative roles of high- versus low-latitude forcing of millennial-scale climate variability are still not well understood. Here we present terrestrial–marine climate profiles from the southwestern Iberian margin, a region particularly affected by precession, that show millennial climate oscillations related to a nonlinear response to the Earth's precession cycle during Marine Isotope Stage (MIS) 19. MIS 19 has been considered the best analogue to our present interglacial from an astronomical point of view due to the reduced eccentricity centred at 785 ka. In our records, seven millennial-scale forest contractions punctuated MIS 19 superimposed to two orbitally-driven Mediterranean forest expansions. In contrast to our present interglacial, we evidence for the first time low latitude-driven 5000-yr cycles of drying and cooling in the western Mediterranean region, along with warmth in the subtropical gyre related to the fourth harmonic of precession. These cycles indicate repeated intensification of North Atlantic meridional moisture transport that along with decrease in boreal summer insolation triggered ice growth and may have contributed to the glacial inception, at ∼774 ka. The freshwater fluxes during MIS 19ab amplified the cooling events in the North Atlantic promoting further cooling and leading to MIS 18 glaciation. The discrepancy between the dominant cyclicity observed during MIS 1, 2500-yr, and that of MIS 19, 5000-yr, challenges the similar duration of the Holocene and MIS 19c interglacials under natural boundary conditions.

A 300 kyr record of upwelling off Oman during the Late Quaternary: evidence of the Asian southwest monsoon

In the northwest Arabian Sea upwelling occurs each summer, driven by the strong SW monsoon winds. Upwelling results in high biological productivity and a distinctive assemblage of plankton species in the surface waters off Oman that are preserved in the sediments along the Oman continental margin, creating a geologic record of monsoon-driven upwelling. Sediments recovered from the Oman continental margin during Ocean Drilling Program leg 117 provide an opportunity to examine how upwelling has varied during the late Quaternary, spanning a longer interval than piston cores recovered prior to the ODP cruise. Variations in foraminifer shell accumulation and in the relative abundance of Globigerina bulloides indicate dominant cycles of variation at 1/100 kyr, the dominant frequency of glacial-interglacial variations, and at 1/23 kyr, the frequency of precessionally driven cycles in seasonal insolation. The strongest monsoon winds (indicated by increased upwelling) occurred during interglacial times when perihelion was aligned with the summer solstice, an orbital change that increased the insolation received during summer in the northern hemisphere. During glacial times upwelling was reduced, and although the precessional cycles were still present their amplitude was smaller. At both frequencies the upwelling cycles are in phase with minimum ice volume, evidence that glacial-interglacial climate changes also include changes to the climate system that influence the low-latitude monsoon. We attribute the decrease in the monsoon winds observed during glacial times to changes in bare land albedo over Asia and/or to changes in the areal extent and seasonal cycle in Asian snow cover that decrease the summer land-sea temperature contrast. Other mechanisms may also be involved. These new upwelling time series differ substantially from previous results, however the previous work relied on cores located farther offshore where upwelling is less intense and other physical mechanisms become important. Our results support the observations derived from atmospheric general circulation models of the atmosphere that indicate that both glacial boundary conditions, and the strength of summer insolation are important variables contributing to cycles in the monsoon winds during the late Quaternary.

Data compilation of ESOP

The work in this sub-project of ESOP focuses on the advective and convective transforma-tion of water masses in the Greenland Sea and its neighbouring areas. It includes observational work on the sub-mesoscale and analysis of hydrographic data up to the gyre-scale. Observations of active convective plumes were made with a towed chain equipped with up to 80 CTD sensors, giving a horizontal and vertical resolution of the hydrographic fields of a few metres. The observed scales of the penetrative convective plumes compare well with those given by theory. On the mesoscale the structure of homogeneous eddies formed as a result of deep convection was observed and the associated mixing and renewal of the intermediate layers quantified. The relative importance and efficiency of thermal and haline penetrative convection in relation to the surface boundary conditions (heat and salt fluxes and ice cover) and the ambient stratification are studied using the multi year time series of hydro-graphic data in the central Greenland Sea. The modification of the water column of the Greenland Sea gyre through advection from and mixing with water at its rim is assessed on longer time scales. The relative contributions are quantified using modern water mass analysis methods based on inverse techniques. Likewise the convective renewal and the spreading of the Arctic Intermediate Water from its formation area is quantified. The aim is to budget the heat and salt content of the water column, in particular of the low salinity surface layer, and to relate its seasonal and interannual variability to the lateral fluxes and the fluxes at the air-sea-ice interface. This will allow to estimate residence times for the different layers of the Greenland Sea gyre, a quantity important for the description of the Polar Ocean carbon cycle.

Surface velocity fields, digital elevation models and ice front positions derived from multi-mission SAR remote sensing data at Sjögren Inlet glaciers, Antarctic Peninsula

Substantial retreat or disintegration of numerous ice shelves have been observed on the Antarctic Peninsula. The ice shelf in the Prince Gustav Channel retreated gradually since the late 1980's and broke-up in 1995. Tributary glaciers reacted with speed-up, surface lowering and increased ice discharge, consequently contributing to sea level rise. We present a detailed long-term study (1993-2014) on the dynamic response of Sjögren Inlet glaciers to the disintegration of Prince Gustav Ice Shelf. We analyzed various remote sensing datasets to observe the reactions of the glaciers to the loss of the buttressing ice shelf. A strong increase in ice surface velocities was observed with maximum flow speeds reaching 2.82±0.48 m/d in 2007 and 1.50±0.32 m/d in 2004 at Sjögren and Boydell glaciers respectively. Subsequently, the flow velocities decelerated, however in late 2014, we still measured about two times the values of our first measurements in 1996. The tributary glaciers retreated 61.7±3.1 km² behind the former grounding line of the ice shelf. In regions below 1000 m a.s.l., a mean surface lowering of -68±10 m (-3.1 m/a) was observed in the period 1993-2014. The lowering rate decreased to -2.2 m/a in recent years. Based on the surface lowering rates, geodetic mass balances of the glaciers were derived for different time steps. High mass loss rate of -1.21±0.36 Gt/a was found in the earliest period (1993-2001). Due to the dynamic adjustments of the glaciers to the new boundary conditions the ice mass loss reduced to -0.59±0.11 Gt/a in the period 2012-2014, resulting in an average mass loss rate of -0.89±0.16 Gt/a (1993-2014). Including the retreat of the ice front and grounding line, a total mass change of -38.5±7.7 Gt and a contribution to sea level rise of 0.061±0.013 mm were computed. Analysis of the ice flux revealed that available bedrock elevation estimates at Sjögren Inlet are too shallow and are the major uncertainty in ice flux computations. This temporally dense time series analysis of Sjögren Inlet glaciers shows that the adjustments of tributary glaciers to ice shelf disintegration are still going on and provides detailed information of the changes in glacier dynamics.

(Table 1) Main constituents of coccolith size fractions expressed as carbonate contribution at Bass River during the PETM

Coccoliths, calcite plates produced by the marine phytoplankton coccolithophores, have previously shown a large array of carbon and oxygen stable isotope fractionations (termed "vital effects"), correlated to cell size and hypothesized to reflect the varying importance of active carbon acquisition strategies. Culture studies show a reduced range of vital effects between large and small coccolithophores under high CO2, consistent with previous observations of a smaller range of interspecific vital effects in Paleocene coccoliths. We present new fossil data examining coccolithophore vital effects over three key Cenozoic intervals reflecting changing climate and atmospheric partial pressure of CO2 (pCO2). Oxygen and carbon stable isotopes of size-separated coccolith fractions dominated by different species from well preserved Paleocene-Eocene thermal maximum (PETM, ~56 Ma) samples show reduced interspecific differences within the greenhouse boundary conditions of the PETM. Conversely, isotope data from the Plio-Pleistocene transition (PPT; 3.5-2 Ma) and the last glacial maximum (LGM; ~22 ka) show persistent vital effects of ~2 per mil. PPT and LGM data show a clear positive trend between coccolith (cell) size and isotopic enrichment in coccolith carbonate, as seen in laboratory cultures. On geological timescales, the degree of expression of vital effects in coccoliths appears to be insensitive topCO2 changes over the range ~350 ppm (Pliocene) to ~180 ppm (LGM). The modern array of coccolith vital effects arose after the PETM but before the late Pliocene and may reflect the operation of more diverse carbon acquisition strategies in coccolithophores in response to decreasing Cenozoic pCO2.

Ba/Ca ratios in benthic foraminifera from recent core tops

The carbon isotope ratio (delta13C) and cadmium content (Cd/Ca) of benthic foraminifera shells have been used to reconstruct deep-water circulation patterns of the glacial oceans. These tracers co-vary with phosphorus in the modern ocean because they are nearly quantitatively regenerated from sinking biological debris in the upper water column. Hence they can be used to reconstruct the distribution of labile nutrients in glacial water masses. Independent constraints on glacial deep-ocean circulation patterns could be provided by a tracer of the distribution of silica and alkalinity, the deeply regenerated constituents of planktonic hard parts. Barium shares key aspects of its behaviour with these refractory nutrients because it is removed from solution in surface waters and incorporated into sinking particles which slowly dissolve deep in the water column and in the sediments. The fractionation of Ba between deep-water masses of the major ocean basins is largely controlled by thermohaline circulation patterns, so Ba conforms to different boundary conditions from Cd and 13C. As Ba substitutes into trigonal carbonates, it is a potential palaeoceano-graphic tracer if the Ba content of foraminifera shells reflects ambient dissolved Ba concentrations. Here we present data from Recent core-top benthic foraminifera which indicate that the Ba content of some recent calcitic benthic foraminifera does co-vary with bottom-water Ba.

Results of calculated surface air temperature anomaly of COSMOS Pliocene experiment 2 in NetCDF format

In this manuscript we describe the experimental procedure employed at the Alfred Wegener Institute in Germany in the preparation of the simulations for the Pliocene Model Intercomparison Project (PlioMIP). We present a description of the utilized Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) and document the procedures that we applied to transfer the Pliocene Research, Interpretation and Synoptic Mapping (PRISM) Project mid-Pliocene reconstruction into model forcing fields. The model setup and spin-up procedure are described for both the paleo- and preindustrial (PI) time slices of PlioMIP experiments 1 and 2, and general results that depict the performance of our model setup for mid-Pliocene conditions are presented. The mid-Pliocene, as simulated with our COSMOS setup and PRISM boundary conditions, is both warmer and wetter in the global mean than the PI. The globally averaged annual mean surface air temperature in the mid-Pliocene standalone atmosphere (fully coupled atmosphere-ocean) simulation is 17.35 °C (17.82 °C), which implies a warming of 2.23 °C (3.40 °C) relative to the respective PI control simulation.

Four lectures on mathematics, delivered at Columbia University in 1911,

I. The defintion of solutions of linear partial differnetial equations by boundary conditions.--II. Contemporary researches in differential equations, integral equations, and integro-differential equations.--III. Analysis situs in connection with correspondences and differential equations.--IV. Elementary solutions of partial differential equations and Green's functions.

Spatial and temporal characteristics of summer precipitation over Central Europe in a suite of high-resolution climate models

High-impact, localized intense rainfall episodes represent a major socio-economic problem for societies worldwide, and at the same time these events are notoriously difficult to simulate properly in climate models. Here, the authors investigate how horizontal resolution and model formulation influence this issue by applying the HARMONIE regional climate model (HCLIM) with three different setups; two using convection parameterization at 15 and 6.25 km horizontal resolution (the latter within the “grey-zone” scale), with lateral boundary conditions provided by ERA-Interim reanalysis and integrated over a pan-European domain, and one with explicit convection at 2 km resolution (HCLIM2) over the Alpine region driven by the 15 km model. Seven summer seasons were sampled and validated against two high-resolution observational data sets. All HCLIM versions underestimate the number of dry days and hours by 20-40%, and overestimate precipitation over the Alpine ridge. Also, only modest added value were found of “grey-zone” resolution. However, the single most important outcome is the substantial added value in HCLIM2 compared to the coarser model versions at sub-daily time scales. It better captures the local-to-regional spatial patterns of precipitation reflecting a more realistic representation of the local and meso-scale dynamics. Further, the duration and spatial frequency of precipitation events, as well as extremes, are closer to observations. These characteristics are key ingredients in heavy rainfall events and associated flash floods, and the outstanding results using HCLIM in convection-permitting setting are convincing and encourage further use of the model to study changes in such events in changing climates.

A Kinetic Vlasov Model for Plasma Simulation Using Discontinuous Galerkin Method on Many-Core Architectures

Thesis (Ph.D.)--University of Washington, 2016-03

Is Mindfulness Beneficial for Job Performance Following Positive and Negative Affective Events? The Attenuating Role of Mindfulness on Self-Serving Bias

Thesis (Ph.D.)--University of Washington, 2016-05