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.

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 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.

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.

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.

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.

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.

Mapping low- and high-density clouds in astrophysical nebulae by imaging forbidden line emission

Emission line ratios have been essential for determining physical parameters such as gas temperature and density in astrophysical gaseous nebulae. With the advent of panoramic spectroscopic devices, images of regions with emission lines related to these physical parameters can, in principle, also be produced. We show that, with observations from modern instruments, it is possible to transform images taken from density-sensitive forbidden lines into images of emission from high- and low-density clouds by applying a transformation matrix. In order to achieve this, images of the pairs of density-sensitive lines as well as the adjacent continuum have to be observed and combined. We have computed the critical densities for a series of pairs of lines in the infrared, optical, ultraviolet and X-rays bands, and calculated the pair line intensity ratios in the high- and low-density limit using a four- and five-level atom approximation. In order to illustrate the method, we applied it to Gemini Multi-Object Spectrograph (GMOS) Integral Field Unit (GMOS-IFU) data of two galactic nuclei. We conclude that this method provides new information of astrophysical interest, especially for mapping low- and high-density clouds; for this reason, we call it `the ld/hd imaging method`.

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.

Comparison Between Single-Diode Low-Level Laser Therapy (LLLT) and LED Multi-Diode (Cluster) Therapy (LEDT) Applications Before High-Intensity Exercise

Background Data and Objective: There is anecdotal evidence that low-level laser therapy (LLLT) may affect the development of muscular fatigue, minor muscle damage, and recovery after heavy exercises. Although manufacturers claim that cluster probes (LEDT) maybe more effective than single-diode lasers in clinical settings, there is a lack of head-to-head comparisons in controlled trials. This study was designed to compare the effect of single-diode LLLT and cluster LEDT before heavy exercise. Materials and Methods: This was a randomized, placebo-controlled, double-blind cross-over study. Young male volleyball players (n = 8) were enrolled and asked to perform three Wingate cycle tests after 4 x 30 sec LLLT or LEDT pretreatment of the rectus femoris muscle with either (1) an active LEDT cluster-probe (660/850 nm, 10/30mW), (2) a placebo cluster-probe with no output, and (3) a single-diode 810-nm 200-mW laser. Results: The active LEDT group had significantly decreased post-exercise creatine kinase (CK) levels (-18.88 +/- 41.48U/L), compared to the placebo cluster group (26.88 +/- 15.18U/L) (p < 0.05) and the active single-diode laser group (43.38 +/- 32.90U/L) (p<0.01). None of the pre-exercise LLLT or LEDT protocols enhanced performance on the Wingate tests or reduced post-exercise blood lactate levels. However, a non-significant tendency toward lower post-exercise blood lactate levels in the treated groups should be explored further. Conclusion: In this experimental set-up, only the active LEDT probe decreased post-exercise CK levels after the Wingate cycle test. Neither performance nor blood lactate levels were significantly affected by this protocol of pre-exercise LEDT or LLLT.

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

Individual-level diet variation in four species of Brazilian frogs

Many natural populations exploiting a wide range of resources are actually composed of relatively specialized individuals. This interindividual variation is thought to be a consequence of the invasion of `empty` niches in depauperate communities, generally in temperate regions. If individual niches are constrained by functional trade-offs, the expansion of the population niche is only achieved by an increase in interindividual variation, consistent with the `niche variation hypothesis`. According to this hypothesis, we should not expect interindividual variation in species belonging to highly diverse, packed communities. In the present study, we measured the degree of interindividual diet variation in four species of frogs of the highly diverse Brazilian Cerrado, using both gut contents and delta(13)C stable isotopes. We found evidence of significant diet variation in the four species, indicating that this phenomenon is not restricted to depauperate communities in temperate regions. The lack of correlations between the frogs` morphology and diet indicate that trade-offs do not depend on the morphological characters measured here and are probably not biomechanical. The nature of the trade-offs remains unknown, but are likely to be cognitive or physiological. Finally, we found a positive correlation between the population niche width and the degree of diet variation, but a null model showed that this correlation can be generated by individuals sampling randomly from a common set of resources. Therefore, albeit consistent with, our results cannot be taken as evidence in favour of the niche variation hypothesis.

Pricing farm-level agricultural insurance: a Bayesian approach

This paper applies Hierarchical Bayesian Models to price farm-level yield insurance contracts. This methodology considers the temporal effect, the spatial dependence and spatio-temporal models. One of the major advantages of this framework is that an estimate of the premium rate is obtained directly from the posterior distribution. These methods were applied to a farm-level data set of soybean in the State of the Parana (Brazil), for the period between 1994 and 2003. The model selection was based on a posterior predictive criterion. This study improves considerably the estimation of the fair premium rates considering the small number of observations.