Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

Stochastic stability for Markovian jump linear systems associated with a finite number of jump times

This paper deals with a stochastic stability concept for discrete-time Markovian jump linear systems. The random jump parameter is associated to changes between the system operation modes due to failures or repairs, which can be well described by an underlying finite-state Markov chain. In the model studied, a fixed number of failures or repairs is allowed, after which, the system is brought to a halt for maintenance or for replacement. The usual concepts of stochastic stability are related to pure infinite horizon problems, and are not appropriate in this scenario. A new stability concept is introduced, named stochastic tau-stability that is tailored to the present setting. Necessary and sufficient conditions to ensure the stochastic tau-stability are provided, and the almost sure stability concept associated with this class of processes is also addressed. The paper also develops equivalences among second order concepts that parallels the results for infinite horizon problems. (C) 2003 Elsevier B.V. All rights reserved.

Modal analysis of anisotropic diffused-channel waveguide by a scalar finite element method

A procedure to model optical diffused-channel waveguides is presented in this work. The dielectric waveguides present anisotropic refractive indexes which are calculated from the proton concentration. The proton concentration inside the channel is calculated by the anisotropic 2D-linear diffusion equation and converted to the refractive indexes using mathematical relations obtained from experimental data, the arbitrary refractive index profile is modeled by a. nodal expansion in the base functions. The TE and TM-like propagation properties (effective index) and the electromagnetic fields for well-annealed proton-exchanged (APE) LiNbO3 waveguides are computed by the finite element method.

MATHEMATICA notebook for computing tetrahedral finite element shape functions and matrices for the Helmholtz equation

A MATHEMATICA notebook to compute the elements of the matrices which arise in the solution of the Helmholtz equation by the finite element method (nodal approximation) for tetrahedral elements of any approximation order is presented. The results of the notebook enable a fast computational implementation of finite element codes for high order simplex 3D elements reducing the overheads due to implementation and test of the complex mathematical expressions obtained from the analytical integrations. These matrices can be used in a large number of applications related to physical phenomena described by the Poisson, Laplace and Schrodinger equations with anisotropic physical properties.

Dynamics of the deconfinement transition of quarks and gluons in a finite volume

We employ finite elements methods for the approximation of solutions of the Ginzburg-Landau equations describing the deconfinement transition in quantum chromodynamics. These methods seem appropriate for situations where the deconfining transition occurs over a finite volume as in relativistic heavy ion collisions. where in addition expansion of the system and flow of matter are important. Simulation results employing finite elements are presented for a Ginzburg-Landau equation based on a model free energy describing the deconfining transition in pure gauge SU(2) theory. Results for finite and infinite system are compared. (C) 2009 Elsevier B.V. All rights reserved.

Finite element analysis of anisotropic optical waveguide with arbitrary index profile

This work presents the application of a scalar finite element formulation for Ex (TE-like) modes in anisotropic planar and channel waveguides with diagonal permittivity tensor, diffused in both transversal directions. This extended formulation considers explicitly both the variations of the refractive index and their spatial derivates inside of each finite element. Dispersion curves for Ex modes in planar and channel waveguides are shown, and the results compared with solutions obtained by other formulations.

Projeto de controladores robustos para aplicações em estruturas inteligentes utilizando desigualdades matriciais lineares

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

P-233-Biodiversity of mucosa associated microbiome of Crohn's disease patients living in middle western São Paulo State, Brazil

Background: The intestinal microbiome (IM) has extensively been studied in the search for a link of bacteria with the cause of Crohn`s disease (CD). The association might result from the action of a specific pathogen and/or an eventual imbalance in bacterial species composition of the gut. The innumerous virulence associated markers and strategies described for adherent and invasive Escherichia coli (AIEC) have made them putative candidate pathogens for CD. IM of CD patients shows dysbiosis, manifested by the proliferation of bacterial groups such as Enterobacteriaceae and reduction of others such as Lactobacillus and Bifidobacterium. The augmented bacterial population comprising of commensal and/or pathogenic organisms super stimulates the immune system, triggering the inflammatory reactions responsible for the clinical manifestations of the disease. Considering the role played by IM in CD and the multiple variables influencing its species composition, resulting in differences among populations, the objective of this study was to determine the bacterial biodiversity in the mucosa associated microbiome of CD patients from a population not previously subject to this analysis, living in the middle west region of Sao Paulo state. Methods: A total of 4 CD patients and 5 controls subjects attending the Botucatu Medical School of the Sao Paulo State University (UNESP) for routine colonoscopy and who signed an informed consent were included in the study. A number of 2 biopsies, one from the ileum and other from any part of the terminal colon, were taken from each subject and immediately frozen at -70[degrees]C until DNA purification. The bacterial biodiversity was assessed by next generation (ion torrent) sequencing of PCR amplicons of the ribosomal DNA 16S V6 region (16S V6 rDNA). The bacterial identification was performed at the genus level, by alignment of the generated DNA sequences with those available at the ribosomal database project (RDP) website. Results: The overall DNA sequence output was based on an average number of 526,427 reads per run, matching 50 bacterial genus 16SrDNA sequences available at the RDB website, and 22 non matching sequences. Over 95% of the sequences corresponded to taxa belonging to the major phyla: Firmicutes, Bacterioidetes, Proteobacteria and Actinobacteria. Irrespective of the intestinal site analyzed, no case-control differences could be observed in the prevalence of Actinobacteria and Firmicutes. The prevalence of Proteobacteria was higher (40%) in the biopsies of control subjects as compared to that of DC patients (16%). For Bacterioidetes, the higher prevalence was observed among DC patients (33% as opposed to 14,5% in controls). The significance for all comparisons considered a p value < 0,05 in a Chi2 test. No mucosal site specific differences could be observed in IM comparisons of CD and control subjects. Conclusions: The rise in the number of Bacterioidetes observed here among CD patients seems to be in agreement with most of studies published thus far. Yet, the reduction in the number of Proteobacteria along with an apparently unaltered population of Actinobacteria and Firmicutes, which include the so called "beneficial" organisms Bifidobacterium and Lactobacillus were rather surprising. These data suggest that the analyses on the role of IM in CD should consider the multiple variables that may influence its species composition.

A fully adaptive multiresolution scheme for shock computations

The scheme is based on Ami Harten's ideas (Harten, 1994), the main tools coming from wavelet theory, in the framework of multiresolution analysis for cell averages. But instead of evolving cell averages on the finest uniform level, we propose to evolve just the cell averages on the grid determined by the significant wavelet coefficients. Typically, there are few cells in each time step, big cells on smooth regions, and smaller ones close to irregularities of the solution. For the numerical flux, we use a simple uniform central finite difference scheme, adapted to the size of each cell. If any of the required neighboring cell averages is not present, it is interpolated from coarser scales. But we switch to ENO scheme in the finest part of the grids. To show the feasibility and efficiency of the method, it is applied to a system arising in polymer-flooding of an oil reservoir. In terms of CPU time and memory requirements, it outperforms Harten's multiresolution algorithm.The proposed method applies to systems of conservation laws in 1Dpartial derivative(t)u(x, t) + partial derivative(x)f(u(x, t)) = 0, u(x, t) is an element of R-m. (1)In the spirit of finite volume methods, we shall consider the explicit schemeupsilon(mu)(n+1) = upsilon(mu)(n) - Deltat/hmu ((f) over bar (mu) - (f) over bar (mu)-) = [Dupsilon(n)](mu), (2)where mu is a point of an irregular grid Gamma, mu(-) is the left neighbor of A in Gamma, upsilon(mu)(n) approximate to 1/mu-mu(-) integral(mu-)(mu) u(x, t(n))dx are approximated cell averages of the solution, (f) over bar (mu) = (f) over bar (mu)(upsilon(n)) are the numerical fluxes, and D is the numerical evolution operator of the scheme.According to the definition of (f) over bar (mu), several schemes of this type have been proposed and successfully applied (LeVeque, 1990). Godunov, Lax-Wendroff, and ENO are some of the popular names. Godunov scheme resolves well the shocks, but accuracy (of first order) is poor in smooth regions. Lax-Wendroff is of second order, but produces dangerous oscillations close to shocks. ENO schemes are good alternatives, with high order and without serious oscillations. But the price is high computational cost.Ami Harten proposed in (Harten, 1994) a simple strategy to save expensive ENO flux calculations. The basic tools come from multiresolution analysis for cell averages on uniform grids, and the principle is that wavelet coefficients can be used for the characterization of local smoothness.. Typically, only few wavelet coefficients are significant. At the finest level, they indicate discontinuity points, where ENO numerical fluxes are computed exactly. Elsewhere, cheaper fluxes can be safely used, or just interpolated from coarser scales. Different applications of this principle have been explored by several authors, see for example (G-Muller and Muller, 1998).Our scheme also uses Ami Harten's ideas. But instead of evolving the cell averages on the finest uniform level, we propose to evolve the cell averages on sparse grids associated with the significant wavelet coefficients. This means that the total number of cells is small, with big cells in smooth regions and smaller ones close to irregularities. This task requires improved new tools, which are described next.

Condição ocular dos indivíduos facectomizados na região centro-oeste do estado de São Paulo: estudo populacional

OBJETIVO: Avaliar as causas de baixa visão e cegueira em indivíduos facectomizados, de amostra da população de cidades da região centrooeste do estado de São Paulo. Métodos: Estudo transversal, observacional, feito em cinco cidades da região centro-oeste do estado de São Paulo, em amostra domiciliar e baseada nos dados do último Censo Demográfico (IBGE, 1995), com escolha sistemática dos domicílios. Foi considerada para o presente estudo uma subamostra de indivíduos facectomizados, dos quais foram obtidos dados de identificação e exame oftalmológico completo. Os dados foram avaliados por estatísticas descritivas, análise de freqüência de ocorrência e proporção de concordância, com intervalo de confiança de 95%. RESULTADOS: Dos indivíduos amostrados, 2,37% haviam sido submetidos à facectomia. Dos 201 olhos operados, 26,9% apresentavam acuidade visual compatível com cegueira ou deficiência visual. Com a melhor correção óptica, a acuidade visual permaneceu <0,3 em 19,0%. O exame refracional proporcionou melhora da acuidade visual para 27,9% dos indivíduos facectomizados. As causas de baixa visão foram os erros refrativos não corrigidos, opacidade de cápsula posterior (19,4%), ceratopatia bolhosa (8,3%) coriorretinite cicatricial (8,3%), afacia (8,3%), degeneração macular relacionada a idade (5,5%), leucoma (5,5%), glaucoma (5,5%), atrofia de papila (5,5,%), descolamento de retina (2,8%), atrofia de epitelio pigmentado da retina (2,8%) e alta miopia (2,8%). CONCLUSÃO: Apesar da catarata ser causa de cegueira que pode ser evitável, mesmo após a correção cirúrgica porcentagem expressiva de indivíduos permanece com baixa visão, em geral, em decorrência de fatores relacionados ao seguimento pós-operatório negligenciado.

Características de portadores de pterígio na região de Botucatu

OBJETIVO: Avaliar as características dos portadores de pterígio na região de Botucatu (SP). MÉTODOS: Portadores de pterígio foram avaliados quanto à idade, sexo, profissão, queixas e características da lesão (primário ou recidivado, tamanho, carnoso ou involutivo). Os dados foram submetidos à avaliação estatística. RESULTADOS: Cerca de metade dos portadores eram do sexo feminino, a maioria com idade superior a 40 anos, que procuraram o tratamento com queixa do efeito anti-estético da lesão. A maioria das lesões era primária (77,08%), grau 2 (69,6%) e do tipo carnoso (86,7%). Pterígio grau 4 esteve presente em 1,4% dos pacientes. CONCLUSÃO: Observou-se maior freqüência de portadores de pterígio entre mulheres, acima dos 40 anos e portadoras de pterígio primário, grau 2 e carnoso. A cegueira por pterígio ainda está presente em nosso meio.

Numerical approximation of the Ginzburg-Landau equation with memory effects in the dynamics of phase transitions

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A REACHABILITY TREE FOR STATECHARTS AND ANALYSIS OF SOME PROPERTIES

Statecharts are an extension to finite state machines with capability for expressing hierarchical decomposition and parallelism. They also have a mechanism called history, to remember the last visit to a superstate. An algorithm to create a reachability tree for statecharts is presented. Also shown is how to use this tree to analyse dynamic properties of statecharts; reachability from any state configuration, usage of transitions, reinitiability, deadlocks, and valid sequence of events. Owing to its powerful notation, building a reachability tree for statecharts presents some difficulties, and we show how these problems were solved in the tree we propose.

Stochastic Stability for Markovian Jump Linear Systems Subject to a Crucial Failure Event

This paper is concerned with the stability of discrete-time linear systems subject to random jumps in the parameters, described by an underlying finite-state Markov chain. In the model studied, a stopping time τ Δ is associated with the occurrence of a crucial failure after which the system is brought to a halt for maintenance. The usual stochastic stability concepts and associated results are not indicated, since they are tailored to pure infinite horizon problems. Using the concept named stochastic τ-stability, equivalent conditions to ensure the stochastic stability of the system until the occurrence of τ Δ is obtained. In addition, an intermediary and mixed case for which τ represents the minimum between the occurrence of a fix number N of failures and the occurrence of a crucial failure τ Δ is also considered. Necessary and sufficient conditions to ensure the stochastic τ-stability are provided in this setting that are auxiliary to the main result.

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. © 2008 IEEE.