A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

A geometric mass-preserving redistancing scheme for the level set function

In this paper we describe and evaluate a geometric mass-preserving redistancing procedure for the level set function on general structured grids. The proposed algorithm is adapted from a recent finite element-based method and preserves the mass by means of a localized mass correction. A salient feature of the scheme is the absence of adjustable parameters. The algorithm is tested in two and three spatial dimensions and compared with the widely used partial differential equation (PDE)-based redistancing method using structured Cartesian grids. Through the use of quantitative error measures of interest in level set methods, we show that the overall performance of the proposed geometric procedure is better than PDE-based reinitialization schemes, since it is more robust with comparable accuracy. We also show that the algorithm is well-suited for the highly stretched curvilinear grids used in CFD simulations. Copyright (C) 2010 John Wiley & Sons, Ltd.

Continuous GRASP with a local active-set method for bound-constrained global optimization

Global optimization seeks a minimum or maximum of a multimodal function over a discrete or continuous domain. In this paper, we propose a hybrid heuristic-based on the CGRASP and GENCAN methods-for finding approximate solutions for continuous global optimization problems subject to box constraints. Experimental results illustrate the relative effectiveness of CGRASP-GENCAN on a set of benchmark multimodal test functions.

Front tracking with moving-least-squares surfaces

The representation of interfaces by means of the algebraic moving-least-squares (AMLS) technique is addressed. This technique, in which the interface is represented by an unconnected set of points, is interesting for evolving fluid interfaces since there is]to surface connectivity. The position of the surface points can thus be updated without concerns about the quality of any surface triangulation. We introduce a novel AMLS technique especially designed for evolving-interfaces applications that we denote RAMLS (for Robust AMLS). The main advantages with respect to previous AMLS techniques are: increased robustness, computational efficiency, and being free of user-tuned parameters. Further, we propose a new front-tracking method based on the Lagrangian advection of the unconnected point set that defines the RAMLS surface. We assume that a background Eulerian grid is defined with some grid spacing h. The advection of the point set makes the surface evolve in time. The point cloud can be regenerated at any time (in particular, we regenerate it each time step) by intersecting the gridlines with the evolved surface, which guarantees that the density of points on the surface is always well balanced. The intersection algorithm is essentially a ray-tracing algorithm, well-studied in computer graphics, in which a line (ray) is traced so as to detect all intersections with a surface. Also, the tracing of each gridline is independent and can thus be performed in parallel. Several tests are reported assessing first the accuracy of the proposed RAMLS technique, and then of the front-tracking method based on it. Comparison with previous Eulerian, Lagrangian and hybrid techniques encourage further development of the proposed method for fluid mechanics applications. (C) 2008 Elsevier Inc. All rights reserved.

A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Qualidade nutricional das refeições servidas em uma unidade de alimentação e nutrição de uma indústria da região metropolitana de São Paulo

OBJETIVO: Avaliar a qualidade nutricional das refeições servidas em uma Unidade de Alimentação e Nutrição de uma fábrica da região metropolitana da cidade de São Paulo. MÉTODOS: Dentre os cardápios praticados no período de um ano (242 dias) na unidade mencionada, foram selecionados 30% por sorteio sistemático, e avaliados utilizando-se o Índice de Qualidade da Refeição, com base nas recomendações da Organização Mundial da Saúde e do Ministério da Saúde brasileiro. Esse índice compõe-se de cinco itens que variam entre zero e 20 pontos cada um: adequação na oferta de hortaliças e frutas; oferta de carboidratos; oferta de gordura total; oferta de gordura saturada e variabilidade do cardápio. No período analisado, foram servidas 367 preparações, agrupadas em 30 categorias, segundo composição e forma de preparo. A correlação de Spearman foi utilizada para investigar a correlação do índice com os nutrientes da refeição. As análises foram realizadas no pacote estatístico STATA, considerando-se o nível de significância de 5%. RESULTADOS: O valor médio do Índice de Qualidade da Refeição foi de 64,60 (DP=21,18) pontos, sendo 44% das refeições classificadas como "refeição que necessita de melhora" e apenas 25% como "adequadas". Além do arroz e do feijão, servidos diariamente, as preparações mais frequentes foram: legumes e frutas (30%), massas e cremes (12%), frituras (9%) e sobremesas com creme (8%). Encontrou-se correlação positiva entre o Índice de Qualidade da Refeição e a vitamina C (r=0,32). CONCLUSÃO: Apesar da presença constante de frutas, legumes e verduras, há a necessidade de adequar a oferta das preparações às recomendações para uma alimentação saudável, que efetivamente colaborem na promoção da saúde.

Exploring the knowledge contained in neuroimages: Statistical discriminant analysis and automatic segmentation of the most significant changes

Objective: The aim of this article is to propose an integrated framework for extracting and describing patterns of disorders from medical images using a combination of linear discriminant analysis and active contour models. Methods: A multivariate statistical methodology was first used to identify the most discriminating hyperplane separating two groups of images (from healthy controls and patients with schizophrenia) contained in the input data. After this, the present work makes explicit the differences found by the multivariate statistical method by subtracting the discriminant models of controls and patients, weighted by the pooled variance between the two groups. A variational level-set technique was used to segment clusters of these differences. We obtain a label of each anatomical change using the Talairach atlas. Results: In this work all the data was analysed simultaneously rather than assuming a priori regions of interest. As a consequence of this, by using active contour models, we were able to obtain regions of interest that were emergent from the data. The results were evaluated using, as gold standard, well-known facts about the neuroanatomical changes related to schizophrenia. Most of the items in the gold standard was covered in our result set. Conclusions: We argue that such investigation provides a suitable framework for characterising the high complexity of magnetic resonance images in schizophrenia as the results obtained indicate a high sensitivity rate with respect to the gold standard. (C) 2010 Elsevier B.V. All rights reserved.

Pulmonary function testing in quadriplegic subjects

Study design: Cross-sectional study. Objective: Pulmonary functional capacity in 23 Brazilian quadriplegic subjects (ASIA A), aged 30 (9.5) years, weight 66 (10.75) kg, height 176 (7) cm, was investigated at 42 ( 64) months postinjury. Setting: University Hospital-UNICAMP, Campinas, Brazil. Method: Subjects performed forced vital capacity ( FVC) and maximal voluntary ventilation (MVV) tests while seated in their standard wheelchairs. Forced Expired Volume after 1 s (FEV1) and FVC/FEV1 ratio were calculated from these tests. Values obtained were compared to three prediction equations from the literature that are used specifically for spinal cord subjects and include different variables in their formulae, such as age, gender, height, postinjury time and injury level. Data are expressed as median (interquartile interval). Differences between values were demonstrated by median confidence interval with significance level set at a 0.05. Results: Obtained data were statistically different from prediction equation results, with FVC 3.11 ( 0.81), 4.46 (0.28), 4.16 (0.33), 4.26 (0.42); FEV1 2.77 (1.03), 3.67 (0.21), 3.66 (0.30), 3.45 (0.39) and MVV 92 (27), 154.2 (11.9), 156.6 (14),157.3 (16.8), where the first value is obtained experimentally and the second, third and fourth values correspond to predicted values. The results obtained from spirometry test in this study differed significantly from the results obtained when prediction equations were used. Conclusion: The use of prediction equations developed to estimate pulmonary function in wheelchair users significantly overestimates pulmonary function of quadriplegic individuals with complete lesions (ASIA group A), in comparison to measured values.

Mercúrio total em cabelo de crianças de uma população costeira, Cananéia, São Paulo, Brasil

Os níveis de mercúrio (Hg) total em cabelos estão diretamente relacionados à alimentação, particularmente ao consumo de peixes por populações costeiras com grande representação caiçara. No presente estudo foram avaliados os níveis de mercúrio total em cabelos de crianças com idade entre 4 e 12 anos, pertencente a três escolas públicas da cidade de Cananéia, São Paulo, Brasil. Os resultados obtidos (mediana e intervalo) para mercúrio total foram de: 0,04mg.kg-1 (0,01-0,77mg.kg-1), 0,39mg.kg-1 (< 0,01-3,33mg.kg-1) e 0,39mg.kg-1 (< 0,01-2,81mg.kg-1) considerando as escolas ES1, ES2 e ES3, respectivamente. Em geral, os valores encontrados estiveram bem abaixo do valor preconizado pela Organização Mundial da Saúde para uma população adulta não exposta ao mercúrio (2,0mg.kg-1). Os baixos valores observados e a inexistência de valores de referência para mercúrio total em cabelos de crianças brasileiras possibilitam a consideração desses valores como possível referência nacional em cabelos de populações costeiras, uma vez que foram obtidos em região de baixo impacto ambiental.

Effects of phonophoresis with Arnica montana onto acute inflammatory process in rat skeletal muscles: An experimental study

This study aimed at verifying the effects of phonophoresis associated with Arnica montana on the acute phase of an inflammatory muscle lesion. Forty Wistar male rats (300 +/- 50 g), of which the Tibialis Anterior muscle was surgically lesioned, were divided into four groups (n = 10 each): control group received no treatment; the ultrasound group (US) was treated in pulsed mode with 1-MHz frequency, 0.5 W/cm(2) intensity (spatial and temporal average - SATA), duty cycle of 1: 2 (2 ms on, 4 ms off, 50%), time of application 3 min per session, one session per day, for 3 days; the phonophoresis or ultrasound plus arnica (US+A) group was treated with arnica with the same US parameters plus arnica gel; and the arnica group (A) was submitted to massage with arnica gel, also for 3 min, once a day, for 3 days. Treatment started 24 h after the surgical lesion. On the 4th day after lesion creation, animals were sacrificed and sections of the lesioned, inflamed muscle were removed for quantitative (mononuclear and polymorphonuclear cell count) and qualitative histological analysis. Collected data from the 4 groups were statistically analyzed and the significance level set at p < 0.05. Results show higher mononuclear cell density in all three treated groups with no significant difference between them, but values were significantly different (p < 0.0001) when compared to control group`s. As to polymorphonuclear cell density, significant differences were found between control group (p = 0.0134) and US, US+A and A groups; the arnica group presented lesser density of polymorphonuclear cells when compared (p = 0.0134) to the other groups. No significant difference was found between US and US+A groups. While the massage with arnica gel proved to be an effective anti-inflammatory on acute muscle lesion in topic use, these results point to ineffectiveness of Arnica montana phonophoresis, US having seemingly checked or minimized its anti-inflammatory effect. (C) 2008 Elsevier B. V. All rights reserved.

Genetic Diversity and a PCR-Based Method for Xanthomonas axonopodis Detection in Passion Fruit

Xanthomonas axonopodis pv. passiflorae causes bacterial spot in passion fruit. It attacks the purple and yellow passion fruit as well as the sweet passion fruit. The diversity of 87 isolates of pv. passiflorae collected from across 22 fruit orchards in Brazil was evaluated using molecular profiles and statistical procedures, including an unweighted pair-group method with arithmetical averages-based dendrogram, analysis of molecular variance (AMOVA), and an assigning test that provides information on genetic structure at the population level. Isolates from another eight pathovars were included in the molecular analyses and all were shown to have a distinct repetitive sequence-based polymerase chain reaction profile. Amplified fragment length polymorphism technique revealed considerable diversity among isolates of pv. passiflorae, and AMOVA showed that most of the variance (49.4%) was due to differences between localities. Cluster analysis revealed that most genotypic clusters were homogeneous and that variance was associated primarily with geographic origin. The disease adversely affects fruit production and may kill infected plants. A method for rapid diagnosis of the pathogen, even before the disease symptoms become evident, has value for producers. Here, a set of primers (Xapas) was designed by exploiting a single-nucleotide polymorphism between the sequences of the intergenic 16S-23S rRNA spacer region of the pathovars. Xapas was shown to effectively detect all pv. passiflorae isolates and is recommended for disease diagnosis in passion fruit orchards.

Hyperplane navigation: A method to set individual scores in fMRI group datasets

Recent studies have demonstrated that spatial patterns of fMRI BOLD activity distribution over the brain may be used to classify different groups or mental states. These studies are based on the application of advanced pattern recognition approaches and multivariate statistical classifiers. Most published articles in this field are focused on improving the accuracy rates and many approaches have been proposed to accomplish this task. Nevertheless, a point inherent to most machine learning methods (and still relatively unexplored in neuroimaging) is how the discriminative information can be used to characterize groups and their differences. In this work, we introduce the Maximum Uncertainty Linear Discrimination Analysis (MLDA) and show how it can be applied to infer groups` patterns by discriminant hyperplane navigation. In addition, we show that it naturally defines a behavioral score, i.e., an index quantifying the distance between the states of a subject from predefined groups. We validate and illustrate this approach using a motor block design fMRI experiment data with 35 subjects. (C) 2008 Elsevier Inc. All rights reserved.

Partial spectral projected gradient method with active-set strategy for linearly constrained optimization

A method for linearly constrained optimization which modifies and generalizes recent box-constraint optimization algorithms is introduced. The new algorithm is based on a relaxed form of Spectral Projected Gradient iterations. Intercalated with these projected steps, internal iterations restricted to faces of the polytope are performed, which enhance the efficiency of the algorithm. Convergence proofs are given and numerical experiments are included and commented. Software supporting this paper is available through the Tango Project web page: http://www.ime.usp.br/similar to egbirgin/tango/.

Two- and three-dimensional shape fabric analysis by the intercept method in grey levels

The count intercept is a robust method for the numerical analysis of fabrics Launeau and Robin (1996). It counts the number of intersections between a set of parallel scan lines and a mineral phase, which must be identified on a digital image. However, the method is only sensitive to boundaries and therefore supposes the user has some knowledge about their significance. The aim of this paper is to show that a proper grey level detection of boundaries along scan lines is sufficient to calculate the two-dimensional anisotropy of grain or crystal distributions without any particular image processing. Populations of grains and crystals usually display elliptical anisotropies in rocks. When confirmed by the intercept analysis, a combination of a minimum of 3 mean length intercept roses, taken on 3 more or less perpendicular sections, allows the calculation of 3-dimensional ellipsoids and the determination of their standard deviation with direction and intensity in 3 dimensions as well. The feasibility of this quick method is attested by numerous examples on theoretical objects deformed by active and passive deformation, on BSE images of synthetic magma flow, on drawing or direct analysis of thin section pictures of sandstones and on digital images of granites directly taken and measured in the field. (C) 2010 Elsevier B.V. All rights reserved.