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.

Lateral spreading and associated slope processes in fractured rock slabs

Landslides of the lateral spreading type, involving brittle geological units overlying ductile terrains, are a common occurrence in the sandstone and limestone plateaux of the northern Apennines of Italy. These instability phenomena can become particularly risky, when historical towns and cultural heritage sites built on the top of them are endangered. Neverthless, the mechanisms controlling the developing of related instabilities, i.e. toppling and rock falls, at the edges of rock plateaux are not fully understood yet. In addition, the groundwater flow path developing at the contact between the more permeable units, i.e. the jointed rock slab, and the relatively impermeable clay-rich units have not been already studied in details, even if they may play a role in this kind of instability processes, acting as eventual predisposing and/or triggering factors. Field survey, Terrestrial Laser Scanner and Close Range Photogrammetry techniques, laboratory tests on the involved materials, hydrogeological monitoring and modelling, displacements evaluation and stability analysis through continuum and discontinuum numerical codes have been performed on the San Leo case study, with the aim to bring further insights for the understanding and the assessment of the slope processes taking place in this geological context. The current research permitted to relate the aquifer behaviour of the rocky slab to slope instability processes. The aquifer hosted in the fractured slab leads to the development of perennial and ephemeral springs at the contact between the two units. The related piping erosion phenomena, together with slope processes in the clay-shales led to the progressive undermining of the slab. The cliff becomes progressively unstable due to undermining and undergoes large-scale landslides due to fall or topple.

Improving gene silencing of siRNAs via tricyclo-DNA modification

;Small interfering RNAs (siRNAs) can be exploited for the selective silencing of disease-related genes via the RNA interference (RNAi) machinery and therefore raise hope for future therapeutic applications. Especially chemically modified siRNAs are of interest as they are expected to convert lead siRNA sequences into effective drugs. To study the potential of tricyclo-DNA (tc-DNA) in this context we systematically incorporated tc-DNA units at various positions in a siRNA duplex targeted to the EGFP gene that was expressed in HeLa cells. Silencing activity was measured by FACS, mRNA levels were determined by RT-PCR and the biostability of the modifed siRNAs was determined in human serum. We found that modifications in the 3'-overhangs in both the sense and antisense strands were compatible with the RNAi machinery leading to similar activities compared to wild type (wt) siRNA. Additional modifications at the 3'-end, the 5'- end and in the center of the sense (passenger) strand were also well tolerated and did not compromise activity. Extensive modifications of the 3'- and the 5'-end in the antisense (guide) strand, however, abolished RNAi activity. Interestingly, modifications in the center of the duplex on both strands, corresponding to the position of the cleavage site by AGO2, increased efficacy relative to wt by a factor of 4 at the lowest concentrations (2 nM) investigated. In all cases, reduction of EGFP fluorescence was accompanied with a reduction of the EGFP mRNA level. Serum stability analysis further showed that 3'-overhang modifications only moderately increased stability while more extensive substitution by tc-DNA residues significantly enhanced biostability.

Cut slope design for the access to an underground copper mine

Slope stability analysis is a major area of research in geotechnical engineering. That being said, very little is written in the geotechnical engineering literature on the design of box-cuts. The goal of this thesis will be to investigate the proper design of a boxcuts, and to design a box-cut for access to an underground copper mine. Issues that need to be considered in the box-cut design include, long term dewatering design, slope stability analysis, and erosion control. The soils at the project site were extremely low permeability, as a result a system of ejectors was designed both to improve the stability of the slopes and prevent flooding. Based on the results of limit equilibrium analysis and finite element analysis, a slope design of two horizontal on one vertical was selection, with a rock fill buttress providing reinforcement. Finally, Michigan DOT standards for seeding were used to provide erosion control

Transformer modeling ATP : internal faults & high-frequency discretization

Power transformers are key components of the power grid and are also one of the most subjected to a variety of power system transients. The failure of a large transformer can cause severe monetary losses to a utility, thus adequate protection schemes are of great importance to avoid transformer damage and maximize the continuity of service. Computer modeling can be used as an efficient tool to improve the reliability of a transformer protective relay application. Unfortunately, transformer models presently available in commercial software lack completeness in the representation of several aspects such as internal winding faults, which is a common cause of transformer failure. It is also important to adequately represent the transformer at frequencies higher than the power frequency for a more accurate simulation of switching transients since these are a well known cause for the unwanted tripping of protective relays. This work develops new capabilities for the Hybrid Transformer Model (XFMR) implemented in ATPDraw to allow the representation of internal winding faults and slow-front transients up to 10 kHz. The new model can be developed using any of two sources of information: 1) test report data and 2) design data. When only test-report data is available, a higher-order leakage inductance matrix is created from standard measurements. If design information is available, a Finite Element Model is created to calculate the leakage parameters for the higher-order model. An analytical model is also implemented as an alternative to FEM modeling. Measurements on 15-kVA 240?/208Y V and 500-kVA 11430Y/235Y V distribution transformers were performed to validate the model. A transformer model that is valid for simulations for frequencies above the power frequency was developed after continuing the division of windings into multiple sections and including a higher-order capacitance matrix. Frequency-scan laboratory measurements were used to benchmark the simulations. Finally, a stability analysis of the higher-order model was made by analyzing the trapezoidal rule for numerical integration as used in ATP. Numerical damping was also added to suppress oscillations locally when discontinuities occurred in the solution. A maximum error magnitude of 7.84% was encountered in the simulated currents for different turn-to-ground and turn-to-turn faults. The FEM approach provided the most accurate means to determine the leakage parameters for the ATP model. The higher-order model was found to reproduce the short-circuit impedance acceptably up to about 10 kHz and the behavior at the first anti-resonant frequency was better matched with the measurements.

Influence of Mechanical and Thermal Boundary Conditions on Stabilizing/Destabilizing Mechanisms in Evaporating Liquid Films

Liquid films, evaporating or non-evaporating, are ubiquitous in nature and technology. The dynamics of evaporating liquid films is a study applicable in several industries such as water recovery, heat exchangers, crystal growth, drug design etc. The theory describing the dynamics of liquid films crosses several fields such as engineering, mathematics, material science, biophysics and volcanology to name a few. Interfacial instabilities typically manifest by the undulation of an interface from a presumed flat state or by the onset of a secondary flow state from a primary quiescent state or both. To study the instabilities affecting liquid films, an evaporating/non-evaporating Newtonian liquid film is subject to a perturbation. Numerical analysis is conducted on configurations of such liquid films being heated on solid surfaces in order to examine the various stabilizing and destabilizing mechanisms that can cause the formation of different convective structures. These convective structures have implications towards heat transfer that occurs via this process. Certain aspects of this research topic have not received attention, as will be obvious from the literature review. Static, horizontal liquid films on solid surfaces are examined for their resistance to long wave type instabilities via linear stability analysis, method of normal modes and finite difference methods. The spatiotemporal evolution equation, available in literature, describing the time evolution of a liquid film heated on a solid surface, is utilized to analyze various stabilizing/destabilizing mechanisms affecting evaporating and non-evaporating liquid films. The impact of these mechanisms on the film stability and structure for both buoyant and non-buoyant films will be examined by the variation of mechanical and thermal boundary conditions. Films evaporating in zero gravity are studied using the evolution equation. It is found that films that are stable to long wave type instabilities in terrestrial gravity are prone to destabilization via long wave instabilities in zero gravity.

Plurality of tree species responses to drought perturbation in Bornean tropical rain forest

Drought perturbation driven by the El Niño Southern Oscillation (ENSO) is a principal stochastic variable determining the dynamics of lowland rain forest in S.E. Asia. Mortality, recruitment and stem growth rates at Danum in Sabah (Malaysian Borneo) were recorded in two 4-ha plots (trees ≥ 10 cm gbh) for two periods, 1986–1996 and 1996–2001. Mortality and growth were also recorded in a sample of subplots for small trees (10 to <50 cm gbh) in two sub-periods, 1996–1999 and 1999–2001. Dynamics variables were employed to build indices of drought response for each of the 34 most abundant plot-level species (22 at the subplot level), these being interval-weighted percentage changes between periods and sub-periods. A significant yet complex effect of the strong 1997/1998 drought at the forest community level was shown by randomization procedures followed by multiple hypothesis testing. Despite a general resistance of the forest to drought, large and significant differences in short-term responses were apparent for several species. Using a diagrammatic form of stability analysis, different species showed immediate or lagged effects, high or low degrees of resilience or even oscillatory dynamics. In the context of the local topographic gradient, species’ responses define the newly termed perturbation response niche. The largest responses, particularly for recruitment and growth, were among the small trees, many of which are members of understorey taxa. The results bring with them a novel approach to understanding community dynamics: the kaleidoscopic complexity of idiosyncratic responses to stochastic perturbations suggests that plurality, rather than neutrality, of responses may be essential to understanding these tropical forests. The basis to the various responses lies with the mechanisms of tree-soil water relations which are physiologically predictable: the timing and intensity of the next drought, however, is not. To date, environmental stochasticity has been insufficiently incorporated into models of tropical forest dynamics, a step that might considerably improve the reality of theories about these globally important ecosystems.

NEW TARGET GENES FOR TUMOR SUPPRESSORS p53 AND p73 IN REGENERATING LIVER

The p53-family of proteins regulates expression of target genes during tissue development and differentiation. Within the p53-family, p53 and p73 have hepatic-specific functions in development and tumor suppression. Despite a growing list of p53/p73 target genes, very few of these have been studied in vivo, and the knowledge regarding functions of p53 and p73 in normal tissues remains limited. p53+/-p73+/- mice develop hepatocellular carcinoma (HCC), whereas overexpression of p53 in human HCC leads to tumor regression. However, the mechanism of p53/p73 function in liver remains poorly characterized. Here, the model of mouse liver regeneration is used to identify new target genes for p53/p73 in normal quiescent vs. proliferating cells. In response to surgical removal of ~2/3 of liver mass (partial hepatectomy, PH), the remaining hepatocytes exit G0 of cell cycle and undergo proliferation to reestablish liver mass. The hypothesis tested in this work is that p53/p73 functions in cell cycle arrest, apoptosis and senescence are repressed during liver regeneration, and reactivated at the end of the regenerative response. Chromatin immunoprecipitation (ChIP), with a p73-antibody, was used to probe arrayed genomic sequences (ChIP-chip) and uncover 158 potential targets of p73-regulation in normal liver. Global microarray analysis of mRNA levels, at T=0-48h following PH, revealed sets of genes that change expression during regeneration. Eighteen p73-bound genes changed expression after PH. Four of these genes, Foxo3, Jak1, Pea15, and Tuba1 have p53 response elements (p53REs), identified in silico within the upstream regulatory region. Forkhead transcription factor Foxo3 is the most responsive gene among transcription factors with altered expression during regenerative, cellular proliferation. p53 and p73 bind a Foxo3 p53RE and maintain active expression in quiescent liver. During liver regeneration, binding of p53 and p73, recruitment of acetyltransferase p300, and an active chromatin structure of Foxo3 are disrupted, alongside loss of Foxo3 expression. These parameters of Foxo3 regulation are reestablished at completion of liver growth and regeneration, supporting a temporary suspension of p53 and p73 regulatory functions in normal cells during tissue regeneration.

Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

Metabolomics reveals trichloroacetate as a major contributor to trichloroethylene-induced metabolic alterations in mouse urine and serum

Trichloroethylene (TCE)-induced liver toxicity and carcinogenesis is believed to be mediated in part by activation of the peroxisome proliferator-activated receptor α (PPARα). However, the contribution of the two TCE metabolites, dichloroacetate (DCA) and trichloroacetate (TCA) to the toxicity of TCE, remains unclear. The aim of the present study was to determine the metabolite profiles in serum and urine upon exposure of mice to TCE, to aid in determining the metabolic response to TCE exposure and the contribution of DCA and TCA to TCE toxicity. C57BL/6 mice were administered TCE, TCA, or DCA, and urine and serum subjected to ultra-performance liquid chromatography coupled with electrospray ionization quadrupole time-of-flight mass spectrometry (UPLC-ESI-QTOFMS)-based global metabolomics analysis. The ions were identified through searching metabolomics databases and by comparison with authentic standards, and quantitated using multiple reactions monitoring. Quantitative polymerase chain reaction of mRNA, biochemical analysis, and liver histology were also performed. TCE exposure resulted in a decrease in urine of metabolites involved in fatty acid metabolism, resulting from altered expression of PPARα target genes. TCE treatment also induced altered phospholipid homeostasis in serum, as revealed by increased serum lysophosphatidylcholine 18:0 and 18:1, and phosphatidylcholine metabolites. TCA administration revealed similar metabolite profiles in urine and serum upon TCE exposure, which correlated with a more robust induction of PPARα target gene expression associated with TCA than DCA treatment. These data show the metabolic response to TCE exposure and demonstrate that TCA is the major contributor to TCE-induced metabolite alterations observed in urine and serum.

What can a unicellular parasite tell us about mitochondrial biogenesis?

In the unicellular parasite Trypanosoma brucei, as in other eukaryotes, more than 95% of all mitochondrial proteins are imported from the cytosol. The recently characterized multisubunit ATOM complex, the functional analogue of the TOM complex of yeast, mediates import of essentially all proteins across the outer mitochondrial membrane in T. brucei. Moreover, an additional protein termed pATOM36, which is loosely associated with the ATOM complex, has been implicated in the import of only a subset of mitochondrial proteins. Here we have investigated more precisely which role pATOM36 plays in mitochondrial protein import. RNAi mediated ablation of pATOM36 specifically depletes a subset of outer mitochondrial membrane proteins including ATOM complex subunits and as a consequence results in the collapse of the ATOM complex as shown by Blue native PAGE. In addition, a SILAC-based global proteomic analysis of uninduced and induced pATOM36 RNAi cells together with in vitro import experiments suggest that pATOM36 might be a novel protein import factor acting on a subset of alpha-helically anchored mitochondrial outer membrane proteins. Identification of pATOM36 interaction partners by co-immunoprecipitation together with immunofluorescence analysis shows that unexpectedly a fraction of the protein is associated with the tripartite attachment complex (TAC). This complex is essential for proper inheritance of the mitochondrial DNA in T. brucei. It forms a physical connection between the single unit mitochondrial DNA and the basal body of the flagellum that is stable throughout the cell cycle. Thus, pATOM36 simultaneously mediates ATOM assembly, and thus protein import, as well as mitochondrial DNA inheritance since it is an essential component of the TAC.

Kinetoplastid-specific pATOM36 mediates membrane insertion of a subset of mitochondrial outer membrane proteins

In trypanosomes, as in other eukaryotes, more than 95% of all mitochondrial proteins are imported into the mitochondrion. The recently characterized multisubunit ATOM complex mediates import of essentially all proteins across the outer mitochondrial membrane in T. brucei. Moreover, an additional protein termed pATOM36, which is loosely associated with the ATOM complex, has been implicated in the import of only a subset of mitochondrial matrix proteins. Here we have investigated more precisely which role pATOM36 plays in mitochondrial protein import. RNAi mediated ablation of pATOM36 specifically depletes a subset of ATOM complex subunits and as a consequence results in the collapse of the ATOM complex as shown by Blue native PAGE. In addition, a SILAC-based global proteomic analysis of uninduced and induced pATOM36 RNAi cells together with in vitro import experiments suggest that pATOM36 might be a novel protein insertase acting on a subset of alpha-helically anchored mitochondrial outer membrane proteins. Identification of pATOM36 interaction partners by co-immunoprecipitation together with immunofluorescence analysis furthermore shows that unexpectedly a fraction of the protein is associated with the tripartite attachment complex (TAC). This complex is essential for proper inheritance of the kDNA as it forms a physical connection between the kDNA and the basal body of the flagellum throughout the cell cycle. Thus, the presence of pATOM36 in the TAC provides an exciting link between mitochondrial protein import and kDNA inheritance.

pATOM36 mediates mitochondrial protein import and mitochondrial DNA inheritance

The multisubunit ATOM complex mediates import of essentially all proteins across the outer mitochondrial membrane in T. brucei. Moreover, an additional protein termed pATOM36, which is loosely associated with the ATOM complex, has been implicated in the import of only a subset of mitochondrial matrix proteins. Here we have investigated more precisely which role pATOM36 plays in mitochondrial protein import. RNAi mediated ablation of pATOM36 specifically depletes a subset of ATOM complex subunits and as a consequence results in the collapse of the ATOM complex as shown by Blue native PAGE. In addition, a SILAC-based global proteomic analysis of uninduced and induced pATOM36 RNAi cells together with in vitro import experiments suggest that pATOM36 might be a novel protein insertase acting on a subset of alpha-helically anchored mitochondrial outer membrane proteins. Identification of pATOM36 interaction partners by co-immunoprecipitation together with immunofluorescence analysis furthermore shows that unexpectedly a fraction of the protein is associated with the tripartite attachment complex (TAC). This complex is essential for proper inheritance of the mtDNA; also called kinetoplast or kDNA; as it forms a physical connection between the kDNA and the basal body of the single flagellum throughout the cell cycle. Thus, the presence of pATOM36 in the TAC provides an exciting link between mitochondrial protein import and kDNA inheritance.

Efecto de un nuevo xanthanólido sesquiterpeno sobre la activación de mastocitos inducida por neurópeptidos pro-inflamatorios

Los mastocitos son células del tejido conectivo que participan en la génesis y modulación de las respuestas inflamatorias celulares. En trabajos previos hemos demostrado que xanthatina (xanthanólido sesquiterpeno aislado de Xanthium cavanillesii Schouw) inhibe la activación de mastocitos inducida por secretagogos experimentales. Sin embargo, se desconoce su efecto sobre la activación de mastocitos inducida por estímulos fisiopatológicos. Estos estímulos incluyen, entre otros, los neuropéptidos pro-inflamatorios sustancia P y neurotensina, responsables de una de las principales vías de inflamación neurogénica. El objetivo del presente trabajo fue estudiar el efecto de xanthatina sobre la activación de mastocitos inducida por sustancia P y neurotensina. Mastocitos peritoneales de rata se incubaron con: 1) PBS (basal); 2) sustancia P (100 Fm); 3) neurotensina (50 Fm); 4) xanthatina (8-320 Fm)+sustancia P; 5) xanthatina (8-320 Fm)+neurotensina. La viabilidad de los mastocitos se evaluó con azul tripán. En las soluciones de incubación se cuantificó serotonina liberada (marcador de activación). En las células se cuantificó serotonina remanente (no liberada) y se analizó la morfología celular por microscopía óptica y electrónica de transmisión. Tratamiento estadístico: ANOVA-1 y Tukey-Kramer. La incubación de mastocitos con xanthatina inhibió (P<0,01), en forma dosisdependiente, la liberación de serotonina inducida por sustancia P y neurotensina, sin modificar la viabilidad celular. Los mastocitos tratados con neuropéptidos mostraron características morfológicas de degranulación, mientras que la morfología de los mastocitos tratados con xanthatina+neuropéptido fue semejante a los basales. En conclusión, xanthatina inhibe la activación de mastocitos inducida por sustancia P y por neurotensina. Este sesquiterpeno podría representar una nueva alternativa en el tratamiento de las inflamaciones neurogénicas.

Grain-size distribution and frictional behaviour of sediment cores from the Nice Slope

The upper shelf of the landslide-prone Ligurian Margin (Western Mediterranean Sea) off Nice well-known for the 1979 Airport Landslide is a natural laboratory to study preconditioning factors and trigger mechanisms for submarine landslides. For this study low-stress ring shear experiments have been carried out on a variety of sediments from >50 gravity cores to characterise the velocity-dependent frictional behaviour. Mean values of the peak coefficient of friction vary from 0.46 for clay-dominated samples (53 % clay, 46 % silt, 1 %) sand up to 0.76 for coarse-grained sediments (26 % clay, 57 % silt, 17 % sand). The majority of the sediments tested show velocity strengthening regardless of the grain size distribution. For clayey sediments the peak and residual cohesive strength increases with increasing normal stress, with values from 1.3 to 10.6 kPa and up to 25 % of all strength supported by cohesive forces in the shallowmost samples. A pseudo-static slope stability analysis reveals that the different lithologies (even clay-rich material with clay content >=50 %) tested are stable up to slope angles <26° under quasi-drained conditions.