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.

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.

Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is known about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they underestimate charge gaps.

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.

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.

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

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.

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.

Particle competition for complex network community detection

In many real situations, randomness is considered to be uncertainty or even confusion which impedes human beings from making a correct decision. Here we study the combined role of randomness and determinism in particle dynamics for complex network community detection. In the proposed model, particles walk in the network and compete with each other in such a way that each of them tries to possess as many nodes as possible. Moreover, we introduce a rule to adjust the level of randomness of particle walking in the network, and we have found that a portion of randomness can largely improve the community detection rate. Computer simulations show that the model has good community detection performance and at the same time presents low computational complexity. (C) 2008 American Institute of Physics.

Impact of Manaus City on the Amazon Green Ocean atmosphere: ozone production, precursor sensitivity and aerosol load

As a contribution to the Large-Scale Biosphere-Atmosphere Experiment in Amazonia - Cooperative LBA Airborne Regional Experiment (LBA-CLAIRE-2001) field campaign in the heart of the Amazon Basin, we analyzed the temporal and spatial dynamics of the urban plume of Manaus City during the wet-to-dry season transition period in July 2001. During the flights, we performed vertical stacks of crosswind transects in the urban outflow downwind of Manaus City, measuring a comprehensive set of trace constituents including O(3), NO, NO(2), CO, VOC, CO(2), and H(2)O. Aerosol loads were characterized by concentrations of total aerosol number (CN) and cloud condensation nuclei (CCN), and by light scattering properties. Measurements over pristine rainforest areas during the campaign showed low levels of pollution from biomass burning or industrial emissions, representative of wet season background conditions. The urban plume of Manaus City was found to be joined by plumes from power plants south of the city, all showing evidence of very strong photochemical ozone formation. One episode is discussed in detail, where a threefold increase in ozone mixing ratios within the atmospheric boundary layer occurred within a 100 km travel distance downwind of Manaus. Observation-based estimates of the ozone production rates in the plume reached 15 ppb h(-1). Within the plume core, aerosol concentrations were strongly enhanced, with Delta CN/Delta CO ratios about one order of magnitude higher than observed in Amazon biomass burning plumes. Delta CN/Delta CO ratios tended to decrease with increasing transport time, indicative of a significant reduction in particle number by coagulation, and without substantial new particle nucleation occurring within the time/space observed. While in the background atmosphere a large fraction of the total particle number served as CCN (about 60-80% at 0.6% supersaturation), the CCN/CN ratios within the plume indicated that only a small fraction (16 +/- 12 %) of the plume particles were CCN. The fresh plume aerosols showed relatively weak light scattering efficiency. The CO-normalized CCN concentrations and light scattering coefficients increased with plume age in most cases, suggesting particle growth by condensation of soluble organic or inorganic species. We used a Single Column Chemistry and Transport Model (SCM) to infer the urban pollution emission fluxes of Manaus City, implying observed mixing ratios of CO, NO(x) and VOC. The model can reproduce the temporal/spatial distribution of ozone enhancements in the Manaus plume, both with and without accounting for the distinct (high NO(x)) contribution by the power plants; this way examining the sensitivity of ozone production to changes in the emission rates of NO(x). The VOC reactivity in the Manaus region was dominated by a high burden of biogenic isoprene from the background rainforest atmosphere, and therefore NO(x) control is assumed to be the most effective ozone abatement strategy. Both observations and models show that the agglomeration of NO(x) emission sources, like power plants, in a well-arranged area can decrease the ozone production efficiency in the near field of the urban populated cores. But on the other hand remote areas downwind of the city then bear the brunt, being exposed to increased ozone production and N-deposition. The simulated maximum stomatal ozone uptake fluxes were 4 nmol m(-2) s(-1) close to Manaus, and decreased only to about 2 nmol m(-2) s(-1) within a travel distance >1500 km downwind from Manaus, clearly exceeding the critical threshold level for broadleaf trees. Likewise, the simulated N deposition close to Manaus was similar to 70 kg N ha(-1) a(-1) decreasing only to about 30 kg N ha(-1) a(-1) after three days of simulation.

Contribution of the second Landau level to the exchange energy of the three-dimensional electron gas in a high magnetic field

We derive a closed analytical expression for the exchange energy of the three-dimensional interacting electron gas in strong magnetic fields, which goes beyond the quantum limit (L=0) by explicitly including the effect of the second, L=1, Landau level and arbitrary spin polarization. The inclusion of the L=1 level brings the fields to which the formula applies closer to the laboratory range, as compared to previous expressions, valid only for L=0 and complete spin polarization. We identify and explain two distinct regimes separated by a critical density n(c). Below n(c), the per particle exchange energy is lowered by the contribution of L=1, whereas above n(c) it is increased. As special cases of our general equation we recover various known more limited results for higher fields, and we identify and correct a few inconsistencies in some of these earlier expressions.