Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.

The reconstruction of skeletal movement: the soft tissue artefact issue

In 3D human movement analysis performed using stereophotogrammetric systems and skin markers, bone pose can only be estimated in an indirect fashion. During a movement, soft tissue deformations make the markers move with respect to the underlying bone generating soft tissue artefact (STA). STA has devastating effects on bone pose estimation and its compensation remains an open question. The aim of this PhD thesis was to contribute to the solution of this crucial issue. Modelling STA using measurable trial-specific variables is a fundamental prerequisite for its removal from marker trajectories. Two STA model architectures are proposed. Initially, a thigh marker-level artefact model is presented. STA was modelled as a linear combination of joint angles involved in the movement. This model was calibrated using ex-vivo and in-vivo STA invasive measures. The considerable number of model parameters led to defining STA approximations. Three definitions were proposed to represent STA as a series of modes: individual marker displacements, marker-cluster geometrical transformations (MCGT), and skin envelope shape variations. Modes were selected using two criteria: one based on modal energy and another on the selection of modes chosen a priori. The MCGT allows to select either rigid or non-rigid STA components. It was also empirically demonstrated that only the rigid component affects joint kinematics, regardless of the non-rigid amplitude. Therefore, a model of thigh and shank STA rigid component at cluster-level was then defined. An acceptable trade-off between STA compensation effectiveness and number of parameters can be obtained, improving joint kinematics accuracy. The obtained results lead to two main potential applications: the proposed models can generate realistic STAs for simulation purposes to compare different skeletal kinematics estimators; and, more importantly, focusing only on the STA rigid component, the model attains a satisfactory STA reconstruction with less parameters, facilitating its incorporation in an pose estimator.

Sensitive and selective detection of free FXIII activation peptide: a potential marker of acute thrombotic events

Coagulation factor XIII (FXIII) stabilizes fibrin fibers and is therefore a major player in the maintenance of hemostasis. FXIII is activated by thrombin resulting in cleavage and release of the FXIII activation peptide (AP-FXIII). The objective of this study was to characterize the released AP-FXIII and determine specific features that may be used for its specific detection. We analyzed the structure of bound AP-FXIII within the FXIII A-subunit and interactions of AP-FXIII by hydrogen bonds with both FXIII A-subunit monomers. We optimized our previously developed AP-FXIII ELISA by using 2 monoclonal antibodies. We determined high binding affinities between the antibodies and free AP-FXIII and demonstrated specific binding by epitope mapping analyses with surface plasmon resonance and enzyme-linked immunosorbent assay. Because the structure of free AP-FXIII had been characterized so far by molecular modeling only, we performed structural analysis by nuclear magnetic resonance. Recombinant AP-FXIII was largely flexible both in plasma and water, differing significantly from the rigid structure in the bound state. We suggest that the recognized epitope is either occluded in the noncleaved form or possesses a structure that does not allow binding to the antibodies. On the basis of our findings, we propose AP-FXIII as a possible new marker for acute thrombotic events.

Segmentation of brain tumor images based on atlas-registration combined with a Markov-Random-Field lesion growth model

We present an automatic method to segment brain tissues from volumetric MRI brain tumor images. The method is based on non-rigid registration of an average atlas in combination with a biomechanically justified tumor growth model to simulate soft-tissue deformations caused by the tumor mass-effect. The tumor growth model, which is formulated as a mesh-free Markov Random Field energy minimization problem, ensures correspondence between the atlas and the patient image, prior to the registration step. The method is non-parametric, simple and fast compared to other approaches while maintaining similar accuracy. It has been evaluated qualitatively and quantitatively with promising results on eight datasets comprising simulated images and real patient data.

Making Research Relevant to Development - Experience from the NCCR North-South

Making research relevant to development is a complex, non-linear and often unpredictable process which requires very particular skills and strategies on the part of researchers. The National Centre of Competence in Research (NCCR) North-South provides financial and technical support for researchers so that they can effectively cooperate with policy-makers and practitioners. An analysis of 10 years of experience translating research into development practise in the NCCR North-South revealed the following four strategies as particularly relevant: a) research orientation towards the needs and interests of partners; b) implementation of promising methods and approaches; c) communication and dissemination of research results; and d) careful analysis of the political context through monitoring and learning approaches. The NCCR North-South experience shows that “doing excellent research” is just one piece of the mosaic. It is equally important to join hands with non-academic partners from the very beginning of a research project, in order to develop and test new pathways for sustainable development. Capacity building – in the North and South – enables researchers to do both: To do excellent research and to make it relevant for development.

From Craton to Rift: Empirically Based Ground-Truth Criteria for Local Events Recorded on Regional Networks

Region-specific empirically based ground-truth (EBGT) criteria used to estimate the epicentral-location accuracy of seismic events have been developed for the Main Ethiopian Rift and the Tibetan plateau. Explosions recorded during the Ethiopia-Afar Geoscientific Lithospheric Experiment (EAGLE), the International Deep Profiling of Tibet, and the Himalaya (INDEPTH III) experiment provided the necessary GT0 reference events. In each case, the local crustal structure is well known and handpicked arrival times were available, facilitating the establishment of the location accuracy criteria through the stochastic forward modeling of arrival times for epicentral locations. In the vicinity of the Main Ethiopian Rift, a seismic event is required to be recorded on at least 8 stations within the local Pg/Pn crossover distance and to yield a network-quality metric of less than 0.43 in order to be classified as EBGT5(95%) (GT5 with 95% confidence). These criteria were subsequently used to identify 10 new GT5 events with magnitudes greater than 2.1 recorded on the Ethiopian Broadband Seismic Experiment (EBSE) network and 24 events with magnitudes greater than 2.4 recorded on the EAGLE broadband network. The criteria for the Tibetan plateau are similar to the Ethiopia criteria, yet slightly less restrictive as the network-quality metric needs to be less than 0.45. Twenty-seven seismic events with magnitudes greater than 2.5 recorded on the INDEPTH III network were identified as GT5 based on the derived criteria. When considered in conjunction with criteria developed previously for the Kaapvaal craton in southern Africa, it is apparent that increasing restrictions on the network-quality metric mirror increases in the complexity of geologic structure from craton to plateau to rift. Accession Number: WOS:000322569200012

Challenges for Participatory Conservation in Times of Global Change: Lessons from a Comparative Analysis and New Developments

Participatory approaches to conservation have been applied worldwide by governments and non-governmental organisations. However, results from a comparative analysis of the impacts of global change on management issues in 13 protected areas in Africa, Latin America, Asia, and Europe show that in many cases the involvement of local people has remained limited, and economic gains for local livelihoods have been limited or non-existent. Viewed from a ‘new institutionalist’ perspective and focusing on power relations and ideologies, the results of this study carried out within the framework of the Swiss National Centre of Competence in Research (NCCR) North-South show that in African cases local people do not feel part of the process and, therefore, become disengaged. In Asia, and even more so in Latin America, local indigenous peoples and their leaders support protected areas as a means to gain political rights over areas threatened by immigration. The European (Swiss) case is the only one where political rights and economic incentives present a context in which participation is of direct interest to local people. Meanwhile, recent debates on new global conservation developments in the context of climate change policy indicate a growing tendency to treat conservation as a commodity. We argue that this can have problematical effects on efforts to devolve power to the local level in the context of conservation.

Non-State United States Support For The Contras During The Reagan Era

Under President Ronald Reagan, the White House pursued a complex foreign policy towards the Contras, rebels in trying to overthrow the Sandinista regime in Nicaragua, in Nicaragua. In 1979, the leftist Sandinista government seized power in Nicaragua. The loss of the previous pro-United States Somoza military dictatorship deeply troubled the conservatives, for whom eradication of communism internationally was a top foreign policy goal. Consequently, the Reagan Administration sought to redress the policy of his predecessor, Jimmy Carter, and assume a hard line stance against leftist regimes in Central America. Reagan and the conservatives within his administration, therefore, supported the Contra through military arms, humanitarian aid, and financial contributions. This intervention in Nicaragua, however, failed to garner popular support from American citizens and Democrats. Consequently, between 1982 and 1984 Congress prohibited further funding to the Contras in a series of legislation called the Boland Amendments. These Amendments barred any military aid from reaching the Contras, including through intelligence agencies. Shortly after their passage, Central Intelligence Agency Director William Casey and influential members of Reagan¿s National Security Council (NSC) including National Security Advisor Robert McFarlane, NSC Aide Oliver North, and Deputy National Security Advisor John Poindexter cooperated to identify and exploit loopholes in the legislation. By recognizing the NSC as a non-intelligence body, these masterminds orchestrated a scheme in which third parties, including foreign countries and private donors, contributed both financially and through arms donations to sustain the Contras independently of Congressional oversight. This thesis explores the mechanism and process of soliciting donations from private individuals, recognizing the forces and actors that created a situation for covert action to continue without detection. Oliver North, the main actor of the state, worked within his role as an NSC bureaucrat to network with influential politicians and private individuals to execute the orders of his superiors and shape foreign policy. Although Reagan articulated his desire for the Contras to remain a military presence in Nicaragua, he delegated the details of policy to his subordinates, which allowed this scheme to flourish. Second, this thesis explores the individual donors, analyzing their role as private citizens in sustaining and encouraging the policy of the Reagan Administration. The Contra movement found non-state support from followers of the New Right, demonstrated through financial and organizational assistance, that allowed the Reagan Administration¿s statistically unpopular policy in Nicaragua to continue. I interpret these donors as politically involved, but politically philanthropic, individuals, donating to their charity of choice to further the principles of American freedom internationally in a Cold War environment. The thesis then proceeds to assess the balance of power between the executive and other political actors in shaping policy, concluding that the executive cannot act alone in the formulation and implementation of foreign policy.

Assessment Of Dynamic Properties Of A Crowd Model For Human-Structure Interaction Modeling

As lightweight and slender structural elements are more frequently used in the design, large scale structures become more flexible and susceptible to excessive vibrations. To ensure the functionality of the structure, dynamic properties of the occupied structure need to be estimated during the design phase. Traditional analysis method models occupants simply as an additional mass; however, research has shown that human occupants could be better modeled as an additional degree-of- freedom. In the United Kingdom, active and passive crowd models are proposed by the Joint Working Group as a result of a series of analytical and experimental research. It is expected that the crowd models would yield a more accurate estimation to the dynamic response of the occupied structure. However, experimental testing recently conducted through a graduate student project at Bucknell University indicated that the proposed passive crowd model might be inaccurate in representing the impact on the structure from the occupants. The objective of this study is to provide an assessment of the validity of the crowd models proposed by JWG through comparing the dynamic properties obtained from experimental testing data and analytical modeling results. The experimental data used in this study was collected by Firman in 2010. The analytical results were obtained by performing a time-history analysis on a finite element model of the occupied structure. The crowd models were created based on the recommendations from the JWG combined with the physical properties of the occupants during the experimental study. During this study, SAP2000 was used to create the finite element models and to implement the analysis; Matlab and ME¿scope were used to obtain the dynamic properties of the structure through processing the time-history analysis results from SAP2000. The result of this study indicates that the active crowd model could quite accurately represent the impact on the structure from occupants standing with bent knees while the passive crowd model could not properly simulate the dynamic response of the structure when occupants were standing straight or sitting on the structure. Future work related to this study involves improving the passive crowd model and evaluating the crowd models with full-scale structure models and operating data.

Statistical shape analysis via principal factor analysis

Statistical shape analysis techniques commonly employed in the medical imaging community, such as active shape models or active appearance models, rely on principal component analysis (PCA) to decompose shape variability into a reduced set of interpretable components. In this paper we propose principal factor analysis (PFA) as an alternative and complementary tool to PCA providing a decomposition into modes of variation that can be more easily interpretable, while still being a linear efficient technique that performs dimensionality reduction (as opposed to independent component analysis, ICA). The key difference between PFA and PCA is that PFA models covariance between variables, rather than the total variance in the data. The added value of PFA is illustrated on 2D landmark data of corpora callosa outlines. Then, a study of the 3D shape variability of the human left femur is performed. Finally, we report results on vector-valued 3D deformation fields resulting from non-rigid registration of ventricles in MRI of the brain.

Fitting 3D morphable models using implicit representations

We consider the problem of approximating the 3D scan of a real object through an affine combination of examples. Common approaches depend either on the explicit estimation of point-to-point correspondences or on 2-dimensional projections of the target mesh; both present drawbacks. We follow an approach similar to [IF03] by representing the target via an implicit function, whose values at the vertices of the approximation are used to define a robust cost function. The problem is approached in two steps, by approximating first a coarse implicit representation of the whole target, and then finer, local ones; the local approximations are then merged together with a Poisson-based method. We report the results of applying our method on a subset of 3D scans from the Face Recognition Grand Challenge v.1.0.

The Benefits of Multilingualism for Statutory and Constitutional Interpretation in South Africa

Taking the South African experience as an example, this article considers the interpretive benefits to be reaped from having access to bi- and multilingual versions of a statutory text. The discussion takes place against the backdrop of a history of statutory bi- and multilingualism in the said jurisdiction as well as, at present, constitutional guarantees of language rights and the “parity of esteem” of eleven official languages. It is argued that, if invoked with due discretion and in a non-rigid way, statutory multilingualism can be a boon to statutory and constitutional interpretation. The South African courts – whose traditional approach to statutory inter-pretation has tended to be literalist, formalistic and formulaic – are, generally speaking, to be commended for their supple use of bilingualism as an aid to interpretation over the years. The advent of constitutional multilingualism and the (potential) availability of statutory texts (and the Constitution) in more than two languages, have moreover created conditions conducive to the further development and refinement of reliance on multilingualism in statutory and constitutional interpretation – certain challenges notwithstanding.

Kyzylkumite, Ti2V3+O5(OH): new structure type, modularity and revised formula

The crystal structure of kyzylkumite, ideally Ti2V3+O5(OH), from the Sludyanka complex in South Baikal, Russia was solved and refined (including the hydrogen atom position) to an agreement index, R1, of 2.34 using X-ray diffraction data collected on a twinned crystal. Kyzylkumite crystallizes in space group P21/c, with a = 8.4787(1), b = 4.5624(1), c = 10.0330(1) Å, β = 93.174(1)°, V = 387.51(1) Å3 and Z = 4. Tivanite, TiV3+O3OH, and kyzylkumite have modular structures based on hexagonal close packing of oxygen, which are made up of rutile TiO2 and montroseite V3+O(OH) slices. In tivanite the rutile:montroseite ratio is 1:1, in kyzylkumite the ratio is 2:1. The montroseite module may be replaced by the isotypic paramontroseite V4+O2 module, which produces a phase with the formula Ti2V4+O6. In the metamorphic rocks of the Sludyanka complex, vanadium can be present as V4+ and V3+ within the same mineral (e.g. in batisivite, schreyerite and berdesinskiite). Kyzylkumite has a flexible composition with respect to the M4+/M3+ ratio. The relationship between kyzylkumite and a closely related Be-bearing kyzylkumite-like mineral with an orthorhombic norbergite-type structure from Byrud mine, Norway is discussed. Both minerals have similar X-ray powder diffraction patterns.

Statistical model-based segmentation of the proximal femur in digital antero-posterior (AP) pelvic radiographs.

PURPOSE    Segmentation of the proximal femur in digital antero-posterior (AP) pelvic radiographs is required to create a three-dimensional model of the hip joint for use in planning and treatment. However, manually extracting the femoral contour is tedious and prone to subjective bias, while automatic segmentation must accommodate poor image quality, anatomical structure overlap, and femur deformity. A new method was developed for femur segmentation in AP pelvic radiographs. METHODS    Using manual annotations on 100 AP pelvic radiographs, a statistical shape model (SSM) and a statistical appearance model (SAM) of the femur contour were constructed. The SSM and SAM were used to segment new AP pelvic radiographs with a three-stage approach. At initialization, the mean SSM model is coarsely registered to the femur in the AP radiograph through a scaled rigid registration. Mahalanobis distance defined on the SAM is employed as the search criteria for each annotated suggested landmark location. Dynamic programming was used to eliminate ambiguities. After all landmarks are assigned, a regularized non-rigid registration method deforms the current mean shape of SSM to produce a new segmentation of proximal femur. The second and third stages are iteratively executed to convergence. RESULTS    A set of 100 clinical AP pelvic radiographs (not used for training) were evaluated. The mean segmentation error was [Formula: see text], requiring [Formula: see text] s per case when implemented with Matlab. The influence of the initialization on segmentation results was tested by six clinicians, demonstrating no significance difference. CONCLUSIONS    A fast, robust and accurate method for femur segmentation in digital AP pelvic radiographs was developed by combining SSM and SAM with dynamic programming. This method can be extended to segmentation of other bony structures such as the pelvis.

Expectation Conditional Maximization-Based Deformable Shape Registration

This paper addresses the issue of matching statistical and non-rigid shapes, and introduces an Expectation Conditional Maximization-based deformable shape registration (ECM-DSR) algorithm. Similar to previous works, we cast the statistical and non-rigid shape registration problem into a missing data framework and handle the unknown correspondences with Gaussian Mixture Models (GMM). The registration problem is then solved by fitting the GMM centroids to the data. But unlike previous works where equal isotropic covariances are used, our new algorithm uses heteroscedastic covariances whose values are iteratively estimated from the data. A previously introduced virtual observation concept is adopted here to simplify the estimation of the registration parameters. Based on this concept, we derive closed-form solutions to estimate parameters for statistical or non-rigid shape registrations in each iteration. Our experiments conducted on synthesized and real data demonstrate that the ECM-DSR algorithm has various advantages over existing algorithms.