Untersuchungen zur Genauigkeit adaptiver unstetiger Galerkin-Simulationen mit Hilfe von Luftblasen-Testfällen

Allgemein erlaubt adaptive Gitterverfeinerung eine Steigerung der Effizienz numerischer Simulationen ohne dabei die Genauigkeit des Ergebnisses signifikant zu verschlechtern. Es ist jedoch noch nicht erforscht, in welchen Bereichen des Rechengebietes die räumliche Auflösung tatsächlich vergröbert werden kann, ohne die Genauigkeit des Ergebnisses signifikant zu beeinflussen. Diese Frage wird hier für ein konkretes Beispiel von trockener atmosphärischer Konvektion untersucht, nämlich der Simulation von warmen Luftblasen. Zu diesem Zweck wird ein neuartiges numerisches Modell entwickelt, das auf diese spezielle Anwendung ausgerichtet ist. Die kompressiblen Euler-Gleichungen werden mit einer unstetigen Galerkin Methode gelöst. Die Zeitintegration geschieht mit einer semi-implizite Methode und die dynamische Adaptivität verwendet raumfüllende Kurven mit Hilfe der Funktionsbibliothek AMATOS. Das numerische Modell wird validiert mit Hilfe einer Konvergenzstudie und fünf Standard-Testfällen. Eine Methode zum Vergleich der Genauigkeit von Simulationen mit verschiedenen Verfeinerungsgebieten wird eingeführt, die ohne das Vorhandensein einer exakten Lösung auskommt. Im Wesentlichen geschieht dies durch den Vergleich von Eigenschaften der Lösung, die stark von der verwendeten räumlichen Auflösung abhängen. Im Fall einer aufsteigenden Warmluftblase ist der zusätzliche numerische Fehler durch die Verwendung der Adaptivität kleiner als 1% des gesamten numerischen Fehlers, wenn die adaptive Simulation mehr als 50% der Elemente einer uniformen hoch-aufgelösten Simulation verwendet. Entsprechend ist die adaptive Simulation fast doppelt so schnell wie die uniforme Simulation.

Estudio y desarrollo de un preproceso basado en la difusión no lineal para la segmentación en imágenes

En esta Tesis Doctoral se aborda la utilización de filtros de difusión no lineal para obtener imágenes constantes a trozos como paso previo al proceso de segmentación. En una primera parte se propone un formulación intrínseca para la ecuación de difusión no lineal que proporcione las condiciones de diseño necesarias sobre los filtros de difusión. A partir del marco teórico propuesto, se proporciona una nueva familia de difusividades; éstas son obtenidas a partir de técnicas de difusión no lineal relacionadas con los procesos de difusión regresivos. El objetivo es descomponer la imagen en regiones cerradas que sean homogéneas en sus niveles de grises sin contornos difusos. Asimismo, se prueba que la función de difusividad propuesta satisface las condiciones de un correcto planteamiento semi-discreto. Esto muestra que mediante el esquema semi-implícito habitualmente utilizado, realmente se hace un proceso de difusión no lineal directa, en lugar de difusión inversa, conectando con proceso de preservación de bordes. Bajo estas condiciones establecidas, se plantea un criterio de parada para el proceso de difusión, para obtener imágenes constantes a trozos con un bajo coste computacional. Una vez aplicado todo el proceso al caso unidimensional, se extienden los resultados teóricos, al caso de imágenes en 2D y 3D. Para el caso en 3D, se detalla el esquema numérico para el problema evolutivo no lineal, con condiciones de contorno Neumann homogéneas. Finalmente, se prueba el filtro propuesto para imágenes reales en 2D y 3D y se ilustran los resultados de la difusividad propuesta como método para obtener imágenes constantes a trozos. En el caso de imágenes 3D, se aborda la problemática del proceso previo a la segmentación del hígado, mediante imágenes reales provenientes de Tomografías Axiales Computarizadas (TAC). En ese caso, se obtienen resultados sobre la estimación de los parámetros de la función de difusividad propuesta. This Ph.D. Thesis deals with the case of using nonlinear diffusion filters to obtain piecewise constant images as a previous process for segmentation techniques. I have first shown an intrinsic formulation for the nonlinear diffusion equation to provide some design conditions on the diffusion filters. According to this theoretical framework, I have proposed a new family of diffusivities; they are obtained from nonlinear diffusion techniques and are related with backward diffusion. Their goal is to split the image in closed contours with a homogenized grey intensity inside and with no blurred edges. It has also proved that the proposed filters satisfy the well-posedness semi-discrete and full discrete scale-space requirements. This shows that by using semi-implicit schemes, a forward nonlinear diffusion equation is solved, instead of a backward nonlinear diffusion equation, connecting with an edgepreserving process. Under the conditions established for the diffusivity and using a stopping criterion I for the diffusion time, I have obtained piecewise constant images with a low computational effort. The whole process in the one-dimensional case is extended to the case where 2D and 3D theoretical results are applied to real images. For 3D, develops in detail the numerical scheme for nonlinear evolutionary problem with homogeneous Neumann boundary conditions. Finally, I have tested the proposed filter with real images for 2D and 3D and I have illustrated the effects of the proposed diffusivity function as a method to get piecewise constant images. For 3D I have developed a preprocess for liver segmentation with real images from CT (Computerized Tomography). In this case, I have obtained results on the estimation of the parameters of the given diffusivity function.

Classifying GABAergic interneurons with semi-supervised projected model-based clustering

Objectives: A recently introduced pragmatic scheme promises to be a useful catalog of interneuron names.We sought to automatically classify digitally reconstructed interneuronal morphologies according tothis scheme. Simultaneously, we sought to discover possible subtypes of these types that might emergeduring automatic classification (clustering). We also investigated which morphometric properties weremost relevant for this classification.Materials and methods: A set of 118 digitally reconstructed interneuronal morphologies classified into thecommon basket (CB), horse-tail (HT), large basket (LB), and Martinotti (MA) interneuron types by 42 of theworld?s leading neuroscientists, quantified by five simple morphometric properties of the axon and fourof the dendrites. We labeled each neuron with the type most commonly assigned to it by the experts. Wethen removed this class information for each type separately, and applied semi-supervised clustering tothose cells (keeping the others? cluster membership fixed), to assess separation from other types and lookfor the formation of new groups (subtypes). We performed this same experiment unlabeling the cells oftwo types at a time, and of half the cells of a single type at a time. The clustering model is a finite mixtureof Gaussians which we adapted for the estimation of local (per-cluster) feature relevance. We performedthe described experiments on three different subsets of the data, formed according to how many expertsagreed on type membership: at least 18 experts (the full data set), at least 21 (73 neurons), and at least26 (47 neurons).Results: Interneurons with more reliable type labels were classified more accurately. We classified HTcells with 100% accuracy, MA cells with 73% accuracy, and CB and LB cells with 56% and 58% accuracy,respectively. We identified three subtypes of the MA type, one subtype of CB and LB types each, andno subtypes of HT (it was a single, homogeneous type). We got maximum (adapted) Silhouette widthand ARI values of 1, 0.83, 0.79, and 0.42, when unlabeling the HT, CB, LB, and MA types, respectively,confirming the quality of the formed cluster solutions. The subtypes identified when unlabeling a singletype also emerged when unlabeling two types at a time, confirming their validity. Axonal morphometricproperties were more relevant that dendritic ones, with the axonal polar histogram length in the [pi, 2pi) angle interval being particularly useful.Conclusions: The applied semi-supervised clustering method can accurately discriminate among CB, HT, LB, and MA interneuron types while discovering potential subtypes, and is therefore useful for neuronal classification. The discovery of potential subtypes suggests that some of these types are more heteroge-neous that previously thought. Finally, axonal variables seem to be more relevant than dendritic ones fordistinguishing among the CB, HT, LB, and MA interneuron types.

Graph-based Semi-Supervised Learning in Acoustic Modeling for Automatic Speech Recognition

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

Semi-Smooth Newton Methods for the Time Optimal Control of Nonautonomous Ordinary Differential Equations

AMS Subj. Classification: 49J15, 49M15

Measurement of carotenoid levels in human serum and a catalog of the lutein conformation populations from semi-empirical calculations

Lutein is a principal constituent of the human macular pigment. This study is composed of two projects. The first studies the conformational geometries of lutein and its potential adaptability in biological systems. The second is a study of the response of human subjects to lutein supplements. Using semi-empirical parametric method 3 (PM3) and density functional theory with the B3LYP/6-31G* basis set, the relative energies of s- cis conformers of lutein were determined. All 512 s-cis conformers were calculated with PM3. A smaller, representative group was also studied using density functional theory. PM3 results were correlated systematically to B3LYP values and this enables the results to be calibrated. The relative energies of the conformers range from 1-30 kcal/mole, and many are dynamically accessible at normal temperatures. Four commercial formulations containing lutein were studied. The serum and macular pigment (MP) responses of human subjects to these lutein supplements with doses of 9 or 20 mg/day were measured, relative to a placebo, over a six month period. In each instance, lutein levels in serum increased and correlated with MP increases. The results demonstrate that responses are significantly dependent upon formulation and that components other than lutein have an important influence serum response.

All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients

In this talk, we propose an all regime Lagrange-Projection like numerical scheme for the gas dynamics equations. By all regime, we mean that the numerical scheme is able to compute accurate approximate solutions with an under-resolved discretization with respect to the Mach number M, i.e. such that the ratio between the Mach number M and the mesh size or the time step is small with respect to 1. The key idea is to decouple acoustic and transport phenomenon and then alter the numerical flux in the acoustic approximation to obtain a uniform truncation error in term of M. This modified scheme is conservative and endowed with good stability properties with respect to the positivity of the density and the internal energy. A discrete entropy inequality under a condition on the modification is obtained thanks to a reinterpretation of the modified scheme in the Harten Lax and van Leer formalism. A natural extension to multi-dimensional problems discretized over unstructured mesh is proposed. Then a simple and efficient semi implicit scheme is also proposed. The resulting scheme is stable under a CFL condition driven by the (slow) material waves and not by the (fast) acoustic waves and so verifies the all regime property. Numerical evidences are proposed and show the ability of the scheme to deal with tests where the flow regime may vary from low to high Mach values.

A new approach for in vitro bioassay to measure tannin biological effects based on a gas production technique

The aim of this paper was to study a method based on gas production technique to measure the biological effects of tannins on rumen fermentation. Six feeds were used as fermentation substrates in a semi-automated gas method: feed A - aroeira (Astronium urundeuva); feed B - jurema preta (Mimosa hostilis), feed C - sorghum grains (Sorghum bicolor); feed D - Tifton-85 (Cynodon sp.); and two others prepared mixing 450 g sorghum leaves, 450 g concentrate (maize and soybean meal) and 100 g either of acacia (Acacia mearnsii) tannin extract (feed E) or quebracho (Schinopsis lorentzii) tannin extract (feed F) per kg (w:w). Three assays were carried out to standardize the bioassay for tannins. The first assay compared two binding agents (polyethylene glycol - PEG - and polyvinyl polypirrolidone - PVPP) to attenuate the tannin effects. The complex formed by PEG and tannins showed to be more stable than PVPP and tannins. Then, in the second assay, PEG was used as binding agent, and this assay was done to evaluate levels of PEG (0, 500, 750, 1000 and 1250 mg/g DM) to minimize the tannin effect. All the tested levels of PEG produced a response to evaluate tannin effects but the best response was for dose of 1000 mg/g DM. Using this dose of PEG, the final assay was carried out to test three compounds (tannic acid, quebracho extract and acacia extract) to establish a curve of biological equivalent effect of tannins. For this, five levels of each compound were added to I g of a standard feed (Lucerne hay). The equivalent effect showed not to be directly related to the chemical analysis for tannins. It was shown that different sources of tannins had different activities or reactivities. The curves of biological equivalence can provide information about tannin reactivity and its use seems to be important as an additional factor for chemical analysis. (C) 2007 Elsevier B.V. All rights reserved.

Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

Predictor-corrector methods of Runge-Kutta type for stochastic differential equations

In this paper we construct predictor-corrector (PC) methods based on the trivial predictor and stochastic implicit Runge-Kutta (RK) correctors for solving stochastic differential equations. Using the colored rooted tree theory and stochastic B-series, the order condition theorem is derived for constructing stochastic RK methods based on PC implementations. We also present detailed order conditions of the PC methods using stochastic implicit RK correctors with strong global order 1.0 and 1.5. A two-stage implicit RK method with strong global order 1.0 and a four-stage implicit RK method with strong global order 1.5 used as the correctors are constructed in this paper. The mean-square stability properties and numerical results of the PC methods based on these two implicit RK correctors are reported.

Nonlinear mixture model for hyperspectral unmixing

Proceedings of International Conference - SPIE 7477, Image and Signal Processing for Remote Sensing XV - 28 September 2009

The role of the efflux mechanisms in multidrug resistance in Mycobacterium tuberculosis

Abstract The emergence of multi and extensively drug resistant tuberculosis (MDRTB and XDRTB) has increased the concern of public health authorities around the world. The World Health Organization has defined MDRTB as tuberculosis (TB) caused by organisms resistant to at least isoniazid and rifampicin, the main first-line drugs used in TB therapy, whereas XDRTB refers to TB resistant not only to isoniazid and rifampicin, but also to a fluoroquinolone and to at least one of the three injectable second-line drugs, kanamycin, amikacin and capreomycin. Resistance in Mycobacterium tuberculosis is mainly due to the occurrence of spontaneous mutations and followed by selection of mutants by subsequent treatment. However, some resistant clinical isolates do not present mutations in any genes associated with resistance to a given antibiotic, which suggests that other mechanism(s) are involved in the development of drug resistance, namely the presence of efflux pump systems that extrude the drug to the exterior of the cell, preventing access to its target. Increased efflux activity can occur in response to prolonged exposure to subinhibitory concentrations of anti-TB drugs, a situation that may result from inadequate TB therapy. The inhibition of efflux activity with a non-antibiotic inhibitor may restore activity of an antibiotic subject to efflux and thus provide a way to enhance the activity of current anti-TB drugs. The work described in this thesis foccus on the study of efflux mechanisms in the development of multidrug resistance in M. tuberculosis and how phenotypic resistance, mediated by efflux pumps, correlates with genetic resistance. In order to accomplish this goal, several experimental protocols were developed using biological models such as Escherichia coli, the fast growing mycobacteria Mycobacterium smegmatis, and Mycobacterium avium, before their application to M. tuberculosis. This approach allowed the study of the mechanisms that result in the physiological adaptation of E. coli to subinhibitory concentrations of tetracycline (Chapter II), the development of a fluorometric method that allows the detection and quantification of efflux of ethidium bromide (Chapter III), the characterization of the ethidium bromide transport in M. smegmatis (Chapter IV) and the contribution of efflux activity to macrolide resistance in Mycobacterium avium complex (Chapter V). Finally, the methods developed allowed the study of the role of efflux pumps in M. tuberculosis strains induced to isoniazid resistance (Chapter VI). By this manner, in Chapter II it was possible to observe that the physiological adaptation of E. coli to tetracycline results from an interplay between events at the genetic level and protein folding that decrease permeability of the cell envelope and increase efflux pump activity. Furthermore, Chapter III describes the development of a semi-automated fluorometric method that allowed the correlation of this efflux activity with the transport kinetics of ethidium bromide (a known efflux pump substrate) in E. coli and the identification of efflux inhibitors. Concerning M. smegmatis, we have compared the wild-type M. smegmatis mc2155 with knockout mutants for LfrA and MspA for their ability to transport ethidium bromide. The results presented in Chapter IV showed that MspA, the major porin in M. smegmatis, plays an important role in the entrance of ethidium bromide and antibiotics into the cell and that efflux via the LfrA pump is involved in low-level resistance to these compounds in M. smegmatis. Chapter V describes the study of the contribution of efflux pumps to macrolide resistance in clinical M. avium complex isolates. It was demonstrated that resistance to clarithromycin was significantly reduced in the presence of efflux inhibitors such as thioridazine, chlorpromazine and verapamil. These same inhibitors decreased efflux of ethidium bromide and increased the retention of [14C]-erythromycin in these isolates. Finaly, the methods developed with the experimental models mentioned above allowed the study of the role of efflux pumps on M. tuberculosis strains induced to isoniazid resistance. This is described in Chapter VI of this Thesis, where it is demonstrated that induced resistance to isoniazid does not involve mutations in any of the genes known to be associated with isoniazid resistance, but an efflux system that is sensitive to efflux inhibitors. These inhibitors decreased the efflux of ethidium bromide and also reduced the minimum inhibitory concentration of isoniazid in these strains. Moreover, expression analysis showed overexpression of genes that code for efflux pumps in the induced strains relatively to the non-induced parental strains. In conclusion, the work described in this thesis demonstrates that efflux pumps play an important role in the development of drug resistance, namely in mycobacteria. A strategy to overcome efflux-mediated resistance may consist on the use of compounds that inhibit efflux activity, restoring the activity of antimicrobials that are efflux pump substrates, a useful approach particularly in TB where the most effective treatment regimens are becoming uneffective due to the increase of MDRTB/XDRTB.

The inhibitory activity of synthetic compounds and ions against transporters of multi-drug resistant bacteria

RESUMO: Sessenta e três derivados de hidantoína foram utilizados para avaliar possíveis efeitos de modulação na actividade das bombas de efluxo (BE) na Salmonella NCTC 13349 utilizando um método fluorimétrico semi-automático. Nenhum dos compostos apresentaram actividade anti-bacteriana até concentrações de 240 mg/L. Entre todos os compostos, SZ-7 demonstrou possuir propriedades de modulação de effluxo na presença de glucose. Para testar esta actividade, estirpes de Salmonella resistentes à ciprofloxacina, induzidas a elevados níveis de resistência com sobre-expressão de BE, foram expostas ao SZ-7. Este derivado afectou a susceptibilidade das estirpes à ciprofloxacina. Uma vez que os 63 compostos estudados apresentaram pouco efeito inibitório /cumulativo, apesar de serem conhecidos pelos seus efeitos moduladores de BE-dependentes de iões em eucariotas, foi questionado o papel dos iões na regulação de BE bacterianas, que poderão influenciar a eficácia de novos compostos. Para este estudo, utilizamos a Escherichia coli AG100 como modelo, devido ao extenso conhecimento no que respeita a estrutura e actividade das BE. Devido à importância de iões de cálcio (Ca2+) nos canais de transporte membranar e na actividade de ATPases, a sua actividade na modulação do efluxo foi investigada. De resultados anteriormente obtidos concluiu-se que a pH 5 o efluxo é independente de energia metabólica; contudo, a pH 8 é absolutamente dependente, sendo que o Ca2+ é indispensável para manter a actividade das ATPases bacterianas. A acumulação/effluxo de EtBr pela E. coli AG100 foi determinada na presença/ausência de Ca2+, clorpromazina (inibidor de ligação de Ca2+ a proteínas), e ácido etilenodiamino tetra-acético (quelante de Ca2+). Acumulação/effluxo aumentou a pH 8, contudo o Ca2+ reverte estes efeitos evidenciando a sua importância no funcionamento das BE bacterianas. Em resumo este trabalho colocou em evidência que muitos aspectos bioquímicos e bioenergéticos devem ser tomados em consideração no estudo da resistência bacteriana mediada por BE.------- ABSTRACT: Sixty-three hydantoin derivatives were evaluated for their modulating effects on efflux pump (EP) activity of Salmonella NCTC 13349 utilizing a semi-automatic fluorometric method. None of the compounds presented antibacterial activities at concentrations as high as 240 mg/L. Among all compounds, SZ-7 showed possible efflux modulating activity in the presence of glucose, indicative of a potential EP inhibitor. To verify its potential effects, ciprofloxacin-resistant Salmonella strains, induced to high level resistance with over-expressing EPs, were exposed to SZ-7. This derivative affected the susceptibility of the ciprofloxacin-resistant strains. Since the 63 compounds studied had very low inhibitory/accumulative effects, even though their known for being efficient in modulating ion-driven eukaryotic EPs, we questioned whether ions had a leading role in regulating bacterial EPs, influencing the effectiveness of new compounds. For this study we used Escherichia coli AG100 as a model, due to the extensive knowledge on its EPs structure and activity. Owing the importance of calcium ions (Ca2+) for membrane transport channels and activity of ATPases, the role of Ca2+ was investigated. From previous results we concluded that at pH 5 efflux is independent of metabolic energy; however, at pH 8 it is entirely dependent of metabolic energy and the Ca2+ ions are essential to maintain the activity of bacterial ATPases. Accumulation and efflux of ethidium bromide (EtBr) by E. coli AG100 was determined in the presence and absence of Ca2+, chlorpromazine (inhibitor of Ca2+-binding to proteins), and ethylenediaminetetraacetic acid (Ca2+ chelator). Accumulation of EtBr increased at pH 8; however Ca2+ reversed these effects providing information as to the importance of this ion in the regulation of bacterial EP systems. Overall this work puts in evidence that many biochemical and bioenergetic aspects related to the strains physiology need to be taken into consideration in bacterial drug resistance mediated by EPs.

Convergence of the mass-transport steepest descent scheme for the sub-critical Patlak-Keller-Segel model

Variational steepest descent approximation schemes for the modified Patlak-Keller-Segel equation with a logarithmic interaction kernel in any dimension are considered. We prove the convergence of the suitably interpolated in time implicit Euler scheme, defined in terms of the Euclidean Wasserstein distance, associated to this equation for sub-critical masses. As a consequence, we recover the recent result about the global in time existence of weak-solutions to the modified Patlak-Keller-Segel equation for the logarithmic interaction kernel in any dimension in the sub-critical case. Moreover, we show how this method performs numerically in one dimension. In this particular case, this numerical scheme corresponds to a standard implicit Euler method for the pseudo-inverse of the cumulative distribution function. We demonstrate its capabilities to reproduce easily without the need of mesh-refinement the blow-up of solutions for super-critical masses.

An objective anisotropic elastic plastic model and Algorithm applicable to bone mechanics

The implicit projection algorithm of isotropic plasticity is extended to an objective anisotropic elastic perfectly plastic model. The recursion formula developed to project the trial stress on the yield surface, is applicable to any non linear elastic law and any plastic yield function.A curvilinear transverse isotropic model based on a quadratic elastic potential and on Hill's quadratic yield criterion is then developed and implemented in a computer program for bone mechanics perspectives. The paper concludes with a numerical study of a schematic bone-prosthesis system to illustrate the potential of the model.