Anteroposterior pelvic radiographs to assess acetabular retroversion: high validity of the "cross-over-sign"

Acetabular retroversion has been proposed to contribute to the development of osteoarthritis of the hip. For the diagnosis of this condition, conventional AP pelvic radiographs may represent a reliable, easily available diagnostic modality as they can be obtained with a reproducible technique allowing the anterior and posterior acetabular rims to be visible for assessment. This study was designed to: (i) determine cranial, central, and caudal anatomic acetabular version (AV) from cadaveric specimens; (ii) establish the validity and reliability of the radiographic measurements of central acetabular anteversion; and (iii) determine the validity and reliability of the radiographic "cross-over-sign" to detect acetabular retroversion. Using 43 desiccated pelvises (86 acetabuli) the anatomic AVs were measured at three different transverse planes (cranially, centrally, and caudally). From these pelvises, standardized AP pelvic radiographs were obtained. To directly measure central AV, a modified radiographic method is introduced for the use of AP pelvic radiographs. The validity and reliability of this radiographic method and of the radiographic cross-over-sign to detect cranial acetabular retroversion were determined. The mean central and caudal anatomic AVs were approximately 20 degrees , and the mean cranial AV was 8 degrees . Cranial retroversion (AV < 0 degrees ) was present in 19 of 86 hips (22%). A linear correlation was found between the central and cranial AV. Below 10 degrees of central AV, all acetabuli were cranially retroverted. Between 10 degrees and 20 degrees , 30% of the acetabuli were cranially retroverted, and above 20 degrees , only 1 of 45 acetabuli was cranially retroverted. The radiographic measurement of the central AV (20.3 +/- 6.5 degrees ) correlated strongly with the anatomic AV (20.1 +/- 6.4 degrees ). The sensitivity of the cross-over-sign to detect a cranial acetabular anteversion of less than 4 degrees was 96%, its specificity 95%, and the positive predictive and negative predictive values 90% and 98%, respectively. Both the modified radiographic anteversion measurements and the cross-over-sign demonstrated substantial inter- and intraobserver reliability. Retroversion is almost exclusively a problem of the cranial acetabulum. The cranial AV is on average 12 degrees lower than the central AV, with the latter directly measurable from AP pelvic radiographs. A central AV of less than 10 degrees was associated with cranial retroversion. The presence of a positive cross-over-sign is a highly reliable indicator of cranial AV of <4 degrees.

Diagnostic workup of patients with acute transverse myelitis: spectrum of clinical presentation, neuroimaging and laboratory findings

STUDY DESIGN: Retrospective 9-year survey. OBJECTIVES: Clinical presentation of acute myelitis syndromes is variable, and neuroimaging and laboratory findings are not specific enough to establish the diagnosis with certainty. We evaluated the spectrum clinical features and paraclinical findings encountered during diagnostic workup and aiding the diagnosis. SETTING: Department of Neurology, Inselspital Bern, Switzerland. MATERIAL: Charts and magnetic resonance imaging (MRI) of 63 patients discharged with the diagnosis of acute transverse myelitis. RESULTS: The diagnosis was supported by abnormal MRI and cerebrospinal fluid (CSF) findings in 52 patients (82.5%) and suspected in the remaining either because of a spinal cord MRI lesion suggestive of myelitis (n=5), or abnormal CSF findings (n=4), or electrophysiological evidence of a spinal cord dysfunction (n=2). Clinical impairment was mild (ASIA D) in the majority. All patients had sensory disturbances, whereas motor deficit and autonomic dysfunction were less frequent. Neurological levels were mainly located in cervical or thoracic dermatomes. Spinal cord lesions were visualized by MRI in 90.4% of the patients and distributed either in the cervical or thoracic cord, or both. Multiple lesions were present in more than half of the patients, and lateral, centromedullary and posterior locations were most common. A high percentage of multiple sclerosis (MS)-typical brain lesions and CSF findings suggested a substantial number of MS-related myelitis in our cohort. CONCLUSION: The diagnostic workup of acute myelitis discloses a broad spectrum of CSF or MRI findings, and may be associated with diagnostic uncertainty due to lack of specific CSF or MRI features, or pathological findings.

High cycle fatigue of AA6082 and AA6063 aluminum extrusions

The high cycle fatigue behavior of hollow extruded AA6082 and AA6063 aluminum extrusions has been studied. Hollow extruded aluminum profiles can be processed into intricate shapes, and may be suitable replacements for fatigue critical automotive applications requiring reduced weight. There are several features inherent in hollow aluminum extrusions, such as seam welds, charge welds, microstructural variations and die lines. The effects of such extrusion variables on high cycle fatigue properties were studied by taking specimens from an actual car bumper extrusion. It appears that extrusion die lines create large anisotropy differences in fatigue properties, while welds themselves have little effect on fatigue lives. Removal of die lines greatly increased fatigue properties of AA6082 specimens taken transverse to the extrusion direction. Without die lines, anisotropy in fatigue properties between AA6082 specimens taken longitudinal and transverse to the extrusion direction, was significantly reduced, and properties associated with the orientation of the microstructure appears to be isotropic. A fibrous microstructure for AA6082 specimens showed great improvements in fatigue behavior. The effects of elevated temperatures and exposure of specimens to NaCl solutions was also studied. Exposure to the salt solution greatly reduced the fatigue lives of specimens, while elevated temperatures showed more moderate reductions in fatigue lives.

Spray characterization under flash boiling conditions of a DI multi-hole injector at high injection pressures in a combustion vessel

Spray characterization under flash boiling conditions was investigated utilizing a symmetric multi-hole injector applicable to the gasoline direct injection (GDI) engine. Tests were performed in a constant volume combustion vessel using a high-speed schlieren and Mie scattering imaging systems. Four fuels including n-heptane, 100% ethanol, pure ethanol blended with 15% iso-octane by volume, and test grade E85 were considered in the study. Experimental conditions included various ambient pressure, fuel temperature, and fuel injection pressure. Visualization of the vaporizing spray development was acquired by utilizing schlieren and laser-based Mie scattering techniques. Time evolved spray tip penetration, spray angle, and the ratio of the vapor to liquid region were analyzed by utilizing digital image processing techniques in MATLAB. This research outlines spray characteristics at flash boiling and non-flash boiling conditions. At flash boiling conditions it was observed that individual plumes merge together, leading to significant contraction in spray angle as compared to non-flash boiling conditions. The results indicate that at flash boiling conditions, spray formation and expansion of vapor region is dependent on momentum exchange offered by the ambient gas. A relation between momentum exchange and liquid spray angle formed was also observed.

NLO thermal dilepton rate at non-zero momentum

The vector channel spectral function and the dilepton production rate from a QCD plasma at a temperature above a few hundred MeV are evaluated up to next-to-leading order (NLO) including their dependence on a non-zero momentum with respect to the heat bath. The invariant mass of the virtual photon is taken to be in the range K2 ~ (πT)2 ~ (1GeV)2, generalizing previous NLO results valid for K2 ≫ (πT)2. In the opposite regime 0 < K2 ≪ (πT)2 the loop expansion breaks down, but agrees nevertheless in order of magnitude with a previous result obtained through resummations. Ways to test the vector spectral function through comparisons with imaginary-time correlators measured on the lattice are discussed.

Thermal 2-loop master spectral function at finite momentum

When considering NLO corrections to thermal particle production in the “relativistic” regime, in which the invariant mass squared of the produced particle is K2 ~ (πT)2, then the production rate can be expressed as a sum of a few universal “master” spectral functions. Taking the most complicated 2-loop master as an example, a general strategy for obtaining a convergent 2-dimensional integral representation is suggested. The analysis applies both to bosonic and fermionic statistics, and shows that for this master the non-relativistic approximation is only accurate for K2 ~(8πT)2, whereas the zero-momentum approximation works surprisingly well. Once the simpler masters have been similarly resolved, NLO results for quantities such as the right-handed neutrino production rate from a Standard Model plasma or the dilepton production rate from a QCD plasma can be assembled for K2 ~ (πT)2.

Assessment of Transverse Isotropy in Clinical-Level CT Images of Trabecular Bone Using the Gradient Structure Tensor

The aim of this study was to develop a GST-based methodology for accurately measuring the degree of transverse isotropy in trabecular bone. Using femoral sub-regions scanned in high-resolution peripheral QCT (HR-pQCT) and clinical-level-resolution QCT, trabecular orientation was evaluated using the mean intercept length (MIL) and the gradient structure tensor (GST) on the HR-pQCT and QCT data, respectively. The influence of local degree of transverse isotropy (DTI) and bone mineral density (BMD) was incorporated into the investigation. In addition, a power based model was derived, rendering a 1:1 relationship between GST and MIL eigenvalues. A specific DTI threshold (DTI thres) was found for each investigated size of region of interest (ROI), above which the estimate of major trabecular direction of the GST deviated no more than 30° from the gold standard MIL in 95% of the remaining ROIs (mean error: 16°). An inverse relationship between ROI size and DTI thres was found for discrete ranges of BMD. A novel methodology has been developed, where transversal isotropic measures of trabecular bone can be obtained from clinical QCT images for a given ROI size, DTI thres and power coefficient. Including DTI may improve future clinical QCT finite-element predictions of bone strength and diagnoses of bone disease.

Factorization and resummation for transverse thrust

We analyze transverse thrust in the framework of Soft Collinear Effective Theory and obtain a factorized expression for the cross section that permits resummation of terms enhanced in the dijet limit to arbitrary accuracy. The factorization theorem for this hadron-collider event-shape variable involves collinear emissions at different virtualities and suffers from a collinear anomaly. We compute all its ingredients at the one-loop order, and show that the two-loop input for next-to-next-to-leading logarithmic accuracy can be extracted numerically, from existing fixed-order codes.

Detection efficiency of microchannel plates for e(-) and pi(-) in the momentum range from 17.5 to 345 MeV/c

High-energy e(-) and pi(-) were measured by the multichannel plate (MCP) detector at the PiM1 beam line of the High Intensity Proton Accelerator Facilities located at the Paul Scherrer Institute, Villigen, Switzerland. The measurements provide the absolute detection efficiencies for these particles: 5.8% +/- 0.5% for electrons in the beam momenta range 17.5-300 MeV/c and 6.0% +/- 1.3% for pions in the beam momenta range 172-345 MeV/c. The pulse height distribution determined from the measurements is close to an exponential function with negative exponent, indicating that the particles penetrated the MCP material before producing the signal somewhere inside the channel. Low charge extraction and nominal gains of the MCP detector observed in this study are consistent with the proposed mechanism of the signal formation by penetrating radiation. A very similar MCP ion detector will be used in the Neutral Ion Mass (NIM) spectrometer designed for the JUICE mission of European Space Agency (ESA) to the Jupiter system, to perform measurements of the chemical composition of the Galilean moon exospheres. The detection efficiency for penetrating radiation determined in the present studies is important for the optimisation of the radiation shielding of the NIM detector against the high-rate and high-energy electrons trapped in Jupiter's magnetic field. Furthermore, the current studies indicate that MCP detectors can be useful to measure high-energy particle beams at high temporal resolution. (C) 2015 AIP Publishing LLC.

Hysteresis phenomena in transverse galloping of triangular cross-section bodies

Transverse galloping is a type of aeroelastic instability characterised by large amplitude, low frequency oscillation of a structure in the direction normal to the mean wind direction. It normally appears in bodies with small stiffness and structural damping, provided the incident flow velocity is high enough. In the simplest approach transverse galloping can be considered as a one-degree-of-freedom oscillator subjected to aerodynamic forces, which in turn can be described by using a quasi-steady description. In this frame it has been demonstrated that hysteresis phenomena in transverse galloping is related to the existence of inflection points in the curve giving the dependence with the angle of attack of the aerodynamic coefficient normal to the incident flow. Aiming at experimentally checking such a relationship between these inflection points and hysteresis, wind tunnel experiments have been conducted. Experiments have been restricted to isosceles triangular cross-section bodies, whose galloping behaviour is well documented. Experimental results show that, according to theoretical predictions, hysteresis takes place at the angles of attack where there are inflection points in the lift coefficient curve, provided that the body is prone to gallop at these angles of attack.

High reynolds numbers turbulent boundary layers

Esta tesis estudia las similitudes y diferencias entre los flujos turbulentos de pared de tipo externo e interno, en régimen incompresible, y a números de Reynolds moderada¬mente altos. Para ello consideramos tanto simulaciones numéricas como experimentos de capas límites con gradiente de presiones nulo y de flujos de canal, ambos a números de Reynolds en el rango δ+ ~ 500 - 2000. Estos flujos de cortadura son objeto de numerosas investigaciones debido a la gran importancia que tienen tanto a nivel tecnológico como a nivel de física fundamental. No obstante, todavía existen muchos interrogantes sobre aspectos básicos tales como la universalidad de los perfiles medios y de fluctuación de las velocidades o de la presión, tanto en la zona cercana a la pared como en la zona logarítmica, el escalado y el efecto del número de Reynolds, o las diferencias entre los flujos internos y externos en la zona exterior. En éste estudio hemos utilizado simulaciones numéricas ya existentes de canales y capas límites a números de Reynolds δ+ ~ 2000 y δ+ ~ 700, respectivamente. Para poder comparar ambos flujos a igual número de Reynolds hemos realizado una nueva simulación directa de capa límite en el rango δ+ ~ 1000-2000. Los resultados de la misma son presentados y analizados en detalle. Los datos sin postprocesar y las estadísticas ya postprocesadas están públicamente disponibles en nuestro sitio web.162 El análisis de las estadísticas usando un único punto confirma la existencia de perfiles logarítmicos para las fluctuaciones de la velocidad transversal w'2+ y de la presión p'2+ en ambos tipos de flujos, pero no para la velocidad normal v'2+ o la velocidad longitudinal u'2+. Para aceptar o rechazar la existencia de un rango logarítmico en u'2+ se requieren números de Reynolds más altos que los considerados en éste trabajo. Una de las conse¬cuencias más importantes de poseer tales perfiles es que el valor máximo de la intensidad, que se alcanza cerca de la pared, depende explícitamente del número de Reynolds. Esto ha sido confirmado tras analizar un gran número de datos experimentales y numéricos, cor¬roborando que el máximo de u'2+, p/2+, y w'2+ aumenta proporcionalmente con el log(δ+). Por otro lado, éste máximo es más intenso en los flujos externos que en los internos. La máxima diferencia ocurre en torno a y/δ ~ 0.3-0.5, siendo esta altura prácticamente independiente del número de Reynolds considerado. Estas diferencias se originan como consecuencia del carácter intermitente de las capas límites, que es inexistente en los flujos internos. La estructura de las fluctuaciones de velocidad y de presión, junto con la de los esfuer¬zos de Reynolds, se han investigado por medio de correlaciones espaciales tridimensionales considerando dos puntos de medida. Hemos obtenido que el tamaño de las mismas es gen¬eralmente mayor en canales que en capas límites, especialmente en el caso de la correlación longitudinal Cuu en la dirección del flujo. Para esta correlación se demuestra que las es¬tructuras débilmente correladas presentan longitudes de hasta 0(75), en el caso de capas límites, y de hasta 0(185) en el caso de canales. Estas longitudes se obtienen respecti-vamente en la zona logarítmica y en la zona exterior. Las longitudes correspondientes en la dirección transversal son significativamente menores en ambos flujos, 0(5 — 25). La organización espacial de las correlaciones es compatible con la de una pareja de rollos casi paralelos con dimensiones que escalan en unidades exteriores. Esta organización se mantiene al menos hasta y ~ 0.65, altura a la cual las capas límites comienzan a organi¬zarse en rollos transversales. Este comportamiento es sin embargo más débil en canales, pudiéndose observar parcialmente a partir de y ~ 0.85. Para estudiar si estas estructuras están onduladas a lo largo de la dirección transver¬sal, hemos calculado las correlaciones condicionadas a eventos intensos de la velocidad transversal w'. Estas correlaciones revelan que la ondulación de la velocidad longitudinal aumenta conforme nos alejamos de la pared, sugiriendo que las estructuras están más alineadas en la zona cercana a la pared que en la zona lejana a ella. El por qué de esta ondulación se encuentra posiblemente en la configuración a lo largo de diagonales que presenta w'. Estas estructuras no sólo están onduladas, sino que también están inclinadas respecto a la pared con ángulos que dependen de la variable considerada, de la altura, y de el contorno de correlación seleccionado. Por encima de la zona tampón e independien¬temente del número de Reynolds y tipo de flujo, Cuu presenta una inclinación máxima de unos 10°, las correlaciones Cvv y Cm son esencialmente verticales, y Cww está inclinada a unos 35°. Summary This thesis studies the similitudes and differences between external and internal in¬compressible wall-bounded turbulent flows at moderately-high Reynolds numbers. We consider numerical and experimental zero-pressure-gradient boundary layers and chan¬nels in the range of δ+ ~ 500 — 2000. These shear flows are subjects of intensive research because of their technological importance and fundamental physical interest. However, there are still open questions regarding basic aspects such as the universality of the mean and fluctuating velocity and pressure profiles at the near-wall and logarithmic regions, their scaling and the effect of the Reynolds numbers, or the differences between internal and external flows at the outer layer, to name but a few. For this study, we made use of available direct numerical simulations of channel and boundary layers reaching δ+ ~ 2000 and δ+ ~ 700, respectively. To fill the gap in the Reynolds number, a new boundary layer simulation in the range δ+ ~ 1000-2000 is presented and discussed. The original raw data and the post-processed statistics are publicly available on our website.162 The analysis of the one-point statistic confirms the existence of logarithmic profiles for the spanwise w'2+ and pressure p'2+ fluctuations for both type of flows, but not for the wall-normal v'2+ or the streamwise u'2+ velocities. To accept or reject the existence of a logarithmic range in u'2+ requires higher Reynolds numbers than the ones considered in this work. An important consequence of having such profiles is that the maximum value of the intensities, reached near the wall, depends on the Reynolds number. This was confirmed after surveying a wide number of experimental and numerical datasets, corrob¬orating that the maximum of ul2+, p'2+, and w'2+ increases proportionally to log(δ+). On the other hand, that maximum is more intense in external flows than in internal ones, differing the most around y/δ ~ 0.3-0.5, and essentially independent of the Reynolds number. We discuss that those differences are originated as a consequence of the inter¬mittent character of boundary layers that is absent in internal flows. The structure of the velocity and pressure fluctuations, together with those of the Reynolds shear stress, were investigated using three-dimensional two-point spatial correlations. We find that the correlations extend over longer distances in channels than in boundary layers, especially in the case of the streamwise correlation Cuu in the flow direc-tion. For weakly correlated structures, the maximum streamwise length of Cuu is O(78) for boundary layers and O(188) for channels, attained at the logarithmic and outer regions respectively. The corresponding lengths for the transverse velocities and for the pressure are shorter, 0(8 — 28), and of the same order for both flows. The spatial organization of the velocity correlations is shown to be consistent with a pair of quasi-streamwise rollers that scales in outer units. That organization is observed until y ~ 0.68, from which boundary layers start to organize into spanwise rollers. This effect is weaker in channels, and it appears at y ~ 0.88. We present correlations conditioned to intense events of the transversal velocity, w', to study if these structures meander along the spanwise direction. The results indicate that the streamwise velocity streaks increase their meandering proportionally to the distance to the wall, suggesting that the structures are more aligned close to the wall than far from it. The reason behind this meandering is probably due to the characteristic organization along diagonals of w'. These structures not only meander along the spanwise direction, but they are also inclined to the wall at angles that depend on the distance from the wall, on the variable being considered, and on the correlation level used to define them. Above the buffer layer and independent of the Reynolds numbers and type of flow, the maximum inclination of Cuu is about 10°, Cvv and Cpp are roughly vertical, and Cww is inclined by 35°.

Aeroacoustics in high speed trains

Este trabajo se centra en el estudio de problemas aeroacústicos en los trenes de alta velocidad. Se han considerado dos escenarios en los que las ondas de presión generadas son críticos para el confort de los pasajeros. Uno es el debido a las ondas de presión que genera el tren cuando entra y sale de un túnel, que a su vez producen saltos de presión de baja frecuencia en el tren (cuando se cruzan con él) y en los alrededores del túnel cuando alcanzan la salida. Se estudia este fenómeno, y se propone un sistema aeroelástico basado en el galope transversal para disminuir la energía de estas ondas, y se analiza la energía extraíble de las ondas utilizando cuerpos con diferentes secciones transversales [Sorribes-Palmer and Sanz-Andres, 2013]. La influencia de la geometría de los portales en la energía radiada hacia el exterior de túnel es analizada experimentalmente, prestando especial atención a las boquillas porosas. Las ondas de presión en el interior del túnel se han analizado mediante el método de las características. Se han realizado ensayos experimentales para estimar la energía reflejada hacia el interior del túnel al alcanzar las ondas de presión el portal de salida del túnel. Se ha estudiado la formación e interacción entre el portal del túnel y la onda de choque generada en los túneles de gran longitud y pequeña fricción. Se propone un método para describir de forma aproximada el ruido radiado al exterior. Por otro lado se ha estudiado el ruido de media y alta frecuencia de origen aerodinámico. Se ha estudiado la influencia del desprendimiento de la capa límite sobre el tren. Se propone una metodología basada en una sección de tren característica para predecir rápidamente el nivel de presión de sonido dentro y fuera del tren para todo el rango de frecuencias. Se han realizado medidas experimentales en vía de los espectros de presión sobre la superficie del tren, y de la transmisibilidad de las uniones entre estructura y revestimiento. Los resultados experimentales se han utilizado en los modelos vibroacústicos. El método de la sección del tren característica es especialmente útil a altas frecuencias cuando todo el tren se puede modelar mediante el ensamblaje de diferentes secciones características utilizando el análisis estadístico de la energía. ABSTRACT This work is focused on the study of aeroacoustic problems in high speed trains. We have considered two scenarios in which the pressure waves generated are critical for passengers comfort. The first one is due to the pressure waves generated by a train entering in a tunnel. These waves generate pressure gauges inside the train (when they find each other) and outside of the tunnel portals. This phenomenon has been studied, and an aeroelastic system based on transverse galloping to reduce the energy of these waves is proposed. The maximum extractable energy by using bodies with different cross-section shapes is analyzed. The influence of the portals geometry in the energy radiated outwards the tunnel is analyzed experimentally, with particular attention to the porous exits. The pressure waves inside the tunnel have been analyzed using the method of characteristics. Experimental tests to estimate the energy reflected into the tunnel when the pressure waves reach the tunnel portal have been performed. We have studied the generation and interaction between the tunnel portal and a shock wave generated in long tunnels with small friction. A method to describe in an approximated way the pressure radiated outside the tunnel is proposed. In the second scenario, middle and high frequency noise generated aerodynamically has been studied, including the influence of the detachment of the boundary layer around the train. A method based on a train section to quickly predict the sound pressure level inside and outside the train has been proposed. Experimental test have been performed on board to evaluate the pressure power spectra on the surface of the train, and the transmissibility of the junctions between the structure and trim. These experimental results have been used in the vibroacoustic models. The low frequency pressure waves generated with the train during the tunnel crossing has been identified in the pressure spectrum. The train characteristic section method is especially useful at high frequencies, when the whole train can be modeled by assembling different sections using the statistical energy analysis. The sound pressure level inside the train is evaluated inside and outside the tunnel.

Ultrathin high-index metasurfaces for shaping focused beams

The volume size of a converging wave, which plays a relevant role in image resolution, is governed by the wavelength of the radiation and the numerical aperture (NA) of the wavefront. We designed an ultrathin (λ/8 width) curved metasurface that is able to transform a focused field into a high-NA optical architecture, thus boosting the transverse and (mainly) on-axis resolution. The elements of the metasurface are metal-insulator subwavelength gratings exhibiting extreme anisotropy with ultrahigh index of refraction for TM polarization. Our results can be applied to nanolithography and optical microscopy.

Reform momentum and its impact on Greek growth. Bruegel Policy Contribution 2015/12, July 2015

Greece has an imperfect track-record of structural reform implementation. However, the poor growth outcome of the Greek programmes is also a consequence of the timing and composition of reforms, which were not optimally geared towards a speedy transition to a new growth model based on the private sector. While the main responsibility for this lies with the Greek authorities, international institutions share the responsibility for the poor growth-enhancing effect of reforms. In the current context, further structural reform efforts should be mainly targeted at supporting Greece's speedy return to solid growth rates. This is not only because poverty and unemployment have reached very high levels, but also for political economy reasons: reforms must quickly be seen to be working in order to buttress the consensus in favour of reform. Further efforts should be made to improve Greece’s business environment and to liberalise product markets, in addition to shifting taxation away from labour and towards consumption. Reforms to improve the quality of institutions should continue and are very much needed in the Greek setting, while taking into account that their demanding implementation might use up administrative capacity and their impact on growth will only be seen over long time horizons.

Quantifying ice-shelf basal melt and ice-stream dynamics using high-resolution DEM and GPS time series

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