Approximation Classes for Adaptive Methods

* This work has been supported by the Office of Naval Research Contract Nr. N0014-91-J1343, the Army Research Office Contract Nr. DAAD 19-02-1-0028, the National Science Foundation grants DMS-0221642 and DMS-0200665, the Deutsche Forschungsgemeinschaft grant SFB 401, the IHP Network “Breaking Complexity” funded by the European Commission and the Alexan- der von Humboldt Foundation.

Theoretical investigations of intercellular communication in the microcirculation: Conducted responses, myoendothelial projections and endothelium derived hyperpolarizing factor

The contractile state of microcirculatory vessels is a major determinant of the blood pressure of the whole systemic circulation. Continuous bi-directional communication exists between the endothelial cells (ECs) and smooth muscle cells (SMCs) that regulates calcium (Ca2+) dynamics in these cells. This study presents theoretical approaches to understand some of the important and currently unresolved microcirculatory phenomena. ^ Agonist induced events at local sites have been shown to spread long distances in the microcirculation. We have developed a multicellular computational model by integrating detailed single EC and SMC models with gap junction and nitric oxide (NO) coupling to understand the mechanisms behind this effect. Simulations suggest that spreading vasodilation mainly occurs through Ca 2+ independent passive conduction of hyperpolarization in RMAs. Model predicts a superior role for intercellular diffusion of inositol (1,4,5)-trisphosphate (IP3) than Ca2+ in modulating the spreading response. ^ Endothelial derived signals are initiated even during vasoconstriction of stimulated SMCs by the movement of Ca2+ and/or IP3 into the EC which provide hyperpolarizing feedback to SMCs to counter the ongoing constriction. Myoendothelial projections (MPs) present in the ECs have been recently proposed to play a role in myoendothelial feedback. We have developed two models using compartmental and 2D finite element methods to examine the role of these MPs by adding a sub compartment in the EC to simulate MP with localization of intermediate conductance calcium activated potassium channels (IKCa) and IP3 receptors (IP 3R). Both models predicted IP3 mediated high Ca2+ gradients in the MP after SMC stimulation with limited global spread. This Ca 2+ transient generated a hyperpolarizing feedback of ∼ 2–3mV. ^ Endothelium derived hyperpolarizing factor (EDHF) is the dominant form of endothelial control of SMC constriction in the microcirculation. A number of factors have been proposed for the role of EDHF but no single pathway is agreed upon. We have examined the potential of myoendothelial gap junctions (MEGJs) and potassium (K+) accumulation as EDHF using two models (compartmental and 2D finite element). An extra compartment is added in SMC to simulate micro domains (MD) which have NaKα2 isoform sodium potassium pumps. Simulations predict that MEGJ coupling is much stronger in producing EDHF than alone K+ accumulation. On the contrary, K+ accumulation can alter other important parameters (EC V m, IKCa current) and inhibit its own release as well as EDHF conduction via MEGJs. The models developed in this study are essential building blocks for future models and provide important insights to the current understanding of myoendothelial feedback and EDHF.^

Algebraic Theory of Minimal Nondeterministic Finite Automata with Applications

Since the 1950s, the theory of deterministic and nondeterministic finite automata (DFAs and NFAs, respectively) has been a cornerstone of theoretical computer science. In this dissertation, our main object of study is minimal NFAs. In contrast with minimal DFAs, minimal NFAs are computationally challenging: first, there can be more than one minimal NFA recognizing a given language; second, the problem of converting an NFA to a minimal equivalent NFA is NP-hard, even for NFAs over a unary alphabet. Our study is based on the development of two main theories, inductive bases and partials, which in combination form the foundation for an incremental algorithm, ibas, to find minimal NFAs. An inductive basis is a collection of languages with the property that it can generate (through union) each of the left quotients of its elements. We prove a fundamental characterization theorem which says that a language can be recognized by an n-state NFA if and only if it can be generated by an n-element inductive basis. A partial is an incompletely-specified language. We say that an NFA recognizes a partial if its language extends the partial, meaning that the NFA's behavior is unconstrained on unspecified strings; it follows that a minimal NFA for a partial is also minimal for its language. We therefore direct our attention to minimal NFAs recognizing a given partial. Combining inductive bases and partials, we generalize our characterization theorem, showing that a partial can be recognized by an n-state NFA if and only if it can be generated by an n-element partial inductive basis. We apply our theory to develop and implement ibas, an incremental algorithm that finds minimal partial inductive bases generating a given partial. In the case of unary languages, ibas can often find minimal NFAs of up to 10 states in about an hour of computing time; with brute-force search this would require many trillions of years.

Symmetric Interior Penalty DG Methods for the Compressible Navier-Stokes Equations I: Method Formulation

In this article we consider the development of discontinuous Galerkin finite element methods for the numerical approximation of the compressible Navier-Stokes equations. For the discretization of the leading order terms, we propose employing the generalization of the symmetric version of the interior penalty method, originally developed for the numerical approximation of linear self-adjoint second-order elliptic partial differential equations. In order to solve the resulting system of nonlinear equations, we exploit a (damped) Newton-GMRES algorithm. Numerical experiments demonstrating the practical performance of the proposed discontinuous Galerkin method with higher-order polynomials are presented.

Discontinuous Galerkin Methods for the Biharmonic Problem

This work is concerned with the design and analysis of hp-version discontinuous Galerkin (DG) finite element methods for boundary-value problems involving the biharmonic operator. The first part extends the unified approach of Arnold, Brezzi, Cockburn & Marini (SIAM J. Numer. Anal. 39, 5 (2001/02), 1749-1779) developed for the Poisson problem, to the design of DG methods via an appropriate choice of numerical flux functions for fourth order problems; as an example we retrieve the interior penalty DG method developed by Suli & Mozolevski (Comput. Methods Appl. Mech. Engrg. 196, 13-16 (2007), 1851-1863). The second part of this work is concerned with a new a-priori error analysis of the hp-version interior penalty DG method, when the error is measured in terms of both the energy-norm and L2-norm, as well certain linear functionals of the solution, for elemental polynomial degrees $p\ge 2$. Also, provided that the solution is piecewise analytic in an open neighbourhood of each element, exponential convergence is also proven for the p-version of the DG method. The sharpness of the theoretical developments is illustrated by numerical experiments.

Population genomics diversity of Plasmodium falciparum in malaria patients attending Okelele Health Centre, Okelele, Ilorin, Kwara State, Nigeria.

Background: Plasmodium falciparum, the most dangerous malaria parasite species to humans remains an important public health concern in Okelele, a rural community in Ilorin, Kwara State, Nigeria. There is however little information about the genetic diversity of Plasmodium falciparum in Nigeria. Objective: To determine the population genomic diversity of Plasmodium falciparum in malaria patients attending Okelele Community Healthcare Centre, Okelele, Ilorin, Kwara State. Methods: In this study, 50 Plasmodium falciparum strains Merozoite Surface Protein 1, Merozoite Surface Protein 2 and Glutamate Rich Protein were analysed from Okelele Health Centre, Okelele, Ilorin, Nigeria. Genetic diversity of P. falciparum isolates were analysed from nested polymerase chain reactions (PCR) of the MSP-1 (K1, MAD 20 and RO33), MSP-2 (FC27 and 3D7) and Glutamate Rich Protein allelic families respectively. Results: Polyclonal infections were more in majority of the patients for MSP-1 allelic families while monoclonal infections were more for MSP-2 allelic families. Multiplicity of infection for MSP-1, MSP-2 and GLURP were 1.7, 1.8 and 2.05 respectively Conclusion: There is high genetic diversity in MSP – 2 and GLURP allelic families of Plasmodium falciparum isolates from Okelele Health Centre, Ilorin, Nigeria.

Validación de la simulación numérica del flujo bifásico hidrodinámico en sistemas de lecho fluido

Este trabajo tiene como objetivo la mejora en la validación de la simulación numérica del flujo bifásico característico del transporte de lecho fluido, mediante la formulación y desarrollo de un modelo numérico combinado Volúmenes Finitos - Elementos Finitos. Para ello se simula numéricamente el flujo de mezcla sólido-gas en una Cámara de Lecho Fluido, bajo implementación en código COMSOL, cuyos resultados son mejores comparativamente a un modelo basado en el método de Elementos Discretos implementado en código abierto MFIX. El problema fundamental de la modelización matemática del fenómeno de lecho fluido es la irregularidad del dominio, el acoplamiento de las variables en espacio y tiempo y, la no linealidad. En esta investigación se reformula apropiadamente las ecuaciones conservativas del fenómeno, tales que permitan obtener un problema variacional equivalente y solucionable numéricamente. Entonces; se define una ecuación de estado en función de la presión hidrodinámica y la fracción volumétrica de sólidos, quedando desacoplado el sistema en tres sub-problemas, garantizando así la existencia de solución del problema general. Una vez aproximados numéricamente ambos modelos, se comparan los resultados de donde se observa que el modelo materia del presente artículo, verifica de forma más eficaz las condiciones de mezcla óptima, reflejada en la calidad del burbujeo y velocidad de mezcla.

Novos registros e potencial epidemiológico de algumas espécies de mosquitos (Diptera, Culicidae), no Estado do Rio Grande do Sul

INTRODUÇÃO: A vigilância entomológica tem se mostrado uma importante estratégia de monitoramento da fauna de culicídeos com vistas a predizer o risco de exposição a espécies vetoras de patógenos. Esse trabalho apresenta uma lista de mosquitos identificados pela primeira vez no Rio Grande do Sul e discute o potencial epidemiológico de algumas espécies ocorrentes no Município de Maquiné com registros em outras regiões do Estado. MÉTODOS: Os mosquitos foram coletados com aspirador de Nasci e armadilhas CDC, entre dezembro de 2006 e dezembro de 2008, em área silvestre, rural e urbana do Município de Maquiné. RESULTADOS: Foram verificadas 55 espécies, das quais 22 são registradas pela primeira vez no estado e 10 são potencialmente vetoras do vírus Saint Louis, Oropouche, Aura, Trocara, Ilhéus, Rocio, Una, West Nile e encefalite equina do leste. CONCLUSÕES: Esses dados demonstram a importância da Vigilância Entomológica como ferramenta de informação e ação para a Vigilância em Saúde.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

A 2-D Delaunay refinement algorithm using an initial prerefinement from the boundary mesh

An improvement to the quality bidimensional Delaunay mesh generation algorithm, which combines the mesh refinement algorithms strategy of Ruppert and Shewchuk is proposed in this research. The developed technique uses diametral lenses criterion, introduced by L. P. Chew, with the purpose of eliminating the extremely obtuse triangles in the boundary mesh. This method splits the boundary segment and obtains an initial prerefinement, and thus reducing the number of necessary iterations to generate a high quality sequential triangulation. Moreover, it decreases the intensity of the communication and synchronization between subdomains in parallel mesh refinement.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

Topology optimization for designing strain-gauge load cells

Load cells are used extensively in engineering fields. This paper describes a novel structural optimization method for single- and multi-axis load cell structures. First, we briefly explain the topology optimization method that uses the solid isotropic material with penalization (SIMP) method. Next, we clarify the mechanical requirements and design specifications of the single- and multi-axis load cell structures, which are formulated as an objective function. In the case of multi-axis load cell structures, a methodology based on singular value decomposition is used. The sensitivities of the objective function with respect to the design variables are then formulated. On the basis of these formulations, an optimization algorithm is constructed using finite element methods and the method of moving asymptotes (MMA). Finally, we examine the characteristics of the optimization formulations and the resultant optimal configurations. We confirm the usefulness of our proposed methodology for the optimization of single- and multi-axis load cell structures.

Predictors of genital warts among women attending two family planning clinics

Objectives: To establish the prevalence and predictors of genital warts among healthy women presenting for contraceptive advice at two family planning clinics, one in a major Australian city and one in a country town in the same state. Methods: Consecutive consenting attendees (n = 1218)at two family planning clinics in Queensland completed a questionnaire and were examined for genital warts. Results: The point prevalence of visible genital warts was 3.3 per cent in the city clinic and 14.4 per cent in the country town. For half of these clients a finding of warts was unexpected, in that the client was unaware of their presence and presentation to the family planning clinic was not specifically for advice about sexually transmitted infections. The major predictor of a finding of warts was client age, with the highest prevalence in 20- to 25-year-olds. Warts were also commoner amongst smokers in the country town but not in Brisbane. However, no analysed sociodemographic variable predicted a finding of warts of which the client was not aware. Conclusions: Genital warts are common among young women presenting for contraceptive advice. Such women are often unaware that they have warts. Examination for genital warts should be a part of any routine examination of sexually active women, and medical practitioners should be aware of appropriate advice for patients who are found to have genital warts on routine examination.

Respiratory effects of tobacco smoking among young adults

Background. Respiratory symptoms associated with smoking habit seem to be age dependent. However, there are few reports about the effect of tobacco in young populations. The objective of this study was to analyze the effect of smoking on respiratory symptoms and lung function in 23- to 25-year-old adults in Brazil. This study had a cross-sectional design and included 2063 young people in the city of Ribeirao Preto, Sao Paulo State. Methods: Subjects completed a questionnaire used by the European Community Respiratory Health Survey and underwent spirometry and bronchial challenge test with methacholine. Multiple logistic regression analysis and multiple linear regression analysis were carried out to assess the association between smoking and respiratory symptoms, bronchial hyperresponsiveness, forced expiratory volume in 1 second (FEV1), and forced vital capacity (FVC), adjusted for confounding variables. Results: Prevalence of smoking habit was 17.2% with consumption (median) of 10 cigarettes per day (interquartile range 3-20). There was a significant association between smoking and respiratory symptoms. Smoking was associated to wheezing with odds ratio (95%Cl) of 6.11 (4.03-9.28) among those smoking :10 cigarettes per day and 3.36 (2.11-5.37) among those smoking <10 cigarettes per day. Associations were found for other respiratory symptoms. Smoking was associated with lower FEV1/FVC ratio. No association was detected between smoking and FEV1 or bronchial hyperresponsiveness. Conclusions: These findings highlight the early health consequences of smoking among young adults. These results prompt the necessity to elaborate urgent programs to reduce tobacco habit in young populations.

Modelo de especificação e controlo de execuções paralelas e distribuídas

O desenvolvimento actual de aplicações paralelas com processamento intensivo (HPC - High Performance Computing) para alojamento em computadores organizados em Cluster baseia-se muito no modelo de passagem de mensagens, do qual é de realçar os esforços de definição de standards, por exemplo, MPI - Message - Passing Interface. Por outro lado, com a generalização do paradigma de programação orientado aos objectos para ambientes distribuídos (Java RMI, .NET Remoting), existe a possibilidade de considerar que a execução de uma aplicação, de processamento paralelo e intensivo, pode ser decomposta em vários fluxos de execução paralela, em que cada fluxo é constituído por uma ou mais tarefas executadas no contexto de objectos distribuídos. Normalmente, em ambientes baseados em objectos distribuídos, a especificação, controlo e sincronização dos vários fluxos de execução paralela, é realizada de forma explicita e codificada num programa principal (hard-coded), dificultando possíveis e necessárias modificações posteriores. No entanto, existem, neste contexto, trabalhos que propõem uma abordagem de decomposição, seguindo o paradigma de workflow com interacções entre as tarefas por, entre outras, data-flow, control-flow, finite - state - machine. Este trabalho consistiu em propor e explorar um modelo de execução, sincronização e controlo de múltiplas tarefas, que permita de forma flexível desenhar aplicações de processamento intensivo, tirando partido da execução paralela de tarefas em diferentes máquinas. O modelo proposto e consequente implementação, num protótipo experimental, permite: especificar aplicações usando fluxos de execução; submeter fluxos para execução e controlar e monitorizar a execução desses fluxos. As tarefas envolvidas nos fluxos de execução podem executar-se num conjunto de recursos distribuídos. As principais características a realçar no modelo proposto, são a expansibilidade e o desacoplamento entre as diferentes componentes envolvidas na execução dos fluxos de execução. São ainda descritos casos de teste que permitiram validar o modelo e o protótipo implementado. Tendo consciência da necessidade de continuar no futuro esta linha de investigação, este trabalho é um contributo para demonstrar que o paradigma de workflow é adequado para expressar e executar, de forma paralela e distribuída, aplicações complexas de processamento intensivo.