Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Ancient Greek harbours used as geo-archives for palaeotsunami research : case studies from Krane (Cefalonia), Lechaion (Gulf of Corinth) and Kyllini (Peloponnese)

Since historical times, coastal areas throughout the eastern Mediterranean are exposed to tsunami hazard. For many decades the knowledge about palaeotsunamis was solely based on historical accounts. However, results from timeline analyses reveal different characteristics affecting the quality of the dataset (i.e. distribution of data, temporal thinning backward of events, local periodization phenomena) that emphasize the fragmentary character of the historical data. As an increasing number of geo-scientific studies give convincing examples of well dated tsunami signatures not reported in catalogues, the non-existing record is a major problem to palaeotsunami research. While the compilation of historical data allows a first approach in the identification of areas vulnerable to tsunamis, it must not be regarded as reliable for hazard assessment. Considering the increasing economic significance of coastal regions (e.g. for mass tourism) and the constantly growing coastal population, our knowledge on the local, regional and supraregional tsunami hazard along Mediterranean coasts has to be improved. For setting up a reliable tsunami risk assessment and developing risk mitigation strategies, it is of major importance (i) to identify areas under risk and (ii) to estimate the intensity and frequency of potential events. This approach is most promising when based on the analysis of palaeotsunami research seeking to detect areas of high palaeotsunami hazard, to calculate recurrence intervals and to document palaeotsunami destructiveness in terms of wave run-up, inundation and long-term coastal change. Within the past few years, geo-scientific studies on palaeotsunami events provided convincing evidence that throughout the Mediterranean ancient harbours were subject to strong tsunami-related disturbance or destruction. Constructed to protect ships from storm and wave activity, harbours provide especially sheltered and quiescent environments and thus turned out to be valuable geo-archives for tsunamigenic high-energy impacts on coastal areas. Directly exposed to the Hellenic Trench and extensive local fault systems, coastal areas in the Ionian Sea and the Gulf of Corinth hold a considerably high risk for tsunami events, respectively.Geo-scientific and geoarcheaological studies carried out in the environs of the ancient harbours of Krane (Cefalonia Island), Lechaion (Corinth, Gulf of Corinth) and Kyllini (western Peloponnese) comprised on-shore and near-shore vibracoring and subsequent sedimentological, geochemical and microfossil analyses of the recovered sediments. Geophysical methods like electrical resistivity tomography and ground penetrating radar were applied in order to detect subsurface structures and to verify stratigraphical patterns derived from vibracores over long distances. The overall geochronological framework of each study area is based on radiocarbon dating of biogenic material and age determination of diagnostic ceramic fragments. Results presented within this study provide distinct evidence of multiple palaeotsunami landfalls for the investigated areas. Tsunami signatures encountered in the environs of Krane, Lechaion and Kyllini include (i) coarse-grained allochthonous marine sediments intersecting silt-dominated quiescent harbour deposits and/or shallow marine environments, (ii) disturbed microfaunal assemblages and/or (iii) distinct geochemical fingerprints as well as (iv) geo-archaeological destruction layers and (v) extensive units of beachrock-type calcarenitic tsunamites. For Krane, geochronological data yielded termini ad or post quem (maximum ages) for tsunami event generations dated to 4150 ± 60 cal BC, ~ 3200 ± 110 cal BC, ~ 650 ± 110 cal BC, and ~ 930 ± 40 cal AD, respectively. Results for Lechaion suggest that the harbour was hit by strong tsunami impacts in the 8th-6th century BC, the 1st-2nd century AD and in the 6th century AD. At Kyllini, the harbour site was affected by tsunami impact in between the late 7th and early 4th cent. BC and between the 4th and 6th cent. AD. In case of Lechaion and Kyllini, the final destruction of the harbour facilities also seems to be related to the tsunami impact. Comparing the tsunami signals obtained for each study areas with geo-scientific data from palaeotsunami events from other sites indicates that the investigated harbour sites represent excellent geo-archives for supra-regional mega-tsunamis.

A new information fusion method for land-use classification using high resolution satellite imagery

Satellite image classification involves designing and developing efficient image classifiers. With satellite image data and image analysis methods multiplying rapidly, selecting the right mix of data sources and data analysis approaches has become critical to the generation of quality land-use maps. In this study, a new postprocessing information fusion algorithm for the extraction and representation of land-use information based on high-resolution satellite imagery is presented. This approach can produce land-use maps with sharp interregional boundaries and homogeneous regions. The proposed approach is conducted in five steps. First, a GIS layer - ATKIS data - was used to generate two coarse homogeneous regions, i.e. urban and rural areas. Second, a thematic (class) map was generated by use of a hybrid spectral classifier combining Gaussian Maximum Likelihood algorithm (GML) and ISODATA classifier. Third, a probabilistic relaxation algorithm was performed on the thematic map, resulting in a smoothed thematic map. Fourth, edge detection and edge thinning techniques were used to generate a contour map with pixel-width interclass boundaries. Fifth, the contour map was superimposed on the thematic map by use of a region-growing algorithm with the contour map and the smoothed thematic map as two constraints. For the operation of the proposed method, a software package is developed using programming language C. This software package comprises the GML algorithm, a probabilistic relaxation algorithm, TBL edge detector, an edge thresholding algorithm, a fast parallel thinning algorithm, and a region-growing information fusion algorithm. The county of Landau of the State Rheinland-Pfalz, Germany was selected as a test site. The high-resolution IRS-1C imagery was used as the principal input data.

Das Paläolithikum der Magdalenahöhle bei Gerolstein

Die vorliegende Arbeit hat das von 1969 – 1972 ergrabene Inventar der Magdalenahöhle bei Gerolstein unter kritischer Berücksichtigung der Originaldokumentation sowie der stratigraphischen und sedimentologischen Beschreibungen erneut untersucht und v.a. hinsichtlich zweier Arbeitshypothesen überprüft. Daneben fanden jedoch auch die Schmuckobjekte und in kursorischer Weise die Tierknochen Betrachtung. Die Elfenbeinobjekte setzen sich aus elf Fragmenten zusammen, die bereits in zerbrochenem Zustand in die Höhle gelangt sein müssen. Sie sind mit mehreren Linienbündeln verziert, die teilweise aus v- oder zickzack-förmigen Motiven bestehen. Auch gestanzte Punktreihen treten auf. In ihrer Größe und Form sind die Elfenbeinobjekte einzigartig. Lediglich aus der Csákvár-Höhle in Ungarn gibt es vergleichbare Stücke, deren genaue Altersstellung jedoch unklar ist. Daneben kommen in der Magdalenahöhle zwei vollständige durchlochte Hirschgrandeln sowie die Fragmente einer durchlochten Grandel sowie eines durchlochten Wolfzahns vor. Diese tragen teilweise Spuren einer Aufhängung bzw. Befestigung. Der Grund für ihre Niederlegung vor Ort ist indes nicht endgültig zu klären. Die überlieferten Tierknochen besitzen verschiedene Grade von Verfärbung und Erhaltung, lassen sich dadurch jedoch nicht verschiedenen Schichten zuordnen. Neben Modifikationen von Carnivoren, darunter v.a. durch Verdauungsprozesse, sind auch an einigen Exemplaren Schnittspuren festgestellt worden. Eine Bärenphalange aus der Fundschicht b1 eröffnet so die Perspektive, die menschliche Belegung erneut mit der Radiokohlenstoffmethode direkt zu datieren. Der Untere paläolithische Fundhorizont besteht aus relativ unspezifischen Quarzartefakten, die von einer opportunistischen Abschlags- und Werkzeugsgewinnung aus lokalen Schottern zeugen. Für den Oberen paläolithischen Fundhorizont zeigt die Steinartefaktanalyse, dass die Abschläge als Herstellungsreste dünner bifazieller Geräte angesprochen werden können. Während dieser Befund alleine auch im Zusammenhang mit den spätmittelpaläolithischen Blattspitzengruppen gesehen werden kann, sprechen die bereits erwähnten vergesellschafteten Schmuckobjekte, der Nachweis eines Klingenabbaus sowie die fast ausschließliche Verwendung exogenen Rohmaterials für einen jungpaläolithischen Kontext, d.h. für eine Affinität zum Solutréen. Die Steinartefakte der Magdalenahöhle zeugen gleichzeitig von einer sehr mobilen Lebensweise, da lediglich eine Phase des Herstellungsprozesses des bifaziellen Geräts vor Ort stattgefunden hat. Gleichzeitig wurden früher an anderem Ort gewonnene Abschläge mitgebracht und als Werkzeug verwendet. Ebenso wie der Kern selbst wurden Abschläge auch wieder abtransportiert. Insgesamt kann daher eine Rohmaterial konservierende Strategie rekonstruiert werden, in der neben dem Kern als Gerät selbst auch die Abschläge Verwendung finden. Da die Magdalenahöhle als östlicher Niederschlag des Solutréen und als Beleg für das maximale Verbreitungsgebiet dieses Technokomplexes gewertet werden muss, reiht sie sich in eine Reihe von Fundstellen und Argumenten ein, wonach das zentrale Mitteleuropa während des LGM s.l. nicht menschenleer war, sondern in sporadischen Exkursionen begangen wurde. Obwohl sich daran noch keine dauerhafte Wiederbesiedlung anschloss, muss vom Bild einer absoluten Siedlungsleere Abstand genommen werden. Weitere Fundstellen und absolutchronologische Datierungen, u.a. der Magdalenahöhle, könnten in Zukunft zu einem noch besseren Verständnis der menschlichen Anpassungsstragien an kaltzeitliche Umwelten beitragen.

Biaxial extensional motion of an inertially driven radially expanding liquid sheet

We consider the inertially driven, time-dependent biaxial extensional motion of inviscid and viscous thinning liquid sheets. We present an analytic solution describing the base flow and examine its linear stability to varicose (symmetric) perturbations within the framework of a long-wave model where transient growth and long-time asymptotic stability are considered. The stability of the system is characterized in terms of the perturbation wavenumber, Weber number, and Reynolds number. We find that the isotropic nature of the base flow yields stability results that are identical for axisymmetric and general two-dimensional perturbations. Transient growth of short-wave perturbations at early to moderate times can have significant and lasting influence on the long-time sheet thickness. For finite Reynolds numbers, a radially expanding sheet is weakly unstable with bounded growth of all perturbations, whereas in the inviscid and Stokes flow limits sheets are unstable to perturbations in the short-wave limit.

Retinal thickness in Parkinson's disease

BACKGROUND: Visual symptoms are common in Parkinson's disease with studies consistently demonstrating reductions in visual acuity, contrast sensitivity, colour and motion perception as well as alterations in electroretinogram latencies and amplitudes. Optical coherence tomography can examine retinal structure non-invasively and retinal thinning has been suggested as a potential biomarker for neurodegeneration in Parkinson's disease. Our aim was to examine the retinal thickness of a cohort of Parkinson's disease subjects (and age-matched controls) to establish the practical utility of optical coherence tomography in a representative older Parkinson's disease group. METHODS: Fifty-one established Parkinson's disease subjects and 25 healthy controls were subjected to ophthalmological assessment and optical coherence tomography (Zeiss Stratus 3000™) of macular thickness and volume and retinal nerve fibre thickness around the optic nerve head. Twenty four percent of control and 20% of Parkinson's disease subjects were excluded from final analysis due to co-morbid ocular pathology. Further data was excluded either due to poor tolerability of optical coherence tomography or poor quality scans. RESULTS: Despite a reduction in both visual acuity and contrast sensitivity in the residual evaluable Parkinson's disease cohort, we did not detect any differences between the two study groups for any measures of retinal thickness, in contrast to previously published work. CONCLUSIONS: In addition to technical problems inherent in the evaluation, the lack of difference between Parkinson's disease and healthy control subjects suggests longitudinal studies, employing newer techniques, will be required to define the role of optical coherence tomography as a potential diagnostic biomarker.

Cortical morphometry in narcolepsy with cataplexy

The sleep-wake disorder narcolepsy with cataplexy is associated with the loss of hypocretin-(orexin-) producing neurons in the lateral hypothalamus. Several studies have reported abnormal cerebral activation in patients with narcolepsy with cataplexy. It remains unclear, however, whether these functional changes are related to structural alterations, particularly at the cortical level. To quantify structural brain changes associated with narcolepsy with cataplexy, we used high-resolution T1-weighted magnetic resonance imaging (MRI) in 12 patients compared with 12 healthy participants matched for age and gender. Subcortical and regional cortical volumes were measured using a method unbiased by non-linear registration. Further whole-brain analyses were conducted, measuring cortical characteristics, such as cortical thickness and gyrification, at thousands of points across each hemisphere using validated algorithms. Statistical analyses accounted for an effect of age and gender. We observed decreased cortical volume in the left paracentral lobule and increased cortical volume in the left caudal part of the middle frontal gyrus in narcoleptic patients compared with controls. Cortical thickness in prefrontal areas was inversely correlated with the severity of narcolepsy. Further, we observed several clusters of cortical thinning in patients with childhood or adolescent onset of narcolepsy compared with patients with adult onset of the disease. Our results suggest that specific anatomical changes may differentiate subgroups of narcolepsy patients with different clinical profiles (such as varying symptom severity or different age at onset). Future studies with larger groups of sleepy patients are required to assess whether distinct patterns of anatomical changes may distinguish narcolepsy from non-hypocretin-deficient hypersomnia disorders.

TALDICE-1 age scale of the Talos Dome deep ice core, East Antarctica

A new deep ice core drilling program, TALDICE, has been successfully handled by a European team at Talos Dome, in the Ross Sea sector of East Antarctica, down to 1620 m depth. Using stratigraphic markers and a new inverse method, we produce the first official chronology of the ice core, called TALDICE-1. We show that it notably improves an a priori chronology resulting from a one-dimensional ice flow model. It is in agreement with a posteriori controls of the resulting accumulation rate and thinning function along the core. An absolute uncertainty of only 300 yr is obtained over the course of the last deglaciation. This uncertainty remains lower than 600 yr over Marine Isotope Stage 3, back to 50 kyr BP. The phasing of the TALDICE ice core climate record with respect to the central East Antarctic plateau and Greenland records can thus be determined with a precision allowing for a discussion of the mechanisms at work at sub-millennial time scales.

Development and remodeling of the vertebrate blood-gas barrier

During vertebrate development, the lung inaugurates as an endodermal bud from the primitive foregut. Dichotomous subdivision of the bud results in arborizing airways that form the prospective gas exchanging chambers, where a thin blood-gas barrier (BGB) is established. In the mammalian lung, this proceeds through conversion of type II cells to type I cells, thinning, and elongation of the cells as well as extrusion of the lamellar bodies. Subsequent diminution of interstitial tissue and apposition of capillaries to the alveolar epithelium establish a thin BGB. In the noncompliant avian lung, attenuation proceeds through cell-cutting processes that result in remarkable thinning of the epithelial layer. A host of morphoregulatory molecules, including transcription factors such as Nkx2.1, GATA, HNF-3, and WNT5a; signaling molecules including FGF, BMP-4, Shh, and TFG- β and extracellular proteins and their receptors have been implicated. During normal physiological function, the BGB may be remodeled in response to alterations in transmural pressures in both blood capillaries and airspaces. Such changes are mitigated through rapid expression of the relevant genes for extracellular matrix proteins and growth factors. While an appreciable amount of information regarding molecular control has been documented in the mammalian lung, very little is available on the avian lung.

Duration of Greenland Stadial 22 and ice-gas Δage from counting of annual layers in Greenland NGRIP ice core

High-resolution measurements of chemical impurities and methane concentrations in Greenland ice core samples from the early glacial period allow the extension of annual-layer counted chronologies and the improvement of gas age-ice age difference (Δage) essential to the synchronization of ice core records. We report high-resolution measurements of a 50 m section of the NorthGRIP ice core and corresponding annual layer thicknesses in order to constrain the duration of the Greenland Stadial 22 (GS-22) between Greenland Interstadials (GIs) 21 and 22, for which inconsistent durations and ages have been reported from Greenland and Antarctic ice core records as well as European speleothems. Depending on the chronology used, GS-22 occurred between approximately 89 (end of GI-22) and 83 kyr b2k (onset of GI-21). From annual layer counting, we find that GS-22 lasted between 2696 and 3092 years and was followed by a GI-21 pre-cursor event lasting between 331 and 369 yr. Our layer-based counting agrees with the duration of stadial 22 as determined from the NALPS speleothem record (3250 ± 526 yr) but not with that of the GICC05modelext chronology (2620 yr) or an alternative chronology based on gas-marker synchronization to EPICA Dronning Maud Land ice core. These results show that GICC05modelext overestimates accumulation and/or underestimates thinning in this early part of the last glacial period. We also revise the possible ranges of NorthGRIP Δdepth (5.49 to 5.85 m) and Δage (498 to 601 yr) at the warming onset of GI-21 as well as the Δage range at the onset of the GI-21 precursor warming (523 to 654 yr), observing that temperature (represented by the δ15N proxy) increases before CH4 concentration by no more than a few decades.

Melting of Northern Greenland during the last interglaciation

Using simulated climate data from the comprehensive coupled climate model IPSL CM4, we simulate the Greenland ice sheet (GrIS) during the Eemian interglaciation with the three-dimensional ice sheet model SICOPOLIS. The Eemian is a period 126 000 yr before present (126 ka) with Arctic temperatures comparable to projections for the end of this century. In our simulation, the northeastern part of the GrIS is unstable and retreats significantly, despite moderate melt rates. This result is found to be robust to perturbations within a wide parameter space of key parameters of the ice sheet model, the choice of initial ice temperature, and has been reproduced with climate forcing from a second coupled climate model, the CCSM3. It is shown that the northeast GrIS is the most vulnerable. Even a small increase in melt removes many years of ice accumulation, giving a large mass imbalance and triggering the strong ice-elevation feedback. Unlike the south and west, melting in the northeast is not compensated by high accumulation. The analogy with modern warming suggests that in coming decades, positive feedbacks could increase the rate of mass loss of the northeastern GrIS, exceeding the recent observed thinning rates in the south.

Estimating the number of motor units using random sums with independently thinned terms

The problem of estimating the numbers of motor units N in a muscle is embedded in a general stochastic model using the notion of thinning from point process theory. In the paper a new moment type estimator for the numbers of motor units in a muscle is denned, which is derived using random sums with independently thinned terms. Asymptotic normality of the estimator is shown and its practical value is demonstrated with bootstrap and approximative confidence intervals for a data set from a 31-year-old healthy right-handed, female volunteer. Moreover simulation results are presented and Monte-Carlo based quantiles, means, and variances are calculated for N in{300,600,1000}.

Soft tissue profile changes following mandibular advancement and setback surgery an average of 12 years postoperatively

PURPOSE: The aim of this study was to assess long-term changes in position of soft tissue landmarks following mandibular advancement and setback surgery. MATERIALS AND METHODS: Twenty-seven patients (14 women, 13 men; mean age, 36 years) who had undergone either mandibular advancement (15 patients) or setback surgery (12 patients), were available for a long-term follow-up an average of 12 years postoperatively. In all of these cases, lateral cephalometric radiographs taken immediately before operation, at 1 week, 14 months, and 12 years postoperatively, were studied. RESULTS: During the 14 months postoperatively, soft tissue chin and mentolabial fold followed its underlying hard tissue in all patients. A continuous skeletal relapse was observable 12 years after mandibular advancement, but soft tissue chin moved more in an anterior direction. After mandibular setback, soft and hard tissue landmarks remained almost unchanged. Over the entire observation period, a thickening of soft tissue at pogonion was generally seen, and particularly a thickening of the whole chin in the setback group. All patients showed a significant lengthening and thinning of the upper lip. In all except 2 males, the patient's body weight increased markedly. CONCLUSION: In contrast to the immediate postoperative stage, soft tissue changes observed an average of 12 years after the primary operation do not directly follow the movements of the underlying skeletal structure. The soft tissue profile changes observed over such a long term seem to be influenced not only by the underlying skeletal structure but also by other factors such as weight gain and aging process.

Synchronisation of the EDML and EDC ice cores for the last 52 kyr by volcanic signature matching

A common time scale for the EPICA ice cores from Dome C (EDC) and Dronning Maud Land (EDML) has been established. Since the EDML core was not drilled on a dome, the development of the EDML1 time scale for the EPICA ice core drilled in Dronning Maud Land was based on the creation of a detailed stratigraphic link between EDML and EDC, which was dated by a simpler 1D ice-flow model. The synchronisation between the two EPICA ice cores was done through the identification of several common volcanic signatures. This paper describes the rigorous method, using the signature of volcanic sulfate, which was employed for the last 52 kyr of the record. We estimated the discrepancies between the modelled EDC and EDML glaciological age scales during the studied period, by evaluating the ratio R of the apparent duration of temporal intervals between pairs of isochrones. On average R ranges between 0.8 and 1.2 corresponding to an uncertainty of up to 20% in the estimate of the time duration in at least one of the two ice cores. Significant deviations of R up to 1.4–1.5 are observed between 18 and 28 kyr before present (BP), where present is defined as 1950. At this stage our approach does not allow us unequivocally to find out which of the models is affected by errors, but assuming that the thinning function at both sites and accumulation history at Dome C (which was drilled on a dome) are correct, this anomaly can be ascribed to a complex spatial accumulation variability (which may be different in the past compared to the present day) upstream of the EDML core.

Federal and state policy influence on woody biomass utilization

A considerable portion of public lands in the United States is at risk of uncharacteristically severe wildfires due to a history of fire suppression. Wildfires already have detrimental impacts on the landscape and on communities in the wildland-urban interface (WUI) due to unnatural and overstocked forests. Strategies to mitigate wildfire risk include mechanical thinning and prescribed burning in areas with high wildfire risk. The material removed is often of little or no economic value. Woody biomass utilization (WBU) could offset the costs of hazardous fuel treatments if removed material could be used for wood products, heat, or electricity production. However, barriers due to transportation costs, removal costs, and physical constraints (such as steep slopes) hinder woody biomass utilization. Various federal and state policies attempt to overcome these barriers. WBU has the potential to aid in wildfire mitigation and meet growing state mandates for renewable energy. This research utilizes interview data from individuals involved with on-the-ground woody biomass removal and utilization to determine how federal and state policies influence woody biomass utilization. Results suggest that there is not one over-arching policy that hinders or promotes woody biomass utilization, but rather woody biomass utilization is hindered by organizational constraints related to time, cost, and quality of land management agencies’ actions. However, the use of stewardship contracting (a hybrid timber sale and service contract) shows promise for increased WBU, especially in states with favorable tax policies and renewable energy mandates. Policy recommendations to promote WBU include renewal of stewardship contracting legislations and a re-evaluation of land cover types suited for WBU. Potential future policies to consider include the indirect role of carbon dioxide emission reduction activities to promote wood energy and future impacts of air quality regulations.